聚醚醚酮手指假体三维有限元模型的建立及生物力学分析

时运 覃祖敢 陈程 程丽佳

摘 要：建立食指三维有限元模型，并将聚醚醚酮材料属性赋值，分析在施加不同载荷下食指各节指骨的应力分布情况.健康男性志愿者1名，对其左手进行电子计算机断层扫描，将扫描数据导入Mimics 21.0软件中得到食指的点云数据，再导入UG 12.0软件中进行逆向处理建立食指模型，最后导入Workbench 2020R1软件中，沿食指指骨纵向分别施加10、20、30、40和50 N载荷，分析食指各节指骨的应力分布情况.结果显示，建立了高精度食指三维有限元模型，模型共计551 410个节点，392 047个四面体单元；随着施加载荷的增加，食指指骨的应力变化范围约0.37～

5.96 MPa，远节指骨的整体应力要明显高于其他指骨，近节指骨所受的应力最小.采用Mimics和Workbench软件可建立精确度较高的食指有限元模型，从而对食指生物力学角度进行定性定量的分析研究，为聚醚醚酮材料作为食指替换材料提供一定的理论依据.

关键词：聚醚醚酮；食指；有限元分析；三维重建；生物力学

中图分类号：R318.08;TQ317

文献标志码：A

0 引 言

临床上高能量损伤导致手部组织严重缺损的情况日益增多，尤其是大范围手部损伤将严重影响手的功能.根据损伤程度轻重选择不同的治疗方式，有的需要截肢，手指截肢患者不仅面临手功能的严重下降，有时甚至还面临着更严重的社会心理障碍^[1-5].因此，寻找适合的手指假体材料显得尤为重要.目前，大多数假体为硅橡胶假体，这种材料会导致关节不稳、假体破裂和硅橡胶滑膜炎等问题.钽和钛作为应用最广的骨植入金属材料，因其存在应力屏蔽与骨结合能力弱等缺点使其应用受限^[6-7].聚醚醚酮以其优异的耐磨性、生物相容性与化学稳定性等优点成为目前最具应用前景的人工骨基体复合材料，并可独立作为人工骨替换材料使用^[8-9].聚醚醚酮的弹性模量为（3.86±0.72）GPa，与人体骨的弹性模量相似，能够避免植入后产生的应力遮挡及松动现象^[10-14].医学三维影像重建及有限元技术成为了骨科领域各个部位生物力学研究的重要手段，但对手部关节方面的有限元分析研究较少.进一步了解手指关节生物力学特点，对于手指关节炎与关节僵硬的诊治具有重要意义.本研究采用了Mimics 软件及Workbench有限元分析软件建立食指关节有限元模型，模拟聚醚醚酮材料在不同载荷条件下食指关节的应力分布，为聚醚醚酮成为食指假体材料提供一定的理论依据.

1 方 法

1.1 手部图像数据采集

健康男性志愿者1名，24岁，体重80 kg，身高181 cm，手部无病史.志愿者知晓及同意试验方法，并将电子计算机断层扫描（CT）资料用于本研究使用.本研究通过成都大学附属医院伦理委员会审核.志愿者俯卧在64排CT（东软医疗系统股份有限公司）床上，左上肢伸于头侧并保持水平前伸.将手部CT图像以DICOM格式保存，CT参数为120～140 kV，

525 mA，层厚0.625 mm.

1.2 建立食指三维模型

将得到的DICOM格式图像导入Mimics 21.0（Materialise公司）软件中，设置阈值155～最大值，使用区域增长（region grow）提取骨骼，再用多层编辑工具（multiple slice edit）在各軸视图层面根据骨骼的形状轮廓提取食指，去除多余的部分.最后，将Mimics处理得到的点云数据导入工程设计软件UG 12.0（Siemens PLM Software公司 ）中进行逆向处理建立食指模型，如图1所示.食指关节三维重建的流程如图2所示.

2 食指关节三维重建的流程

1.3 建立食指有限元模型

将食指关节STP模型导入Workbench 2020R1（Ansys公司）软件中进行网格化（见图3）和材料属性赋值.经检索国内外文献研究^[15-18]，食指关节对应的聚醚醚酮材料弹性模量为4 GPa，泊松比为0.38.网格采用四面体单元，且网格大小为0.5 mm.

1.4 设置边界条件和施加载荷

将食指的近节指骨、中节指骨与远节指骨的指骨底作为固定边界，约束限制6个方向的自由度.为了模拟食指各指骨结构的受力情况，在每个指骨的前端沿着指骨径向分别施加10、20、30、40和50 N载荷^[13].

2 结 果

对食指的近节指骨、中节指骨与远节指骨分别施加10、20、30、40和50 N的不同载荷，得出食指各节指骨对应的冯·米塞斯（Von Mises）应力云图（选取50 N为例）及应力峰值，如图4和图5所示.

