韧带损伤影响腰椎稳定性的生物力学特征

doi:10.12307/2024.070

中国组织工程研究 ›› 2024, Vol. 28 ›› Issue (21): 3286-3292.doi: 10.12307/2024.070

• 骨与关节生物力学 bone and joint biomechanics • 上一篇下一篇

韧带损伤影响腰椎稳定性的生物力学特征

李银倩1，吕杰1，丁立军2，王多多3，郭潘靖3，曹金凤4，周楠5，吕强5

上海健康医学院，1康复学院，2医疗器械学院，上海市 201318；3 上海理工大学健康科学与工程学院，上海市 200093；4上海大学期刊社，上海市 200444；5上海中医药大学附属岳阳中西医结合医院，上海市 200437

收稿日期:2023-03-23 接受日期:2023-05-10 出版日期:2024-07-28 发布日期:2023-09-27
通讯作者: 吕杰，博士，副教授，硕士生导师，上海健康医学院康复学院，上海市 201318
作者简介:李银倩，女，1998年生，山东省淄博市人，汉族，2023年上海理工大学(上海健康医学院联合培养)毕业，硕士，主要从事生物力学方面的研究。
基金资助:
国家自然科学基金面上项目(82074571)，项目负责人：周楠

Biomechanical characteristics of ligament injury affecting lumbar spine stability

Li Yinqian1, Lyu Jie1, Ding Lijun2, Wang Duoduo3, Guo Panjing3, Cao Jinfeng4, Zhou Nan5, Lyu Qiang5

1College of Rehabilitation Sciences, 2College of Medical Instrumentation, Shanghai University of Medicine & Health Sciences, Shanghai 201318, China; 3School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China; 4Periodicals Agency of Shanghai University, Shanghai 200444, China; 5Yueyang Hospital of Integrated Traditional Chinese and Western Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai 200437, China

Received:2023-03-23 Accepted:2023-05-10 Online:2024-07-28 Published:2023-09-27
Contact: Lyu Jie, MD, Associate professor, Master’s supervisor, College of Rehabilitation Sciences, Shanghai University of Medicine & Health Sciences, Shanghai 201318, China
About author:Li Yinqian, Master, College of Rehabilitation Sciences, Shanghai University of Medicine & Health Sciences, Shanghai 201318, China
Supported by:
National Natural Science Foundation of China (General Program), No. 82074571 (to ZN)

摘要/Abstract

摘要：

文题释义：

腰椎失稳：由于腰椎活动超出正常生理活动范围的现象即为腰椎失稳，这将引起下腰痛、脊柱畸形和椎间盘病变等一些列症状。所有的腰椎软组织都有助于稳定腰椎，因此任何组织结构的退行性、创伤性或破坏性损伤都会导致一定程度的腰椎不稳定性。
有限元分析法：利用数学近似的方法对现实物理系统进行模拟，即将连续求解区域离散为有限个单元，对单元进行力学分析，最后整合。它可以很好地应用于骨科生物力学研究，是在计算机辅助下模拟脊柱病理最有效的工具。

背景：韧带是维持腰椎稳定的重要结构，由于产生的机械应力这些结构容易退化，而对于韧带损伤的研究较少。

目的：基于三维有限元模型模拟韧带损伤，确定各个韧带在运动状态下的活动度及应力方面的作用，系统探讨韧带受损对腰椎稳定性的影响和生物力学意义。
方法：建立腰椎L4-L5有限元模型，约束L5下表面所有自由度，施加5 N·m力矩模拟腰椎的各个运动状态，通过改变韧带弹性模量模拟韧带渐进式受损。

结果与结论：①大部分韧带损伤伸展活动度都出现了明显的增大；相比之下侧弯活动度无明显变化。②关节囊韧带受损后前屈、后伸和侧弯的活动度明显增大；前纵韧带受损后伸展活动度明显增大；横突韧带损伤导致侧弯活动度明显增大。③单韧带损伤后，伸展时活动度变化最为显著；单韧带受损后剩余韧带应力增大，尤其是关节囊韧带；整体来看棘间韧带和黄韧带应力变化最不明显；韧带应力在侧弯状态下变化最小，而扭转状态下变化最大。④韧带受损并未导致椎间盘压力发生重大变化，提示韧带损伤导致腰椎活动度增大，影响腰椎稳定性。⑤关节囊韧带稳定前屈、后伸和侧弯；前纵韧带稳定后伸；横突韧带稳定侧弯。⑥后伸运动对韧带损伤敏感，可通过后伸运动检查韧带病理变化。⑦单韧带损伤后剩余韧带给予应力补偿以维持腰椎稳定性。⑧棘间韧带和黄韧带损伤对周围韧带影响最小，关节囊韧带损伤对周围韧带影响最大。⑨韧带损伤对侧弯运动下的剩余韧带应力影响最小，而对扭转运动下的韧带应力影响最大，韧带损伤患者应避免扭转运动；韧带损伤不影响椎间盘压力。

https://orcid.org/0000-0001-6342-2735 (李银倩)

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

关键词: 有限元, 腰椎, 韧带, 渐进性损伤, 生物力学

