肘关节有限元模型模拟碰撞损伤的生物力学特性分析

来源：

唐鹏程 ，朱 英 ，赵敏珠 ，代亚磊 ，张 丽  ，李红卫 ，李剑波 ， 400016 重庆，重庆医科大学：基础医学 院法 医学教研 室 ，重庆市刑事侦查工程技术研究 中心 ；400021 重庆 ，重庆市公安局刑警总队技术处

[摘要]

目的 利用有限元法建立肘关节的三维有限元模型,动态模拟肘关节以不同角度撞击地面,分析肘关节的生物力学特性｡

方法 选取1名健康青年男性志愿者,肘关节屈曲90°行CT扫描,应用Mimics､Geomagic､UG､Hypermesh等软件对所得数据进行三维重建,建立肘关节的三维有限元模型;在Abaqus软件中模拟肘关节以不同角度(前臂与地面呈30°､60°､90°)撞击地面时,计算骨皮质破裂的临界速度及显著破裂时的速度,分析应力传导及其大小､分布等的变化规律｡

结果 肘关节屈曲90°状态发生正面碰撞时应力传导主要集中在肱骨髁､肱骨滑车及尺骨鹰嘴等部位,而桡骨头应力传导较小;当肘关节与地面撞击(前臂与地面分别呈30°､60°､90°)时,骨皮质破裂临界速度分别为11.4､9.3､13.8nv/s;前臂与地面分别呈30°､60°撞击,撞击速度分别达15.0､13.0m/s时,尺骨鹰嘴呈粉碎性骨折;前臂与地面呈90°撞击,撞击速度达18.0m/s时,肱骨远端呈明显横行骨折｡

结论 应用有限元建模分析软件构建的人体肘关节有限元模型可用于外力作用下肘关节损伤机制的研究,也可用于肘关节损伤时暴力的大小､速度､方向的推断｡

[关键词]肘关节;有限元分析;损伤;生物力学

肘关节由肱骨下端和尺骨、桡骨上端构成,包括3个关节,即肱尺关节肱桡关节和桡尺近侧关节。肘关节可作前屈、后伸运动也参与了前臂的旋前和旋后运动,是人日常生活中最重要的运动关节之一。肘关节因关节囊的保护,其损伤相较于其他部分并不常见。主要发生在交通事故、暴力斗殴等意外事故中。以往对肘关节的研究集中在其解剖结构及韧带的形态功能上,而对肘关节损伤的机制缺乏系统、深入的研究。同时,由于关节的解结构较为复杂,在医学实践中对其损伤为直接还是间接暴力形成,损伤的发生与受力方向、方式和力的作用面的关系,以及骨折线走向与受力方向的关联程度等的判断仍是难点和重点。因此针对关节损伤机制的研究具有重要的法医学和临床医学价值。近年来,随着有限元技术的快速发展,有限元模拟实验已广泛应用于医学的各个领域48。有限元技术可以将复杂的实体模型简化为数字模型并对这些模型进行计算机软件分析以替代复杂昂贵乃至无法实现的实验。因此,本研究应用有限元法结合高分辨率CT及仿真模拟软件建立完善的肘关节模型模拟时关节以不同角度碰撞地面,探讨肘关节应力和应变等生物力学性能以期更深远地了解肘关节的生物力学特性研究肘关节碰撞时损伤形成的机制。

1材料与方法

1.1研究对象

青年男性志愿者1名26岁,身高171cm,体质量

68kg,既往无肘关节前臂等损伤及相关疾病史,并通过肘关节及前臂的X线CT及MRI排除相关疾病。本研究经重庆医科大学伦理学委员会批准(2017)研究对象已签署知情同意书。

1.2设备和软件PhilipsMX8000型64排螺旋CT扫描机;HPZ600工作站:Inter(R)XeonCPU32G(十二核)48G内存。建模及计算分析软件:Mimics16.0(比利时Mater- alize公司)Geomagic(美国3DSystem公司)UG(Sie- mens PLM Software公司)Hypermesh(美国Altair公司)Abaqus(美国Simulia公司)。

1.3模型建立

用64排螺旋CT对志愿者的左肘关节屈曲90°进行CT扫描,层厚05mm,将得到的数据导入Mimics16.0进行三维重建建立肘关节三维数字模型。在 CT断层图像中,不同组织的CT值不同因此以600灰度值为界限严格区分皮质骨、松质骨。再根据尸体解剖样本肘关节屈曲90°状态描绘韧带骨间膜等不能显影的软组织,其他部分予以去除同时参照解剖图谱,运用Geomagic软件对模型进行处理,建立肘关节内、外侧副韧带桡骨头环状韧带及前臂骨间膜及关节软骨。最后运用UG软件对模型进行修饰、加工,剔除几何模型缺陷得到肘关节有限元模型(图1)。

