三维有限元生物力学分析后交叉韧带断裂后膝关节内翻及外旋时的应力

doi:10.3969/j.issn.2095-4344.2017.07.010

摘要/Abstract

摘要：

文章快速阅读：

文题释义：

膝关节韧带损伤：①内侧副韧带损伤：为膝外翻暴力所致，即膝伸直位，膝或腿部外侧受强大暴力打击或重压，使膝过度外展，内侧副韧带可发生部分或完全断裂；②外侧副韧带损伤：主要为膝内翻暴力所致，即膝或腿部内侧受暴力打击或重压，使膝过度内收，外侧副韧带可发生部分或完全断裂。在严重创伤时，侧副韧带、十字韧带和半月板可同时损伤；③前交叉韧带损伤：膝关节伸直位下内翻损伤和膝关节屈曲位下外翻损伤都可以使前交叉韧带断裂；④后交叉韧带损伤：无论膝关节属于屈曲位或伸直位，来自前方的使胫骨上端后移的暴力都可以使后交叉韧带断裂。

后交叉韧带：后交叉韧带起自股骨内侧髁的外侧面，斜向后下方，止于胫骨髁间隆起的后部和外侧半月板的后角。当膝关节活动时，后交叉韧带能防止腿骨向后移位外，还可限制膝关节的过伸、过屈及旋转活动，交叉韧带损伤常与胫侧副韧带或半月板损伤同时发生。

摘要

背景：后交叉韧带损伤合并膝关节周围其他韧带损伤的发生率相对较高，然而这些韧带之间内在动态应力变化情况的报道目前尚属罕见。

目的：模拟正常及后交叉韧带完全断裂(Ⅲ度损伤)的膝关节三维有限元模型，计算并分析正常膝关节模型和后交叉韧带完全断裂模型中内侧副韧带、外侧副韧带和前交叉韧带的应力分布及变化情况。

方法：采集正常志愿者右膝关节的CT及MRI 图像，利用E-feature Biomedical Modeler，ANSYS等软件建立正常膝关节三维有限元模型，模拟正常膝关节屈曲0°、30°、60°时施加内外翻及内外旋力矩，分析后交叉韧带、外侧副韧带、前交叉韧带及内侧副韧带的应力，模拟后交叉韧带完全断裂时分析外侧副韧带、前交叉韧带和内侧副韧带的应力。

结果与结论：三维有限元生物力学分析数据显示后交叉韧带断裂后外侧副韧带在膝关节内翻及外旋运动时应力明显增大，前交叉韧带及内侧副韧带的应力变化较小。提示后交叉韧带断裂可能造成膝关节其他韧带的继发损伤，进一步支持了恢复膝关节动态稳定的重要性。

中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程
ORCID: 0000-0002-0616-9347(谢强)

关键词: 骨科植入物, 数字化骨科, 膝关节, 后交叉韧带, 前交叉韧带, 外侧副韧带, 内侧副韧带, 三维有限元, 泊松比, 生物力学

Abstract:

BACKGROUND: The incidence of posterior cruciate ligament injury combined with other knee ligament injuries is relatively high, but studies on the internal dynamic stress change of these ligaments are still rare.

OBJECTIVE: To simulate three-dimensional finite element model of the normal human knee and the complete posterior cruciate ligament rupture (III injury), calculate and analyze the stress distribution and changes of the medial collateral ligament, the lateral collateral ligament and the anterior cruciate ligament in the normal knee model and the complete posterior cruciate ligament rupture model.

METHODS: We collected CT and MRI images of the healthy male right knee, and established three-dimensional finite element model of the normal human knee by using E-feature Biomedical Modeler, and ANSYS software. The simulation of a flexion moment at 0°, 30°, 60°, varus/valgus and internal/external rotation torque was applied on the established three-dimensional finite element model of the normal knee. The stress of posterior cruciate ligament, lateral collateral ligament, anterior cruciate ligament and medial collateral ligament was analyzed. During simulation of complete posterior cruciate ligament rupture, the stress of lateral collateral ligament, anterior cruciate ligament and medial collateral ligament was analyzed.

