Ansys应用：临床骨科中应用的有限元分析法-新理论与新进展

2017-03-23 by:CAE仿真在线来源:互联网

目前,有限元分析在骨生物力学研究、内固定材料研究、骨科疾病病理机制研究等方面的应用已越来越广泛。

几何模型的建立

将CT、核磁共振成像等二维的断层影像数据,导入Mimics中,通过阈值分割、区域增长、空腔填充和三维建模等步骤获取人体几何模型。

再通过geomagic等软件对模型进行修改,使分析模型尽可能的是由光滑面或体构成。

划分有限元网格

有限元网格划分就是将整体划分为个体进行分析,最后将个体集合为整体进行研究。其中有限元网格划分的质量直接影响分析结果的精度, 主要表现在以下3方面:①几何模型结构细节处理。 ②网格划分密度。③网格过渡。

一般会把重建三维分析模型划分成有限个四面体或六面体,这样就会造成分离出来的正面体 (正四面体或正六面体)变形,其后果即是越发加强了对有限元网格划分这一操作的质量要求。

骨骼材料属性定义

现在对骨组织的属性赋值主要是依据其CT值计算出骨密度值,再对骨组织进行划分有限元网格单元,根据每个网格单位CT值来对骨组织的属性赋值,从而描述骨的生物力学特性。

确定边界条件

人体骨骼的承载和运动是和其周围的肌肉等相关组织分不开的,由于这些结构之间的配合比较复杂,要建立一个完整的有限元分析模型是比较困难的。所以,为了便于分析求解往往会简化分析模型。待简化之后明确正确的边界条件也是有限元分析的重要环节之一。

分析与结果验证

分析作为有限元分析仿真的执行环节,主要工作是确定一种适合分析模型本身的有限元算法,为了证实分析所得到的结果是否正确,需要利用所掌握的相关知识结合临床综合分析,或者进行可行性实验进行反复验证并分析,得到相对客观数据,从而归纳并分析得出其一般规律,提高所得数据的可信度。

有限元分析在脊柱骨折研究领域的应用

利用有限元分析法可以对脊柱骨折的生物力学功能和致病机制等进行研究及评价。

Silve等通过有限元法证实进行椎弓根螺钉置入时,皮质骨提供的强度占椎体强度的10%,在相对较高的压力载荷下,终板中下方的松质骨发生骨折时,椎体的强度明显下降,而骨内压则随之升高。

Cao等建立了三椎体腰椎有限元分析模型,分析结果标明34%的椎体压缩骨折发生在上终板或下终板附近,而63%的椎体压缩骨折发生在上下终板之间。结果认为骨折的发生与负荷类型没有关系, 而与皮质骨的弹性模量成正比例关系,与松质骨的弹性模量成反比例关系。

Dall Ara等通过非线性解剖特异的定量CT数据位依据建立的有限元分析模型,预测了体内椎体的力量。随后应用离体实验验证了该模型的正确性。

Matsuura等利用标本研究得出有限元分析可能是预测椎体成形术后带来的原位椎体骨折应力负荷的有用工具。

有限元分析在关节置换中的应用

有限元分析技术为了解关节置换患者术前关节功能状态、术中假体的选择、术后预测提供了一种全新的方法。

Chang等通过对正常股骨头外形进行有限元分析后,建议将假体股骨柄的股骨颈设计成相对较细,距股骨头顶端相对较短的形态,其理论依据是这样可以适当减少对周围骨质的破坏,从而减少骨结构重塑,同时假体稳定性可以保障。

Patil等通过有限元法分析全髋关节置换后髋关节的应力分布表明,当髋臼杯外展角由35° 增加到55°时,髋关节内应力会逐步上升,其磨损程度也会增大;而随着髋臼杯的前倾角的增大,髋关节内应力会逐渐减小。从而论证得出髋臼杯置入位置和其磨损存在必然的联系。

Sakagoshi等利用有限元分析法技术来研究股骨头坏死面积、术中假体安放位置与髋关节置换后生物力学的变化三者之间的关系, 其实验出现随着股骨头坏死面积的增大,骨-骨水泥界面的应力会逐渐增大;当假体内倾角处于内翻位时,骨水泥的最大应力增加。内倾角110°逐渐增大至 140°时,假体的峰值应力逐步减小,从而论证出股骨头坏死面积大节假体的三维有限元模型,在股骨头中心加大关节面负荷量,采用非线性接触弹塑性分析聚乙烯内衬关节面的接触应力。从而证实表明,选择股骨头直径相对较大的髋臼假体,并严格控制假体置入在安全范围内,术后减少上小与假体位置不当均是影响髋关节置换后生物力学环境的重要因素。

有限元分析在四肢与骨盆骨折中的应用

有限元分析法可以在四肢骨折固定方式选择以及术后对固定效果的评价提供相应依据。

Bresina等在1995年用有限元分析评估点状内固定(PC-Fix)系统的内固定效果;

