作者:李志远，申勇张英泽，王林峰，刘艳兵，武刚，张宇宸，任虎
[目的]建立人工颈椎间盘置换术前、术后(C4~7)的三维有限元模型，分析颈椎间盘置换术的特点，并探讨术后对邻近节段生物力学的影响。方法]选取1名行人工颈椎间盘置换的男性病人的术前、术后CT片，范围为C4~7(包括椎体、邻近椎体和邻近椎间盘)分别建立术前、术后的三维脊柱功能单元的有限元模型。对模型加载2N·m的力矩和10N·m的转矩，分析其在垂直载荷、前屈、后伸、侧屈及旋转的情况下，手术前后的椎体的变化以及椎体皮质骨、松质骨、椎间盘、终板的应力改变，以及活动度的变化。[结果]手术后轴向负荷力线后移，人工间盘邻近椎体上下终板所受应力，以及邻近椎间盘所受应力有所增加，但所受压力、张力不明显。节段活动度的变化减小不明显。「结论]人工颈椎间盘置换术有效的改善了脊柱功能单位的生物力学性能，但需谨慎选择人工椎间盘的大小，并且最大限度保留终板骨性组织，以避免出现术后假体下沉、移位、脱出、节段后凸等不良影响。
人工假体;颈椎间盘;生物力学;有限元分析
随着颈椎的发病率逐年增加，前路减压植骨融合技术广泛应用于颈
椎病患者的手术治疗。然而，临床研究显示融合会导致颈椎活动度的减少，
加速临近节段的退行性病变。为解决异常应力作用于融合后邻近节段的问题，避免前路减压植骨融合术的诸多并发症，作为“非融合固定技术”的代表-颈椎间盘置换手术在临床上得到广泛应用。然而关于人工颈椎间盘置换术对邻近节段影响的报道较少，因此，作者建立颈椎间盘置换手术术前、术后两个脊柱功能单元的三维有限元模型。为了能充分反映出人体的生物学特性，区别于以往基于尸体建模的模式，采用基于患者的实体建模，加强对椎间盘、皮质骨、松质骨、韧带及颈后纵韧带等后部结构的模拟，使其更接近于真实结构[1]，观察模型在正常范围内的垂直压缩、前屈后伸、侧弯情况下，手术前后活动度，以及邻近椎体、邻近椎间盘和后部结构的生物力学变化，研究其对邻近节段的生物力学影响。
1资料与方法
1.1CT扫描
选取1名行“颈椎人工间盘置入术”后6个月的志愿者，男性，48岁，身高cm，体重75kg，未发现骨性异常与脊柱畸形，采用CT自上而下行螺旋扫描，扫描范围包括C4~7全部骨性结构、韧带、椎间盘及人工椎间盘假体。
1.2建模方法
运用Photoshop语言数字图像处理程序对C4~7骨性结构和椎间盘区域进行边缘检测，提取边界坐标，并整理成点云文件。将点云文件输入逆向工程软件Geomagic，生成NURBS曲面后，将其实体模型通过IGES文件导入到有限元分析软件Ansys90中去，对椎体皮质骨、松质骨、椎间盘用三维十节点四面体结构实体单元进行网格划分，上下终板和纤维环被简化成一体，中间髓核部分被区分出来。根据CT值选取钛合金盖板轮廓，同骨组织类似在逆向工程思路下构建实体模型。高分子结构髓核部分在CT图像中显示并不清晰，在操作
中按确髓轴线位置，并根据人工椎间盘结构确定半径尺寸，然后建立在轴向方向略大于真实尺寸的初始髓核实体模型。通过布尔操作，用上下2个钛合金盖板实体模型减去初始髓核实体模型中多余的部分，即得到完整的人工椎间盘实体模型。在这个模型中，包括盖板、钛铢和前部凸缘在内的金属结构被统一建立为上下两个钛结构模型，在计算应用中这些简化对力学响应基本不构成影响。壳和髓核之间采用无摩擦的接触分析，这与Rousseau等[1]的假设不同，但更接近真实情况。韧带结构使用仅受拉的二节点杆单元模拟，包括前纵韧带、后纵韧带、*韧带、关节囊韧带、棘上韧带和棘间韧带。所有结构材料参数、韧带截面积都按文献报道选取，数值及单元划分情况如表1所示。
13计算与分析
对模型C7下表面全自由度约束，在C4椎体上表面分别加载均为2N·m前屈、后伸、侧弯和轴向旋转的弯矩。主要计算分析在各种状态下脊柱功能单元(functionalspinalunitFSU)生物力学变化。
对于一个ESU.下椎体是相对于其自身的下一椎体旋转或绕自身的旋转轴旋转的，在分析之初，必须把该旋转连带造成的上一椎体的运动量减去。在下一椎体相对固定的条件下，在上一相对活动的椎体上表面选择
两个节点，加载前后节点的三维坐标可以方便取得。对任意节点加载前后
坐标的连线可以在三维空间内做其垂直平分面。由于有限元分析中椎体的
弹性模量远大于其他组织的弹性模量，椎体的应变相对较小，所以可以近
似认为任意两个节点所得到的HAM都是一致的。但不恰当地选择两个节点
也可能造成较大的计算误差，这在其他文献中都没有提到。本文选择6对节点，将计算得到的HAM参数均值作为结果，减小了计算误差。模型结果包括:所完成的C4~7有限元模型共计包含个shell单元，个Solid单元79个缆绳单元。
2结果
21有限元模型结果表1模型构建及各材料系数组织单元类型单元数弹性模量术后脊柱功能单元模型如图1图2。
22模型分析计算结果如表2，图3所示。
2.3结果
图1颈椎术后模型(C4~7)侧位
图2颈椎术后模型(C4~7)正位(1)在垂直荷载压力下，C5下终板和6上终板表面应力与术前相比，应力分布有所变化，局部应力增加，但低于相关实验研究所报道的终板受损变形的最低值。
(2)与术前相比，C5、6椎间盘和C6、7椎间盘所受应力变化不大，但应力分布比术前集中。
(3)Bryan植入后该节段的活动度较术前略有减少，但在正常范围之内。
(4)本实验研究结果与王岩等[2]的临床研究结果相近，并且运动趋势与Goffin[3]相同。
表2垂直载荷下，术前和术后各组织vonMises应力比较(MPa)C5
下终板C5椎体C6上终板C45椎间盘C67椎间盘人工间盘术前5..43...97-术后.97..85图3模型脊柱功能单元在不同状态下的生物力学情况(MPa)
图3a前屈图3b后伸图3c侧弯图3d轴向旋转
3讨论
近年来,诸多临床研究报道[4]，相比颈前路融合术后而言，人工颈椎间盘置换术有着保留颈椎运动功能，术后对邻近节段影响较小等优势。而本实验的研究结果也再次印证颈椎间盘置换术后，颈椎运动趋势近似于术前，并且节段活动度也接近正常。这一结果与Galbusera等[4]学者的研究一致。因此，再次说明了人工颈椎间盘置换术在恢复颈椎稳定性的同时，也较好地保留了颈椎的运动功能，这一点得到了DuggalN等[5]影像学分析结果的支持，并且与HaSK[6]所报道的生物力学理论吻合。然而，值得
