不同骨强度下全膝置换过程中发生股骨前皮质切迹的三维有限元分析

doi:10.12307/2025.144

中国组织工程研究 ›› 2025, Vol. 29 ›› Issue (9): 1775-1782.doi: 10.12307/2025.144

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

不同骨强度下全膝置换过程中发生股骨前皮质切迹的三维有限元分析

周金海1，李江伟2，王序全1，庄颖1，赵瑛1，杨渝勇1，王嘉嘉1，杨阳1，周仕炼1

1贵黔国际医院骨科，贵州省贵阳市 550024；2贵州医科大学附属医院骨科，贵州省贵阳市 550004

收稿日期:2023-12-19 接受日期:2024-02-29 出版日期:2025-03-28 发布日期:2024-10-09
通讯作者: 李江伟，男，主任医师，贵州医科大学附属医院骨科，贵州省贵阳市 550004
作者简介:周金海，男，1994年生，贵州省赫章县人，汉族，硕士，主治医师，主要从事骨与关节损伤方面的研究。
基金资助:
贵州省卫健委科学基金项目(gzwkj2021-507)，项目负责人：庄颖

Three-dimensional finite element analysis of anterior femoral notching during total knee arthroplasty at different bone strengths

Zhou Jinhai1, Li Jiangwei2, Wang Xuquan1, Zhuang Ying1, Zhao Ying1, Yang Yuyong1, Wang Jiajia1, Yang Yang1, Zhou Shilian1

1Department of Orthopedics, Guiqian International Hospital, Guiyang 550024, Guizhou Province, China; 2Department of Orthopedics, Affiliated Hospital of Guizhou Medical University, Guiyang 550004, Guizhou Province, China

Received:2023-12-19 Accepted:2024-02-29 Online:2025-03-28 Published:2024-10-09
Contact: Li Jiangwei, Chief physician, Department of Orthopedics, Affiliated Hospital of Guizhou Medical University, Guiyang 550004, Guizhou Province, China
About author:Zhou Jinhai, Master, Attending physician, Department of Orthopedics, Guiqian International Hospital, Guiyang 550024, Guizhou Province, China
Supported by:
Guizhou Provincial Health Commission Science Fund Project, No. gzwkj2021-507 (to ZY)

摘要/Abstract

摘要：

文题释义：
股骨前皮质切迹：是指全膝关节置换过程中股骨前髁截骨不良导致的股骨前方皮质受损现象，可能是截骨板位置放置不良、后参考技术的应用、后髁角增大、股骨前皮质解剖形态变异、股骨假体尺寸选择性缩小、后移入口点的位置等原因导致的；其与膝关节股骨假体周围骨折存在高风险性。切迹的深度是引起骨折风险的直接因素之一。
膝关节股骨假体周围骨折：指全膝关节置换后出现膝关节线15 cm范围内或者距离假体近端5 cm范围内的股骨远端骨折。近年来膝关节股骨假体周围骨折的发生率越来越高，导致其发生的原因与术中的失误操作关系密切，有研究表明股骨前皮质切迹是主要原因之一。

摘要
背景：全膝关节置换后发生膝关节股骨假体周围骨折是常见并发症之一，目前对不同骨质强度条件下发生膝关节股骨假体周围骨折的生物力学研究欠缺，三维有限元分析可为临床提供生物力学基础。
目的：探讨不同骨质强度下全膝关节置换过程中发生股骨前皮质切迹的生物力学变化，为临床预防膝关节置换后发生膝关节股骨假体周围骨折提供力学理论基础。
方法：获取健康成年人的股骨CT数据，运用Mimics、Geomagic studio、Solidworks软件建立膝关节股骨侧置换三维模型，然后构建不同深度的股骨前皮质切迹模型，将模型导入到ANSYS软件并分析不同骨强度、不同股骨前皮质切迹深度对股骨髁上生物应力的影响，并分析股骨前皮质切迹骨水泥填补前、后的股骨前髁截面应力变化。
结果与结论：①任何骨质强度下，髁上应力都随股骨前皮质切迹深度加深而增加；在正常骨质情况下，当股骨前皮质切迹深度在3 mm和4 mm间有一个应力突变点；在骨质疏松情况下，当股骨前皮质切迹深度在2 mm和3 mm间有一个应力突变点；②膝关节置换过程中发生股骨前皮质切迹且深度超过骨皮质厚度时，随着骨质强度降低，股骨髁上应力逐渐增大；③使用骨水泥填补股骨前皮质切迹深度为3 mm的模型，股骨前髁截面应力下降；④结果显示，膝关节置换过程中应该避免股骨前皮质切迹出现，特别是骨质疏松患者；如术中出现股骨前皮质切迹且深度超过骨皮质厚度时，可以采用骨水泥均匀填充股骨前皮质切迹以降低股骨髁上应力，降低膝关节股骨假体周围骨折的发生率。

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

关键词: 全膝关节置换, 股骨前皮质切迹, 骨质疏松, 膝关节股骨假体周围骨折, 生物力学, 三维有限元分析

