膝关节屈伸运动形态的有限元分析

doi:10.3969/j.issn.2095-4344.2017.03.013

中国组织工程研究 ›› 2017, Vol. 21 ›› Issue (3): 396-400.doi: 10.3969/j.issn.2095-4344.2017.03.013

• 数字化骨科 digital orthopedics • 上一篇下一篇

膝关节屈伸运动形态的有限元分析

张立超¹，张立敏²，吕永明¹，王智慧¹，杨阳¹，徐飞¹，戴海峰¹，李嘉¹，曹向宇¹，武丽珠³

¹承德医学院附属医院骨科，河北省承德市 067000；²解放军第266医院，河北省承德市 067000；³承德市中心医院急诊科，河北省承德市 067000

修回日期:2016-11-17 出版日期:2017-01-28 发布日期:2017-03-14
作者简介:张立超，男，1981年生，河北省宽城县人，满族，2005年承德医学院毕业，硕士，主治医师，主要从事骨科方面的临床研究。
基金资助:
河北省2014年度医学科学研究重点课题计划项目(Zl20140191)

Finite element analysis of knee flexion and extension movement

Zhang Li-chao¹, Zhang Li-min², Lv Yong-ming¹, Wang Zhi-hui¹, Yang Yang¹, Xu Fei¹, Dai Hai-feng¹, Li Jia¹, Cao Xiang-yu¹, Wu Li-zhu³

¹Department of Orthopedics, Affiliated Hospital, Chengde Medical College, Chengde 067000, Hebei Province, China; ²the 266 Hospital of Chinese PLA, Chengde 067000, Hebei Province, China; ³Department of Emergency, Chengde Central Hospital, Chengde 067000, Hebei Province, China

Revised:2016-11-17 Online:2017-01-28 Published:2017-03-14
About author:Zhang Li-chao, Master, Attending physician, Department of Orthopedics, Affiliated Hospital, Chengde Medical College, Chengde 067000, Hebei Province, China
Supported by:
the Key Medical Science Plan Program of Hebei Province in 2014, No. Zl20140191

摘要/Abstract

摘要：

文章快速阅读：

文题释义：

膝关节炎：是一种以退行性病理改变为基础的疾患。多见于中老年人群，其症状多表现为膝盖红肿痛、上下楼梯痛、坐起立行时膝部酸痛不适等。也会有患者表现肿胀、弹响、积液等，如不及时治疗，则会引起关节畸形。在膝关节部位还常患有膝关节滑膜炎、韧带损伤、半月板损伤、膝关节游离体、腘窝囊肿、髌骨软化、鹅足滑囊炎、膝内/外翻等关节疾病。

人工膝关节：人工膝关节置换是在近代人工髋关节成功应用于患者后逐渐发展起来的一种治疗膝关节疾病的新技术。人工膝关节完全参照了正常人膝关节的解剖形状，是一种仿生设计制品。模仿人体膝关节的结构及活动方式，人工膝关节由4个部件组成：股骨部分、胫骨部分、髌骨部分以及聚乙垫。

摘要

背景：三维有限元分析在骨科中被广泛应用，在膝关节屈曲运动中的研究不多。

目的：对膝关节屈曲运动的生物力学特性进行三维有限元分析。

方法：建立膝关节的三维有限元模型以及膝关节屈曲30°，60°以及120°时的三维有限元模型，分析膝关节不同屈曲位股骨前移情况、股骨内移情况、胫骨内旋情况和胫骨内翻情况。

