有限元分析在步态生物力学研究中应用进展.docVIP

有限元分析在步态生物力学研究中应用进展 　　摘要：有限元法应用于步态生物力学研究后不断获得发展，已成为当今生物力学研究的重要方法。综述了该方法在计算机建模、生物材料学、临床问题和体育医疗器械开发等领域的应用进展，并根据其发展脉络作出了趋势展望。 　　关键词：生物力学；有限元分析；步态研究 　　中图分类号：G804.01文献标识码：B文章编号：1007-3612(2007)08-1080-03 　　 　　有限元法(finite element method,FEM)是一种在工程科学技术中广泛应用的数学物理方法，用于模拟并解决各种工程力学、热学、电磁学等物理场问题。有限元法应用于生命科学的定量研究,已取得了长足进展。在步态生物力学研究领域，尽管足底压力测量、影像分析和表面肌电等技术被广泛应用，但足部骨骼、软组织等结构内部的应力传递机理始终无法得到力学解释。当FEM应用于足部生物力学研究时，复杂的骨胳几何结构、边界条件和材料的不均匀性等问题便找到了可能的解决途径。 　　 　　1步态的数值建模几何参数 　　 　　1.1X光照片从不同方位的多幅?X?光照片获得几何数据重建足部的三维模型是一种方便、经济的方式。Jocob 和Patil［1,2］通过某健康者和某患有Hansen氏病的主体的X光照片获得的足部几何数据，进行了有限元对比分析。但是该方法所得到的模型相对粗糙, 只能给出足部的大致轮廓, 难以细致分辨足部各骨及其组织。 　　1.2 CT扫描图像的三维重建CT主要为横断层成像，但带有越来越多的图像后处理系统，如多维断层重建、三维图像重建、扫描后再次放大、立体模型与几何模型测量法等，可以清晰地反映出物体内部的不同结构和组织，对研究及建模带来很大方便。Camacho等［3］分析了间隔为1mm 的人体足部CT断层扫描照片(图1)，为进行有限元分析建立了一套较精确的模型。 　　1.3MRI图像重建MRI为核磁共振成像，该技术既具有成像清晰、解析度高的特点，又能清晰显示人体结构的组织学差图1人体足部CT断层扫描建立模型异和生化变化，在不改变体位的情况下，可直接获取横、矢、冠、斜4种断层图象。研究人员可以得到十分细致的足部几何模型。Gefen 等［4］从MRI图像中重建出正常人足部骨??的几何构形。Zhang Ming等［5］通过Coronal MRI技术，间隔2 mm建立了糖尿病足的模型(图2)。 　　图2糖尿病足的模型通过上述医学成像系统获得的足部骨骼的几何学数据是比较精确的，但对于韧带、肌肉和足底腱膜等其它组织，则还需要从解剖学作更深入的研究，利用图像后处理分析（由输出的通用图像数据文件，存储在计算机中；采用直方图修正法来改善图像质量，或高斯滤波进行图像平滑以减少噪声，Sobel梯度算子检测边缘）来获得可靠的模型几何数据。常用的影像分析软件有mimics、3d-doctor等。通常方法是从上述软件得到各组织的边界数据，再从制图软件（solidwork、pro/e、ug等）建立实体模型，这样才有利于有限元分析软件（ansys、abaqus、femlab等）更好的处理。 　　 　　2步态的数值建模材料力学参数 　　 　　Huiskes［6］研究发现在准静态荷载下，骨和软骨可视为各向同性的线弹性材料。但是，Clift［7］对关节软骨做了有限元分析，证实关节软骨的应力―应变关系随时间的变化呈现出一定的非线性，故把关节软骨看作一种非均质、各向异性的材料更合理。Keyak［8］等把骨的应变张量看成力学激励，这样，应变张量中的各元素就是骨再造方程中的系数，然而，系数之多又给计算机模拟数值解造成了一定困难。Koosltra［9］等假设骨能够根据其应力―应变状态来调整小梁骨的方向性和内部密度分布，用应变能作为力学激励，由于它是标量，使方程的系数为常数，可以保证数值解的取得。因此，骨单元的划分一般采用六面体块单元或楔形单元。考虑到关节软骨在足部运动中起的重要作用，如果把关节软骨都简化为线弹性材料，计算结果必将有很大误差。因此，从模型的整体出发，在建模过程中一方面要注意把一些对足部运动影响不大的关节软骨(如跗骨间关节)加以简化，以减小模型的复杂度；另一方面，为保证计算的精确性，对于一些大关节软骨(如 距下关节)则要考虑为非线性，最好结合影像分析，获取其运动参数，以提供较为准确的边界条件。如何定义足部各关节的运动方式，并建立相应的足部模型约束条件是足部有限元分析的难点。 　　韧带，足底皮肤及肌肉等软组织属于非线性材料，所以可划分为簧单元。Race and Amis［10］对正常人体下肢韧带进行了单轴拉伸实验，得出韧带的材料参数可按下式取值： 　　 　　3步态有限元分析的应用 　　 　　3.1医学的临床诊断应用Jocob