论文部分内容阅读
目的:踝关节是人体承重关节中最末端的关节,由于其特定的解剖结构和复杂的运动机制,造成了该关节是运动损伤的常见部位。本研究通过CT数据的采集,对人体踝部主要结构进行三维重建,从而建立有限元数字模型,在不同的载荷下,对踝关节各组成部位的受力进行有限元分析。以此方法来探讨踝关节各组成骨的应力分布规律及发生形变情况。并通过与临床统计数据和踝足系统其他有限元模型的生物力学研究对比,来验证本模型的有效性,从生物力学的角度分析踝关节损伤的原因,为临床上研究踝关节损伤的机制提供理论支持。
方法:本实验首先选取一名中国成年女子(志愿者),34岁。身高164cm,体重60kg,踝部无畸形,无明显外伤史。对其右踝足部行螺旋CT扫描。范围自胫腓骨中下段1/3至足跟部,总长度16.625cm。采用仪器为美国GE公司的lightspeed16排螺旋CT。扫描参数为:120kv,240mA,螺旋层厚1.25mm,曝光时间12.053s。共得到133幅分辨率为512×512像素的二维断层图像。将CT扫描数据以DICOM格式存储入可读写光盘。为简化模型,本实验选取包含胫骨下段、腓骨下段、距骨在内的共82幅横断扫描图像,建立踝关节的数字模型。将DICOM文件导入MIMICS10.0软件,对原始图像进行蒙皮、分割、光滑等处理,生成踝部的三维模型。并以点云的形式输入逆向工程软件Geomagic9.0,再生成实体模型。将得到的实体模型经过装配,导入到有限元分析软件Ansys10.0中去,格式为IGES。采用十节点的四面体结构实体单元SOLID92模拟骨组织,三维仅受压杆单元LINK8模拟软骨,最终生成踝关节的三维有限元模型。
对模型施加载荷,进行求解试算。对胫腓骨进行全方位约束,在距骨下端加载纵向向上压力300N(相当于静止站立时,单踝关节所承受的压力),运算求解得出正常生理状态下,静止中立位相时胫腓骨中下段和距骨的形变图、Von Mises应力云图以及位移图。然后在相同的约束条件下,增加载荷至600N、900N、1200N、1500N、1800N,(资料表明,踝关节最大能承受人体5倍体重的压力)模拟坠落伤时,踝关节中立位相时,各组成骨所承受的应力和位移情况。通过观察所得踝关节各部的形变图、应力云图和位移图变化,判断哪些部位最容易发生骨折。并将其与临床统计数据和已有的生物力学理论对比,验证模型的有效性。
结果:本实验成功建立了包括胫腓骨中下1/3段和距骨的踝关节数字有限元模型。采用10节点四面体solid92模拟骨组织,仅受压杆单元link8模拟关节软骨;网格划分尺寸为7mm。建立了共包含11603个单元,17882个节点的踝关节的数字模型。得到了正常生理状态下,中立位相时静止站立及高空坠落情况下,踝关节各组成骨所发生的形变图、应力云图和位移图,并以形象、直观的彩图形式展示结果。踝关节在正常生理情况下,中立位静止站立时(只承受自身重力时),最大应力区集中在胫骨中下段前内侧皮质、胫骨远端关节面外侧部和髁间线前部、距骨滑车外侧部、内踝;最大位移在距骨头;逐渐增加载荷,至相当正常人体6倍体重时,踝关节各部所受应力均逐渐增加;其中胫骨中下段最大应力集中区由前侧皮质移动至内侧皮质。这与汤欣的四柱分型理论相符合。
结论:
1、利用CT扫描原始数据建立了踝关节三维可视化模型。
2、建立了包括胫腓骨中下1/3及距骨的踝关节三维有限元模型。
3、分别模拟中立位相静止站立以及高空坠落的工况,对踝关节各组成骨进行有限元分析,得到了踝关节各组成骨的形变图、Von Mises应力云图以及位移图。从图上可以看出,中立位相时,踝关节各组成骨所受应力主要分布在胫骨远端关节面外侧、髁间线前部、距骨滑车的外侧、内踝和胫骨中下段前缘。这完全符合汤欣的四柱理论。临床统计数据也表明,踝关节损伤中以胫骨远端关节面外侧、内踝最多见。与本实验结果相符。从而间接证实了本有限元模型的有效性,为踝关节损伤病因学研究提供一定的理论支持。也为后续的生物力学相关研究提供了一个二次开发的数字平台。