大動脈弁FSI解析
理论与物理
主动脉瓣随着心脏每次收缩而开闭，每天承受约10万次载荷循环。瓣叶的开闭行为由血流与瓣膜组织的强相互作用决定，因此仅进行流体或仅进行结构分析是不够的。这是人工瓣膜（机械瓣、生物瓣、TAVI瓣）设计优化、瓣膜病进展预测、手术规划中不可或缺的技术。
瓣膜开闭时的力学是怎样的呢？
收缩期左心室压力超过主动脉压力时瓣膜打开，血液射出。舒张期在反向压力梯度下瓣膜关闭。瓣叶是厚度约0.5mm的薄组织，反复经历大变形。瓣膜后流在Valsalva窦内形成涡流，达芬奇早已观察到这个涡流有助于瓣膜关闭。
控制方程
瓣膜的FSI需要解哪些方程呢？
流体侧是不可压缩Navier-Stokes方程。收缩期峰值雷诺数 $Re \approx 5000$〜$8000$，有时需要考虑湍流转捩。
结构侧将瓣叶建模为超弹性壳。常使用考虑纤维增强的Fung型或Lee-Sack型本构关系。
其中 $E_{11}$ 是纤维方向（周向），$E_{22}$ 是正交方向的Green-Lagrange应变。
瓣叶的接触是如何处理的呢？
瓣膜关闭时三片瓣叶相互接触（对合）。这是主动脉瓣FSI最大的技术难点。结构侧需要处理接触力学，流体侧需要处理间隙趋近于零的极限情况。浸入边界法的优点在于无需改变流体网格拓扑即可处理接触。
主动脉瓣每天承受10万次开闭
主动脉瓣的小叶（瓣叶）是厚度仅约0.5mm的组织，在收缩期承受约120mmHg（16kPa）的压力差，每天开闭10万次。设计寿命约30年。支撑这种惊人耐久性的是胶原纤维交叉排列的各向异性结构。要在FSI理论中再现这种纤维结构，简单的各向同性弹性体是不够的，需要Fung型超弹性模型等。理论的复杂性，正是生物组织卓越性的体现。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\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%以下为目标 |
数值解法与实现
主要有三种方法。
| 方法 | 优点 | 缺点 |
|---|---|---|
| ALE-FEM | 界面精度高 | 瓣膜关闭时网格崩溃 |
| Immersed Boundary (IB) | 对大变形·接触处理能力强 | 界面处存在涂抹效应 |
| Immersogeometric | IGA精度 + IB灵活性 | 实现复杂 |
Immersogeometric听起来很陌生？
这是Kamensky, Hsu, Bazilevs (2015)提出的方法，将基于NURBS的瓣叶模型嵌入固定的流体网格。它通过Nitsche法的界面条件克服了IB法的弱点（delta函数的涂抹效应）。这是德克萨斯大学奥斯汀分校Bazilevs研究室引领的领域。
IBAMR/IBFE的实现
具体的软件是？
IBAMR（Immersed Boundary Adaptive Mesh Refinement）是开源软件的代表。由Griffith教授（北卡罗来纳大学）开发，使用IBFE（Immersed Boundary Finite Element）法将瓣叶的FEM结构嵌入基于SAMRAI自适应网格的流体中。
流体使用带惩罚项的IB法求解Navier-Stokes方程。
瓣叶的弹性力 $\mathbf{F}$ 由FEM计算，通过delta函数分散到流体网格上。
时间步长与CFL条件
需要多大的时间步长呢？
瓣膜开闭大约发生在30ms内（从收缩期开始到完全开放）。根据CFL条件 $\Delta t \leq h/|\mathbf{u}_{max}|$，峰值流速1.5 m/s，最小网格宽度0.1mm时 $\Delta t \leq 67\mu$s。实际计算中采用 $\Delta t = 10$〜$50\mu$s。一次心跳（0.8s）的模拟需要16,000〜80,000步。
计算成本真庞大啊。
是的。必须使用AMR（自适应网格细化）仅在瓣膜附近进行细化。在IBAMR中，将瓣膜表面数毫米内设为最细网格，远处设为粗网格，可将计算量降至1/10以下。
IBM与ALE的选用——主动脉瓣分析的关键判断点
在主动脉瓣FSI的数值解法中，常被讨论的选择是“浸入边界法（IBM）还是ALE法”。处理瓣膜大变形·接触时IBM有利，但边界附近精度易下降。ALE法精度高，但需要应对完全关闭时的网格畸变。在实际应用中，有研究小组采取“开放期用ALE，关闭期切换为IBM”的混合策略。选择哪种方法，代码的复杂性会有很大不同。
整体式方法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
独立求解各物理场，在界面交换数据。实现容易，可活用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（守恒性好）、RBF插值（对网格不一致鲁棒性强）。需要在守恒性与精度间取得平衡。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的一致性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止过松弛导致发散，是加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，反弹又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，根据前一次的修正量自动调整下一次修正量的自适应方法。
实践指南
典型的流程是这样的。
1. 瓣膜形状模型创建: 从参数化形状（瓣叶高度、对合高度、瓣环直径）创建CAD，或从心脏超声/CT重建患者特异性形状
2. 主动脉窦（Valsalva窦）模型: 构建从瓣环到升主动脉的流体区域
3. 网格生成: 瓣叶（壳或实体）与流体区域。IB法下瓣膜网格独立于流体网格
4. 材料参数设定: 瓣叶的超弹性常数。根据双轴拉伸试验数据标定
5. 边界条件: 入口施加左心室压力波形，出口施加主动脉压力波形或三元件Windkessel模型
6. 计算执行: 3次心跳以上（消除初始瞬态）
7. 后处理: 瓣口面积（EOA）、压降、反流量、瓣叶应力
EOA（有效瓣口面积）是什么？
是Effective Orifice Area的缩写，作为人工瓣膜最重要的性能指标。由Gorlin公式定义。
$Q_{rms}$ 是RMS流量（mL/s），$\Delta p_{mean}$ 是平均压降（mmHg）。对于TAVI瓣膜，EOA > 1.0 cm² 被认为是良好的。
材料参数标定
瓣叶的材料常数如何确定？
将双轴拉伸试验数据拟合到Fung型或Lee-Sacks模型。天然瓣膜典型值为 $c_0 = 2$〜$10$ kPa，纤维方向的刚度参数 $c_1$ 通常是正交方向 $c_2$ 的5〜20倍。
| 瓣膜类型 | $c_0$ (kPa) | $c_1$ | $c_2$ | 出处 |
|---|---|---|---|---|
| 天然主动脉瓣 | 2〜10 | 10〜50 | 1〜5 | Billiar & Sacks (2000) |
| 牛心包（TAVI瓣膜） | 5〜20 | 30〜80 | 5〜15 | 因产品而异 |
| 猪瓣膜（生物瓣） | 3〜15 | 15〜60 | 2〜10 | Stella & Sacks (2007) |
钙化的瓣膜如何处理？
主动脉瓣狭窄症中瓣叶会发生钙化。通常采用从CT中提取钙化区域，并将该区域的弹性模量设为10〜100倍的方法。钙化程度直接关系到TAVI瓣膜植入后的瓣周反流（PVL），因此精确建模非常重要。
从CT图像重建瓣膜形状的实际困难
在主动脉瓣FSI分析的实际工作中，最耗时的是“患者特异性几何形状重建”。从CT或MRI图像中分割出瓣叶形状是项艰巨任务，即使经验丰富的研究者，一个病例也可能花费数天。近年来出现了基于深度学习的自动分割技术，时间缩短至数小时。即便如此，仍有人说“分析结果的好坏，80%取决于几何形状的质量”。
分析流程的比喻
你吹过气球吗？那个瞬间，其实发生了高级的流体-结构耦合。内部空气压力（流体）推动橡胶壁（结构）扩张→扩张的壁改变内部压力分布→改变的压力进一步使壁变形……这种“投接球”在每个计算步中反复进行，就是FSI分析。
相关主题
なった
詳しく
報告