基于abaqus的血管支架有限元分析

寂静回声

血管支架是一种薄壁管状结构，利用其可扩张的特性，借以支撑血管狭窄部位来治疗冠状动脉。　　本文利用solidworks绘制了血管支架及相关模型，并利用abaqus分析了某型号血管支架植入人体过程中的静态受力分析。对血管支架的研究设计以及血管支架的医学植入具有一定的理论指导意义。血管支架治疗的过程为：在植入前，将血管支架放在气囊外侧，并将其运输到动脉血管狭窄部位处。通过给气囊充气，致使血管支架支撑狭窄部位的血管壁。为了全面仿真某型号血管支架的植入过程，本文建立了气囊、血管支架、血小板、血管四个物理模型。
　　1.1 气囊
　　气囊的图形如下图所示。两边为半球，中间为圆柱。气囊的初始直径为1.4mm，厚度为0.02mm，长度为10mm。由于气囊很薄，模型中使用薄膜单元。通过SolidWorks的abaqus插件将实体无缝连接到abaqus。

血管支架由主筋和连接筋组成。主筋的草绘图形如下图2所示，连接筋的示意图如下图所示。




2.1 气囊
　　气囊的单元采用显示动力学缩减积分四节点薄膜单元M3D4R。网格划分中单元的尺寸为0.1mm，气囊的材料模型选用弹性模型。
　　2.2 血管支架
　　血管支架的单元采用显示动力学缩减积分八节点单元C3D8R。为了使血管支架可以进行扫描划分方法，需要对从creo导入进来的模型进行拓扑合并以及相关的切割操作。网格划分的单元尺寸为0.02mm。主筋的材料模型选用Ti-Ni形状记忆合金。相关材料的物理性能参数采用东北大学裴丽丽硕士论文《生物医用TiNi形状记忆合金的制备及性能研究》中相关性能参数。
　　2.3 血小板
　　血小板的单元采用显示动力学缩减积分八节点单元C3D8R。网格划分的单元尺寸为0.4mm。血小板的材料模型选用6项式超弹性模型，其中材料模型中D取0，材料具有不可压缩性。
　　2.4 血管
　　血管的单元采用显示动力学缩减积分八节点单元C3D8R。网格划分的单元尺寸为0.5mm。血小板的材料模型选用Neo-Hookean超弹性模型。整个模型划分网格的节点总数量为213067，单元的数量为146814。

3.1 载荷
　　气囊进入血管壁紧缩部位，对气囊进行充气，到一定压力后保持，进而泄气。因此，本文采用如下的加载方式：在0-0.03s之间，给气囊内表面施加3.6Mp的压力；在0.03~0.05s之间保持压力3.6Mpa；从0.05~0.06s之间，压力逐渐将为0。选用abaqus中的光滑幅值加载曲线。
　　3.2 边界条件
　　固定血管的两端，使其U1，U2，U3等于0。同时血管由于模型的对称，在血管中间部位只存在径向位移，因此需要约束血管中间面的轴向位移，为了加载这一边界条件，在模型的建立中，需要对模型进行分割处理。
　　3.3 接触的设置
　　在载荷分析的初始，设置主筋和连接筋为Tie连接，同时血小板和血管也为Tie连接。在第一步分析时，设置气囊和血管支架、血管支架和血小板、气囊和血小板为surface-to-surface 连接。算法采用罚刚度算法，滑移选用有限滑移。接触类型中，正向压力采用硬接触，切向压力中，摩擦系数选为0.2。同时，气囊的外表面设置为自接触。
　　3.4 模型的求解
　　由于本文模型运用的是abaqus的显示动力学求解准静态模型，因此，需要采用一定的方法来加速问题的模拟。在abaqus中准静态加速分析的方法包含质量放大，加大加载速率等方式。
　　Abaqus explicit中稳定的时间增量与单元的特征长度成正比，与材料的膨胀速度成反比；而材料的膨胀速度与弹性模量的根方成正比，与材料的密度根方成反比。增加abaqus explicit 求解时的稳定增加时间增量，以加快求解速度。因此，在模型的求解中通过提供材料的密度进而对物体的质量进行放大，以实现求解速度的加快，经过多次尝试之后，本文采用的质量放大系数为1000。