Abstract: BACKGROUND: Ligaments are important structures in maintaining the stability of the lumbar spine, and these structures are prone to degradation due to the generated mechanical stress. However, there are few studies on ligament injuries.
OBJECTIVE: To determine the range of motion and stress of each ligament in the state of motion based on a three-dimensional finite element model, and to systematically explore the influence of ligament injury on the stability of the lumbar spine and its biomechanical significance.
METHODS: The L4-L5 lumbar finite element model was established. All free forces on the lower surface of L5 were constrained, and a torque of 5 N·m was applied to simulate the motion states of the lumbar spine. Progressive ligament damage was simulated by changing Young’s modulus of the ligament.
RESULTS AND CONCLUSION: (1) The extension range of motion of most ligament injuries increased significantly. In contrast, there was no significant change in the lateral bending range of motion. (2) The range of motion of capsular ligament injury increased significantly in flexion, extension and lateral bending. The extension range of motion increased significantly after the anterior longitudinal ligament injury. Intertransverse ligament injury resulted in a significant increase in the lateral bending range of motion. (3) After a single ligament injury, the most significant change in range of motion was observed during extension. After a single ligament injury, the stress of the remaining ligaments increased, especially the capsular ligament. The stress changes of the interspinous ligament and ligamentum flavum were the least obvious. Ligament stress changes least in lateral bending and most in torsion. (4) Ligament damage did not result in major changes in disc pressure, indicating that ligament injury leads to an increased range of motion of the lumbar spine and affects the stability of the lumbar spine. (5) Capsular ligament was stable in flexion, extension and lateral bending. The anterior longitudinal ligament showed a stable extension; the intertransverse ligament stabilizes the lateral curve. (6) Extension exercise is sensitive to a ligament injury, and the pathological changes of the ligament can be examined by extension exercise. (7) Stress compensation was given to the remaining ligaments to maintain the stability of the lumbar spine after a single ligament injury. (8) Interspinous ligament and ligamentum flavum injuries have the least impact on the peripheral ligaments, while capsular ligament injuries have the greatest impact on the peripheral ligaments. (9) The ligament injury has the least effect on the residual ligament stress during lateral bending exercise, while it has the greatest effect on the ligament stress during the twisting exercise. Patients with ligament injury should avoid twisting exercises. Ligament injuries do not affect disc pressure.

Key words: finite element, lumbar vertebrae, ligament, progressive injury, biomechanics

中图分类号:

李银倩, 吕杰, 丁立军, 王多多, 郭潘靖, 曹金凤, 周楠, 吕强. 韧带损伤影响腰椎稳定性的生物力学特征[J]. 中国组织工程研究, 2024, 28(21): 3286-3292.

Li Yinqian, Lyu Jie, Ding Lijun, Wang Duoduo, Guo Panjing, Cao Jinfeng, Zhou Nan, Lyu Qiang. Biomechanical characteristics of ligament injury affecting lumbar spine stability[J]. Chinese Journal of Tissue Engineering Research, 2024, 28(21): 3286-3292.

图/表（结果） 3

2.1 椎骨活动度结果韧带的损伤通常会导致腰椎活动度的增加。图5显示了140个活动度结果，在“T型图”上，表示了图表中每个分支的加载模式(侧弯和旋转的左右运动取平均值)。列出了不同程度上由于不同韧带受损而导致的腰椎(L4-L5)活动度的变化百分比。每个单元格显示了某运动状态下该韧带损伤一定程度时的腰椎节段活动度。比如，在前屈状态下，关节囊韧带损伤50%时腰椎节段活动度变化百分比为9.65%。

此外，这些单元格中的文本表示4种运动状态下的活动度随韧带受损程度的加深而增大。其中，后伸时活动度变化趋势最为明显，大部分韧带损伤后活动度都出现了明显的增大；相比之下侧弯活动度无明显变化，大部分韧带受损后活动度几乎不变，整体来看侧弯最不容易受韧带损伤的影响。关节囊韧带受损导致前屈的最大变化的发生，相比健康韧带，关节囊韧带受损程度50%的情况下活动度增加9.65%。对于后伸而言，关节囊韧带和前纵韧带受损对活动度影响最大，其变化百分比分别为32.37%和8.95%。扭转变化最为明显只发生在关节囊韧带受损后，活动度增大10.16%。侧弯活动度的变化最不明显，最大变化发生在关节囊韧带受损后，为6.15%。
2.2 韧带应力结果韧带的损伤通常会导致剩余韧带最大应力增加。图6显了每组模拟结果，受损韧带显示在每个“T形图”的中心。在T型图上，表示了图表中每个分支的加载模式(横向弯曲和轴向旋转的左右运动取平均值)，列出了不同程度上由于韧带受损而导致的剩余韧带最大应力的变化百分比。每个单元格显示了某运动状态下该韧带损伤一定程度时的剩余韧带最大应力变化百分比。比如，T型图中标记为“XALL”，这意味着这一组实验中受损韧带为前纵韧带。这些单元格中的文本表明，在前屈状态下，前纵韧带损伤10%-50%时关节囊韧带最大应力变化百分比为0.55%-2.98%。
结果显示一种韧带受损将导致剩余韧带应力增加。任何一种韧带受损都会在扭转期间导致关节囊韧带最大应力明显增大。关节囊韧带受损会在扭转期间增大棘上韧带、棘间韧带和黄韧带的应力(24.56%，15.39%，12.99%)；在前屈和后伸期间则增大所有韧带应力。韧带受损后棘间韧带在侧弯运动状态下出现了不同程度的应力减小，然而棘间韧带应力明显增大发生在横突韧带受损后的侧弯状态下。棘上韧带在后纵韧带受损50%时的扭转状态下有一个明显应力的突变，由-0.31%增加到20.32%。整体来看韧带应力在侧弯状态下最不容易受影响，最大应力变化百分比大部分都不超过1%。而扭转状态下的剩余韧带最大应力变化最大。

