[1]方斌,何伟,展磊,等.不同坏死范围下股骨头坏死区应力分布的有限元分析[J].中医正骨,2012,24(10):10-15.
　FANG Bin*,HE Wei,ZHAN Lei,et al.Finite element analysis of stress distribution over femoral head necrosis zones in different necrosis areas[J].The Journal of Traditional Chinese Orthopedics and Traumatology,2012,24(10):10-15.
点击复制

不同坏死范围下股骨头坏死区应力分布的有限元分析()
分享到：

《中医正骨》[ISSN:1001-6015/CN:41-1162/R]

卷:
第24卷
期数:
2012年10期
页码:
10-15
栏目:
基础研究
出版日期:
2012-10-20

文章信息/Info

Title:
Finite element analysis of stress distribution over femoral head necrosis zones in different necrosis areas
作者:
方斌1; 何伟1; 展磊2; 张庆文1; 魏秋实3; 陈镇秋1; 陈鹏1; 张美超4
1.广州中医药大学第一附属医院,广东 广州 510407;
2.广东药学院附属新市医院, 广东 广州 510410;
3.广州军区广州总医院,广东 广州 510010;
4.南方医科大学,广东 广州 510515
Author(s):
FANG Bin*; HE Wei; ZHAN Lei; ZHANG Qing-wen; WEI Qiu-shi; CHEN Zhen-qiu; CHEN Peng; ZHANG Mei-chao.
*The First Affiliated Hospital of Guangzhou University of Traditional Chinese Medicine,Guangzhou 510407,Guangdong,China
关键词:
股骨头坏死 坏死范围 应力分布 塌陷 有限元分析
Keywords:
Femur head necrosis;  Necrosis area;  Stress distribution;  Collapse;  Finite element analysis
摘要:
目的:探讨股骨头坏死范围对坏死区应力分布的影响。方法:采用中国虚拟人男1号数据集,利用Mimics软件分别建立坏死组织体积分数为30%、50%和60%的股骨头坏死有限元模型,并利用Ansys软件进行分析。对上述股骨头坏死模型分别施以1 440 N、2 400 N、4 200 N载荷,根据不同载荷及坏死范围,共设置9个股骨头坏死模型,分别测定坏死区表面和底部的应力分布情况、应力峰值,并分析股骨头塌陷的可能性。结果:①股骨头坏死区表面应力分布情况。相同坏死范围下,股骨头在不同载荷作用下,坏死区表面应力分布基本一致,载荷越大,相应区域的应力值越大; 坏死组织体积分数为30%时,应力分布越接近活骨组织其应力值越大,表面峰值应力区位于死骨与活骨交界边缘的外后方; 坏死组织体积分数为50%、60%时,股骨头坏死表面范围超过了股骨头承重面,二者应力分布类似,股骨头峰值应力区位于股骨头后外侧,坏死股骨头表面后方出现多点、多区的应力集中,前外侧则出现1条相对于周围应力增高的应力带。②股骨头坏死区底部应力分布情况。相同坏死范围下,股骨头在不同载荷作用下,坏死区底部的应力分布基本一致,载荷越大,相应区域的应力值越大; 坏死组织体积分数为30%时,坏死区底部应力分布以外侧、后侧为主,越接近活骨组织其应力值越大,峰值应力区位于死骨与活骨交界边缘的后方; 坏死组织体积分数为50%、60%时不同载荷下股骨头坏死区底部应力分布类似,外侧及后侧存在明显的多点、多区应力增高,峰值应力区位于股骨头坏死区基底部的外侧,基底部前方存在散点应力增高区。③股骨头坏死区表面与底部的峰值应力。股骨头坏死区底部的应力明显较股骨头坏死区表面高; 在4 200 N载荷作用下不同坏死范围的3种模型及2 400 N载荷作用下坏死组织体积分数为30%的坏死模型的坏死区底部的峰值应力(0.992 MPa、0.685 MPa、0.692 MPa、0.567 MPa)均超过坏死骨组织的临界应力(0.55 MPa)。④不同坏死范围下股骨头塌陷的可能性。9个坏死股骨头模型仅坏死区底部峰值应力点及附近小部分区域超过临界应力,其余大部分区域均未超过临界应力,即所有模型中坏死组织仅发生应力代偿不全性微骨折。结论:坏死范围影响股骨头坏死区的应力分布; 坏死区底部的应力高于坏死区表面; 应力峰值超过临界值时出现软骨下深部松质骨微骨折,是导致股骨头塌陷的直接原因。
Abstract:
Objective:To explore the effect of femoral head necrosis areas on the stress distribution over the femoral head necrosis zones.Methods:Data set of Virtual Chinese Human male No.1 was adopted,and the finite element models of femoral head necrosis with necrotic tissue volume fractions of 30%,50% and 60% were established through Mimics software respectively,and then they were analysed by Ansys software.These models of femoral head necrosis were performed with loads of 1 440 N,2 400 N and 4 200 N respectively,and 9 models of femoral head necrosis were set in total according to the different loads and different necrosis areas.Then the stress distribution and peak stress on surface and bottom of necrosis zones were measured respectively,and the possibility of femoral head collapse was analysed.Results:Under the same necrosis areas,stress distribution over the surface of femoral head necrosis zones were basically in the same although femoral head bore different loads.The larger the load was,the greater the stress value was in the corresponding areas.When necrotic tissue volume fraction reached 30%,the closer the stress distribution came to that of live bone tissue,the greater the corresponding stress value was,and the surface peak stress zone was located at the posterolateral approach to the edge of boundary of dead bones and live bones.When necrotic tissue volume fraction reached 50% or 60%,the similar stress distribution were shown and the scope of surface of femoral head necrosis was larger than that of weight-bearing surface of femoral head.The peak stressed zone of femoral head was in the posterolateral approach to femoral head,also the multipoint and multiregion stress concentration appeared in the rear surface of necrotic femoral head,while one stress zone with relative higher stress compared with that of ambient stress was shown on anterolateral femoral head.Under the same necrosis areas,stress distribution over the bottom of femoral head necrosis zones were basically in the same although femoral head bore different loads.The larger the load was,the greater the stress value was in the corresponding areas.When necrotic tissue volume fraction reached 30%,the stress was mainly distributed over the lateral side and rear side of bottom of necrosis zone,and the closer it came to live bone tissue,the greater its stress value was,and the peak stress zone was located behind the edge of boundary of dead bones and live bones.When necrotic tissue volume fraction reached 50% or 60%,the similar stress distribution over the bottom of femoral head necrosis zone were shown under different loads,and the stress increased significantly in multipoint and multiregion of the lateral and rear side,while peak stress zone was located at the lateral side of basilar part of femoral head necrosis zones and sporadic stress concentrated area was located at the front of basilar part.The stress on bottom of femoral head necrosis zone was significantly higher than that on the surface of femoral head necrosis zone.For the three kinds of models under 4 200 N load in different necrosis areas and necrosis model with necrotic tissue volume fraction of 30% under 2 400 N load,the peak stress(0.992 MPa,0.685 MPa,0.692 MPa,0.567 MPa)on bottom of the necrosis zones were all larger than the critical stress(0.55 MPa)of necrotic bone tissue.For the 9 models of necrotic femoral head,only the peak stress on the bottom of necrosis zones and nearby small areas were larger than the critical stress,while the peak stress on most of the rest areas were all smaller than the critical stress,namely,for all the models,the tissue necrosis was only presented as stress decompensation microfracture.Conclusion:Necrosis area can affect the stress distribution over femoral head necrosis zone and stress on the bottom of necrosis zone is higher than that on the surface of necrosis zone.The microfracture of subchondral cancellous bone in the deep appear when peak stress is larger than the critical stress,which directly leads to the collapse of femoral head.

