有限元法预测运动护膝在不同运动状态下对膝关节韧带的影响

doi:10.12307/2023.733

中国组织工程研究 ›› 2023, Vol. 27 ›› Issue (36): 5771-5777.doi: 10.12307/2023.733

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

有限元法预测运动护膝在不同运动状态下对膝关节韧带的影响

曾露露，谢红

上海工程技术大学纺织服装学院，上海市 201620

收稿日期:2022-10-17 接受日期:2022-11-25 出版日期:2023-12-28 发布日期:2023-03-24
通讯作者: 谢红，博士，教授，上海工程技术大学纺织服装学院，上海市 201620
作者简介:曾露露，女，1997年生，贵州省绥阳县人，汉族，上海工程技术大学在读硕士，主要从事运动生物力学和功能性运动服装的研究。
基金资助:
国家重点研发计划主动健康和老龄化科技应对重点专项(2018YFC2000900) ，项目负责人：谢红

Finite element method for predicting the effect of sports knee brace on knee ligaments under different sports conditions

Zeng Lulu, Xie Hong

School of Textile Fashion, Shanghai University of Engineering Science, Shanghai 201620, China

Received:2022-10-17 Accepted:2022-11-25 Online:2023-12-28 Published:2023-03-24
Contact: Xie Hong, PhD, Professor, School of Textile Fashion, Shanghai University of Engineering Science, Shanghai 201620, China
About author:Zeng Lulu, Master candidate, School of Textile Fashion, Shanghai University of Engineering Science, Shanghai 201620, China
Supported by:
the National Key Special Research and Development Plan, Active Health and Aging Science and Technology Response Project of China, No. 2018YFC2000900 (to XH)

摘要/Abstract

摘要：

文题释义：
服装压力：无论在静止还是运动状态下，人体表面与织物之间均存在接触关系，穿着服装时产生的正压力导致肌肤变形，刺激肌肤深处的压觉点，从而使人体感受到服装压力。服装压力是探究服装舒适性的重要组成部分，也是人体-织物力学模型研究的重要指标，常用于压缩服、医疗袜、功能性材料及智能可穿戴监测服装等研究领域。该研究采用弹性套筒式运动护膝进行实验，其与人体表面形成紧密接触，因此以服装压力实验验证膝关节-运动护膝模型的有效性。
Von Mises应力：是一种基于剪切应变能、考虑了第一、二、三主应力的等效应力，其遵循材料力学的第四强度理论，即畸变能密度理论，可用于疲劳、破坏等评价，判断材料的屈服性能。在相同工况下，Von Mises应力值越大，其计算云图分布越不均匀，应力集中会加大损伤的概率。该研究的目的是借助有限元法，探讨在不同运动状态下，不同材料的运动护膝对膝关节韧带的影响，因此以Von Mises等效应力为主要观察指标。

背景：膝关节韧带是维持膝关节运动稳定性的重要组成部分，在人体运动过程中极易受到损伤，运动护膝常用于预防膝关节运动损伤，但其防护性能未能明确。
目的：建立健康成年人的膝关节有限元模型及运动护膝模型，运用有限元法预测运动护膝在不同运动状态下对膝关节韧带的影响，以Von Mises等效应力为观察指标，探讨不同材料运动护膝的防护性能。
方法：以1名男性健康志愿者的CT影像为数据来源，利用Mimics、Solidworks、Abaqus等软件获得膝关节有限元模型及运动护膝模型；以股骨内外髁中点为参考点施加后向134 N集中力，模拟临床前抽屉实验，获取股骨相对位移及主要韧带的生物力学响应，验证膝关节模型的有效性；给运动护膝施加位移载荷，模拟膝关节直立位穿着护膝的状态，并与服装压力测试结果进行对比，验证膝关节-护膝模型的有效性；分别在无护膝和两种不同材料护膝作用下模拟膝关节屈曲0°，30°，60°运动，分析不同载荷下前交叉韧带、后交叉韧带、内侧副韧带及外侧副韧带的应力应变情况。
结果与结论：①在0°，30°，60°屈曲角度下，并且施加320 N垂直压缩力及134 N股骨后向推力后，相较于裸膝，穿着两种不同材料的运动护膝使得前交叉韧带的峰值应力在0°和30°有所降低，但在屈曲60°状态下峰值应力有所上升；后交叉韧带在屈曲0°时的峰值应力有所上升，在30°和60°状态下的峰值应力有所降低；内侧副韧带和外侧副韧带在0°，30°，60°的峰值应力均有所降低。②结论：不同材料的运动护膝能在一定程度上对膝关节韧带起到防护作用，但在特定角度下会对前交叉韧带和后交叉韧带产生额外的负荷，该研究结果对膝关节韧带的生物力学研究、运动损伤防护及运动护膝的设计与生产等具有积极意义。
https://orcid.org/0000-0003-0169-2806 (曾露露)

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

