中国组织工程研究 ›› 2025, Vol. 29 ›› Issue (33): 7223-7230.doi: 10.12307/2025.914

• 骨与关节有限元分析Finite element analysis of bones and joints • 上一篇    下一篇

内、外踝韧带缺损对邻近核心肌腱生物力学影响的有限元分析

徐天宇1,陈默迪1,谢明茹1,叶兴华2,潘朝晖2   

  1. 1山东第二医科大学临床医学院,山东省潍坊市  261053;2陆军第八十集团军医院骨科研究所,山东省潍坊市  261053
  • 收稿日期:2024-09-19 接受日期:2024-11-19 出版日期:2025-11-28 发布日期:2025-04-12
  • 通讯作者: 潘朝晖,博士,主任医师,硕士生导师,陆军第八十集团军医院骨科研究所,山东省潍坊市 261053
  • 作者简介:徐天宇,男,1998年生,山东第二医科大学在读硕士,主要从事显微重建外科及数字骨科研究。
  • 基金资助:
    潍坊市卫健委科研项目(WFWSJK-2023-040),项目负责人:潘朝晖

Finite element analysis of biomechanical effect of medial or lateral malleolar ligament defects on its neighboring core tendons

Xu Tianyu1, Chen Modi1, Xie Mingru1, Ye Xinghua2, Pan Zhaohui2   

  1. 1School of Clinical Medicine, Shandong Second Medical University, Weifang 261053, Shandong Province, China; 2Institute of Orthopedics and Traumatology of Chinese PLA, Hospital of Chinese PLA 80th Group Army, Weifang 261053, Shandong Province, China
  • Received:2024-09-19 Accepted:2024-11-19 Online:2025-11-28 Published:2025-04-12
  • Contact: Pan Zhaohui, MD, Chief physician, Master’s supervisor, Institute of Orthopedics and Traumatology of Chinese PLA, Hospital of Chinese PLA 80th Group Army, Weifang 261053, Shandong Province, China
  • About author:Xu Tianyu, Master candidate, School of Clinical Medicine, Shandong Second Medical University, Weifang 261053, Shandong Province, China
  • Supported by:
    Weifang Municipal Health Commission Research Project, No. WFWSJK-2023-040 (to PZH)

摘要:


文题释义:

有限元分析:以有限数量的相互作用单元模拟复杂系统,将求解域划分为多个小域,各域内采用近似解,通过求解整体条件获得系统解。
邻近核心肌腱:在空间位置上相互靠近或接近,对于维持身体的平衡、稳定性和姿势控制至关重要的肌腱。


背景:高能量损伤导致的内、外踝韧带断裂、缺损是足踝外科常见难题之一,邻近肌腱常作为移植材料重建缺损,探究韧带缺损状态下邻近核心肌腱的力学特性对于选择肌腱供区具有重要意义。

目的:探究踝关节活动过程中,内、外踝韧带不同程度缺损时,其各自邻近核心肌腱间的相互作用及生物力学变化。
方法:基于一名50岁健康男性左足的CT影像数据,通过MIMICS软件提取并构建表面立体模型,利用Geomagic Wrap进行拟合曲面,并在SOLIDWORKS中构建包含不同程度韧带缺损的骨-软骨-韧带-肌腱足踝复合体模型。采用Ansys Workbench进行有限元分析,通过前抽屉实验验证模型有效后,模拟踝关节内、外翻及内、外旋条件,观察腓骨短、长肌腱与胫前、后肌腱的Von-Mises应力最大值变化及分布特征。