如图7所示，韧带受损并未导致椎间盘的压力发生重大变化。在负载情况下椎间盘通常承受80%-82%的载荷。尽管周围组织中的载荷分布发生了变化，但椎间盘内的压力保持稳定。

参考文献

[1] SHARMA M, LANGRANA NA, RODRIGUEZ J. Role of ligaments and facets in lumbar spinal stability. Spine (Phila Pa 1976). 1995;20(8): 887-900.
[2] GILLESPIE KA, DICKEY JP. Biomechanical role of lumbar spine ligaments in flexion and extension: determination using a parallel linkage robot and a porcine model. Spine (Phila Pa 1976). 2004; 29(11):1208-1216.
[3] 杨雄.基于韧带损伤腰椎力学失稳模型的建立与应用[D].重庆:重庆大学,2019.
[4] HEUER F, SCHMIDT H, KLEZL Z, et al. Stepwise reduction of functional spinal structures increase range of motion and change lordosis angle. J Biomech. 2007;40:271-280.
[5] HEUER F, SCHMIDT H, CLAES L, et al. Stepwise reduction of functional spinal structures increase vertebral translation and intradiscal pressure. J Biomech. 2007;40:795-803.
[6] 吴聪聪. 胸腰段脊柱后方韧带复合体序贯损伤对其稳定性的有限元分析[D]. 温州:温州医科大学,2019.
[7] NG HW, TEO EC, ZHANG Q. Influence of cervical disc degeneration after posterior surgical techniques in combined flexion-extension - a nonlinear analytical study. J Biomech Eng. 2005;127(1):186-192.
[8] MACKERLE J. Finite element modeling and simulations in orthopedics: a bibliography 1998-2005. Comput Methods Biomech Biomed Eng. 2006;9(3):149-199.
[9] SCHMIDT H, HEUER F, DRUMM J, et al. Application of a calibration method provides more realistic results for a finite element model of a lumbar spinal segment. Clin Biomech (Bristol, Avon). 2007;22(4):377-384.
[10] BOWDEN AE, GUERIN HL, VILLARRAGA ML, et al. Quality of motion considerations in numerical analysis of motion restoring implants of the spine. Clin Biomech (Bristol, Avon). 2008;23(5):536-544.
[11] AYTURK UM, PUTTLITZ CM. Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine. Comput Methods Biomech Biomed Eng. 2011;14(8):695-705.
[12] UENO K, LIU YK. A three dimensional nonlinear finite element model of lumbar interverbral joint in torsion. J Biomech(Eng). 1987;109:200-209.
[13] GUDAVALLI MR, TRIANO JJ. An analytical model of lumbar motion segment in flexion. J Manipulative Physiol Ther. 1999; 22(4):201-208.
[14] BERMEL EA, BAROCAS VH, ELLINGSON AM. The role of the facet capsular ligament in providing spinal stability. Comput Methods Biomech Biomed Engin. 2018;21(13):712-721.
[15] ALEXANDER T, NIKOALOS M, SAVVAS S, et al. A finite element model technique to determine the mechanical response of a lumbar spine segment under complex loads. Journal of Applied Biomechanics. 2012;28:448-456.
[16] PANJABI MM, GOEL VK, TAKATA K. Physiologic strains in the lumbar spinal ligaments. An in vitro biomechanical study 1981 volvo award in biomechanics. Spine. 1982;7(3):192-203.
[17] ZANDER T, ROHMANN A, BERGMANN G. Analysis of simulated single ligament transection on the mechanical behaviour of a lumbar functional spinal unit. Biomed Tech. 2004;49(1-2):27-32.
[18] HOPKINS TJ, WHITE AA, 3RD. Rehabilitation of athletes following spine injury. Clin Sports Med. 1993;12(3):603-619.
[19] GUO LX, LI WJ. A biomechanical investigation of thoracolumbar burst fracture under vertical impact loads using finite element method. Clin Biomech. 2019;68:29-36.
[20] AFFOLTER C, KEDZIERSKA J, VIELMA T, et al. Estimating lumbar passive stiffness behaviour from subject-specific finite element models and in vivo 6DOF kinematics. J Biomech . 2020;102:109681.
[21] HSIEH YY, KUO YJ, CHEN CH, et al. Biomechanical Assessment of Vertebroplasty Combined with Cement-Augmented Screw Fixation for Lumbar Burst Fractures: A Finite Element Analysis. Appl Sci Basel. 2020; 10(6):14.

[22] 张文,王兰,施勤,等.腰椎行椎间孔入路椎间融合术固定的有限元分析[J].医用生物力学,2014,29(5):405-410.
[23] RASTEGAR S, ARNOUX PJ, WANG XY, et al. Biomechanical analysis of segmental lumbar lordosis and risk of cage subsidence with different cage heights and alternative placements in transforaminal lumbar interbody fusion. Comput Methods Biomech Biomed Eng. 2020;23(9): 456-466.
[24] GILLESPIE KA, DICKEY JP. Biomechanical role of lumbar spine ligaments in flexion and extension: Determination using a parallel linkage robot and a porcine model. Spine. 2004;29(11):1208-1216.
[25] 方新果,赵改平,王晨曦,等. 基于CT图像腰椎L4~L5节段有限元模型建立与分析[J].中国生物医学工程学报,2014,33(4):487-492.
[26] 王宏卫,刘新宇,万熠.人体腰椎L4-L5段有限元模型建立及力学有效性验证[J].医学与哲学(B),2017,38(5):50-53.
[27] 颜文涛,赵改平,方新果,等.人体腰椎L(4~5)节段有限元建模及分析[J].生物医学工程学杂志,2014,31(3):612-618.
[28] 蒋强,丁宇,刘金玉,等.有限元模拟全内镜下精准椎板开窗减压术及生物力学分析[J]. 中国组织工程研究,2020,24(12):1891-1896.
[29] OXLAND TR, LIN RM, PANJABI MM. Three‐dimensional mechanical properties of the thoracolumbar junction. J Orthop Res. 1992;10:573-580.
[30] QIU TX, TEO EC, ZHANG QH. Effect of bilateral facetectomy of thoracolumbar spine T11–L1 on spinal stability. Med Biol Eng Comput. 2006;44:363-370.
[31] 庞胤,尹帅,张海峰,等. 退变对腰椎4-5节段椎间盘生物力学的影响[J]. 山东医药,2021,61(4):55-57.
[32] ALAPAN Y, DEMIR C, KANER T, et al. Instantaneous center of rotation behavior of the lumbar spine with ligament failure. J Neurosurg Spine. 2013;18:617-626.
[33] 牛陆,李娜,柳茵,等.后纵韧带硬化对颈椎力学影响的三维有限元分析[J]. 北京生物医学工程,2018,37(1):40-50.
[34] WIDMER J, ETH MSC, CORNAZ F, et al. Biomechanical contribution of spinal structures to stability of the lumbar spine-novel biomechanical insights. Spine. 2020;20(10):1705-1716.
[35] ELLINGSON AM, SHAW MN, GIAMBINI H, et al. Comparative role of disc degeneration and ligament failure on functional mechanics of the lumbar spine. Comput Methods Biomech Biomed Engin. 2016; 19(9):1009-1018.
[36] NASERKHAKI S, ARJMAND N, SHIRAZI-ADL A, et al. Effects of eight different ligament property datasets on biomechanics of a lumbar L4-L5 finite element model. J Biomech. 2018;70:33-42.
[37] GREGORY A, FORELL V, ANTON E, et al. Biomechanical implications of lumbar spinal ligament transection.Comput Methods Biomech Biomed Eng. 2014; 17(15): 1685-1695.
[38] 付方达,邵晶奕,童培建,等.脊柱轴线力学失稳诱导椎间盘退变小鼠模型的建立及评价[J].中国中医骨伤科杂志,2020,28(1):1-7.
[39] 黄河,李宁宁,胡朝晖,等.纤维环穿刺法与腰椎失稳法建立椎问盘退变模型[J].中国现代医生,2011,49(16):22-24.

引言

腰椎是由各种主动和被动结构组成的复杂结构，允许直立姿势及复杂的运动，并在日常活动中承受载荷。腰椎的主动结构主要是一些肌肉系统，腰椎的被动结构则是一些腰椎间关节及韧带。当主动神经肌肉成分控制运动时，腰椎的被动结构提供稳定性、限制和定义运动范围，并保护椎管内的神经结构[1-2]。除椎骨外，主要的解剖结构是连接椎体的椎间盘、覆盖椎体腹段的前纵韧带、附着于椎体背部的后纵韧带、连接脊柱椎板的黄韧带、横突韧带连接横突、棘间韧带和棘上韧带连接背部棘突、关节囊韧带包裹小关节。不同韧带的力学性质及功能也有所不同，前纵韧带很弱地附着于椎间盘的前面，抵抗椎间盘突出的能力有限；后纵韧带附着于椎间盘的后面和椎体的邻近缘，防止椎间盘向后移动；黄韧带具有较高比例的弹性蛋白，保证了在伸展运动时椎管不会过度弯曲；横突韧带加强椎体间连接的稳定性；棘间韧带和棘上韧带位于运动中心最远端[3]；关节囊韧带包裹在椎间关节的关节囊外部，以便于脊柱运动。由于其解剖位置和特定组成，每个被动结构都有助于在特定运动和特定载荷条件下保持稳定性。由于运动产生的机械应力，这些结构容易退化，从而改变其有效性并导致疼痛或故障，因此，了解每个脊柱结构的特殊功能并理解退行性的病理机制至关重要。在这方面，过去进行了几项横断面研究[4-6]，均在体外对脊柱节段施加弯矩，依次移除韧带结构，并测量活动度的变化。
有限元方法允许改变各种参数，可以系统了解每个结构和条件对脊柱生物力学响应的作用。因此腰椎有限元模型成为预测腰椎正常力学和病理力学的标准工具，它们提供了在不同临床条件下重复检测同一脊柱节段的机会，消除了临床检测和尸体检测中发现的高度受试者特异性变异。先前许多研究者利用腰椎有限元模型进行实验[7-12]。GUDAVALLI等[13]使用有限元模型检查韧带横断导致的屈曲载荷分担和韧带应变。BERMEL等[14]通过有限元模型研究了小关节囊韧带在脊柱稳定性中的作用，随着关节囊韧带的移除，活动度增加，弯曲刚度降低，并观察到运动轴的改变。
正常模型可以反映人体正常生理载荷下的力学性质，但是一旦某个部位发生病变，则腰椎系统内结构的力学性质会随之改变，研究者越来越聚焦于具有临床意义的病理模拟研究。这些研究通过建立病态有限元模型来模拟人体病变，在模型上进行多次重复实验预测不同情况对解剖结构的影响，这对于腰椎疾病的诊断和治疗具有重要意义。ALEXANDER 等[15]利用L3-L4有限元模型研究了前屈状态下韧带响应，发现前屈中的棘上韧带和棘间韧带以及伸展中的关节囊韧带存在失效风险。HEUER等[4-5]建立L4-L5有限元模型，通过逐步减少法依次去除韧带得到腰椎平移和运动范围增加的结论。PANJABI等[16]利用腰椎模型施加载荷计算韧带长度变化，研究了各韧带在各个运动状态下的应变，发现前屈过程中棘间韧带和棘上韧带应变最大；伸展过程中前纵韧带应变最大；侧弯过程中横突韧带、黄韧带和关节囊韧带张力最大；旋转过程中关节囊韧带是张力最大的韧带。ZANDER等[17]使用有限元模型模拟韧带横断导致的节段旋转和韧带载荷分担的变化。GUDAVALLI等[13]使用有限元模型预测韧带横断导致的前屈状态下载荷分担和韧带应变，实验结果表明棘间韧带承受载荷最大，其次是黄韧带、关节囊韧带、横突韧带和棘上韧带。韧带应变表明棘间韧带应变增加最大，其次是棘上韧带、黄韧带、关节囊韧带和横突韧带。然而过去往往直接横断韧带，对于韧带损伤的模拟研究尚且不足。并且，更多学者侧重于研究某种运动模式下的生物力学特性，然而腰椎系统作为一个整体共同分担载荷变化，因此应该进行更为系统的研究。
因此，此次研究的目的有：①开发并验证腰椎(L4-L5)的三维有限元模型；②系统地改变韧带特性，模拟韧带受损；③确定各个韧带在前屈、后伸、侧弯及扭转运动状态下的活动度及应力方面的作用，系统地探讨韧带受损对腰椎稳定性影响和生物力学意义。
此次研究的创新点和特色为：①以往软组织的研究大多集中在椎间盘，或者是单元运动状态和单一韧带，此次研究各运动状态下每个韧带的生物力学特性，是更为系统的研究；②此次研究有限元模型中的韧带建立为六面体单元韧带，更接近人体腰椎解剖学结构，更符合真实人体情况；③对于韧带病理研究，以往研究大多数直接横断韧带，此文则采用改变韧带参数的方法来模拟韧带渐进式损伤，更符合实际。作者假设韧带受损将导致活动度增大，韧带应力增大，韧带受损将影响腰椎稳定性。此外，由于超生理负荷(如运动项目和车辆碰撞)导致的腰椎韧带损伤很少进行手术治疗[18]。脊柱手术的侵入性，加上脊椎疼痛的非特异性使得修复脊椎韧带损伤的手术技术不够成熟。然而，治疗的缺乏并不意味着对韧带损伤后果的忽视，对腰椎韧带损伤后果的进一步了解可以阐明腰椎韧带在临床和腰痛中的作用，可为临床的干预治疗提供合理建议和有力依据。因而从根源上了解韧带受损的生物力学特性具有鲜明的临床意义。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

材料方法

1.1 设计建立腰椎L4-L5有限元模型，约束L5下表面所有自由的，施加5 N·m力矩模拟腰椎各运动状态，改变韧带弹性模量模拟韧带渐进式受损。
1.2 时间及地点实验于2022年9月在上海健康医学院医学影像学院完成。

1.3 材料选择1名合格的健康志愿者的L4/5节段作为观察对象，男性，22岁，对志愿者进行影像检查排除脊椎病变及损伤。为确保模型完整性，对志愿者进行L3-S1节段CT横断面扫描(64排，层厚1 mm，西门子公司)，获取DICOM格式图像(图1)。

此次研究经上海健康医学院医学伦理委员会批准(医院伦理批件号：2022-SWJWLC-03-310109198608301517，审批时间：2022-08-27)，志愿者知情同意并自愿参与此次研究。
台式计算机辅助建模平台(电脑配置：处理器：Intel(R) Core(TM) i7-10700F CPU@2.90 GHz 2.90 GHz；RAM 48.0 G；Windows 10)；Mimics 20.0建模软件进行原始模型的提取与转化；3-matic 20.0建模软件进行部分附属结构建立、网格生成及模型优化；HyperMesh 14.0进行零件结构的组装、网格生成、韧带结构建立，Ansys Workbench 2020 软件添加材料属性、边界条件、坐标及载荷设定及进行生物力学分析。
1.4 方法
1.4.1 有限元模型的建立使用Mimics进行图像预处理，选择默认骨骼灰度阈值初步分割椎体。再利用编辑蒙版(Edit Masks)命令对图层进行逐层分离，分别对L4/L5椎骨与毗邻软组织进行准确分割。选择恰当的分离工具圆形(Circle)、拉索(Lasso)及方形(Square)等，提取过程中不同的解剖部位选用不同的分离工具，以创建脊椎和所有骨骼特征的实体面罩。应用形态学操作和区域增长的方式，生成L4-L5节段椎体3D模型。将得到的三维几何边界数据以stl格式文件导出。

