血管支架有限元优化设计
(完整版)血管支架设计
膨胀
3.支架优化设计
分析结果
膨胀
3.支架优化设计
分析结果
径向支撑力
3.支架优化设计
分析结果
压握
3.支架优化设计
分析结果
压握
3.支架优化设计
分析结果
柔顺
3.支架优化设计
分析结果
柔顺
3.支架优化设计
分析结果
柔顺
3.支架优化设计
2.血管支架设计
参数计算
ac b
参数 金属覆盖率 （Φ4mm） 轴向短缩率 径向回弹率 支架壁厚 显影效果 支架压握直径 均匀扩张程度 径向支撑力 柔顺性
数值 10%
≤ 2% ≤ 2% 0.004in 优秀 <1.1mm 优秀 良好 良好
参数 标称直径 长度 波形数 节长 节数
数值 4mm 20mm 7 5mm 4
2.血管支架设计
二维模型建立
AutoCAD
三维模型建立
SolidWorks
2.血管支架设计
支架建模II-SolidWorks
2.血管支架设计
支架建模II-SolidWorks
2.血管支架设计
支架建模III
SolidWorks
2.血管支架设计
支架建模III
Pro/E
2.血管支架设计
支架建模Ⅳ
Inventor
2.血管支架设计
支架建模Ⅳ
Inventor
3.支架优化设计
优化设计过程
径向支撑力 径向反弹比 轴向缩短率 轴向柔顺性
支架模型
有限元模型
求解
结果分析 N
设计 要求
利用有限元技术优化支架设计并应用
1序言某汽车厂研发的一款纯电动微型轿车，后桥为拖曳臂式半独立悬架，在样车路试阶段，汽车制动时后轮存在异响、抖动等问题，急需解决。
2原因分析汽车厂开发的后桥总成如图1所示。
经认真检查分析，在排除了轮毂、制动盘及钳体总成的原因后，我们认为制动时产生异响、抖动，可能有两个原因：一是钳体总成安装支架（见图2，以下简称支架）的刚度不够，制动时支架变形；二是支架的两耳面不在同一平面上，平行度超差。
这些原因造成钳体总成安装后，制动块与制动盘不平行，从而导致制动时出现抖动、异响。
支架的平行度经过测量为0.25mm，确实超差比较严重。
针对支架的刚度变形量，我们采用有限元技术分析刚性变形量的大小。
图1后桥总成图2安装支架3有限元CAE分析支架由6.5mm厚的Q235钢板弯曲而成，属于悬臂梁结构。
首先利用有限元CAE 技术对支架进行受力分析，支架用材料属性见表1。
表1支架用材料属性有限元CAE分析的主要步骤为：①对支架数模指派材料为Q235。
②划分网格，本零件划分网格类型为CTETRA(10)，单元大小为2.5mm，如图3所示。
③添加边界条件，约束及载荷。
④利用Nastran求解器进行计算，查看位移大小，最大变形位移量为0.1069mm，支架的弹性变形量严重超差，如图4所示。
图3划分网格图4位移量计算结果4产品优化设计为了解决支架弹性变形较大的问题，必须增加支架的刚度，刚开始提出的方法是在支架两边各焊接一个加强板，如图5所示。
经有限元受力分析，刚度变形最大位移量为0.0856mm，效果并不明显，问题没有得到有效解决，此方法行不通。
图5焊接加强板示意根据以上探索及分析，要想解决根本问题，必须重新设计支架来满足要求。
第一，重新设计支架外形，结合后桥总成的空间要求，给支架的外轮廓加上翻边，增加其整体刚度。
第二，增加材料的厚度，但厚度不能无限制增加，要考虑到产品要求、工艺要求以及成本等因素。
经过综合考虑，将支架的材料厚度增加到8mm。
有限元法在血管支架生物力学性能评价中的应用_刘磊
▲图 5
轴向短缩率：
fs = l0 － ls / l0
（ 2）
扩张不均匀率 Ｒ = dmax － dmin / dave2
（ 3）
弹性回缩率 κ = dave1 － dave2 / dave2
（ 4）
式中 fs 为轴向短缩率，ls 为扩张后长度，l0 为初始长度； Ｒ 为
Key words： finite element； endovascular stents； biomechanical performance； expansion simulation
血管支架成形术是将支架植入体内用来支撑血管和压 缩血管腔内斑块，防止血管回弹及内膜斑块向血管内突出的 一种介入手术。支架经输送系统送达病变部位后扩张，保持 血流畅通。血管支架作为人体植入物，必须具备优异的生物 力学性，包括： 扩张性、支撑性、柔顺性等，其中扩张性是重要 的指标，直接影响支架的植入过程以及术后的支撑性和血管 内再狭窄率。
向扩张，分析了 3 种截面、4 种尺寸的支架在扩张过程中主要 部位的残余应力； 周永恒等［2］建立球囊、支架模型，通过对球 囊内部加压扩张支架，分析了支架扩张过程的径向和轴向位 移、非线性回复量； 倪中华等［3］ 建立支架、球囊耦合扩张模 型，通过四个分析阶段： 支架外表面施加压强压缩支架、球囊 内表面加压扩张支架、压强保持、压强卸载，对支架变形机理 进行了研究，并通过实验验证了有限元模拟结果与支架实际 扩张后形貌的一致性； 王文雯等［4］将变形过程分为： 对压握 壳施加径向位移压缩支架，对球囊施加径向位移扩张支架， 撤销球囊三个阶段，对三种结构支架进行了扩张模拟，分析 其生物力学性能； 李红霞等［5］对球囊施加气压使支架扩张后 再卸载，通过提取扩张后模型的节点信息重新建立变形后的 血管、血栓、支架及流体的耦合模型，进行了支架内血流动力 学分析。国外学者在进行血管支架的有限元分析中，同样做 了大量的扩张模拟工作。Dumoulin 和 Cochelin、
人体血管支架有限元分析与结构拓扑优化共3篇
人体血管支架有限元分析与结构拓扑优化共3篇人体血管支架有限元分析与结构拓扑优化1人体血管支架有限元分析与结构拓扑优化随着现代医学的发展，血管支架已成为血管疾病治疗中不可或缺的一种工具。
血管支架可以通过膨胀和固定在动脉内部，从而恢复狭窄或闭塞部位的血流通畅。
因此，如何提高血管支架的稳定性和生物相容性已成为关注的焦点。
本文旨在探讨如何通过有限元分析和结构拓扑优化，提高人体血管支架的性能。
有限元分析的基本原理是将一个复杂的结构模型分解为小的单元，在每个单元内进行力学分析。
这种分析可以模拟不同的荷载状态和材料性质，从而评估结构的行为和性能。
在血管支架模型的有限元分析中，一个主要的问题是如何精确模拟支架材料和血管组织的非线性应力应变行为。
此外，由于支架植入后会受到血流和动脉脉动的影响，因此在分析中必须考虑这些因素的复杂效应。
一种有效的方法是使用仿真软件，在计算机中模拟血管支架的力学行为。
这种方法可以显示支架在不同荷载状态下的应力和变形，从而评估支架的性能。
这些结果可以用于优化支架的设计，以提高其性能和生物相容性。
例如，通过有限元分析，可以确定支架的形状、大小、横截面积和壁厚等参数，以最大限度地减少支架内部的应力和扭曲，从而提高其稳定性。
然而，即使在最优化的设计下，支架材料也可能不足以承受日常使用和长期暴露的逆境。
在这种情况下，我们可以采用结构拓扑优化的方法进一步优化支架的性能。
结构拓扑优化是一种在已有结构中寻找最优分布的方法，以最大限度地减少材料的使用量并提高结构的性能。
这种方法在血管支架设计中有广泛的应用，因为它可以减少支架内部的应力和材料冲击，并提高支架材料的生物相容性。
例如，我们可以使用结构拓扑优化来探索支架材料的排列，在保证结构稳定性的同时尽可能减少材料的使用量。
我们还可以使用形状和参数优化技术来优化支架的设计，以最大限度地减少支架内的应力和变形。
这些技术可以进一步提高支架的性能，使其适用于更广泛的应用场景。
基于生物力学性能的血管内支架弧梁单元结构优化设计_劳永华
中国组织工程研究与临床康复 第14卷 第30期 2010–07–23出版Journal of Clinical Rehabilitative Tissue Engineering Research July 23, 2010 Vol.14, No.30ISSN 1673-8225 CN 21-1539/R CODEN: ZLKHAH5581School of Bioscience and Engineering, South China University of Technology, Guangzhou510006, Guangdong Province, ChinaLao Yong-hua ☆, Studying fordoctorate, Lecturer, School of Bioscience and Engineering, South China University of Technology, Guangzhou510006, Guangdong Province, China yhlao@Correspondence to: Huang Yue-shan, Doctor, Associate professor, School of Bioscience and Engineering, South China University of Technology, Guangzhou510006, Guangdong Province, China bmbyshan@ Supported by: the Basic ScientificResearch of Central High Institutes, South China University of Technology, No. 2009ZM0075*Received: 2010-02-21Accepted: 2010-05-25基于生物力学性能的血管内支架弧梁单元结构优化设计*☆劳永华，支晓兴，林泽枫，李小您，岑人经，黄岳山Optimal design of intravascular stent arc-strut element structure based on biomechanical performanceLao Yong-hua, Zhi Xiao-xing, Lin Ze-feng, Li Xiao-nin, Cen Ren-jing, Huang Yue-shanAbstractBACKGROUND: Support the expansion of blood vessels in the implanted process of requires stent structure with complex biomechanical properties. The structure can achieve optimal overall performance by integrated multi-objective optimization. OBJECTIVE: Based on designs of various intravascular stent structures, the features parameters were extracted for optimization analysis of biomechanics.METHODS: Based on the analysis of biomechanical behavior: metal-based coverage, the maximal stent stress of expansion, maximal stent radius of expansion, axial recoil rate, radial recoil rate, finite element analysis software ADINA was used to establish 81 stent structures of closed, parallel, and open-arc strut element with different circumstantial unit number, arc-strut width and axial length of the stent unit element; the non-dimensional method and multi-objective optimum ranking was used to analyze the stent biomechanical behavior.RESULTS AND CONCLUSION: Results showed that mixture structure of stent with closed, parallel, and open-arc strut element can be designed, and that the structure possesses mainly large axial length of the stent’s closed, parallel strut element and small arc-strut width uniformly. It is significantly important for the stent structural design and mechanical behavior evaluation system.Lao YH, Zhi XX, Lin ZF, Li XN, Cen RJ, Huang YS. Optimal design of intravascular stent arc-strut element structure based on biomechanical performance.Zhongguo Zuzhi Gongcheng Yanjiu yu Linchuang Kangfu. 2010;14(30): 5581-5585. [ ]摘要背景：在扩张置入体内支撑血管的过程中血管内支架的结构需具备多种生物力学性能，支架结构可以通过综合多目标优化达到其整体性能的最佳状况。
血管支架结构有限元模拟与分析
大连理工大学硕士学位论文血管支架结构有限元模拟与分析姓名：王丽申请学位级别：硕士专业：计算机技术指导教师：吴迪20060611血管支架结构有限元模拟与分析ｍ。
所以，我们要改变既往“我国心血管病发病率和死亡率低于发达国家”的认识误区。
并且要清醒地认识到以冠心病为主的心血管病已成为我国一个重要的公共卫生问题，对冠心病的治疗方法和相关器械的研究已成为相关科学家们研究的热点之一。
图１．２粥样硬化斑块使冠状动脉腔变狭窄的示意图Ｆｉｇ．１．２Ｓｃｈｅｍａｔｉｃｄｉａｇｒａｍｏｆｃｏｒｏｎａｒｙａｒｔｅｒｙ’ｓｎａｒｒｏｗｉｎｇｃａｕｓｅｄｂｙａｔｈｅｒｏｓｃｌｅｒｏｔｉｃｂｌｏｃｋ１．２冠心病的治疗方法目前冠心病的治疗方法有药物治疗，外科手术治疗和介入性治疗三种。
（１）药物治疗：是指用硝酸酯类药、Ｂ受体阻滞剂、钙离子拮抗剂、转换酶抑制剂、调节血脂药、抗凝药物和中药等药物的治疗。
这种方法只能起到预防或短期的疗效。
（２）外科手术治疗：是指冠状动脉旁路移植术，即冠状动脉搭桥术（ＣｏｒｏｎａｒｙＡｒｔｅｒｙＢｙｐａｓｓＧｒａｆｔｉｎｇ，ＣＡＢＧ）。
这种手术一般比较复杂，需开胸，创伤大，恢复慢。
（３）介入性治疗：目前应用最广泛的是经皮穿刺冠状动脉腔内成形术（ＰｅｒｃｕｔａｎｅｏｕｓＴｒａｎｓｌｕｒｎｉｎａｌＣｏｒｏｎａｒｙＡｎｇｉｏｐｌａｓｔｙ，ＰＴＣＡ）和冠状动脉内支架术等。
ＰＴＣＡ也称为微创伤疗法。
这种治疗技术是在医学影像设备的导引下，利用穿刺针、导丝、导管和支架等器械所进行的诊断与治疗操作。
微创伤治疗技术是目前西方国家治疗各种威胁人类健康疾病的最流行手段，在所有的手术中微创伤介入治疗所占的比例已经超过８５％，微创伤治疗技术代表着国际上最高手术水平，亦是全世界医疗界发展的必然趋势。
与传统外科手术相比，进行微创伤介入治疗，无须开刀，只需局部麻醉，具有出血少、创伤小、并发症少、安全可靠、术后恢复快等优点，大大减轻了病人所承受的痛苦，降低了手术者的操作难度，手术时间及住院时间显著缩短，费用亦明显降低。
冠脉支架体外扩张过程的有限元分析
研究选用的支架材料型号分别为镁合金 ＷＥ ３（ ４ 泊松 比为 其施加反 向位移 ， 支架发生小距离的压缩 回弹 ； 扩张阶段对球囊 ０ ７ ， １Ｌ ． ）３ ６ 不锈钢和 Ｌ ０ ２ ６ ５钻铬合金（ 泊松 比为 ０ ）材料模 型采 施加径 向位移 ， ．， ３ 使支架 内径扩 至 Ｏ３ ｍ 然后对球囊施 加反向 ． ｍ， ０ 用 Ｖ ｎＭ ｓｓ ｏ ｉ 屈服准则和各 向同性等 向强化 准则来分析支架 的 ｅ 位移 ， 支架又发生小距离的扩张回弹。
的球囊结构 ， 使其膨胀时 ， 支架变形结果 良好 。一些科研人 员『 ５ Ｊ 定 到材料 、 几何和接触等复杂 的非线性 问题 ， 而支架结构的轴向又
量分 析了球囊扩张式 Ｂ ｅ ｃｔ支架扩张变形及其卸载后轴向 具有周期性 ， ＸＶ ｌ ｉ ｏ ｙ 因此为了减小计算量 ， 只选取 了支架 的一个几何单
１ ｄｐｅ ｅ ｃｍ Ｉ ｉｌｅ ｎｏ ｒ ｐ ￣ ｓａ ｍ ｔｆｉ
ｓ ｅｌ ｈｌ
大变形结果。 通过单轴拉伸试验 ， 获得各种材料的特性参数 ， 如表
１ 所示 。
Ｉ Ｉ－ ｉｐａ ｅ ｎ ｆ ａｌ０ －Ｉ ｄ ｓ ｌｃ ｍｅ ｔ ｂ ｌ “ － ｏ ０
郭飞飞 冯海全 江旭 东 韩 青松
（ 内蒙古 工业 大学 机械 工程学 院 ， 和浩 特 ００ ５ ） 呼 １０ １
Ｔｈ ｉ ｉ ｅｍｅ ｔ ａｌｓｓ ｏ ｒ ｎ ｒ ｔ ｎ ｔｏ Ｅ ｐ ｎ ｉ ｎ Ｐｒ ｃ ｓ ｅ Ｆ ｎｔ Ｅｌ ｅ ｎ Ａｎ ｙ ｉ ｆ Ｃｏ ｏ ａ ｙ Ｓ ｅ ｔ ｎ Ｖｉ ｘ ａ ｓｏ ｏ ｅ ｓ ｉ ｒ
的基础上发展起来的又一种治疗 冠心病 的手段 ， 使用 时 ， 在球囊
缩短和径 向回弹效果 ， 并与体外实验数据进行 了对 比。一些科研 人员 建立 了球囊扩 张式冠脉支 同建模策 略对支架扩张变形的影响。 一些科研人 员 通过分析球囊扩张式支架的一个重复结构单元 比较 了不 同连 接形式 的三种支架在扩张过程中发生的变化 ， 出了小单元分析 提
血管内支架设计和检测中必不可少的助手一有限元分析法
。
每 年大 约 有 7 0 fL 美 元 的 全 球 市 场 而 且 有 不 断 增 加 的
，
趋势
。
从 发 展 历 史 来 看 研 究 者将 主 要 的精 力 放 在 了提
，
支 架 的概 念是D o t t e r e
第 第
l 9 6 9 年首先提 出的
Ⅲ
；
，
用来描
P a lm
”。
高支 架 的生 物相容性 机 械力学性 能和植 人 后 的再 狭
．
宋金 子
．
国 家食 品 药 品 监 督管 理 局 济 南 医 疗 器 械 质 量 监 督检 验 中心 ( 山 东
【 摘
济南
．
2 5 0 10 1 )
。
要 】 介 绍 了 国 内外 应 用 有 限 元 分 析 方 法
．
在 不 同 释 放 方 式 支架
、
支 架耐 疲 劳 性 能 及 支 架 柔 顺 性 等 方 面 研 究 的进 展
维普资讯
… 0
+
u
‘
e
$
’
Ⅻ
r 1
r 1
文章 编 号
：
1 6 7 1 7 1 0 4 (2 0 0 8 )0 4 0 2 8 4 0 3
—
— —
血 管 内 支 架 设 计和 检 测 中必 不 可 少 的助 手 有 限 元分析 法
一
【 作
者】 张博
．
万敏
．
姚秀君
s
，
f a t ig
du
t f i
，
ra n
b ilit y
it e
a n
d th e f l
血管支架结构有限元模拟与分析
随着我国经济水平的发展、人民生活水平的提高、饮食结构的改变及人口迅速老龄 化，心血管病的发病率和死亡率呈上升趋势，是全球上升较快的国家。据统计，我国心 血管病和脑卒中的总发病率和死亡率已超过许多发达国家。１９９９年我国农村和城市男性 ３５～７４岁人群中心血管病总死亡率分别为４１３／１０万和３８９／１０万，冠心病死亡率分别为 ６４／１０ＴＹ和１０６／１０ＴＹ，脑卒中死亡率分别为２４３／１０万和２１７／１０万；而同期美国同年龄段男 性人群心血管病总死亡率为３６０１１０万，冠心病和脑卒中死亡率分别为２３０１１０万和４１／１０ＴＪ
ＣＡＥ ｓｏｆｔｗａｒｅ ＡＮＳＹＳ ｉｓ ｕｓｅｄ ｔｏ ａｎａｌｙｚｅ ｔｈｅ ｄｉｌａｔａｔｉｏｎ ｏｆ ｔｈｅ ＩｌｅＷ ｓｔｅｎｔ．Ａｆｔｅｒ ｓｅｖｅｒａｌ
ｍｏｄｉｆｉｃａｔｉｏｎｓ，ａ ｎｅｗ ｔｙｐｅ ｏｆ ｓｔｅｎｔ ｗｉｔｈ ｎｏ－ｄｏｇｂｏｎｉｎｇ ｉｓ ｄｅｓｉｇｎｅｄ．
Ｋｅｙ Ｗｏｒｄｓ：Ｅｎｄｏｖａｓｃｕｌａｒ Ｓｔｅｎｔ；Ｓｔｒｕｃｔｕｒｅ Ｄｅｓｉｇｎｉｎｇ；Ｆｉｎｉｔｅ Ｅｌｅｍｅｎｔ Ｍｅｔｈｏｄ
Ａｂｓｔｒａｃｔ
Ｃｏｒｏｎａｒｙ ｓｔｅｎｔｓ ａｒｅ ｇｅｎｅｒａｌｌｙ ｔｕｂｕｌａｒ－ｓｈａｐｅｄ ｅｘｐａｎｄａｂｌｅ ｄｅｖｉｃｅｓ ｗｈｉｃｈ ｆｕｎｃｔｉｏｎ ｔｏ ｈｏｌｄ
ｏｐｅｎ ａ ｓｅｇｍｅｎｔ ｏｆ ａ ｂｌｏｏｄ ｖｅｓｓｅｌ．Ｔｈｅｙ ａｒｅ ｐａｒｔｉｃｕｌａｒｌｙ ｓｕｉｔａｂｌｅ ｆｏｒ ｕｓｅ ｔｏ ｓｕｐｐｏｒｔ ａｎｄ ｈｏｌｄ
血管支架结构有限元模拟与分析
ｍ。所以，我们要改变既往“我国心血管病发病率和死亡率低于发达国家”的认识误区。 并且要清醒地认识到以冠心病为主的心血管病已成为我国一个重要的公共卫生问题，对 冠心病的治疗方法和相关器械的研究已成为相关科学家们研究的热点之一。
基于HyperStudy的心血管支架结构优化
基于HyperStudy的心血管支架结构优化Optimization of Stent Structure based on HyperStudy张雯,刘祥坤1高超2(1上海微创医疗器械（集团）有限公司，上海，201203)(2上海惠骋实业发展有限公司，上海，200436)摘要：血管内支架的力学性能研究对解决在支架植入后的破坏问题非常重要。
使用常规方法研究支架性能会有很大困难，成本昂贵。
随着有限元法的发展，在支架设计和性能检测中有限元已成为不可或缺的工具。
为了提高支架设计的性能和效率，有限元技术及专业CAE软件结合，对冠脉血管支架进行多学科结构优化设计。
本文将支架几何模型转化为参数化有限元模型，建立morph后的形状变量模型，提交HyperStudy 采用多种优化策略进行针对前两种工况的形状优化，并最终进行验证。
关键词：血管支架，多学科结构优化，HyperStudyAbstract:Mechanical properties of intravascular stent is very important in solving the problem of damage after stent implantation. The study of stents properties will be difficult and expensive using traditional methods. With the development of finite element method, the new method has become an indispensable tool in the stent design and performance testing. In order to improve the performance and efficiency of the stent design, multidisciplinary structural optimization with FEM and professional CAE software should be used in coronary stent design. This paper will convert the geometry models into parametric FEM model, establish a shape-variable model, submit to HyperStudy for shape optimization of the first two working conditions, and validate the result.Key words:stents, multidisciplinary structural optimization, HyperStudy1 引言血管内支架放置术已广泛应用于心血管疾病介入性治疗，但支架在植入后的破坏问题还需要研究解决，因此支架的力学性能研究非常重要。
基于有限元分析的综合支架结构设计与优选
今天， 对工程结构 的分析与计算采用近似的数值解
已成 为现 实又非 常 有效 的手 段 ， 中有 限 单值计算方法 。通过有 ２ ］ 限元法设计综合支架能够快速、 高效 的满足设计要 求。本文采用 了 Ｓｌ ｗ ｒｓ０６软件 中的 Ｃ ｓ ｏ ｏｄ ｏｋ２０ ｉ ｏｍ ． ｓｏ ｓ ｗ ｒ 模块 ， 中进行有限元建模时只对综合支架 ｋ 其 中安装部分的标准件作 了简化， 整体支架保持原始 状态， 这使得计算结果更加接近实际的受力情况 ， 分 析也更加准确可信 。
Ｆ ｂ． ｏ ８ ｅ ２ ０
基于有限元分析的综合支架结构设计与优选
王 艳 彩
（ 廊坊职业技术学院 机械 工程 系，河北 廊坊 ０ ５ ０ ） ６０ ０ 摘 要 ：本文在分析海洋平 台防爆加热器橇块综 合支架结构特 点 的基础 上 ，通过有限元法对 综合支架进行 强度
能力分析 ，找到该支架结构设计上存在的问题并作 了相应 的改 进设计 ，使得综 合支架 的结构 在满足强 度和刚度
Ｏ 引 言
１ 防爆 电加热器橇 块 的综合 支架结构 的设计 方法
１ １ 综合 支架结构组 成 ．
防爆电加热器是海上平台中不可缺少的部分 ， 防爆电加热器橇块则是海上平台装备中的一个机械 模块。海上环境条件恶劣并且复杂多变， 橇块 中作
为管线 支撑结 构 的管 道支架 必须要 有足够 的强度 和
的同时减轻 了整体的重量 ，在材料 的利用上趋 于优化 的设计 ；通过对 综合 支架 的有限元模型施加不 同的惯性力 ， 得 出不 同载荷下综合支架应力变化规律 ，从 而确定 出综合支架所能承受的最大惯性力 。 关键词 ：综合支架 ；有 限元法 ；优化 ；承载能力
血管内支架结构的设计优化
Ｌ ｂｏ ａｏｒ ｏｏｎｃＩｆ ｒ ａｉｎ Ｔｅ ｈ ｏ ｙ， ｕｔ ｉａ Ｎｏｒ ｌＵｎｉｅｒｉ Ｇ ｕ ｇ ｈ ａ ｒ ｔ ｙｏｆ Ｐｈ ｔ ｉ ｎｏ ｍ ｔ ｃ ｏｌｇ Ｓｏ ｈ Ｃｈｎ ｏ ｍａ ｖ ｓｔ ｙ， ａｎ ｚ ｏｕ，１ ０６ ５ ００
表 １三 种 支 架 设 计 方 案
Ｔａ １ Ｔｈ ｅ ｅｓｇｎ ｈ ｔ ｎ ｂ． ｒ ｅ ｄ ｉ ｓｏｆｔ ｅ ｓ ｅ ｔ
编号
ａ
沟槽 长度（ ｍ
２Ｏ ．
１８ ． １５ ．
４０ ．
４２ ． ４５ ：
４Ｏ ．
４０ ． ４０ ．
的关系ｐ】 一 这要求有获得 支架结构特征及变形行为 的 。
方法及优化支架形状结构的科学理论 。
支架实体的二分之一（ 支架网梁的厚度为０ Ｉ， 图２ ｏ ．１Ｔ ９Ｔ Ｉ １
由于血管内支架 的价格高及其制备技术复杂性 ， 采用试错实验确定支架设计参数需消耗大量资源『】 ５。 ． 近 ６ 年来，人们研究 出用有限元分析法计算支架的力学行 为， 并指导支架的形状设计， 模拟结果与真实应用情 况 很接近ｌ】 ７，已被美国 Ｆ Ａ作为支架产品可靠性 的评判 Ｉ ８ Ｄ 依据之一 。 本文用有限元法模拟了球囊扩张式不锈钢血 管内支架在内压力作用下的变形过程， 并通过计算机仿 真显示出支架尺寸及应力大小的变化， 这种仿真模拟方 法可使血管 内支架的结构设计方案更具科学性。
及变形应 力。 模拟结 果表 明通过改变支架结构 的尺寸参数可获得均匀扩张的设计方案 ． 并探索了降低支架变形应力大小的 方法。 【 关 键 词】 血管内支架：有限元分析：结构设计
面向个性化需求的医用血管支架数字化设计平台及其关键技术的研究(可编辑)
面向个性化需求的医用血管支架数字化设计平台及其关键技术的研究东南大学博士学位论文面向个性化需求的医用血管支架数字化设计平台及其关键技术的研究姓名:陈功申请学位级别:博士专业:机械制造及其自动化指导教师:易红20060704摘要摘要随着人民生活水平的提高,心血管疾病发病率越来越高,心血管狭窄引起的各种血管疾病已经成为严重危害人类身体健康的疾病。
在冠心病的各种治疗方法中,支架植入术因其创伤小、效果好,成为当前治疗心血管狭窄的新型方法。
本文以血管支架作为研究对象,研究面向个性化需求的血管支架数字化设计平台及其关键技术,重点研究血管狭窄部位个性化数据的获取、虚拟环境可视化和工程分析集成技术,开发应用于血管支架的数字化设计原型系统。
本文将如何快速、低成本地设计并制造出高质量、面向具体病人个性化数据的血管支架作为研究的着眼点,把网络化、数字化手段引入血管支架产品设计和制造过程中。
本文的具体研究内容、成果和创新点如下:研究并建立了面向病人血管个性化数据的血管支架数字化设计平台。
应用网络化、数字化设计技术改造血管支架的生产流程,通过一个面向个性化制造任务的信息共享平台,将病人、医院和生产厂家以及原材料供应商等多部门通过网络技术和信息技术整合在一起,以满足血管支架个性化、单品种、变批量生产的需求;在分析其制造特点和需求的基础上,提出并实现了面向个性化需求的血管支架数字化设计平台的平台体系结构,包括网络拓扑结构、功能框架结构、技术体系结构以及工作流程。
基于高斯灰度模板的血管边界分割方法,研究并开发了基于血管影像的支架参数获取工具。
提出了一种针对图像的基于高斯灰度模板的血管边界分割方法进行血管的边缘检测和跟踪;对血管狭窄程度的评价体系进行了研究并实现了血管直线狭窄率测量工具;提出了一种基于平面图形主惯性轴求法的血管轮廓长、短轴获取方法,进行了血管面积狭窄率的评估和沿轴向的分析,如管腔面积、斑块面积、偏心度、圆形度、最小外接矩形等;研究并实现了基于的血管直径和病变长度的测量工具;研究了协议的实现方法,提取其中的高精度信息,用于血管测量和建模。
缺血性脑卒中取栓支架优化设计与有限元分析
Ｃｏｒｒｅｓｐｏｎｄｉｎｇ ａｕｔｈｏｒ： ＣＨＥＮＧ Ｙｕｎｚｈａｎｇ （ Ｅ⁃ｍａｉｌ： ｃｙｚ２００８＠ ｕｓｓｔ ｅｄｕ ｃｎ）
【 Ａｂｓｔｒａｃｔ】 Ｏｂｊｅｃｔｉｖｅ Ｍｅｃｈａｎｉｃａｌ ｔｈｒｏｍｂｅｃｔｏｍｙ （ ＭＴ） ｉｓ ａｎ ｅｆｆｅｃｔｉｖｅ ｍｅｔｈｏｄ ｔｏ ｔｒｅａｔ ａｃｕｔｅ ｉｓｃｈｅｍｉｃ
同筋厚数值的支架在压握和释放过程中均是安全有效的，支架支撑单元的弯曲处最易发生支架疲劳失
效。 支架的径向支撑力和最大主应变峰值随着筋厚的增加而增加，支架筋厚数值对血管表面 Ｍｉｓｅｓ 应
力影响不大。 研究结果在理论上可为取栓支架的优化设计提供一定的依据。
关键词 急性缺血性脑卒中；取栓支架；有限元分析；径向支撑力；最大主应变
ｓｔｒｏｋｅ （ ＡＩＳ） ． Ｔｈｅ ｃｌｉｎｉｃａｌ ｕｓｅ ｏｆ ｓｔｅｎｔ ｒｅｔｒｉｅｖｅｒ ｉｓ ａ ｃｏｍｐｌｉｃａｔｅｄ ｐｒｏｃｅｓｓ ｏｆ ｍｅｃｈａｎｉｃａｌ ｃｈａｎｇｅｓ． Ｃｅｒｔａｉｎ ｓｔｕｄｉｅｓ
ｈａｖｅ ｆｏｕｎｄ ｔｈａｔ ｔｈｅ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｏｆ ｔｈｅ ｓｔｅｎｔ ｒｅｔｒｉｅｖｅｒ ｈａｖｅ ａｎ ｉｍｐｏｒｔａｎｔ ｉｍｐａｃｔ ｏｎ ｔｈｅ
Ｒｅｓｕｌｔｓ Ｔｈｅ ｍａｘｉｍｕｍ ｐｒｉｎｃｉｐａｌ ｓｔｒａｉｎ ｏｆ ｔｈｅ ｓｔｅｎｔｓ ｗｉｔｈ ０ ０４ ｍｍ，０ ０５ ｍｍ ａｎｄ ０ ０６ ｍｍ ｔｈｉｃｋｎｅｓｓ ａｒｅ ｌｏｃａｔｅｄ
ａｔ ｔｈｅ ｂｅｎｄｉｎｇ ｏｆ ｔｈｅ ｂｒａｃｋｅｔ ｓｕｐｐｏｒｔ ｕｎｉｔ，ｗｈｉｃｈ ａｒｅ ２ ８２７％，３ ４０３％ ａｎｄ ４ ７６８％，ｆａｒ ｌｅｓｓ ｔｈａｎ ｔｈｅ ｆｒａｃｔｕｒｅ
人体血管支架有限元分析
可降解支架的缺陷
1，高分子可降解 支架材料的弹性模 量低。支架体积臃 肿，容易产生不良 的血液涡流，进而 可能引起再狭窄的 发生
2，纯铁支架的降 解速度太慢
性能研究意义
3，镁合金支架在 降解过程中的支撑 能力有待于进一步 改善
为了对可降解支架的性能特点进行深入了解，并提出更好的
支架设计，有必要从生物力学的角度进一步研究可降解支架
2021/6/16
优化后有限元模型
优化前有限元模型
24
柔度分析和支撑性能对比
柔度分析表明，Sopt与Sori比较，降低了14%的柔度
2021/6/16
Sopt
25
Sori
优化后的Conor支架设计
2021/6/16
拓扑优化可以对新型的 储库式药物洗脱支架平台进 行结构优化。在不影响支架 平台载药性能的前提下，经 拓扑优化后的单个支撑筋的 结构刚度提高了14％。
2021/6/16
4
介入治疗的背景
PTCA+支架植入术
2021/6/16
5
支架植入术应用的领域
心脑血管
外周血管
非血管
冠颅 状内 动动 脉脉
颈肾胸肢静髂
动动腹体脉动
脉脉主动
脉
动脉
脉
肝胆食气尿 脏道道管道
2021/6/16
6
主要内容
• 介入治疗的背景 • 有限元法与验证实验 • 模拟测试 • 问题和未来应用
数，实现多次而快速的设
计调整；
Ansys具有与Pro/E的专用
接口，可以使得用Pro/E建
的几何模型能够无损地传
输到Ansys中，从而确保
后续有限元模型的精确性
。
可降解聚合物血管支架结构优化设计
可降解聚合物血管支架结构优化设计赵丹阳;刘韬;李红霞;王敏杰【摘要】Due to the low stiffness of polymers, polymeric stent has lower radial support capability compared to metallic stent. Therefore, the width and thickness of the stent are usually increased to improve its radial support capability, which can not only reduce the flexible performance of the stent and the area of the vascular lumen, but also increase the surface coverage and thus increase the risk of in-stent restenosis. In order to design polymeric stent with smaller strut width and thickness and improve its radial support capability, an optimization method combining with Kriging surrogate model and finite element method was used to optimize the geometries of stent. Kriging surrogate model was used to construct the approximate function relationship between design objectives and design variables. Optimized Latin Hypercube Sampling method was used to select the initial sample points. EI function was used to balance global and local search and tend to find the global optimal solution. As an example, ART18Z polymeric stent was studied in this paper. Firstly, the strut width and thickness of the stent were respectively reduced by 0.02 mm, and then the optimization method was used to optimize the key geometric parameters of ART18Z stent. The numerical results show that the overall service performance of ART18Z stent was improved after optimization and the proposed optimization method can be effectively applied to the optimal design of the polymeric stent.%聚合物血管支架由于材料刚度较低导致其径向支撑能力相对于金属血管支架较弱,通常采用增大支架筋宽和厚度的方式来提高其径向支撑能力,但这不仅会降低支架的柔顺性能,减小血管管腔获得面积,还会增大表面覆盖率,从而增大支架内再狭窄的风险.为了设计出具有较小筋宽和厚度的聚合物血管支架,提高其径向支撑能力,本文采用一种将Kriging 代理模型和有限元方法相结合的优化方法来优化支架的结构.采用Kriging代理模型建立设计目标和设计变量之间近似的函数关系,采用优化拉丁超立方抽样方法选取初始样本点,采用EI函数平衡局部和全局搜索,以便获得全局最优解.选取ART18Z聚合物支架作为算例,首先将支架的筋宽和厚度各减小0.02 mm,然后采用优化方法优化ART18Z支架的几何结构参数.数值结果表明,优化后ART18Z支架的综合服役性能得到改善,文中提出的优化方法能有效地应用于聚合物血管支架的优化设计.【期刊名称】《力学学报》【年(卷),期】2017(049)006【总页数】9页(P1409-1417)【关键词】可降解聚合物,血管支架,力学性能,Kriging;代理模型,优化设计【作者】赵丹阳;刘韬;李红霞;王敏杰【作者单位】大连理工大学机械工程学院,大连 116023;大连理工大学机械工程学院,大连 116023;大连理工大学机械工程学院,大连 116023;郑州大学微纳成型技术国家级国际联合研究中心及河南省微成型重点实验室,郑州 450001;大连理工大学机械工程学院,大连 116023【正文语种】中文【中图分类】R318.01随着微创介入治疗技术的发展，血管支架植入术已经发展成为治疗心血管狭窄率疾病最有效的方法之一.可降解聚合物血管支架是血管支架的重要发展方向，它避免了永久性金属支架作为永久异物所导致的炎症反应、平滑肌增生和迁移、内膜增生以及细胞外基质的聚集等病变引起的支架内再狭窄和远期安全性问题[1-8].然而，聚合物血管支架由于聚合物材料本身的弹性模量较低导致支架径向支撑能力相对于金属血管支架较弱，在球囊扩张并卸载后支架产生较大的径向弹性回缩，减小了支架内管腔获得面积和血管内血流通过量，并影响血管管腔重塑[9].目前，通常采用增大支架筋宽和厚度的方式来提高聚合物支架的径向支撑能力，如Ideal支架的厚度是0.2mm，Igaki-Tamai支架和ART18Z支架的厚度均为0.17mm，BVS支架和DeSolve支架的厚度都是0.15mm[10-11].然而，增大支架的筋宽和厚度不仅会降低支架的柔顺性能，减小血管管腔获得面积，而且较大的表面覆盖率还会增大新生内膜增生几率，从而增大支架内再狭窄的风险.所幸的是，支架的结构也会对支架的力学性能产生影响，可以通过改进支架的结构设计来改善支架的力学性能. 目前，血管支架结构设计的研究主要依赖于实验和临床，存在许多的困难和局限性，不仅费时，难度大而且成本较高，同时，某些重要参数，如支架结构内部的应力和应变分布等也较难通过实验获取.随着计算机科学技术的发展，数值计算技术模拟分析血管支架行为的研究解决了仅依靠实验存在的困难和局限，近年来被广泛采用.目前基于数值计算的支架结构设计主要是结构参数的对比分析，即对比分析几款不同结构的支架设计，或同款结构几组不同尺寸参数支架设计的力学性能，从中选取力学性能“最好”的一个设计作为支架的“最佳”设计[12-21].这种方法易于研究影响支架设计的因素，但只能在有限个设计中选取“最佳”设计，不能在设计区间内找出全局最优设计.并且，支架的有限元计算涉及到大变形和弹塑形共存等非线性问题，单个计算量较大，而要找到较好的支架设计，需要做大量的对比分析，其计算量将十分巨大.与对比分析法相比，代理模型的方法能建立设计变量和设计目标之间难以显性表达的近似函数关系，将复杂的工程问题转化为可求解的数学问题，大大降低计算成本.将数值计算与代理模型相结合的优化算法能更简便快捷地对支架进行优化.Atherton和Bate[22]用代理模型的方法研究了支架对血管的壁面剪切应力的影响，优化了支架结构.Harewood等[23]用代理模型研究了支架一个单环的径向刚度.Clune等[24]利用响应面技术优化了支架的抗疲劳特性.在众多代理模型技术中，Kriging代理模型在对未知信息的模拟精度上具有明显的优势，并且其计算模型更加简洁，能够更加方便地应用于工程问题的研究.Kriging代理模型包含参数模型和非参数随机过程两个因素，因此，比仅具有参数模型的插值方法[25-26]在具体应用中更加灵活，比非参数化模型具有更强的全局预测能力.Li等[27-28]利用基于Kriging代理模型的自适应优化算法分别优化了支架结构和支架涂层，消除了支架扩张过程中的狗骨头效应，延长了药物释放的有效时间.李红霞等[29]采用与有限元技术相结合的Kriging代理模型对基于不同扩张模拟方式的支架进行了优化.Pant等[30]采用Kriging代理模型优化了支架力学性能、药物释放效果和柔顺性能.Srinivas等[31]对简化的支架二维稳态流动模型进行了优化.利用代理模型技术，使得血管支架更为复杂和细致的研究成为可能.综上所述，针对目前可降解聚合物血管支架普遍存在的径向支撑能力较弱的问题，为了设计出具有较小筋宽和厚度，同时具有较好支撑性能的聚合物血管支架，本文以ART18Z聚合物支架为研究对象，先将其筋宽和厚度各减小 0.02mm，即支架的筋宽减小15.38%，厚度减小13.33%，然后采用改进的Kriging代理模型的优化方法，以减小支架扩张后的径向弹性回缩为设计目标，对其结构进行优化设计，以达到改善ART18Z聚合物支架服役性能的目的.以长 13.75mm，厚 0.17mm，外直径 3.36mm的ART18Z支架为研究对象.带有直桥连接单元的ART18Z聚合物支架的三维模型如图1所示.采用血管、血栓、支架和球囊的三维分析模型进行有限元分析，利用ANSYS软件包分析聚合物支架在狭窄血管内的扩张过程，其三维有限元模型如图2所示.由于整个模型具有对称性，故在此研究中选取了整个模型的1/12(周向1/6、轴向1/2)进行模拟计算.在支架扩张过程计算的有限元模型中，血管长7.9mm，厚度保持一致，均为0.15mm；血栓斑块近端的厚度为0.48mm，远端的厚度为0.08mm，在扩张开始的时候血栓斑块并没有与支架接触；球囊长7.75mm，厚0.05mm. PLLA具有良好的生物相容性和可控的降解速率，在可降解聚合物血管支架的制备上具有广阔的应用前景.ART18Z支架材料为PLLA，其弹性模量为3363MPa，泊松比0.45，屈服强度40MPa，等效为双线性、各向同性弹塑性材料.血管和血栓弹性模量分别为1.75MPa和2.19MPa，泊松比均为0.499，等效为线性、各向同性的近似不可压缩材料.球囊的弹性模量参数为C10=1.0688MPa，C01=0.710918MPa，泊松比0.495，等效为超弹性橡胶材料.在进行有限元分析中，血管、血栓和支架均采用8节点Solid 185实体单元，球囊则采用4节点Shell 181壳单元.模拟过程中的材料性能参数来源于文献[32].材料属性如表1所示.支架扩张是一个非线性大变形过程，本文采用具有较高计算精度的六面体网格划分模型.考虑球囊和支架的接触，以及支架和血栓的接触.为模拟支架在球囊膨胀作用下扩张和球囊卸载后支架的变形过程，在球囊上施加两步径向位移载荷：第一步，将支架的外直径扩张到血管的内部直径；第二步，将球囊卸载到最初的状态.血管、血栓、支架和球囊的对称部分施加对称约束，球囊远端约束轴向刚体位移，血管、血栓和支架的远端自由.在扩张过程中，支架在球囊膨胀的作用下逐渐被扩张，将狭窄血管撑开.当球囊卸载后，支架由于塑性变形仍然处于扩张状态，从而支撑狭窄血管.图3描述的是聚合物支架在狭窄血管内的膨胀过程，其中(a)∼(c)指的是支架从原始状态逐渐扩张到最大程度的过程，(d)∼(e)描述的是球囊卸载后的过程.径向回缩率(radial recoil,RR)是血管支架扩张性能的重要技术指标.由于PLLA材料刚度较低，导致PLLA支架的径向支撑能力较弱，支架在球囊卸载之后会出现较大的径向弹性回缩，这将减小血管内血流通过量，并影响血管管腔重塑.此外，支架的覆盖面积较大还会造成血管壁的内膜增生，因此通常将其限制在20%以下.本文先将支架的筋宽和厚度各减小0.02mm，然后以支架的关键结构尺寸为设计变量，以支架扩张后的径向弹性回缩为设计目标对支架进行优化设计，优化问题可定义为式中，Rloading是支架完全扩张时的半径，Runloading是卸载后支架的半径，X 是支架的设计变量，包含支架的几何参数a和b，如图1所示，图中a和b分别是菱形孔的长和宽；Xu和Xd分别是设计变量的上下限，且2.2mm6 a 62.6mm，1.4mm6 b 61.6mm.在工程优化设计中，采用代理模型建立设计变量与设计目标之间难以显性表达的函数关系，可将复杂工程问题转化为可求解的数学问题，大大降低计算成本.采用基于改进的Kriging代理模型的优化算法对血管支架进行优化设计过程的流程如图4所示，其主要过程为：(a)定义优化问题，即定义设计变量、设计目标和约束条件；(b)通过优化拉丁超立方取样方法获取初始样本点，并通过有限元分析计算每个设计样本点，得出每个样本点的响应值(即径向回缩率)；(c)基于样本点信息，Kriging代理模型采用二次多项式响应面形式来构建设计目标和设计变量之间近似的函数关系；(d)选择样本点中响应值最小的点作为优化程序执行的初始点，并选择优化算法来获得基于近似函数关系的最优设计；(e)检查收敛条件，当满足收敛条件时，优化程序停止.(1)Kriging代理模型Kriging代理模型是一种具有“统计性”的近似技术，包括回归部分和非参数部分，由多项式和随机分布组成式中，β为回归系数，f(x)为x的多项式，提供模拟的全局近似；z(x)为随机分布的误差，提供对模拟局部偏差的近似，具有如下统计特性式中，xi和xj是训练样本中的任意两个点，R(θ,xi,xj)是带有参数θ的相关函数，表征训练样本点之间的空间相关性.这里采用高斯相关函数式中，n为设计变量的数量，和为训练样本点的第k个分量.给定训练样本S=[x1,x2,···,xni]和对应的响应Y=[y1,y2,···,yns]，在式 (2)假定的基础上，任意一个待测点x*的响应值可由已知训练样本响应值Y的线性组合来表达则误差为把式(2)代入式(6)得式中，F=[f1,f2,···,fns],Z=[z1,z2,···,zns].为保证模拟过程的无偏性，误差的均值应该为0，即可得，FTc−f=0.则式(7)的预测方差为式中代表新样本点x*和各样本点的空间相关性.通过最小化预测值的预测方差来确定方差系数c，可以得到至此，可以根据式(12)得到任意新点x*处的预测值(2)取样方法采用优化拉丁超立方抽样方法取样，该方法由拉丁超立方抽样方法改进而来的.拉丁超立方取样的主要优点是确保选出的样本点可以代表样本空间中的所有部分，其次是这种取样方式获取的样本点数目不受变量维数限制，可以是任意整数，使用起来方便.但因为拉丁超立方取样具有随机性，会使得获取的样本点时好时差，因此，本文采用的是由Joseph等[33]提出的改进的拉丁超立方取样，即优化拉丁超立方取样.它是在拉丁超立方取样的基础上，以最小化样本点中点与点之间的两两相关性，并同时最大化点与点之间的距离为标准，对样本点的空间分布重新进行优化.图5为拉丁超立方取样及优化拉丁超立方取样图，通过比较可以看出优化拉丁超立方取样产生的样本点效果更好.(3)EI函数最大化期望提高是考虑预测值和预测方差加权的一种加点方法.期望提高(expected improvement，EI)是指计算在一个给定点的响应值目标提高的概率.对于任意一个设计点 x,在计算其响应值 y(x)之前，y(x)是未知的，但kriging能够预测它的均值和均方差σ2(x).假设当前最优设计的响应值为Ymin，则该点的响应值目标的提高I(x)=Ymin−y(x),它服从正态分布，则其概率密度函数为那么响应值目标提高的期望值为分部积分可得式中，Φ和ϕ分别为正则化概率分布函数和概率密度函数.式(16)右端包括两项之和，第1项为预测值与当前最优值Ymin之差乘以提高的概率，如果要求该项达到最大值，就需要求预测值和预测方差均很小，即该项要求寻找预测值小且预测准确的点；第2项为预测方差与概率密度函数的乘积，当预测的不确定性较大时该项较大，而概率密度函数又限制了预测值不可远离当前最优解.因此，最大化期望提高就是寻找具有以下特点的点：其预测值比当前最优响应值小，并且或者存在着较大的预测的不确定性.(4)收敛条件式中，Ymax和Ymin分别是样本点中最大和最小的响应值，fk和fk−1分别是通过ANSYS计算得到的第k步和k−1步的响应值，表示Kriging在第k步的预测值，ε1，ε2和ε3分别为给定的收敛精度.由于可降解聚合物血管支架优化设计中存在模型复杂、计算量大和高度的非线性等问题，本文采用改进的Kriging代理模型优化设计方法对ART18Z支架的关键几何尺寸进行优化设计.在该优化算法中，采用Kriging代理模型建立设计变量和设计目标之间近似的函数关系，运用优化拉丁超立方抽样方法获取最初样本点，并利用期望提高(EI)函数平衡局部和全局搜索，防止得到局部最优解.以支架的菱形孔的长和宽为设计变量，采用优化拉丁超立方抽样方法在设计空间内选取30个初始样本点，经过21步迭代后得到最优结果.表2对比了优化支架、原始设计支架和参照支架(在支架原始设计的基础上，仅将筋宽和厚度减小0.02mm，其他条件不变)的结果.与原始设计支架相比，优化支架与原始设计支架的径向回缩率差别很小，这说明即使优化支架的筋宽减小了15.38%，厚度减小了13.33%，但优化支架具有与原始设计支架相近的径向支撑能力；且优化支架筋宽和厚度较小，使得优化支架柔顺性优于原始支架，更利于支架在血管中的输送.此外，优化支架的表面覆盖率比原始设计支架的小，从而减小了新生内膜增生几率，降低了支架内再狭窄的风险；优化支架的血管管腔获得面积也比原始支架的大，血流流经血管更加通畅，有利于保证支架内血流通过量.总的来说，优化支架在明显减小筋宽和厚度的前提下，具有与原始设计支架相近的径向支撑能力，提高了支架的柔顺性，增大了血管管腔获得面积，降低了支架表面覆盖率，从而改善了支架的综合服役性能.与单纯减小支架筋宽和厚度的参照支架相比，尽管优化支架与其有相同的宽度和厚度，但是优化支架的径向回缩率比参照支架的小2.58%，这表明优化支架具有更好的径向支撑能力，而优化支架的血管管腔获得面积也比参照支架的大6.21%，更利于保证支架内血流通过量，且优化支架表面覆盖率与参照支架很相近.综上所述，虽然减小支架的筋宽和厚度有利于降低血管支架的表面覆盖率、增大管腔面积并提高支架的柔顺性，但单纯减小支架的筋宽和厚度，必然会降低支架的力学性能，尤其是径向支撑能力.然而，在减小支架筋宽和厚度的前提下，通过对支架结构进行优化设计，提高支架的力学性能，使其与筋宽和厚度较大的支架具有相似的力学特性，同时还能增大管腔面积，减小覆盖率，提高其综合服役性能，从而有利于降低再狭窄的风险.图6显示的是其他三个设计变量为最优值时，每个设计变量对径向弹性回缩的影响.支架的径向弹性回缩率随着支架筋宽和厚度的增大而减小，这是因为随着支架宽度和厚度的增大，支架的径向刚度也随之增大，从而使得支架的径向弹性回缩率减小.菱形孔的长度在设计范围内有一个最优值使得支架的径向弹性回缩最小，而菱形孔的宽度在设计范围有一个值使得支架的径向弹性回缩最大.这是因为研究的支架结构中，既有菱形也有六边形，而且他们共边.菱形的结构影响了六边形的几何结构，从而对支架的径向弹性回缩产生了影响.此外，支架的径向弹性回缩不仅与支架结构相关，而且与聚合物支架的材料和膨胀过程相关.本文针对可降解聚合物血管支架存在的径向刚度较弱，以及血管支架优化设计中存在的多重非线性耦合，传统的梯度类优化算法易落于局部最优，智能算法存在超大规模计算量的问题，提出了可降解聚合物血管支架的代理模型优化设计方法.并选取典型的ART18Z支架作为算例对其进行结构尺寸优化设计.经过优化后改善了该支架的综合服役性能.数值结果表明：该优化设计方法能有效地应用于可降解聚合物血管支架的结构优化设计.Kriging代理模型的优化算法具有较好的精确度和稳健性，对于多学科优化问题，尤其是对于解决支架优化这类模型复杂、计算量大、高度的非线性的问题具有很大优势.除了本文所提出的优化支架径向回缩的问题，此优化算法还能解决其他支架优化的问题，如改善支架的扩张性能、提高支架服役期的疲劳寿命、降低支架轴向缩短率等.Kriging代理模型的优化算法能有效地运用于血管支架的优化设计，建立设计目标和设计变量之间的近似函数关系，代替优化过程中费时的有限元计算，解决了冠脉支架设计中的设计目标与设计变量之间关系复杂且不能显式表达的问题，为血管支架的优化设计提供新的思路和方法.作为可降解聚合物血管支架优化设计的阶段性研究，本文在优化设计的目标函数中仅将径向弹性回缩作为设计目标，而将覆盖率作为约束来考虑，在后期工作中可在目标函数中适当考虑多个设计目标；此外，在优化过程中，对样本点的计算需要耗费大量计算量，在后期工作中可考虑血管支架优化设计的并行算法，以提高优化效率.【相关文献】1 Rivard A,Andr´es V.Vascular smooth muscle cell proliferation in the pathogenesis of atherosclerotic cardiovascular diseases.Histology and Histopathology,2000,15(2):557-5712 葛均波.支架内再狭窄 2014：当前问题.医心评论，2014(Ge Junbo.Intraventricular restenosis 2014:Current problem.Medical Heart Review,2014(in Chinese))3 葛均波,葛雷,黄榕翀.慢性完全闭塞病变介入治疗进展与展望.心血管病学进展,2007,28(2):165-167(Ge Junbo,Ge Lei,Huang Rongzhong.Advances in interventional therapy of chronic total occlusion.Advances in Cardiovascular Diseases,2007,28(2):165-167(in Chinese))4 高润霖.药物洗脱支架研究现状及进展.中国实用内科杂志,2006,26(15):1121-1123(Gao Runlin.Research status and development of drug-eluting stents.Chinese Journal of Practical Internal Medicine,2006,26(15):1121-1123(in Chinese))5 易勇,陈玉成,曾智.冠脉支架内再狭窄防治研究进展.心脏杂志,2003,15(6):560-562(YiYong,ChenYucheng,ZengZhi.Advances in prevention and treatment of coronary stent restenosis.Heart Magazine,2003,15(6):560-562(in Chinese))6 Sousa JE,Serruys PW,Costa MA.New frontiers in cardiology drugeluting stents:PartI.Circulation,2003,107(17):2274-22797 Jr DJ.Role of adjunct pharmacologic therapy in the era of drugelutingstents.Atherosclerosis Supplements,2005,6(4):47-528 刘赵淼,南斯琦,史艺.中等严重程度冠状动脉病变模型的血流动力学参数分析.力学学报,2015,47(6):1058-1064(Liu Zhaomiao,Nan Siqi,Shi Yi.Hemodynamic parameters analysis for coronary artery stenosis of intermediate severity model.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):1058-1064(in Chinese))9 Bresslo ffNW.Multi-Objective Design of a Biodegradable Coronary ArteryStent.Cardiovascular and Cardiac Therapeutic Devices:Springer Berlin Heidelberg,2013:1-2810 Zhang YJ,Bourantas CV,Farooq V,et al.Bioresorbable sca ff olds in the treatment of coronary artery disease.Medical Devices Evidence&Research,2013,6(1):37-4811 Kraak RP,Grundeken MJ,Koch KT,et al.Bioresorbable sca ff olds for the treatment of coronary artery disease:Current status and future perspective.Expert Review of Medical Devices,2014,11(5):467-48012 Pauck RG,Reddy putational analysis of the radial mechanical performance of PLLA coronary artery stents.Medical Engineering&Physics,2014,37(1):7-1213 Rogers C,Tseng DY,Squire JC,et al.Balloon-artery interactions during stent placement:A fi nite element analysis approach to pressure,compliance,and stent design as contributors to vascular injury.Circulation Research,1999,84(4):378-38314 Etave F,Finet G,Boivin M,et al.Mechanical properties of coronary stents determined by using fi nite element analysis.Journal of Biomechanics,2001,34(8):1065-107515 Migliavacca F,Petrini L,Colombo M,et al.Mechanical behavior of coronary stents investigated through the fi nite element method.Journal of Biomechanics,2002,35(6):803-81116 Mori K,Saito T.E ff ects of stent structure on stent fl exibility measurements.Annals of Biomedical Engineering,2005,33(6):733-742.17 Wu W,Yang DZ,Huang YY,et al.Topology optimization of a novel stentplatformwithdrugreservoirs.MedicalEngineering&Physics,2008,30(9):1177-118518 DeBM,VanIR,VerheggheB,etal.Finiteelementanalysisandstent design:Reduction of dogboning.Technology&Health Care O fficial Journal of the European Society for Engineering&Medicine,2006,14(4-5):233-24119 Conway C,Sharif F,Mcgarry JP,et al.A computational test-bed to assess coronary stent implantation mechanics using a populationspeci fi c approach.Cardiovascular Engineering and Technology,2012,3(4):374-38720 Wang WQ,Liang DK,Yang DZ,et al.Analysis of the transient expansion behavior and design optimization of coronary stents by finite element method.Journal of Biomechanics,2006,39(1):21-3221 王伟强,王丽,杨大智等.血管支架有限元优化设计.生物医学工程学杂志,2008,25(2):372-377(Wang Weiqiang,Wang Li,Yang Dazhi,et al.Optimization of vascular stent fi nite element method.Biomedical Engineering,2008,25(2):372-377(in Chinese))22 Atherton MA,Bates RA.Robust optimization of cardiovascular stents:A comparison of methods.Engineering Optimization,2003,36(2):1-1123 Harewood F,Thornton R,Ireland M.Step Change in Design:Exploring Sixty Stent Design Variations Overnight C.Altair Product Design Workshop,201124 Clune R,Kelliher D,Robinson JC,et al.NURBS modeling and structural shape optimization of cardiovascular stents.Structural&MultidisciplinaryOptimization,2014,50(1):1-1025 杜超凡,章定国.光滑节点插值法:计算固有频率下界值的新方法.力学学报,2015,47(5):839-847(Du Chaofan,Zhang Dingguo.Node-based smoothed point interpolation method:A new method for computing lower bound of natural frequency.Chinese Journal of Theoretical and Applied Mechanics,2015,47(5):839-847(in Chinese))26 杜超凡,章定国,洪嘉振.径向基点插值法在旋转柔性梁动力学中的应用.力学学报,2015,47(2):279-288(Du Chaofan,Zhang Dingguo,Hong Jiazhen.A meshfree method based on radial point interpolation method for the dynamic analysis of rotating fl exible beams.Chinese Journal of Theoretical and Applied Mechanics,2015,47(2):279-288(in Chinese))27 Li HX,Wang XY.Design optimization of balloon-expandable coronarystent.Structural&Multidisciplinary Optimization,2013,48(4):837-84728 Li HX,Zhang YH,Zhu B,et al.Drug release analysis and optimization for drug-eluting stents.Scienti fi c World Journal,2013,2013(4):827-83929 李红霞,王希诚.基于不同扩张模拟方式的支架优化设计.哈尔滨工业大学学报,2011(s1):267-272(Li Hongxia,Wang Xicheng.The stent optimization design based on di ff erent expansion simulation.Journal of Harbin Institute of Technology,2011(s1):267-272(in Chinese))30 Pant S,Bresslo ffNW,Limbert G.Geometry parameterization and multidisciplinary constrained optimization of coronary stents.Biomechanics and Modeling in Mechanobiology,2012,11(1):61-8231 Srinivas K,Nakayama T,Ohta M,et al.Studies on Design Optimization of Coronary Stents.Journal of Medical Devices,2008,2(1):121-13632 Chua SND,Macdonald BJ,Hashmi MSJ.Finite element simulation of slottedtube(stent)with the presence of plaque and artery by balloon expansion.Journal of Materials Processing Technology,2004,155-156(6):1772-177933 Joseph VR,Hung Y.Orthogonal-maximin latin hypercube designs.StatisticaSinica,2012,18(1):171-186。
CAE-心血管支架模拟
心血管支架模拟1介绍下图一为心血管支架作用原理示意图，心血管支架用以张开硬化的动脉。
图心血管支架作用原理示意图支架的成功在很大程度上取决于如何将支架在动脉中支撑起来。
在支架设计过程和临床手术过程，利用有限元分析技术作为工具进行支架动脉相互作用研究已经成为一个公认的趋势。
在支架支撑的过程中涉及到非常多的非线性技术，如生物组织特性、大动脉壁变形、支架和动脉壁的滑动接触。
例题中均会涉及到这些。
2问题描述建立的冠状动脉支架和严重阻塞的冠状动脉模型。
动脉被简化为两层的直圆柱体，用一层代表动脉壁和另一个代表钙化斑块。
下图为动脉和支架的尺寸。
图动脉与支架尺寸模拟支架置入术分为三个步骤：➢扩大动脉升高压力（球囊血管成形术）。
➢放置的支架。
➢收缩动脉，使用平均血压和使得支架和动脉壁接触。
3模型3.1支架模型采用Beam189单元建立支架模型并划分了1760个网格，如图所示。
图心血管支架模型梁单元模型的计算效率优于实体单元。
支架为3.5毫米的直径，长度15毫米。
构建支架的金属丝是圆形截面，外直径为0.1毫米。
镍钛合金材料是常用作支架材料，镍钛合金的非线性材料行为需要单独讨论。
因此对于这个问题，材料采用线弹性的316L不锈钢材料来代替。
3.2动脉和斑块模型采用实体单元模拟动脉血管和斑块，如图所示：图动脉血管和斑块（钙化层）有限元模型动脉层共划分了9000个单元，采用的是SOLID185单元，采用增强应变公式（KEYOPT（2）=3）和混合U-P公式（KEYOPT（6）=1），这是由于生物材料为不可压缩材料，采用上述公式可以克服体积锁定。
钙化层采用SOLID185单元，有9000个网格。
单元采用全积分算法，材料采用线性材料。
动脉和钙化层之间采用共节点的方式连接。
基于圣维南原理，动脉和钙化层都作了延长，减少边界对计算结果的影响。
接近边界的地方采用密的网格，增强收敛性。
3.3支架和斑块的接触模型支架和动脉之间的接触采用线面接触类型，支架采用CONTAC177接触单元，接触方法采用拉格朗日乘子法，接触行为采用标准接触行为，打开步长二分，动脉内部面采用TARGE170单元，摩擦系数为0。
有限元圆角L形支架分析及对本专业作用
有限元圆角L形支架分析及对本专业作用支架是什么呢？支架是一种微小的网状装置，放置在动脉、血管或其他管道内，保持结构扩张以疏通血管。
支架用于哪里？通常用于治疗由于体内血管阻塞或损伤引发的
疾病，包括冠状动脉疾病、外周动脉疾病、颈动脉疾病、肾血管性高血压，以及腹主动脉瘤等。
其他使用支架的还有保持扩张阻塞或受损的输尿管、保持胆汁在阻塞的胆管内流动、帮助气道阻塞的病人呼吸等。
支架是如何工作？收缩的支架通过一个小切口穿入，安装到受影响的导管内。
支架扩张使导管变宽，并提供支撑以防止导管变窄。
通常用的支架有两大类型：1）气球胀开支架：其由由塑性变形的材料制成，通过气球胀开进而扩张，气球放气后，除了稍微回弹，支架仍保持扩张状态。
2）自展式支架：由具有较大弹性应变的材料制成，消除约束传递系统后，支架通过弹性变形恢复至扩张状态。
支架使用的材料：
1）不锈钢广泛应用于气球胀开支架；气球胀开支架的替代材料有钽、铂合金、铌合金、钴合金等；
2）镍钛诺超弹性合金用于自展式支架；
3）较为前沿的新型支架材料有形状记忆聚合物，生物可降解聚合物和生物可降解金属。
支架的制造加工工艺：
1）大多数冠状动脉支架是由激光切割制造，激光束切割出复杂
的设计形态；大多数支架是从金属管切割而来，有些是由金属板切割，之后卷成管状；
2）其他制造方法包括卷绕、编织、针织、水切割和光化学蚀刻等。
支架的几何设计有连续环、单独环和编织等。
其中，占据大部分市场的连续环是由一系列可扩展的Z形或S形单元组成；单独环是由单个Z形或S形环构成，用于支撑移植物；而编织设计，具有卓越的覆盖性能，但在扩张期间会大幅度缩短。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
胀，即两端先翘起的“狗骨头”状膨胀现象，是造成支架端部急性动脉损伤的主要原因之一，如何最大程度地减少支架瞬时膨胀过程中的“狗骨头”现象是减少这种损伤的关键因素。
因此，本研究预利用有限元技术来对血管支架瞬时膨胀过程中的力学行为进行分析，以期通过支架结构上的优化设计达到减少或消除其“狗骨头”膨胀现象的目的。
２模型与方法
２．１有限元模型
２．１．１几何模型球囊伎架系统在血管内的安放过程是一个十分复杂的生物力学过程，再加上对斑块、动脉等这些软组织以及它们之间的交互作用的力学模拟通常具有高度的非线性，因此，为了简化计算，本研究只考虑了球囊与支架这两个组成部分以及它们之间的相互作用。
所模拟的支架原型采用专有支架设计，支架的原始外径为１．５６ｍｍ，支柱横截面近似为０．１瑚ｍ×０．１ｍｍ的矩形，支架的原始长度为１１．９ｍｍ。
这里采用比支架在每端长出来Ｏ．５ｍｍ的球囊与上述支架匹配，组成球囊伎架输
送系统，图１为其组合图，根据模型在长轴方向上的对称性，图中显示的组合模型为总模型的１／２。
球囊伎架系统的三维几何模型先由ＰＩ｝０厄ｎｇｉｎｅｅｒ软件包建立，然后输入到如ｙｓ通用有限元软件包进行模拟分析。
２．１．２材料属性支架材料选用３１６Ｌ医用不锈钢，取其弹性模量为２０１ＧＰａ，泊松比为０．３。
因支架的几何大变形主要靠材料的局部塑性变形提供，材料特性采用ｖｏｎＭｉｓｅＳ屈服准则和各向同性多线性强化准则。
球囊在支架的大部分变形阶段中所受的分力很小，如果模拟真实的、折叠的球囊，必须进行动态模拟，且计算机性能要求很高，同时还要考虑球囊自身的接触等因素。
所以为简化计算而又不影响支架变形效果，用一厚为０．１８ｍｍ的圆柱壳虚拟材料模型来模拟球囊。
此处把球囊作为各向同性线弹性材料，采用弹性模量为１０ＭＰａ。
同时考虑其不可压缩性，泊松比选取为０．４９９９。
２．１．３网格划分和约束设置根据支架的结构特点，选用２０节点ｓｏｌｉｄｌ８６单元进行网格划分。
针对球囊的壳状特征，选用４节点Ｓｈｅｌｌｌ８１单元进行网格划分［８Ｊ。
综合考虑支架、球囊各自的求解精度和整体求解效率的因素，球囊网格划分得相对粗糙，支架则相对精细。
最终网格划分结果如图２所示，支架筋的厚度和宽度方向上都至少有２层单元。
因球囊和支架问摩擦的参数未知，同时也为了进一步简化模型，在球囊和支架间建立了无摩擦接触模型。
以柱坐标为基础建立了球囊伎架系统的边界约束条件。
因为支架绕轴１２０。
旋转对称，所以。
实际计算中取球囊伎架在周向上的１／３作为研究对象。
对于球囊，为了保证其只发生径向上的位移，约束其两端在轴向上的位移，约束其两侧节点在周向上的位移，并偏移耦合其两侧节点的剩余自由度。
对于支架，约束其中间对称面在轴向上位移，约束其两侧循环对称面上在周向上的位移，循环对称面上其他两个方向的位移进行偏移耦合，支架端部则不受约束。
这样即保证了支架在模拟过程中不会产生刚性移动，同时又不阻碍其在径向和轴向上的自由变形。
圈１球囊／支架组合模型图
ＦｉｇｌＢａｌｌ０吼／如ｎｔ舾嘲ＩｌｂｌｙｌＩＨ埘
图２球囊伎架模型网格划分示例图
ｎｇ２
ｎｍ叫眵眦ｉ佃ｏｆｌＴ翳ｈｅｄｈ¨鲫ｌ／§ｔｅＩｎａｓｓ哪ｂｌｙ
２．１．４求解控制球囊伎架系统的膨胀，靠给球囊内壁加压来实现，所加压力值为１．０Ⅷａ（约１０ａｔｍ）。
系统从最初原始状态开始向外膨胀。
因为球囊被当作线弹性圆柱壳材料来处理，在实际模拟分析中发现，对于单纯的球囊壳体，其在持续的受压膨胀过程中。
会不可避免地出现膨胀失稳的状态，亦即在球囊受压膨胀时，会出现直径持续增加，而压力不增加，反而减小的情况，这在数值求解上将会表现出刚度矩阵的奇异，对于这种情况，如果仍然采用常规的Ｎｅ叭ｏｎ—Ｒａｐｈｓｏｎ迭代数值求解方法会导致严
重的收敛问题。
为此，此处采用弧长法进行静力求
图４球囊／支架模型在内压力作用下．量大膨胀直径达到２．８硼左右时的膨胀结果。
彩色云图代表着球囊伎架系统在径向上位移的变化
（ａ）模型Ａ；（ｂ）模型Ｂ；（ｃ）模型Ｃ；（ｄ）模型Ｄ
Ｆｌｇ４Ｓｉｍｍａｔｉ蛐ｒ荡ｌｌＩｔｓｏｆｔｌ埒舟ｄｉｍｍｓｐｌ粼科咖协ｏｆｆ幽丌ｎ－０ｄｅｋ．ＴＩ埒ｓｉｍｕｈｔｊ硼ｏｆｎ蟹∞驰璐ｉ蛐ｏｆｂ柚咖／ｓｔｅＩＩｔ掣ｓ咖惴ｋｎＩｌｉ珊ｔｅｄｗｈｍｔｌ虻眦ｘｉｍ唧出卸—ｅ协ｒｏｆ２．８删ｎｗ罄ｒ目Ｋ捌
（ａ）ｍｏｄｅｌＡ；（ｂ）ｍｏｄｅｌＢ；（ｃ）ｍｏｄｅｌＣ；（ｄ）ｍ３ｄｅｌＤ
作者早期的一篇关于支架防“狗骨头”设计的文章中曾指出，球囊伎架系统在瞬时膨胀过程中的“狗骨头”指标不仅与支架端部的结构强度相关，而且还与其球囊的过长量紧密相关。
通过减少球囊的过长量和增加支架端部网丝的宽度均可以减弱支架的“狗骨头”状膨胀现象，但是只有两种方法共同应用才能起到更加明显的效果【９Ｊ。
现在看来，通过增加支架端部支撑体筋的宽度这一途径来减少支架“狗骨头”现象，有可能会对支架与血管的柔性匹配（Ｃ０ｎｆｏｒｎｎａｂｉｌｉｔｙ，指支架撑开后顺应血管原始形状的能力，它和支架撑开后的轴向弯曲性能成正相关）产生影响，因为过度的增加支架端部支撑体筋的宽度，一方面会使得撑开后的支架端部弯曲刚性增大，不能很好地匹配血管柔性波动，影响到血液动力学，另一方面有可能增加支架端部血管壁的应力集中，促使血管内膜过度增生，从而引发后期的支架内再狭窄［１０】。
而采用本文的这种由中间到端部逐渐地减小环状支撑体宽度的方法来进行支架防“狗骨头”的设计，较单纯的增加支架端部网丝尺寸的方法相比，一方面可以增加支架瞬时膨胀的均匀性，另一方面则不会过多地影响支架端部的弯曲刚性，也就是说，能够更好地与血管进行柔性匹配，从而降低支架端部的再狭窄几率，不过，这一结论还有待进一步的流体动力学和相应的实验验证。
４结论
（１）在支架的结构参数中环状支撑体的宽度对
支架膨胀压的影响非常显著，环状支撑体宽度小的结构所需要的临界以及后续膨胀压比较高。
在支架使用直径确定的情况下。
对于这种典型的“Ｚ”字形或正弦波形环状支撑体结构的设计，应避免选取过窄的参数值。
（２）通过本文有限元模拟结果可知，在设计血管支架时，其环状支撑体结构的非对称设计，即环状支
撑体的宽度从中部到端部依次逐渐递减变化，可以有效抑制了球囊伎架系统在瞬时膨胀过程中的“狗骨头”现象。
从而有可能避免单纯靠增加支架端部筋宽来设计防“狗骨头”支架时所带来的弊端。
（３）本文关于血管支架的有限元分析结果，可以为其它类似血管支架的结构优化设计提供有益的帮助。
参考文献
ｌ“代，ＨｕＤＹ．ＴｈｅｐＩＤｇ嗍ｏｆｉｎ眦岫ａｒｔｅｒｙ吼饥ｔ．伽腿ＭｅｄｉｃａＩＤ酣略Ｉｎｆ姗ｔｉ∞，２０００；６（２）：６［李田昌，胡大一．
冠状动脉内支架进展．中国医疗器械信息．２０００；６（２）：６］
２ＥｒｌｄＲ，地ｕｄｅＭ，ｍｐｐＨＷ，矗４ｆ＿Ｃｏ∞ｎａ驴矾ｅｒｙ吼朗ｔ啦ｏｏｍｐａ捌ｗｉｔｈｂａｌｌ饯ｍａｎｇ；ｏｐｌａｓｃｙｆｏｆｒｅｓｔ唰５ａｆｔ口ｉｌｌｉｔ柚ｂａｌｌ０。
ｒＩ锄ｇｉｏｐＩａｓ哆．Ｒ吲舶曲ＳｔａｎｔＳｃｕｄｙＧ∞印．ＮＥｎｇｌＪ
Ｍｅｄ，１９９８；。