基于CT灰度值赋值的股骨头坏死有限元模型对比

doi:10.3969/j.issn.2095-4344.2414

中国组织工程研究 ›› 2020, Vol. 24 ›› Issue (3): 395-400.doi: 10.3969/j.issn.2095-4344.2414

• 骨与关节图像与影像 bone and joint imaging • 上一篇下一篇

基于CT灰度值赋值的股骨头坏死有限元模型对比

薛志鹏，李泰贤，李䶮，何海军，黄泽青，孙继高，陈卫衡

中国中医科学院望京医院，北京市 100102

收稿日期:2019-06-11 修回日期:2019-06-12 接受日期:2019-07-20 出版日期:2020-01-28 发布日期:2019-12-26
通讯作者: 陈卫衡，博士，教授，博士生导师，中国中医科学院望京医院骨关节三科，北京市 100102
作者简介:薛志鹏，男，1991年生，内蒙古自治区赤峰市人，中国中医科学院在读博士研究生，医师，主要从事骨关节疾病的临床与基础研究。
基金资助:
国家自然科学基金(81373801)；北京市自然科学基金(7182186)

Comparison of finite element models of osteonecrosis of the femoral head based on CT gray-assigned method

Xue Zhipeng, Li Taixian, Li Yan, He Haijun, Huang Zeqing, Sun Jigao, Chen Weiheng

Wangjing Hospital of China Academy of Chinese Medical Sciences, Beijing 100102, China

Received:2019-06-11 Revised:2019-06-12 Accepted:2019-07-20 Online:2020-01-28 Published:2019-12-26
Contact: Chen Weiheng, MD, Professor, Doctoral supervisor, Wangjing Hospital of China Academy of Chinese Medical Sciences, Beijing 100102, China
About author:Xue Zhipeng, Doctoral candidate, Physician, Wangjing Hospital of China Academy of Chinese Medical Sciences, Beijing 100102, China
Supported by:
the National Natural Science Foundation of China, No. 81373801; the Natural Science Foundation of Beijing, No. 7182186

摘要/Abstract

摘要：

文题释义：
CT灰度值赋值：此赋值法利用CT灰度的不同，将组织划分为不同的等级材料，利用数学公式将CT灰度值与弹性模量进行转换运算，将不同骨组织材料进行相应赋值，能够较真实的体现骨组织内部的复杂情况，使模型更加准确。
股骨头坏死应力集中：正常股骨头内应力延压力性骨小梁均匀分布在股骨头内，坚硬的股骨距承受大部分应力。股骨头坏死的应力因坏死区骨小梁的断裂、硬化带及肉芽带的形成无法均匀传导下行，股骨头负重区、坏死区应力分布不均匀出现多个应力集中区域，与正常股骨头不同的是，应力分散于股骨颈内外侧。

背景：目前研究股骨头有限元模型多采用单例或较少样本建模用于特定的生物力学研究，对于模型稳定性研究不多。

目的：以正常股骨头和股骨头坏死建模进行多样本的模型对照，通过应力分布规律和力学参数的对比，分析模型的准确性和稳定性，为股骨头坏死塌陷防治提供力学基础。

方法：选择20髋经1年非手术治疗稳定未塌陷的股骨头为实验组，20髋单侧股骨头坏死患者的健侧设为正常组。采集股骨头CT数据分别建立有限元模型，观察正常股骨头与股骨头坏死应力分布、股骨头负重区最大等效应力与最大总形变量，对比并统计分析。研究方案经中国中医科学院望京医院医学伦理委员会批准，患者签署知情同意书。