由于椎间盘的灰度值与周围软组织极其相近，不宜采用区分灰度阈值的方法建立其几何模型。将stl格式文件导入3-matic软件对椎间盘进行初步建模。利用抽壳(Hollow)功能区分骨松质及骨皮质，椎体定义为被0.5 mm厚皮质骨外壳包围的松质骨核心。划分椎体四面体网格，生成igs格式文件。单元的划分类型直接影响计算精度，在相同情况下，六面体网格密度较四面体网格密度小，且更易收敛。此文有限元模型中，对特别关注的韧带及椎间盘采用六面体单元。将igs格式文件导入Hypermesh划分六面体纤维环基质及髓核，通过rods下的linear soild功能建立纤维环纤维，其中髓核占整个椎间盘面积的50%。腰椎模型韧带为厚度1 mm的前纵韧带、后纵韧带、横韧带、黄韧带、棘间韧带、棘上韧带和关节囊韧带(图2)[19-21]。

将模型的cdb文件导入Ansys进行材料赋予、接触设置及边界属性加载等，并进行计算及结果查看。所有材料均简化为各向同性材料，用弹性模量和泊松比2个参数描述。各个部分的材料参数由既往文献获得[22-23](表1)。关节囊韧带与L4关节突为不分离接触，椎骨与椎间盘、椎骨与其他韧带为绑定接触。

1.4.2 模型有效性验证约束L5椎体下表面所有自由度，对特定节点施加相应力矩，模拟腰椎的前屈、后伸、侧弯和扭转4种运动状态。对健康腰椎模型施加10 N·m力矩，以确保与参考文献载荷条件一致，通过活动度验证腰椎模型的有效性。其中，侧弯和扭转分为左右两侧，数据取平均值。由验证结果得，前屈、后伸、侧弯和扭转下的活动度数据分别为6.08°，3.55°，4.19°及2.71°(图3)。通过将有限元分析结果与文献中的实验数据进行定量比较[16，24-25]，侧弯数据略小于文献，其他数据均在文献数据范围内。

在椎骨L4上表面施加500，1 000，1 500，2 000 N的均匀轴向载荷，得到L4的轴向位移分别为0.29，0.58，0.86，1.15 mm。与已通过验证的腰椎模型进行比较[25-27]，验证结果显示，此模型数据均在合理的范围内(图4)。所建模型可以用于生物力学研究。

1.4.3 测试条件限制下椎体所有自由度，在L4椎体特定位置施加5 N·m力矩，模拟前屈、后伸、侧弯及扭转4种运动状态。蒋强等[28]通过改变弹性模量来模拟韧带受损程度，此文参照这种方式来模拟韧带受损，如将弹性模量降低10%即该韧带损伤10%，此文模拟韧带健康及受损10%-50%的情况。改变5次L4-L5之间的前纵韧带弹性模量，分别执行4种运动状态测试，并记录椎间盘压力、活动度及韧带应力结果以进行比较。对后纵韧带、黄韧带、棘上韧带、棘间韧带、关节囊韧带及横突韧带重复同样的过程。每模拟一次韧带损伤状况就进行4种运动状态的测试并记录数据。在任何给定的测试中，只有一种韧带的一个受损程度，一共完成了1 368个数据的记录。
1.5 主要观察指标椎骨的活动度以及韧带的应力，椎骨活动度值通过运动前后空间坐标计算获得，应力指标由软件结果观察所得。

讨论