按照本研究材料聚醚醚酮的材料参数，在施加不同载荷过程中，食指应力峰值变化范围约为0.37～

5.96 MPa，其中，近节指骨的应力峰值变化范围在0.37～1.84 MPa，中节指骨的应力峰值变化范围在0.56～2.81 MPa，远节指骨的应力峰值变化范围在1.19～5.96 MPa.

随着施加载荷的增加，食指各节指骨所承受的应力呈现不同程度的增大.在不同载荷作用下，远节指骨的整体应力要明显高于其他指骨，其中，近节指骨所受的应力最小，且各节指骨的最大应力主要集中在指骨的中段.

3 讨 论

有限元分析法是一种用于分析复杂固体力学和结构力学系统的方法，能够对人体骨骼的生物力学特性进行研究，因此被广泛应用于骨科生物力学研究领域.金波等^[19]采用Mimics结合Ansys有限元分析软件对膝关节半月板有限元模型进行研究，结果显示，有限元分析结果与实体基本一致，因此，该方法可以作为临床治疗半月板损伤的有效手段.张帅等^[20]利用SolidWorks软件建立种植体在人类下颌骨的三维有限元模型，对种植体及下颌骨各部位的应力分布进行研究，结果显示，基台有缓冲装置可降低种植体及各部位所受应力.

手指的生物力學是由骨、关节、韧带、肌肉和肌腱的机械力学作用的复杂组合^[21].食指是仅次于手的最重要的手指，因为其不仅具有外展和内收的能力，还具有弯曲和伸展的能力，并且靠近拇指.多项研究指出，食指在日常生活中使用的重要性，失去食指会造成捏、握与旋转等力量下降约20%^[22-23].聚醚醚酮作为一种骨科植入物材料，其弹性模量与人骨相当，具有优异的生物力学性能、生物相容性和放射线透过性，植入人体后可有效降低应力屏蔽效应，且磁共振成像（MRI）不产生伪影^[24].随着3D打印技术趋于成熟，可结合CT和MRI等医学扫描技术个性化定制聚醚醚酮植入物.目前，3D打印聚醚醚酮逐渐被应用于临床，例如，口腔、人造骨关节、颌骨及颅骨等缺损修复^[25-26].聚醚醚酮材料正逐渐成为骨与软骨替代材料中的一类重要的生物材料，并应用于医学的各种领域.

本研究采用CT影像重建技术对食指整段进行三维重建，并在有限元中进行聚醚醚酮材料属性赋值，对食指有限元模型沿着指骨径向施加不同载荷，分析食指应力的变化范围和分布情况，从生物力学角度对食指的结构受力展开定性定量的分析研究.通过分析得出，聚醚醚酮食指在沿指骨径向分别施加不同载荷时，其应力峰值分布范围在0.37～5.96 MPa，聚醚醚酮食指内的最大应力是其屈服应力的4.0%，不会因受压力过大而导致变形.因此，聚醚醚酮食指种植体材料符合一定的机械力学和生物力学性能要求，为聚醚醚酮材料作为食指替换材料提供一定的理论依据.此外，关于食指的生物力学性能方面参考研究较少，因此，其实际生物力学性能还有待进一步研究.

参考文献：

[1]薛金伟.面向前臂截肢患者的虚拟现实康复方法研究[D].深圳：中国科学院大学（中国科学院深圳先进技术研究院），2020.

[2]Kong Y K，Seo M T，Kang H S，et al.Effect of finger interaction on individual finger：Index finger[J].Work，2012，41（1）：5696-5698.

[3]李金虎，杨新蕾，米增法.早期整形手术治疗手部深度烧伤患者的临床效果[J].中外医学研究，2022，20（1）：124-127.

[4]徐胜，刘志诚，涂鑫涛，等.基于手势识别的手部康复智能评估算法研究[J].自动化仪表，2021，42（12）：23-27.

[5]耿鹏鹏，许路，田洁.手部重度烧伤患者伤残接受度与应对方式的相关性分析[J].现代诊断与治疗，2021，32（20）：3289-3290.

[6]曾春，金柯含，李凯，等.聚醚醚酮与钛网在颅骨修补术中的应用分析[J].中华神经外科杂志，2021，37（12）：1219-1222.

[7]赵丙辉，李沅，李文娇.3D打印技术制备骨科植入材料聚醚醚酮的力学性能研究及氨基改性对材料成骨性能的影响[J].中国医疗器械信息，2021，27（13）：10-12.

[8]Oladapo B I，Zahedi S A，Ismail S O，et al.3D printing of PEEK and its composite to increase biointerfaces as a biomedical material-A review[J].Colloid Surface B，2021，203（3）：111726-1-111726-9.

[9]张勇，吴正敏，陈溯.聚醚醚酮支架对All-on-4远中倾斜种植体应力影响的三维有限元分析[J].北京口腔医学，2021，29（5）：309-312.