结果与结论：①建立了正常股骨近端、无坏死股骨近端和坏死骨有限元模型，单元数和节点数分别为 502 568±114 196、692 608±154 678，449 954±125 824、623 311±171 401，19 133±13 167、27 577± 19 131；②模拟单足站立位进行载荷设定的云图显示：当施以体质量2.5倍压强作用于股骨头负重区表面时，正常股骨头负重区表面应力均匀，应力延压力性骨小梁均匀分布在股骨头内，股骨距承受大部分应力；股骨头坏死负重区表面及坏死区出现应力集中区域，头内应力分散于股骨头颈交界，股骨颈内外侧应力相当，且股骨头坏死产生的形变多于正常股骨头；③股骨头坏死与正常股骨头的负重区最大总形变量分别为(4.14±1.31)，(1.36±0.22) mm，最大等效应力分别为(1.94±0.77)，(0.75±0.19)MPa，差异有统计学意义(P < 0.05)，且两组数据偏于集中，模型稳定性较好。通过多样本的正常股骨头与股骨头坏死比较，证明了基于CT灰度值直接赋值反映了股骨头坏死的实际力学特性，具有较好的准确性和稳定性。

ORCID: 0000-0001-5227-532X(薛志鹏)；0000-0003-0831-1201(陈卫衡)

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

关键词: CT灰度值赋值, 股骨头坏死, 有限元模型, 单元数, 节点数, 应力分布, 最大等效应力, 最大总形变量

Abstract:

BACKGROUND: Previous studies have reported that femoral head finite element models are mostly modeled with single or few samples for specific biomechanical research, but there is little research on model stability.

OBJECTIVE: To compare the models of normal femoral head and osteonecrosis of the femoral head with multiple samples, and to analyze the accuracy and stability of the models through the comparison of stress distribution and mechanical parameters, so as to provide mechanical basis for prevention and treatment of collapse of osteonecrosis of the femoral head.

METHODS: Totally 20 sides uncollapsed of osteonecrosis of the femoral head one year of non-surgical treatment were selected as the experimental group, and the healthy side of 20 patients with unilateral osteonecrosis of the femoral head were set as the normal group. The CT data of the femoral head were collected to establish the finite element model. The stress distribution of normal femoral head and osteonecrosis of the femoral head, the maximum equivalent stress and the maximum total deformation at the weight-bearing area of the femoral head were observed and compared. This study was approved by the Medical Ethics Committee of Wangjing Hospital of China Academy of Chinese Medical Sciences. Patients signed the informed consent.

RESULTS AND CONCLUSION: (1) The finite element models of normal proximal femur, non-necrotic proximal femur and necrotic bone were established. The number of elements and nodes were 502 568±114 196, 692 608±154 678; 449 954±125 824, 623 311±171 401; 19 133±13 167, 27 577±19 131, respectively. (2) When the load was set by simulating one-foot standing position, the cloud image showed that when 2.5 times body weight applied to the weight-bearing area of the femoral head; the surface stress of the weight-bearing area of the normal femoral head was uniform. The stress was uniformly distributed in the femoral head along the stress trabeculae, and the calcar femorale bears the most. The stress concentration area appeared on the surface of the weight-bearing area and the necrotic area of osteonecrosis of the femoral head. The stress was scattered and distributed on the inner and outer sides of the femoral neck and the femoral head of osteonecrosis of the femoral head produced more deformation than the normal femoral head. (3) The maximum total deformation of the weight-bearing area of the osteonecrosis of the femoral head and the normal femoral head was (4.14±1.31) mm and (1.36±0.22) mm and the maximum equivalent stress was (1.94±0.77) MPa and (0.75±0.19) MPa, respectively, and with statistically significance (P < 0.05). Moreover, two groups of data tend to be concentrated and the models are stable. Through the comparison of multi-sample normal femoral head and osteonecrosis of the femoral head, the CT gray-assigned method reflects the actual mechanical properties of osteonecrosis of the femoral head, and has good accuracy and stability.

Key words: CT gray-assigned method, osteonecrosis of femoral head, finite element model, number of elements, number of nodes, stress distribution, maximum equivalent stress, maximum total deformation

中图分类号:

薛志鹏, 李泰贤, 李䶮, 何海军, 黄泽青, 孙继高, 陈卫衡. 基于CT灰度值赋值的股骨头坏死有限元模型对比[J]. 中国组织工程研究, 2020, 24(3): 395-400.