3.1 单韧带损伤下的腰椎运动分析此次研究证实了作者的初步假设，渐进性韧带损伤会导致腰椎活动度增加。先前的研究已经调查了韧带断裂在单个L4-L5脊柱单元上的作用，并且类似地报道了随着周围稳定结构的丧失，运动增加[29-32]。牛陆等[33]通过增大韧带参数模拟韧带硬化，其利用三维有限元分析探讨后纵韧带硬化后对颈椎生物力学的影响，分析探讨其可能引起的颈部稳定性改变；发现后纵韧带硬化后，前屈及旋转运动使后纵韧带承载更多的力，可能加重已有病变或引起继发损伤；仰伸运动会使最大受力部位转移到小关节，导致继发性小关节损伤的可能性增加。蒋强等[28]降低弹性模量参数模拟退变腰椎，利用退变腰椎有限元模型评估分析全内镜下不同分区精准椎板开窗减压手术对腰椎活动度及应力分布的影响。
此文图5结果表明，渐进式韧带损伤导致所有运动方向的活动度增加。屈曲的最大变化发生在关节囊韧带损伤后，而关节囊韧带和前纵韧带的损伤导致了后伸时活动度的明显增加，这表明它们在矢状面运动(屈曲和伸展)中各自的重要性。屈曲和伸展提供了更多的运动，并且因为它与脊柱健康有关，是许多临床检查的焦点。先前BERMEL等[14]通过三维模型模拟研究发现前纵韧带及关节囊韧带在腰椎运动中的重要作用：关节囊韧带在前屈、后伸和侧弯中的主要作用和前纵韧带在后伸运动中的重要贡献，关节囊韧带的横断不仅导致活动度增大，而且使得弯曲刚度降低、旋转轴发生改变，关节囊韧带功能下降使得其大幅度拉伸以保持腰椎稳定。近来WIDMER等[34]针对腰椎韧带在腰椎中的作用进行了体外实验，通过位移控制的逐步复位研究，对22具人类尸体的腰椎节段进行生物力学测试，实验得腰椎间盘在侧弯运动中承担主要载荷(67%)，横突韧带、棘上韧带和棘间韧带在所有载荷方向上的贡献都很小(分别小于2%和小于6%)，前纵韧带对稳定后伸运动贡献最大，关节囊韧带为主要运动载荷承担者，关节囊韧带对扭转运动起到稳定作用，与此次实验结果有良好的一致性。此外，该研究还发现除椎间盘外，黄韧带和后纵韧带主要在屈曲中起作用(分别为22%和16%)。
此次研究表明后伸对韧带损伤敏感，因此通过此运动更容易发现韧带细小活动度变化。韧带损伤后，侧弯只发生很小的活动度变化，因此侧弯受韧带损伤的影响最小，则冠状面运动(侧弯)对韧带病理变化不敏感。ELLINGSON等[35]有相似的研究结果，除此之外他们的研究提出通常韧带退化和椎间盘病理变性伴随发生。ELLINGSON 等[35]发现后伸运动对韧带损伤敏感，此外侧弯运动对椎间盘病理变化敏感，更容易检测到与退化相关的细微变化，而不会被周围稳定结构的相关变化所掩盖。在韧带损伤的治疗措施中，诊断和康复最为重要。此文研究单韧带损伤后腰椎运动，得到各韧带损伤对腰椎稳定性影响的直观数据，得到每种韧带在运动中的作用，这将为医生在临床诊断和康复治疗中起到重要参考意义。
3.2 单韧带损伤下的生物力学分析韧带的损伤通常会导致该节段剩余韧带的最大应力增加，韧带损伤扰乱了腰椎内部载荷分布。为了保持腰椎稳定，单一韧带受损时其他韧带协同补偿。图6结果得出单韧带损伤后剩余韧带应力发生不同程度的增加，其中关节囊韧带在单韧带受损时总是给出显著力学补偿，表明关节囊韧带受损对周围韧带影响最大，对于维持腰椎运动稳定性极为重要。BERMEL等[14]和NASERKHAKI等[36]就此做出验证，其研究发现横断任何韧带都导致关节囊韧带给出最大力学补偿，随着关节囊韧带的移除，弯曲刚度降低，螺旋轴模式发生改变。此外，在关节囊韧带机械能力下降的情况下，关节囊韧带拉伸大幅增加。另外，横断关节囊韧带并没有影响前屈运动，可能是因为棘间韧带、棘上韧带和横突韧带给出力学补偿。
此外，此次研究结果显示关节囊韧带受损后棘间韧带、棘上韧带和横突韧带在扭转期间最大补偿了力学改变。整体来看棘间韧带和黄韧带应力变化最不明显。后纵韧带受损后的扭转状态下棘上韧带应力发生突变，后纵韧带损伤程度为40%-50%，棘上韧带应力增长百分比从-0.31%上升到20.32%。总之，以上情况都将导致关节囊韧带、棘上韧带、棘间韧带、棘上韧带和横突韧带在扭转期间被过度使用，因此，韧带受损患者应格外注意避免扭转运动。当其余单韧带受损时，剩余韧带同样也会在其他运动状态下给予应力补偿。横突韧带受损后侧弯期间有且仅有棘间韧带给出明显应力补偿，棘间韧带应力增大30.72%。GREGORY等[37]发现，横断韧带、剩余韧带最大应力增加，所有生物力学变化均发生在同一水平，而对相邻节段影响较小。同时，在L3-L4腰椎节段后纵韧带横断对周围韧带应力的影响最小。
韧带损伤后，剩余韧带应力增加，但由图7的结果可知椎间盘压力几乎不变。GREGORY等[37]发现韧带的横断增加其他韧带的应变、所有其他韧带的负荷以及椎间盘上的力矩，但不会显著影响椎间盘的压缩载荷。究其原因可能是韧带横断时，椎间盘力矩受到影响，韧带位置和力臂可以使旋转中心不受引起垂直压缩和水平剪切的载荷垂直分量的影响，从而产生前向平移。虽然韧带损伤不会改变椎间盘压力，但是长期韧带损伤造成腰椎失稳会引起椎间盘损伤[38-39]。韧带损伤对椎间盘应力无明显改变。尽管整个周围组织的剩余负荷分布发生了变化，但椎间盘的压力保持相当稳定。
大多数研究采取直接横断韧带的方法，此文认为渐进式损伤更符合韧带退化或受损现象。而关于渐进式损伤的研究比较少，因而此次研究在韧带相关疾病的检查、预防及康复等方面具有重要参考价值。正如所有的研究一样，几何结构的简化可能导致理想化的结果。此外，使用了简化的材料特性，在未来的研究中将包括更合适的材料特性。
3.3 结论渐进性韧带损伤将影响腰椎稳定性，改变韧带力学分布，剩余韧带应力增大以进行协同补偿。其中，关键性的补偿机制如下：关节囊韧带给予每个韧带损伤时扭转运动下的力学补偿；剩余韧带给予关节囊韧带损伤时各个运动方向上的力学补偿，其中棘间韧带和棘上韧带在扭转期间给予最大力学补偿。后纵韧带受损后的扭转状态下由棘上韧带给出最大应力补偿，横突韧带受损后侧弯期间由棘间韧带给出明显应力补偿。棘间韧带和黄韧带损伤对周围韧带影响最小，关节囊韧带损伤对周围韧带影响最大。
关节囊韧带稳定屈曲和伸展，前纵韧带稳定后伸运动。后伸对韧带病理变化敏感，通过后伸运动更容易发现韧带细小病变。韧带损伤对侧弯运动下的剩余韧带应力影响最小，而对扭转运动下的韧带应力影响最大，韧带损伤患者应尽量避免扭转运动，以防过度使用韧带。韧带受损对椎间盘压力无明显影响，椎间盘压力保持稳定。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

