成人颈关节结构与微有限元力学论文范文

［摘要］目的分析成人颈椎关节突不同区域、C2～C7椎序间、左右侧别微结构的变化规律及其与关节突微有限元力学反应之间的关系，预测及评价由关节突微观形态及力学载荷的改变导致颈椎失稳、骨质疏松、关节突骨折等引起关节突损伤等颈椎疾患的风险。方法取由内蒙古医科大学解剖学教研室提供的36个成人新鲜颈椎上、下关节突，行Micro-CT扫描，手动勾选上、下关节突2个感兴趣区域(VOI)，并计算关节突微结构参数;用环钻从关节突取体积为8mm×8mm×8mm的立方体进行微有限元力学分析，对比分析关节突微结构与力学参数之间的关系。结果主效应分析上、下关节突区域间BV/TV、BS/BV、Tb．Sp参数之间差异有统计学意义(P＜0．05);关节突左、右侧别间，仅Tb．Th参数值差异有统计学意义(P＜0．05);不同椎序间关节突BV/TV、Tb．N、BS/BV、Tb．Th、Tb．Sp、Tb．Pf参数值差异均有统计学意义(P＜0．05);交互效应仅在不同椎序间与上、下关节突区域间交互效应Tb．N参数值差异有统计学意义(P＜0．05)。微有限元力学表现为关节突骨小梁微有限元模型所受的六个方向的应力分量为σx、σy、σz、σxy、σxz、σyz，分别对应的应力值为S11、S22、S33、S12、S13、S23，其中S33即在σz轴向应力值最大，符合关节突加载载荷的方向(沿Z轴方向加载)，S12即σxy二维平面所受应力最小，表明颈椎在运动时关节突所受切应力小于主应力。结论C3关节突BV/TV、Tb．N参数值最大，BV/TV、Tb．Th与关节突应力呈正相关，与骨折程度、受力集中程度、内部剪切力呈负相关，因此，C3关节突载荷能力较好，发生骨折等损伤风险相对较小。C2关节突BS/BV、Tb．Sp、Tb．Pf参数值最大，Tb．Pf、Tb．Sp与应力呈负相关，与骨折、骨质疏松、增生、退变等呈正相关，说明C2关节突发生骨折、骨质疏松、增生、退变等颈椎疾病风险较其他关节突大。关节突骨小梁有限元模型在整个变形过程中塑性应变累计结果为零，即在此载荷条件下未发生塑性变形，没有发生屈服，说明在成人正常生理载荷条件下关节突骨小梁有限元模型未发生骨折。

［关键词］颈椎;关节突;微结构;微有限元力学;Micro-CT

颈椎关节突增生、退变是引起颈椎不稳和神经根受压的重要原因，被喻为慢性颈痛的发动机。关节突不规则的增生和肥大等因素可使关节突形态改变，最终导致关节突结构形态和受力状态改变，颈椎稳定性丧失，因此颈椎关节突区域微结构与微有限元力学研究对临床研究和治疗相关疾病具有重要意义。但鲜有报道成人颈椎关节突区域微结构参数预测表观力学、致病属性的对比研究。关节突的解剖特征是其发生退变与创伤性关节病变的解剖学基础。为适应机体机械力学环境，承重骨和非承重骨在结构和微结构上存在差异，骨质量和结构总是以适应力学载荷的形式存在，椎骨区域微结构与其应力存在相关性(屈服准则)［1-3］。运动员整体胫骨皮质骨和小梁骨的密度和微结构值在应力断裂和非应力断裂之间无显著差异，区域分析显示两者之间存在显著差异［4］。胫骨前、后、左、右区域微结构因承力不同而异［5］。椎体背侧骨小梁体积、骨小梁体积分数、骨小梁连接度均大于腹侧，且背侧骨小梁分离度小于腹侧，椎体区域间BV/TV、连接密度和屈服强度之间有很强的相关性［6］。基于以上文献报道采用同一骨质分区方法探讨同一部位不同区域微结构及力学变化的方法，本研究通过对比分析关节突区域微结构参数(BV/TV、Tb．N等)和微有限元力学参数之间的相关性，探讨关节突2个VOI微结构参数与微有限元力学参数的关系，为有效评估颈椎关节突因节段应力过于集中导致退行性改变及损伤的风险，现报告如下。

1材料与方法

1．1材料

选取由内蒙古医科大学解剖教研室提供的新鲜成人颈椎上、下关节突共36个标本(均为男性，平均年龄52岁)，剔除肌肉、横突、棘突保留双侧完整关节突(C2～C7左、右侧关节突)行Micro-CT扫描(内蒙古医科大学附属医院核医学科)。数据来源于内蒙古医科大学数字医学中心;扫描设备依托内蒙古医科大学附属医院核医学科Micro-CT扫描仪(SiemensInveonPET/CT，德国)，扫描参数为80kV设备电压，500mA电流，螺旋扫描模式，X射线发射源为可变焦发射源，焦点直径小于6μm，最大分辨率为16．75μm，扫描层数为1024pixels，扫描图像以CAT格式储存。