1.4肘关节网格模型的划分及组织材料属性赋予

运用Hypermesh软件对肘关节模型进行高质量的网格划分,对于应力集中或相互接触的地方,增加网格密度,以提高接触区域的计算精度。网格划分后各模型参数:尺侧副韧带单元数28379节点数7627,桡侧副韧带单元数37526节点数9900环状韧带单元数18296,节点数5273骨间膜单元数94086节点数22863软骨单元数93753节点数32064骨皮质单元数207267节点数50610骨松质单元数129996.节点数32559地面单元数418节点数420。最后将肘关节模型中的各种组织材料全部简化为各向同性的均质性弹性材料,查阅文献[10]对肘关节各组织材料进行属性(表1)赋值。经测算中等身材成人的上肢质量约45kg而有限元建模时只能计算得到骨骼的质量,因此运用Abaqus的装配模块，计算得出上臂及前臂软组织质量并将其赋予相应部位。

1.5肘关节的装配及定位

运用Abaqus对装配好的肘关节模型进行旋转或平移等操作最终得到装配好的肘关节30°60°及90°模型(图2)。

1.6定义接触关系、边界条件及分析步长

采用通用接触算法肘关节各组织之间的接触滑动系数为0.01。按照肘关节运动特点设定边界条件，对地面进行固支约束。肘关节与地面的撞击时间定义为100ms.

1.7 计算

在Abaqus6.14中分别对肘关节模型施加一定的初始撞击速度,应用Abaqus/Explicit中“brittle cracking准则对骨皮质破裂行为进行判定,通过不断“试凑”“迭代”计算出肘关节撞击地面后骨皮质破裂的临界速度并得到骨皮质断裂裂口处的最大应力。在此速度基础上增加撞击速度可得到骨皮质显著破裂的速度。

2结果

21前臂与地面呈30°撞击地面时的计算结果

撞击前后3个时间点肘关节位移云图见图3。当前臂与地面呈30°撞击地面时骨皮质破裂的临界速度为11.4m/s。在该临界速度下尺骨骨皮质破裂;此时，骨折处最大应力为135.1MPa(图4A)。当肘关节初速度达15m/s,则尺骨呈显著破裂(图4B)。

2.2前臂与地面呈60°撞击地面时的计算结果

撞击前后3个时间点肘关节位移云图见图5。当前臂与地面呈60°撞击地面时,骨皮质破裂的临界速度为9.3m/s。在该临界速度下,尺骨骨皮质破裂;此时,骨折处最大应力为2121MPa(图6A)。当肘关节初速度达13m/s,则尺骨呈显著破裂(图6B)。

2.31前臂与地面呈90°撞击地面时的计算结果

撞击前后3个时间点肘关节位移云图见图7。当前臂与地面呈90°撞击地面时,骨皮质破裂的临界速度为13.8m/s。在该临界速度下,肱骨骨皮质破裂:此时,骨折处最大应力为1485MPa(图8A)。当肘关节初速度达18m/s时,肱骨呈显著破裂(图8B)。

2.4结果分析

肘关节与地面呈不同角度撞击时骨皮质破裂的临界速度及显著破裂时的速度见表2。当前臂与地面呈30°60°时在尺骨鹰嘴及肱骨滑车等部位出现应力集中现象(图46);当前臂与地面呈90°时在肱骨外上髁及肱骨内上髁等部位出现应力集中现象(图8);各个工况中桡骨头应力传导较小,且单元应力峰值远低于桡骨的屈服强度。

3 讨论

肘关节由肱骨下端和尺骨、桡骨上端构成,是由3个关节共同包裹在1个关节囊内组成的复关节。肘关节损伤受多结构的影响,外力作用瞬间关节内诸骨及关节辅助结构所受应力复杂多变,为肘关节的生物力学分析带来挑战。由于有限元显式动力学(Abaqus/ Explicit)能迅速地分析许多独立物体相互作用体系的复杂接触问题,因此特别适合于分析受冲击载荷和在复杂接触条件下结构内部相互作用的瞬间动力反应。相比隐式分析程序(Abaqus/Standard)中材料的失效和破坏常常导致严重的收敛困难,Abaqus/Explicit能很好地模拟材料的失效行为。近年来,已有研究表明，利用有限元法进行人体各组织生物力学分析科学、可行12]。因此,本课题采用Abaqus/Explicit模块对肘关节模型碰撞地面进行损伤分析,模拟肘关节碰撞地面时的情况以及生物力学性能。

