正常与骨质疏松髋关节模型的建立及有限元分析

doi:10.12307/2023.911

中国组织工程研究 ›› 2024, Vol. 28 ›› Issue (9): 1342-1346.doi: 10.12307/2023.911

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

正常与骨质疏松髋关节模型的建立及有限元分析

郭苏童1，冯德宏1，郭宇1，王凌1，丁育健1，刘仪1，钱正瑛2，李明洋3

南京医科大学附属无锡人民医院，1骨科，2医学工程处，江苏省无锡市 214000；3北京中诺恒康生物科技有限公司，北京市 102299

收稿日期:2022-12-09 接受日期:2023-02-04 出版日期:2024-03-28 发布日期:2023-07-25
通讯作者: 冯德宏，主任医师，教授，南京医科大学附属无锡人民医院骨科，江苏省无锡市 214000
作者简介:郭苏童，男，1996年生，江苏省连云港市人，汉族，南京医科大学在读硕士，主要从事关节外科、生物材料与骨组织工程研究。
基金资助:
无锡市社会发展科技示范(医疗卫生类) 项目(N20192006)，项目负责人：冯德宏；无锡市卫生健康委科研项目(重大项目)(Z202010)，项目负责人：冯德宏

Construction and finite element analysis of normal and osteoporotic hip models

Guo Sutong1, Feng Dehong1, Guo Yu1, Wang Ling1, Ding Yujian1, Liu Yi1, Qian Zhengying2, Li Mingyang3

1Department of Orthopedics, 2Medical Engineering Division, Wuxi People’s Hospital, Nanjing Medical University, Wuxi 214000, Jiangsu Province, China; 3Beijing Zhongnuo Hengkang Biotechnology Co., Ltd., Beijing 102299, China

Received:2022-12-09 Accepted:2023-02-04 Online:2024-03-28 Published:2023-07-25
Contact: Feng Dehong, Chief physician, Professor, Department of Orthopedics, Wuxi People’s Hospital, Nanjing Medical University, Wuxi 214000, Jiangsu Province, China
About author:Guo Sutong, Master candidate, Department of Orthopedics, Wuxi People’s Hospital, Nanjing Medical University, Wuxi 214000, Jiangsu Province, China
Supported by:
Wuxi Social Development Science and Technology Demonstration (Medical and Health) Project, No. N20192006 (to FDH); Scientific Research Project of Wuxi Health Commission (Major Project), No. Z202010 (to FDH)

摘要/Abstract

摘要：

文题释义：

髋关节：由股骨头及髋臼构成，属球窝关节，是典型的结构稳定的杵臼关节。
有限元分析：利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。利用简单且又相互作用的元素(即单元)，就可以用有限数量的未知量去逼近无限未知量的真实系统。

背景：骨密度是临床上判断骨骼强度的金标准，但骨密度对骨量变化的敏感较低，只有骨量明显降低时骨密度才会出现大幅变化，故骨密度对骨强度变化和骨折危险度的预测能力有限。

目的：建立正常与骨质疏松髋关节模型，分析单腿站立工况下正常及骨质疏松患者髋部应力及形变情况。
方法：选择1例健康成年女性志愿者为研究对象，年龄36岁，获得该志愿者的髋部CT数据并以DICOM格式保存。对髋关节模型进行三维重建，通过灰度赋值法赋予材料属性，按照经验公式，获得正常与骨质疏松髋关节模型。设定相同的边界条件和载荷，模拟单腿站立位状态下正常与骨质疏松髋关节应力及形变情况。

结果与结论：①在正常及骨质疏松髋关节有限元模型中，股骨颈内侧区域应力分布较为集中；②在髋骨中，应力分布主要集中于髋臼上部；③正常髋关节模型比骨质疏松性髋关节模型在股骨颈内侧、髋臼上部的应力峰值大，可能是由于骨质疏松性骨骼骨强度降低导致；④正常及骨质疏松髋关节模型的Von Mises峰值都集中于股骨颈内侧，髋骨Von Mises峰值较小，说明骨质疏松对髋骨受力整体影响相对较小；⑤单腿站立位下形变方面，正常髋关节模型最大形变位于髋臼与股骨头处，骨质疏松髋关节模型最大形变位于股骨大转子上部；⑥提示有限元分析法模拟骨质疏松症骨组织的相关参数，可能会提高临床上对骨质疏松患者骨强度变化的监测和骨折风险的预测能力，从生物力学角度解释了股骨转子间、股骨颈是骨质疏松性髋部骨折的好发部位。