关键词: 运动护膝, 韧带, 防护性, 有限元模拟, 生物力学

Abstract: BACKGROUND: Knee ligaments are an important part of maintaining the stability of knee joint movement, which is very vulnerable to injury during exercise, and sports knee braces are often used to prevent knee injuries, but their protective performance is not clear.
OBJECTIVE: A finite element model of the knee joint of a healthy adult and a sports knee brace model were established. The finite element method was used to predict the effect of the sports knee brace on the knee joint ligament under different movement states so as to observe the protective performance of different knee braces using Von Mises equivalent stress as the observation index.
METHODS: Using CT medical images of one male health volunteer as the data source, Mimics, Solidworks, and Abaqus software were used to obtain the finite element model of the knee joint and the sports knee brace model. The median epicenter of the femur and the medial condyle was set as the reference point to apply backward 134 N concentration force to the point, to simulate the preclinical drawer experiment, to obtain the biomechanical response of the femur relative displacement and major ligaments, and to verify the effectiveness of the knee model. The displacement load was applied to the sports knee brace, to simulate the state of wearing the knee brace in an upright position and to compare with the clothing stress test results to verify the effectiveness of the knee joint-knee brace model. The 0°, 30° and 60° movements of knee joint flexion were simulated without knee braces and knee braces of two different materials. The stress strains of the anterior cruciate ligament, posterior cruciate ligament, medial collateral ligament and lateral collateral ligament were analyzed under different loads.
RESULTS AND CONCLUSION: (1) The peak stress of the four ligaments was changed at 0°, 30°, and 60° angles, and after applying 320 N vertical compressive force and 134 N femur backward thrust, compared with bare knee, the peak stress of the anterior cruciate ligament was reduced at 0° and 30°, but the peak stress in the 60° flexion increased. The peak stress of the posterior cruciate ligament increased at 0° flexion and decreased in the 30° and 60° flexion. Peak stress of the medial collateral ligament and lateral collateral ligament was reduced at 0°, 30°, and 60° angles. (2) Conclusion: Different materials of sports knee braces can protect the knee ligaments to a certain extent, but at specific angles will cause additional loads on the anterior cruciate ligament and posterior cruciate ligament, and the results of this study have the positive significance for the biomechanical research of knee ligaments, sports injury protection and the design and production of sports knee braces.

Key words: sports knee braces, ligament, protective performance, finite element simulation, biomechanics

中图分类号:

曾露露, 谢红. 有限元法预测运动护膝在不同运动状态下对膝关节韧带的影响[J]. 中国组织工程研究, 2023, 27(36): 5771-5777.

Zeng Lulu, Xie Hong. Finite element method for predicting the effect of sports knee brace on knee ligaments under different sports conditions[J]. Chinese Journal of Tissue Engineering Research, 2023, 27(36): 5771-5777.

图/表（结果） 9

2.1 膝关节模型验证膝关节在受到134 N后向推力时，股骨在前后方向上的位移为3.92 mm，与前人研究结果相似，其中GABRIEL等[20]为尸体实验，SONG等[21]、万超等[22]和何川等[15]均是有限元模型；前交叉韧带、后交叉韧带、内侧副韧带及外侧副韧带的峰值应力分别为19.95，8.84，4.0及1.15 MPa，从韧带受力分析可知，前交叉韧带的受力载荷大于其他主要韧带，各条韧带的峰值应力均出现在胫骨、股骨附着区，这与临床诊断前交叉韧带是防止胫骨前移的主要作用，且膝关节韧带损伤主要集中在骨骼止点部位的临床现象相一致[21，23]。
2.2 膝关节-护膝模型验证在模拟膝关节分别穿上两种护膝后，提取膝关节测量点的接触压力，并与服装压力测试数据进行对比。首先对模型测量点进行设定，与上述服装压测试操作一致，设置髌骨中心为参考点，依次选定9个测量点，每个测量点与气囊式压力传感器面积大小一致，最终计算其接触压力平均值。对9个压力测量点的模拟值与实验值进行相对误差分析，两种材料9个测量点的相对误差均小于20%，与多数服装压力研究中的误差范围一致[24-25]，因此认为该膝关节-护膝模型具有可靠性，各测量点对比见图2。