结果与结论：①股骨前移：在膝关节屈曲30°，60°和120°时股骨均出现后移，膝关节屈曲60°时的股骨后移值大于屈曲30°和屈曲120°(P < 0.05)；膝关节屈曲120°股骨后移值大于屈曲30°，但差异无显著性意义(P > 0.05)；②股骨内移：在膝关节屈曲30°时股骨外移，在膝关节屈曲120°时膝关节内移，膝关节屈曲120°时股骨内移值明显大于屈曲30°和屈曲60°时(P < 0.05)；膝关节屈曲30°和屈曲60°时外移值比较差异无显著性意义(P > 0.05)；③胫骨内旋：在膝关节屈曲30°，60°和120°时胫骨均出现内旋，3者比较差异均无显著性意义(P > 0.05)；④胫骨内翻：在膝关节屈曲30°，60°和120°时胫骨均出现内翻，其中屈曲120°时胫骨内翻不明显，膝关节屈曲60°时胫骨内翻度数最大，明显高于屈曲120°时(P < 0.05)；膝关节屈曲30°时胫骨内翻度数大于屈曲120°时，但差异无显著性意义(P > 0.05)；⑤结果表明，膝关节在屈曲60°时股骨后移最明显，膝关节屈曲120°时股骨内移值最大，膝关节屈曲30°和屈曲120°时胫骨内旋明显，膝关节屈曲60°时胫骨明显内翻。

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

关键词: 骨科植入物, 数字化骨科, 膝关节, 力学模型, 屈曲运动, 有限元分析, 生物学特性

Abstract:

BACKGROUND: Three dimensional finite element analysis is widely used in orthopedics, but research in the knee flexion movement is not much.

OBJECTIVE: To analyze biomechanical properties of knee flexion using finite element analysis.

METHODS: Three dimensional finite element models and models of knee flexion at 30 degrees, 60 degrees and 120 degrees were established. Femoral forward movement, femur inward movement, tibia internal rotation and tibia vara were analyzed at different flexion angles.

RESULTS AND CONCLUSION: (1) When knee flexion was at 30 degrees, 60 degrees and 120 degrees, the femur had backward movement. The backward movement value was larger at 60 degrees than at 30 and 120 degrees (P < 0.05). The backward movement value was larger at 120 degrees than at 30 degrees (P > 0.05). (2) The femur had outward movement at 30 degrees, and inward movement at 120 degrees. The inward movement value was significantly larger at 120 degrees than at 30 and 60 degrees (P < 0.05). No significant difference in outward movement was determined between the 30 and 60 degrees (P > 0.05). (3) Tibia internal rotation was not significantly different at 30, 60 and 120 degrees of flexion (P > 0.05). (4) Tibia vara was found at 30, 60 and 120 degrees of flexion, but the tibia vara was not obvious at 120 degrees. The tibia vara was maximum at 60 degrees, and significantly higher than that at 120 degrees (P < 0.05). The tibia vara was larger at 30 degrees than at 120 degrees (P > 0.05). (5) These findings verify that backward movement was most obvious at 60 degrees. Femur inward movement was most large at 120 degrees. Tibia internal rotation was noticeable at 30 and 120 degrees. Tibia vara was remarkable at 60 degrees.

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

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

中图分类号:

R318

张立超，张立敏，吕永明，王智慧，杨阳，徐飞，戴海峰，李嘉，曹向宇，武丽珠. 膝关节屈伸运动形态的有限元分析[J]. 中国组织工程研究, 2017, 21(3): 396-400.

Zhang Li-chao, Zhang Li-min, Lv Yong-ming, Wang Zhi-hui, Yang Yang, Xu Fei, Dai Hai-feng, Li Jia, Cao Xiang-yu, Wu Li-zhu. Finite element analysis of knee flexion and extension movement[J]. Chinese Journal of Tissue Engineering Research, 2017, 21(3): 396-400.

图/表（结果） 2

参考文献

[1] 翟云,高根德,徐守宇. 膝关节骨关节炎的基础研究进展[J]. 中国骨伤,2012,28(1):83-87.

[2] 徐小敏.江苏省高校高水平篮球运动员膝关节损伤现状及对策研究[D].苏州大学,2010.

[3] 黄伟华.人体膝关节有限元建模及损伤生物力学探索性研究[D].南方医科大学,2014.

[4] 樊星.等速肌力测定分析膝关节骨关节炎对膝关节屈伸肌群影响[D].中南大学,2013.

[5] 徐小会.全膝关节置换术治疗晚期类风湿性关节炎的中期随访研究[D].第四军医大学,2014.