Xue Zhipeng, Li Taixian, Li Yan, He Haijun, Huang Zeqing, Sun Jigao, Chen Weiheng. Comparison of finite element models of osteonecrosis of the femoral head based on CT gray-assigned method[J]. Chinese Journal of Tissue Engineering Research, 2020, 24(3): 395-400.

图/表（结果） 5

2.1 一般资料研究共纳入患者35例，其中男性20例，女性15例；双侧发病14例，左侧10例，右侧11例。纳入20髋经1年治疗稳定未塌陷的股骨头为实验组，20髋单侧股骨头坏死患者的健侧为正常组进行研究。两组患者年龄、侧别、体质量指数等一般临床资料比较见表1，两组患者年龄、侧别、体质量指数差异无显著性意义(P > 0.05)。

2.2 模型建立结果分别建立了20个正常股骨近端、不包括坏死骨的无坏死股骨近端和坏死骨的有限元模型，且不同模型的单元数与节点数均有不同，同一模型的数值较集中，有利于后续有限元分析的运算，具体见表2。

2.3 股骨头坏死力学传导分析结果对正常及坏死的股骨头模拟单足站立位进行应力分布(图4)及形变分析(图5)。云图显示：当施以体质量2.5倍压强作用于股骨头负重区表面时，正常股骨头负重区股骨头表面应力均匀，股骨头内应力延压力性骨小梁均匀分布在股骨头内；股骨头坏死显示股骨头负重区、坏死骨出现多个应力集中区域，股骨头表面应力分布不均匀，股骨头内应力分散于股骨头颈交界处，内外侧应力相当，且股骨头坏死产生的形变多于正常股骨头。

2.4 最大等效应力、最大总形变量比较三维有限元模型结果显示，股骨头坏死负重区较正常股骨头负重区比较最大等效应力与最大总形变量明显增大，差异有显著性意义(P < 0.05)，说明力学参数较真实的反映了股骨头的力学特性，且两组标准差偏于集中，数值稳定性较好，具体结果见表3。

参考文献