2.3 膝关节运动预测模拟结果
2.3.1 不同护膝材料对前交叉韧带的应力分布影响在0°，30°，60°屈曲角度下，并且施加320 N垂直压缩力及134 N股骨后向推力后，膝关节无或穿戴护膝状态下，前交叉韧带峰值应力见表5，应力均集中在韧带前端，是较易产生损伤的部位。在穿着两种不同材料护膝运动后，前交叉韧带的峰值应力在0°和30°状态下均有所降低，材料2的效果较好，使其分别降低了7.9 MPa和1.09 MPa，材料1使其分别降低了6.38 MPa和0.26 MPa，但在屈曲60°状态下，韧带的峰值应力相较于无护膝时均有所上升，其中材料1使其上升10.82 MPa，材料2使其上升15.73 MPa，前交叉韧带在0°，30°，60°状态下的应力情况见图3。

2.3.2 不同护膝材料对后交叉韧带的应力分布影响在0°，30°，60°屈曲角度下，并且施加320 N垂直压缩力及134 N股骨后向推力后，膝关节无或穿戴护膝状态下，后交叉韧带的峰值应力见表6，峰值应力在屈曲过程中从股骨、胫骨连接处转移到后前侧，应力呈下降趋势。在穿着两种不同材料护膝运动后，后交叉韧带在屈曲0° 状态下均有所上升，材料1、材料2分别使其上升1.62 MPa和2.72 MPa；在屈曲30°与屈曲60°状态下，后交叉韧带的应力均有所降低，材料1使其降低0.47 MPa和0.45 MPa，材料2使其降低1.01 MPa和0.36 MPa，后交叉韧带在0°，30°，60°状态下的应力情况见图4。

2.3.3 不同护膝材料对内侧副韧带的应力分布影响在0°，30°，60°屈曲角度下，并且施加320 N垂直压缩力及134 N股骨后向推力后，膝关节无或穿戴护膝状态下，内侧副韧带的峰值应力见表7，应力均集中在股骨、胫骨附着区。在穿着两种不同材料护膝后，内侧副韧带在各个角度的峰值应力均有所降低，其中材料1的运动护膝使其在0°，30°和60°屈曲角度下分别减少了0.63，0.28和 0.94 MPa，材料2的运动护膝使其分别减少了0.76，0.01和0.64 MPa，内侧副韧带在0°，30°，60°状态下的应力情况见图5。

2.3.4 不同护膝材料对外侧副韧带的应力分布影响在0°，30°，60°屈曲角度下，并且施加320 N垂直压缩力及134 N股骨后向推力后，膝关节无或穿戴护膝状态下，外侧副韧带的峰值应力见表8，应力均集中在股骨、腓骨附着区。在穿着两种不同材料护膝后，外侧副韧带在各个角度的峰值应力均有所降低，其中材料1的运动护膝使其在0°，30°和60°屈曲角度下分别减少了0.19，0.09和 0.82 MPa，材料2的运动护膝使其分别减少了0.01，1.43和0.78 MPa，外侧副韧带在0°，30°，60°状态下的应力情况见图6。

参考文献

