双向流体结构耦合

分类：流体分析（CFD） | 综合版 2026-04-06

CAE visualization for two way fsi theory - technical simulation diagram

双向流体结构耦合

双向流体结构耦合的理论基础

概要

🧑🎓

老师！今天是关于双向流体结构耦合的讲座吧？这是怎样的东西呢？

🎓

变形→流动→力的双向迭代。强耦合/弱耦合。

控制方程

$$ \mathbf{u}_f|_{\Gamma}=\dot{\mathbf{d}}_s|_{\Gamma} $$

$$ \boldsymbol{\sigma}_f\cdot\mathbf{n}|_{\Gamma}=\boldsymbol{\sigma}_s\cdot\mathbf{n}|_{\Gamma} $$

🧑🎓

老师的解释很清楚！双向流体结构耦合的困惑已经消除了。

离散化方法

🧑🎓

这个方程，在计算机上实际上如何求解呢？

🎓

使用有限元法（FEM）进行空间离散化。建立单元刚度矩阵，构建整体刚度方程。

🎓

转换为弱形式（变分形式），使用试验函数和形状函数进行Galerkin方法的定式化。单元类型的选择（低阶单元 vs. 高阶单元、完全积分 vs. 低减积分）与解的精度和计算成本的权衡直接相关。

矩阵求解算法

🧑🎓

矩阵求解算法，具体是什么意思呢？

🎓

使用直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程。对于大规模问题，预处理迭代法很有效。

求解法	分类	内存使用	适用规模
LU分解	直接法	O(n²)	小到中等规模
Cholesky分解	直接法（对称正定）	O(n²)	小到中等规模
PCG法	迭代法	O(n)	大规模
GMRES法	迭代法	O(n·m)	大规模·非对称
AMG预处理	预处理	O(n)	超大规模

🧑🎓

也就是说，在有限元法的地方偷工减料，之后就会吃苦头吧。我铭记在心！

商业工具中的实现

🧑🎓

那么，要做双向流体结构耦合需要什么样的软件呢？

工具名称	开发者/现在	主要文件格式
Ansys Fluent	Ansys Inc.	.cas, .dat, .msh, .jou
Simcenter STAR-CCM+	Siemens Digital Industries Software	.sim, .java, .csv
Ansys CFX	Ansys Inc.	.cfx, .def, .res, .ccl
OpenFOAM	开源（OpenCFD/ESI、OpenFOAM Foundation）	字典文件（blockMeshDict等）, .foam
COMSOL Multiphysics	COMSOL AB	.mph

厂商系统和产品整合历史

🧑🎓

各个软件的来历，有没有什么戏剧化的故事呢？

Ansys Fluent

🧑🎓

接下来是Ansys Fluent的讲座吧。内容是什么呢？

🎓

Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。

现在的所属：Ansys Inc.

Simcenter STAR-CCM+

🧑🎓

接下来是Simcenter STAR的讲座吧。内容是什么呢？

🎓

CD-adapco开发。2016年被西门子收购并整合到Simcenter品牌。多面体网格是特点。

现在的所属：Siemens Digital Industries Software

Ansys CFX

🧑🎓

"Ansys CFX"请告诉我！

🎓

AEA Technology (UK) 开发的CFX。2003年被Ansys收购。耦合求解器方式是特点。

现在的所属：Ansys Inc.

🧑🎓

哦〜，开发历史的讲座，太有趣了！再多讲一些吧。

文件格式和互操作性

🧑🎓

在不同软件之间传递数据时有什么注意事项吗？

格式	扩展名	种类	概述
CGNS	.cgns	CFD数据	CFD通用记号系统。CFD结果的标准交换格式。
VTK	.vtk/.vtu	可视化	可视化工具包格式。在ParaView等中使用。

🎓

在不同求解器之间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是高阶单元和特殊单元（内聚单元、用户定义单元等）通常不能在求解器之间直接转换。

🧑🎓

原来，格式看似简单，但实际上深度很大呢。

实务中的注意事项

🧑🎓

有没有教科书里没有的"现场知识"？

🎓

网格收敛性确认、边界条件的妥当性验证、材料参数的灵敏度分析非常重要。

🎓

网格依赖性验证：用至少3个网格密度等级确认收敛性

边界条件的妥当性：设置物理上有意义的约束条件

