不同材料赋值方法下踝关节三维有限元模型的应力及位移变化
|
摘要
背景:人体三维有限元仿真分析由于其自身的优越性,被广泛应用于人体生物力学研究,但对骨骼材料属性的赋值并未有确切的定论。目的:对踝关节进行力学有限元分析,探索出更符合临床实际情况的材料属性赋值方法。方法:对1例正常男性志愿者踝关节CT断层数据进行三维重建并进行曲面拟合,然后进行有限元网格划分,分别运用均一赋值法与灰度值赋值法2种方式对骨骼进行材料属性的赋予,模拟双腿站立、单足站立、正常步态、上楼梯、下楼梯等5种工况下踝关节的生物力学响应。结果与结论:①在5种工况下胫骨关节面、距骨体的应力峰值均一赋值组明显高于灰度值赋值组,2种赋值方法下腓骨关节面上的应力峰值相差不大;②均一赋值组在距骨颈处应力分布较为集中且分布区域面积较小,灰度值赋值组在距骨颈处应力分布较为均匀且分布区域面积覆盖整个距骨颈处;③均一赋值组形变量胫骨大于腓骨,灰度值赋值组腓骨大于胫骨;④均一赋值组胫骨与腓骨形变量的差值大于灰度值赋值组;⑤说明在对踝关节进行生物力学有限元分析中,灰度值赋值法更符合临床实际情况,具有更高的准确性且能在一定程度上解决由于个体差异性带来的建模误差。
BACKGROUND: Three-dimensional finite element simulation analysis of human body is widely used in human biomechanical research because of its own superiority, but the assignment of bone material properties has not been definitive.OBJECTIVE: To explore a more practical method of material attribute assignment by conducting finite element analysis of ankle joint biomechanics.METHODS: Three-dimensional reconstruction of the ankle joint CT data from a normal male volunteer was performed and underwent surface fitting. Then the finite element model was meshed, and the two methods of uniform assignment and gray value assignment were used to give the material properties of bones. The biomechanical response of the ankle joint in five working conditions: double leg standing, single foot standing, normal gait, upstairs and downstairs were simulated.RESULTS AND CONCLUSION:(1) In the five working conditions, the peak stress of body of talus in the uniformity group was obviously higher than that in the gray value assignment group. The stress value of the fibular joint surface was similar in the two methods.(2) In the uniform assignment group, the stress distribution was concentrated and the distribution area was narrow in the talar neck. The gray value assignment group had a uniform stress distribution at the talar neck and the distribution area covered the entire talar neck.(3) The tibial deformation was greater than fibula in the uniform assignment group, and the fibula deformation was greater than tibia in the gray value assignment group.(4) The uniform assignment group was larger than the gray value assignment group in the difference between the deformation of the tibia and the fibula.(5) In summary, in the finite element analysis of the ankle joint biomechanics, the gray value assignment method is more in line with the clinical actual situation, which has high accuracy and can solve the modeling error caused by the individual difference to some extent.
BACKGROUND: Three-dimensional finite element simulation analysis of human body is widely used in human biomechanical research because of its own superiority, but the assignment of bone material properties has not been definitive.OBJECTIVE: To explore a more practical method of material attribute assignment by conducting finite element analysis of ankle joint biomechanics.METHODS: Three-dimensional reconstruction of the ankle joint CT data from a normal male volunteer was performed and underwent surface fitting. Then the finite element model was meshed, and the two methods of uniform assignment and gray value assignment were used to give the material properties of bones. The biomechanical response of the ankle joint in five working conditions: double leg standing, single foot standing, normal gait, upstairs and downstairs were simulated.RESULTS AND CONCLUSION:(1) In the five working conditions, the peak stress of body of talus in the uniformity group was obviously higher than that in the gray value assignment group. The stress value of the fibular joint surface was similar in the two methods.(2) In the uniform assignment group, the stress distribution was concentrated and the distribution area was narrow in the talar neck. The gray value assignment group had a uniform stress distribution at the talar neck and the distribution area covered the entire talar neck.(3) The tibial deformation was greater than fibula in the uniform assignment group, and the fibula deformation was greater than tibia in the gray value assignment group.(4) The uniform assignment group was larger than the gray value assignment group in the difference between the deformation of the tibia and the fibula.(5) In summary, in the finite element analysis of the ankle joint biomechanics, the gray value assignment method is more in line with the clinical actual situation, which has high accuracy and can solve the modeling error caused by the individual difference to some extent.
引文
[1]马信龙.踝关节骨折的损伤机制[J].中华骨科杂志,2013,33(4):429-432.
[2]关国平,孟强,牛文鑫,等.踝关节扭伤防护及其康复研究进展[J].医用生物力学,2016,31(1):78-82.
[3]曾茂兰,王静.踝关节运动损伤应急处理及其康复治疗[J].甘肃医药,2017,36(1):25-26.
[4]张琳琳.人体上肢生物力学建模和典型运动的生物力学研究[D].上海:上海交通大学,2009.
[5]唐刚.人体典型运动生物力学仿真分析[D].上海:上海交通大学,2011.
[6]王博文.标本含湿率及加载方式对骨微观力学特性测试的影响研究[D].西安:第四军医大学,2014.
[7]郑海峰.髂骨新月形骨折伴骶髂关节脱位两种内固定方法的生物力学研究[D].大连:大连医科大学,2016.
[8]何旅洋,张志强,郑百林.外骨骼框架生物力学设计及强度分析[J].医用生物力学,2014,29(6):504-510.
[9]尹峰,王晓东,梁炜,等.建立锁骨个性化锁定接骨板模型及有限元分析[J].中国组织工程研究,2016,20(35):5244-5249.
[10]Zhou EC,Tang P,Zhu CY,et al.Finite element analysis on the effect of lateral wedge insole intervention on the contact characteristics of the subtalar joint.Chin J Orthop Traumatol.2017;30(1):38.
[11]Wang X,Sanyal A,Cawthon PM,et al.Prediction of new clinical vertebral fractures in elderly men using finite element analysis of CT scans.J Bone Mineral Res.2012;27(4):808-816.
[12]欧阳汉斌,谢普生,邓羽平,等.基于CT灰度赋值的有限元模型建立及其在截骨矫形中的应用[J].南方医科大学学报,2016,36(7):979-983.
[13]Dall'Ara E,Luisier B,Schmidt R,et al.A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro.Bone.2013;52(1):27-38.
[14]李瑜,邓双城,杜程,等.中立位踝关节生物力学有限元仿真[J].北京石油化工学院学报,2017,25(4):51-58.
[15]张宏宁,朱永展,李雪,等.踝关节骨折中外踝移位对胫距关节接触影响的三维有限元分析[J].中华创伤骨科杂志,2017,19(9):769-774.
[16]郭国新,郭继涛,李伟,等.基于有限元模型的踝关节生物力学分析[J].中国组织工程研究,2012,16(17):3056-3060.
[17]刘清华,马建兵,姚建锋.不同载荷对正常踝关节影响的有限元研究[J].中华创伤骨科杂志,2013,15(8):704-708.
[18]李瑜,邓双城,杜程,等.中立位踝关节生物力学有限元仿真中立位踝关节生物力学有限元仿真[J].北京石油化工学院学报,2017,25(4):51-58.
[19]Chen G,Wu FY,Liu ZC,et al.Comparisons of node-based and element-based approaches of assigning bone material properties onto subject-specific finite element models.Med Eng Phys.2015;37(8):808.
[20]徐志才,胡广洪,黄振宇,等.胫骨模型对膝关节有限元分析结果影响的探讨[J].中国数字医学,2014,9(4):69-72.
[21]谢亮文.不同属性建模方法对股骨近端有限元分析影响的研究[D].福州:福建医科大学,2015.
[22]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.
[23]秦林林,马海波,张卫,等.中国北方汉族健康人骨密度正常值[J].中国骨质疏松杂志,2002,8(2):110-111.
[24]Cheung JT,An KN,Zhang M.Consequences of partial and total plantar fascia release:a finite element study.Foot Ankle Int.2006;27(2):125.
[25]Simon SR,Paul IL,Mansour J,et al.Peak dynamic force in human gait.J Biomech.1981;14(12):817,821-819,822.
[26]Stacoff A,Diezi C,Luder G,et al.Ground reaction forces on stairs:effects of stair nclination and age.Gait Posture.2005;21(1):24.
[27]Anderson DD,Goldsworthy JK,Li W,et al.Physical validation of a patient-specific contact finite element model of the ankle.JBiomech.2007;40(8):1662-1669.
[28]唐勇涛,魏思奇,程刚,等.基于两种不同建模方式的胫骨有限元对比分析[J].中国骨科临床与基础研究杂志,2016,8(1):26-32.
[29]尚禹,彭亮,刘加成,等.均匀材料和基于CT灰度值材料的股骨、胫骨有限元分析[J].中国生物医学工程学报,2008,27(5):722-727.
[2]关国平,孟强,牛文鑫,等.踝关节扭伤防护及其康复研究进展[J].医用生物力学,2016,31(1):78-82.
[3]曾茂兰,王静.踝关节运动损伤应急处理及其康复治疗[J].甘肃医药,2017,36(1):25-26.
[4]张琳琳.人体上肢生物力学建模和典型运动的生物力学研究[D].上海:上海交通大学,2009.
[5]唐刚.人体典型运动生物力学仿真分析[D].上海:上海交通大学,2011.
[6]王博文.标本含湿率及加载方式对骨微观力学特性测试的影响研究[D].西安:第四军医大学,2014.
[7]郑海峰.髂骨新月形骨折伴骶髂关节脱位两种内固定方法的生物力学研究[D].大连:大连医科大学,2016.
[8]何旅洋,张志强,郑百林.外骨骼框架生物力学设计及强度分析[J].医用生物力学,2014,29(6):504-510.
[9]尹峰,王晓东,梁炜,等.建立锁骨个性化锁定接骨板模型及有限元分析[J].中国组织工程研究,2016,20(35):5244-5249.
[10]Zhou EC,Tang P,Zhu CY,et al.Finite element analysis on the effect of lateral wedge insole intervention on the contact characteristics of the subtalar joint.Chin J Orthop Traumatol.2017;30(1):38.
[11]Wang X,Sanyal A,Cawthon PM,et al.Prediction of new clinical vertebral fractures in elderly men using finite element analysis of CT scans.J Bone Mineral Res.2012;27(4):808-816.
[12]欧阳汉斌,谢普生,邓羽平,等.基于CT灰度赋值的有限元模型建立及其在截骨矫形中的应用[J].南方医科大学学报,2016,36(7):979-983.
[13]Dall'Ara E,Luisier B,Schmidt R,et al.A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro.Bone.2013;52(1):27-38.
[14]李瑜,邓双城,杜程,等.中立位踝关节生物力学有限元仿真[J].北京石油化工学院学报,2017,25(4):51-58.
[15]张宏宁,朱永展,李雪,等.踝关节骨折中外踝移位对胫距关节接触影响的三维有限元分析[J].中华创伤骨科杂志,2017,19(9):769-774.
[16]郭国新,郭继涛,李伟,等.基于有限元模型的踝关节生物力学分析[J].中国组织工程研究,2012,16(17):3056-3060.
[17]刘清华,马建兵,姚建锋.不同载荷对正常踝关节影响的有限元研究[J].中华创伤骨科杂志,2013,15(8):704-708.
[18]李瑜,邓双城,杜程,等.中立位踝关节生物力学有限元仿真中立位踝关节生物力学有限元仿真[J].北京石油化工学院学报,2017,25(4):51-58.
[19]Chen G,Wu FY,Liu ZC,et al.Comparisons of node-based and element-based approaches of assigning bone material properties onto subject-specific finite element models.Med Eng Phys.2015;37(8):808.
[20]徐志才,胡广洪,黄振宇,等.胫骨模型对膝关节有限元分析结果影响的探讨[J].中国数字医学,2014,9(4):69-72.
[21]谢亮文.不同属性建模方法对股骨近端有限元分析影响的研究[D].福州:福建医科大学,2015.
[22]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.
[23]秦林林,马海波,张卫,等.中国北方汉族健康人骨密度正常值[J].中国骨质疏松杂志,2002,8(2):110-111.
[24]Cheung JT,An KN,Zhang M.Consequences of partial and total plantar fascia release:a finite element study.Foot Ankle Int.2006;27(2):125.
[25]Simon SR,Paul IL,Mansour J,et al.Peak dynamic force in human gait.J Biomech.1981;14(12):817,821-819,822.
[26]Stacoff A,Diezi C,Luder G,et al.Ground reaction forces on stairs:effects of stair nclination and age.Gait Posture.2005;21(1):24.
[27]Anderson DD,Goldsworthy JK,Li W,et al.Physical validation of a patient-specific contact finite element model of the ankle.JBiomech.2007;40(8):1662-1669.
[28]唐勇涛,魏思奇,程刚,等.基于两种不同建模方式的胫骨有限元对比分析[J].中国骨科临床与基础研究杂志,2016,8(1):26-32.
[29]尚禹,彭亮,刘加成,等.均匀材料和基于CT灰度值材料的股骨、胫骨有限元分析[J].中国生物医学工程学报,2008,27(5):722-727.