Abstract: BACKGROUND: Periprosthetic fracture of the femoral of the knee after total knee arthroplasty is one of the common complications, and there is a lack of biomechanical research on the periprosthetic fractures of the femoral of the knee under different bone strength conditions. The three-dimensional finite element analysis can provide a biomechanical basis for clinical practice.
OBJECTIVE: To investigate the biomechanical changes of anterior femoral notching after total knee arthroplasty under different bone strengths, and to provide a mechanical basis for the clinical prevention of supracondylar femoral periprosthetic fractures after knee arthroplasty.
METHODS: The femoral CT data of healthy adults were obtained, and the three-dimensional model of femoral lateral replacement of the knee joint was established by Mimics, Geomagic studio, and Solidworks software. Anterior femoral notching models of different depths were constructed, and the models were imported into ANSYS software to analyze the changes of biological stress on the femoral condyle with different bone strengths and different anterior femoral notching depths. The stress changes of the femoral anterior condyle section after and before the filling of anterior femoral notching with bone cement were analyzed.

RESULTS AND CONCLUSION: (1) Under any bone strength, the supracondylar stress increased with the depth of anterior femoral notching. In normal bone conditions, there was a stress abrupt change point when the anterior femoral notching depth was between 3 mm and 4 mm. In the case of osteoporosis, there was a stress abrupt point when the anterior femoral notching depth was between 2 mm and 3 mm. (2) When anterior femoral notching occurred during knee arthroplasty and the depth exceeded the thickness of the bone cortex, the supracondylar stress of the femoral gradually increased as the bone strength decreased. (3) The stress of the anterior femoral condyle section decreased when the model with an anterior femoral notching depth of 3 mm was filled with bone cement. (4) The results show that anterior femoral notching should be avoided during knee arthroplasty, especially in patients with osteoporosis. If anterior femoral notching occurs during surgery, bone cement can be used to evenly fill the anterior femoral notching to reduce the supracondylar stress of the femur and reduce the incidence of periprosthetic fractures of the femoral joint

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

Key words: total knee arthroplasty, anterior femoral notching, osteoporosis, periprosthetic fracture of knee joint and femur, biomechanics, three-dimensional finite element analysis

中图分类号:

周金海, 李江伟, 王序全, 庄颖, 赵瑛, 杨渝勇, 王嘉嘉, 杨阳, 周仕炼. 不同骨强度下全膝置换过程中发生股骨前皮质切迹的三维有限元分析[J]. 中国组织工程研究, 2025, 29(9): 1775-1782.

Zhou Jinhai, Li Jiangwei, Wang Xuquan, Zhuang Ying, Zhao Ying, Yang Yuyong, Wang Jiajia, Yang Yang, Zhou Shilian. Three-dimensional finite element analysis of anterior femoral notching during total knee arthroplasty at different bone strengths[J]. Chinese Journal of Tissue Engineering Research, 2025, 29(9): 1775-1782.

图/表（结果） 7

2.1 股骨髁上应力分析云图 在ANSYS 2023 R1 软件中对模型进行网格化并赋予不同材料值，给予外力载荷后得到应力云图如图8-10所示。

2.2 股骨髁上应力特点
2.2.1 实验结果汇总由表4、图11可看出，任何骨质强度下，髁上应力都随AFN深度加深而增加。同一骨质条件下，随着骨质条件下降，股骨髁上应力也逐渐增大；在正常骨质情况下，当AFN深度在3 mm和4 mm间有一个应力突变点；在骨质疏松情况下，当AFN深度在2 mm和3 mm间有一个应力突变点。

2.2.2 同种AFN深度情况下骨质强度对股骨髁上应力的影响分析将每一种AFN深度下骨质条件对髁上应力的影响分别进行分析，将3种骨质条件的弹性模量作为自变量设为X，股骨髁上最大等效应力作为因变量Y，选择单元线性回归统计分析，分析自变量X与因变量Y之间的关系。将每种AFN深度情况下的数据导入SPSS 25.0进行单元线性回归分析，得到分析数据如表5所示。

由表5可以看出，在AFN深度为3 mm时，P值为0.044 < 0.05，具有统计学意义；回归方程系数为-0.03，表明切迹深度超过骨皮质(2.96 mm)厚度时，骨质强度与髁上最大等效应力存在负相关关联；即骨质强度越低，股骨髁上应力越大。
2.2.3 同种骨质强度下AFN深度对髁上应力的影响分析将每一种骨质强度下的AFN深度对髁上应力的影响单独进行分析，将AFN深度H作为自变量设为X，股骨髁上最大等效应力作为因变量Y，选择单元线性回归统计分析，分析自变量X与因变量Y之间的关系。将数据导入SPSS 25.0进行单元线性回归分析，得到分析数据如表6所示。

由表6可看出，每种骨质条件下的预测模型拟合度均> 0.9，拟合度极佳；回归方程系数均为正值，说明此回归方程为正相关方程，自变量与因变量之间存在正相关关联；3种情况P值均< 0.05，存在统计学意义。总结上述可得出：不同骨质强度下，AFN深度明显影响股骨髁上最大等效应力，且随着AFN深度加深，股骨髁上最大等效应力随之增加。
2.3 骨水泥填补AFN对股骨髁上最大等效应力的影响 AFN深度为0，3 mm及骨水泥填补深度为3 mm时髁上最大等效应力结果如表7所示。将上述3组数据进行单因素方差分析，结果如下：

