扁平足站立中期足底压力与踝关节分析.docxVIP

?

?

扁平足站立中期足底压力与踝关节分析

?

?

王光正石光林吴声翔高聪梁潇帝商搏世

摘?要：为给扁平足生物力学研究以及矫形鞋垫设计提供数字化平台，构建并验证了扁平足三维有限元模型.首先基于扁平足患者的右脚CT图像数据，依次利用Mimics、Geomagic、UG和ABAQUS进行处理并生成扁平足足部有限元模型;然后通过模拟双足站立中期得到的足底压力数据与设备采集的数据进行对比，验证模型的有效性.仿真结果显示：其压力峰值与正常足一样，主要集中在跖骨头区和足跟区，依次是0.20MPa和0.24MPa，而设备测量值依次是0.21MPa和0.26MPa，误差分别为4.8%和7.7%;在足跟区，扁平足的压力峰值与压力分布主要集中在足跟内侧，与正常足不同;踝关节应力最大值是在距骨颈部位置，形变最大值是在腓骨处.

关键词：扁平足;有限元分析;足底压力;踝关节;生物力学

：TH122∶R318.01?DOI：10.16375/45-1395/t.2021.02.007

0引言

近年来，有限元分析法凭借低成本和高效率的特点在医学领域得到迅速发展，尤其是在骨科方面的研究.由于传统研究主要依靠尸体来进行试验，要通过相关的实验设备来获取数据，所以传统方法进行试验比较困难，数据并不理想.更多学者开始采用有限元方法进行生物力学研究，并证明了这种方法可行且有效，这也促使学者对人体各部位进行生物力学有限元分析的相关研究.金乾坤等[1]建立正常足有限元模型并分析足底压力以及跖骨骨骼应力;张明等[2]通过模拟足部平衡站立时的仿真数据得出：跟腱力为总体重的75%的时候足底压力分布与测量结果接近;徐鉴等[3]通过有限元分析研究不同类型距下关节制动器治疗Ⅱa期成人获得性平足效果的差异;何晓宇等[4]建立正常骨骼和常见疾病有限元模型并进行生物力学分析，比较它们之间的骨骼应力分布情况.但是关于扁平足的有限元分析比较有限，尤其是在足底压力方面.

本文拟通过构建扁平足有限元模型并进行足底压力分布和踝关节研究，为扁平足足部生物力学研究以及矫形鞋垫设计提供参考依据.

1方法

1.1??足部数据的提取以及三维模型的构建

基于一名男性扁平足患者（19岁，体重72kg，足部有内侧足弓塌陷，拇外翻，后脚跟有一定程度的外翻症状）的右下脚CT断层图像，采用医学建模软件Mimics21.0，根据解剖结构对DICOM格式的CT图像分别对足部软组织和骨骼进行选择性分割，利用区域增长填补小的空隙和修整轮廓线，使整个模型外轮廓线光滑而连续，通过计算模型选项分别建立足部软组织和骨骼的三维模型.构建好的模型信息以STL格式进行保存并将STL格式文件导入逆向工程软件Geomagic13.0进行曲面优化和重构处理.去除模型表面非特征性不光滑的地方，光滑松弛表面，然后用曲面去拟合模型表面，最后生成表面光滑且连续的几何模型.将构建好的几何模型以Step格式存储并导入三维建模软件UG10进行组合装配，然后参考CT图像足底筋膜和足部韧带的起止点及走向信息，建立三維曲线模拟足部韧带，构建相对完整的足部模型（见图1）.

1.2??网格划分

将完整的足部模型导入ABAQUS6.13，分别对骨骼、韧带、足底筋膜以及地面支撑板进行网格划分，划分后单元总数为164233，网格节点总数为245785.其中足底筋膜、足部韧带均采用Truss单元，骨骼采用四面体C3D8R单元，地面支撑板采用六面体C3D8R单元，软组织采用四面体C3D4单元.

1.3??材料属性

整个足踝模型由趾骨、跖骨、楔骨、骰骨、足舟骨、距骨、跟骨、胫骨以及腓骨这几部分构成.除了软组织之外，骨骼、足底筋膜和韧带这几部分均视为均匀、各向同性、线性弹性，其弹性模量、泊松比、横截面积以及单位类型参数见表1[5-8].外围软组织采用超弹性材料，具体的超弹性材料系数见表2[2].

1.4??边界条件及接触设定

为了真实模拟患者双足站立中期的情况，外围软组织的上表面、胫骨和腓骨上端面进行完全固定约束.由于扁平足患者双足站立时，单足承受人体一半的重力，所以在地面支撑板下表面中心施加方向向上的360N集中力（见图2）.

在整个分析过程中关节间、软组织与骨骼接触以及足底与支撑板之间出现相互接触的情况.关节间设定为无摩擦接触，外围软组织与支撑板设定摩擦接触，其摩擦系数为0.6.另外足部骨骼与外围软组织之间设定为绑定接触.

1.5??足底压力测试

通过德国Go-tec公司的足底压力测量系统，对扁平足患者进行足底压力数据的采集.采用测量系统在静态下采集到的患者足底压力数据，观察整体压力分布以及压力峰值情况.

2结果与分析

2.1??足底压力验证

从有限元仿真的足底压力云图（图3）可知（从俯视图方向观察），足跟区域的压力峰值是最大