[1] MAJEWSKI M, SUSANNE H, KLAUS S. Epidemiology of athletic knee injuries: A 10-year study. Knee. 2006;13(3):184-188.
[2] POSCH M, SCHRANZ A, LENER M, et al. In recreational alpine skiing, the ACL is predominantly injured in all knee injuries needing hospitalisation. Knee Surg Sports Traumatol Arthrosc. 2021;29(6):1790-1796.
[3] BROPHY RH, WOJTYS EM, MACK CD, et al. Factors associated with the Mechanism of ACL tears in the National Football League: A Video-Based Analysis. Orthop J Sports Med. 2021;9(11):23259671211053301.
[4] MANIAR N, VERHAGEN E, BRYANT AL, et al. Trends in Australian knee injury rates: An epidemiological analysis of 228,344 knee injuries over 20 years. Lancet Reg Health West Pac. 2022;21:100409.
[5] 陈凌娴,李俊,王敏.护膝防护性能及其功能设计研究进展[J].毛纺科技,2022,50(3):117-123.
[6] MOON J, KIM H, LEE J, et al. Effect of wearing a knee brace or sleeve on the knee joint and anterior cruciate ligament force during drop jumps: A clinical intervention study. Knee. 2018;25(6):1009-1015.
[7] 荆晓伟, 母应秀. 硅网支具对前交叉韧带损伤及膝关节稳定性的影响[J]. 中国组织工程研究,2020,24(4):505-510.
[8] ALTHOMALI OW. Influence of using knee sleeve and lateral wedge insole on knee loading among healthy individuals during stair negotiation. Gait Posture. 2022;92:103-109.
[9] VALLDECABRES R, DE BENITO AM, LITTLER G, et al. An exploration of the effect of proprioceptive knee bracing on biomechanics during a badminton lunge to the net, and the implications to injury mechanisms. Peer J. 2018;6:e6033.
[10] JIANG T, TIAN S, FAN X, et al. Kinematics and kinetics of lower-extremity joints in parachuting landing with backpack and knee brace. Med Eng Phys. 2020;86:1-7.
[11] BOTTONI G, HERTEN A, KOFLER P, et al. The effect of knee brace and knee sleeve on the proprioception of the knee in young non-professional healthy sportsmen. Knee. 2013;20(6):490-492.
[12] KRAMI M, QIAN Z, ZOU Z, et al. Subject-specific finite element modelling of the human foot complex during walking: sensitivity analysis of material properties, boundary and loading conditions. Biomech Model Mechanobiol. 2018;17(2):559-576.
[13] PENA E, CALVO B, MARTINEZ 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.
[14] 万超杰.人体大臂软组织活体力学特性分析[D].哈尔滨:哈尔滨工业大学,2017.
[15] 何川,李彦林,张振光,等.不同屈曲状态下膝关节韧带生物力学的有限元分析[J].中国运动医学杂志,2015,34(7):662-669.
[16] 孙亚博. 基于有限元模拟的针织服装穿着压力分析[D].天津:天津工业大学,2021.
[17] GROOD ES, SUNTAY WJ. A joint coordinate system for the clinical description of three-dimensional motions: application to the knee. J Biomech Eng. 1983;105(2):136-144.
[18] SKINNER H. Current diagnosis and treatment in orthopedics.USA: McGraw-Hill Professional Publishing. 2006.
[19] 黎蒙. 女性羽毛球运动员头顶球单腿落地动作的膝关节生物力学有限元分析[D].上海:上海体育学院,2020.
[20] GABRIEL MT, WONG EK, WOO SLY, et al. Distribution of in situ forces in the anterior cruciate ligament in response to rotatory loads. J Orthop Res. 2004;22(1):85-89.
[21] SONG YH, 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.
[22] 万超,郝智秀,温诗铸.前交叉韧带力学特性差异对膝关节有限元仿真结果的影响[J].医用生物力学,2012,27(4):375-380.
[23] SHERMAN MF, LIEBER L, BONAMO JR, et al. The long-term followup of primary anterior cruciate ligament repair. Defining a rationale for augmentation. Am J Sports Med. 1991;19(3):243-255.
[24] YU A, YICK KL, NG SP, et al. Numerical simulation of pressure therapy glove by using Finite Element Method. Burns. 2016;42(1):141-151.
[25] 于洁.基于人体臂部模型的服装压测试与数值模拟[D].天津:天津工业大学,2019.
[26] CHEW KT, LEW HL, DATE E, et al. Current evidence and clinical applications of therapeutic knee braces. Am J Phys Med Rehabil. 2007; 86(8):678-686.
[27] COUDEYRE E, NGUYEN C, CHABAUD A, et al. A decision-making tool to prescribe knee orthoses in daily practice for patients with osteoarthritis. Ann Phys Rehabil Med. 2018;61(2):92-98.
[28] 庞欣.高弹护膝纺织品的开发与压力性能分析[D].杭州:浙江理工大学,2018.
[29] 沈妙音.抗菌热湿舒适性服用机织物的研究与开发[D].杭州:浙江理工大学,2021.
[30] 匡新谋,苏敏茹,李浩,等.预处理工艺对回收涤纶/锦纶/氨纶混合纤维改性聚丙烯力学性能影响[J].工程塑料应用,2020,48(9): 78-82.
[31] WIGGINS AJ, GRANDHI RK, SCHNEIDER DK, et al. Risk of Secondary Injury in Younger Athletes After Anterior Cruciate Ligament Reconstruction: A Systematic Review and Meta-analysis. Am J Sports Med. 2016;44(7):1861-1876.
[32] ANDRIACCHI TP, MÜNDERMANN A, SMITH RL, et al. A framework for the in vivo pathomechanics of osteoarthritis at the knee. Ann Biomed Eng. 2004;32(3):447-457.
[33] FOCKE A, STEINGREBE H, MÖHLER F, et al. Effect of Different Knee Braces in ACL-Deficient Patients. Front Bioeng Biotechnol. 2020;8:964.
[34] HOLDEN MA, CALLAGHAN M, FELSON D, et al. Clinical and cost-effectiveness of bracing in symptomatic knee osteoarthritis management: protocol for a multicenter, primary care, randomized, parallel-group, superiority trial. BMJ Open. 2021;11(3):e048196.
[35] PIERRAT B, OULLION R, MOLIMARD J,et al. Characterisation of in-vivo mechanical action of knee braces regarding their anti-drawer effect.Knee. 2015;22(2):80-87.
[36] MOHD SHARIF NA, USMAN J, WAN SAFWANI WKZ, et al. Effects of simple knee sleeves on pain and knee adduction moment in early unilateral knee osteoarthritis. Proc Inst Mech Eng H. 2019;233(11): 1132-1140.
[37] SINCLAIR J, TAYLOR PJ. Effects of a Prophylactic Knee Sleeve on Anterior Cruciate Ligament Loading During Sport-Specific Movements. J Sport Rehabil. 2019;28(1):1-7.
[38] BALTACI G, AKTAS G, CAMCI E, et al. The effect of prophylactic knee bracing on performance: balance, proprioception, coordination, and muscular power.Knee Surg Sports Traumatol Arthrosc. 2011;19(10): 1722-1728.
[39] GIOTIS D, TSIARAS V, RISTANIS S, et al. Knee braces can decrease tibial rotation during pivoting that occurs in high demanding activities. Knee Surg Sports Traumatol Arthrosc. 2011;19:1347-1354.
[40] FLEMING BC, RENSTROM PA, BEYNNON BD, et al. The influence of functional knee bracing on the anterior cruciate ligament strain biomechanics in weightbearing and nonweightbearing knees. Am J Sports Med. 2000;28(6):815-824.