2.3.1 AFN深度为0 mm和3 mm的比较如图12所示，从表7中可看出，正常没有AFN组和AFN深度为3 mm组股骨髁上最大等效应力分别为(35.99±1.62)及(57.76±6.90) MPa，差异近12 MPa，差异有显著性意义(P=0.02 < 0.05)。
2.3.2 AFN深度为3 mm和采用骨水泥填充AFN比较如图13所示，AFN深度为3 mm+骨水泥、AFN深度为3 mm时股骨髁上最大等效应力分别为(42.40±2.37) MPa及(57.76±6.90) MPa，用骨水泥填补AFN后髁上最大等效应力降低，差异有显著性意义(P=0.018 < 0.05)。
2.3.3 正常没有AFN和AFN深度为3 mm+骨水泥的比较如图14所示，没有AFN组和骨水泥填补3 mm深度AFN组的股骨髁最大等效应力分别为(35.99±1.62) MPa和(42.40±2.37) MPa，相差不到7 MPa，在3组数据的单因素方差分析的多重比较中，两个样本差异无显著性意义(P=0.119 > 0.05)。

参考文献

[1] WANG JC, PIPLE AS,HILL WJ, et al. Computer-Navigated and Robotic-Assisted Total Knee Arthroplasty: Increasing in Popularity Without Increasing Complications. J Arthroplasty. 2022;37(12):2358-2364.
[2] LOMBARDO DJ, SILJANDER MP, SOBH A, et al. Periprosthetic fractures about total knee arthroplasty. Musculoskelet Surg. 2020;104(2):135-143.
[3] CARON É, GABRION A, EHLINGER M, et al. French society of orthopedic surgery and traumatology (SOFCOT); Complications and failures of non-tumoral hinged total knee arthroplasty in primary and aseptic revision surgery: A review of 290 cases. Orthop Traumatol Surge Res. 2021;107(3):102875.
[4] MA LL, YU XR, WENG XS, et al. Possible Risk Factors for Severe Complications Occurring after Primary Total Knee Arthroplasty. Chin Med Sci J. 2022;37(4):303-308.
[5] DELSMANN MM, SCHMIDT C, MÜHLENFELD M, et al. Prevalence of osteoporosis and osteopenia in elderly patients scheduled for total knee arthroplasty. Arch Orthop Trauma Surg. 2022;142(12):3957-3964.
[6] SKIBICKI HE, PONZIO DY, BRUSTEIN JA, et al. A cautionary case: osteoporotic femur fracture after robotic-assisted total knee arthroplasty. Osteoporos Int . 2021;32(10):2125-2129.
[7] STAMIRIS D, GKEKAS NK, ASTERIADIS K, et al. Anterior femoral notching ≥ 3 mm is associated with increased risk for supracondylar periprosthetic femoral fracture after total knee arthroplasty: a systematic review and meta-analysis. Eur J Orthop Surg Traumatol. 2022;32(3):383-393.
[8] ZHANG JC, ZHANG LS, ZHOU H, et al. Stress distribution patterns during the gait cycle in patients with anterior femoral notching following total knee replacement. BMC Musculoskelet Disord. 2022;23(1):718.
[9] CHATURONG P, AKRAPORN S, PAKPOOM R, et al. Risk factors of early periprosthetic femoral fracture after total knee arthroplasty. BMC Musculoskelet Disord. 2021;22(1):1009.
[10] PURANIK HG, MUKARTIHAL R, PATIL SS, et al. Does Femoral Notching During Total Knee Arthroplasty Influence Periprosthetic Fracture. A Prospective Study. J Arthroplasty. 2019;34(6):1244-1249.
[11] WANG DG, LI Y, YIN HL, et al. Three-dimensional finite element analysis of optimal distribution model of vertebroplasty. Ann Palliat Med. 2020; 9(3):1062-1072.
[12] XIE XJ, CAO SL, TONG K, et al. Three-dimensional finite element analysis with different internal fixation methods through the anterior approach. World J Clin Cases. 2021;9(8):1814-1826.
[13] 李志帅,张红倩,李丽,等.有膝关节有限元分析的研究热点及趋势[J].中国组织工程研究,2022,26(15):2412-2418.
[14] FALCINELLI C, WHYNE C. Image-based finite-element modeling of the human femur. Comput Methods Biomech Biomed Engin. 2020;23(14): 1138-1161.
[15] KANAIZUMI A, SUZUKI D, NAGOYA S, et al. Patient-specific three-dimensional evaluation of interface micromotion in two different short stem designs in cementless total hip arthroplasty: a finite element analysis. J Orthop Surg Res. 2022;17(1):437.
[16] 熊恒恒,聂伟志.三维有限元分析对骨关节和相关软组织损伤时应力状态的精准模拟[J].中国组织工程研究,2022,26(36):5875-5880.
[17] RHO JY, HOBATHO MC, ASHMAN RB. Relations of mechanical properties to density and CT numbers in human bone. Med Eng Phys. 1995;17(5): 347-355.
[18] HOU YJ. Investigation in the dependency of stiffness of cancellous bone on apparent density—based on the combination model of rod-rod structure and perforated plate structure. Sheng Wu Yi Xue Gong Cheng Xue Za Zhi. 2006;23(1):78-81.
[19] JETHANANDANI R, PATWARY MB, SHELLITO AD, et al. Biomechanical consequences of anterior femoral notching in cruciate-retaining versus posterior-stabilized total knee arthroplasty. Am J Orthop. 2016; 45(5):E268.
[20] WIJAYATHUNGA VN, JONES AC, OAKLAND RJ, et al. Development of specimen-specific finite element models of human vertebrae for the analysis of vertebroplasty. Proc Inst Mech Eng H. 2008;222(2):221-228.
[21] COMPLETO A1, SIMÕES JA, FONSECA F. Revision total knee arthroplasty: The influence of femoral stems in load sharing and stability. Knee. 2009; 16(4):275-279.
[22] WU WD, HAN ZH, HU BI, et al. A graphical guide for constructing a finite element model of the cervical spine with digital orthopedic software. Ann Transl Med. 2021;9(2):169.
[23] KAZEMI M, DABIRI Y, LI LP. Recent advances in computational mechanics of the human knee joint.Computational and mathematical methods in medicine. 2013;2013:718423.
[24] KIM JG, KANG KT, WANG JW. Biomechanical Difference between Conventional Transtibial Single-Bundle and Anatomical Transportal Double-Bundle Anterior Cruciate Ligament Reconstruction Using Three-Dimensional Finite Element Model Analysis. J Clin Med. 2021; 10(8):1625.
[25] MEHBOOB H, TARLOCHAN F, MEHBOOB A, et al. Finite element modelling and characterization of 3D cellular microstructures for the design of a cementless biomimetic porous hip stem. Mater Design. 2018;149: 101-112.
[26] ACEVEDO C, SYLVIA M, SCHAIBLE E, et al. Contributions of material properties and structure to increased bone fragility for a given bone mass in the UCD‐T2DM rat model of type 2 diabetes. J Bone Miner Res. 2018;33(6):1066-1075.
[27] GAUTAM A, CALLEJAS MA, ACHARYYA A, et al. Shape-memory-alloy-based smart knee spacer for total knee arthroplasty: 3D CAD modelling and a computational study. Med Eng Phys. 2018;55:43-51.
[28] LI Y, GAO YH, LU D, et al. Analysis of the effect of tibial torsion on tibial osteotomy in knee arthroplasty using a three-dimensional computed tomography-based modelling technique. BMC Musculoskelet Disord . 2019;20(1):361.
[29] PHILIPPE M, SARA C, INES K, et al. Physiological joint line total knee arthroplasty designs are especially sensitive to rotational placement – A finite element analysis. Plos One. 2018;13(2):e0192225.
[30] JEFFREY NK, KAETLYN RA, RICHARD FL. Diagnosis and Treatment of Hip and Knee Osteoarthritis: A Review. JAMA. 2021;325(6):568-578.
[31] BRYAN SB, JEREMIAH JM, SAMUEL FT, et al. Complication Rates in Total Knee Arthroplasty Performed for Osteoarthritis and Post-Traumatic Arthritis: A Comparison Study. J Arthroplasty. 2020;35(2):371-374.
[32] ZALZAL P, BACKSTEIN D, GROSS AE, et al. Notching of the anterior femoral cortex during total knee arthroplasty characteristics that increase local stresses. J Arthroplasty. 2006;21:737-743.
[33] CULP RW, SCHMIDT RG, HANKS G, et al. Supracondylar fracture of the femur following prosthetic knee arthroplasty. Clin Orthop Relat Res. 1987;(222):212-222.
[34] HOENIG T, KATHRYN EA, BELINDA RB, et al. Bone stress injuries. Nature reviews. Dis Primers. 2022;8(1):26.
[35] CEN HP, JIA YM, WU XG, et al. Effects of the microcrack’s shape, size and direction on the poroelasticbehaviors of a single osteon: a finite element study. Acta Bioeng Biomech. 2016;18(1):3-10.
[36] STOCKLEY I, MCAULEY JP, GROSS AE. Allograft reconstruction in total knee arthroplasty. Bone Joint J. 1992;74(3):393-397.
[37] RAND JA. Bone deficiency in total knee arthroplasty. Use of metal wedge augmentation. Clin Orthop Relat Res. 1991;271(271):63.
[38] RÖHNER E, HEINECKE M, MATZIOLIS G. Bone defect management in revision knee arthroplasty. Der Orthopade. 2021;50(12):1004-1010.
[39] AGGARWAL AK, BABURAJ V. Managing bone defects in primary total knee arthroplasty: options and current trends. Musculoskelet Surg. 2021;105(1):31-38.
[40] ZOHRE Z, MAJEDEH N, MAJID N, et al. Effect of Cement Type and Cementation Technique on the Retention of Implant-Supported Restorations. J Long Term Eff Med Implants. 2020;30(1):61-67.
[41] AJAY K, RAJESH G. Fracture Toughness of Acrylic PMMA Bone Cement: A Mini-Review. Indian J Orthop. 2021;55(5):1208-1214.

