股骨颈骨折空心钉内固定后股骨近端骨质疏松的有限元分析

doi:10.12307/2023.909

中国组织工程研究 ›› 2024, Vol. 28 ›› Issue (6): 862-867.doi: 10.12307/2023.909

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

股骨颈骨折空心钉内固定后股骨近端骨质疏松的有限元分析

薛晓峰1，魏永康2，乔晓红1，杜玉勇1，牛建军1，任立新1，杨慧峰1，张治民1，郭媛2，陈维毅2

1山西省吕梁市人民医院，山西省吕梁市 033000；2太原理工大学生物医学工程学院，山西省太原市 030000

收稿日期:2022-11-23 接受日期:2023-01-29 出版日期:2024-02-28 发布日期:2023-07-11
通讯作者: 乔晓红，博士，主任医师，山西省吕梁市人民医院，山西省吕梁市 033000
作者简介:薛晓峰，男，1979年生，山西省临县人，汉族，2008年山西医科大学毕业，硕士，副主任医师，主要从事骨科创伤、关节和运动医学研究。魏永康，男，1995年生，河南省商丘市人，汉族，太原理工大学在读硕士，主要从事肌骨系统生物力学研究。
基金资助:
山西省科学技术厅基金资助(201903D321150)，项目负责人：乔晓红

Finite element analysis of osteoporosis in proximal femur after cannulated screw fixation for femoral neck fracture

Xue Xiaofeng1, Wei Yongkang2, Qiao Xiaohong1, Du Yuyong1, Niu Jianjun1, Ren Lixin1, Yang Huifeng1, Zhang Zhimin1, Guo Yuan2, Chen Weiyi2

1Lvliang People’s Hospital, Lvliang 033000, Shanxi Province, China; 2College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030000, Shanxi Province, China

Received:2022-11-23 Accepted:2023-01-29 Online:2024-02-28 Published:2023-07-11
Contact: Qiao Xiaohong, MD, Chief physician, Lvliang People’s Hospital, Lvliang 033000, Shanxi Province, China
About author:Xue Xiaofeng, Master, Associate chief physician, Lvliang People’s Hospital, Lvliang 033000, Shanxi Province, China Wei Yongkang, Master candidate, College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030000, Shanxi Province, China
Supported by:
a grant from Shanxi Provincial Department of Science and Technology, No. 201903D321150 (to QXH)

摘要/Abstract

摘要：

文题释义：

股骨颈骨折空心钉内固定术：是指采用3枚空心钉呈正三角形或倒三角形排列、平行分散内固定治疗股骨颈骨折的一种手术方式，能够使股骨颈骨折端保持相对稳定。空心钉固定治疗股骨颈骨折的优点有创伤小、可经过皮肤直接操作、固定效果较好、手术操作难度较小、空心结构能够显著降低关节腔内压力，并具有动力加压作用，有利于骨折顺利愈合。
有限元法：是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解，每个子区域都成为简单的部分，这种简单部分就称作有限元。该方法在经过近几十年的高速发展已经较为成熟，被广泛应用于肌骨系统生物力学相关分析中。

背景：股骨颈骨折空心钉内固定术后由于患肢短期内常常不能负重，且高刚度的内植物对骨折断端存在应力遮挡效应，易导致患肢出现骨质疏松，股骨近端生物力学分布发生变化，术后股骨头坏死发病率较高，目前关于股骨颈骨折术后股骨近端骨质疏松对股骨近端及空心钉生物力学影响的研究较少。

目的：通过有限元分析探讨股骨颈骨折术后发生骨质疏松对空心钉内固定治疗的生物力学影响，探究生物力学因素在股骨头坏死进程中的作用。
方法：获取1例股骨颈骨折患者股骨CT扫描数据，利用Mimics 19.0、3-Matic、UG 11.0、Hypermesh 14.0、Abaqus软件建立空心钉治疗股骨颈骨折的股骨近端模型，利用Abaqus软件分析1种术后股骨近端无骨质疏松、3种术后股骨近端骨质疏松的有限元模型，测量分析4种模型不同部件的应力、接触压力、位移峰值及云图，对比分析股骨头内部应力变化及分布情况。

结果与结论：股骨头及下前空心钉的应力、接触压力随骨质疏松程度变化较大，4种模型的位移峰值随着骨质疏松程度加重而缓慢增长。通过单因素方差分析，结果显示骨质疏松程度对不同部件的应力、接触压力、位移峰值无显著性影响，股骨头内部应力分布随骨质疏松发生相应变化，股骨近端的生物力学环境变化对股骨头坏死有着重要影响。

https://orcid.org/0000-0002-8344-6715 (薛晓峰)

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

关键词: 股骨, 股骨颈骨折, 空心钉, 内固定术, 骨质疏松, 股骨头坏死, 生物力学, 有限元