[1] LESPASIO MJ, SODHI N, MONT MA. steonecrosis of the Hip: A Primer.Perm J. 2019;23. pii: 18-100.
[2] 李子荣.股骨头坏死:早期诊断与个体化治疗[J].中国矫形外科杂志,2013,21(19):1909-1911.
[3] 叶奕亨,陈凯,金可可,等.塌陷前期股骨头坏死的保头手术治疗的进展[J].中国骨伤,2017,30(3): 287-292.
[4] 中国医师协会骨科医师分会显微修复工作委员会,中国修复重建外科专业委员会骨缺损及骨坏死学组,中华医学会骨科分会显微修复学组.成人股骨头坏死临床诊疗指南(2016)[J].中华骨科杂志,2016, 36(15):945-954.
[5] MOTOMURA G, YAMAMOTO T, YAMAGUCHI R, et al. Morphological analysis of collapsed regions in osteonecrosis of the femoral head.J Bone Joint Surg Br.2011;93-B(2):184-187.
[6] KARASUYAMA K, YAMAMOTO T, MOTOMURA G, et al. The role of sclerotic changes in the starting mechanisms of collapse: a histomorphometric and FEM study on the femoral head of osteonecrosis. Bone.2015;81:644-648.
[7] BROWN TD, PEDERSEN DR, BAKER KJ, et al. Mechanical consequences of core drilling and bone-grafting on osteonecrosis of the femoral head. J Bone Joint Surg Am.1993;75(9):1358.
[8] 何海军,陈卫衡,鲁超,等.基于CT图像建立可供动态力学模拟分析的股骨头坏死有限元模型[J].中国中医骨伤科杂志,2016, 24(9):5-10.
[9] 谢亮文.不同属性建模方法对股骨近端有限元分析影响的研究[D].福州:福建医科大学,2015. .
[10] 杨宾宾,刘耀升,刘蜀彬,等.多种髓芯减压术治疗股骨头坏死的有限元研究[J].中华损伤与修复杂志(电子版),2017,12(1):39-45.
[11] MONT MA, HUNGERFORD DS. Non-tranmatie avascular necrosis of the femoral head. J Bone Joint Surg(Am).1995; 77:459.
[12] GARDENIERS JWM.The ARCO perspective for reaching one uniform staging system of osteonecrosis//Bone circulation and vascularization in normal and pathological conditions. Springer,Boston,MA,1993:375-380.
[13] TANG Q, HE HJ, CHEN WH, et al. Mechanical evaluation of osteonecrosis of femoral head preand post-medical treatment: patient-specific computations.Neuroscience and Biomedical Engineering.2017;5(2):94-100.
[14] 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.
[15] PENG L, BAI J, ZENG X, et al. Comparison of isotropic and orthotropic material property assignments on femoral finite element models under two loading conditions. Med Eng Phys. 2006;28(3):227-233.
[16] 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.
[17] BROWN TD, WAY ME, FERGUSON AB JR. Mechanical characteristics of bone in femoral capital aseptic necrosis. Clin Orthop Relat Res.1981;(156):240-247.
[18] LI B, HU R, SUN L, et al. A CARE-compliant article: Biomechanics of treating early-stage femoral-head osteonecrosis by using a β-tricalcium phosphate bioceramic rod system: a 3-dimensional finite-element analysis. Medicine (Baltimore).2018;97(25):e10808.
[19] READ L, NIGG BM. An analysis of hip joint loading during walking, running, and skiing. Med Sci Sports Exerc.1999; 31(1):131-142.
[20] 李子荣.晚期股骨头坏死人工关节置换术的选择[J].中华关节外科杂志(电子版)[J].2008.2(1):51-52.
[21] ÜNAL MB, CANSU E, PARMAKSIZOĞLU F, et al. Treatment of osteonecrosis of the femoral head with free vascularized fibular grafting: results of 7.6-year follow-up. Acta Orthop Traumatol Turc. 2016;50(5):501-506.
[22] 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.
[23] BROWN TD, BAKER KJ, BRAND RA. Structural consequences of subchondral bone involvement in segmental osteonecrosis of the femoral head. J Orthop Res.1992;10(1): 79-87.
[24] SHI J, CHEN J, WU J, et al. Evaluation of the 3D finite element method using a tantalum rod for osteonecrosis of the femoral head. Med Sci Monit.2014;20:2556-2564.
[25] 庞智晖,魏秋实,周广全,等.个体股骨头坏死三维有限元模型的建立与应用[J].生物医学工程学杂志. 2012(2):251-255.
[26] 张海峰,宋翠荣,庞胤,等.基于CT数据建立髋关节三维有限元模型的研究[J].中国临床医学影像杂志, 2015,26(9):667-670.
[27] 彭军,胡冬发,邓立柳,等.股骨头坏死钽棒支撑植入术模拟及有限元模型的建立[J].中国数字医学,2014,9(4):66-68.
[28] TAYLOR WR, ROLAND E, PLOEG H, et al. Determination of orthotropic bone elastic constants using FEA and modal analysis. J Biomech.2002;35(6):767-773.
[29] 杨立峰,肖东民,彭春雷,等.自制股骨头内撑器治疗早期股骨头坏死的有限元分析[J].中国组织工程研究,2018,22(7):1096-1101.
[30] 李孝林,吕志鹏,刘迎军.股骨头坏死有限元模型的建立[J].中国组织工程研究与临床康复,2011,15(30): 5522-5525.
[31] CHEN L, HONG G, FANG B, et al. Predicting the collapse of the femoral head due to osteonecrosis: From basic methods to application prospects. J Orthop Translat.2017;11:62-72.
[32] KNAPIK JJ, REYNOLDS K, HOEDEBECKE KL. Stress Fractures: Etiology, Epidemiology, Diagnosis, Treatment, and Prevention. J Spec Oper Med. Summer.2017;17(2):120-130.
[33] MALEKIPOUR F, OETOMO D, LEE PV. Equine subchondral bone failure threshold under impact compression applied through articular cartilage. J Biomech. 2016;49(10): 2053-2059.
[34] ESCUDIER JC, OLLIVIER M, DONNEZ M, et al. Superimposition of maximal stress and necrosis areas at the top of the femoral head in hip aseptic osteonecrosis. Orthop Traumatol Surg Res. 2018;104(3):353-358.
[35] MEDIOUNI M, SCHLATTERER DR, KHOURY A, et al. Optimal parameters to avoid thermal necrosis during bone drilling: A finite element analysis. J Orthop Res.2017;35(11): 2386-2391.
[36] 史风雷,张美超.建立股骨头坏死有限元模型的方法[J].中国组织工程研究与临床康复,2008,12(30):5890-5891.
[37] BAE JY, KWAK DS, PARK KS, et al. Finite element analysis of the multiple drilling technique for early osteonecrosis of the femoral head. Ann Biomed Eng. 2013;41(12):2528-2537.
[38] 袁高翔,张伟滨.有限元分析在骨骼肌肉系统模型材料特性研究中的应用[J].国际骨科学杂志,2011, 32(6):352-355+366.
[39] 冀勇,胡永成,张春秋,等.股骨近端骨干组配式假体在下蹲步态下的三维有限元分析[J].天津理工大学学报,2018,34(3):12-16.
[40] 刘欣伟,闫寒,刘中洋,等.包含肌肉､韧带组织的骨盆､髋臼3D有限元模型的构建[J].临床军医杂志, 2014,42(4):331-335.
[41] 李泰贤. 健脾活骨方治疗非创伤性股骨头坏死生物力学机理与塌陷危险因素分析[D].北京:中国中医科学院,2018.