RESULTS AND CONCLUSION: Three-dimensional finite element biomechanical analysis data revealed that after posterior cruciate ligament rupture, the stress of lateral collateral ligament during inversion and external rotation of the knee joint obviously increased, but the stress of anterior cruciate ligament and medial collateral ligament slightly altered. These findings suggested that posterior cruciate ligament rupture may cause secondary injury to other ligaments of the knee joint, and further support the importance of the dynamic stability of the knee joint.

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

Key words: Knee Joint, Ligaments, Finite Element Analysis, Tissue Engineering

中图分类号:

R318

谢强. 三维有限元生物力学分析后交叉韧带断裂后膝关节内翻及外旋时的应力[J]. 中国组织工程研究, 2017, 21(7): 1036-1040.

Xie Qiang. Three-dimensional finite element model for biomechanical analysis of stress in knee inversion and external rotation after posterior cruciate ligament rupture[J]. Chinese Journal of Tissue Engineering Research, 2017, 21(7): 1036-1040.

图/表（结果） 1

参考文献

[1] Brekelmans WA, Poort HW, Slooff TJ. A new method to analyse the mechanical behaviour of skeletal parts. Acta Orthop Scand. 1972;43(5):301-317.

[2] Rybicki EF, Simonen FA, Weis EB Jr. On the mathematical analysis of stress in the human femur. J Biomech. 1972; 5(2):203-215.

[3] Zelle J, Van der Zanden AC, De Waal Malefijt M, et al. Biomechanical analysis of posterior cruciate ligament retaining high-flexion total knee arthroplasty. Clin Biomech (Bristol, Avon). 2009;24(10):842-849.

[4] Beillas P, Lee SW, Tashman S, et al. Sensitivity of the tibio-femoral response to finite element modeling parameters. Comput Methods Biomech Biomed Engin. 2007;10(3): 209-221.

[5] Mesfar W, Shirazi-Adl A. Computational biomechanics of knee joint in open kinetic chain extension exercises. Comput Methods Biomech Biomed Engin. 2008;11(1):55-61.

[6] Chantarapanich N, Nanakorn P, Chernchujit B, et al. A finite element study of stress distributions in normal and osteoarthritic knee joints. J Med Assoc Thai. 2009;92 Suppl 6:S97-103.

[7] Guess TM, Thiagarajan G, Kia M, et al. A subject specific multibody model of the knee with menisci. Med Eng Phys. 2010;32(5):505-515.

[8] Gíslason MK, Stansfield B, Nash DH. Finite element model creation and stability considerations of complex biological articulation: The human wrist joint. Med Eng Phys. 2010; 32(5):523-531.

[9] Frehill B, Crocombe A, Cirovic S, et al. Initial stability of type-2 tibial defect treatments. Proc Inst Mech Eng H. 2010;224(1): 77-85.

[10] 刘安庆,张银光,王春生,等.人股骨生物力学特性的三维有限元分析[J].西安医科大学学报,2001,22(3):242-244.

[11] 邹磊,詹朝双,史本龙,等.国人正常膝关节几何学参数测量及临床意义[J].解剖与临床,2010,15(4):243-246.

[12] Ramaniraka NA, Saunier P, Siegrist O, et al. Biomechanical evaluation of intra-articular and extra-articular procedures in anterior cruciate ligament reconstruction: a finite element analysis. Clin Biomech (Bristol, Avon). 2007;22(3):336-343.

[13] LeRoux MA, Setton LA. Experimental and biphasic FEM determinations of the material properties and hydraulic permeability of the meniscus in tension. J Biomech Eng. 2002;124(3):315-321.

[14] Li G, Lopez O, Rubash H. Variability of a three-dimensional finite element model constructed using magnetic resonance images of a knee for joint contact stress analysis. J Biomech Eng. 2001;123(4):341-346.