https://orcid.org/0000-0001-5582-7249 (郭苏童)

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

关键词: 髋关节, 骨质疏松, 灰度值赋值法, 应力, 形变, 有限元

Abstract: BACKGROUND: Bone mineral density is the clinical gold standard for determining bone strength, but bone mineral density is less sensitive to changes in bone mass, with large changes in bone mineral density only occurring when bone mass is significantly reduced, so bone mineral density has limited ability to predict changes in bone strength and fracture risk.
OBJECTIVE: A model of the normal and osteoporotic hip joint was developed to analyze the stresses and deformation in the hip of normal and osteoporotic patients under single-leg standing conditions.
METHODS: A healthy adult female volunteer at the age of 36 years was selected as the study subject. The CT data of the hip joint of this volunteer were obtained and saved in DICOM format. The hip joint model was reconstructed in three dimensions, and the material properties were assigned by the gray value assignment method to obtain the normal and osteoporotic hip joint models according to the empirical formula. The same boundary conditions and loads were set to simulate the stresses and deformation in the normal and osteoporotic hip joints in the single-leg standing position.
RESULTS AND CONCLUSION: (1) In the finite element model of the normal and osteoporotic hip, the stress distribution was more concentrated in the medial region of the femoral neck. (2) In the hip bone, the stress distribution was mainly concentrated in the upper part of the acetabulum. (3) The stress peaks in the medial femoral neck and upper acetabulum were larger in the normal hip model than in the osteoporotic hip model, probably due to the reduced bone strength of the osteoporotic bone. (4) The peak Von Mises of both normal and osteoporotic hip models were concentrated on the medial femoral neck, and the peak Von Mises of the hip bone was smaller, indicating that the overall effect of osteoporosis on hip bone stresses was relatively small. (5) In terms of deformation in the single-leg standing position, the maximum deformation in the normal hip model was located at the acetabulum and femoral head, and the maximum deformation was located at the upper part of the greater trochanter of the femur. (6) It is suggested that the finite element analysis method to model the values of parameters related to bone tissue in osteoporosis may improve clinical prediction of bone strength changes and fracture risk. It is explained from the biomechanical view that the intertrochanteric femur and femoral neck are good sites for osteoporotic hip fractures.

Key words: hip joint, osteoporosis, gray value assignment method, stress, deformation, finite element

中图分类号:

郭苏童, 冯德宏, 郭宇, 王凌, 丁育健, 刘仪, 钱正瑛, 李明洋. 正常与骨质疏松髋关节模型的建立及有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1342-1346.

Guo Sutong, Feng Dehong, Guo Yu, Wang Ling, Ding Yujian, Liu Yi, Qian Zhengying, Li Mingyang. Construction and finite element analysis of normal and osteoporotic hip models[J]. Chinese Journal of Tissue Engineering Research, 2024, 28(9): 1342-1346.

图/表（结果） 4

参考文献