Abstract: BACKGROUND: After the internal fixation of cannulated screws in femoral neck fractures, because the affected limb is often unable to bear weight in the short term and the implants with high stiffness have a stress shielding effect on the fracture end, it is easy to cause osteoporosis of the affected limb and changes in the biomechanical distribution of the proximal femur, the incidence of osteonecrosis of the femoral head is high after surgery. At present, few studies have been conducted on the biomechanical effects of osteoporosis at the proximal end of the femur occurring after femoral neck fracture surgery on femoral neck fracture treated with cannulated screws.
OBJECTIVE: Using finite element analysis, to investigate the biomechanical effects of osteoporosis occurring after femoral neck fracture surgery on femoral neck fracture treated with cannulated screws and explore the role of biomechanical factors in osteonecrosis of the femoral head.
METHODS: Based on the obtained CT scan data of the femur in a patient with a femoral neck fracture, a proximal femoral model for internal fixation for femoral neck fracture was established by Mimics 19.0, 3-Matic, UG 11.0, Hypermesh 14.0, and Abaqus software. One finite element model of the proximal femur without osteoporosis and three finite element models of the proximal femur with osteoporosis were analyzed using Abaqus software. The stress, contact pressure, displacement peak and cloud map under different components of the four models were measured and analyzed, and the internal stress changes and distribution of the femoral head were compared and analyzed.

RESULTS AND CONCLUSION: The stresses and contact pressures of the femoral head and lower anterior cannulated screws varied more with the degree of osteoporosis. The peak displacement of the four models increased slowly with the degree of osteoporosis. By one-way analysis of variance, there was no significant effect of the degree of osteoporosis on the peak stress, contact pressure, and displacement of the different components. The internal stress distribution of the femoral head changed with the degree of osteoporosis. Changes in the biomechanical environment of the proximal femur have an important impact on osteonecrosis of the femoral head.

Key words: femur, femoral neck fracture, cannulated screw, internal fixation, osteoporosis, femoral head necrosis, biomechanics, finite element

中图分类号:

薛晓峰, 魏永康, 乔晓红, 杜玉勇, 牛建军, 任立新, 杨慧峰, 张治民, 郭媛, 陈维毅. 股骨颈骨折空心钉内固定后股骨近端骨质疏松的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 862-867.

Xue Xiaofeng, Wei Yongkang, Qiao Xiaohong, Du Yuyong, Niu Jianjun, Ren Lixin, Yang Huifeng, Zhang Zhimin, Guo Yuan, Chen Weiyi. Finite element analysis of osteoporosis in proximal femur after cannulated screw fixation for femoral neck fracture[J]. Chinese Journal of Tissue Engineering Research, 2024, 28(6): 862-867.

图/表（结果） 6

2.1 模型验证将未发生骨质疏松的股骨近端模型进行有限元分析，其具体工况：材料属性、边界条件、载荷等设定同上所述。计算结果表明：未发生骨质疏松的股骨近端模型的应力峰值为43.8 MPa，应力由股骨颈下缘经股骨内侧向股骨干传递。对比TAYLOR 等[26]的研究结果，该研究模型股骨近端的应力分布与其模型股骨近端应力分布基本相同，故该研究所建立的有限元模型有效，见图4。

2.2 不同模型计算结果的对比分析 4种不同模型的有限元计算结果(1种股骨颈骨折空心钉内固定术后股骨近端无骨质疏松模型、3种股骨颈骨折空心钉内固定术后股骨近端骨质疏松模型)，其应力都主要分布在空心钉附近，当不同刚度的材料共同承担外部载荷时，刚度大的材料将会承载较高的载荷，刚度低的材料承担较小的载荷，该应力结果符合应力遮挡效应。
股骨近端无骨质疏松模型与骨质疏松模型的等效应力分析：随着股骨近端骨质疏松程度不断加重，其股骨头上所承担的载荷越来越小，应力峰值从267.9 MPa降低到78.9 MPa。而相对于下前空心钉，当股骨近端发生轻微骨质疏松时，其应力峰值相对股骨近端无骨质疏松模型略微减小，随着股骨近端骨质疏松程度不断加重，其等效应力峰值发生较大的增高，从120.6 MPa增高到189.0 MPa，而4种有限元模型的其他部件：上空心钉、下后空心钉、股骨干的应力峰值变化不太明显。4种模型的应力云图，见图5和图6，不同部件的应力峰值见表3。