韧带损伤影响腰椎稳定性的生物力学特征

Biomechanical characteristics of ligament injury affecting lumbar spine stability

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 3

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	李智斐, 杨尹, 陈华龙, 梁钦秋, 钟远鸣, 张翼升. 颈椎前路椎体次全切除减压融合钛笼倾斜角度与术后钛笼下沉的有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1313-1319.
[2]	欧阳北平, 马向阳, 罗春山, 邹小宝, 陆廷盛, 陈啟鸰. 新型钉尾横连在后路寰枢椎内固定中的三维有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1320-1324.
[3]	陈萌萌, 包利, 陈浩, 贾璞, 冯飞, 侍管, 唐海. 新型棘突间撑开融合装置BacFuse修复腰椎退行性病的生物力学特征[J]. 中国组织工程研究, 2024, 28(9): 1325-1329.
[4]	何心愉, 周红海, 姜宏, 马智佳, 苏少亭, 林泽宏, 田君明, 陈龙豪, 刘柏杰. L5/S1椎间盘突出重吸收的腰骶矢状面参数改变[J]. 中国组织工程研究, 2024, 28(9): 1330-1335.
[5]	梁成, 张麟麒, 王冠, 李文, 段可, 李忠, 鲁晓波, 卓乃强. 不同内植物修复单侧不稳定骨盆后环损伤的有限元及生物力学分析[J]. 中国组织工程研究, 2024, 28(9): 1336-1341.
[6]	郭苏童, 冯德宏, 郭宇, 王凌, 丁育健, 刘仪, 钱正瑛, 李明洋. 正常与骨质疏松髋关节模型的建立及有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1342-1346.
[7]	杨骏良, 路坦, 徐彪, 姜亚琼, 王富成. 前交叉韧带部分断裂对膝关节应力影响的三维有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1347-1353.
[8]	李超杰, 居来提•买提肉孜, 艾尔西丁•阿不来提, 郑辉, 涂浒第. 不同屈曲角度时内侧髌股韧带重建的有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1359-1364.
[9]	席金涛, 鲁齐林, 汪洋, 王小娟, 吕鹏, 陈龙, 石震, 谢维, 竺义亮, 李绪贵. 经椎间孔腰椎椎间融合后融合器后移的危险因素分析[J]. 中国组织工程研究, 2024, 28(9): 1394-1398.
[10]	余伟杰, 刘爱峰, 陈继鑫, 郭天赐, 贾易臻, 冯汇川, 杨家麟. 机器学习在腰椎间盘突出症诊治中的优势和应用策略[J]. 中国组织工程研究, 2024, 28(9): 1426-1435.
[11]	翁汭, 林东鑫, 郭海威, 张文胜, 宋雨珂, 林宏衡, 李文超, 叶林强. 椎间盘结构异常类型及相关机械载荷与生物力学因素[J]. 中国组织工程研究, 2024, 28(9): 1436-1442.
[12]	白晨, 杨文骞, 孟志超, 王宇泽. 损伤前交叉韧带修复及促进移植物愈合的策略[J]. 中国组织工程研究, 2024, 28(9): 1457-1463.
[13]	左心玮, 刘港, 白惠中, 徐林, 赵毅, 任敬佩, 胡传宇, 穆晓红. 痉挛型脑性瘫痪儿童腰椎与髋关节发育的关系[J]. 中国组织工程研究, 2024, 28(8): 1247-1252.
[14]	张熙辉, 李峥嵘, 李仕能, 邢增宇, 王蛟. 前交叉韧带重建后Pro-kin平衡系统训练患膝本体感觉和平衡功能的评价[J]. 中国组织工程研究, 2024, 28(8): 1259-1264.
[15]	霞黑达·依拉尔江, 尼加提·吐尔逊, 热依拉·库尔班, 白布加甫·叶力思, 陈欣. 三维有限元分析不同特征骨组织中应力的分布规律[J]. 中国组织工程研究, 2024, 28(8): 1277-1282.