[1] YAVROPOULOU MP, MAKRAS P, ANASTASILAKIS AD. Bazedoxifene for the treatment of osteoporosis. Expert Opin Pharmacother. 2019;20(10):1201-1210.
[2] NELLANS KW, KOWALSKI E, CHUNG KC. The epidemiology of distal radius fractures. Hand Clin. 2012;28(2):113-125.
[3] 李祖国, 扈佐鸿, 刘浩, 等.单侧经皮椎体成形术治疗骨质疏松性上胸椎骨折的临床疗效观察[J].中华临床医师杂志(电子版),2015,9(15): 2820-2823.
[4] 潘盛, 杨佳生, 赖以毅, 等.CT与MRI在骨质疏松性与非骨质疏松性腰椎骨折中的临床应用[J].医学影像学杂志,2022,32(9):1638-1641.
[5] GUZON-ILLESCAS O, PEREZ FERNANDEZ E, CRESPÍ VILLARIAS N, et al. Mortality after osteoporotic hip fracture: incidence, trends, and associated factors. J Orthop Surg Res. 2019;14(1):203.
[6] KARRES J, KIEVIET N, EERENBERG JP, et al. Predicting Early Mortality After Hip Fracture Surgery: The Hip Fracture Estimator of Mortality Amsterdam. J Orthop Trauma. 2018;32(1):27-33.
[7] QASEEM A, FORCIEA MA, MCLEAN RM, et al. Treatment of Low Bone Density or Osteoporosis to Prevent Fractures in Men and Women: A Clinical Practice Guideline Update From the American College of Physicians. Ann Intern Med. 2017;166(11):818-839.
[8] SCHREIBER JJ, KAMAL RN, YAO J. Simple Assessment of Global Bone Density and Osteoporosis Screening Using Standard Radiographs of the Hand. J Hand Surg Am. 2017;42(4):244-249.
[9] WU S, TODO M, UMEBAYASHI D, et al.Risk assessment of vertebral compressive fracture using bone mass index and strength predicted by computed tomography image based finite element analysis.Clin Biomech (Bristol, Avon). 2021;85:105365.
[10] 郭文文, 刘静, 曹慧, 等.正常与骨质疏松肱骨的三维重建及有限元分析[J].中国医疗设备,2018,33(4):34-37.
[11] 杨锐敏, 吴文正, 郑永泽, 等.不同松质骨体积分数影响股骨近端表观力学响应的有限元分析[J].中国组织工程研究,2021,25(36):5765-5770.
[12] 秦大平, 张晓刚, 权祯, 等.不同方法治疗骨质疏松性胸腰椎压缩骨折椎体力学稳定性变化差异的有限元分析[J].中华中医药杂志,2021,36(8): 4886-4895.
[13] Pustoc’h A, Cheze L.Normal and osteoarthritic hip joint mechanical behaviour: a comparison study. Med Biol Eng Comput. 2009;47(4):375-383.
[14] 陈凯宁, 农明善, 叶青, 等.钙化软骨层硬度变化对骨软骨结构应力影响的有限元分析[J].中国组织工程研究,2019,23(12):1881-1886.
[15] WRONSKI S, KAMINSKI J, WIT A, et al. Anisotropic bone response based on FEM simulation and real micro computed tomography of bovine bone. Comp Methods Biomech Biomed Eng. 2019;22:S465 - S467.
[16] 郭苏童, 郭宇, 王凌, 等.发育性髋关节发育不良患者不同髋关节外展角度下髋部应力分布的有限元分析[J].中国组织工程研究,2023,27(27): 4265-4270.
[17] YANG P, LIN TY, XU JL, et al.Finite element modeling of proximal femur with quantifiable weight-bearing area in standing position. J Orthop Surg Res. 2020;15(1):384.
[18] 王汝良, 张凯旋, 胡霖霖, 等.基于CT的正常与骨质疏松股骨模态分析[J].软件,2019,40(4):77-80.
[19] GRECU D, PUCALEV I, NEGRU M, et al. Numerical simulations of the 3D virtual model of the human hip joint, using finite element method. Rom J Morphol Embryol. 2010;51(1):151-155.
[20] 张海峰, 尹爱华, 董毅, 等.有限元法分析不同负荷下髋臼区的应力分布[J].中国组织工程研究,2016,20(39):5867-5872.
[21] 丁海, 朱振安.髋臼的解剖形态及生物力学研究进展[J].医用生物力学, 2008,23(5):411-414.
[22] 张馨元, 丁晓红, 段朋云, 等.基于骨骼材料分布特征的髋关节有限元模型建立及其应力分析[J].北京生物医学工程,2020,39(5):462-469.
[23] 陈国栋, 罗羽婕,王锐英.有限元分析在股骨生物力学研究中的应用[J].实用医学杂志,2011,27(2):334-336.
[24] 陈诺, 邱兴, 汪帝, 等.基于CT灰度值赋予材料属性的马拉松运动员髋关节有限元模型建立及其应力分析[J].中华医学杂志,2022,102(9):679-682.
[25] Honig S, Chang G. Osteoporosis: an update. Bull NYU Hosp Jt Dis. 2012; 70(3):140-144.
[26] 石少辉, 潘伟, 吴国平, 等. 60岁及以上髋部和椎体骨折患者血清25-羟基维生素D水平分析[J].中华全科医师杂志,2020,19(3):233-237.
[27] 唐佩福.老年髋部骨折的诊治现状与进展[J].中华创伤骨科杂志,2020,22(3):197-199.
[28] 原发性骨质疏松症诊疗指南(2017).中华骨质疏松和骨矿盐疾病杂志, 2017,10(5):413-444.
[29] WAINWRIGHT SA, MARSHALL LM, ENSRUD KE, et al. Hip fracture in women without osteoporosis. J Clin Endocrinol Metab. 2005;90(5):2787-2793.
[30] BOUXSEIN ML, KARASIK D. Bone geometry and skeletal fragility. Curr Osteoporos Rep. 2006;4(2):49-56.
[31] HOLLENSTEINER M, SANDRIESSER S, BLIVEN E, et al. Biomechanics of Osteoporotic Fracture Fixation. Curr Osteoporos Rep. 2019;17(6):363-374.
[32] CHRISTEN D, WEBSTER DJ, MÜLLER R. Multiscale modelling and nonlinear finite element analysis as clinical tools for the assessment of fracture risk.Philos Trans A Math Phys Eng Sci. 2010;368(1920):2653-2668.
[33] KOPPERDAHL DL, ASPELUND T, HOFFMANN PF, et al.Assessment of incident spine and hip fractures in women and men using finite element analysis of CT scans. J Bone Miner Res. 2014;29(3):570-580.
[34] 秦大平, 张晓刚, 宋敏, 等.有限元分析在骨质疏松性椎体压缩骨折脊柱力学动态变化中的应用[J].中华中医药杂志,2019,34(1):206-211.
[35] 杜根发. 基于断裂力学探讨骨密度影响老年股骨颈骨折的有限元分析[D].广州:广州中医药大学,2017.
[36] 郑利钦, 林梓凌, 何祥鑫, 等.动态载荷下股骨转子间区域皮质骨厚度对骨折类型影响的有限元分析[J].医学研究生学报,2018,31(10):1043-1046.
[37] 郑利钦, 林梓凌, 陈心敏, 等.载荷速率对股骨颈骨折裂纹扩展影响的有限元分析[J].中国组织工程研究,2019,23(20):3148-3152.
[38] 丁海, 朱振安, 薛晶, 等.骨质疏松症对松质骨骨小梁应力与微损伤关系的影响[J].医用生物力学,2015,30(1):68-73.
[39] 高耀东, 段宇星, 郭鹏年, 等.基于髋关节CT图像三维实体髋骨重建及髋骨承受力数据有限元分析[J]. 中国组织工程研究,2019,23(28):4564-4569.
[40] HAO Z, WAN C, GAO X, et al.The effect of boundary condition on the biomechanics of a human pelvic joint under an axial compressive load: a three-dimensional finite element model. J Biomech Eng. 2011;133(10):101006.
[41] ARKUSZ K, KLEKIEL T, NIEZGODA N, et al. The influence of osteoporotic bone structures of the pelvic-hip complex on stress distribution under impact load. Acta Bioeng Biomech. 2018;20(1):29-38.
[42] 刘欣伟, 闫寒, 刘中洋, 等.包含肌肉、韧带组织的骨盆、髋臼3D有限元模型的构建[J].临床军医杂志,2014,42(4):331-335.
[43] 骆健, 王立华, 王涛, 等.不同材料赋值方法下踝关节三维有限元模型的应力及位移变化[J].中国组织工程研究,2019,23(18):2822-2826.