结果验证：与理论解、实验数据、已知基准问题的比较

🧑🎓

哎呀，双向流体结构耦合真是奥深呢… 不过有了老师的讲解，已经整理得差不多了！

🎓

恩，很好的状态呢！实际上动手去做才是最好的学习。有不明白的地方随时来问我。

Coffee Break 闲话

FDA将双向FSI规定为"必须"的那一天

在2000年代初期，机械式人工心脏瓣膜的疲劳破损事故频发。当时的分析主要采用将流体和结构分开求解的单向FSI，完全忽视了"瓣膜关闭时刻流体反过来再击打瓣膜"的耦合效应。FDA由此吸取教训，在体内植入医疗器械的认可审查中要求提交双向耦合分析，使得双向FSI一夜之间从学术界走向工程实务的标准。当涉及生命安全时，计算精度要求的标准也会因此改变呢。

双向流体结构耦合的数值计算方法

数值方法详解

🧑🎓

具体用什么算法来求解双向流体结构耦合呢？

离散化的公式化

🎓

用形状函数 $N_i$ 逼近未知量：

$$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$

🎓

用公式表达就是这样。

$$ K_e = \int_{\Omega_e} B^T \, D \, B \, d\Omega \approx \sum_{g=1}^{n_g} w_g \, B^T(\xi_g) \, D \, B(\xi_g) \, |J(\xi_g)| $$

基本方程式的离散形式

🎓

用公式表达就是这样。

$$ \mathbf{u}_f|_{\Gamma}=\dot{\mathbf{d}}_s|_{\Gamma} $$

$$ \boldsymbol{\sigma}_f\cdot\mathbf{n}|_{\Gamma}=\boldsymbol{\sigma}_s\cdot\mathbf{n}|_{\Gamma} $$

🧑🎓

嗯，光是看公式不太能理解呢… 这表达的是什么呢？

🎓

将连续体的控制方程离散化后，得到如下代数方程组：

$$ [K]\{u\} = \{F\} $$

🎓

这里 $[K]$ 是总体刚度矩阵（或等效的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量呢。

🧑🎓

啊，这样啊！连续体的控制方程原来是这样离散化的呀。

单元技术

🧑🎓

"单元技术"这个词听过，但可能还没真正理解…

单元类型	次数	节点数(3D)	精度	计算成本
四面体1阶	线性	4	低（剪切锁定）	低
四面体2阶	二次	10	高	中
六面体1阶	线性	8	中	中
六面体2阶	二次	20	非常高	高
棱柱	线性/二次	6/15	中～高	中

积分方案

🧑🎓

积分方案，具体是什么意思呢？

🎓

完全积分：精确积分所有项。刚度过度评估趋势（锁定）

低减积分

：减少积分点数。提高计算效率，但存在沙漏模式发生的风险
选择低减积分 (B-bar法)：分离体积项和偏差项进行积分。回避锁定

🧑🎓
听到这儿，总算明白单元类型为何重要了！

收敛性和稳定性

🧑🎓
不收敛了的话，最先检查什么？

🎓

h-细化：细分网格（减小单元大小 h）以提高精度

p-细化：提高单元多项式次数以改进精度

hp-细化：同时优化 h 和 p

🎓
收敛速度：二阶单元以 $O(h^2)$ 阶数降低误差（光滑解的情况）

🧑🎓
原来，细分网格看似简单，但实际上深度很大呢。

求解器设置的建议

🧑🎓
具体用什么算法来求解双向流体结构耦合呢？

参数推荐值备注

迭代法收敛判定 $10^{-6}$ 残差范数基准

预处理方法 ILU(0) or AMG 取决于问题规模

最大迭代次数 1000 非收敛时应重新审视设置

内存模式 In-core 尽可能使用

风上差分（Upwind）

一阶风上：数值扩散大但稳定。二阶风上：精度提高但振荡风险。高Reynolds数流动中必需。

中心差分（Central Differencing）

二阶精度，但当Pe数 > 2时会产生数值振荡。低Reynolds数扩散支配流动时适用。

TVD方案（MUSCL、QUICK等）

使用限制函数抑制数值振荡的同时维持高精度。对冲击波和陡峭梯度的捕捉有效。

有限体积法 vs 有限元法