1．2方法

利用Micro-CT扫描仪自带的InveonResearchWorkplace三维重建软件(由内蒙古医科大学附属医院核医学科提供)进行后处理。将Micro-CT连续扫描的关节突断层原始数据以“无破损及压缩”的方式导入InveonResearchWorkplace，进入主界面选择“Gener-alAnalysis”选项，显示出关节突横断面、冠状面、矢状面三个平面的二维图像，在三个平面图像中均选择“512pixels”像素大小，然后在主界面选择Segmenta-tion＞Orientation＞XZPlane”，在关节突横断面二维图像中央区上均匀选取2个(Superior、Inferior)感兴趣区间体积(tissuevolume，TV)，TV厚度为100pixels(pix-el=16．75μm)，横断面TV面积为300×300pixels，最终获得TV为300×300×100pixels的立方体(TV=36．864mm3)，测量感兴趣区间内骨小梁相关参数，观察关节突的显微形态学结构，做出分析描述(图1)。骨小梁各结构参数及其意义见表1。

1．3统计学处理

采用SPSS16．0软件对数据进行统计学分析，应用Bartlett检验法进行方差齐性检验。关节突椎序间、左右侧别、四个区域间对微结构参数的影响用多因素方差分析(两两多重比较)，确立检验水准为α=0．05，以P＜0．05为差异有统计学意义。

2结果

2．1微结构参数结果

由上、下关节突不同区域、左右侧别、C2～C7不同椎序间多重比较微结构参数进行对比主效应及交互效应分析如下，其中主效应:①上、下关节突区域间BV/TV、BS/BV、Tb．Sp参数之间差异有统计学意义;②关节突左、右侧别间，仅Tb．Th参数值差异有统计学意义;③不同椎序间关节突BV/TV、Tb．N、BS/BV、Tb．Th、Tb．Sp、Tb．Pf参数值差异均有统计学意义(P＜0．05)。分别对①②、①③、②③三者间两两交互效应分析，结果显示，左右两侧的上、下关节突区域内关节突微结构参数值交互效应差异均无统计学意义(P＞0．05);不同椎序间与上、下关节突区域间交互效应仅Tb．N参数值差异有统计学意义(P＜0．05);不同椎序间与左、右侧别交互效应各微结构参数值差异均无统计学意义(P＞0．05)，见表2。其中，BV/TV、Tb．N参数值C3关节突最大，BV/TV、Tb．Th与关节突应力呈正相关，与骨折程度、受力集中程度，内部剪切力呈负相关，因此，C3关节突负载能力最大，发生骨折等损伤风险较小。BS/BV、Tb．Sp、Tb．Pf参数值C2关节突最大。Tb．Pf、Tb．Sp与应力呈负相关，与骨折、骨质疏松、增生、退变等呈正相关，说明C2关节突发生骨折风险、骨质疏松、增生、退变等颈椎疾病概率较其他关节突大。

2．2微有限元力学结果

通过对微有限元输出结果分析，关节突75%Mises最大应力值为1．497×104Pa，最小值为1．957×103Pa(图2a～c)，说明此关节突骨小梁模型所能承受的最大主应力为1．497×104Pa。若应力超过最大主应力，在法线方向就会超过抗压强度而开裂，开裂的方向就是法线方向(力学加载方向即Z轴方向)，以此预测关节突骨折的风险及骨折方向;关节突骨小梁模型所受的六个方向的应力分量为σx、σy、σz、σxy、σxz、σyz，其中σx、σy、σz分别为坐标系X、Y、Z轴方向的应力分量，分别对应的应力值为S11、S22、S33;σxy、σxz、σyz分别为坐标系XY、XZ、YZ平面方向的应力分量，分别对应的应力值为S12、S13、S23(图2d～i)，其中S33即在σz轴向应力值最大，符合关节突加载方向的载荷(沿Z轴方向加载)，S12即σxy二维平面所受应力最小，表明颈椎在运动时关节突主应力大于其他方向的切应力。图2j为最大弹性应变应力云图，在最小弹性应力为1．487×104～3．060×101Pa时，关节突骨小梁模型进入弹性应变状态，应力超过最大弹性应力1．760×104Pa，关节突骨小梁模型变为进入塑性应变状态(图2k、l)。等效塑性应力应变云图显示等效塑性应变(PEEQ)值均为零，表明在主应力间的最大差值为2．497×104Pa时，整个关节突骨小梁有限元模型在整个变形过程中塑性应变累计结果为零，即在此载荷条件下未发生塑性应变，没有发生屈服，即在此生理载荷条件下关节突骨小梁有限元模型未发生骨折现象(图2m、n)。