引言

骨质疏松症是一种以单位体积内骨量减少、骨脆性增加、骨密度降低、骨皮质变薄及骨显微结构改变为主要特征的全身性骨代谢性相关疾病[1]。由于骨质疏松导致骨结构易损，患者极易发生脆性骨折，如桡骨骨折[2]、胸椎骨折[3]、腰椎骨折等[4]。在所有的脆性骨折及并发症中，髋部骨折是一种常见且严重的类型[5]。近年来，国内人口老龄化趋势日益明显，骨质疏松症发病率逐年升高，骨质疏松髋部骨折的发生率也随之升高[6]。骨密度是临床上判断骨骼强度的金标准，其水平降低预示骨质疏松骨折发生的风险增加，但骨密度对骨量变化的敏感较低，只有骨量明显降低时骨密度值才会出现大幅变化[7-8]。随着计算机技术和有限元理论的发展，通过计算机建立更接近真实状态的有限元模型评估骨折风险和治疗效果的方法为骨质疏松的诊治提供了重要的参考。WU等[9]基于CT图像预测骨密度和骨强度评估椎体压缩性骨折的风险，并进行了有限元分析，结果表明骨密度值不能精准预测骨质疏松导致的椎体骨折风险，通过引入基于骨密度和基于骨强度的指标，可以更精确地评估椎体骨折风险。郭文文等[10]建立了正常肱骨和骨质疏松肱骨三维有限元模型，并分析了肱骨骨折的机制。杨锐敏等[11]按梯度建立了不同松质骨体积分数的股骨近端三维模型，探究松质骨体积分数对股骨近端表观力学响应的影响，为髋部骨折发生机制研究提供参考。秦大平等[12]通过有限元方法分析不同微创方法治疗骨质疏松性椎体压缩骨折伤椎高度恢复不全状态下，不同工况下脊柱及相邻椎体附属结构的应力变化，为骨质疏松性椎体压缩骨折筛选治疗策略提供依据。
此次研究通过计算机软件系统对正常髋关节CT数据进行处理，利用Mimics 软件、3-matic软件基于灰度值赋值法建立正常与骨质疏松症髋关节的三维有限元模型，并利用Ansys Workbench 软件进行分析，比较单腿站立工况下正常与骨质疏松症髋关节有限元模型的应力及形变情况。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