引言

材料方法

1.1 设计有限元法预测运动护膝的防护性能。
1.2 时间及地点实验于2021年11月至2022年8月在上海工程技术大学人体工效与功能服装实验室完成。
1.3 对象选择1名中国男性健康志愿者，年龄24岁，身高175 cm，体质量65 kg；既往无膝关节创伤史，无骨关节炎史。使用CT对研究对象的膝关节直立位进行螺旋式扫描，扫描参数为层厚0.5 mm，层间距1 mm，分辨率512×512像素。CT扫描共计成像707张，将其保存为医学数字成像(digital imaging and communications in medicine，DICOM)格式存储在光盘中。受试者对试验知情同意。
该研究方案的实施符合上海市第一人民医院的相关伦理要求。
1.4 材料
1.4.1 实验设备 HP Z4 G4工作站(CPU：i9-10900X；内存：64 GB；显卡：NVIDIA Quadro P2200)、数字式织物厚度仪(YG141D型)、电子天平(LE104E)、动态热机械分析仪(DMA850)、气囊式接触压力测试仪(AMI3037-10-II)由上海工程技术大学纺织服装学院提供。
1.4.2 软件设备 Mimics医学影像控制软件、3-Matic Research(Materialise公司，比利时)；Geomagic Studio逆向工程软件(Geomagic 公司，美国)；Solidworks三维建模软件(Dassault Systemes，美国)；Hypermesh网格前处理软件(Altair公司，美国)；Abaqus有限元分析软件(ABAQUS公司，美国)。
1.5 方法
1.5.1 三维模型的建立把DICOM格式的膝关节二维图片导入Mimics中，通过区域增长、计算实体等操作构建模型，以STL(Standard Template Library)格式导出，再将其导入逆向工程软件Geomagic Studio中，对部件进行多边形和精确曲面处理，以STEP格式导出，再导入Solidworks软件中进行装配、布尔操作等，处理完后以STEP格式导出膝关节整体模型。将从Solidworks软件中导出的模型再导入Hypermesh软件中进行网格处理，将模型以INP格式导出。模型包括：软组织、股骨、胫骨、腓骨、髌骨、前交叉韧带、后交叉韧带、内侧副韧带、外侧副韧带。在3-Matic Research软件中对膝关节软组织进行抽壳处理后以STL格式导出，再导入Geomagic Studio中进行平滑、轮廓线勾勒处理以STEP格式导出，再将其导入Solidworks中以X_T格式导出，将膝关节模型及护膝模型导入Abaqus 2019中进行整体装配。
1.5.2 运动护膝材料测试研究选择针对健康人体运动时常用的弹性套筒式护膝为主要研究类型，选择氨纶、聚酯纤维及锦纶为原料的运动护膝产品。通过预实验挑选两款力学参数差异的复合面料，分别记作材料1、材料2。利用电子天平称量面积3 cm×3 cm 织物的质量，利用YG141D型数字式织物厚度仪测量材料的织物厚度，通过物理计算，即密度为质量与体积之比，分别求得2种运动护膝材料的密度值；利用 DMA850 动态热机械分析仪来对面料试样的弹性模量进行测试，选择薄膜拉伸夹具，将拉伸速率设置为100 mm/min，根据应力和应变计算其弹性模量，用最小二乘法将数据拟合为应力-应变直线，计算出其相关系数，若大于 0.97，则认为该直线的斜率即为弹性模量。每种试样测量3次取平均值，材料属性见表1。