3讨论

1983年，Denis提出了脊柱的三柱理论，后柱承担颈椎运动载荷的20%，关节突属于后柱的主要应力结构，为颈椎损伤的治疗奠定了生物力学基础［7-8］。本研究结果表明颈椎关节突作为颈椎仅次于椎体的主要承重部位，上下关节突、不同椎序间、左右关节突对整个颈椎生理活动过程中的载荷贡献都不尽相同。C3关节突负载能力最大，发生骨折等损伤风险较小，其原因可能是C3作为上、下颈椎过渡椎体，骨结构形态总是适应其功能活动，从力学角度分析，其关节突起到承上启下的作用，既要承载头颅及上颈椎的载荷，还要维持下颈椎的生理曲度及稳定性，因此C3关节突承载能力相对较好。Tb．Sp、Tb．Pf参数值C2关节突最大。Tb．Pf、Tb．Sp与应力呈负相关，与骨折、骨质疏松、增生、退变等呈正相关，说明C2关节突发生骨折风险、骨质疏松、增生、退变等颈椎疾病概率较其他关节突大，C2作为上颈椎结构特殊椎体，寰齿关节主要完成上颈椎旋转屈曲等活动，活动范围较大，容易受到剪切应力而损伤，发生齿突骨折、寰齿关节脱位、关节突骨折等损伤。颈椎在运动时关节突在主应力大于其他方向的切应力，且当颈椎关节突累计受力或瞬间暴力超过弹性变形就会进入塑性变形阶段，在塑性阶段应力值超过关节突所承受应力阈值后就会加大关节突损伤风险。颈椎关节突形态结构的改变通常使颈椎三柱遭到破坏，严重丧失稳定性，并伴有不同程度的脊髓损伤，严重者可导致全身瘫痪甚至危及生命。颈椎的各关节突关节承受着颈椎一定范围内的前屈、后伸、压缩、牵拉、剪切、扭转等负荷，此结构特征是颈椎的运动枢纽，对维持颈椎稳定及正常生理活动起着至关重要的作用［9］。在颈椎脊髓运动节段(C2～C7)，关节突受压缩、剪切和轴向载荷的影响，关节突损伤可能与椎间盘退变有关，关节突、韧带及其囊的完整性对颈椎屈曲、伸展、侧弯和旋转的稳定性至关重要，颈椎上、下关节突关节面曲度、深度、关节突骨质与其承载能力相关。本研究经Micro-CT观测成人颈椎关节突不同侧别、椎序、上下关节突微结构的变化不尽相同，加载关节突微有限元模型，分析关节突的受力情况及应力阈值以评估颈椎关节突损伤的风险。颈椎关节突的异常改变包括先天性或创伤性、单侧或双侧的关节突脱位、关节突绞锁、关节突伴脱位外伤性骨折、关节突不对称、关节突炎等，可造成关节突软骨下骨及骨髓强化、重塑等微结构改变及力学能力的改变［10］。依据Wolff定律，骨质量和结构总是以适应力学载荷的形式存在，在力学试验和微有限元分析中，通常以加载条件、应力应变来评估小梁结构［11-12］。然而，骨小梁模型可以用来预测小梁骨的表观屈服行为［13］。虽然本研究没有模拟屈服后的破坏过程，但对屈服阶段(即用于评估骨强度的第一阶段)进行了模拟。已有研究采用双线性组织本构模型描述骨组织材料的非线性，其组织压缩应变和拉伸屈服应变分别为0．8%和0．48%，且多数脆性骨折是直接由骨质疏松引起，通常发生在骨小梁(松质骨)占优势的部位［14-15］。小梁骨在负荷传递和能量交换中起着重要的作用，主要以膝关节、髋关节和脊柱分布较丰富，其在椎体中承载超过75%的载荷，除了骨体积分数(骨组织体积与整体体积之比)外，还有详细的显微结构，包括小梁形态(棒状和板状骨小梁)和小梁方向［16-17］。骨小梁连通性对于控制骨小梁的机械性能也是重要的，基于骨小梁的形态参数，如棒状骨体积分数(BV/TV)和板状骨小梁体积分数与骨载荷能力、弹性模量和屈服强度高度相关［18-19］。从关节突微结构参数对比分析到微有限元力学分析颈椎关节突的结构功能及力学反应鲜有报道，本研究结合成人颈椎关节突微结构的观测与微有限元力学的分析，多重比较C2～C7上、下关节突骨小梁的分布与受力，有效评价颈椎关节突对整个颈椎运动的结构载荷与力学贡献，为人体颈椎生理活动防护、临床诊疗提供重要参考资料。

作者：冯会梅 王星 马渊 李筱贺 李志军 张少杰 单位：内蒙古医科大学基础医学院解剖学教研室