材料方法

1.1 设计建立正常及骨质疏松髋关节模型，并进行有限元分析。
1.2 时间及地点实验于2022年8-10月在在南京医科大附属无锡人民医院完成。
1.3 对象选择1例健康成年女性志愿者为研究对象，36岁，身高160 cm，体质量50 kg，既往无髋关节骨折、肿瘤、畸形、骨代谢性疾病及肝肾功能严重异常等病史，既往未服用如双膦酸盐类和激素类及其他影响骨代谢的药物。获得该志愿者的髋关节CT数据后以DICOM格式保存。
志愿者及家属对实验知情同意，并与志愿者签署“知情同意书”。此次研究的实施符合《赫尔辛基宣言》和南京医科大学附属无锡人民医院对研究的相关伦理规定，伦理批件号为：(2017)伦审第(01)号。
1.4 材料
1.4.1 建模设备 CT扫描机(南京医科大学附属无锡人民医院提供)；华硕品牌笔记本电脑(硬件配置：CPU为intel i5-4200U 处理器，硬盘750 G，内存8 G)。
1.4.2 操作系统 Windows 7 旗舰版64位软件操作系统(Microsoft公司，美国)。
1.4.3 图像处理及分析软件医学影像数据处理软件Mimics 21.0(Materialise公司，比利时)；网格划分软件3-matic 13.0(Materialise公司，比利时)；约束条件界定的有限元分析专用软件Ansys workbench 18.2(Ansys公司，美国)。
1.5 方法
1.5.1 有限元模型的构建志愿者取仰卧位，双膝关节并拢伸直，髋关节放松，处中立位，位于扫描视野中心，纵轴方向保持不动。采用像素512×512，层厚1.0 mm的髋关节CT数据，获取包括矢状面、横断面、冠状面的图像信息，同时以DICOM格式存储。
将DICOM格式的图像数据导入Mimics软件后，确定图像方向与人体坐标轴的方向一致。选定合适的灰度值，根据不同组织成像的差异性，显示骨骼与软组织间的区别，提取出髋部骨骼的轮廓。对图像进行阈值分割、蒙版编辑、区域增长等操作后，三维重建，得到最初的髋关节骨骼三维模型。将Mimics软件中生成的髋关节三维有限元模型导入3-matic软件中，检查三维有限元模型后进一步优化等处理。进行优化时需注意力度，保持模型的基本特征。髋关节软骨的主要作用是润滑和缓震，为模拟真实情况下髋部骨与骨之间力的传导，添加软骨结构十分有必要。结合既往研究及此模型关节实际间隙，通过3-matic软件建立厚度1.0 mm的股骨头软骨和厚度为1.0 mm的髋臼软骨，此方法建立的软骨与骨骼贴合较好[13-14]。

有限元模型建立过程复杂，模型处理须精确，否则将产生实验误差，甚至无法得出正确结果。在骨骼建模中，多数研究者采用逆向建模技术重建出三维模型，但不同软件之间的切换、导入，易导致模型的信息丢失，且建模费时费力，建模结果的准确性会影响后续的计算结果。此次研究基于Mimics软件对髋关节CT数据进行重建，利用3-matic软件生成高质量的有限元模型，在保证模型精度的同时减少了时间成本。建模结果如图1所示。

