股骨有限元模型构建及钝力性损伤的研究.docVIP

股骨有限元模型构建及钝力性损伤的研究 　　摘 要：对股骨进行多层螺旋CT扫描，利用Mimics软件将CT医学影像数据转化为有效的股骨有限元模型。将模型在Workbench显示动力学Ls-Dyna模块中进行网格划分、材料赋值及其它前处理，用碰撞物以不同的速度对其正面不同位置进行模拟打击，得到股骨被撞击结果，从而为不同股骨骨折受伤源头提供依据。实验结果表明，碰撞物的速度不同、被撞击的位置不同，骨折发生时间和位置也不一样。 　　关键词：股骨有限元；CT扫描；显示动力学 　　DOIDOI：10.11907/rjdk.171231 　　中图分类号：TP319 　　文献标识码：A 文章编号：1672-7800（2017）007-0137-03 　　0 引言 　　钝力性损伤主要包括运动的制伤物作用于相对静止的人体和运动的人体撞击静止的物体，同类型钝力性损伤程度和形态是推断暴力来源和致伤特征分析的主要依据。随着计算机技术的发展，有限元方法和CT扫描技术结合应用在生物工程学上，以解决一些人体损伤的问题[1]。利用软件建立和人体类似的三维模型，用钝性物体以不同的速度对人体不同位置进行撞击，分析应力、应变、位移、速度、加速度等应力响应参数变化，为钝力性损伤鉴定提供依据[2]。本文拟用钢管以不同的速度撞击股骨的不同位置，得到股骨在不同条件下各时间点的应力、应变和位移情况，为相关问题研究提供参考。 　　1 有限元模型创建及研究 　　1.1 有限元模型创建 　　利用CT机及AW图像后处理工作站对尸体股骨头上方至膝关节处进行扫描，将扫描得到的数据导入Mimics中，利用阈值分割工具将股骨和股骨中段的软组织提取出来，然后使用蒙罩编辑和区域增长等工具分离出股骨和软组织，之后用填充工具将空洞填充和修补，将提取出来的蒙罩三维重建，得到股骨及其中段软组织的三维模型。 　　1.2 研究方法 　　根据相关文献和实验得到的数据，将股骨简化为与应变率无关的塑性材料模型，密度ρ=1 990kg/m3，弹性模量E=14 635MPa，泊松比μ=0.343 5，屈服应力δs=133MPa，切线模量1 024MPa。失效准则采用最大主应变准则[3]，包括它的拉伸阈值（εTmax=0.073）和压缩阈值（εCmax=0.010 4）。中段软组织采用Mooney-Rivlin材料模型建立，材料相关属性[4]如下：密度ρ=3 600kg/m3，两个Mooney-Rivlin常数为C01=4.1、C10=0.41。碰撞的钢管材料采用常见的结构钢，密度ρ=7 850kg/m3，弹性模量E=200 000MPa，泊松比μ=0.3，体积模量K=166 670MPa，剪切模量G=76 923MPa。 　　因为是模拟股骨被钝物撞击，涉及到瞬态、大应变、大型变、材料的破坏、材料的完全失效或伴随复杂接触的结构问题，所以通过Ansys显示动力学求解。Ansys显示动力学分析模块包括Explicit Dynamics、ANSYS AUTODYN及Workbench 16.0的Ls-Dyna。 　　Ls-Dyna是世界上最著名的通用显式非线性有限元分析程序，能够模拟真实世界的各种复杂问题，特别适合求解各种二维、三维非线性结构碰撞、金属成型等非线性动力冲击问题，它以Lagrange算法为主，兼有ALE和Euler算法；以显式求解为主，兼有隐式求解功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算），是通用的结构分析非线性有限元程序[5]，操作和?Y果查看都十分快捷，本文采用Ls-Dyna模块进行股骨模拟碰撞[6]。 　　为了模拟人体股骨受打击情况，给股骨模型上端施加一个垂直的力，考虑人体情况，拟添加一个250N的压力，股骨的下端固定，用直径为40mm、壁厚为5mm、长度为1.5m的结构钢钢管正面打击，分别以20m/s、10m/s的速度去撞击股骨模型中段，然后以10m/s的速度撞击股骨模型靠近膝关节的位置，用Workbench计算股骨和钢管的应力变化云图，得到股骨的von mises应力变化表和股骨在不同时间点的应变、形变情况[7]，得到股骨中段被钢管以不同速度撞击和股骨中段、下段被相同速度撞击的生物力学响应，从而判断股骨骨折情况，为司法鉴定中已知股骨的受伤情况来反推受伤源头提供一定的科学依据。 　　2 分析过程及结果 　　2.1 添加材料 　　使用Worbench16.0进行分析。打开Workbench16.0，启动Ls-Dyna模块。在Engineering Data中选择新建一种材料，这种材料用于股骨的各项弹塑性材料模型。在Tool Box中依次选择密度、各项同性弹塑性材料、双线性各项同性硬化和主要的应变失效，密度中填入1