[15] Ramaniraka NA, Terrier A, Theumann N, et al. Effects of the posterior cruciate ligament reconstruction on the biomechanics of the knee joint: a finite element analysis. Clin Biomech (Bristol, Avon). 2005;20(4):434-442.

[16] Peña E, Calvo B, Martinez MA, et al. Influence of the tunnel angle in ACL reconstructions on the biomechanics of the knee joint. Clin Biomech (Bristol, Avon). 2006;21(5):508-516.

[17] Moglo KE, Shirazi-Adl A. Biomechanics of passive knee joint in drawer: load transmission in intact and ACL-deficient joints. Knee. 2003;10(3):265-276.

[18] Peña E, Calvo B, Martínez MA, et al. A three-dimensional finite element analysis of the combined behavior of ligaments and menisci in the healthy human knee joint. J Biomech. 2006; 39(9):1686-1701.

[19] LaPrade RF, Muench C, Wentorf F, et al. The effect of injury to the posterolateral structures of the knee on force in a posterior cruciate ligament graft: a biomechanical study. Am J Sports Med. 2002;30(2):233-238.

[20] Song Y, Debski RE, Musahl V, et al. A three-dimensional finite element model of the human anterior cruciate ligament: a computational analysis with experimental validation. J Biomech. 2004;37(3):383-390.

引言

材料方法

1.1 设计 三维有限元分析。

1.2 时间及地点 试验于2013年7月在中南大学湘雅医学院完成。

1.3 对象 选取1名健康的成年男性志愿者，年龄25岁，身高173 cm，体质量60 kg，体型标准，无膝关节疼痛及外伤史，符合中国正常人膝关节解剖数据^[11]。试验已经获得参与者知情同意。

1.4 方法

1.4.1 数据采集利用Philip Tomoscans R7000 64 SCT扫描机和Siemens 1.5T超导磁共振采集志愿者膝关节CT图像及MRI图像。通过CT扫描获骨骼的边界，通过MRI扫描获得韧带、半月板和软骨边界，将数据导出并刻盘存为DICOM格式，来获得模型重建的数据来源。

1.4.2 建立正常人膝关节三维有限元模型

膝关节各结构三维实体模型的建立：将Dicom格式的 CT和MRI扫描数据分别导入 E-feature Biomedical Modeler软件(辉擎信息科技公司)，利用系统中的图像分割和组织自动提取功能，进行精确三维分割与模型重建从CT扫描数据重建出股骨、胫腓骨及髌骨三维模型(图1A)；从MRI扫描数据中重建出主要韧带、半月板、和软组织表面三维模型(图1B)。

利用E-feature Biomedical Modeler软件中的表面网格编辑工具对各部分三维模型，特别是从MRI扫描数据重建的三维模型进行必要的编辑修改，使模型更光滑，柔顺，得到高质量的表面模型(图2)。

膝关节三维有限元模型网格化：在 E-feature Biomedical Modeler中将模型生成为高质量体网格，并将各个部位集合成为膝关节三维有限元模型，再将Biomedical Modeler系统中膝关节三维有限元模型导入ANSYS 12.0有限元分析软件，将线性的四面体单元全部转换为非线性的二次10节点四面体和6节点壳单元，以提高数值计算的精度。

设置材料属性：根据实验要求及参照文献[12-14]定义骨组织为各向同性的弹性材料，关节软骨、半月板为均匀、连续、各向同性的弹性材料，韧带为均匀、连续、各向同性的粘弹性材料，并设置各材料的参数(表1)。

边界设定及定义接触区：设定关节软骨内表面与骨组织的表面相固定，半月板的前角、后角以及内侧半月板外周边与胫骨平台的边缘相固定，以此来模拟半月板在胫骨平台的附着。用80个线性弹簧单元来模拟股四头肌的作用，使其与股骨走行一致。在对模型进行加载后，模型各结构中会出现接触区域，在进行分析前有必要对接触区域进行假设，将此区域均设定为非线性无摩擦接触，设定14处可能的接触区：髌股关节1处，前交叉韧带和后交叉韧带之间1处，膝关节内外侧间室各2处(半月板-胫骨软骨、股骨软骨-半月板)，膝关节内/外侧副韧带、前/后交叉韧带与股骨及胫骨(或腓骨)间各4处。