[6] 郑瑜峰.人工膝关节置换前后髌骨高度对膝关节功能的影响[D].第四军医大学,2014.

[7] Ip D.Does addition of low-level laser therapy (LLLT) in conservative care of knee arthritis successfully postpone the need for joint replacement? Lasers Med Sci. 2015;30(9): 2335-2339.

[8] Barlow T, Dunbar M, Sprowson A, et al.Development of an outcome prediction tool for patients considering a total knee replacement—the Knee Outcome Prediction Study (KOPS).BMC Musculoskelet Disord. 2014;15:451. doi: 10.1186/1471-2474-15-451.

[9] 王国伟. 老年骨关节炎患者全膝关节置换术后疗效影响因素分析[D].山东大学,2008.

[10] 王鑫瑞,王涛,张志强.全膝关节置换围置换期氨甲环酸的临床应用[J].中国组织工程研究,2015,19(44):7188-7192.

[11] 胡勇,谢辉,杨述华.三维有限元分析在脊柱生物力学中应用研究[J].医用生物力学,2006,21(3):246-250.

[12] 王一民. 跟骨骨折的有限元分析[D].南方医科大学,2013.

[13] Razfar N, Reeves JM, Langohr DG,et al.Comparison of proximal humeral bone stresses between stemless, short stem, and standard stem length: a finite element analysis[J].J Shoulder Elbow Surg. 2016 Jan 22. pii: S1058-2746(15) 00625-4. doi: 10.1016/j.jse.2015.11.011.

[14] Allred JJ, Flores-Hernandez C, Hoenecke HR Jr,et al. Posterior augmented glenoid implants require less bone removal and generate lower stresses: afinite element analysis.J Shoulder Elbow Surg. 2016 Jan 14. pii: S1058-2746(15)00574-1. doi: 10.1016/j.jse.2015.10.003.

[15] Harith H, Schmutz B, Malekani J,et al.Can we safely deform a plate to fit every bone? Population-based fit assessment and finiteelement deformation of a distal tibial plate[J].Med Eng Phys. 2015 Dec 28. pii: S1350-4533(15)00270-2. doi: 10.1016/j.medengphy.2015.11.012.

[16] Wang B,Li L,Chen A.Research on method of rule-based finite element quick modeling//Networking and Digital Society (ICNDS),2010 2nd International Conference on IEEE,2010;2: 540-543.

[17] 郑锋,林海滨.人骨三维有限元分析的研究进展[J]. 医学综述, 2011,17(11):1689-1691.

[18] 周伊顗.三维建模与有限元分析在眼科和骨科的应用[D].清华大学,2005.

[19] 何川.前交叉韧带虚拟重建系统的应用及术后膝关节内结构的三维有限元研究[D].昆明医科大学,2014.

[20] 马雪梅,张春秋,高丽兰,等. 骑行时不同屈曲角度膝关节软骨受力分析[J].中国生物医学工程学报,2016,35(1):119-123.

[21] 齐玮. 高屈曲状态下膝关节生物力学研究[D].中国人民解放军医学院,2013.

[22] 王西十,白瑞蒲.关于人膝关节生物力学模型的研究现状[J]. 力学进展,1999,29(2):101-107.

[23] 孙永森.人体膝关节生物力学模型研究[D].大连理工大学,2006.

[24] 许玉林,孙允高,张春林. 膝关节生物力学模型研究现状[J]. 医用生物力学,2003,18(4):250-255.

[25] 周吉彬,郝智秀,杨义勇,等.人体膝关节韧带生物力学研究进展[J]. 生物医学工程学杂志,2006,23(4):903-906.

[26] 刘爱峰,马剑雄,魏强,等.膝关节生物力学模型研究进展[J]. 中国矫形外科杂志,2012,20(19):1774-1776.

[27] 孙永森. 人体膝关节生物力学模型研究[D].大连理工大学,2006.

[28] Gokkus K, Atmaca H, U?ur L,et al.The relationship between medial meniscal subluxation and stress distribution pattern of the knee joint: Finite element analysis.J Orthop Sci. 2016; 21(1):32-37.