股骨近端无骨质疏松模型与骨质疏松模型的接触压力分析：随着股骨近端骨质疏松程度不断加重，其股骨头及下前空心钉上的接触压力峰值越来越小，股骨头接触压力峰值从379.3 MPa降低至114.6 MPa，下前空心钉接触压力峰值从141.0 MPa降低至88.9 MPa。4种有限元模型的其他部件：上空心钉、下后空心钉、股骨干接触压力峰值仅有细微的变化，见表4。随着股骨近端骨质疏松程度的不断加重，其最大位移也在不断增大。股骨颈骨折空心钉内固定术后股骨近端无骨质疏松模型与轻微骨质疏松、中度骨质疏松、重度骨质疏松模型的最大位移发生在股骨头处，分别为0.80，0.81，0.83，0.88 mm。3枚空心钉的最大位移也略有增长，而股骨干的位移基本没有变化。利用单因素方差分析法对数据进行处理分析，在显著性水平α=0.05条件下，骨质疏松程度对不同部件的应力、接触压力、位移峰值无显著性影响，其F值分别为F=0.03，F=0.12，F=0.14均小于F0.05(3，16)=3.24。

该研究还重点分析了股骨头内部的应力分布情况，探究了股骨近端生物力学环境变化对股骨头坏死的影响。随着股骨近端骨质疏松程度的加重，股骨头内部应力分布也发生相应变化。股骨近端无骨质疏松模型在股骨头上外侧区域附近的应力峰值为13.8 MPa，轻微骨质疏松时，该区域应力峰值降低为10.8 MPa，随着骨质疏松程度的进一步加重，应力峰值位置发生改变，转移至股骨前上内侧附近，中度骨质疏松时该区域应力峰值为7.9 MPa，重度骨质疏松时该区域应力峰值增长至8.7 MPa，并与临床实际坏死区做对比，该病例的CT片示：股骨头坏死区域位于股骨头的前上内侧，与该研究中的有限元模型发生骨质疏松后股骨头内部应力峰值转移位置基本相符，见图7。研究表明，股骨颈骨折空心钉内固定术后股骨近端发生骨质疏松，会改变股骨头内部应力分布，应力峰值位置的转移可能是该位置坏死的生物力学因素之一。

参考文献