加载条件：设定胫腓骨捆绑并固定，约束其在X、Y、Z三个方向的活动自由度和旋转，将股骨设为承载面，胫骨下端固定，设置为与下表面牢固结合，无移动。根据正常膝关节生理运动状态下会出现内外翻及内外旋的运动，遂向模型施加弯曲力矩以模拟膝关节屈曲0°、30°、60°的同时再向模型分别施加20 N•m的内外翻力矩、8 N•m的内外旋力矩来模拟胫骨内外翻和内外旋的生理运动^[15-18]。

1.4.3 建立后交叉韧带断裂模型在已建立的正常膝关节三维有限元模型(图3)的基础上除去后交叉韧带，模拟后交叉韧带完全断裂(Ⅲ度损伤，见图4)，并模拟膝关节分别屈曲0°、30°、60°。

1.4.4 应力分析

模拟正常膝关节模型：在屈曲0°、30°、60°时，分别施加20 N•m的内外翻力矩、8 N•m的内外旋力矩，记录后交叉韧带、外侧副韧带、前交叉韧带及内侧副韧带的平均应力值及等效应力分布图。

模拟后交叉韧带完全断裂：在屈曲0°、30°、60°时，分别施加20 N•m的内外翻力矩、8 N•m的内外旋力矩，记录并分析外侧副韧带、前交叉韧带及内侧副韧带的平均应力值及等效应力分布图。

1.5 主要观察指标 正常膝关节模型与后交叉韧带完全断裂模型在施加不同外力时外侧副韧带、前交叉韧带、内侧副韧带的应力变化及等效应力分布图。

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

讨论

3.1 有限元分析结果 在后交叉韧带断裂后膝关节在做内翻运动时屈曲到60°位外侧副韧带的应力变化更为明显。LaPrade等^[19]报道外侧副韧带在对抗膝关节内翻运动时与后交叉韧带协同发挥作用，在屈曲60°时其应力达最大。在后交叉韧带断裂后膝关节在做外旋运动时屈曲到30°位外侧副韧带的应力变化更为明显，这与Song等^[20]的研究结果相一致。在这2种应力情况下外侧副韧带的股骨端应力值较高，可能成为后交叉韧带断裂后继发外侧副韧带损伤的最主要部位。而20 N•m内外翻力矩和8 N•m内外旋力矩作用下，在后交叉韧带断裂模型中前交叉韧带及内侧副韧带在不同屈曲角度时的等效平均应力值与正常膝关节模型差别很小。

3.2 有限元分析在骨科领域应用的不足之处 作者所建立的模型在进行生物力学分析时还存在一定的不足之处：①所建立的模型中并没有包括肌肉、关节囊和后方稳定结构等膝关节的其他解剖结构；②模型只屈曲到60°，没有再模拟更多屈曲角度的活动；③为了简化建立模型过程，模型中所设定各个结构的材料属性相对单一；④只模拟了膝关节在不同屈曲角度的静态力学情况，对于膝关节动态的屈曲活动还不能真实反映。

3.3 总结 作者成功建立正常人膝关节三维有限元模型的，去除后交叉韧带用来模拟后交叉韧带完全断裂，在正常膝关节模型和后交叉韧带完全断裂模型上，模拟膝关节分别屈曲 0°、30°、60°，同时分别施加20 N•m的内外翻力矩、8 N•m的内外旋力矩，通过数据分析发现在后交叉韧带完全断裂模型中外侧副韧带在施加内翻及外旋应力时，其平均应力值与正常膝关节模型比较明显增大，前交叉韧带及内侧副韧带在膝关节运动时应力变化较小。这说明后交叉韧带完全断裂后会影响膝关节内在动态稳定性，造成膝关节内其他结构的慢性损害，这也为证明预防后交叉韧带断裂后继发损伤及恢复重建膝关节的稳定的重要性提供理论依据。

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