[29] Zhu GD, Guo WS, Zhang QD, et al.Finite Element Analysis of Mobile-bearing Unicompartmental Knee Arthroplasty: The Influence of Tibial Component Coronal Alignment.Chin Med J (Engl). 2015;128(21):2873-2878.

[30] Tanska P, Mononen ME, Korhonen RK. A multi-scale finite element model for investigation of chondrocyte mechanics in normal and medial meniscectomy human knee joint during walking. J Biomech. 2015;48(8):1397-1406.

[31] Feikes JD, O'Connor JJ, Zavatsky AB. A constraint-based approach to modelling the mobility of the human knee joint. J Biomech. 2003;36(1):125-129.

[32] Martelli S, Ellis RE, Marcacci M, et al. Total knee arthroplasty kinematics: Computer simulation and intraoperative evaluation. J Arthroplasty. 1998;13(2):145-155.

[33] Wilson DR, Feikes JD, O'Connor JJ. Ligaments and articular contact guide passive knee flexion. J Biomech. 1998;31(12): 1127-1136.

[34] Galbusera F, Freutel M, Dürselen L, et al.Material models and properties in the finite element analysis of knee ligaments: a literature review.Front Bioeng Biotechnol. 2014;2:54. doi: 10.3389/fbioe.2014.00054.

[35] Sun ZH, Liu YJ, Li H.Femoral stress and strain changes post-hip, -knee and -ipsilateral hip/knee arthroplasties: a finite element analysis.Orthop Surg. 2014;6(2):137-144. doi: 10.1111/os.12105.

[36] Krouskop TA,Trono R,Adamski AJ.The effects of pylon shape on bone-pylon interface performance in direct skeletal attachment. J Biomed Mater Res.1976;10(3):345-369.

[37] Bendjaballah MZ,Shirazi-Adl A,Zukor DJ.Biomechanics of the human knee joint in compression:reconstruction,mesh generation and finite element analysis. Knee.1995;2(2): 69-79.

引言

膝关节是下肢活动枢纽，负重大，运动多，是主要的下肢关节，膝关节是人体关节中功能和结构最复杂的关节，由髌骨、股骨、胫骨和腓骨，以及周围的韧带和肌肉等组成，是结构完备的滑膜关节，具有关节囊、关节腔和关节面滑膜关节必备结构，以及软骨、半月板、滑囊、韧带和脂肪垫等辅助结构，膝关节的运动是各个关节的相互协调运动。因为膝关节的运动量比较大，结构比较复杂，因此膝关节的损伤位居所有关节的首位，研究膝关节生物力学特点，对于预防膝关节炎、膝关节损伤具有重要意义^[1-5]。人工膝关节置换术是晚期关节炎常用的一种比较有效的方法，能够减轻膝关节疼痛，恢复膝关节功能，恢复关节力线，手术成功率比较高^[6-7]，但术后有出现假体过度磨损、假体松动、功能性失败等风险，膝关节置换术失败的原因主要有假体因素、手术操作因素以及病理性的因素，研究膝关节的生物力学特性对于设计和制造膝关节假体也具有重要意义^[8-10]。三维有限元分析操作简单，可重复性好，在骨科生物力学研究领域被广泛应用^[11-18]。何川^[19]采用有限元法对膝关节屈曲0°，30°，60°和90°时膝关节韧带的生物力学特点进行分析，发现膝关节韧带在膝关节不同屈曲角度时的应力大小和应力分布不同，发生损伤的风险也不同；马雪梅等^[20]对骑行时膝关节不同屈曲角度时膝关节软骨的受力情况进行有限元分析，发现膝关节屈曲角度大时应力比较大，髌骨软骨上侧和胫骨平台软骨后侧承受的应力比较大。何川和马雪梅等的研究对不同屈曲度数时膝关节韧带和膝关节软骨的力学特点进行了研究，但关于膝关节不同屈曲位时股骨和胫骨的生物力学特点没有涉及，文章对膝关节不同屈曲位时股骨和胫骨的生物力学特点进行研究。有研究将膝关节的屈曲度分为低屈曲度、中屈曲度和高屈曲度^[21]，屈曲度< 30°为低屈曲度，屈曲度在30°-120°为中屈曲度，屈曲度> 120°为高屈曲度，结合何川^[19]的膝关节屈曲度划分方法，文章对膝关节屈曲30°，60°，120°时股骨和胫骨的生物力学特点进行研究。