FVM：自然满足守恒律。CFD主流。FEM：复杂形状·多物理场有利。SPH等无网格法也在发展。

CFL条件（库朗数）

显式法：CFL ≤ 1为稳定条件。隐式法：CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：1个时间步内信息不超过1个网格。

残差监控

连续方程、动量、能量各残差下降3～4个数量级时判定为收敛。质量守恒的残差尤为重要。

松弛系数

压力：0.2～0.3，速度：0.5～0.7为一般初值。发散时降低松弛系数。收敛后提高以加速。

非定常计算的内部迭代

各时间步内迭代至定常解。内部迭代数：5～20次为目标。残差在时间步间波动时应重新检视时间步长。

参数	推荐值	备注
迭代法收敛判定	$10^{-6}$	残差范数基准
预处理方法	ILU(0) or AMG	取决于问题规模
最大迭代次数	1000	非收敛时应重新审视设置
内存模式	In-core	尽可能使用

双向流体结构耦合的实务应用

实践指南

🧑🎓
老师，请讲讲"实践指南"！

🎓
解说双向流体结构耦合的实务分析流程和注意事项。

🧑🎓
听到这儿，总算明白为什么双向流体结构耦合很重要了！

分析流程

🧑🎓
从第一步开始就教我！该从什么开始呢？

🎓
1. 前处理 (Pre-processing)

导入CAD数据并简化形状

定义材料特性

网格生成（确定单元类型和大小）

设置边界条件和荷载条件

🎓
2. 求解 (Solving)

求解器设置（方法、收敛标准、输出控制）

投入作业并执行计算

收敛监控

🎓
3. 后处理 (Post-processing)

结果可视化（位移、应力、其他物理量）

结果验证和妥当性确认

报告编制

网格生成的最佳实践

🧑🎓
网格的好坏怎样判断呢？

单元品质指标

🧑🎓
请讲讲"单元品质指标"！

指标理想值允许范围影响

纵横比 1.0 < 5.0 精度下降

雅可比比 1.0 > 0.3 单元退化

翘曲 0° < 15° 精度下降

歪斜度 0° < 45° 收敛性恶化

锥角比 0 < 0.5 精度下降

网格密度的确定

🧑🎓
网格密度的确定，具体是什么意思呢？

🎓

应力集中部：至少配置3层以上的单元

应力梯度大的区域：单元大小改为周围1/3～1/5

荷载施加点附近：局部细分

远场区域：使用粗网格以确保计算效率

边界条件的设置指南

🧑🎓
听说边界条件弄错的话全部都完蛋呢…

🎓

注意过约束：刚体移动的约束只需6个自由度

活用对称条件：减少计算规模

荷载的等价分配：集中荷载 vs. 分布荷载的选择

🧑🎓
啊，这样啊！过约束注意就是这样的机制呢。

按商业工具的实现步骤

🧑🎓
有各种各样的软件吧？各自的特点请教教我！

工具名称开发者/现在主要文件格式

Ansys Fluent Ansys Inc. .cas, .dat, .msh, .jou

Simcenter STAR-CCM+ Siemens Digital Industries Software .sim, .java, .csv

Ansys CFX Ansys Inc. .cfx, .def, .res, .ccl

OpenFOAM 开源（OpenCFD/ESI、OpenFOAM Foundation）字典文件（blockMeshDict等）, .foam

COMSOL Multiphysics COMSOL AB .mph

Ansys Fluent

🧑🎓
接下来是Ansys Fluent的讲座吧。内容是什么呢？

🎓
Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。

现在的所属：Ansys Inc.

Simcenter STAR-CCM+

🧑🎓
接下来是Simcenter STAR的讲座吧。内容是什么呢？

🎓
CD-adapco开发。2016年被西门子收购并整合到Simcenter品牌。多面体网格是特点。

现在的所属：Siemens Digital Industries Software

🧑🎓
老师的解释很清楚！工具名的困惑已经消除了。

常见失败和对策

🧑🎓
初学者容易犯什么样的失败呢？想事先知道！

症状原因对策

计算不收敛网格品质不良、不适当的边界条件网格改善、约束条件重新审视

应力异常大应力奇点、网格依赖奇点回避、局部网格细分

位移非现实材料常数错误、单位系不一致确认输入数据