引言

材料方法

1.1 设计多样本计算机模拟仿真生物力学实验。

1.2 时间及地点于2017年6月至2019年1月在中国中医科学院望京医院及北京理工大学生命科学学院完成。

1.3 对象

1.3.1 病例来源选择在中国中医科学院望京医院骨关节三科就诊并接受系统治疗的成人股骨头坏死患者，共35例，选择20髋经1年治疗稳定未塌陷的股骨头为实验组，20髋单侧股骨头坏死患者的健侧为正常组。患者知晓全部过程自愿参加，并签署知情同意书。该研究方案经中国中医科学院望京医院医学伦理委员会批准。

1.3.2 诊断标准参照MONT等^[11]提出的股骨头坏死的诊断标准，分期参考国际骨循环学会(Association Research Circulation Osseous，ARCO)的标准^[12]。

1.3.3 纳入标准至少有一侧为ARCO II期的股骨头坏死；随访1年股骨头无塌陷者。

1.3.4 排除标准合并其他髋关节疾病；CT图像有噪点或无法拷贝。

1.4 材料软硬件设备：GE公司Bright speed 16排螺旋CT；联想(Lenovo ThinkStation)P510图形塔式工作站；图像处理软件为Mimics 16.0；逆向工程软件：Geomagic Studio 2012(64 bit)；有限元分析软件：Ansys 14.0。

1.5 方法

1.5.1 原始资料收集采用GE公司Bright speed 16排螺旋CT进行双髋关节骨组织窗扫描，层厚1.25 mm，层间隔为0.625 mm，得到二维CT图像，将获得的图像原始数据以DICOM 3.0标准保存。

1.5.2 建立模型髋关节CT DICOM 3.0图像的bim格式导入Mimics 16.0，基于CT灰度值进行阈值分割，界定阈值在226-1 672 Thresholdinig单位，修剪蒙板，擦除噪点，填充目标区域，建立股骨近端和坏死骨区域三维模型(图1)。此步骤需2名从事股骨头坏死研究的高年资医师共同依据X射线片、CT和MRI确认坏死区域。