[1] SWART E, ROULETTE P, LEAS D, et al. ORIF or Arthroplasty for Displaced Femoral Neck Fractures in Patients Younger Than 65 Years Old: An Economic Decision Analysis. J Bone Joint Surg Am. 2017;99(1):65-75.
[2] TIAN FM, ZHANG L, ZHAO HY, et al. An increase in the incidence of hip fractures in Tangshan, China. Osteoporos Int. 2014;25(4):1321-1325.
[3] 张长青,张英泽,余斌,等.成人股骨颈骨折诊治指南[J].中华创伤骨科杂志,2018,20(11):921-928.
[4] 徐凯航,纪方.青壮年股骨颈骨折的治疗进展[J].中华创伤骨科杂志,2020,22(6):549-552.
[5] SLOBOGEAN GP, SPRAGUE SA, SCOTT T, et al. Complications following young femoral neck fractures. Injury. 2015;46(3):484-491.
[6] STOCKTON DJ, O’HARA LM, O’HARA NN, et al. High rate of reoperation and conversion to total hip arthroplasty after internal fixation of young femoral neck fractures: a population-based study of 796 patients. Acta Orthop. 2019;90(1):21-25.
[7] PATTERSON JT, ISHII K, TORNETTA P 3RD, et al. Open Reduction Is Associated With Greater Hazard of Early Reoperation After Internal Fixation of Displaced Femoral Neck Fractures in Adults 18-65 Years. J Orthop Trauma. 2020;34(6):294-301.
[8] SLOBOGEAN GP, STOCKTON DJ, ZENG B, et al. Femoral Neck Fractures in Adults Treated With Internal Fixation: A Prospective Multicenter Chinese Cohort. J Am Acad Orthop Surg. 2017;25(4):297-303.
[9] LI M, COLE PA. Anatomical considerations in adult femoral neck fractures: how anatomy influences the treatment issues? Injury. 2015; 46(3):453-458.
[10] UEO T, TSUTSUMI S, YAMAMURO T, et al. Biomechanical aspects of the development of aseptic necrosis of the femoral head. Arch Orthop Trauma Surg (1978). 1985;104(3):145-149.
[11] 张永飞,张义修.股骨颈骨折术后股骨头坏死的力学因素[J].骨与关节损伤杂志,2001,16(4):270-272.
[12] KIM YM, LEE SH, LEE FY, et al. Morphologic and biomechanical study of avascular necrosis of the femoral head. Orthopedics. 1991;14(10): 1111-1116.
[13] WANG Y, MA JX, YIN T, et al. Correlation Between Reduction Quality of Femoral Neck Fracture and Femoral Head Necrosis Based on Biomechanics. Orthop Surg. 2019;11(2):318-324.
[14] SUN W, WANG BL, LI ZR, et al. Chinese specialist consensus on diagnosis and treatment of osteonecrosis of the femoral head. Orthop Surg. 2011;3(2):131-137.
[15] YUE Y, YANG H, LI Y, et al. Combining ultrasonic and computed tomography scanning to characterize mechanical properties of cancellous bone in necrotic human femoral heads. Med Eng Phys. 2019;66:12-17.
[16] MA JX, HE WW, ZHAO J, et al. Bone Microarchitecture and Biomechanics of the Necrotic Femoral Head. Sci Rep. 2017;7(1):13345.
[17] CIARELLI MJ, GOLDSTEIN SA, KUHN JL, et al. Evaluation of orthogonal mechanical properties and density of human trabecular bone from the major metaphyseal regions with materials testing and computed tomography. J Orthop Res. 1991;9(5):674-682.
[18] HARP JH, ARONSON J, HOLLIS M. Noninvasive determination of bone stiffness for distraction osteogenesis by quantitative computed tomography scans. Clin Orthop Relat Res. 1994;(301):42-48.
[19] HOBATHO MC, RHO JY, ASHMAN RB. Anatomical variation of human cancellous bone mechanical properties in vitro. Stud Health Technol Inform. 1997;40:157-173.
[20] 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.
[21] 张国栋,廖维靖,陶圣祥,等.股骨有限元分析赋材料属性的方法[J].中国组织工程研究与临床康复,2009,13(43):8436-8441.
[22] LENGSFELD M, SCHMITT J, ALTER P, et al. Comparison of geometry-based and CT voxel-based finite element modelling and experimental validation. Med Eng Phys. 1998;20(7):515-522.
[23] MEI J, LIU S, JIA G, et al. Finite element analysis of the effect of cannulated screw placement and drilling frequency on femoral neck fracture fixation. Injury. 2014;45(12):2045-2050.
[24] ZENG W, LIU Y, HOU X. Biomechanical evaluation of internal fixation implants for femoral neck fractures: A comparative finite element analysis. Comput Methods Programs Biomed. 2020;196:105714.
[25] EBERLE S, GERBER C, VON OLDENBURG G, et al. A biomechanical evaluation of orthopaedic implants for hip fractures by finite element analysis and in-vitro tests. Proc Inst Mech Eng H. 2010;224(10):1141-1152.
[26] TAYLOR ME, TANNER KE, FREEMAN MA, et al. Stress and strain distribution within the intact femur: compression or bending? Med Eng Phys. 1996;18(2):122-131.
[27] TIANYE L, PENG Y, JINGLI X, et al. Finite element analysis of different internal fixation methods for the treatment of Pauwels type III femoral neck fracture. Biomed Pharmacother. 2019;112:108658.
[28] LI J, ZHAO Z, YIN P, et al. Comparison of three different internal fixation implants in treatment of femoral neck fracture-a finite element analysis. J Orthop Surg Res. 2019;14(1):76.
[29] FREITAS A, TOLEDO JÚNIOR JV, FERREIRA DOS SANTOS A, et al. Biomechanical study of different internal fixations in Pauwels type III femoral neck fracture - A finite elements analysis. J Clin Orthop Trauma. 2020 ;14:145-150.
[30] 赵耀,许新忠,程文丹,等.股骨颈动力交叉钉系统内固定治疗股骨颈骨折的早期疗效分析[J].中国骨与关节损伤杂志,2022,37(4): 389-391.
[31] 孙玉生.不同数量空心螺钉内固定治疗股骨颈骨折的有限元分析及生物力学测定[D].蚌埠:蚌埠医学院,2020.
[32] JIANG D, ZHAN S, WANG L, et al. Biomechanical comparison of five cannulated screw fixation strategies for young vertical femoral neck fractures. J Orthop Res. 2021;39(8):1669-1680.
[33] 梁跃闯,李艳军,于鹤童,等.正三角构型, 倒三角构型空心螺钉内固定及股骨颈动力交叉钉系统固定治疗股骨颈骨折的有效性分析[J].创伤外科杂志,2022,24(9):659-664.
[34] 刘俊俊,付建国,柳威,等.三角稳定固定系统对股骨颈骨折生物力学特性的有限元分析[J].中国医学物理学杂志,2021,38(7):893-897.
[35] BHANDARI M, TORNETTA P 3RD, HANSON B, et al. Optimal internal fixation for femoral neck fractures: multiple screws or sliding hip screws? J Orthop Trauma. 2009;23(6):403-407.
[36] TSAO AK, ROBERSON JR, CHRISTIE MJ, et al. Biomechanical and clinical evaluations of a porous tantalum implant for the treatment of early-stage osteonecrosis. J Bone Joint Surg Am. 2005;87 Suppl 2:22-27.
[37] 邱士超,刘利,吴迪.股骨头坏死的生物力学及有限元分析相关研究[J].临床医药文献电子杂志,2019,6(28):43-44.
[38] XIAO D, YE M, LI X, et al. Biomechanics analysis for the treatment of ischemic necrosis of the femoral head by using an interior supporting system. Int J Clin Exp Med. 2015;8(3):4551-4556.
[39] ZHANG Y, TIAN L, YAN Y, et al. Biomechanical evaluation of the expansive cannulated screw for fixation of femoral neck fractures. Injury. 2011;42(11):1372-1376.
[40] 赵春伶,贾少薇,李剑,等.基于 3D 打印多孔支架和植入体的结构设计研究进展[J].医用生物力学,2019,34(4):446-452.
[41] 褚楷,张兴琳,鲁兴,等.股骨颈骨折内固定物取出后股骨头内部微骨折发生风险研究[J].中国修复重建外科杂志,2020,34(9):1091-1095.
[42] 潘显明,胡修德,谭映军,等.82 例青壮年股骨颈骨折治疗结果的评价[J].中华创伤杂志,2000,16(3):145-147.
[43] 王荣升,韩雪,王静,等.股骨颈骨折行闭合复位空心钉内固定后股骨头坏死[J].中国老年学杂志,2020,40(23):4984-4986.
[44] 马剑雄,何伟伟,赵杰,等.股骨头坏死发病机制研究的最新进展[J].中国组织工程研究,2017,21(27):4397.
[45] SEBESTYÉN A, MESTER S, VOKÓ Z, et al. Wintertime surgery increases the risk of conversion to hip arthroplasty after internal fixation of femoral neck fracture. Osteoporos Int. 2015;26(3):1109-1117.
[46] 王义生,李劲峰.股骨头坏死发病机制的研究现状与展望[J].中华实验外科杂志,2020,37(6):1001-1010.
[47] 申洪全,陆慧,张孝华,等.计算机辅助 2 枚空心钉置入内固定治疗股骨颈骨折的生物力学分析[J].中国组织工程研究,2020,24(6): 882-887.
[48] 余霄,张迪峰,宋蒙胜,等.骨质疏松性股骨颈骨折空心钉内固定术后股骨颈短缩对髋关节生物力学影响的有限元分析[J].中华医学杂志,2020,100(33):2628-2632.
[49] PEDERSEN DR, BRAND RA, DAVY DT. Pelvic muscle and acetabular contact forces during gait. J Biomech. 1997;30(9):959-965.