计算时间过大不必要的细分、低效求解法网格优化、并行计算

质量保证检查清单

🧑🎓
有没有教科书里没有的"现场知识"？

🎓

用3个以上的网格密度等级确认收敛性了吗

验证了力的平衡（反力合计）吗

确认结果在物理合理的范围内了吗

与已知理论解或基准问题做过比较吗

🧑🎓
哎呀，双向流体结构耦合真是奥深呢… 不过有了老师的讲解，已经整理得差不多了！

🎓
恩，很好的状态呢！实际上动手去做才是最好的学习。有不明白的地方随时来问我。

Coffee Break 闲话

从双向FSI首次理解"血管能主动振动"的事
在用双向FSI评估腹主动脉瘤（血管局部膨胀）破裂风险的过程中，发现了有趣的事实。血管壁不是单纯的弹性管，而是根据血流脉动主动变形，以此来降低流体阻力。单向FSI由于"将壁视为不动"来分析，完全忽视了这个效应。通过双向FSI，血管的"聪慧性"第一次得到了定量评估，这也推进了人工血管设计的改进。

双向流体结构耦合的软件比较

商业工具比较

🧑🎓
有各种各样的软件吧？各自的特点请教教我！

🎓
解说支持双向流体结构耦合的主要商业CAE工具的功能比较，以及各产品的历史背景。

适用工具列表

🧑🎓
那么，要做双向流体结构耦合需要什么样的软件呢？

工具名称开发者/现在主要文件格式

Ansys Fluent Ansys Inc. .cas, .dat, .msh, .jou

Simcenter STAR-CCM+ Siemens Digital Industries Software .sim, .java, .csv

Ansys CFX Ansys Inc. .cfx, .def, .res, .ccl

OpenFOAM 开源（OpenCFD/ESI、OpenFOAM Foundation）字典文件（blockMeshDict等）, .foam

COMSOL Multiphysics COMSOL AB .mph

Ansys Fluent

🧑🎓
接下来是Ansys Fluent的讲座吧。内容是什么呢？

🎓
Fluent Inc.开发。2006年被Ansys收购。基于非结构化网格的通用CFD求解器。

现在的所属：Ansys Inc.

Simcenter STAR-CCM+

🧑🎓
接下来是Simcenter STAR的讲座吧。内容是什么呢？

🎓
CD-adapco开发。2016年被西门子收购并整合到Simcenter品牌。多面体网格是特点。

现在的所属：Siemens Digital Industries Software

🧑🎓
听到这儿，总算明白为什么开发历史很重要了！

Ansys CFX

🧑🎓
"Ansys CFX"请告诉我！

🎓
AEA Technology (UK) 开发的CFX。2003年被Ansys收购。耦合求解器方式是特点。

现在的所属：Ansys Inc.

OpenFOAM

🧑🎓
OpenFOAM，具体是什么意思呢？

🎓
Imperial College London发源的开源CFD。OpenCFD Ltd（ESI Group旗下）和The OpenFOAM Foundation并行开发。

现在的所属：开源（OpenCFD/ESI、OpenFOAM Foundation）

🧑🎓
啊，这样啊！开发机制就是这样呢。

功能比较矩阵

🧑🎓
预算和时间都有限，成本最强的是哪个？

功能 Fluent Star-CCM+ CFX OpenFOAM

基本功能 ○ ○ ○ ○

高级功能 ○ ○ ○ △

自动化/脚本 ○ ○ ○ ○

并行计算 ○ ○ ○ ○

GPU支持 △ △ △ ○

转换时的风险

🧑🎓
转换时的风险，具体是什么意思呢？

🎓

单元类型的不兼容：求解器特有单元无法用中立格式表现

材料模型的差异：同名但内部实现不同的情况存在

边界条件的重定义：大多情况下需要手动重设

结果数据的比较：输出变量的定义（节点值 vs. 单元值、积分点值）有差异

🧑🎓
啊，这样啊！不同工具之间的数据转换就是这样的机制呢。

许可证形式

🧑🎓
听说过"许可证形式"，但可能还没真正理解…

工具许可证特点

商用FEA 节点锁定/浮动高价但有官方支持

OpenFOAM GPL 免费但支持有偿

COMSOL 节点锁定/浮动按模块购买

