基于断裂力学模拟3种不同跌倒状态下股骨近端骨折的有限元分析
|
摘要
目的运用断裂力学分析方法,探究3种不同边界条件下模拟跌倒的最佳有限元模型。方法将1例老年股骨颈骨折健侧髋关节以及大腿螺旋CT的DICOM格式影像资料导入Mimics中三维重建出股骨模型,将模型导入Hypermesh中设置材料属性,定义3种不同的仿真跌倒状态下的边界条件:股骨长轴与水平线成30°同时股骨颈施加垂直载荷(A)、股骨长轴与水平线成10°同时股骨颈施加垂直载荷(B)、股骨长轴与水平线成25°同时垂直股骨颈施加水平载荷(C),采用LS-DYNA软件对3种模型运算分析。结果相同大小的载荷下,断裂力学分析A组模型出现股骨颈骨折、B组股骨大结节处骨折、C组股骨干上段骨折;Von Mises应力云图分析:A组模型的应力主要集中在股骨颈部位,股骨转子间以及股骨干应力很小;B组模型中主要应力集中于股骨大结节处(也是有限元模型中固定限制的节点),C组模型应力集中在股骨颈以及股骨干部位。结论基于Hypermesh/LS-DYNA断裂力学软件环境下,A组边界条件的模型仿真跌倒更具有可信度。
Objective To explore the optimal finite element model of simulated fall under three different boundary conditions by using the fracture mechanics analysis method.Methods A DICOM format image data of an elderly femoral neck fracture and thigh spiral CT was imported into Mimics to reconstruct the femor model three-dimensionally.And thenthe model was introduced into Hypermesh to define the material properties and the three different boundary conditions.The details were as follows:the longitudinal axis of the femor was 30°while the femoral neck was applied with a vertical load(A),the femoral long axis was at 10° to the horizontal line and the femoral neck was applied with a vertical load(B),the femoral long axis was at 25° and the vertical femoral neck was applied horizontally Load(C).The three models were analyzed by using the LS-DYNA software.Results Under the same magnitude force of the load,the fracture mechanics analysis of group A model turned to be femoral neck fracture,group B showed femoral nodular fracture and group C was femoral shaft fracture;Von Mises stress cloud analysis:The stress of group A model was mainly concentrated in the femoral neck since femoral intertrochanter and femoral shaft stress is very small;the stress ofgroup B model mainlywasconcentrated in the femor nodules(also fixed limit nodule in finite element model);the stress of the group C was concentrated in the femoral neck and the femoral stem.Conclusion The simulation of the boundary condition of group A is more credible,based on the Hypermesh/LS-DYNA fracture mechanics software environment analysis.
Objective To explore the optimal finite element model of simulated fall under three different boundary conditions by using the fracture mechanics analysis method.Methods A DICOM format image data of an elderly femoral neck fracture and thigh spiral CT was imported into Mimics to reconstruct the femor model three-dimensionally.And thenthe model was introduced into Hypermesh to define the material properties and the three different boundary conditions.The details were as follows:the longitudinal axis of the femor was 30°while the femoral neck was applied with a vertical load(A),the femoral long axis was at 10° to the horizontal line and the femoral neck was applied with a vertical load(B),the femoral long axis was at 25° and the vertical femoral neck was applied horizontally Load(C).The three models were analyzed by using the LS-DYNA software.Results Under the same magnitude force of the load,the fracture mechanics analysis of group A model turned to be femoral neck fracture,group B showed femoral nodular fracture and group C was femoral shaft fracture;Von Mises stress cloud analysis:The stress of group A model was mainly concentrated in the femoral neck since femoral intertrochanter and femoral shaft stress is very small;the stress ofgroup B model mainlywasconcentrated in the femor nodules(also fixed limit nodule in finite element model);the stress of the group C was concentrated in the femoral neck and the femoral stem.Conclusion The simulation of the boundary condition of group A is more credible,based on the Hypermesh/LS-DYNA fracture mechanics software environment analysis.
引文
[1]GIBSON A A,HAY A W,RAY D C.Patients with hip fracture admitted to critical care:epidemiology,interventions and outcome[J].Injury,2014,45(7):1066-1070.
[2]范璐,林华,陈新,等.跌倒和髋部骨密度对绝经后妇女骨折风险的影响[J].中华骨质疏松和骨矿盐疾病杂志,2011,4(4):241-246.
[3]JARVINEN T L,SIEVANEN H,KHAN K M,et al.Shifting the focus in fracture prevention from osteoporosis to falls[J].BMJ,2008,336(7636):124-126.
[4]BESSHO M,OHNISHI I,OKAZAKI H,et al.Prediction of the strength and fracture location of the femoral neck by CT based finite-element method:apreliminary study on patients with hip fracture[J].J Orthop Sci,2004,9(6):545-550.
[5]SUTTER E G,WALL S J,MEARS S C,et al.The effect of cement placement on augmentation of the osteoporotic proximal femur[J].Geriatr Orthop Surg Rehabil,2010,1(1):22-26.
[6]CODY DD,GROSS GJ,HOU FJ,et al.Femoral strength is better predicted by finite element models than QCT and DXA[J].J Biomech,1999,32(10):1013-1020.
[7]李鹏飞,杜根发,林梓凌,等.基于LS-DYNA模拟老年股骨颈骨折的有限元分析[J].中国组织工程研究,2016,20(44):6606-6611.
[8]韩勇.车辆与行人碰撞中下肢动力学响应和损伤生物力学分析[D].长沙:湖南大学,2006.
[9]HOBATHO M C,RHO J Y,ASHMAN R B.Anatomical variation of human cancellous bone mechanical properties in vitro[J].Stud Health Technol Inform,1997,40(40):157-173.
[10]张国栋,廖维靖,陶圣祥,等.股骨有限元分析赋材料属性的方法[J].中国组织工程研究与临床康复,2009,13(43):8436-8441.
[11]孙培栋.侧方跌倒高度及髋保护器对髋部冲击影响的实验及有限元分析[D].广州:南方医科大学,2012.
[12]康颖安.断裂力学的发展与研究现状[J].科技致富向导,2014,16(24):39-42.
[13]薛付霞.断裂力学在边坡稳定分析中的应用研究[D].重庆:重庆交通大学,2014.
[14]薛明琛,赵永生,汤美安.断裂力学和损伤力学在混凝土中的应用[J].山西建筑,2007,33(36):19-20.
[15]范小宁,徐格宁,杨瑞刚.基于损伤-断裂力学理论的起重机疲劳寿命估算方法[J].中国安全科学学报,2011,21(9):58-63.
[16]胡建军,许洪斌,金艳,等.基于有限元计算的金属断裂准则的应用与分析[J].锻压技术,2007,32(3):100-104.
[17]李平.不同跌倒姿势与髋部应力变化关系的有限元分析[D].广州:广州中医药大学,2014.
[2]范璐,林华,陈新,等.跌倒和髋部骨密度对绝经后妇女骨折风险的影响[J].中华骨质疏松和骨矿盐疾病杂志,2011,4(4):241-246.
[3]JARVINEN T L,SIEVANEN H,KHAN K M,et al.Shifting the focus in fracture prevention from osteoporosis to falls[J].BMJ,2008,336(7636):124-126.
[4]BESSHO M,OHNISHI I,OKAZAKI H,et al.Prediction of the strength and fracture location of the femoral neck by CT based finite-element method:apreliminary study on patients with hip fracture[J].J Orthop Sci,2004,9(6):545-550.
[5]SUTTER E G,WALL S J,MEARS S C,et al.The effect of cement placement on augmentation of the osteoporotic proximal femur[J].Geriatr Orthop Surg Rehabil,2010,1(1):22-26.
[6]CODY DD,GROSS GJ,HOU FJ,et al.Femoral strength is better predicted by finite element models than QCT and DXA[J].J Biomech,1999,32(10):1013-1020.
[7]李鹏飞,杜根发,林梓凌,等.基于LS-DYNA模拟老年股骨颈骨折的有限元分析[J].中国组织工程研究,2016,20(44):6606-6611.
[8]韩勇.车辆与行人碰撞中下肢动力学响应和损伤生物力学分析[D].长沙:湖南大学,2006.
[9]HOBATHO M C,RHO J Y,ASHMAN R B.Anatomical variation of human cancellous bone mechanical properties in vitro[J].Stud Health Technol Inform,1997,40(40):157-173.
[10]张国栋,廖维靖,陶圣祥,等.股骨有限元分析赋材料属性的方法[J].中国组织工程研究与临床康复,2009,13(43):8436-8441.
[11]孙培栋.侧方跌倒高度及髋保护器对髋部冲击影响的实验及有限元分析[D].广州:南方医科大学,2012.
[12]康颖安.断裂力学的发展与研究现状[J].科技致富向导,2014,16(24):39-42.
[13]薛付霞.断裂力学在边坡稳定分析中的应用研究[D].重庆:重庆交通大学,2014.
[14]薛明琛,赵永生,汤美安.断裂力学和损伤力学在混凝土中的应用[J].山西建筑,2007,33(36):19-20.
[15]范小宁,徐格宁,杨瑞刚.基于损伤-断裂力学理论的起重机疲劳寿命估算方法[J].中国安全科学学报,2011,21(9):58-63.
[16]胡建军,许洪斌,金艳,等.基于有限元计算的金属断裂准则的应用与分析[J].锻压技术,2007,32(3):100-104.
[17]李平.不同跌倒姿势与髋部应力变化关系的有限元分析[D].广州:广州中医药大学,2014.