结果与结论:①前抽屉实验中距骨最大位移5.208 5 mm,与文献报道的数据相似,验证模型有效;②4种工况下,踝关节单侧单束韧带缺损对邻近核心肌腱的应力最大值影响较小,单侧多束韧带缺损则显著增加肌腱应力最大值;除胫前肌腱各段应力载荷均偏高外,腓骨短、长肌腱及胫后肌腱的高应力区域均集中在肌腱远端;③外踝韧带一旦缺损,内翻、内旋时腓骨短肌腱的应力最大值及变化量均高于腓骨长肌腱;距腓前韧带缺损时内翻,腓骨短肌腱应力最大值较正常相比增加0.951 2 MPa,腓骨长肌腱应力最大值降低0.065 1 MPa;跟腓韧带缺损时内旋,腓骨短肌腱应力最大值增加2.352 9 MPa,腓骨长肌腱应力最大值降低0.269 2 MPa;④提示内踝韧带缺损时外翻、外旋,胫前、后肌腱的应力最大值变化复杂,取决于缺损韧带类型,尤其全层韧带缺损将显著增加二者应力最大值。

https://orcid.org/0009-0003-5851-2818(徐天宇)

中国组织工程研究杂志出版内容重点:人工关节;骨植入物;脊柱;骨折;内固定;数字化骨科;组织工程

关键词: 踝关节, 肌腱, 韧带缺损, 邻近核心肌腱, 有限元分析, 生物力学

Abstract: BACKGROUND: Ligament rupture and defect of the lateral or medial malleolus caused by high-energy injuries are common challenges in foot and ankle surgery. Their neighboring core tendons are often used as grafts to reconstruct the deficient ligaments. It is of paramount importance to investigate the mechanical properties of such tendons in the context of ligament defects to provide a suitable donor tendon.
OBJECTIVE: To investigate the interactive dynamics and biomechanical alterations among their core tendons during ankle joint motions under varying degrees of lateral or medial malleolar ligament defect.
METHODS: Based on CT imaging data of the left foot of a 50-year-old healthy male, a surface stereolithography model was extracted and constructed using MIMICS software. After Geomagic Wrap software was employed to fit the surfaces, a bone-cartilage-ligament-tendon ankle complex model incorporating varying degrees of ligament deficiencies was created within SOLIDWORKS software. Finite element analysis was then conducted using Ansys Workbench software, and the model’s validity was verified through a simulated anterior drawer test. Following validation, the mechanical response of the ankle under the conditions of internal and external rotation, as well as inversion and eversion, was simulated. The variation and distribution patterns of the maximum Von-Mises stress in the peroneus brevis and longus tendons, as well as the anterior and posterior tibial tendons, were observed.
RESULTS AND CONCLUSION: (1) In the anterior drawer test, the maximum talar displacement reached 5.208 5 mm, which was similar to the data in the previous literature, thereby the effectiveness of the model was validated. (2) Under four loading conditions, the defect of unilateral single-bundle ligaments exerted minimal influence on the maximum stress of adjacent core tendons, whereas the defect of unilateral multi-bundle ligament significantly increased the maximum tendon stress. Except for the consistently high stress across segments of the anterior tibial tendon, the high-stress regions in the long and short peroneal tendons and posterior tibial tendon were concentrated at their distal ends near the insertions. (3) Regarding the defect of the lateral malleolar ligament, the maximum stress and its variation in the peroneus brevis tendon during inversion and internal rotation were higher than those in the peroneus longus tendon. During inversion under the condition of the defect of the anterior talofibular ligament, the maximum stress in the short peroneal tendon increased by 0.951 2 MPa compared to that of normal condition, while that in the long peroneal tendon decreased by 0.065 1 MPa. Under the condition of the defect of the calcaneofibular ligament during internal rotation, the maximum stress in the short peroneal tendon increased by 2.352 9 MPa, while the maximum stress in the long peroneal tendon decreased by 0.269 2 MPa. (4) During eversion and external rotation under the defect of medial malleolar ligament, the variations in the maximum stress of the anterior and posterior tibial tendons were complex and depended on the type of ligament defect. Notably, full-thickness ligament defect significantly augmented the maximum stress in both tendons.

Key words: ankle joint, tendon, ligament defect, neighboring core tendon, finite element analysis, biomechanics

中图分类号: