1. 引言

血管疾病是如今社会最常见的一种疾病，该疾病的防治工作在取得初步成效的同时，又面临新的严峻挑战。作为心血管介入手术中常用的医疗器械之一，血管医疗支架 [1] 通常被定义为一种应用于植入型外科手术的管状金属结构，该结构的主要功能是对阻塞的冠状动脉或病变的血管提供机械支持，恢复血液在管道内的正常流通。

近几年对编织结构支架的研究逐渐兴起，相较于传统连接筋式的医疗支架 [2] [3] [4]，因其结构的稳定性和顺应性，具有很强的研究意义。计算机仿真技术 [5] 作为一种发展迅速的研究方法，能够有效模拟支架置于病变的血管组织并与血管内壁相互作用的过程。如今国内外学者对编织支架的有限元仿真研究较少，主要原因在于支架的编织特性 [6] [7] [8]，采用传统的建模方式很难实现。另一方面模拟支架受到血管径向压缩的过程中，在编织交错位置处存在大量的非线性接触，这些接触导致ANSYS在运行仿真时计算量繁琐且不易得到收敛。

综上所述，本文的研究数据通过可视化软件Constent实现不同几何结构参数支架模型的输出。该输出文件包含了支架编织路径中心线上所有点的柱坐标，基于Python语言编写脚本在建模软件AnsysSpaceclaim中读取坐标值并自动生成编织支架的三维模型。对比解析模型，主要从两个方面对编织支架的有限元模型进行验证分析：几何结构方面，提取支架与血管内壁作用过程中直径与整体伸缩长度的变化关系。力学性能方面，分析支架在输送完成后能否为血管内壁提供足够的径向支持，以及预测支架在受到血管壁作用时应力集中的部位。

2. 解析模型

解析模型 [9] [10] 由Jeweb和Clerk提出并且经过了实验验证，描述了支架在受到轴向拉伸或径向压缩时几何结构和力学性能的变化。该模型假设编织支架由多个独立的螺旋编织线构成，并且这些螺旋线只发生弹性变形。编织线丝的直径为d，不受外力作用时，支架横截面可简化如图1所示：

支架的平均直径D₀和外侧直径D_e0的关系如下：

$D_0=D_{e0}-2d$ (1.1)

如图2所示支架的表面沿着编织方向切割可转化为平面四边形，编织参数主要包括编织角度 $\alpha$ 、长度L和平均直径D。当支架受到轴向拉伸或径向压缩时，支架长度从初始值L₀伸长为L，直径从初始值D₀缩小为D，编织角度初始值从 ${\alpha }_{\text{0}}$ 减小为 $\alpha$。

支架长度L与直径D的变化关系如下所示：

$L=\frac{L_0}{\cos {\alpha }_0}\sqrt{1-{\left(\frac{D\sin {\alpha }_0}{D_0}\right)}^2}$ (1.2)

等效应力 ${\sigma }_e$ 与支架长度L的关系可由式(1.3)代入式(1.4)中得到：

$\alpha =\arccos \left(\frac{L\cos {\alpha }_0}{L_0}\right)$ (1.3)

${\sigma }_e=\frac{16\sin {\alpha }_0}{\pi d^3{D}_0}\sqrt{{\left(4EI\left(\sin \alpha -\sin {\alpha }_0\right)\right)}^2+3{\left(G{I}_p\left(\cos \alpha -\cos {\alpha }_0\right)\right)}^2}$ (1.4)

支架的径向力 $R_F$ 与直径D的关系可由式(1.5)代入式(1.6)中得到：

$\alpha =\arcsin \left(\frac{D\sin {\alpha }_0}{D_0}\right)$ (1.5)

$R_F=4\pi n\frac{L_0}{p_0}\left[\frac{G{I}_p}{C_3}\left(\frac{2\cos \alpha }{C_3}-C_1\right)-\frac{EI}{C_3\tan \alpha }\left(\frac{2\sin \alpha }{C_3}-C_2\right)\right]\tan \alpha$ (1.6)

其中n为编织线丝数量，系数 $C_1$ 、 $C_2$ 和 $C_3$，惯性矩I，极惯性矩 $I_p$，剪切模量G分别如下所示：

$C_1=\frac{\sin 2{\alpha }_0}{D_0}$ 、 $C_2=\frac{2{\sin }^2{\alpha }_0}{D_0}$ 、 $C_3=\frac{D_0}{\sin {\alpha }_0}$ (1.7-1.9)

$I=\frac{\pi d^4}{64}$ 、 $I_p=\frac{\pi d^4}{32}$ 、 $G=\frac{E}{2\left(1+\mu \right)}$ (1.10-1.12)

3. 有限元模型建立

3.1. 编织支架三维模型的建立及网格划分