引言

全膝关节置换是膝骨关节炎的终末期治疗手段，随着全球老龄化及膝关节假体的不断更新换代，其手术量也逐年上升[1]；同时手术带来的并发症也逐年上升。膝关节股骨假体周围骨折是全膝关节置换术后常见并发症之一，不仅增加了患者的创伤及住院费用，同时导致膝关节功能丢失，是临床医师和患者需高度重视的问题[2-4]。相关研究表明导致膝关节股骨假体周围骨折发生的2个主要危险因素为股骨骨质强度和股骨前皮质切迹(anterior femoral
notching，AFN)[5-8]；然而一些研究提示正常骨质强度患者发生AFN与股骨髁上假体周围骨折之间无相关性[9-10]。骨质疏松患者全膝关节置换术中发生AFN的研究尚少，而临床中行膝关节置换患者大多为骨质疏松患者[5]，故骨质疏松膝关节置换术中发生AFN相关生物力学研究迫在眉睫。综上所述，目前AFN与膝关节股骨假体周围骨折的关系尚存争议，且对出现AFN情况的补救措施的研究尚缺乏，尤其是在对不同骨质强度发生AFN的“安全深度”缺乏足够的认识，从而可能增加膝关节股骨假体周围骨折的风险性。随着骨科数字化的迅速发展，三维有限元分析成为骨科相关力学研究中重要的工具之一。相关研究表明，有限元分析手段相比于实体力学研究存在可重复性高、更简洁化、更准确的优点，逐渐发展为骨科力学研究的主流工具，对骨科关节置换领域做出了巨大贡献[11-16]。
此次研究旨在通过三维有限元方法研究不同骨质强度下全膝关节置换术中发生AFN的生物力学变化及使用骨水泥填补AFN后的髁上应力变化，希望为临床预防全膝关节置换后发生膝关节股骨假体周围骨折和指导术中出现AFN的处理方案提供力学理论基础。
中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