Rogge等通过对桡骨远端骨折的三维有限元模型采用不同内固定方式所产生的不同生物力学特征进行分析,从而以寻找出最优固定方式。

Qian等在尸体实验中研究发现厚度、极限载荷和极限应变在更大的近端部分伸肌,机械性能的变化可能会影响预测力的传输,术前建模时应考虑伸肌装置的作用予以建模。既往骨盆骨折的治疗是基于广泛的临床经验。

Fan 等在有限元模型的基础上,对3种固定系统的生物力学通过有效强度水平、压力分布、力转移,沿骨折线位移进行评估,发现所有3个固定系统可以用来获得有效的功能结果。前柱钢板结合四边形螺钉是最优方法丁字形的髋臼的骨折。这个固定系统可以减少很多与其他固定系统有关风险和限制。

有限元分析在病理机制方面的研究应用

Chen等在 2003年,通过有限元分析对于糖尿病患者足底应力分布改变进行了研究,这对于矫形支具的设计和材料选择都有很大的帮助。

Goto等运用压缩屈曲理论, 通过有限元仿真分析模拟出脊柱侧凸形成的过程,研究其病理机制。

Campoli等研究股骨的形状和密度分布及其振动固有频率之间的关系,得出形状和密度可能导致大型计算固有频率的变化,不同的有限元模型需要准确估计骨的自然频率。

有限元分析在骨应力分析中的应用

采用有限元分析技术对正常步态下足踝部的解剖结构和力学接触特点进行研究,从而获取足踝部在各种受力情况下各关节面之间接触应力情况、应力分布部位以及应力值。也可对踝关节骨折术后导致创伤性关节炎、退行性病变等并发症采用建立有限元分析模型来进行力学分析。

黄诸侯等对跟骨的应力分析表明,在踝关节中立位,足背伸20°时,跟骨应力最大点为跟距外侧的关节面与由内后斜向前外方向的相交处,除此之外后距关节面和跟腱之间的部位承受应力也相对较大。

有限元分析在骨科中的发展方向

过去对骨的力学分析中,由于人体骨骼周围环境的复杂,其周围的肌肉、骨连接等,对骨骼受力均有影响,对其进行数据分析时会出现数据量过大处理困难的问题。故通常在分析时都会采用简化模型的方法来处理, 或者把一些影响因素转换成载荷或约束来进行分析,这样必然会影响结果的精准性。

随着计算机技术的飞速发展,并行运算技术的广泛应用,和GUP 图形加速卡的出现使得个人超级计算机成为可能,有限元分析运用的领域也越来越广泛,它不再局限于对单个实体进行处理,而是可以处理包括多部件类似装配体形式的分析模型。

Pang等对跟骨接骨板、螺钉和骨等形成的装配体进行研究发现在相同的载荷作用下的模型。跟骨的最大主应力集中在皮质骨支撑斜面,且跟骨的最大主应力和内固定系统与支撑螺钉固定模型中是小于无需支撑螺钉固定模型的,支撑螺钉固定模型中的应力更分散。这些为处理复杂非线性接触问题创造了条件,同时为研制和选择符合个体化生物力学的内(外)固定材料,并为一些退行性疾病的发生发展机制的研究提供更为精准参考。

Rotariu等在对复杂下肢截肢的患者步态周期进行分析研究得出有限元分析有助于选择合适假肢并为设计出合适的假肢也提供了良好的条件。

Choi等应用有限元分析研究纳米医学中单个细胞的力学。

有限元分析在骨科研究领域的优势及劣势

精确的诊断和正规的治疗计划对获得良好功能及较好远期疗效是很重要的,有限元分析通过计算机可自动提供三维图像,同时可给予不同的实验条件,进行反复模拟实验分析,快速获取完整多样的结果数据, 从而在术前发现潜在的问题,研究模型表面和内部区域应力-应变状态的变化。这些均可以为选择合适的植入物,减少偏差,从而为术中因植入物选择等问题节省部分手术时间,同时在手术结束后可以对术后效果进行有效评价。

有限元分析大多数为静态分析其应力-应变分布的状态,有限元方法在运动中的生物力学分析稍显不足,这使仿真分析的准确度和效率降低。在分析仿真的过程中,其结果本身的精确度受到很多因数的影响。例如,有限单元的划分对仿真计算结果的精度影响就很大。合理的单元形状和尺寸能极大提高分析结果的精度和可行度,反之则影响精度,甚至造成分析仿真过程失败。其次,因为骨折的骨折线是随机变化,因此不可能找到一个通用骨折固定系统,而有限元分析是针对个体案例的生物力学分析,故在应用上不易普及。

结论

有限元分析是一门成熟的学科,它在骨科中的应用从过去的静态分析发展到现在的动态仿真;

从单个简单模型发展到现在的多模型、复杂模型的接触分析;

从线性分析发展到现在的非线性接触或材料分析;

从单物理场发展到多物理场耦合分析。

开放分享：优质有限元技术文章,助你自学成才