由于支架的编织特性 [8]，需要从参数化的角度进行模型构建。基于已开发的可视化软件Constent，对编织支架进行离散化处理，导出含有支架编织路径上所有柱坐标点的集合，该导出文件后缀为txt。此外，该点集的密度可在软件Constent进行设置。如图3所示，在软件ANSYS Spaceclaim中，基于Python语言编写脚读取包含支架编织路径上的所有柱坐标点的导出文件，自动生成编织支架的三维模型，最终将三维模型导入ANSYS Workbench。

如图4所示，该编织支架的初始直径为15.4 mm，初始长度为30 mm，初始编织角度58.76˚，编织线丝的直径为0.2 mm，末端环式结构的曲率半径为1 mm。经过网格划分后的编织支架三维模型，采用梁单元结构，网格尺寸大小为0.01 mm，网格共4944个单元。

3.2. 编织支架的接触定义

由于编织支架在交错位置存在大量的非线性接触，使用连接器对编织线之间的接触进行定义将产生严重的接触渗透和过约束问题 [11]。与此同时，由于求解器需要迭代计算以达到收敛，所以一个非线性问题所需的求解时间将远远超过一个线性问题的求解时间，为了描述编织导线之间的数量庞大的接触以及减少接触面之间的互相渗透，本节采用ANSYS中自动接触检测的设定方法取代了传统的手工方式接触定义法，其中自动接触检测的接触模拟方法采用广义拉格朗日算法，编织线丝之间的接触类型为无摩擦接触，法向刚度设置为0.01。

3.3. 编织支架的材料模型

编织支架的材料为镍钛合金，如表1所示，该材料模型 [12] 在ANSYS软件中设置为经典弹性模型。

3.4. 压缩管模型的建立及网格划分

如图5所示，压缩管的长度为80 mm，外侧直径为16 mm，内侧直径为15.9 mm，管壁厚度为0.1 mm。网格划分采用线性和二次单元阶混合模式，网格尺寸为0.01 mm。

3.5. 边界条件

在径向压缩的过程中，压缩管的内壁与支架的外侧发生接触，该接触定义为无摩擦接触，压缩管沿着径向压缩3 mm的同时，限制支架的轴向移动和周向旋转。此时，支架作用于压缩管的径向力为与施加在压缩管的径向位移相反的作用力。

3.6. 支架径向压缩过程的模拟

如图6所示，当支架受到压缩管的作用径向压缩3 mm时，可以观察到支架的两端沿着轴向滑移并收缩，应力集中主要位于支架两端环式结构的弯折位置。支架的初始直径为15.4 mm，受到压缩后支架中间区域的直径为9.4 mm。为了与解析模型进行对比分析，该有限元仿真模型主要提取以下两个方面的数据：1) 几何结构方面：支架长度随直径的变化，其中支架长度的变化通过比较两端节点间的绝对长度与初始长度的差值定义；2) 力学性能方面：应力随支架长度的变化关系，其中应力为从支架的中部均匀选取8个节点提取的平均应力值；以及径向力与支架直径的变化关系。

4. 结果对比分析

如图7所示，当支架通过与压缩管作用被径向压缩3 mm时，支架的直径从15.4 mm压缩至9.4 mm，长度有30 mm增加至50 mm，编织角度从58.9˚减小至30.69˚。由解析模型得到的支架长度与直径的变化关系与仿真分析结果具有较好的一致性。

如图8所示，对比解析模型与仿真分析结果，等效应力与长度的变化关系表现出基本的一致性。当支架受到径向压缩3 mm时，其初始长度30 mm伸长至50 mm，最大等效应力达到54 MPa，该数值小于镍钛合金材料的屈服强度。

如图9所示，对于解析模型和有限元模型，在支架直径受到压缩从15.4 mm至9.4 mm的过程中，支架作用于压缩管的径向力逐渐升高。但是仿真得到的径向力要明显大于由解析模型计算得到的径向力，主要原因在于编织支架模型末端的环式结构，该结构相比末端开式结构具有更高的稳定性。

5. 结论

本文运用有限元方法，对具有编织特性的支架进行了三维建模，模拟了编织支架受到压缩管作用并产生径向压缩3 mm的过程，通过与经实验验证的解析模型对比分析，得出以下结论：

1) 该编织支架的有限元模型能够有效预测支架长度与直径的变化关系、等效应力与支架长度的变化关系，以及支架径向力与支架直径的变化趋势。

2) 支架在受到压缩管作用时，可以观察到支架的末端沿着轴向滑移并收缩，应力主要集中于末端环式结构的弯折位置附近。

3) 参考编织支架受到压缩管作用并产生径向压缩3 mm的过程，可以预测当编织支架受到狭窄血管壁的作用时，相比于末端开放式结构的编织支架，具有末端环式结构的支架能够提供更多的径向机械支撑。

4) 实际上在支架与血管壁的弹性限制之间取得平衡关系从而贴近血管壁的过程中，编织线丝之间的接触以及支架与血管壁的接触都存在一定程度的摩擦，本文还需进一步进行径向压缩实验评估摩擦因素对仿真环境的影响。

基金项目

国家自然科学基金资助项目(51405300)。

NOTES

^*通讯作者。

参考文献