文章通过建立膝关节屈曲30°，60°，120°时的三维有限元模型，分析膝关节不同屈曲位股骨前移、内移情况以及胫骨内旋、胫骨内翻情况。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

材料方法

1.1 设计 三维有限元分析实验。

1.2 时间及地点 于2014年10月至2015年11月在承德医学院附属医院实验室完成。

1.3 材料 Mimics 15.0软件，购自于比利时Materialise公司。

1.4 对象 选取健康男性39岁志愿者1名，体质量78 kg，身高176 cm，志愿者无内外科疾病史，没有肢体残疾和外伤史，CT扫描和拍X射线片示膝关节正常。

1.5 方法

获得志愿者膝关节影像资料：志愿者取仰卧位，CT扫描机按照0.6 mm的层厚进行志愿者膝关节扫描，扫描范围：膝关节股骨下端1/3至膝关节胫骨上段1/3。将扫描的膝关节图像以DICOM格式保存。

建立膝关节三维模型：将CT扫描机扫描的DICOM格式的膝关节图像导入Mimics软件，将膝关节图像进过无损模型转换，填补图像空洞，生成膝关节骨骼蒙罩，对蒙罩进行编辑，生成三维模型，对三维模型进行包裹、平滑，勾画出膝关节骨骼以及韧带组织蒙罩，生成膝关节的三维有限元模型(见图1)，以Binary STL格式保存。

获得志愿者膝关节不同屈曲位置影像资料：CT扫描机按照0.6 mm的层厚扫描志愿者膝关节屈曲30°，60°和120°时的图像，扫描范围：膝关节股骨下端1/3至膝关节胫骨上段1/3。将扫描的膝关节图像以DICOM格式保存。

建立膝关节屈曲三维模型：将CT扫描机扫描的DICOM格式的膝关节屈曲图像导入Mimics软件，将膝关节屈曲图像进过无损模型转换，填补图像空洞，生成膝关节屈曲骨骼蒙罩，对蒙罩进行编辑，生成膝关节屈曲三维模型，对膝关节屈曲三维模型进行包裹、平滑，勾画出膝关节屈曲骨骼以及韧带组织蒙罩，生成膝关节屈曲的三维有限元模型(见图1)，以Binary STL格式保存。

网络划分：将建立的膝关节模型进行四面体网络划分。

膝关节材料赋值设置：股骨的弹性模量为14 357 MPa，松比设置为0.34；髌骨的弹性模量为14 357 MPa，松比设置为0.34；胫骨的弹性模量为18 456 MPa，松比设置为0.31；软骨的弹性模量为109 MPa，松比设置为0.19；前后交叉韧带的弹性模量为302 MPa，松比设置为0.29；股-髌韧带的弹性模量为447 MPa，松比设置为0.29；内外侧韧带的弹性模量为300 MPa，松比设置为0.29；周围肌肉的弹性模量为413 MPa，松比设置为0.29。

1.6 主要观察指标 膝关节不同屈曲位股骨前移、股骨内移、胫骨内旋和胫骨内翻情况。

1.7 统计学分析 采用SPSS软件进行分析，均数两两之间比较采用t 检验，取P=0.05为显著标准。

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

讨论