Code_Aster GPL EDF开发的OSS求解器

选择指南

🧑🎓
到底要选哪个，评判标准能教教我吗？

🎓
在双向流体结构耦合的工具选择中应考虑：

🎓

分析规模：数万～数亿DOF的可扩展性

物理模型：所需构成则单元类型的对应情况

工作流：与CAD的联动、自动化的容易度

成本：初期投资 + 年度保守 + 教育成本

支持：技术支持的质量和响应

🧑🎓
哎呀，双向流体结构耦合真是奥深呢… 不过有了老师的讲解，已经整理得差不多了！

🎓
恩，很好的状态呢！实际上动手去做才是最好的学习。有不明白的地方随时来问我。

Coffee Break 闲话

ANSYS SystemCoupling普及了"工具间联动"文化

在双向FSI商业工具的普及中，ANSYS SystemCoupling的角色很大。此前Fluent（流体）和Mechanical（结构）是分开的软件，数据联动需要自制脚本。SystemCoupling在2010年代初期成熟，充当两者的"中介"时，非专业程序员也能运转双向FSI的环境整备完毕。"向非专家也开放多物理场"的思想，如今已经波及其他公司的工具。

双向流体结构耦合的先进研究

先进课题和研究趋势

🧑🎓
双向流体结构耦合的领域，今后会怎样发展呢？

🎓
来看双向流体结构耦合的最新研究动向和先进方法。

最新数值方法

🧑🎓
接下来是最新数值方法的讲座吧。内容是什么呢？

🧑🎓
嗯，光是看公式不太能理解呢… 这表达的是什么呢？

🎓

等几何分析 (IGA)：直接使用NURBS基函数，实现CAD-CAE间的无缝联动

粒子法 (SPH, MPM)：无网格方法追踪大变形、破壊

位相场法 (Phase-Field)：用隐式表现界面，复杂界面追踪

机器学习支援：代理模型、物理约束神经网络 (PINN)

高性能计算 (HPC) 适配

并行化方法概述适用求解器

MPI (领域分割) 分布内存型。大规模问题的标准全主要求解器

OpenMP 共享内存型。节点内并行多数求解器

GPU (CUDA/OpenCL) GPGPU活用。特别在显式法有效 LS-DYNA, Fluent等

混合 MPI+OpenMP 节点间+节点内并行大规模HPC环境

双向流体结构耦合的故障排除

故障排除

常见错误和对策

🧑🎓
老师也因为双向流体结构耦合而通宵调试过吗？（笑）

1. 收敛失败

🧑🎓
收敛失败，具体是什么意思呢？

🎓
症状：求解器未在指定迭代数内收敛而异常终止

🎓
可能的原因：

网格品质不足（过度歪斜的单元）

材料参数设置不当

不当的初始条件

非线性性过强（

指标	理想值	允许范围	影响
纵横比	1.0	< 5.0	精度下降
雅可比比	1.0	> 0.3	单元退化
翘曲	0°	< 15°	精度下降
歪斜度	0°	< 45°	收敛性恶化
锥角比	0	< 0.5	精度下降

症状	原因	对策
计算不收敛	网格品质不良、不适当的边界条件	网格改善、约束条件重新审视
应力异常大	应力奇点、网格依赖	奇点回避、局部网格细分
位移非现实	材料常数错误、单位系不一致	确认输入数据
计算时间过大	不必要的细分、低效求解法	网格优化、并行计算

功能	Fluent	Star-CCM+	CFX	OpenFOAM
基本功能	○	○	○	○
高级功能	○	○	○	△
自动化/脚本	○	○	○	○
并行计算	○	○	○	○
GPU支持	△	△	△	○

工具	许可证	特点
商用FEA	节点锁定/浮动	高价但有官方支持
OpenFOAM	GPL	免费但支持有偿
COMSOL	节点锁定/浮动	按模块购买
Code_Aster	GPL	EDF开发的OSS求解器

并行化方法	概述	适用求解器
MPI (领域分割)	分布内存型。大规模问题的标准	全主要求解器
OpenMP	共享内存型。节点内并行	多数求解器
GPU (CUDA/OpenCL)	GPGPU活用。特别在显式法有效	LS-DYNA, Fluent等
混合 MPI+OpenMP	节点间+节点内并行	大规模HPC环境