引言

股骨颈骨折是临床常见的创伤性骨折，约为全身骨折发病率的3.61%[1]，占髋部骨折的54%[2]。临床上，年轻患者(年龄< 65岁被定义为“年轻患者”[3])股骨颈骨折多由高能量创伤所致，对其治疗仍是一项挑战。目前，临床上对于青壮年股骨颈骨折患者，强调早期手术、解剖复位，尽可能早期恢复股骨头血运循环[4]，减少坏死率。虽然，闭合复位空心钉内固定技术的应用将股骨颈骨折愈合率提升到90%以
上[5]，极大改善了患者预后，但是青壮年股骨颈骨折术后股骨头坏死率至今仍无明显改善，股骨头坏死的发生率仍然高达40%-80%[6-7]，其临床预后结果仍不理想，尤其是股骨头坏死率较老年人高[8]。因此，青壮年股骨颈骨折被称为“未解决骨折”[9]。关于股骨颈骨折术后股骨头坏死的机制尚不明确，有研究表明部分患者在股骨颈周围血运良好的状况下，其仍发生了股骨头坏死[10]，因此一些学者开始意识到生物力学因素在股骨头坏死进程中的作用。股骨头坏死多发生于骨骼的塑形期内，股骨内部应力改变对内固定术后股骨头坏死发生及坏死进程的发展有着重要意义，股骨头坏死不仅是血供导致的，而且生物力学在其中也起到重要作用，故坏死的发生是生物学与生物力学共同作用的结果[11-12]。有研究表明股骨颈骨折复位质量差是股骨头坏死的危险因素，随股骨头空间位置的变化，股骨头受力方向与骨小梁方向不一致，导致股骨头内部应力分布发生变化，长期的不良刺激引发骨的重建与吸收，进而改变骨小梁结构[13]，发生股骨头坏死[14]。坏死硬化带区与非硬化带区密度弹性模量关系不同[15]，股骨头坏死区生物力学性能明显减弱，而硬化区生物力学性能增强，说明股骨头骨结构与功能的不协调可能在股骨头坏死的发病机制中起到了一定作用[16]。
股骨颈骨折空心钉内固定术后患肢短期内常常不能负重，导致患肢出现骨质疏松，有研究指出50岁以上的正常成年人，如果卧床超过1周，其可能会丢失全身骨量的1%，而三四个月后可丢失全身骨量的20%，而股骨颈骨折患者内固定术后往往需要不负重超过3个月，甚至半年，因此患肢骨量明显较健侧下降。
该研究拟通过股骨颈骨折术后股骨近端发生骨质疏松的有限元分析，来探讨股骨颈骨折术后骨质疏松对空心钉内固定治疗的生物力学影响，探究生物力学因素在股骨头坏死进程中的作用。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

材料方法