材料方法

1.1 设计 单一样本的计算机模拟实验。
1.2 时间及地点 试验于2021-05-01/2023-05-01在贵黔国际医院骨科完成。
1.3 对象 选择1名骨质条件正常的青年男性，年龄30岁，身高175 cm，体质量70 kg。既往无膝关节创伤病史、膝关节无内外翻畸形、膝关节CT检查无骨折、畸形、肿瘤及无骨质疏松病史。另外选择1名男性膝关节置换术后志愿者膝关节CT数据，65岁，身高165 cm，体质量75 kg。
此次研究经贵黔国际医院医学伦理委员会批准(医院伦理批件号：贵黔2021伦审第(07)号，审批时间：2021-04-13)，志愿者均知情同意并自愿参与此次研究。
1.4 仪器设备 ①电脑型号：神舟 CP65S 笔记本电脑；②图像采集设备：128排双源螺旋CT；③图像编辑软件：Mimics 19.0(Materialises interactive medical images control system，Materiali，Leuve，Belgium)；Geomagic studio 2013(Raindrop Company，Marble Hill USA)；Solidworks 2017(Dassault Systemes Company Massachusetts USA)；
ANSYS 2023 R1(ANSYS Company，PA USA)。
1.5 实验方法

1.5.1 在Mimics进行股骨三维模型的建立及提取假体模型志愿者常规进行右下肢全长CT扫描，获取医学二维图像，将该数据以DICOM格式导出储存。在定义扫描解剖姿势(前后、左右、上下)方位一致后，交互式医学影像控制系统Mimics 19.0软件将自动输入DICOM格式的医学文件；采用剖面线的阂值分割方法(Profile line)，建立mask区分骨骼和其他组织结构，将其分割为独立的完整的右侧股骨体、股骨皮质骨3D模型、股骨松质骨3D模型(图1)；该志愿者股骨骨皮质厚度为2.96 mm。同法将患者膝关节置换术后CT数据重建成三维有限元模型，在Mimics 19.0软件里用区域增长把股骨髁假体从膝关节置换术后模型中把股骨假体STL模型提取出来(图2)。

1.5.2 处理假体、骨皮质和骨松质模型将假体、骨皮质、骨松质的STL文件格式分别导入Geomagic studio 2013软件，进行点云数据处理，最后分别进行假体数据模型、骨皮质、骨松质的多边形处理曲面拟合，见表1，为后续进行ANSYS精确分析打下基础。最后并以STP文件格式输出(图3)。

1.5.3 股骨侧假体置换不同深度股骨前皮质切迹模型的建立以CT扫描坐标系为原始工作坐标系，现将Geomagic studio 2013中Step格式的皮质骨、松质骨模型导入至Solidworks 2017中装配成一个完整的股骨模型，利用该软件模拟出股骨头的点云中心，与膝关节中心相连构成股骨机械轴。按照全膝关节置换的常规操作流程，进行膝关节股骨侧虚拟截骨，建立膝关节股骨置换模型(图4)。利用Solidworks软件分割功能制作不同深度的AFN模型(图5)，最后将各个模型以X-T格式进行保存，输出导入至ANSYS 2023 R1软件进行力学分析，得出股骨髁上区域的应力分布特点。

1.5.4 模型的有限元分析 ①将几何模型导入到ANSYS 2023 R1有限元分析软件中，查阅既往文献，国内外学者基于定量CT测定骨密度做出了分类标准；骨质的体积密度在240 mg/cm³以上者为正常骨质，骨量减少在25%以上为骨质疏松(根据诊治要求又可以分为轻、中度)，减少在37%以上为严重骨质疏松。此次实验采用定量CT测量的骨密度作为骨质疏松的诊断标准，以方便骨密度与不同骨质条件下骨松质、骨皮质的弹性模量换算，换算经验公式为E=1 780ρ1.88 [17-18]。在分析材料库中分别建立正常皮质骨、松质骨；轻中度骨质疏松皮质骨、松质骨；重度骨质疏松皮质骨、松质骨；假体的材料属性参数及Ⅲ型丙烯酸树酯骨水泥[19-20](表2)，进入Mechanical工作界面，在Geometry中对股骨皮质骨、松质骨、假体等模型赋予相关的材料。②接下来对模型进行网格划分，在ANSYS 2023 R1 目录树中，利用Mesh功能对从 Solidworks中导出的X-T文件模型进行四面体的网格单元划分，尺寸统一大小为5 mm(图6)，为了保证计算的精度达到分析的要求，对网格的类型和网格大小进行控制(表3)。③设置边界条件和载荷：为尽量施加膝关节生理情况的受力载荷，此次研究采用45%步态周期支撑相中足趾离地时的静力学分析，将股骨模型近端所有节点的 X、Y、Z 轴方向上的自由度均被约束为0 mm；然后予以施加胫股关节轴向压缩负荷2 100 N(内侧髁60%+外侧髁40%)，内外旋转扭矩为7 Nm，髌股关节面的压缩负荷为225 N，股骨髁部远端的剪切力为220 N[21] (图7)。

