有限元在心脏瓣膜支架设计中应用.docVIP

有限元在心脏瓣膜支架设计中应用 　　摘要： 通过SolidWorks建立某心脏瓣膜支架的三维模型，用Abaqus对瓣膜支架压握入鞘过程与疲劳情况进行分析，得到支架的应变分布和疲劳性能.结果可以判断支架设计的优劣，为支架优化设计提供参考. 　　关键词： 心脏； 瓣膜支架； 介入治疗； 有限元分析 　　中图分类号： R318.11；TB115.1文献标志码： B 　　引言 　　主动脉瓣膜狭窄是一种常见的心血管疾病，但对病人的生活造成严重影响.经导管主动瓣膜植入术（Transcatheter Aortic Valve Implantation，TAVI）治疗主动脉瓣膜狭窄已在世界范围内得到广泛应用.自2006年起，国内开始TAVI的基础研究，研制出球囊扩张式、自扩张式等不同类型的主动脉瓣膜支架.本文针对自扩张式的支架展开研究. 　　本文采用有限元法（Finite Element Method，FEM）计算心脏瓣膜支架压握入鞘和疲劳载荷下支架的应变分布与安全性能，为支架的优化设计提供参考. 　　1分析方法与模型 　　1.1分析步骤 　　根据支架生产加工到植入人体的过程，瓣膜支架有限元分析分为以下3个步骤：（1）模拟支架扩张定型过程；（2）模拟支架被压握至输送器的过程，即利用压握工具将支架压到指定尺寸；（3）模拟人体中瓣环收缩、扩张的循环对支架的疲劳寿命影响. 　　1.2几何模型 　　通过SolidWorks建立瓣膜支架的三维几何模型，见图1.瓣膜支架结构具有一定的对称性，在保证结果准确的基础上，考虑到计算时间和效率，选择1/6个支架进行分析.图 1瓣膜支架三维几何模型 　　压握工具和扩张工具均为理想的圆柱面，采用SFM3D4R单元建模. 　　1.3分析过程 　　利用Abaqus/Standard求解器求解计算，并进行相应的前、后处理. 　　1.3.1边界条件 　　在第一步的模拟中，利用一个圆柱形刚性管代替扩张工具.对刚性管施加位移载荷，使支架从原始尺寸开始扩张，一直到支架所要求的给定尺寸；然后去除刚性管，利用2个刚性面作为支架定型的模具，与支架接触，使支架最终定型. 　　在第二步的模拟中，利用一个圆柱形刚性管代替压握工具.对刚性管施加位移载荷，使支架从原始尺寸进行压握，一直到支架所要求的压握尺寸，模拟压握. 　　在第三步的模拟中???在支架瓣环处与支架升主段建立2个刚性面，代替与支架接触的血管.对刚性面施加位移载荷，模拟血管收缩与舒张对支架的作用. 　　1.3.2材料属性 　　瓣膜支架由NiTi形状记忆合金制成.NiTi形状记忆合金具有优良的形状记忆和超弹性性能，以及良好的抗腐蚀性和生物相容性.[1]在过去的30年里，NiTi合金在牙科、骨科以及介入治疗等医学领域得到广泛应用.由NiTi合金制成的支架具有良好的柔顺性，可以与复杂的人体内管道很好地贴合.本文分析采用Abaqus内建的NiTi材料超弹性模型，NiTi材料模型见图2. 　　图 2NiTi材料模型 　　1.3.3单元类型 　　在分析过程中，单元类型选择六面体非协调模式C3D8I（见图3）.其优点主要有：（1）克服剪切自锁问题，在单元扭曲较小情况下得到的应力应变结果很精确；（2）在弯曲过程中，厚度方向只需很少的单元就可得到与二次单元相当的结果，而计算成本却明显降低；（3）使用增强变形梯度的非协调模式，单元交界处不会重叠或开洞，因此很容易扩展到非线性、有限应变的位移. 　　（a）立体图 　　（b）平面图 　　图 3八节点六面体单元示意 　　1.3.4疲劳分析判断依据 　　支架在扩张后瓣环处、流出道附近分别与血管壁贴合，受到血管壁的作用力.在疲劳分析时，血管收缩与舒张对支架的影响可通过对代替血管的刚性面施加交替变化的位移进行模拟，由此引起的支架应变变化示意见图4.由图4可知，交变应变εa=εmax-εmin2 （1）平均应变εm=εmax+εmin2 （2）图 4应变变化示意 　　Pelton图[2]常用于评估NiTi材料疲劳性能，可以清晰地表达出测试点应变的分布状态，见图5.可知，若测试点分布在疲劳临界线下方，说明在若干次的疲劳试验后，支架发生断裂的可能性不大；反之，若测试点分布在其上方，则断裂极有可能发生，风险较大. 　　图 5Pelton示意 　　2分析结果 　　NiTi材料的破坏主要由拉伸应变造成，因此，扩张定型和压握过程将主要关注最大主应变；而在疲劳分析中，关注的则是最大交变应变. 　　2.1扩张定型 　　扩张定型支架应变分布见图6.在扩张定型过程中，支架整体应变水平较低，最大应变出现在网格节点内侧，与实际情况相符.NiTi材料的最大主应变不应超过12%，扩张定型中支架最大主应变为6.696%.因此，在