血管支架展开流体结构相互作用
概述
先生、血管ステントの展開シミュレーションでFSIが必要な理由は何ですか？
理论与物理
支架是一种在血管内扩张金属网状结构以拓宽狭窄血管的医疗器械。支架展开后的血流模式强烈依赖于支架形状和血管壁变形，并直接关系到再狭窄（neointimal hyperplasia）的风险。已知在WSS（壁面剪切应力）较低的区域更容易发生再狭窄，因此通过FSI准确预测血流场和壁面应力非常重要。
支架有哪些类型？
主要有三种类型。球囊扩张型（BES：冠状动脉支架）、自膨型（SES：镍钛合金，颈动脉/外周血管用）、药物洗脱型（DES：涂覆再狭窄抑制药物）。它们的力学行为各不相同，因此结构模型也会有所不同。
控制方程
支架展开的FSI需要求解什么样的方程组？
这是一个三个结构体（支架、球囊、血管壁）与血流耦合的问题。
支架使用SUS316L或镍钛合金（超弹性）建模。镍钛合金的超弹性通常采用Auricchio-Taylor模型作为标准。
相变应变 $\boldsymbol{\varepsilon}^{tr}$ 在奥氏体-马氏体相变时产生。
血管壁采用Holzapfel-Gasser-Ogden模型，血流采用不可压缩Navier-Stokes方程。存在斑块时，使用线弹性或塑性模型来表现硬质层。
球囊是如何建模的？
球囊由聚酰胺等聚合物薄膜制成，使用膜单元（membrane element）建模。通过逐步增加内压，将支架从压握直径扩张到目标直径。有时也会再现球囊的折叠形状，但这需要在计算成本之间权衡。
镍钛合金——"形状记忆"引发的支架革命
现代血管支架中广泛使用的镍钛合金（Ni-Ti合金）具有形状记忆效应和超弹性两大特性。在体温（37℃）附近转变为奥氏体相，使以收缩状态插入的支架在血管内自然展开。要在FSI理论中处理镍钛合金，需要超弹性本构关系（如Brinson法则等），而非通常的线弹性体。这种材料模型的实现，是提升支架FSI分析难度的最大因素之一。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常示例】夏天铁轨伸长导致间隙变窄——温度升高→热膨胀→产生应力的典型例子。电子基板在焊接后翘曲也是由于不同材料热膨胀率差异所致。发动机气缸体因高温部和低温部的温差产生热应力，最坏情况下会导致裂纹。
- 流体-结构耦合（FSI）项：流体压力/剪切力使结构变形，结构变形又改变流体区域，是双向相互作用。【日常示例】强风导致悬索桥缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风流，进而放大振动。心脏血流与血管壁的弹性变形、飞机机翼的颤振（气动弹性不稳定性）也是典型的FSI问题。有时仅需单向耦合即可，但变形较大时双向耦合是必须的。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度升高，温度变化改变电阻，形成反馈回路。【日常示例】电暖炉的镍铬丝通电发热（焦耳热）变红——温度升高导致电阻变化，电流分布也随之改变。IH电磁炉的涡流发热、输电线路温度升高导致的垂度增加也是这种耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常示例】天气预报中结合"气温数据"和"风力数据"计算体感温度时，若观测点不同就需要插值——CAE的耦合分析中，结构网格和CFD网格通常也不一致，因此界面处的数据传递（插值）精度直接关系到结果的可靠性。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性较强时
- 时间尺度分离：各物理场特征时间差异较大时，可通过子循环提高效率
- 界面条件一致性：需确认耦合界面处能量/动量守恒在数值上得到满足
- 不适用的情形：三个以上物理场同时强耦合时，有时可能需要整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢：约12×10⁻⁶，铝：约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或N（集中力） | 确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度依赖于网格密度比。5%以下为目标 |
数值解法与实现
通常分为两个阶段。
第1阶段: 支架展开（仅结构）
1. 支架的压握（缩径）
2. 球囊加压扩张
3. 球囊减压后的回弹（回弹）
第2阶段: 展开后的FSI
1. 将展开后的变形形状设为初始条件
2. 导入残余应力
3. 执行脉动血流的FSI分析
导入残余应力很重要呢。
展开后的支架存在较大的残余应力，忽略此应力会大幅降低壁面应力的预测精度。可以使用Abaqus的*IMPORT功能或Ansys的ICTRL/RESUME命令来继承前次分析的应力场。
接触算法
支架与血管壁的接触是如何处理的？
存在三个接触对：支架-球囊、支架-血管壁、球囊-血管壁。
支架与血管壁的接触是如何处理的？
存在三个接触对：支架-球囊、支架-血管壁、球囊-血管壁。
| 接触对 | 方法 | 摩擦系数 |
|---|---|---|
| 支架-球囊 | 面-面接触 | 0.05〜0.1 |
| 支架-血管壁 | 面-面接触 | 0.1〜0.2 |
| 支架-斑块 | 面-面接触 | 0.2〜0.3 |
Abaqus中使用*CONTACT PAIR或General Contact，LS-DYNA中使用*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE。增广拉格朗日法比罚函数法更能抑制穿透。
网格策略
支架的细长筋杆网格划分看起来很难。
筋杆截面（典型尺寸80〜120μm×80〜120μm）用六面体单元分割为3×3以上。整个支架将达到50万〜200万个单元。
血管壁在支架筋杆附近进行局部细化，单元尺寸设为筋杆宽度的1/3以下。流体网格需要在筋杆间隙也布置单元，整体通常达到100万〜500万个单元。
那种规模的计算时间大概要多久？
展开分析（准静态）通常需要8〜24小时，FSI（数个心动周期）需要24〜72小时。这是基于隐式解法、16〜32核并行计算的典型预估时间。
支架展开的两步分析——分开求解折叠与展开
血管支架展开的FSI分析，通常分为"折叠工序的结构分析"→"血管内展开的FSI分析"两步进行。关键点在于将折叠时产生的残余应力作为展开分析的初始条件继承下来。如果省略这一步，直接从展开后的支架形状开始分析，实际展开力及对血管的负荷可能会被低估20〜30%。这种两步法在医疗器械的注册申报资料中也逐渐被要求提供。
整体式法
将所有物理场作为一个联立方程组同时求解。对于强耦合问题稳定，但实现复杂且内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。需要在守恒性和精度之间取得平衡。
子迭代
在每个耦合步内进行充分的迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止过度松弛导致发散，是一种加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于"平衡跷跷板"。一方推得太用力，另一方就会弹起，反作用力又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，这是一种根据上次修正量自动调整下次修正量的自适应方法。
实践指南
1. 支架形状创建: 从CAD中创建展开形状的筋杆图案，并卷绕成圆柱体
2. 血管形状获取: 从CT/MRI获取患者特异性形状，或使用IDEAL形状（直管+狭窄部）
3. 材料定义: 支架（SUS316L: 弹塑性 / 镍钛合金: Auricchio超弹性）、血管壁（Holzapfel-Gasser-Ogden）、斑块（弹塑性）
4. 压握〜展开分析: Abaqus/Explicit 或 Abaqus/Standard
5. 展开后形状的FSI分析: Ansys System Coupling 或 Abaqus co-simulation
6. 后处理: 支架的von Mises应力（疲劳评估）、WSS、OSI（再狭窄风险评估）
也要做疲劳评估吗？
FDA（美国食品药品监督管理局）要求对支架进行疲劳寿命评估。需要满足4亿次循环（相当于10年）的疲劳耐久性要求。在Goodman图上绘制平均应力和交变应力，确认其低于疲劳极限。
安全系数SF ≥ 2 是普遍标准。
常见陷阱
实际工作中容易陷入的难点在哪里？
| 问题 | 原因 | 对策 |
|---|---|---|
| 支架未与血管壁贴合 | 回弹被高估 | 设置球囊过度扩张率（10〜20%） |
| 筋杆断裂判定 | 网格依赖的应力集中 | 使用子模型法进行局部精细化 |
| 不现实的狗骨现象 | 球囊加压顺序 | 采用中央先行加压的分步加压 |
| WSS分布左右不对称 | 血管弯曲、分叉的影响 | 入口条件设置足够长的入口段 |
狗骨现象是什么？
球囊扩张型支架中，两端先扩张而中央延迟扩张的现象。这会导致支架端部对血管壁施加过大应力，成为再狭窄的原因。需要通过球囊顺应性和支架的单元设计来控
なった
詳しく
報告