1.5.5 AFN深度为3 mm使用骨水泥填补模型的有限元分析根据上述实验步骤得出应力云图及应力值，发现当切迹深度由2 mm增加到3 mm时应力增大出现骤变，故在Solidworks软件中建立用骨水泥填补深度为3 mm的AFN模型，X-T格式进行保存输出导入至ANSYS 17.0软件进行有限元网格化，有限元网格划分后节点为23 907，单元数为12 616。予以材料赋值后使用上诉方法进行力学分析。

1.6 统计学分析 各组数据用统计表和统计图表示，不同骨质强度之间的比较和不同AFN深度之间的比较分别运用单元线性回归分析；而骨水泥填补AFN深度为3 mm、AFN深度为3 mm、未发生AFN之间运用单因素方差分析，其检验水准均为 α=0.05。所有数据处理都在SPSS 25.0软件上完成。文章统计学方法已经贵州医科大学、贵黔国际医院生物统计学专家审核。

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

讨论

3.1 结论 以上结果显示，任何骨质强度下，随着AFN深度加深，股骨髁上应力增大，并集中于切迹处。在正常骨质情况下，当AFN深度在3 mm和4 mm间有一个应力突变点；在骨质疏松情况下，当AFN深度在2 mm和3 mm间有一个应力突变点，再一次证明AFN与膝关节股骨假体周围骨折存在相关性。使用骨水泥填补AFN深度为3 mm的模型，股骨前髁截面应力下降。综上所述，全膝关节置换术中应该避免AFN出现，特别是骨质疏松患者；如术中出现AFN且深度超过骨皮质厚度时，可以采用骨水泥均匀填充AFN以降低股骨髁上应力，降低膝关节股骨假体周围骨折的发生率。
3.2 正常膝关节股骨置换有限元模型的建立 随着计算机技术的飞速发展，有限元分析被广泛用于各大领域；有限元技术在1972年应用于骨科生物力学后[22]，数字医学在骨科领域得到迅速的发展[23]。与传统生物建模方法相比，三维有限元模型不仅可模拟各种传统生物力学实验无法完成的实验研究，而且在建模速度和时间上有着明显的优势[24]。
使用有限元方法进行生物力学研究的关键步骤之一是建立稳定的、准确的、更接近实际情况的模型。对于骨科三维建模领域来说，对不同组织赋值至关重要，通过国内外的大量研究，对骨组织的赋值方法有均值赋值法、灰度赋值及皮-松质骨赋值法[25]；每种赋值方法各有优缺点，均值赋值法是为了简化实验，但与真实情况存在较大差异；灰度赋值及皮-松质骨赋值法更接近实际情况[26]，两者无明显差异，此次实验CT数据为正常骨质条件志愿者提供，故选择皮-松质骨赋值法，更符合真实生物力学特征。虽然三维有限元模型的建立是基于实体数据和真实情况，但构建好模型后需要进行模型有效性的验证，验证方法有与实体生物力学实验结果相比较与已验证的计算机模型做比较[27]；第一种方法成本极高且很难获取良好的实体标本，临床中大都使用第二种方法进行验证。此次实验模型与既往向相关文献研究结果相近[28-29]，可用于有限元的分析，对临床的应用具有一定的参考性和指导性意义。
3.3 不同骨质强度发生AFN对股骨髁上应力的影响 随着关节外科的迅速发展，全膝关节置换术能有效缓解晚期膝关节的疼痛，恢复患肢力线，而达到良好的膝关节功能，提高了患者生活质量，已经成为膝骨关节炎、类风湿关节炎、膝关节创伤性关节炎，以及其他一些非化脓性关节病后期的有效治疗措施，被广泛应用于临床[30]。但其仍存在众多并发症，诸如感染、血栓、假体周围骨折、慢性疼痛等并发症给患者及医生带来不小麻烦[31]；随着目前膝关节置换患者高龄化，膝关节置换术后假体周围骨折可导致较高的致死率。假体周围骨折中发生率最高的是膝关节股骨假体周围骨折，导致其发生的原因很多，相关研究表明骨质疏松、AFN是其中两个重要原因。
大量研究表明AFN将导致膝关节股骨假体周围骨折[7-8]，但仍有少数学者认为微小的股骨前皮质与股骨髁上骨折无关[9-10]，目前尚存在争议，此次实验研究表明任何骨质强度下，随着AFN深度的加深，股骨髁上最大等效应力越大；在正常骨质条件下，AFN深度在超过3 mm时；而骨质疏松情况下，当AFN深度超过2 mm，应力出现骤变，且出现应力集中，应力出现骤变后随深度的加深应力增大幅度较前减少。此次实验骨皮质厚度约为2.96 mm，出现应力骤变且集中的AFN深度即为骨皮质完全切除，骨松质完全暴露；这与既往研究的结果显示一致，AFN深度超过3 mm时，股骨髁上应力集中，很有可能导致膝关节股骨假体周围骨折[32-33]。这是因为正常情况下，骨单元作为骨皮质的基本单位，在受到外界生理载荷情况下在间质骨交界处可适当诱发显微裂纹和释放能量，避免应力集中[34]，在病理载荷下，骨单元可以利用最外层间骨板尽量阻止或偏转显微裂纹的迅速生长与扩展，防止宏观裂纹甚至应力性骨折的出现[35]。故当皮质骨完全切除时，髁上应力出现应力骤变且集中，许多研究表明骨质疏松与假体周围骨折相关联。骨质疏松以骨皮质变薄、骨小梁减少为特点，从而骨组织的强度及刚度降低，受到同样的外力更容易导致骨折。此次研究受同样的外力，不同骨质强度下膝关节股骨置换后股骨髁上的最大等效应力与骨质强度呈现负相关，在AFN深度超过骨皮质厚度时，随着骨质强度降低，股骨髁上应力明显增大，说明发生AFN时，骨质疏松患者更容易出现股骨髁上骨折。综上所述，骨质疏松患者行膝关节置换更容易出现假体周围骨折，围术期需行抗骨质疏松治疗，术前需作充分的术前计划，术中需仔细操作，避免AFN的出现，术后需个体化降低患者功能锻炼强度。而在正常骨质情况下，同样需避免AFN的发生。
3.4 膝关节置换术中骨缺损的处理 全膝关节置换术中的膝关节骨缺损是让手术医师头疼的问题，处理不好将明显提高关节翻修率。明确其分型是处理的前提，根据不同的基础参照有很多分型，临床中常用的分型为Stokley分型[36]，根据缺损是否累及周围皮质骨分为包容性和非包容性；另一种常用的临床分型Rand分型是根据缺损深度进行分型[37]，A型为缺损深度< 5 mm；B型为缺损深度为5-10 mm；C型为缺损> 10 mm；D型为空隙型骨缺损，周围骨皮质完整；E型为D型的基础上外周骨皮质不完整。其他的分型与上述两种分型差异不大，且临床中应用不多，在此不一一阐述。对于膝关节置换术中骨缺损的处理，临床中常用方法有金属垫块、自体骨移植、异体骨移植、骨水泥填充[38]。对于Rand A型缺损，使用骨水泥进行填补缺损是一种快速、简单、经济、稳定的处理方案[39]。临床中常用的骨水泥分为磷酸钙骨水泥和丙烯酸酯类骨水泥两种[40]，各有其优缺点。磷酸钙骨水泥生物组织相容性好，对周围软组织损伤较小，常用于抗生素和骨诱导因子的载体；但其机械强度不高，不能用于负重区域的骨缺损修复。丙烯酸酯类骨水泥虽然其生物相容性较差，但其由固、液双组分材料组成，具有及时塑形、机械强度良好等优点，且大量学者正在研究避其短板，发挥其优点[41]。
目前一些学者将聚甲基丙烯酸丁酯微球与脆性聚甲基丙烯酸甲酩共混制成低模量聚甲基丙烯酸甲酯骨水泥，将其性能进行了改良，临床中常用于人体密度高组织的置换或修复。此次实验结合上述骨缺损分型及处理方案，选择丙烯酸酯类骨水泥填补AFN，研究其最大应力的变化进而探索骨水泥的修复作用。当使用骨水泥填补AFN后股骨髁上应力减少，再一次说明骨水泥对修复生物体坚硬组织具有一定效果。
3.5 此次研究的不足之处及对未来的展望 此次实验通过构建不同骨质条件下AFN的膝关节股骨置换模型，探索了不同骨质强度下AFN深度对股骨髁上应力的影响，但骨组织是一种极为复杂的生物复合材料，故其生物力学是十分复杂的，骨组织周围组织对其力学特点也存在影响，因此以下方面是此次实验的不足之处及对今后的展望：①膝关节的组成复杂，周围软组织对其运动起到非常重要的作用，此次实验基于实验条件及技术，简化了实验模型，今后的研究需全面、更加真实地构建模型；②有限元分析法模型的单元划分、节点选择及材料赋值在一定程度上是被假定的，并不能真实代表人体骨组织，不能准确代表骨组织内不同区域的弹性模量变化；③此次实验进行了既往文献中常用的静力学分析，而膝关节的运动是极为复杂的，上下楼梯、平地行走、游泳、爬山时其所受的应力是不一样的，今后应探索膝关节各方面运动的动力学研究，更加真实、全面地体现膝关节的生物力学特点；④此次实验中AFN深度分组较少，数据样本不大，未进行AFN长度和横截面积的影响研究；今后应研究AFN的长度、横截面积对股骨髁上应力的影响；⑤皮质骨、松质骨是骨组织的宏观组成，此次研究基于宏观概念，不能准确体现骨组织的生物复杂性，今后应在骨组织的微观结构层面作生物力学研究。
中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