1.5.2 网格划分有限元模型网格的形状、大小、疏密及分布与结果的准确性密切相关。划分网格不仅要考虑髋关节的结构特点，还需考虑运算的收敛性及分析的精确性。建立好髋关节模型后，进入3-matic软件自动生成面网格，过滤细小的三角片、减少网格数量、网格自动优化、填补空洞、检查相交网格并删除，最后创建体网格。此次研究共生成体网格单元数为154 561个，节点数为277 913个。经网格检验，此次研究无不合格网格。网格划分结果见图2。

1.5.3 骨骼材料属性赋予当前建模方法多采用均质材料属性，然而在实际生理状态下，人体骨骼是由多种各向异性的材料构成[15]，不同个体的骨组织在生物力学表现上不尽相同[16]，不能反映出个体的骨质特性，故建立个体化材料属性的模型对髋关节有限元分析十分重要。将3-matic中建立的髋关节模型体网格导入Mimics 21.0软件中，通过材料赋值模块可以得到CT数据上每一点的灰度值，通过公式计算对每个网格单元进行材料属性赋予。依据下面文献中的经验公式ρ=-13.4+1 017Gv，E=-388.8+5 925ρ(ρ为密度，Gv为灰度值，E为弹性模量)，泊松比设为0.3[17]，对正常髋关节有限元模型进行材料参数赋值。以正常髋关节骨骼的弹性模量减少66%、泊松比保持不变为条件，赋予髋关节骨骼材料属性[10，18]，生成骨质疏松髋关节骨骼有限元模型。灰度值材料属性分配结果见图3。赋予髋臼软骨及股骨头软骨弹性模量为10.5 MPa，泊松比为0.45[19]。

1.5.4 设置边界条件及负载分别将2种模型以CDB格式导入ANSYS Workbench软件，在股骨上段的下端横断面施加垂直向上的载荷[20]。因单腿站立状态下髋关节承受的负荷约为体质量的81%，故施加大小为患者体质量81%的载荷，为405 N[21]。对髋骨近端及耻骨联合处进行约束，在X、Y、Z三个方向上均无位移和转位。
1.6 主要观察指标通过模拟加载进行有限元仿真分析，观察单腿站立工况下正常及骨质疏松髋关节有限元模型的应力及形变情况。

讨论