1.1 设计基于CT扫描数据的有限元建模及有限元分析。
1.2 时间及地点 2021年4月至2022年9月在太原理工大学生物医学工程学院完成。
1.3 材料
1.3.1 实验设备及软件东芝320排CT Aquilion ONE (日本Toshiba公司)；Mimics 19.0及3-Matic 19.0软件(比利时Materialise公司)；UG NX 11.0软件(德国Siemens PLM Software公司)；Hypermesh 14.0软件(美国Altair公司)；Abaqus 17.0软件(美国Dassault公司)。
1.3.2 志愿者资料选取1名64岁，身高170 cm，体质量67 kg，左侧股骨颈骨折Pauwels角为47°的男性志愿者，骨质良好，术前无其他疾病，术后6年患侧股骨头前上内侧发生坏死，已详细告知相关实验内容及实验目的，志愿者知情并签署知情同意书。该研究方案的实施符合《赫尔辛基宣言》和山西省吕梁市人民医院及太原理工大学对研究的相关伦理要求(医院伦理批件号：2020LYLL-03)。
1.4 实验方法

1.4.1 空心钉内固定股骨颈骨折的几何建模将志愿者股骨CT扫描数据以Dicom格式导入Mimics 19.0软件中，利用该软件的阈值分割、空腔填充、动态区域增长等命令构建出股骨近端的三维模型，再将该模型导入3-Matic 19.0软件中对模型进行去噪、平滑、光滑处理，并对股骨颈中心面处进行分割，以建立股骨颈骨折Pauwels角为50°的骨折模型。在UG 11.0软件中建立空心钉模型(直径4 mm，内部空心圆柱直径2 mm，其中2枚空心钉长为85 mm、1枚空心钉长为90 mm，空心钉螺纹部分简化为直径为6.3 mm的圆柱以简化计算)。将处理后的股骨颈骨折模型导入UG 11.0软件中，按实际手术置钉位置(正三角形内固定)，完成股骨颈骨折模型与空心钉模型的装配，得到空心钉内固定股骨颈骨折的几何模型，见图1。

1.4.2 有限元分析
(1)网格划分：在Hypermesh 14.0软件中对空心钉内固定骨折模型进行网格划分，生成四面体网格。有限元计算结果的精度受网格划分质量的影响较大，网格尺寸越小则生成的单元数目及节点数越多，有限元计算结果的精度越高。但随着单元数及节点数的增高，有限元计算量也会相应增大。在不影响有限元计算结果可靠性的前提下，该研究将股骨近端模型网格大小设置为1 mm，空心钉模型网格大小为0.5 mm，见图2。股骨头、股骨干、上空心钉、下后空心钉、下前空心钉的节点与单元数，见表1。

(2)定义材料属性：该研究采用的是患者术后(6年)复诊的CT数据，为了避免空心钉及坏死区对骨结构材料赋值的影响，故将患肢赋予了对侧健康股骨的材料属性，并定义为无骨质疏松模型。材料属性的赋值分为2个部分：①骨结构的材料赋值：采用Mimics提供的经验公式对股骨模型进行赋值[17-20]：
ρ=-13.4+1 017·Gv
E=-388.8+5 925·ρ

ρ(kg/m3)为骨的密度，E(Pa)为骨的弹性模量。在对骨结构材料赋值时，骨材料属性分类越多，骨的材料赋值越接近实际生理状态，有限元计算结果越接近真实情况，但相应的有限元计算量也会显著增加，该研究将股骨模型的材料属性分为10种[21]，泊松比为0.3[22]，此种材料属性分类既能满足有限元计算分析精度的需求，又能减少相关计算量。②内植物的材料赋值：空心钉的弹性模量为206 000 MPa，泊松比为0.3[23]。股骨、空心钉模型材料属性均为各向同性，见图3。为模拟3种股骨颈骨折空心钉内固定术后股骨近端发生骨质疏松，该研究将这3种模型股骨近端的弹性模量分别设为无骨质疏松股骨近端弹性模量的76%(轻微骨质疏松)、52%(中度骨质疏松)、28%(重度骨质疏松)，以模拟股骨近端骨质疏松程度不断加重的渐变过程。

(3)相互作用：股骨颈的骨折面完全断裂，设置为相互接触，摩擦系数为0.46[24]。空心钉具有螺纹的前端部分与股骨头接触面设置为绑定关系，空心钉的光滑体部分与股骨接触面设为接触，摩擦系数为0.3[25]。

(4)边界条件及载荷：定义模型底面及向上3 cm内所有节点为完全固定。该研究探究的工况为患者单足站立，故参考TAYLOR等[26]的研究，单腿站立时的髋关节反作用力大小约为体质量的4.28倍，作用力方向与矢状面呈20°夹角，与冠状面呈3°夹角，作用在股骨头与髋臼窝的实际接触面上。外展肌力作用方向沿肌肉收缩方向与矢状面呈20°夹角。髂腰肌力作用方向沿肌肉收缩方向与矢状面呈9°夹角，与冠状面呈47°夹角。髂胫束力作用方向为沿重力方向竖直向下。外展肌、髂胫束和髂腰肌3个主要肌群作用力的大小按体质量倍数关系计算得出，具体数值见表2。