1.5.3 服装压力测试服装压力是验证膝关节-护膝模型有效性的重要方式。研究采用AMI3037-10-II气囊式接触压力测试系统进行测量，传感器量程为0-35 kPa，精度为±0.1 kPa，
受试者与模型人体的身高、体型相近。在膝关节中，髌骨结构凸出，是膝关节伸膝功能的重要组成部分，能较好反映服装压的变化，且便于膝关节模型测量点的选择，因此以髌骨中心为参考点。首先在皮肤上手动勾勒出髌骨结构范围，确定大致中点位置；其次测量膝关节模型高度及运动护膝的长度，由受试者穿戴两款护膝进行预实验，确定两款护膝覆盖人体的范围大小，经对比膝关节模型及受试者穿戴护膝覆盖范围，选择最佳的测量点。对比后，在髌骨参考点向上量取6 cm为第一个测量点，测量该点所在水平腿围，平均量取其余3个测量点，向下同理，总共9个测量点。测量时受试者静止站立，两款护膝各测量5次取平均值，具体测量值见表2。

1.5.4 材料属性的设置在Abaqus中设置模型的材料属性。运动护膝近似于壳体，其材料属性见表1；为了对膝关节模型进行有效性验证，膝关节内各结构属性均采自以往文献，该模型中的骨骼设为各向同性线弹体材料[12]，材料属性见表3；软组织及韧带定义为超弹体材料[13-14]，选用Neo-hooker模型，Neo-hooke模型函数如下：

软组织及各韧带的材料属性见表4。

1.5.5 网格划分膝关节软组织、骨骼及韧带均在Hypermesh软件中进行网格划分、网格类型选择及网格质量检查，其中骨骼及韧带的网格类型为四面体、C3D10M类型，股骨、胫骨及腓骨的网格大小均为5 mm，髌骨为3 mm，韧带均为1 mm；软组织为四面体、C3D4类型，网格大小为6 mm；护膝模型在Abaqus软件中进行网格划分，网格类型为四面体、S4R类型，网格大小为2 mm，网格质量均良好，膝关节、护膝网格模型如图1所示。

1.5.6 接触定义在该膝关节模型中，将除骨骼、韧带以外的半月板、关节软骨、肌肉及皮肤等组织简化为整体软组织。韧带与骨骼、骨骼与软组织内部设置为绑定接触[15]，护膝内表面与软组织表面设为面与面接触，摩擦系数为0.2，接触方向属性为硬接触[16]。
1.5.7 边界条件、载荷及验证基于GROOD等[17]的研究建立膝关节局部坐标系，以股骨内外髁连线中点为坐标系的原点，X、Y、Z轴定义了膝关节在6个自由度上的屈曲、内外翻、内外旋与平移运动，并将股骨髁中点设置为参考点，耦合股骨表面所有节点于该点上，有助于后续力的加载及数据提取，且均采用显示动力学进行计算。
(1)膝关节模型的有效性验证：前抽屉实验是重要的临床诊断、评价方式之一，用于检查膝关节韧带的松弛度及健康程度[18]，也是有限元模型有效性验证的常用方式，此次研究采用该方法验证膝关节有限元模型的有效性，固定胫腓骨远端，对股骨参考点施加134 N的后向推力，模拟前抽屉实验。
(2)膝关节-护膝模型的有效性验证：人体静止站立时，膝关节所受负荷约为人体质量的0.43倍[19]，此次研究对象体质量65 kg，经计算后以320 N为模型单腿重力载荷，固定胫腓骨远端，对股骨参考点施加320 N垂直向下的压缩力，模拟人体静止站立单腿受力状态，经预计算后，给运动护膝施加160 mm的位移载荷，分别模拟膝关节直立位穿戴两种不同材料运动护膝的状态，提取膝关节表面测量点的接触压力，并与服装压力测试数据进行对比。
(3)膝关节屈曲运动的模拟预测：在人体屈曲运动中，膝关节韧带起着牵引、稳定及力学传递的重要作用，而前抽屉实验可以检验韧带的健康程度，为了更好地突出关节韧带的力学变化，以重力载荷及前抽屉力为该膝关节模型屈曲受力状态，固定胫腓骨远端，以股骨内外髁连线为旋转轴，在股骨参考点施加一定的角速度，模拟膝关节屈曲0°，30°，60°的运动状态，并且在不同屈曲角度下对股骨参考点施加320 N垂直向下的压缩力及134 N股骨后向推力。
1.6 主要观察指标 ①膝关节模型验证结果；②膝关节-护膝模型验证结果；③膝关节运动预测模拟结果。

讨论