骨质疏松性髋部脆性骨折是在无创伤或轻、中度创伤下出现的局部骨组织病变，是髋部骨强度下降的体现，是骨质疏松症的最终结果之一[25]。髋部骨折后，因活动减少会出现失用性骨丢失，这将进一步加剧骨质疏松症的严重程度。严重骨质疏松髋部骨折患者易发生医源性骨折、骨缺损及内固定失效、假体松动等并发症，而且骨质疏松骨折后二次骨折的风险明显升高[26]。且股骨头血供特殊，骨折部位血供减少，还易引发股骨头缺血性坏死。故老年髋部骨折患者功能恢复不佳，生活质量降低，甚至可能引发死亡[27]。
骨质疏松症是老年髋部骨折不容忽视的重要危险因素，但骨质疏松症缺少典型症状，往往发生骨折才开始临床干预。X射线片、CT、MRI和核医学等影像检查结果以及双能X射线骨密度测量、定量CT骨密度测量结果是骨质疏松症诊断、风险预测和疗效评价的主要依据[28]。作为诊断骨质疏松的金标准，骨密度测量间接反映骨量高低和骨骼强度从而预测骨质疏松骨折的发生风险[11]。然而，单纯依靠骨密度进行诊断具有一定的局限性，准确度较低，较易漏诊。约54%的骨质疏松髋部骨折患者骨密度值未达到骨质疏松症的诊断标准[29]。除了骨密度值外，骨骼的几何特点、密度空间分布和骨材料特性也对骨强度有影响[30]。骨质疏松症评价的最佳途径是测量骨强度的方法[31]，但是临床上鲜有对骨强度的无创测量手段。有研究表明，有限元分析比骨密度可更有效地评估骨折风险[32-33]。对于骨质疏松的有限元研究主要集中在肱骨[10]、股骨[11]、腰椎等部位[34]，鲜有使用有限元方法对骨质疏松性髋关节的研究。随着科学技术的发展，有限元法以其准确、快速地评价骨强度而得到广泛应用。有限元分析法模拟骨质疏松症骨组织的相关参数，可能会提高临床上对骨质疏松患者骨强度变化的监测和骨折风险的预测能力[29]。
松质骨承担主要弹性形变，骨松质中骨小梁的断裂、减少与骨皮质的变薄、微损伤在脆性骨折的发生起着重要作用[35]。有研究表明转子间皮质骨越薄，应力范围越趋向分散于转子间部位，最终可能形成转子间骨折[36]。大部分髋关节有限元研究基于静态条件下，认为应力集中区域为骨折起始部位，这种方法预测的骨骼应力情况、骨强度与实体实验的关联性较强。但仅依靠应力云图常难以预测骨折发生的准确部位[37]。施加载荷后，骨质发生应变，若载荷超过应变范围也将发生骨折[38]。
在正常髋关节有限元模型中，股骨颈的内侧区域应力最大，Von Mises峰值为4.263 MPa；在髋骨中，髋臼上部应力分布较为集中，Von Mises峰值为4.257 MPa。在骨质疏松髋关节有限元模型中，股骨颈内侧区域应力分布较为集中，Von Mises峰值为4.249 MPa；在髋骨中，应力分布主要于集中髋臼上部，Von Mises峰值为4.209 MPa。正常髋关节模型比骨质疏松性髋关节模型在股骨颈内侧、髋臼上部的应力峰值大，可能是由于骨质疏松性骨骼骨强度降低导致。正常及骨质疏松的Von Mises峰值都集中于股骨颈内侧，髋骨Von Mises峰值较小，说明骨质疏松对髋骨受力整体影响相对较小。单腿站立位下形变方面，正常髋关节模型最大形变位于髋臼与股骨头处，为0.643 mm；骨质疏松髋关节模型最大形变位于股骨大转子上部，为0.600 mm。可能因为正常站立位应力下松质骨具有弹性稳定作用，随着骨质疏松患者股骨大转子处骨松质减少，应变量上升[11]。这些都从生物力学的角度上说明股骨转子间、股骨颈是骨质疏松性髋部骨折的好发部位。
髋关节的解剖复杂，模拟其应力及应变分布较为困难，而既往研究仅建立了股骨或髋骨三维模型并分析其应力，结果或存在较大误差[11，39]。此次研究构建了含髋骨、股骨近端、髋臼软骨、股骨头软骨的有限元模型，或可在一定程度上解决这个问题。有研究构建了无骶髂韧带髋关节、无股骨近端骨盆髋关节及完整髋关节的有限元模型[40]，也有研究的建模方式不仅涉及了股骨近端、髋骨、骶骨，而且考虑到了韧带等软组织结构对应力分布的影响[41]，这些研究均表明有限元分析受到边界条件设置的影响较大。但也有研究表明构建髋骨周围的肌肉及韧带作为有限元分析的边界条件，不但费时费力，而且可能与实际情况有出入[42]。
此次研究利用Mimics软件基于髋关节CT影像数据进行三维重建，利用3-matic软件对模型进行优化、光顺、网格划分，生成高质量的有限元模型。后基于灰度值属性赋予正常及骨质疏松髋关节模型的材料分布特性，既减少了因不同软件的文件导入造成信息缺失而降低准确性，又保证了模型的网格划分质量。再将所建模型通过有限元分析软件ANSYS Workbench设定边界条件，模拟加载负荷，计算单腿站立位工况下正常与骨质疏松患者髋关节的应力及形变分布情况，结果真实可信。此外，此次研究模型可通过施加不同负荷等进一步研究其他运动状态下髋关节的应力分布[20，43]，还可根据研究内容模拟骨折、骨质增生、股骨头坏死等状态，为患者疾病的预防、诊断、治疗、康复等提供科学的参考[16]。此次研究所建髋关节个体化有限元模型与实际高度符合，能得到更精确、更接近实际的结果，所提供方法或可成为一种预测骨质疏松性骨折的方法，为髋部骨折发生机制的研究及疾病的预防、诊治等提供参考和借鉴。
然而，此次研究存在以下不足和局限性：第一，此模型只涉及骨骼和软骨，后续研究可建立包含肌肉、韧带等软组织的有限元模型并进行力学分析，以期获得更接近人体实际的分析结果；第二，此次研究样本量小，所建髋关节模型是根据单一受试者的几何信息建立的，后续需进行更多深入且完善的研究；第三，不同研究所建立的髋关节模型的建模方法、载荷及边界条件常不同，难以定量比较；第四，髋关节的材料参数与实际或存在误差，后续应将有限元与实体试验结合并进行对比分析，以求得到更为精确的结果。中国组织工程研究杂志出版内容重点：人工关节；骨植入物；脊柱；骨折；内固定；数字化骨科；组织工程