膝关节活动范围比较广泛，是人体最复杂的关节，由3个骨平台、2个软骨平面、5个关节韧带以及滑液和结缔组织构成，为双关节结构，具有平面关节的特性和球窝关节功能。人体膝关节损伤的原因主要有伸屈关节运动，剪切运动，跳跃或负重物时受力不均匀，应力的突然急剧变化，机械碰撞或摔跌伤等，膝关节的损伤在生物力学上主要是直接或间接暴力作用于膝关节骨组织或软骨组织时引起的膝关节微小结构裂伤或折断。由于人体结构比较复杂，对人体活体的研究具有一定的局限性，因此直接对人体膝关节的生物力学特点进行研究比较困难，在尸体上对膝关节进行研究又缺乏人体运动时的生理环境，因此建立膝关节的力学和解剖学模型是研究膝关节生物力学特点的有效方法。关于膝关节生物力学研究以往常用的方法有动物膝关节模型、物理膝关节模型、数字膝关节模型以及虚拟膝关节模型等^[22-26]，但由于上述模型均具有一定的局限性，没有对骨架和关节软组织进行综合考虑，和膝关节的解剖结构和生物状态参数有比较大的差距，不能正确的反应膝关节的生物力学特点，膝关节三维有限元分析能够克服上述方法的不足，在膝关节的研究中被广泛应用^[27-31]。

随着技术的不断发展，膝关节力学模型也不断发展。物理模型是比较早的膝关节力学模型，最常见的膝关节物理模型为将膝关节股经关节运动视为交叉四连杆结构模型^[32]，在胫骨平台和股骨髁接触时，交叉韧带长度基本不变，由交叉韧带长度和交叉韧带在股骨以及胫骨上的附着位置计算股骨髁与胫骨平台的几何形状，将膝关节简化为股胫关节二维模型^[33]，膝关节交叉四连杆模型在膝关节运动力学研究的早期具有重要意义。随后出现了膝关节股胫关节四连杆结构的三维模型，在膝关节股胫关节四连杆结构的三维模型中将膝关节屈曲过程中，前后交叉韧带和内侧副韧带设定为不能伸缩的线约束，胫骨髁设立为刚性平面，股骨后髁设定为刚性球面，在屈膝过程中股胫关节面处于接触状态，整个膝关节结构为5个约束把刚性胫骨和刚性股骨相连起来，膝关节股胫关节四连杆结构三维模型和膝关节股胫关节四连杆结构二维模型共同存在的缺陷是不能描述股胫关节在伸膝和屈膝时的内外翻运动和内外旋转运动，考虑到交叉韧带的中间纤维长度在高度屈膝时能够发生改变，因此该模型在高度屈膝时不适用。国内关于膝关节力学模型的研究比较少，主要集中的膝关节的二维模型的研究上。随着生物力学和计算机技术的不断进步，有限元分析发也不断完善，有限元分析法能够节约研究时间，降低研究成本，能够对活体组织进行生物力学分析，在膝关节生物力学的研究中的应用也不断成熟^[34-36]。何川^[19]通过建立膝关节有限元模型，对膝关节不同屈曲角度时的膝关节韧带的生物力学进行研究，发现膝关节屈曲度数不同，膝关节韧带的应力大小和应力分布不同；马雪梅等^[20]通过建立膝关节有限元模型，对骑行时膝关节不同屈曲角度时膝关节软骨的受力情况进行有限元分析，发现膝关节屈曲角度不同膝关节软骨的受力也不同。何川和马雪梅等的研究没有涉及膝关节不同屈曲位时股骨和胫骨的生物力学特点，本文中对膝关节不同屈曲位时股骨和胫骨的生物力学特点进行研究。

