目的 利用有限元方法建立椎体前缘不同压缩状态模型，并分析椎体前缘压缩程度与邻椎终板应力的相关性。
方法 基于正常成年人志愿者的胸腰段CT数据，导入Mimics，经图像分割、修补及三维网格化获取胸腰椎T₁₁～L₁3个节段数据，网格赋材质属性后导入ABAQUS软件，获取逼真的三维有限元模型。对模型进行生理载荷加载，验证正常模型的有效性。将L1模型底部全自由度固定，在T₁₁模型上方施加前屈、后伸及轴向压缩载荷，分别将椎体前沿压缩至90%、80%、70%......10%9个状态，提取T11下终板及L1上终板的MISES应力，获得了压缩状态与邻近终板应力关系曲线。
结果 该模型高度逼真，能真实反映实际受力状态；T₁₁椎体下位终板和L1椎体上位终板的应力值与T₁₂椎体（研究对象）前沿的压缩程度正相关。
结论 应力增加可能导致终板骨折可能性增加，进而增大相邻椎体的骨折风险，有限元分析结果能够为临床实践提供依据。
关键词：胸腰椎；有限元分析；ABAQUS；三维模型；压缩状态； 终板 应力
脊柱压缩性骨折是老年性骨质疏松常见的并发症，主要发生在胸腰段椎体，以T₁₂椎体最为多见。椎体压缩性骨折由于椎体前缘的高度丢失，可导致局部后凸畸形、生理曲线发生改变，继而使骨折椎体前方应力集中及骨折椎体周围应力增大，从而导致邻近椎体再骨折的发生。有研究表明，椎体压缩性骨折即使接受了骨水泥强化治疗，邻近椎体的再骨折发生率仍可高达23%，可以设想椎体压缩骨折后若脊柱生物力学改变无法纠正，仍可增大邻近骨折的风险。本研究通过有限元方法[2]，探讨T₁₂椎体前缘不同压缩程度对邻近椎体终板应力及其与邻近椎体再骨折的相关性，为临床评估病情及指导治疗提供参考。
1资料与方法
1.1建立胸腰椎T₁₁～L1节段三维有限元模型
1.1.1材料选择
1例成年健康男性志愿者，身高173cm，体重72kg。拍摄X线检查并排除病变、创伤、畸形等情况。采用螺旋CT，连续扫描胸腰椎T10～L2节段，获取连续断层图像，层厚0.625mm，导出*.dcm格式文件，并保存。采用Materialise上海公司提供的Mimics16.0医学图像处理软件，导入上述扫描的断层图像，通过阀值区分，区域增长，模型修补等一系列必要操作，重建胸腰椎T₁₁～L1三维模型。采用PRO/E5.0结构设计软件对胸腰椎T₁₁～L1三维模型进行组装。将三维模型导入达索公司的ABAQU6.13-1有限元分析软件，通过材质属性赋值，网格划分，模型组装，完成三维有限元模型。
1.1.2方法
系统结构：椎骨区分皮质骨与松质骨，皮质骨厚1.5mm；模型包括黄韧带、翼状韧带、后纵韧带、小关节和椎间盘等。骨质模型采用高阶四面体单元，韧带等简化为采用truss绳索单元，椎间盘采用高阶六面体单元。所有结构均采用均质、连续、线弹性有限元材料。椎骨与韧带采用共节点约束，小关节之间采用contact接触关系，接触摩擦系数为0.1。胸腰椎T₁₁～L1三维模型结构的具体材质属性、单元类型以及单元数量见附表。
1.2正常胸腰椎节段有限元模型有效性验证
依据正常生理条件下胸腰椎所承受的载荷，对1.1节中建立的正常胸腰椎T₁₁～L1有限元模型进行载荷施加以及边界约束，载荷条件为在T11椎骨上表面一点施加15N·m扭矩载荷，T1底面全自由度固定约束。进行前屈、后伸、左旋、右旋、左侧屈、右侧屈6个方向运动，验证胸腰椎运动自由度数据，并与BROWN[6]、MARKOLF[7]、TENCER[8]及SCHLTZ等[9]的离体胸腰椎生物力学测试结果做对比。
1.3胸腰椎T₁₁～L1节段压缩有限元模型分析
单纯研究T₁₂椎体前缘压缩程度与邻近椎体终板应力的相关性，获取T₁₂椎体前缘压缩程度-邻椎终板应力峰值曲线图。用以分析老年性骨折疏松性T12椎体楔形压缩性骨折后，T₁₂椎体前缘压缩程度是否为邻近椎体再骨折的危险因素。以T12椎体前缘高度/后缘高度的数值建立不同前缘压缩状态模型，每10%为1个状态，依次建立90%、80%、70%......10%9个前缘压缩状态模型；在三维有限元模型上对L1下表面进行全自由度固定约束，在T₁₁上缘施加1.2节中的载荷；9个状态模型分别模拟实体生物力学进行生理载荷下前屈、后伸及轴向压缩运动；获取T₁₂椎体下终板及L1椎体上终板应力数据；获取T₁₂椎体前缘压缩程度与邻近椎体终板应力的相关性。
2结果
2.1正常胸腰椎T11～L1节段三维有限元模型有效性验证
正常胸腰椎T11～L1节段有限元模型包括前后纵韧带、棘间韧带、小关节等解剖结构，与真实胸腰椎模型具有较高的拟合度（见图1）。
通过1.2节中的载荷、边界，该正常模型在前屈、后伸、左右旋转、左右侧屈的结果与BROWN[6]、MARKOLF[7]、TENCER[8]及SCHLTZ等[9]的离体实验结果相吻合。证明该模型可以用于下一阶段的胸腰椎T11～L1压缩状态的生物力学有限元分析。该模型的分析结果见图2、3。
计算结果表明，本文研究的正常模型在轴向位移、伸直力矩、侧屈力矩以及扭转力矩等载荷条件下，与文献中的结果基本相似，随着载荷的增大，均呈现增大趋势，结果趋势与文献一致。
2.2胸腰椎T₁₁～L1节段压缩三维有限元模型分析
在正常胸腰椎T₁₁～L1节段椎骨有限元模型验证的基础上，在ABAQUS6.12-1中，以T12椎体前缘高度/后缘高度的数值建立不同前缘压缩状态模型，每10%为1个状态，建立出了9个工况的前缘压缩状态模型，提取出了前屈、后伸及旋转工况下T₁₁下终板及L1上终板的应力大小及分布情况，分别列出各状态终板的应力对比数据，见图4～9。
前屈载荷工况下，9种状态结果表明，当压缩比例从10%增加到90%时，T₁₁下终板应力从0.812MPa提高到4.019MPa；L1上终板应力从0.6141MPa提高到3.056MPa。后伸载荷工况下，9种状态结果表明，当压缩比例从10%增加到90%时，T₁₁下终板应力从1.151MPa提高到5.827MPa；L1上终板应力从0.6051MPa提高到2.908MPa。
轴向压缩载荷工况下，9种状态结果表明，当压缩比例从10%增加到90%时，T₁₁下终板应力从1.431MPa提高到6.333MPa；L1上终板应力从0.847MPa提高到3.535MPa。胸腰椎T₁₁～L1系统有限元模型中，椎骨与椎间盘之间，椎骨之间韧带连接等，模拟真实脊柱内的各部分连接状况，真实反映椎间盘的受力状态。提取各载荷条件下9种状态的T₁₁下终板及L1上终板的MISES应力。前屈90%状态最大MISES应力T₁₁下终板为4.019MPa，L1上终板为3.056MPa；后伸90%状态最大MISES应力T₁₁下终板为5.827MPa，L1上终板为2.908MPa；轴向压缩90%状态最大MISES应力T₁₁下终板为6.333MPa，L1上终板为3.535MPa。随着压缩状态的增大，椎间盘受到的压迫力更大，其椎间盘表现出来的应力水平随着压缩程度的增大而增大。
计算结果表明，随着压缩状态的增大，椎间盘受到的压迫力更大，其椎间盘表现出来的应力水平随着压缩程度的增大而增大。T₁₁下终板应力及L1上终板应力表现出的趋势相同。见图10。
3讨论
3.1正常胸腰椎T11～L1有限元模型有效性验证
本文研究胸腰椎T₁₁～L1中计算结果表明，随着压缩状态的增大，椎间盘受到的压迫力更大，其椎间盘表现出来的应力水平随着压缩程度的增大而增大。T₁₁下终板应力及L1上终板应力表现出的趋势相同系统三维有限元模型，该模型包括胸腰椎T₁₁～L1椎骨、椎间盘终板、髓核、纤维环以及各部分牵连韧带。依据韧带的相关特性以及材质属性，对于韧带在有限元中，采用非线性TRUSS单元来模拟韧带的特性，只承受拉力，不承受压力，实现了韧带的功能；椎骨与椎间盘之间设定为绑定关系[10-12]。要进行有限元分析，需先保证模型准确，建模方法可靠。在此基础上才能保证由该建模方法获得的模型是可靠的，为分析结果的可靠性提供最基础的条件。参考国外学者研究方法[13-14]，对正常胸腰椎T₁₁～L1有限元模型的有效性验证方法是：将模型定义和加载与标本实验相同的边界条件和载荷，分别进行了前屈、后伸、侧屈以及轴向旋转等运动模拟，然后将有限元结果与实验结果BROWN等[6]的研究做比较，以验证模型有效性。
结果表明，本文建立的正常胸腰椎T₁₁～L1节段椎骨及椎间盘、韧带系统三维有限元模型的有效性验证结果与参考文献结果相吻合，能够准确无误地呈现胸腰椎的受力状态，充分证明能够用于下一步的胸腰椎压缩状态有限元分析。
3.2压缩状态胸腰椎T₁₁～L1系统三维有限元模型分析
通过对椎体有限元模型前屈、后伸及轴向载荷的分析研究，获取的T₁₁下终板及L1上终板的MIS-ES应力大小及分布，应力结果显示，椎体前缘压缩程度越高，其终板应力越大。通过研究T₁₂椎体前缘压缩程度与邻椎终板应力的相关性，可以通过数据分析得知，老年性骨折疏松性T₁₂椎体楔形压缩性骨折后，椎体压缩程度为邻近椎体再骨折的危险因素，对指导临床治疗具有较大的意义。T₁₂椎体前缘压缩程度与邻近椎体终板应力的具有正相关性，结果表明，T₁₂椎体前缘压缩程度为邻近椎体再骨折的危险因素；通过精确建模手段及三维有限元计算方法，对于指导临床，降低治疗后邻近椎体再骨折的发生机率具有重要意义[15]。
本文有限元结果表明，在90%～60%区间，椎间盘能够较好的保持其形态，但在50%～10%，随着压缩程度的增加，椎间盘产生剧烈变形，并有变形溢出现象，这是由于椎间盘髓核等类似近似不可压缩流体，力学强度很低，而纤维环也属于高泊松比物质，在椎骨*力暴**压缩强度条件下，将会产生大变形，本研究针对此种现象，对椎间盘进行高阶六面体单元划分，该类型单元能够支持材料大变形计算，较好地模拟了椎间盘的变形状态[16]。
本研究首次以T₁₂椎体前缘不同压缩状态三维有限元模型分析邻近椎体终板的应力情况；以健康成人脊柱胸腰段三维有限元模型为基础，排除骨质疏松、椎体骨水泥填充、脊柱退变等干扰性因素，单纯研究T₁₂椎体前缘压缩程度与邻近椎体终板应力的相关性。轴向压缩状态下终板应力值最大，其次是前屈压缩状态及后伸压缩状态，分析原因为达到相同的压缩程度，轴向压缩条件下，整个椎间盘均受到压迫，产生的压应力最大，而前屈和后伸载荷条件下，椎间盘只有局部受压迫，其压应力小于轴向压缩状态。
综合研究内容，本文所研究的胸腰椎T₁₂～L1节段三维有限元模型包含椎间盘纤维环、髓核，韧带，小关节等特征，具有精确的几何外形，能够准确模拟实际临床力学环境，力学研究结果表明，通过有限元算法能够为临床提供良好的术前指导。