参考文献/References:

[1] Ha YC,Jung WH,Kim JR,et al.Prediction of collapse in femoral head osteonecrosis:a modified Kerboul method with use of magnetic resonance images[J].J Bone Joint Surg Am,2006,88(Suppl 3):35-40.
[2] Hernigou P,Lambotte JC.Volumetric analysis of osteonecrosis of the femur.Anatomical correlation using MRI[J].J Bone Joint Surg Br,200l,83(5):672-675.
[3] Nishii T,Sugano N,Ohzono K,et al.Significance of lesion size and location in the prediction of collapse of osteonecrosis of the femoral head:a new three-dimensional quantification using magnetic resonance imaging[J].J Orthop Res,2002,20(1):130-136.
[4] Yang JW,Koo KH,Lee MC,et al.Mechanics of femoral head osteonecrosis using three-dimensional finite element method[J].Arch Orthop Trauma Surg,2002,122(2):88-92.
[5] Etienne G,Mont MA,Ragland PS.The diagnosis and treatment of nontraumatic osteonecrosis of the femoral head[J].Instr Course Lect,2004,53:67-85.
[6] Steinberg DR,Steinberg ME,Garino JP,et al.Determining lesion size in osteonecrosis of the femoral head[J].J Bone Joint Surg Am,2006,88(Suppl 3):27-34.
[7] Lee MS,Tai CL,Senan V,et al.The effect of necrotic lesion size and rotational degree on the stress reduction in Transtrochanteric rotational osteotomy for femoral head osteonecrosis a three-dimensional finite-element simulation[J].Clin Biomech(Bristol,Avon),2006,21(9):969-976.
[8] Floerkemeier T,Lutz A,Nackenhorst U,et al.Core decompression and osteonecrosis intervention rod in osteonecrosis of the femoral head:clinical outcome and finite element analysis[J].Int Orthop,2011,35(10):1461-1466.
[9] 史风雷,张美超.建立股骨头坏死有限元模型的方法[J].中国组织工程研究与临床康复,2008,12(30):5890-5891.
[10] Brown TD,Ferguson AB Jr.Mechanical property distributions in the cancellous bone of the human proximal femur[J].Acta Orthop Scand,1980,51(3):429-437.
[11] Brown TD,Pedersen DR,Baker KJ,et al.Mechanical consequences of core drilling and bone-grafting on osteonecrosis of the femoral head[J].J Bone Joint Surg Am,1993,75(9):1358-1367.
[12] Chen WP,Tai CL,Tan CF,et al.The degrees to which Iranstrochanteric rotational estcotomy moves the region of osteonecrotic femoral head out of the weight-bearing area as evaluated by computer simulation[J].Clin Biomech(Bristol,Avoll),2005,20(1):63-69.
[13] Yoshida H,Faust A,Wilckens J,el al.Three-dimensional dynamic hip contact area and pressure distribution during activities of daily living[J].J Biomech,2006,39(11):l996-2004.
[14] Langlais F,Fourastier J.Rotation esteotomies for ogleonecresis of the femoral head[J].Clin Orthop Relat Res,1997,(343):110-123.

备注/Memo

备注/Memo:
通讯作者:何伟 E-mail:hewei1123@21cn.com
更新日期/Last Update: 2012-10-20