1.5.3 生成网格模型 Geomagic Studio 2012软件导入model文件，利用软件工具将坏死骨三维模型与股骨近端做表面润滑操作。润滑操作需注意力度，以保持原型为原则。再次，利用布尔运算(Boolean Command)，将整个股骨头与坏死区作减法运算，得到不包括坏死骨的无坏死股骨近端三维模型与坏死骨三维模型。最后，应用ANSYS Mechanical APDL选择单元类型，进行网格划分，保存nodes节点数、elements单元数和prep7文件(图2)。

1.5.4 材料赋值与力学加载参考前期研究^[13-16]，根据CT灰度值-表观密度-计算弹性模量值，公式：ρ[kg/m³]= 0.046 4·HU+1 000，E[Pa]=3 790·ρ3 (ρ表示密度，E表示杨氏模量)。其中坏死骨赋予10级材料属性，股骨近端赋予30级材料属性^[17-18]。参考BROWN等^[7]研究结果设置股骨头的负重区，此区域位于股骨内外方向对向股骨头中心的弧度为40°，前后方向对向股骨头中心的弧度为80°，模拟单足站立位进行载荷设定和静态分析，施加压力以体质量的2.5倍进行模拟计算^[19]，见图3。

1.6 主要观察指标有限元模型节点数与单元数；实验组与正常组的应力分布情况及股骨头负重区的最大等效应力与最大总形变量的变化。

1.7 统计学分析采用SPSS 20.0对数据进行统计学分析，其中计量资料符合正态分布则将数据表示为x(_)±s，采用t 检验统计分析；不符合正态分布则用中位数代表平均水平，采用秩和检验；计数资料采用卡方检验或者秩和检验进行统计分析。P < 0.05则代表差异有显著性意义。

讨论

3.1 股骨头坏死有限元模型现状与对比股骨头坏死是致残率极高的骨伤科疑难病，未经治疗的患者80%股骨头会在1-4年内塌陷，不得不通过人工关节置换恢复关节功能，给社会、家庭带来了沉重的负担^[20]，因此，预防或延缓塌陷具有重要的研究意义。股骨头塌陷是生物学和生物力学共同作用的结果^[21-22]，除生物学研究之外，股骨头作为负重关节髋关节的组成部分，其生物力学的认识和研究对髋关节生理功能的体现尤为重要，因此有必要进行股骨头的生物力学研究，模拟真实股骨头坏死的受力情况，为股骨头坏死塌陷的防治提供数据基础。

有限元分析是分析股骨头坏死生物力学的主要研究方法之一，建立真实准确、稳定性好的有限元模型是其关键步骤^[23-24]。目前的研究中，不同研究者采用的方法和工具各不相同，包括利用几何模型、力学参数或是利用真实股骨头标本逆向建模，尽可能的还原股骨头的组织特性和形态不规则性。在影像资料选择方面，现阶段的有限元模型主要采用CT图像来提供几何形态信息和参数，股骨头区域分割包括区分皮质骨、松质骨赋值，不区分皮质骨、松质骨赋值，其中，区分皮质骨和松质骨的建模方法需准确区分皮质骨、松质骨和坏死核的范围；不区分皮质骨、松质骨的建模采用整体建模方法。有研究将股骨头分为软骨、软骨下坏死区、肉芽带、硬化带和正常骨小梁进行人为划分建模^[25]，但是人为划分区域在技术上显得相对复杂，人为误差不可避免，且需要大量时间操作，难以模拟多个坏死区及不同坏死区交叉的真实情况。

在研究对象的选择方面也存在差异。有研究选择健康成人股骨头作为研究对象^[18^，^26]，通过赋值及力学加载模拟股骨头坏死受力；有研究选择成人股骨头坏死患者^[27]；有研究利用尸体的股骨头解剖标本进行生物力学研究^[28]；有采用人工关节置换术后的新鲜股骨头标本作为对象进行CT扫描^[29]，其样本多是单例或较少样本。