文章快阅

延伸阅读

全膝关节置换术(total knee arthroplasty,TKA)是膝关节骨性关节炎的终末期治疗手段，随着全球老龄化及膝关节假体的不断更新换代，其手术量也逐年上升；同时其带来的并发症也逐年上升。膝关节股骨假体周围骨折是TKA术后常见并发症之一，不仅增加了患者的创伤及住院费用，同时导致膝关节功能丢失，是临床医师和患者需高度重视的问题。相关研究表明导致其发生的两个主要危险因素为骨质强度和股骨前皮质切迹(Anterior femoral notching, AFN)；然而一些研究提示正常骨质强度患者发生AFN与股骨髁上假体周围骨折之间无相关性。骨质疏松患者TKA术中发生AFN研究尚少，而临床中行膝关节置换患者大多为骨质疏松患者，故骨质疏松膝关节置换术中发生AFN相关生物力学研究迫在眉睫。综上所述，目前骨质疏松患者TKA术中发生AFN研究尚缺，AFN与膝关节股骨假体周围骨折尚存争议，且目前对出现AFN情况的补救措施的研究尚缺，尤其是在对不同骨质强度发生AFN的“安全深度”缺乏足够的认识，从而可能增加膝关节股骨假体周围骨折的风险性。随着骨科数字化的迅速发展，三维有限元成为骨科相关力学研究中重要的工具之一。相关研究表明，有限元分析手段相比于实体力学研究存在以下优点：可重复性高，更简洁化，更准确，逐渐发展为骨科力学研究的主流工具，对骨科关节置换领域作出了巨大贡献。本研究旨在通过三维有限元方法研究不同骨质强度下TKA术中发生AFN的生物力学变化及使用骨水泥填补AFN后的髁上应力变化。为临床预防TKA术后发生膝关节股骨假体周围骨折和指导术中出现AFN的处理方案提供力学理论基础。