膝关节的有限元模型最初只有简单的骨性连接模型，只包括胫骨平台和股骨远端，不包括髌骨在内^[37]，随后开始测量和分析对膝关节局部受力情况，以及测量骨松质和骨皮质的负载情况，数据主要由CT扫描得到，随着计算机的发展，膝关节模型中加入了髌骨，随后在膝关节模型中加入了半月板、胫骨软骨和股骨软骨，研究膝关节伸展拉力，随着影像技术的不断发展，磁共振获得的膝关节数据比CT获得的膝关节数据更加精确，有研究开始用磁共振扫描获得膝关节的数据资料。近年来在膝关节模型中加入了膝关节周围韧带和肌肉等组织，对膝关节的稳定性进行模拟。建立膝关节有限元模型的建模方法主要有：①人工建模法：在20世纪没有将二维影像学图像转化为三维图像的软件，三维建模软件和影像图像软件是分开的，三维模型的建立主要靠肉眼分辨，并手动描绘三维轮廓，在进行网络设置和划分，由于人工建模法比较费时，现在基本被淘汰；②半自动建模法：随着技术的发展，三维建模软件出现多功能一体化，能够支持多种图像格式，使二维影像学图像向三维图像转换更加方便，更加快捷，能够节约大量人工时间，并且门槛比较低，比较常用的软件用Mimics构膜软件，Geomagic Studio影像三维检测软件，以及ABAQOUS三维有限元处理软件，半自动建模法是目前膝关节三维有限元分析的主流；③全自动建模法：全自动建模法是一种新的技术，全自动建模法只需患者的影像学资料，和数据库比对，根据和数据库的差异，将数据库中的标准模型根据患者的影像学资料进行修改，自动生成患者的三维有限元模型，改建模方法仍在测试阶段。

膝关节在不同位置屈曲过程中生物力学均发生变化，在膝关节伸直位时，膝关节的稳定性主要由前交叉韧带前内侧束维持，内外侧副韧带发挥次要作用，后交叉韧带基本不发挥作用，在屈曲度数不断增加时，膝关节的稳定性主要有前交叉韧带前内侧束河后交叉韧带维持，前交叉韧带后外侧束基本没有太大作用，在屈曲120°时膝关节的稳定主要由后交叉韧带和前交叉韧带前内侧束维持，内外侧副韧带基本不发挥作用。在膝关节不同屈曲位时，股骨和胫骨的位置也发生相应变化，为研究膝关节不同屈曲过程中股骨和胫骨的生物力学变化，本实验对膝关节屈曲运动的生物力学特性进行三维有限元分析，采用半自动建模法建立膝关节的三维有限元模型以及膝关节屈曲30°、屈曲60°以及屈曲120°时的三维有限元模型，分析膝关节不同屈曲位股骨前移情况、股骨内移情况、胫骨内旋情况和胫骨内翻情况。发现在膝关节屈曲30°、屈曲60°和屈曲120°时股骨均出现后移，膝关节屈曲60°时的股骨后移值大于屈曲30°和屈曲120°，膝关节屈曲120°股骨后移值大于屈曲30°，但差异无显著性意义。在膝关节屈曲30°时股骨外移，在膝关节屈曲120°时膝关节内移，膝关节屈曲120°时股骨内移值明显大于屈曲30°和屈曲60°时，膝关节屈曲30°和屈曲60°时外移值比较差异无显著性意义。在膝关节屈曲30°、屈曲60°和屈曲120°时胫骨均出现内旋，胫骨内旋度数在膝关节屈曲30°、屈曲60°和屈曲120°时比较差异均无显著性意义。在膝关节屈曲30°、屈曲60°和屈曲120°时胫骨均出现内翻，其中屈曲120°时胫骨内翻不明显，膝关节屈曲60°时胫骨内翻度数最大，明显高于屈曲120°时，膝关节屈曲30°时胫骨内翻度数大于屈曲120°时，但差异无显著性意义。由此可见，膝关节在屈曲60°时股骨后移最明显，膝关节屈曲120°时股骨内移值最大，膝关节屈曲30°和屈曲120°胫骨内旋明显，膝关节屈曲60°时胫骨明显内翻。膝关节在30°，60°，120°不同屈曲度数运动过程中，股骨和胫骨在不同屈曲度数时的相对运动特点不同，股骨和胫骨的相对运动特点为膝关节假体的设计以及膝关节损伤后的功能康复锻炼提供了一定的生物力学依据。

文章的不足之处在于：①没有建立半月板模型，未考虑半月板对股骨和胫骨生物力学的影响；②没有将人种、性别和体质量的差异对研究的影响；③是对膝关节准静态屈曲运动的研究。在后续的研究中将考虑把半月板模型设计到研究中去，并注意人种、性别和体质量的差异，对膝关节动态屈曲运动过程中的生物力学进行研究分析。

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

