5岁儿童上颈椎三维有限元模型的建立及对照分析.docxVIP

? ? 5岁儿童上颈椎三维有限元模型的建立及对照分析 ? ? 刘振杰，解英静，侯玲玲 儿童上颈椎尚未发育成熟，与成人相比，在解剖结构及力学稳定性方面存在很大差异，在遭遇外伤、上呼吸道感染时，容易诱发寰枢关节半脱位。近些年来，为评价各种病理状态下儿童颈椎、颅脑的损伤情况，学者们建立了不同年龄段的颈椎有限元模型，如吕文乐 等[1]建立了6岁儿童的全颈椎有限元模型；曹立波 等[2]对3岁儿童C4～5颈椎有限元模型进行了拉伸、弯曲试验验证。本研究通过建立5岁儿童上颈椎及颅底的三维有限元模型，并与成人模型进行对照研究，分析两者在旋转、前屈、后伸、侧屈运动时上颈椎活动范围及稳定性的差异，为临床正确认识该年龄段儿童上颈椎病变提供参考。 1 材料与方法 1.1 正常上颈椎几何模型的构建排除标准：3个月内有上呼吸道感染病史、颈椎及颅脑发育畸形、颈部外伤及手术史。分别选取1名5岁儿童及1名45岁成年男性为研究对象，采用螺旋CT连续扫描颅骨至C5椎体，层厚0.625 mm。将CT图像以DICOM格式导入Mimics 14.0软件中，根据灰度值不同，分割出图像中骨质部分，然后调用其中的Edit Masks功能对颅骨、C1～4椎体横断面图像进行修补和擦除，去除之间交叉像素，以方便建立各部分骨骼的3D模型，通过Smoothing功能光顺3D模型，去除表面的骨突及非特征性孔隙，以点云格式导入逆向工程软件Geomagic Studio 10.0进行封装、区域划分，拟合出各个骨骼的NURBS曲面模型。 1.2 正常上颈椎有限元模型的构建NURBS曲面模型导入有限元前处理软件Hypermesh 13.0中，通过2D模块的automesh功能将所有模型划分为2D单元，选取椎体上下终板、小关节突等部位的单元进行复制、偏移等操作，建立椎间盘、关节软骨等结构的3D模型，并根据解剖参数[3-4]建立寰枕前膜、寰枕后膜、翼状韧带、横韧带、齿突尖韧带及项韧带等其他相关韧带，见表1。 表1 各解剖结构的单元类型、力学参数[3-4] 在原始CT图像基础上建立的模型，其冠状面、矢状面、横断面并不与全局标准坐标系定位平面平行，从而影响载荷添加的准确性。为此，我们通过调整模型局部坐标系的定位，将其与全局标准坐标系一致，其中坐标原点定位在C4椎体前下缘中点，X-Z为冠状面、Y-Z为矢状面、X-Y平面为横断面。将边界条件设在C4椎体下表面，约束各节点的所有自由度。根据枕骨内板曲面弧度计算颅脑的几何中心作为加载质点，利用MPC-tie约束质点与颅底内壁各节点的连接关系，在颅骨几何中心点上分别添加20 N及40 N的垂直载荷模拟颅脑重量，围绕坐标系X、Y、Z轴添加1.5 N·m的扭矩载荷模拟头部前屈、后伸、旋转及侧屈运动。将上述完整的有限元综合模型导入Abaqus 6.14软件中进行计算及后处理。见图1。 图1 上颈椎CT图像及有限元模型图 A.儿童上颈椎CT横断面；B.儿童上颈椎CT三维重建；C.儿童上颈椎有限元模型；D.成人上颈椎CT横断面；E.成人上颈椎CT三维重建；F.成人上颈椎有限元模型 2 结果 2.1 模型有效性验证本研究在样本骨骼点云数据的基础上应用数字建模，所建有限元模型符合人体解剖学特征[5]。儿童模型包括17 440个节点和93 410个单元，成人模型包括29 700个节点和144 828个单元。加载工况为在颅骨几何中心添加1.5 N·m力矩模拟颅骨在3个坐标轴线方向的旋转运动，结果显示成人模型中C0～1、C1～2的活动范围与寨旭 等[6]研究的在体活动度一致，从而验证了模型的有效性。 2.2 上颈椎各节段活动范围分析见表2。上颈椎(C0～1、C1～2、C2～3)的旋转、前屈、后伸活动范围儿童模型显著大于成人模型，差异均有统计学意义(P0.05)；C0～1侧屈活动范围儿童模型大于成人模型，差异有统计学意义(P0.01)。 2.3 齿状突及关节突部位的平均Mises应力分析见表3、图2。成人和儿童模型齿状突及关节突部位的平均Mises应力大小依次为后伸前屈旋转侧屈。各运动方向的平均Mises应力儿童模型C0～1、C1～2关节突 [C0～1为1.53～4.22(3.08±0.94) MPa，C1～2为1.18～4.16(3.11±0.76) MPa]显著大于C2～3关节突[0.95～3.82(2.37±1.15) MPa]，差异均有统计学意义(P0.05)。在各运动方向上，C0～1、C1～2关节突的平均Mises应力儿童模型均显著大于成人模型，C2～3关节突成人模型均显著大于儿童模型，差异均有统计学意义(P0.05)。 表2 上颈椎有限元模型各节段活动范围比较 表3 上颈椎有限元模型齿状突、关节突的平均Mises应力比较 图2 两组模型前屈及后伸运动中Mises应力分