主动脉瓣FSI解析
主动脉瓣FSI的理论基础
主动脉瓣随心脏收缩而开闭，每天承受约10万次负荷循环。瓣叶的开闭行为由血流与瓣膜组织的强相互作用决定，仅用流体或结构解析是不够的。是人工瓣膜（机械瓣、生体瓣、TAVI瓣）设计优化、瓣膜病进展预测、手术计划的必需技术。
瓣膜开闭时的力学是怎样的呢？
收缩期左室压超过主动脉压时瓣膜开启，血液被驱出。舒张期由于逆压力梯度，瓣膜关闭。瓣叶是厚度约0.5mm的薄组织，反复承受大变形。瓣膜后流的Valsalva窦内形成涡流，这个涡流辅助瓣膜关闭，这是达·芬奇早就观察到的现象。
支配方程
瓣膜的FSI解析中求解什么方程呢？
流体侧是不可压缩Navier-Stokes方程。收缩期峰值处雷诺数约为 $Re \approx 5000$〜$8000$，应考虑湍流转变。
结构侧将瓣叶建模为超弹性壳。采用考虑纤维增强的Fung型或Lee-Sacks型本构关系。
其中 $E_{11}$ 是纤维方向（周向）的Green-Lagrange应变，$E_{22}$ 是垂直方向的应变。
瓣叶的接触怎样处理呢？
瓣膜关闭时3片瓣叶相互接触（coaptation）。这是主动脉瓣FSI的最大技术难点。结构侧需处理接触力学，流体侧需处理间隙趋于零的极限。Immersed Boundary法的优势在于无需改变流体网格拓扑就能处理接触。
主动脉瓣每天10万次开闭承受住了
主动脉瓣的小叶（瓣叶）厚度仅约0.5mm，在收缩期承受约120mmHg（16kPa）的压力差，每天开闭10万次。设计寿命约30年。支撑这种惊人耐久性的是胶原蛋白纤维交叉排列的各向异性结构。FSI理论中要再现这种纤维结构，简单的等向性弹性体是不够的，必须用Fung型超弹性模型等。理论的复杂性反映了生体组织的精妙之处。
主动脉瓣FSI的数值计算手法
主要有3种方法。
| 手法 | 优点 | 缺点 |
|---|---|---|
| ALE-FEM | 界面精度高 | 瓣膜关闭时网格破坏 |
| Immersed Boundary (IB) | 大变形·接触强 | 界面处产生模糊 |
| Immersogeometric | IGA精度+IB灵活性 | 实现复杂 |
Immersogeometric是什么呢？
Kamensky, Hsu, Bazilevs (2015)提出的手法，将基于NURBS的瓣叶模型嵌入固定流体网格中。用Nitsche法处理界面条件来克服IB法的弱点（δ函数的模糊）。是由德克萨斯大学奥斯汀分校Bazilevs研究室主导的研究领域。
IBAMR/IBFE的实现
具体的软件是什么呢？
IBAMR（Immersed Boundary Adaptive Mesh Refinement）是开源的代表。由北卡罗来纳大学Griffith教授开发，用IBFE（Immersed Boundary Finite Element）法在SAMRAI自适应网格上将瓣叶的FEM结构嵌入流体中。
流体用带罚函数法的IB法求解Navier-Stokes方程。
瓣叶的弹性力 $\mathbf{F}$ 由FEM计算，用δ函数分散到流体网格。
时间步长与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。1个心跳（0.8s）的仿真需要16,000〜80,000步。
计算成本太大了。
是啊。因此必须用AMR（自适应网格细分）。IB AM R用瓣膜周围数mm内的最细网格，远处用粗网格，能把计算量抑制在1/10以下。
IBM与ALE的选择——主动脉瓣解析的判断要点
主动脉瓣FSI的数值方法中，"隐式边界法（IBM）还是ALE法"的选择一直是讨论的焦点。瓣膜的大变形·接触处理用IBM有利，但边界附近的精度易下降。ALE法精度高，但完全关闭时需要对网格粉碎进行处理。实际工作中，"开口期用ALE，关闭期切换到IBM"的混合策略有研究小组采用。选择哪一种会很大程度上影响代码的复杂度。
主动脉瓣FSI的实务应用
典型的流程是这样的。
1. 瓣膜形状模型创建: 从参数化形状（瓣叶高度、对合高度、瓣膜直径）创建CAD，或从心脏超声/CT重建患者特异性形状
2. 主动脉窦（Valsalva窦）模型: 从瓣膜环到升主动脉的流体区域构建
3. 网格生成: 瓣叶（壳或实体）与流体区域。IB法中瓣叶网格与流体网格独立
4. 材料参数设定: 瓣叶的超弹性常数。从双轴拉伸试验数据同定
5. 边界条件: 入口施加左室压力波形，出口施加主动脉压力波形或3单元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图像中分割出瓣膜小叶形状是难题，经验丰富的研究人员处理1个病例也要花数天。近年来深度学习自动分割出现，时间缩短到数小时。但即便如此，"解析结果的好坏8成取决于几何质量"这一论断在该领域仍成立。
主动脉瓣FSI的软件比较
FEBio初次听说。
TAVI瓣的留置模拟
TAVI瓣的设计中怎样利用FSI呢？
TAVI（经导管主动脉瓣置换）的留置过程分阶段仿真。
1. 支架的压缩: 用Abaqus对自膨胀型/球囊扩张型支架压缩
2. 导管运输: 支架在主动脉弓内的变形
3. 瓣膜展开: 支架在钙化天然瓣膜内的扩张
4. FSI评估: 展开后瓣叶的开闭行为、EOA、PVL（瓣膜周围逆流）的预测
Medtronic和Edwards Lifesciences也用FSI吗？
FDA提交资料中有越来越多包含FSI结果的案例。但企业内的具体手法往往未公开。从学术层面，Dasi教授（乔治亚理工）和Sotiropoulos教授（明尼苏达大学）的小组发表了领先研究。
SimVascular和OpenFOAM——开源力量开创心脏瓣膜研究
主动脉瓣FSI解析中商用工具以ANSYS Fluent和Abaqus出名，但研究界OpenFOAM的SimVascular（斯坦福大学发）正迅速普及。从患者特异性血管形状的导入、边界条件设置到一整套流程均受支持，免费可用。与商用工具的验证论文也在增加，"预算受限的大学医院也能进行本格瓣膜解析"的环境正在形成。
主动脉瓣FSI的前沿研究
收缩期峰值流速达1〜1.5 m/s（狭窄瓣超过4 m/s），瓣后流呈湍流转变。雷诺数 $Re = 5000$〜$10000$ 时层流假设不适用。
LES（Large Eddy Simulation）或DNS（Direct Numerical Simulation）理想，但壁面分解LES需要 $\Delta x \sim 10\mu$m，计算量巨大。Wall-modeled LES或DES（Detached Eddy Simulation）是现实选择。
湍流对瓣膜的影响是什么呢？
湍流剪切应力（turbulent shear stress）参与瓣叶损伤和红细胞溶血。粘性剪切应力（VSS）超过150 Pa时有溶血风险。通过FSI精确求解湍流场，能定量评估人工瓣膜的溶血风险。
数字孪生
瓣膜的数字孪生可能吗？
将患者特异性瓣膜形状、血行动力学数据、材料参数整合，构建"心脏瓣膜数字孪生"的设想正在推进。实时状态更新需要ROM（Reduced Order Model），POD-Galerkin法和Deep Operator Network（DeepONet）等正在研究。
飞利浦公司已用HeartModel AI实用化从心脏超声自动提取瓣膜形状的技术。将其与FSI基的ROM结合，实现瓣膜功能自动评估是未来目标。
生物可吸收瓣
下一代人工瓣膜有什么技术呢？
Tissue Engineered Heart Valve（TEHV）受关注。将生物可吸收支架播种细胞，体内重塑为接近天然瓣膜的组织。FSI需考虑支架分解伴随的材料性质时间变化，需要生长·重塑连成。
组织密度 $\rho_s$ 随应力与flow-induced刺激生长（$k_g$），支架分解（$k_d$）。目前还在动物实验阶段，但FSI在设计中的角色核心。
钙化瓣膜的FSI——刚硬不均匀引起的应力集中
主动脉瓣狭窄症中瓣叶沉积钙，局部弹性率变为正常组织的10〜100倍。这种硬质钙化部位与柔软正常组织的边界产生应力集中，增加瓣叶破损风险。前沿FSI研究中，用μCT（微型CT）获取钙化的3D分布，作为不均质材料模型纳入解析的手法正在尝试。用于优化手术时机的研究。
主动脉瓣FSI的故障排查
瓣膜FSI特有的难点有几个。
1. 瓣膜关闭时的网格破坏（ALE法）
症状: 瓣叶关闭瞬间出现negative volume错误停止。
原因: 3片瓣叶接近·接触时流体单元被挤压而破坏。
对策:
- 用Overset（Chimera）网格，瓣叶周围设置组件网格
- 切换到IB法或Immersogeometric法（回避接触拓扑变化）
- ALE法中提高网格更新频率，配合重网格化
2. 瓣叶的穿透（接触失败）
症状: 瓣叶互相穿透，得到非物理结果。
原因: 接触检出的时间步长粗糙，或罚函数系数不足。
对策:
- 增加接触罚函数系数（但过大会导致振动）
- 切换到增强拉格朗日法
- 在瓣膜关闭前后将时间步长减少到1/5〜1/10
3. 非生理的瓣口面积
瓣膜开不了，或开得过大，这是为什么呢？
瓣膜开不了的情况:
- 瓣叶刚性参数过高
- 入口流量/压力波形峰值不足
- 初应力（预应力）设定不足
瓣膜开过大的情况:
- 材料刚性过低（检查拟合数据）
- 胶原纤维方向设定错误（周向与径向反向）
4. Valsalva窦内的涡流不形成
症状: 瓣膜后流无涡流观察，瓣膜关闭延迟。
对策:
- 精确建模Valsalva窦的形状（窦深和宽度很重要）
- 窦内流体网格充分细化
- 出口边界离窦足够远（最少5个瓣膜直径以下流）
| 验证项目 | 期望值 | 允许范围 |
|---|---|---|
| EOA（天然瓣） | 3.0〜4.0 cm² | ±0.5 cm² |
| 峰值流速（正常瓣） | 1.0〜1.5 m/s | ±0.3 m/s |
| 逆流率（正常瓣） | < 5% | < 10% |
| 瓣膜关闭时间 | 30〜50 ms | ±20 ms |
瓣膜FSI的验证是与心脏超声测值比较呢。临床数据的核对是品质保证的关键。
瓣膜关闭时计算发散——接触处理的陷阱
主动脉瓣FSI解析最常见的故障是"瓣膜关闭瞬间计算发散"。3片小叶在中心接触的瞬间，接触力急增导致时间步破坏。对策是"逐步提高罚函数刚性""提高接触检出频率""关闭前后缩小Δt"的组合。某研究小组实现了只在关闭前后自动调整Δt为1/10的自适应时间步，解决了这个问题。
关联主题
价值
细节
错误