1.5 主要观察指标该研究进行了4种不同模型的有限元分析，1种股骨颈骨折空心钉内固定术后股骨近端无骨质疏松模型、3种股骨颈骨折空心钉内固定术后股骨近端骨质疏松模型。主要对比分析了不同模型下各部件的应力、接触压力、位移等参数。
1.6 统计学分析采用SPSS 21.0 统计软件对实验数据进行了分析整理，利用单因素方差分析研究骨质疏松程度对股骨及空心钉的应力峰值、接触压力、位移的影响，显著性水平α=0.05，当F值≥F0.05(3，16)=3.24时，认为骨质疏松程度对计算结果有显著性影响。

讨论

有限元分析已广泛应用于肌骨系统生物力学中，大量学者对股骨颈骨折内固定进行了相关研究，探究不同手术类型内固定：空心钉联合钢板内固定，加压螺钉组合防旋螺钉内固定，空心钉内固定治疗股骨颈骨折等[27-30]。有研究表明不同空心钉数量[31]、位置结构[32-34]、角度及不同加压种类螺钉对内固定均有影响[35]，其中空心钉内固定术后股骨骨质疏松及股骨头坏死的相关研究一直是研究热点。有学者通过有限元探究了骨质疏松性股骨颈骨折术后股骨颈缩短对骨折愈合的影响，骨质疏松患者更易发生股骨颈缩短，随着股骨颈的缩短，股骨及内固定物的应力分布发生改变[32]。目前股骨头坏死的机制尚无明确定论，国内外针对此发病机制提出了很多说法：酗酒、使用类固醇药物、创伤、血液疾病[36]、骨内压力值增高、血运障碍、脂类代谢出现紊乱、免疫机制紊乱等因素，在股骨头坏死的治疗过程中，生物力学因素有着重要地位[37]。对于股骨头坏死的治疗，有研究从生物力学角度发现多孔钽植入物可以在骨缺损的情况下为软骨下骨提供较好支撑，且多孔钽与骨植入物在股骨头内的应力和应变相似，内部的支撑系统有助于防止股骨头塌陷[38-39]。而内固定术后高弹性模量的内植物会产生应力遮挡，降低对骨折部位的应力刺激，影响骨重塑，可能导致骨折不愈合、骨吸收、骨质疏松及假体松动等问题[39-40]，而且股骨头坏死更易出现在骨的重塑期。有学者发现取出内固定物后，骨密度对股骨头内部出现微骨折有着密切的相关性[41]。潘显明等[42]认为股骨颈骨折术后股骨头发生坏死与患者开始负重的时间密切相关，应避免术后10个月内过早的负重，以减少股骨头坏死的发生。患者术后康复存在各类争议，沃尔夫定律指出外部压力对骨骼的生长起到调节作用，合理的外力刺激会促进骨骼生长，青壮年患者骨折后因长时间不能负重，其骨小梁结构易发生变化，骨小梁变细、变脆，进而影响股骨近端生物力学分布，导致发生骨质疏松及股骨头坏死，且年龄越大，越容易出现骨质疏松，同时由于其股骨颈周围的微循环较差，发生股骨头坏死的概率越大[43]。股骨头坏死的发病机制至今仍是一大难题，股骨颈骨折患者发生股骨头坏死多出现在负重生活后的1-3年内，此时大多骨折部位已经愈合，生物力学因素可能对发生股骨头坏死起着更为重要的影响[44]。但目前普遍认为，术后均应当避免早期的完全负重和功能锻炼[45]，以保证骨折愈合阶段股骨近端稳定的生物力学环境，为骨愈合及血供重建提供稳定力学环境。
该研究通过股骨颈骨折术后股骨近端发生骨质疏松的有限元分析，来探讨股骨颈骨折术后发生骨质疏松对空心钉治疗股骨颈骨折的生物力学影响，探究生物力学因素对股骨头坏死的作用。对比分析了青壮年股骨颈骨折空心钉内固定术后早期与术后远期股骨近端发生不同程度骨质疏松后的应力、接触压力、位移等数据，结果表明随着骨质疏松程度的加重，股骨近端模型内部的应力均发生了相应的变化，在股骨头上外侧区域的应力峰值转移至股骨下内侧区域，这可能是股骨头近端坏死的生物力学因素之一。骨质疏松及负重学说都表明股骨头的负重过载可能引发股骨头坏死，股骨头内部坏死区与健康骨区负载不一致时会加剧股骨头的坏死[46]。空心钉内固定股骨颈骨折术后约60%的患者会在5年内发生股骨颈缩短[47]，术后股骨颈缩短也会影响髋关节内部生物力学环境，股骨颈缩短在股骨头坏死发生中也可能会起到一定作用[48]。
该研究也存在局限性，首先髋关节具有复杂的生理结构，简化的空心钉内固定骨折模型中，并未建立相应的髋臼、软骨、肌肉、韧带等组织；另外，一个完整的步态周期中，不同阶段髋关节内部接触力是不同的[49]，单一的边界条件及载荷仅能模拟患者单足站时股骨的受力分析。作者认为，该研究虽存在一些不足，但其研究成果仍能为临床治疗提供一定的参考依据。针对上述局限性，后续将建立包含髋臼、软骨、肌肉、韧带等组织的髋关节模型，并建立更多的工况以模拟完整步态周期不同阶段，通过对比分析，更深层次研究骨质疏松对内固定术后股骨近端的生物力学影响。该研究探讨了术后股骨头发生坏死的可能因素，并指出当患者术后发生骨质疏松时，应尽早介入治疗，以避免可能发生的股骨头坏死，为实际临床治疗提供生物力学依据。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

