一种股骨材料仿真软件的开发及可行性分析-论文.pdfVIP

中国科技论文在线 一种股骨材料仿真软件的开发及可行性分析 重庆大学机械传动国家重点实验室，机械工程学院，重庆（400044 ） E-mail ：meadowba@ 摘 要：针对多体几何关系模型中骨材料赋值困难的问题，通过两种方法建立股骨材料仿真 的有限元结构，并对比其生物力学特性来探讨自主提出的骨材料仿真方法的可行性。先由股 骨CT 图像数据建立其三维实体模型，利用CT 图像灰度值与骨材料弹性模量间的关系，借 助ANSYS 软件参数设计语言APDL 与VC++平台结合开发骨材料属性赋值软件，从而建立 起基于几何关系的股骨各向异性非均匀有限元结构A；再由CT 图像数据直接利用商业专用 软件 Simpleware 得到基于象素的近似真实股骨材料特性的有限元参考结构 B 。通过计算分 析得到，结构A 与参考结构B 的应力应变分布相近，且比较点的相对误差在5%左右，仿真 结果表明这种自主开发的骨材料属性赋值软件能够较好的实现骨组织材料仿真。 关键词：股骨；材料仿真；各向异性；有限元 中图分类号：R68; TP391 1 引言 为解决植入假体存在的松动和下沉等并发症，在设计定制型人工髋关节股骨内置假体过 程中，需要研究定制假体对股骨力学性能的影响，从而优化设计人工股骨假体。探讨假体对 股骨的生物力学特性，首先要研究股骨的非均质材料特性。目前，在骨材质研究领域，很多 学者把股骨仿真模型的材料被看成是各向同性的，这就不能真实反映骨组织微构造的本构关 系[1] ；还有些学者[2-3]基于 CT 图像结合不同软件建立了股骨有限元结构，但都只是单纯的股 骨有限元结构，且对骨材料赋值的研究也不够深入。本文将深入研究多体几何关系模型中骨 材料赋值问题。利用 CT 数据，有两种建立股骨三维有限元结构的方法[4]，一种是基于几何 关系的，另一种是基于象素的。基于象素的有限元结构可以获得较好材料特性，但给植入假 体对股骨的响应特性及其接触分析研究带来困难；利用基于几何关系的模型可以有效地利用 强大的辅助软件 CAD、CAE 来深入研究骨—假体之间的响应特性，从而优化设计假体，但 其不足之处就是较难模拟多体结构中骨的真实材质。为此，本文就是要利用基于几何关系模 型能便于分析多体相互作用的优点，解决其骨材质真实性模拟的困难。建立股骨精确的几何 模型，尽可能模拟骨材料特性的真实值，尤其是基于几何关系建立正交各向异性的骨有限元 [1，5] 结构，是当前该领域的一个重要研究方向，也是难点课题 。 本文首先从 CT 数据出发，根据灰度值与骨材料的转换关系，开发材料赋值软件，实现 股骨材料各向异性非均匀模拟，建立起基于几何关系的股骨材料仿真有限元结构。其次借助 于商业专用软件 Simpleware 建立基于象素的股骨有限元结构。然后对上述两种不同方法建 立起的有限元结构在步态单肢负重状态下进行计算其应力、应变，通过分析比较，来检验自 主开发的材料赋值软件的可行性，进而利用该材料赋值软件来帮助深入研究骨组织材料的生 物力学特性。 2 材料与方法 研究股骨标本取自成年健康男性尸体的左侧股骨。利用第三军医大学西南医院德国产西 门子SOMATOM SENSATION16 排螺旋 CT 机扫描整个股骨，扫描技术参数为：间距为 本课题得到国家高等学校博士学科点专项科研基金（20060611014 ）；重庆大学研究生科技创新基金 （200801A1A0090268 ）资助。 -1- 中国科技论文在线 1mm，层厚为 1mm，窗宽 500，窗位 40 ，一共得到 443 张 CT 图像。编写计算程序来精细 地仿真股骨表观密度的区域性分布，通过灰度值与弹性模量之间的转换，完成将股骨 CT 灰 度图像数据向骨的材料属性的转换，对各单元的各材料主方向上弹性模量具体化，实现股骨 材料的非均匀及各向异性描述，建立起更接近真实股骨的计算力学模型。 2.1 CT 值与灰度的变换 按 CT 图片格式，