该研究利用通过CT灰度值直接赋值的方法，区别于以往，利用20个正常股骨头与20个股骨头坏死的有限元模型进行多样本的对比研究，严格的区域分割、模拟真实情况进行力学加载，证明了有限元模型结果的准确性与稳定性。

3.2 建立有限元模型的关键环节有限元模型是研究股骨头坏死生物力学的重要基础，模型的精确性和稳定性是后续各种实验研究的关键，因此，建模过程应充分考虑模型的准确性和稳定性。

区域分割尽量减少人为因素参与^[30]。人为因素的参与为有限元模型带来了一定的偏倚，因此应在保证模型质量前提下尽量减少人为因素参与，最大程度的将人工处理与软件处理相结合。此次研究在前期分割皮质骨与松质骨建模的基础上，进一步减少人工处理，采用CT灰度值直接赋值，仅分割坏死区和股骨近端，其中坏死区的划分由2名从事股骨头坏死研究的高年资医生共同依据X射线片、CT、MRI检查综合决定。坏死区的组织分割细化尽管更精确的还原了其生物学特性，但增加了人为操作且技术复杂，需要大量时间，不利于临床多样本的重复研究，故此次研究未进一步细化坏死区。

建立具有良好的几何和物理相似性的有限元模型是进行精确力学分析的前提^[31]。此次研究在润滑过程中，严格控制润滑力度，多次尝试，以保证股骨头原始几何形态为准。另一方面，研究集中对比了20例正常股骨头与股骨头坏死的模型力学传导情况，正常股骨头负重区股骨头表面应力均匀，应力延压力性骨小梁均匀分布在股骨头内，股骨距承受大部分应力；而股骨头坏死在负重区、坏死骨区域出现应力集中区域，股骨头内应力分散于股骨颈内外侧。同时，实验组与正常组的云图显示股骨头坏死的形变更大，过大的形变超过一定范围得不到及时修复就会出现塌陷，这与股骨头坏死应力集中密切相关^[32-35]。

模型的稳定性是股骨头坏死的后续各项力学研究的基本保证^[36]。而大部分有限元模型的建立多是采用单样本或较少样本，为进行某种特定的力学研究建模，或研究生物支撑材料的力学效能，或不同手术间的模拟测试，尽管取得了一定的结果，但却无法保证模型的稳定性。此次研究将正常股骨头与坏死股骨头给予体质量2.5倍的压力进行模拟时，提取最大等效应力与最大总形变量参数，股骨头坏死与正常股骨头之间均具有统计学差异(P < 0.0.1)，且力学参数集中，模型稳定性较好，并与BAE等^[37]发表的研究结论接近。
3.3 实验设计的局限性为了保证模型的准确性和稳定性，建模过程的每一步进行了充分的文献研究与实践尝试，并通过对髋关节真实受力情况的模拟，对比了正常股骨头与股骨头坏死的应力分布，但此次研究也存在着一定的局限性。髋关节是负重关节，日常步态包括慢走、快走、慢跑等速度的变化，髋关节周围肌肉群的参与程度也不同，但此次研究只模拟了单足站立位进行载荷设定和静态分析。建模过程中，控制人为因素参与的程度，精确的分配赋值材料数目，髋臼与肌肉等组织结构的建模，多种步态的力学模拟^[38-41]，都是完善模型建立的必要研究环节。此次研究只针对坏死骨和股骨近端进行建模，尚未增加髋臼及肌肉等组织结构的建模；其次，根据赋值公式直接定义了坏死区材料属性的赋值，并未进一步分割细化。在下一步工作中，将细化建模步骤，在保证模型的准确性前提下，对各个关键环节进行更为细致的相关因素分析。

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

基于CT灰度值赋值的股骨头坏死有限元模型对比

Comparison of finite element models of osteonecrosis of the femoral head based on CT gray-assigned method

PDF

可视化

摘要/Abstract

引用本文

使用本文

图/表（结果） 5