正常与骨质疏松髋关节模型的建立及有限元分析

Construction and finite element analysis of normal and osteoporotic hip models

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 4

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	李智斐, 杨尹, 陈华龙, 梁钦秋, 钟远鸣, 张翼升. 颈椎前路椎体次全切除减压融合钛笼倾斜角度与术后钛笼下沉的有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1313-1319.
[2]	欧阳北平, 马向阳, 罗春山, 邹小宝, 陆廷盛, 陈啟鸰. 新型钉尾横连在后路寰枢椎内固定中的三维有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1320-1324.
[3]	梁成, 张麟麒, 王冠, 李文, 段可, 李忠, 鲁晓波, 卓乃强. 不同内植物修复单侧不稳定骨盆后环损伤的有限元及生物力学分析[J]. 中国组织工程研究, 2024, 28(9): 1336-1341.
[4]	杨骏良, 路坦, 徐彪, 姜亚琼, 王富成. 前交叉韧带部分断裂对膝关节应力影响的三维有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1347-1353.
[5]	李超杰, 居来提•买提肉孜, 艾尔西丁•阿不来提, 郑辉, 涂浒第. 不同屈曲角度时内侧髌股韧带重建的有限元分析[J]. 中国组织工程研究, 2024, 28(9): 1359-1364.
[6]	闵美鹏, 吴进, URBA RAFI, 张文杰, 高嘉, 王云华, 何斌, 范磊. 全髋关节置换前人工智能(AI)规划的作用与意义[J]. 中国组织工程研究, 2024, 28(9): 1372-1377.
[7]	杨策凯, 蔡卓延, 陈明, 刘昊, 翁汭, 崔健超, 张顺聪, 姚珍松. 绝经后女性椎旁肌退化与经皮穿刺椎体成形后再骨折的相关性[J]. 中国组织工程研究, 2024, 28(9): 1414-1419.
[8]	霞黑达·依拉尔江, 尼加提·吐尔逊, 热依拉·库尔班, 白布加甫·叶力思, 陈欣. 三维有限元分析不同特征骨组织中应力的分布规律[J]. 中国组织工程研究, 2024, 28(8): 1277-1282.
[9]	章晓云, 刘桦, 柴源, 陈锋, 曾浩, 高振罡, 黄有荣. 益肾固疏方干预老年性骨质疏松症患者骨代谢标志物的变化及临床疗效[J]. 中国组织工程研究, 2024, 28(8): 1155-1160.
[10]	代越星, 郑利钦, 吴敏辉, 李志鸿, 李少彬, 郑德声, 林梓凌. 血管数量对小血管网计算流体力学的影响[J]. 中国组织工程研究, 2024, 28(8): 1206-1210.
[11]	童奕博 , 李明辉. 骨质疏松性椎体骨折患者椎体成形后邻近椎体再发骨折的影响因素[J]. 中国组织工程研究, 2024, 28(8): 1241-1246.
[12]	魏源标, 林展, 陈燕梅, 杨腾辉, 赵晓, 陈扬声, 周燕辉, 杨敏超, 黄飞麒. 颈椎矢状位端提手法对椎间盘及小关节影响的有限元分析[J]. 中国组织工程研究, 2024, 28(6): 827-832.
[13]	张瑞, 王琨, 沈子聪, 毛路, 吴小涛. 内镜下椎间孔成形和椎板间成形对椎间盘、椎体峡部生物力学特性的影响[J]. 中国组织工程研究, 2024, 28(6): 833-839.
[14]	康志杰, 曹振华, 许阳阳, 张云凤, 金凤, 苏宝科, 王利东, 仝铃, 刘清华, 方源, 沙丽蓉, 亮亮, 李濛濛, 杜奕霏, 林琳, 王海燕, 李筱贺, 李志军. 强直性脊柱炎腰椎-骶骨椎间盘有限元模型建立及应力分析[J]. 中国组织工程研究, 2024, 28(6): 840-846.
[15]	张敏, 彭婧, 张强, 陈德旺. 有限元法分析老年骨质疏松患者L3/4椎板减压椎间融合的力学性能[J]. 中国组织工程研究, 2024, 28(6): 847-851.