正常步态下距骨三维有限元模型的建立.docx

《正常步态下距骨三维有限元模型的建立.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、正常步态下距骨三维有限元模型的建立摘要：目的：建立一个具有高度几何相似性的足踝部三维有限元模型，并用此模型静态地分析人体不同步态相时距骨的生物力学特性，量化距骨的应力-应变状况。方法：基于志愿者足的三维CT扫描数据，利用Mimics13.0、Geomagic10.0软件对足踝相关组织进行几何重建，导入Hypermesh10.0软件中进行网格划分，赋予材料属性，最后导入相）时的受力状况，进行有限元分析。结果：建立距骨及周围结构的三维有限元模型，共21865个节点、73440个单元，具有较好的几何相似性。正常步态中从落地相到离地相中等效应力峰值在距

2、骨滑车分别为3.0MPa、4.3MPa、4.8MPa；在距骨颈分别为1.3MPa、1.9MPa、2.8MPa；在距舟关节分别为2.8MPa、3.0MPa、3.4MPa；在距下关节分别为2.2MPa、1.8MPa、1.5MPa。结论：本研究中创建的三维有限元模型，经验证是一个正确可靠的模型，可以帮助临床医生和其他研究人员更好的理解正常步态下距骨的许多生物力学特性。关键词：距骨；生物力学；有限元分析；步态引言：距骨具有独特的解剖结构和功能，无肌肉附着，大部分骨质被关节软骨面包围，是人体重力传至足部的枢纽，故有“骨半月板”之称［1］。人体在行走时，足

3、踝的生物力学变化比较复杂，因此研究正常步态下距骨生物力学变化较为困难。全面了解距骨在正常步态下的应力分布可以为我们研究距骨生理学和病理学提供十分有用的信息，理解距骨在正常步态下力的传递对足的损伤及治疗有重要意义。活体组织研究能得到较为精确数据，但因伦理学的原因受到限制。随着计算机技术和有限元理论的不断发展，人们开始大量使用数值模型和有限元法分析复杂的结构。与以往的生物力学实验相比，有限元方法可以建立高度几何相似及物理相似的有限元模型，既可以反映区域性力学特性，有可以反映整体信息；既可以进行精确的数字分析，有可以进行形象的、直观的定性研究。有限元

4、分析方法的最大特点是对生物体的无损伤，可以模拟活体组织内部的生物力学行为，完成其他研究方法所不能实现的加载方式及约束条件，得到客观实体实验所难以得到的研究结果。可以通过改变载荷加载方式、改变材料特性等方法进行个体化受力分析，研究省时快捷，费用低廉，应用面广。本研究基于CT数据建立距骨及周围结构的三维有限元模型，分析不同步态下距骨结构的应力变化，以期为距骨生物力学研究提供新的研究手段。1材料与方法1.1数据收集：男性志愿者1名，25岁，身高170cm，体质量60kg。先行X线检查，排除足部肿瘤、畸形等其他病变，进行右足踝部64排螺旋CT平扫，层厚

5、0.45mm，数据以DICOM格式输出保存。1.2有限元模型建立CT扫描的数据导入三维重建软件Mimics13.0，通过阈值分开骨组织及软组织，建立距骨及周围结构的几何模型，输出为STL文件。然后导入逆向工程软件GeomagicStudio10.0中，对模型进行除噪点、平滑，根据各关节面的几何形状，在各骨面上划分软骨边界，最后拟合曲面，输出为Iges格式。再导入有限元前处理软件Hyperwork10.0中，按照上述生成的软骨边界，生成软骨；依据韧带解剖数据，通过五条纤维束三维排列连接韧带附着点建立韧带模型，最后对实体模型进行网格划分、设置材料属

6、性等处理。最终导入ABAQUS6.9软件进行计算分析等后处理。有限元模型见图1。1.3材料参数：骨性结构和软骨均模拟为各向同性的线弹性材料，韧带使用非线性单轴连接单元来模拟其只受拉、不受压的特性。骨、软骨、韧带的材料特性由实验文献确定［2-6］，列于表1、表2。1.4边界条件及载荷：对步态过程中位相的划分方式很多，但通常将其分为落地相、中立相、离地相3个位相。正常步态踝关节运动方式由距骨体滑车关节面确定，侧面观，距骨体滑车形状近似圆锥一部分，腓侧曲率半径大于胫侧，踝关节在矢状面屈伸运动轴不是水平的，自内上向外下倾斜，并且也不是恒定不变的，但是在

7、正常步态中，踝关节屈伸运动的瞬时旋转中心位于距骨体内并且十分靠近，因此可以近似认为是一个点。本研究屈伸运动轴及旋转中心的确定根据以前学者的研究方法，确定距骨滑车胫、腓侧穹弧的圆心，连接两个圆心，即为旋转轴，两圆心的中点为旋转中心［5,7-8］。踝关节在步态运动中的受力状况研究采取将其作为相对静止加载的近似方法，加载在胫腓骨截面上。关节面之间建立接触对，摩擦系数为0.01［5］。不同位相受力状况参考文献数据［9］，见表3。2结果：建立距骨及周围结构三维有限元模型，其形态还原性好，重建效果较理想，能任意旋转，可获得详细、满意的三维信息。模型分析结果

8、与以往研究具有相似性［4-5,10］，认为模型有效，可以用来进一步分析计算。足踝部不同位相VonMises应力分布如图2、3。距骨不同区域最高等效应力