护膝按不同性能、材料、结构及用途均有不同的分类[5，26-27]，此次研究选择针对健康人体运动时常用的弹性套筒式护膝为主要研究类型，其可以在不影响人体正常运动的前提下，尽可能减少膝关节的活动幅度和范围，一方面避免不科学的运动方式而引起的关节、肌肉及韧带损伤；另一方面有利于减少肌肉振动、增强肌肉本体感受等。在材料方面，弹性纤维及面料紧度可以很好地延缓应力衰变，增加弹性护膝的稳定性[5]，其中氨纶弹性纤维具有优良的延伸性、较低的弹性模量及高弹回复率等特点，其伸长率达400%-800%，是弹性织物中常用的原料之一[28]；聚酯纤维具有良好的保形性、抗皱性、弹性回复性及较高的强度[29]，锦纶具有良好的耐磨性及弹性，复合织物能加强运动护膝的弹性、延伸性、耐磨性及防护性等[30]，因而此次试验选择氨纶、聚酯纤维及锦纶为原料的运动护膝产品。
此次研究的目的在于利用有限元模拟方法，预测健康人体在裸膝及穿着两种不同材料护膝的状态下进行不同运动时膝关节韧带的生物力学响应，以Von Mises等效应力为观察指标，探讨了运动护膝对膝关节韧带损伤的影响。研究结果表明，运动护膝对膝关节韧带有一定的防护作用，但在某些特定角度下会加大前交叉韧带和后交叉韧带的负荷。研究结果显示，在0°，30°，60°屈曲角度下，并且施加320 N垂直压缩力及134 N股骨后向推力后，相较于裸膝状态，膝关节韧带的峰值应力均有所变化，前交叉韧带的峰值应力在0°及30°状态下均有所降低，此时材料2的防护效果较好，但在60°状态时，穿着两种材料护膝后前交叉韧带的峰值应力均有所上升，这增加了韧带受损的风险；后交叉韧带的峰值应力在0°状态下均有所上升，在30°及60°状态下均有所降低；内侧副韧带及外侧副韧带在各个角度下的峰值应力均有所降低，材料1与材料2的防护效果差异不大。
膝关节韧带是重要的稳定结构，在人体运动过程中极易受到损伤，如韧带撕裂、拉伤等不仅会影响膝关节运动及神经肌肉变化，还会降低膝关节本体感受能力，这进一步加大了继发性损伤及慢性疾病如骨关节炎的风险[4，31-32]。膝关节护具常被用于预防运动损伤、矫正畸形及损伤康复，从而改善膝关节运动功能及延缓疾病进展，其类型众多，功能各异，美国矫形外科学会从传统上定义了4种膝关节护具的分类[26]，分别为预防型、功能型、康复型及髌股型；而根据材料、结构及刚度等不同，又可分为套筒式、铰链式及卸荷式护具[27]。这些膝关节护具大多用于膝骨关节炎治疗、前交叉韧带重建治疗及术后康复等用途，且部分研究已证实其功能及辅助作用[33-36]。其中PIERRAT等[35]选择了健康人体及膝关节损伤患者，通过关节测量仪评估了套筒式护膝及铰链式护膝的静态抗抽屉能力，结果表明，无论是健康人体还是膝关节受损患者，两种类型的护膝均会有效减小前交叉韧带受力，铰链式护膝的效果更突出，可能是由于刚性支撑能吸收更多力，从而达到更好的效果，在该研究中膝关节是小屈曲状态下进行实验，因此与此研究中健康人体在0°及30°角度下，穿戴套筒式护膝后前交叉韧带受力下降的趋势一致。在日常生活中，正常健康人体使用护膝的目的在于预防运动损伤，降低运动过程中所受损伤的概率，针对竞技运动员，其有提高肌力及改善成绩的作用，部分研究表明，运动护膝可以有效减小膝关节运动范围、改善肌肉本体感受、减小前交叉韧带应力及胫骨旋转运动[37-40]，多数研究利用三维运动测量系统、角度传感器等工具，得出人体外部的动力学及运动学数据，但较难获得人体内部组织的变化。此次研究选择运动人群常用的弹性套筒式运动护膝为主要研究类型，其由不同比例的锦纶、聚酯纤维及氨纶原料制成，利用有限元法，预测了不同材料的运动护膝对健康人体膝关节韧带的防护效果。结果表明，弹性套筒式运动护膝对健康膝关节韧带的防护效果有限，相同原料但比例不同的护膝在性能上并没有太大的差异，对膝关节韧带并不能达到全面防护的预期设想，应基于不同运动状态下，对护膝的结构、材料及功能进行设计与研究，如运动护膝的分区设计、刚性或半刚性设计、多功能面料选择等，以达到对膝关节韧带针对性预防及防护的目的。
此次研究基于有限元法，探讨了弹性套筒式运动护膝对正常膝关节韧带的影响，对正常膝及病理膝的研究均有借鉴意义。但研究也具有一定的局限性：首先，在构建模型时，将半月板及关节软骨等组织简化为包覆骨骼及韧带的整体软组织，忽略了部分组织的力学性能；其次，膝关节在屈曲过程中受力复杂，为了突出韧带力学变化，将其受力简化为重力载荷及前抽屉力，模型并不能完全体现膝关节完整的受力状态；最后，研究只针对膝关节屈曲0°，30°及60°下的运动状态，且运动护膝只涉及弹性套筒式护膝，类型较单一，未来研究应更加注重模型的完整性及精确性。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