股骨颈骨折空心钉内固定后股骨近端骨质疏松的有限元分析

Finite element analysis of osteoporosis in proximal femur after cannulated screw fixation for femoral neck fracture

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 6

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	王强, 李世芸, 熊鹰, 李甜甜. 不同颈椎前路椎间融合系统内固定颈椎的生物力学变化[J]. 中国组织工程研究, 2024, 28(6): 821-826.
[2]	魏源标, 林展, 陈燕梅, 杨腾辉, 赵晓, 陈扬声, 周燕辉, 杨敏超, 黄飞麒. 颈椎矢状位端提手法对椎间盘及小关节影响的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 827-832.
[3]	张瑞, 王琨, 沈子聪, 毛路, 吴小涛. 内镜下椎间孔成形和椎板间成形对椎间盘、椎体峡部生物力学特性的影响[J]. 中国组织工程研究, 2024, 28(6): 833-839.
[4]	康志杰, 曹振华, 许阳阳, 张云凤, 金凤, 苏宝科, 王利东, 仝铃, 刘清华, 方源, 沙丽蓉, 亮亮, 李濛濛, 杜奕霏, 林琳, 王海燕, 李筱贺, 李志军. 强直性脊柱炎腰椎-骶骨椎间盘有限元模型建立及应力分析[J]. 中国组织工程研究, 2024, 28(6): 840-846.
[5]	张敏, 彭婧, 张强, 陈德旺. 有限元法分析老年骨质疏松患者L3/4椎板减压椎间融合的力学性能[J]. 中国组织工程研究, 2024, 28(6): 847-851.
[6]	黄培镇, 董航, 蔡群斌, 林梓凌, 黄枫. 顺行和逆行髓内钉治疗不同部位股骨干骨折的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 868-872.
[7]	谭能贤, 吴文正, 郑楚荣, 罗列良, 古鹏, 欧阳崇志, 郑晓辉. 不同空心加压螺钉固定方式对股骨颈垂直型骨折的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 873-878.
[8]	王明明, 张中, 孙建华, 赵刚, 宋华, 颜华东, 吕彬. 胫骨远端骨折伴软组织损伤3种不同微创固定方式的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 879-885.
[9]	侯泽欣, 许本柯, 戴媛, 何川, 张朝驹, 李孝林. 老年人跌倒腕关节背伸损伤机制的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 886-890.
[10]	宁天亮, 王坤, 王领彪, 韩鹏飞. 基于三点力原理对足内翻矫形器矫正效果的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 891-899.
[11]	章丽莉, 张兴来, 张杰, 郑洁皎. 视觉反馈对慢性踝关节不稳者落地过程中生物力学的影响[J]. 中国组织工程研究, 2024, 28(6): 900-904.
[12]	吴钟汉, 王景坤, 李涛, 许新忠, 余水生, 程里, 田大胜, 汤健, 荆珏华. 股骨近端防旋髓内钉治疗外侧壁完整型和外侧壁危险型股骨转子间骨折[J]. 中国组织工程研究, 2024, 28(6): 911-916.
[13]	赵汝顺, 郝阳泉, 许鹏, 郑鑫, 姜永宏, 张玉婷, 王孟飞, 鲁超. 不同坏死病灶位置对非创伤性股骨头坏死自然病程的影响[J]. 中国组织工程研究, 2024, 28(6): 917-921.
[14]	凯依塞尔•阿布都克力木, 麦麦提敏•阿卜力米提, 李磊, 杨晓凯, 张玉坤, 刘帅. 女性腰椎退行性病变患者腰椎CT值对骨质疏松症的诊断作用[J]. 中国组织工程研究, 2024, 28(6): 945-949.
[15]	王力平, 连天星, 胡永荣, 杨红胜, 曾智谋, 刘浩, 屈波. 胸部CT椎体HU值在2型糖尿病骨质疏松症机会性筛查中的价值[J]. 中国组织工程研究, 2024, 28(6): 950-954.