人体脊椎腰骶段的三维有限元仿真模拟0724.docVIP

( 概述 有限元方法是根据变分法原理求解数学上可描述的物理问题的一种数值计算方法，其基本是用简单问题代替复杂问题再求解它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的近似解，然后推导求解这个域总的满足条件，从而得到问题的解 由于人体脊椎在结构形状、材料特性以及承载能力等方面都比较复杂，以往的研究手段都难以获得全域性信息，而有限元方法可用任意形状的网格来分割区域，布置节点，尤其是对脊椎这样复杂的结构有较好的适应性。本文将在Simpleware软件的辅助下建立人体脊椎腰骶段L3-L4段的三维有限元网格模型，然后给出实际情况中有限元应力应变的计算方法，最后利用ANSYS软件对得到的有限元网格模型进行受力仿真模拟，并比较验证了腰椎间盘突出情况下脊椎的应力应变分布。 有限元网格模型的建立 本文采用的数据是65帧连续的DICOM格式的人体脊椎矢状位图像，其中图像序列间隔为1.25mm，像素间隔也是1.25mm，分辨率为256x256，如图1所示。 首先将DICOM图像导入到专业的医学图像处理 软件Simpleware中，然后对二维DICOM数据进行手动 分割，得到人体脊椎腰骶L3-L4段椎体外部、椎体内部、椎间盘、间盘内部的二维轮廓，通过三维重建得到与真实椎体高度相似的三维模型，最后对模型进行网格划分，进而得到人体脊椎腰骶段L3-L4段的三维有限元体网格模型，为下一步的有限元力学分析提供仿真模型，图2和图3分别是三维有限元模型的单元和节点表示。 图1 人体脊椎矢状位图像 图2 单元表示 图3 节点表示 在有限元网格模型输出时，需要赋予各结构材料不同的科学参数：弹性模量（杨氏模量）和泊松比，从而将椎骨各部分简化为各项同性的弹性材料。其中椎体外表面单元用皮质骨描述，内部单元用松质骨描述，椎间盘纤维环基质也可以简化为各向同性弹性材料，现将其设置为可伸缩的弹簧单元，而将髓核假定为一种不可压缩的固体材料。 图3中浅蓝色节点表示的是脊椎的L3段皮质骨，红色节点表示的是脊椎的L4段皮质骨，深蓝色节点表示的是脊椎L3、L4段内的松质骨，黄色节点表示的是两段椎骨间的纤维环，紫色节点则表示的是纤维环内的髓核。 表1给出了输出的有限元网格模型各结构材料所包含的节点、单元数量和材料参数值。 表1 各结构单元数量和材料参数 组成部分 节点数 六面体数 四面体数 单元总数 弹性模量 泊松比 皮质骨L3 20901 6847 52991 59838 12000 0.3 皮质骨L4 22625 7762 55949 63711 12000 0.3 松质骨 12553 6418 22785 29203 100 0.2 纤维环 6616 2658 13385 16043 92 0.45 髓核 1835 915 2442 3357 1.0 0.499 有限元的应力应变计算 将得到的网格模型导入到专业有限元分析软件ANSYS中，为了简化计算本文仅采用SOLID45这个八节点六面体模型来模拟所有的四面体和六面体计算单元。 3.1 八节点六面体单元的应力应变计算 在实际的计算中，需要将八节点六面体单元映射为正方体单元，二者之间的坐标变换关系可以用局部坐标表示的形函数来表示：，其中 ， 表示实际单元的节点坐标。用u，v，w分别表示节点沿x轴，y轴和z轴方向的位移变化，若令 ，，则八节点六面体单元的应变位移为：。 将节点的应力状态用六个分量表示：， 其中是正应力，是剪应力；节点的应变状态也用六个分量表示：，其中 是正应变，是剪应变，则节点应力可以通过由弹性模量E和泊松比构成的三维弹性矩阵[D]和节点应变的乘积得到：，其中 节点的应变和位移的关系为： ，，，，，，则有，其中 根据虚位移原理就可以计算出三维八节点等参单元的刚度矩阵，即：，其中，表示单元的体积。 对于负荷矩阵的计算，因为所采用的是集中力载荷，所以只需要把该负荷放在某个适当的节点，并使它的负荷分量与节点的坐标方向保持某种一致，就可以获