有限元法预测运动护膝在不同运动状态下对膝关节韧带的影响

Finite element method for predicting the effect of sports knee brace on knee ligaments under different sports conditions

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 9

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	钟毅征, 黄培镇, 蔡群斌, 郑利钦, 何兴鹏, 董航. 影响骨小梁微有限元模型最大应力的骨微结构指标[J]. 中国组织工程研究, 2023, 27(9): 1313-1318.
[2]	吴韬光, 聂少波, 陈华, 朱正国, 祁麟, 唐佩福. 新型多维交叉锁定钢板固定股骨转子下骨不连的生物力学特征[J]. 中国组织工程研究, 2023, 27(9): 1330-1334.
[3]	彭志鑫, 闫文刚, 王坤, 张振江. 3D打印前臂外固定支具的有限元分析与结构优化设计[J]. 中国组织工程研究, 2023, 27(9): 1340-1345.
[4]	吴天亮, 陶秀霞, 徐宏光. 三维有限元法分析不同骨密度对单纯斜外侧腰椎椎间融合后融合器下沉的影响[J]. 中国组织工程研究, 2023, 27(9): 1352-1358.
[5]	刘金玉, 张晗硕, 崔洪鹏, 潘灵之, 赵博然, 李菲, 丁宇. 脊柱微创治疗脊髓型颈椎病的有限元生物力学分析及精准化运动康复方案[J]. 中国组织工程研究, 2023, 27(9): 1359-1364.
[6]	温星花, 丁焕文, 成凯, 闫晓楠, 彭元昊, 王宇宁, 刘康, 张挥武. 比格犬股骨大段骨缺损髓内钉固定方案设计的三维有限元建模分析[J]. 中国组织工程研究, 2023, 27(9): 1371-1376.
[7]	贺垠皓, 李晓声, 陈宏文, 陈铁柱. 3D打印多孔钽金属治疗发育性髋关节发育不良：现状及应用前景[J]. 中国组织工程研究, 2023, 27(9): 1455-1461.
[8]	刘佳辛, 贾鹏, 门玉涛, 刘璐, 王烨明, 叶金铎. 基于三周期极小曲面骨小梁结构的设计及优化[J]. 中国组织工程研究, 2023, 27(7): 992-997.
[9]	李惠琴, 王春娟, 王禹, 王伟峰, 陈定根, 李娜. 无托槽隐形矫治器治疗下颌不同形态扭转牙的三维有限元分析[J]. 中国组织工程研究, 2023, 27(7): 1050-1054.
[10]	刘广銮, 郭宗磊, 戈进, 黄栋, 王业华. 前交叉韧带损伤伴发内侧半月板后根撕裂：解剖学的危险因素分析[J]. 中国组织工程研究, 2023, 27(5): 663-668.
[11]	熊波涵, 余洋, 卢晓君, 王旭, 杨腾云, 张瑶璋, 廖欣宇, 周晓翔, 何璐, 李彦林. 促进前交叉韧带重建腱骨愈合的研究与进展[J]. 中国组织工程研究, 2023, 27(5): 779-786.
[12]	李雅萍, 刘红, 高真, 陈晓琳, 黄武杰, 江征. 太极拳练习伴膝关节疼痛者下肢生物力学表现的三维运动分析[J]. 中国组织工程研究, 2023, 27(4): 520-526.
[13]	李世浩, 李奇, 李震, 张媛媛, 刘苗苗, 欧阳懿, 徐卫国. 前交叉韧带损伤与重建术后患者的足底压力及步态分析[J]. 中国组织工程研究, 2023, 27(4): 626-631.
[14]	陈鹏, 王玲, 董诗雨, 丁月, 贾绍辉, 寇现娟, 郑成. 全身振动训练对前交叉韧带重建后影响的Meta分析[J]. 中国组织工程研究, 2023, 27(36): 5875-5883.
[15]	文鹏飞, 李亚宁, 路玉峰, 郝林杰, 王亚康, 马涛, 宋伟, 张育民. 腰椎-骨盆-髋关节有限元模型建立及生物力学分析[J]. 中国组织工程研究, 2023, 27(36): 5741-5746.