参考文献

相关文章 15

引言

材料方法

讨论

文章快阅

延伸阅读

编辑推荐

Metrics

本文评价

[1]	李中峰, 陈明海, 凡一诺, 魏秋实, 何伟, 陈镇秋. 右归饮治疗激素性股骨头坏死作用机制的网络药理学分析[J]. 中国组织工程研究, 2021, 25(8): 1256-1263.
[2]	刘立华, 孙伟, 王云亭, 高福强, 程立明, 李子荣, 王江宁. 头颈部开窗减压治疗L1型激素性股骨头坏死：单中心前瞻性临床研究[J]. 中国组织工程研究, 2021, 25(6): 906-911.
[3]	刘钊, 徐西林, 申意伟, 张晓峰, 吕航, 赵军, 王政春, 刘旭卓, 王海涛. 塌陷预测方法联合分期分型对股骨头坏死治疗的指导作用与前景[J]. 中国组织工程研究, 2021, 25(6): 929-934.
[4]	李时斌, 赖渝, 周毅, 廖建钊, 章晓云, 张璇. 激素性股骨头坏死发病机制及相关信号通路的靶点效应[J]. 中国组织工程研究, 2021, 25(6): 935-941.
[5]	郑小龙, 何晓铭, 龚水帝, 庞凤祥, 杨帆, 何伟, 刘少军, 魏秋实. 酒精性股骨头坏死患者的骨转换特点[J]. 中国组织工程研究, 2021, 25(5): 657-661.
[6]	张国梅, 祝军, 胡杨, 焦红卫. E-Max瓷嵌体三维有限元模型粘接界面应力分析[J]. 中国组织工程研究, 2021, 25(4): 537-541.
[7]	曾祥洪, 梁博伟. 股骨头坏死保髋治疗的新策略[J]. 中国组织工程研究, 2021, 25(3): 431-437.
[8]	王岩, 董本超, 王颖, 孙磊, 卢斌, 柏豪豪, 田爱现, 马剑雄, 马信龙. 股骨头坏死动物模型的造模方法及特点[J]. 中国组织工程研究, 2021, 25(2): 292-297.
[9]	周毅, 陈跃平, 章晓云, 赖渝, 廖建钊, 李时斌. 复方中药圣愈汤治疗股骨头坏死机制的网络药理学分析[J]. 中国组织工程研究, 2021, 25(17): 2687-2696.
[10]	辛鹏飞, 柯梦楠 , 张海涛, 乌日莎娜, 李子祺, 庄至坤, 魏秋实, 何伟. 活血化瘀中药治疗股骨头坏死共同作用机制的网络药理学数据[J]. 中国组织工程研究, 2021, 25(17): 2727-2733.
[11]	刘蓬然, 焦瑞, 陶金, 陈辉, 戴纪杭, 颜连启. 不同界面假体全髋关节置换治疗老年髋关节疾病的比较[J]. 中国组织工程研究, 2021, 25(15): 2347-2351.
[12]	苏敬阳, 张耀杰, 曹斌, 郝学伟, 李晓, 韩永台. 多孔钽棒治疗ARCO Ⅰ-Ⅱ期非创伤性股骨头坏死后塌陷及相关因素分析[J]. 中国组织工程研究, 2021, 25(15): 2400-2404.
[13]	张煦坚, 赵振群, 刘万林. 激素性股骨头缺血坏死发病机制中的内质网应激[J]. 中国组织工程研究, 2021, 25(11): 1759-1765.
[14]	刘锌, 杜斌, 高丽丽, 孙光权, 王旭, 袁鹏, 林煊烨. 右归饮干预激素性股骨头坏死大鼠骨髓间充质干细胞的自噬及命运[J]. 中国组织工程研究, 2021, 25(1): 20-25.
[15]	王田田, 王建忠. 骨髓间充质干细胞治疗早期股骨头坏死的应用与前景 [J]. 中国组织工程研究, 2020, 24(7): 1117-1122.