[10]吴水华，陈朝晖，范双石，等.个性化聚醚醚酮植入物在儿童颅骨修补术中的临床应用及分析[J].中华神经创伤外科电子杂志，2021，7（4）：252-255.

[11]Theivendran K，Arshad F，Hanif U K，et al.Carbon fibre reinforced PEEK versus traditional metallic implants for orthopaedic trauma surgery：A systematic review[J].J Clin Orthop Trauma，2021，23：101674-1-101674-10.

[12]Hosameldin A，Osman A，Hussein M，et al.Three dimensional custom-made PEEK cranioplasty[J].Surg Neurol Int，2021，12（30）：587-1-587-8.

[13]康永强，吴永伟，马运宏，等.正常掌指关节三维有限元模型的建立及生物力学分析[J].中华手外科杂志，2019，35（5）：364-367.

[14]康永强.2-5掌指关节屈伸活动中有效关节接触面积的三维可视化研究[D].苏州：苏州大学，2017.

[15]Shrivastava S P，Dable R，Raj A P N，et al.Comparison of mechanical properties of PEEK and PMMA：An in vitro study[J].J Contemp Dent Pract，2021，22（2）：179-183.

[16]Manzoor F，Golbang A，Jindal S，et al.3D printed PEEK/HA composites for bone tissue engineering applications：Effect of material formulation on mechanical performance and bioactive potential.J Mech Behav Biomed Mater[J].2021，121（15）：104601-1-104601-11.

[17]Baek I，Kwon O，Lim C M，et al.3D PEEK objects fabricated by fused filament fabrication （FFF）[J].Materials，2022，15（3）：898.

[18]Honigmann P，Sharma N，Schumacher R，et al.In-hospital 3D printed scaphoid prosthesis using medical-grade polyetheretherketone （PEEK） biomaterial[J].Biomed Res Int，2021，2021（12）：1301028-1-1301028-8.

[19]金波，胡云根，韓磊.半月板三维有限元模型建立及力学分析[J].中国骨伤，2020，33（8）：766-770.

[20]张帅，吕川，李进红，等.基台颈部缓冲层对种植体应力分布影响的三维有限元分析[J].华西口腔医学杂志，2020，38（5）：537-540.

[21]Duncan S F M，Saracevic C E，Kakinoki R.Biomechanics of the Hand[J].Hand Clin，2013，29（4）：483-492.

[22]李洋，张开治，杨寿君，等.足部和手部有限元分析模型的建立及应用进展[J].吉林医学，2007，28（7）：957-960.

[23]林柳兰，周建勇.3D打印聚醚醚酮及其复合材料修复骨缺损的应用研究[J].中国组织工程研究，2020，24（10）：1622-1628.

[24]赵靖，王笛，刘继全，等.3D打印技术在医学领域应用的现状及问题[J].中国现代医学杂志，2017，27（12）：71-74.

[25]Wong K C.3D-printed patient-specific applications in orthopedics[J].Orthop Res Rev，2016，8（15）：57-66.

[26]贺昆，张永红，牛金亮，等.基于MRI的人体指骨间关节三维有限元模型的建立及生物力学分析[J].中国医疗前沿，2010，5（9）：6-7.

（责任编辑：伍利华）

Establishment of Three-Dimensional Finite Element Model and Biomechanical Analysis of Index Finger by Polyetheretherketone

SHI Yun ¹，QIN Zugan ¹，CHEN Cheng²，CHENG Lijia^1，2

（1.School of Mechanical Engineering，Chengdu University，Chengdu 610106，China; 2.School of Preclinical Medicine，Chengdu University，Chengdu 610106，China）

Abstract：The three-dimensional finite element model of index finger was established，and the attributes of polyetheretherketone material were assigned to analyze the stress distribution of each phalanx of index finger under different loads.A healthy male volunteer was selected to perform CT scanning on his left hand.The scan data was imported into Mimics 21.0 software to obtain the point cloud data of index finger，and then imported into UG 12.0 software for reversed processing to establish index finger model.Finally，the data were put into the Workbench 2020R1 software.Ten，twenty，thirty，forty and fifty Newton loads were applied along the longitudinal direction of phalanges to analyze the stress distribution of index finger bones.The results show that a high-precision three-dimensional finite element model of index finger has been established，with a total of 551 410 nodes and 392 047 tetrahedral elements.With the increase of load，the peak is stress of index finger bones ranges from 0.37 to 5.96 MPa，and the overall stress of the distal phalanx was significantly higher than that of other bones，of which the proximal phalanx suffers the smallest amount of stress.It can be concluded that Mimics and Workbench software can be used to establish the index finger finite element model with high accuracy so as to conduct qualitative and quantitative analysis of index finger biomechanics，providing a theoretical basis for polyether ether ketone material as index finger replacement material.

Key words：polyetheretherketone;index finger;finite element analysis;three-dimensional reconstruction;biomechanics