标准步态下骨盆三维有限元模型构建及其生物力学意义

标准步态下骨盆三维有限元模型构建及其生物力学意义 标准步态下骨盆三维有限元模型构建及其生物力学意义 Threedimensional finite element pelvic modeling during standard gait and its biomechanical significances 【Abstract】 AIM: To construct a threedimensional finite element pelvic model hich includes both 3D pelvic muscle and acetabular contact forces during standard gait. METHODS: A normal pelvis from an adult male subject as scanned by CT and the images of every crosssection ere obtained. A pelvic model including both threedimensional pelvic muscle and acetabular contact forces during standard gait as constructed ith PATRAN 2005R2 softare. RESULTS: The constructed 3D finite element pelvic model clearly reflected the real pelvic anatomy and biomechanical behavior, especially the muscle force and hip contact force. The model as divided into 113 028 nodes and 137 524 units. CONCLUSION: The construction of 3D finite element pelvic model provides basic data hich are critical for accurately modeling either normal loads or stresses and strains, or the effects of abnormal conditions. 【Keyords】 pelvis; threedimensional finite element; acetabular contact force; pelvic muscle force 【摘要】 目的: 构建包含骨盆肌肉及头臼作用力的标准步态下骨 盆三维有限元模型. 方法: 选择标准成年男性志愿者行骨盆CT扫描 成像得到骨盆每层横截面图像, 运用Mimics软件行三维重建， 利用 有限元分析软件PATRAN 2005R2构建髋臼三维有限元模型,并将标准 步态中肌肉收缩力及头臼作用力也设计到该模型中. 结果: 所构建 髋臼模型共划分为113 028个结点、137 524个单元, 模型模拟了骨 盆皮质骨、松质骨、关节软骨以及股骨头结构的特性. 所建模型 结构完整，空间结构测量准确度高，单元划分精细，重点突出，尤其 肌肉及头臼作用力加载模拟形象、逼真，客观反映髋臼真实解剖形态 结构和生物力学特点. 结论: 构建的髋臼三维有限元模型为正常髋 臼或异常髋臼力学研究提供可循模型. 【关键词】 骨盆;三维有限元;髋臼作用力;骨盆肌肉收缩力 0引言 对于髋关节的生物力学的研究，有必要了解跨髋关节肌肉的收缩力的大小方向及部位，完整地了解髋臼接触力与骨盆肌肉收缩力力学机制，有助于进一步模拟髋关节力学行为、解剖学与手术后改变的影响，虽然有文献,1-3,讨论过骨盆髋臼生物力学行为，然而这些先前的模型或没包括髋臼接触力或不包括肌肉收缩力. 因此其结果已受到大多学者的质疑. 作为实验研究骨盆与髋臼负载的需要，本实验试图报导一个既包含了肌肉又包含了髋臼接触力的三维有限元模型，这些数据将对于准确模拟髋关节与骨盆的正常负载或异常负载将尤为重要. 1材料和方法 1.1材料? 成年男性志愿者骨盆(男性，40岁);? 西门子SOMATOM Volume Zoom CT 机;? DELL工作站;? Mimics软件;? Patran 2005 r2 软件;? Windos XP操作系统. 1.2方法 1.2.1骨盆三维有限元模型构建具体可以分为3个步骤? 步骤一: 对一成年男性志愿者(40岁)的骨盆进行CT断层成像;在CT成像过程中要志愿者在骨盆纵轴方向保持不动，每隔1 mm层厚扫描一次. 所得图像以DIMCOM格式存入CT机，刻录光盘，从而得到示骨盆每层横截面的图像. ? 步骤二: 将上述以DIMCOM格式存储的图像导入Mimics软件，设定门槛下界灰度值为270，上界不限. 通过Region groing 选择感兴趣区域，进行骨盆三维重建. ? 有限元分析与计算: 三维重建图像以out格式导出，并导入Patran 2005 r2软件，并重新进行网格划分，建立皮质骨的外表面与内表面，其内外表面都相互连接，分别建立骨盆的皮质骨及松质体模. 与骨盆骨相互作用的股骨头也被建模，此外关节软骨也被建模，以保证一个光滑且真实的髋关节作用力的加载,4,，除了髋关节头臼作用力，22块附着在骨盆上的肌肉也被设计到该模型中，肌肉的生理附着点被绘制到有限元网格中，其相关参数皆来自于Pedersen1997年的研究结果,5,. 图1髋臼皮质骨与松质骨有限元模型 略 表1步态循环中支撑相时负载受力参数 略 1.2.2设置单元属性根 据骨盆的解剖结构特点，骨性结构模拟材料为皮质骨、松质骨、软骨下骨及关节软骨，材料参数的选择参考文献,6-7,(表3). 1.2.3加载情况固定骶髂关节处的节点模拟骶骨的支持，只对左半侧骨盆进行加载，在耻骨联合处模拟弹性边界条件. 根据步态循环周期中左支撑相的第4个子步相的髋关节的头臼之间的作用力,从股骨头的中心点处进行加载. 此外肌肉收缩力也被考虑，收缩力的方向是根据远近端的附着点获得,8,，并且步态中骨盆肌肉与股骨的相对位置关系也被考虑进去. 本实验中设定股骨头内收15?，同时在前后位上骨盆与股骨头的角度是可以变化的. 肌肉的收缩力以分布载荷运用到模型的表面，肌肉的收缩力的方向是由其远近端的生理附着点的连线，并同时兼顾髋关节的屈伸角度对远端生理附着点的坐标影响. 假设条件:本实验所涉及的生物材料均假定为均质、连续和各向同性. 受力时模型各截面不产生相互滑动，各单元有足够的稳定性;材料受力变形为小变形. 表25步态载荷中肌肉收缩力的大小 表3骨盆三维有限元模型单元属性表 略 2结果 2.1模型建立本研究中，应用CT图像专用的DICOM医学数字图像通讯标准，CT专用的Mimics图像识别分割建模软件和PATRAN有限元专用软件，将薄层CT横断面图像进行分割识别，建立了适合生物力学研究的髋臼三维有限元模型，肌肉的生理附着点被绘制到有限元网格中，并由肌肉的远近端的附着点获得肌肉收缩力的方向，并兼顾髋关节的屈伸角度对远端生理附着点的坐标影响. 髋臼三维有限元模型共划分为113028个单元，137524个节点. 股骨头采用四节点tet10单元，共划分为37090个单元50392个节点. 松质骨模型在髂骨中央有椭圆形孔，及髂前下棘处有缺失，既该两处为皮质骨充满. 旋转观察模型形态与骨盆解剖形态具有满意的相似性，可任意切割和调整几何及材料参数以模拟不同临床与实验状态. 表4模型的有限元如下表 略 2.2模型的受力分析模拟步态4髋关节的头臼之间的作用力及附着在骨盆上的22块肌肉的收缩力对模型进行加载. 得到该情况下髋臼 的力场分布. 结果显示:主应力集中在髋臼软骨的臼顶区，传递至软骨下骨后，沿弓状线及臼顶的外上皮质至骶髂关节. 3讨论 3.1骨盆有限元模型的应用1995年Dalstra等首次详细报告了骨盆三维有限元模型. 2000年Garcia等应用三维有限元对内固定后的骨盆进行应力分析，研究骨盆骨折后不同固定对其稳定性的影响;2006年Ye运用有限元分析影响髋关节假体置换稳定性的因素,9,. 国内张春才等,10,报道了应用三维有限元法对骨盆骨折损伤机制的分析. 在骨折内固定的应力分布、人工关节的设计优化研究，骨折危险性预测、骨折愈合、预测骨的塑型以及骨质疏松症研究方面有限元模型已是不可缺少的方法,7,. 图3单腿站立位时皮质骨,松质骨和髋臼软骨受力分布 略 本研究应用16排螺旋扫描 CT扫描，并用MIMICS软件直接读取DICOMR的CT文件，选取兴趣区域行3D计算，获得髋臼的三维模型，并导入MSC公司PATRAN软件，建立了髋臼的三维有限元模型，并对重点部位进行细划分. 与标本的相似程度高，形态结构准确完整. 操作过程使图像失真减少到最小. 以往,4-5,建立骨盆三维有限元模型采用的 CT胶片扫描或照相技术, 在对图像反复采样过程中, 图像的细节和信息在采样过程中丢失较多, 而且工作繁重而导致建模差异较大, 出现建模质量的问题,影响到计算结果的准确性. 本实验三维建模过程为全数据化处理，避免了数据收集过程中关键信息的丢失，最大限度上保证了建模的准确性和精确性. 3.2模型的力学相似性本研究结合了结合头臼之间的作用力与肌肉的收缩力，不当在外形上最大限度上保证了建模的准确性和精确性. 并且在力学加载上更与实体相似程度高，更加准确完整. 本研究头臼作用力由股骨头中心点处进行加载，同时参予髋关节运动有关的22块肌肉也被考虑，肌肉收缩力的方向由远近端附着点决定，并同时兼顾了髋关节屈伸角度对肌肉远端附着点的影响;肌肉收缩力大小与方向及头臼作用力大小与方向参考相关文献,5-6,,并转换为本研究的坐标系，更加真实地模拟实际步行情况. 3.3髋臼受力的分析生理状况下，骨骼的结构和功能在很大程度上 依赖于其所处的力学环境. 在髋臼骨折治疗中, 施加的力学环境, 对治疗效率、准确的愈合是至关重要的. 本研究尝试对髋