膝关节屈伸运动形态的有限元分析

Finite element analysis of knee flexion and extension movement

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 2

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	王建平, 张晓辉, 余进伟, 魏绍亮, 张新民, 许幸新, 曲海军. 三维图像配准及坐标变换膝关节运动分析在机构学中的应用[J]. 中国组织工程研究, 2022, 26(在线): 1-5.
[2]	魏国强, 李云峰, 王一, 牛晓芬, 车丽芳, 王海燕, 李志军, 史国鹏, 白灵, 莫凯, 张晨晨, 许阳阳, 李筱贺. 非均匀材料股骨在不同负荷情况下的生物力学分析[J]. 中国组织工程研究, 2022, 26(9): 1318-1322.
[3]	吕倩忆, 陈芯仪, 郑慧娥, 何灏龙, 李棋龙, 陈楚淘, 田浩梅. 不同角度肘按法对正常人腰椎椎体及后部结构的应力和位移分析[J]. 中国组织工程研究, 2022, 26(9): 1346-1350.
[4]	张玉芳, 吕蒙, 梅钊. 青少年脊柱侧弯全脊柱生物力学模型的构建及验证[J]. 中国组织工程研究, 2022, 26(9): 1351-1356.
[5]	张吉超, 董跃福, 牟志芳, 张震, 李冰言, 徐祥钧, 李佳意, 任梦, 董万鹏. 骨关节炎患者在不同步态角度下膝关节内部生物力学变化的有限元分析[J]. 中国组织工程研究, 2022, 26(9): 1357-1361.
[6]	吴冰霜, 汪志, 唐懿, 唐晓俞, 李棋. 前交叉韧带重建：从腱骨止点到腱骨愈合的研究进展[J]. 中国组织工程研究, 2022, 26(8): 1293-1298.
[7]	郭嘉, 丁琼桦, 刘泽, 吕思懿, 周泉程, 高玉花, 白春雨. 间充质干细胞来源外泌体的生物学特性及免疫调控作用[J]. 中国组织工程研究, 2022, 26(7): 1093-1101.
[8]	刘冬铖, 赵继军, 周子红, 吴沼锋, 俞颖豪, 陈宇浩, 冯德宏. 开放楔形胫骨高位截骨手术不同力线矫正参考方法的比较[J]. 中国组织工程研究, 2022, 26(6): 827-831.
[9]	邵洋洋, 张峻霞, 姜美姣, 刘泽龙, 高昆, 于淑晗. 青年男性佩戴护膝跑步优势侧和非优势侧下肢关节的运动学特征[J]. 中国组织工程研究, 2022, 26(6): 832-837.
[10]	温明韬, 梁学振, 李嘉程, 许波, 李刚 . 两种方式固定Sanders Ⅱ型跟骨骨折后的力学稳定性[J]. 中国组织工程研究, 2022, 26(6): 838-842.
[11]	王海龙, 李龙, 买合木提·亚库甫, 陈洪涛, 刘旭, 伊力哈木·托合提. 髋臼假体在Lewinnek 安全区内应力分布的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 843-847.
[12]	黄浩, 洪嵩, 瓦庆德. 全膝关节置换过程中股骨假体旋转对髌股关节接触压影响的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 848-852.
[13]	魏兵, 常山. 脊柱骨折矢状面不同角度置钉方式的有限元分析[J]. 中国组织工程研究, 2022, 26(6): 864-869.
[14]	袁景, 孙效虎, 陈慧, 乔永杰, 王立新. 成人继发性膝外翻畸形股骨远端的数字化测量与分析[J]. 中国组织工程研究, 2022, 26(6): 881-885.
[15]	周建国, 刘世伟, 袁长红, 毕声荣, 杨国平, 胡伟全, 刘辉, 钱锐. 后交叉韧带保留型假体全膝关节置换治疗膝骨性关节炎合并膝外翻畸形[J]. 中国组织工程研究, 2022, 26(6): 892-897.