曾智13]杨运平等[14]将肘关节伸直位有限元模型运用于生物力学有限元分析中,通过三维CT数据建立了成人肘关节及前臂三维有限元模型并对比分析桡骨头应力的传导规律与肘关节实体模型数据验证了模型的有效性:张建新等建立了肘关节伸直位三维有限元模型并对桡骨头切除术后并发症进行三维有限元分析:谭军16通过三维CT重建研究了活体肘关节侧副韧带,探讨肘关节屈曲过程中,肘关节侧副韧带的生物力学变化情况。人体肘关节功能位为屈曲90°左右,同时由于人体的自我保护意识.因抵抗或碰撞引起损伤时肘关节往往呈屈曲状态,因此.本研究在前期研究的基础上建立肘关节屈曲90°状态时有限元模型。HAMEL等7及杨资洋等18进行的有限元分析分别采用模型与撞击面呈45°90°及0°45°的夹角撞击,发现撞击角度是骨折发生及损伤严重程度的影响因素。因此,本研究通过模拟肘关节有限元模型分别与地面呈30°60°90°撞击地面能够更加全面准确地分析肘关节碰撞损伤时的生物力学特点。同时，本研究充分考虑了肘关节内、外侧副韧带桡骨头环状韧带,关节软骨及骨间膜等关节附属结构对肘关节碰撞损伤时的影响19-21结合影像学技术,运用 GeomagicUG等软件准确地在骨性有限元模型上构建了上述关节附属结构精确了肘关节力学模型,使其能更准确地用于分析肘关节碰撞生物力学研究。肘关节碰撞是一个动态的过程,在碰撞时肱骨、尺桡骨会以肘关节为中心点做伸屈及旋转等运动。本研究按照肘关节运动特点设定边界条件,综合考虑肘关节在实际碰撞过程中肱骨、尺桡骨的运动情况。通过对肘关节模型进行损伤分析及计算得到肘关节屈曲90°状态下分别以不同角度撞击地面骨皮质破裂时的临界速度以及显著破裂时的速度。在各个工况中尺骨鹰嘴、肱骨髁及肱骨滑车等部位出现了高应力的红色区域,表明这几个部位在正面撞击下易出现应力集中现象,也就是骨折多发部位。此外,本研究模拟结果发现桡骨头因解剖位置关系,在正面撞击中不易与撞击物直接发生碰撞,而是通过移位碰撞肱骨后产生应力传导,因此肘关节屈曲90°状态发生正面碰撞时应力传导主要集中在肱骨与尺骨,而桡骨头与肱骨和尺骨相比应力传导较小。这与MORREY等22发现“肘关节随屈曲角度增大,经桡骨传导的应力逐渐减小”的研究结果基本一致。

本研究通过建立人体肘关节三维有限元模型对肘关节撞击伤的损伤特点以及相关生物力学特点进行了初步的分析进一步明确了肘关节直接碰撞损伤时力的大小、速度及方向等，为生物力学应用于肘关节损伤研究奠定了基础。但由于人体组织结构的复杂性，加之个体特异性用模拟实验反应损伤过程仍有一定的局限性,临床仿真度及指导性有待进一步提高。本课题组将在后续开展的实验中采用尸体标本进行验证和对比，以期为生物力学在肘关节损伤的医学实践中应用提供依据。

来源：

唐鹏程 ，朱 英 ，赵敏珠 ，代亚磊 ，张 丽  ，李红卫 ，李剑波 ， 400016 重庆，重庆医科大学：基础医学 院法 医学教研 室 ，重庆市刑事侦查工程技术研究 中心 ；400021 重庆 ，重庆市公安局刑警总队技术处

文章仅作为学术交流，版权归原作者所有。如涉及版权，请联系删除！

PART  01

针对骨学、关节外科、普外科、囗腔科等提供医学力学有限元分析仿真、培训、临床手术模拟分析等；

代做有限元：脊椎、腰椎、颈椎、上下肢、假肢、前臂、血管支架、牙齿、骨关节等；

服务对象：各省市、自治区从事运动生物力学、生物医学工程、基础医学、临床医学、囗腔力学、骨骼力学的教学、研究与应用的教师、科研工作者、各级教练员等方面的相关人员；国内各重点大学、科研院所相关研究领域的博士、硕士相关研究生和学者等。

相关项目代做或课程培训可主页私信小编