#br#

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

不同骨强度下全膝置换过程中发生股骨前皮质切迹的三维有限元分析

Three-dimensional finite element analysis of anterior femoral notching during total knee arthroplasty at different bone strengths

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 7

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	许浩, 丁璐, 李潇. 中央螺丝机械磨损对莫氏锥度种植系统力学影响的三维有限元分析[J]. 中国组织工程研究, 2025, 29(在线): 1-9.
[2]	周佳俊, 马飞, 冷叶波, 徐世财, 何宝强, 李洋, 廖烨晖, 唐强, 唐超, 王清, 钟德君. 全脊柱MRI评估骨质疏松患者椎体骨折的分布特点及临床意义[J]. 中国组织工程研究, 2025, 29(9): 1883-1889.
[3]	孙晓君, 王华明, 张德宏, 宋学文, 黄晋, 张辰, 裴生太. 有限元法在儿童发育性髋关节发育不良及治疗中的作用[J]. 中国组织工程研究, 2025, 29(9): 1897-1904.
[4]	赵济宇, 王少伟. 叉头框转录因子O1信号通路与骨代谢[J]. 中国组织工程研究, 2025, 29(9): 1923-1930.
[5]	李良奎, 黄永灿, 王鹏, 于滨生. 颈椎前路椎体骨化物可控前移融合对后纵韧带骨化物和内植物影响的有限元分析[J]. 中国组织工程研究, 2025, 29(9): 1761-1767.
[6]	徐彪, 路坦, 姜亚琼, 阴玉娇. 有限元分析不同程度冈上肌断裂对肩关节应力的影响[J]. 中国组织工程研究, 2025, 29(9): 1768-1774.
[7]	付恩洪, 杨行, 梁成, 张小刚, 张亚丽, 靳忠民. 基于OpenSim预测青少年跖屈肌无力的下肢生物力学行为[J]. 中国组织工程研究, 2025, 29(9): 1789-1795.
[8]	赵嘉诚, 任诗齐, 祝秦, 刘佳佳, 朱翔, 杨洋. 原发性骨质疏松潜在生物标志物的生物信息学分析[J]. 中国组织工程研究, 2025, 29(8): 1741-1750.
[9]	陈帅, 金杰, 韩化伟, 田宁晟, 李志伟. 两样本孟德尔随机化分析循环炎症细胞因子与骨密度的因果关联[J]. 中国组织工程研究, 2025, 29(8): 1556-1564.
[10]	章镇宇, 梁秋健, 杨军, 韦相宇, 蒋捷, 黄林科, 谭桢. 新橙皮苷治疗骨质疏松症的靶点及对骨髓间充质干细胞成骨分化的作用[J]. 中国组织工程研究, 2025, 29(7): 1437-1447.
[11]	芦劼明, 李亚静, 杜培洁, 徐冬青. 人工和天然草地对青年女性跳跃落地时下肢生物力学表现的影响 [J]. 中国组织工程研究, 2025, 29(6): 1101-1107.
[12]	张帅, 李子春, 徐亦豪, 谢晓峰, 郭忠圣, 赵清扬. 经颅磁声电刺激强度对小鼠前额叶皮质网络可塑性的影响[J]. 中国组织工程研究, 2025, 29(6): 1108-1117.
[13]	李岳尧, 张民, 杨家驹. 肉苁蓉苷A通过JNK/MAPK通路抑制破骨细胞活性[J]. 中国组织工程研究, 2025, 29(6): 1144-1151.
[14]	郑琳, 靳文君, 罗珊珊, 黄芮, 王杰, 程余婷, 安哲庆, 熊玥, 巩仔鹏, 廖健. 杜仲促进去势大鼠牙槽骨成骨的作用[J]. 中国组织工程研究, 2025, 29(6): 1159-1167.
[15]	黄小彬, 葛继荣, 李生强, 谢丽华, 黄景文, 何艳艳, 薛立鹏. 不同滋阴补肾法干预去势大鼠破骨通路的机制[J]. 中国组织工程研究, 2025, 29(6): 1214-1219.