有限元分析在股骨股骨生物力学研究中的应用.docVIP

有限元法(finite elementmethod,FEM)是一种工程力学研究方法,它是用较简单的问题代替复杂问题后再求解，计算过程中将求解域看成是由许多称为有限元的小的互连子域组成，有限元是那些集合在一起能够表示实际连接域的离散单元，对每一单元假定一个合适的(较简单的)近似解，然后退到求解这个域总的满足条件，如结构的平衡的条件，从而得到问题的解。有限元法正逐渐被应用到矫形外科领域, 极大的促进了生物力学的发展。早期FEM主要应用于简单的应力分析。随着三维非线性有限元软件的开发, FEM的应用取得长足进展,目前主要从微观水平和定量测定两个方面研究骨各种时间依赖的生物力学特性。股骨联接髋关节和膝关节，是人体最大的骨头，是人体正常运动不可或缺的支撑和活动单元。作者主要从股骨力学性能研究，应力分析，预测骨折发生的危险度及其好发部位，指导骨折固定物的设计和改进及手术评价与术式选择方面综述了近年来有限元法在股骨的一些研究进展。 股骨力学性能研究。 股骨力学性能研究已有百余年的历史.自从1892年Wolff提出Wolff定律以来，人们对股骨的生物与力学特征进行了广泛研究.例如，Koch证实骨小梁按主应力迹线排布，当荷载变化时，骨小梁将按新的应力迹线重新取向[[3,4]. Carte:等[[5]研究提出股骨小梁体系受力并不遵从Wb班的主应力迹线定律.文献[6,7】详细地研究了骨小梁结构与力学性能的关系.Mullende:等[8]提出骨结构优化的自组织控制模型，并结合有限元方法模拟了股骨结构的受力过程.Keyak等[9]通过试验和有限元方法模拟研究了静载和冲击荷载下股骨的破坏行为.Cody等[lo]提出应用有限元方法预测股骨颈的断裂强度比应用QCT和DXA方法更为有效.Keyak等[ll]利用CT扫描技术构建了股骨结构有限元模型，并利用该模型定量分析了股骨的破坏荷载.Bayraktar等[12]系统地总结了股骨组织力学性质的实验数据和测试方法，通过试验详细地测量骨皮质的弹性模量、屈服应变，为准确预测骨组织的破坏性质提供了参考.王远亮等[4]系统地总结和分析了骨组织生物力学的研究.王以进等[13]研究了扭矩对股骨的破坏作用.吴维才等[14]通过试验研究了扭转荷载下股骨的力学性质.赵宝林等[15]通过试验研究了股骨下端松质骨的拉伸、扭转、剪切、弯曲等力学性能.李德源等[16]利用数值模拟方法分析了撞击荷载作用下具有横观各向同性松质骨的动力响应特征. 2股骨应力分析 应力分析是FEM的主要用途之一。松质骨应力特性依赖于骨密度和骨小梁结构的各向异性,传统的实验方法根本无法全面、精确描述其应力特性。三维有限元分析( finite element analysis, FEA)可以很好地表现骨关节等组织应力分布的总体趋势[1]。李苏皖等【1】应用有限元法，通过建立股骨近端模型，分析人正常站立股骨的应力分析取得很好的结果。刘杨等【33】建立了完整的正常人股骨三维有限元模型，并应用于正常站立位生物力学分析，验证了股骨骨折多发生于股骨体中下段的理论，并证明了老年人多发股骨颈骨折，并且骨折难以愈合的事实，从生物力学角度说明股骨颈骨折后应避免活动。Gallagher等(1982年采用FEM对股骨的应力作了详尽的研究,他通过对比骨小梁分布和应力分布状态,发现骨小梁主空间构架分布与主应力方向一致。在单足站立时,股骨各截面上以压缩应力为主,且分布恒定,而非弯曲变形;同时可观察到股骨的压缩应力分散,并发现这种应力分散与股骨的截面几何形状有关[2]。Key-ak等[3]分析了引起股骨骨折的各种应力,发现扭转和剪切应力在股骨骨折发生过程中起重要作用。 3. 预测股骨骨折发生的危险度及其好发部位 临床用DEXA或QCT等骨量/骨密度测定的方法来预测某一部位的骨折危险性,有一定的局限性,因为骨质疏松性骨折的发生与骨量、骨强度和骨的结构均有关。对骨质疏松患者骨折危险性的预测较好和较准确的指标是骨的力学特性。Cody等[7]通过比较QCT、DEXA和三维FEA作为股骨颈骨折的预测指标后认为,三维FEA是预测骨折载荷的最好方法,并能用于在体骨强度变化的监测。二维FEA基于一种简单、非侵入性、低花费的图像技术来获得与骨强度特性有关的资料,并且计算量很少,可用于大样本骨质疏松人群中股骨颈骨折的预测研究[8]。而仅依靠数字X线片通过代数重建技术(algebraic reconstruction technique,ART)来建立三维有限元模型,有利于提高FEA在预测骨折发生风险方面的作用[9]。鞠杨等【34】建立了一个针对高危女性人群股骨结构的三维有限元模型，分析了运动荷载作用下此类股骨的应力和应变分布以及股骨破坏的可能部位.结果表明:运动荷载使股骨产生较大的横向位移和弯曲变形;股骨颈处的应力