单向流体结构耦合
单向流体结构耦合理论基础
单向FSI的概念
老师，单向（One-Way）FSI与双向（Two-Way）FSI有什么区别？
在单向FSI中，流体分析得到的压力和剪切力作为结构分析的荷载条件转移，假设结构变形对流体没有影响。换言之，信息流是"流体→结构"的单向流动。
这个假设在以下条件下是合理的：
1. 结构的位移相对于特征长度非常小（$\delta/L \ll 1$）
2. 结构的固有振动频率与流体激励频率充分分离
3. 结构密度相对于流体密度足够大（$\rho_s / \rho_f \gg 1$）
具体可以用在什么问题上？
以下场景可以通过单向FSI获得充分的精度。
| 应用 | 流体 | 结构 | 评估项 |
|---|---|---|---|
| 汽车外装面板 | 行驶风 | 钢板面板 | 风压变形、应力 |
| 管道流体力 | 内部流动 | 管道结构 | 压力荷载引起的应力 |
| 建筑风荷载 | 风 | 钢结构/混凝土 | 风压分布、层间位移 |
| 涡轮机器叶片 | 定常气动力 | 叶片 | 离心力+气动静应力 |
| 电子设备热变形 | 冷却气流 | 基板 | 温度分布→热应力 |
也就是说，变形很小且不需要回传到流动场的情况下就用单向。计算成本也比双向便宜得多？
便宜多了。只需各解一遍CFD和FEA，而双向FSI（反复多次求解）的计算时间是1/10到1/100的量级。设计初期的筛选或强度粗略估算来说效率很高。
桥梁设计中"省略"单向FSI的代价
单向FSI是"流体对结构的影响较小时"的有效简化，但如果判断错误会造成严重后果。实际上有一座吊桥的设计中，将流体力和结构分析完全独立进行，用静态换算处理流体荷载的动态分量，结果投入使用后强风时的振动大幅超过设计预测。这是"分开解答应该和在一起解答结果相同"这种假设造成的失败。在选择单向FSI前，必须定量验证其假设的有效性。
单向流体结构耦合的数值计算方法
数据转移方法
将CFD结果转到FEA时，网格不同的话怎么转移数据呢？
需要从CFD网格（通常较细）到FEA网格（通常较粗）的数据映射。主要方法如下。
| 方法 | 原理 | 精度 | 守恒性 |
|---|---|---|---|
| 最近邻点法 | 复制最近点的值 | 低 | 无 |
| 距离倒数加权法 | 按距离倒数加权 | 中 | 无 |
| 守恒映射法 | 按面积/体积加权 | 高 | 有（力的总和守恒） |
| 轮廓保留法 | 基于形状函数的插值 | 高 | 无（分布形状守恒） |
| RBF插值 | 径向基函数插值 | 非常高 | 取决于设置 |
在力（压力和剪切力）转移中，守恒映射法很重要。它保证CFD面的力的合计与FEA面的力合计严格一致。Ansys System Coupling和STAR-CCM+的Co-Simulation提供了这个选项。
没有守恒性会怎样？
比如CFD结果升力100 N转移到FEA时，FEA上的积分值可能变成95 N或105 N。结构整体的力平衡破裂了，变形和反力结果会出现误差。特别是在整体荷载是设计判断标准的情况下（比如叶片根部应力），守恒映射是必须的。
时间序列数据转移
非定常CFD结果转到结构时怎么处理？
有两个方法。
1. 瞬时值转移：将CFD每个时间步的压力分布输入到FEA的瞬态分析。对所有时间步进行映射。精度高但数据量巨大。
2. 统计量转移：转移时间平均压力和压力变动的RMS（或PSD）。定常荷载用时间平均处理，疲劳荷载用RMS方法评估。数据少，更实用。
PSD（功率谱密度）怎样转换成结构的疲劳荷载？
用随机振动分析（Nastran SOL 111/112），以压力PSD为输入，计算应力的PSD，然后用Dirlik法或Steinberg法推算疲劳寿命。火箭整流罩音响荷载分析等广泛应用这个方法。
基于文件 vs 软件耦合
单向FSI的数据交换可以通过文件吗？
单向FSI用文件方式足够实用。
| 方法 | 优点 | 缺点 |
|---|---|---|
| 基于文件（CSV、CGNS） | 简洁、软件无依赖 | 手工映射、无守恒性保证 |
| Ansys Workbench内部 | GUI完整、自动映射 | 仅限Ansys产品 |
| System Coupling | 守恒性保证、自动化 | 需要许可证 |
| Python脚本 | 灵活、可自动化 | 需要自行实现 |
文件方式的话，从Fluent输出CGNS/EnSight格式的壁面压力，导入到Ansys Mechanical或Abaqus作为导入压力。映射精度取决于软件的导入功能，转移后一定要验证力的总和。
"数据映射"的精度掌控着单向FSI的成败
单向FSI中最容易被忽视的是接口数据映射的精度。CFD网格与结构网格通常使用不同的形状和密度，节点不一致。力或压力的补间映射若草率地用最近邻点法，会产生非平凡的误差，使荷载积分值偏差数个百分点。这可能导致结构分析的峰值应力误差达10~20%。"数据转移是简单的工作"这种认识，往往是精度劣化的隐因。
单向流体结构耦合的实际应用
Ansys Workbench中的单向FSI
请详细说说Ansys Workbench执行单向FSI的具体步骤。
在Workbench项目示意图中按如下方式配置。
1. 放置Fluent Analysis（System A），完成CFD分析
2. 放置Static Structural（System B）
3. 将System A的Solution单元拖拽到System B的Setup单元
4. Mechanical中自动生成导入压力
5. 验证映射，执行结构分析
就拖拽一下就能连接，很方便呢。
Workbench的优势。但也要注意几点。
- CFD和FEA网格应独立优化（CFD重在边界层分辨，FEA重在应力集中区）
- 确认导入压力的映射目标面与CFD壁面空间对齐
- 转移后在Mechanical中验证力反应值与CFD合力一致
STAR-CCM+ + Abaqus Co-Simulation
从STAR-CCM+转到Abaqus怎么做？
STAR-CCM+有导出功能可把壁面数据输出为Abaqus INP格式。步骤如下。
1. STAR-CCM+完成CFD分析
2. 壁面零件→右键→导出→选择Abaqus Surface
3. 在Abaqus中将输出的.INP文件作为*DSLOAD（分布荷载）读入
4. 执行结构分析
非定常数据也能这样转吗？
非定常的话逐步输出INP，在Abaqus中定义为振幅表格数据。或者用Co-Simulation Engine自动同步时间。
单向FSI的合理性验证
怎样判断单向FSI的结果精度是否足够？
有简单的检查方法。从单向FSI得到结构的位移 $\delta$，检查
的话，单向FSI假设基本合理。超过5-10%就需要双向FSI。
还有一个指标是无量纲位移参数 $\Pi$。
这里 $q_{\infty}$ 是动压，$E$ 是杨氏模量，$L/t$ 是特征长度/厚度比。$\Pi \ll 1$ 时单向足够，$\Pi \sim O(1)$ 就需要双向。
实际上是先用单向跑一遍，位移太大就换成双向，这样比较实用对吧。
完全同意。单向FSI计算成本低，先做筛选，再判断是否需要双向，这样很高效。
电子基板热应力分析中单向FSI每天都在用
单向FSI在最平凡但用途最广的地方是电子设备的热应力分析。先用CFD（或热传导分析）求基板温度分布，再转到FEM结构分析计算焊点热应力和疲劳。智能手机、EV功率模块等数十万至数百万件产品中，这个单向FSI工作流是可靠性保证的核心。不起眼，但可能是全球计算最频繁的FSI分析。
单向流体结构耦合的软件比对
热-结构的单向耦合
把温度分布转到结构算热应力，这也是单向FSI的一种吗？
完全是。CHT（共轭传热）+热应力分析就是单向FSI的典范。工作流如下。
用Fluent算CHT得到固体域温度分布，转到Ansys Mechanical进行热应力分析。Workbench的话温度数据转移自动化了。
能同时转温度和压力吗？
可以。从同一个Fluent分析，Mechanical可以同时转温度和压力。导入压力和导入体温度两个荷载会被设置。在涡轮叶片分析中，离心力+气动+热应力3个荷载同时考虑是标准做法。
应用案例
单向FSI的产业应用案例有哪些？
举几个代表性例子。
| 领域 | 对象 | CFD荷载 | 结构评估 | 工具 |
|---|---|---|---|---|
| 汽车 | 后视镜 | 行驶风压 | 振动、变形 | Fluent→Mechanical |
| 航空 | 发动机短舱 | 巡航气动力 | 疲劳寿命 | STAR-CCM+→Nastran |
| 土木 | 高层建筑 | 风荷载 | 层间位移 | CFD→ETABS |
| 电力 | 冷却塔 | 内外表面风压 | 壳体应力 | Fluent→ABAQUS |
| 半导体 | 晶圆 | 冷却气流 | 热变形 | Fluent→Mechanical |
后视镜的振动分析很接近生活呢。
汽车行驶中后视镜的振动直接关系到驾驶视野清晰度，是重要评价指标。用CFD算压力变动的PSD，再用FEA的随机振动分析预测振幅。目标是让镜子的固有频率远离风压变动的卓越频率。
Python自动化
能用Python把单向FSI工作流自动化吗？
用Ansys PyFluent和Ansys PyMechanical可以完全用Python自动化。
大体流程是这样的。
1. 用PyFluent执行CFD分析，壁面压力输出CGNS
2. 用Python把压力数据映射到FEA网格上（scipy.interpolate.RBFInterpolator很方便）
3. 用PyMechanical施加荷载执行结构分析
4. 结果评估和报告生成一气呵成
参数研究也很容易吧。
流速5个条件、迎角3个条件，共15个工况用for循环跑一遍。每个工况的最大应力、位移汇总到CSV，用Pandas/Matplotlib可视化设计空间。效率很高。
热单向FSI改变了EV动力总成设计
电动汽车的电机和逆变器中，用CFD计算冷却水路的热传递系数，转到结构FEM算热应力和热变形，这种热单向FSI的应用急速普及。传统燃油车时代，冷却分析和结构分析通常是分开的，但EV的发热密度局部极高，决定了哪里会变热，直接影响可靠性，所以连成评估变成了必需。EV的转型改变了CAE工作流。
单向流体结构耦合的前沿研究
声学荷载引起的结构响应
火箭整流罩上的音响荷载分析也是单向FSI吗？
完全是。火箭发射时产生的高强度声环境（约140~160 dB）给整流罩内搭载设备造成振动荷载。这个分析称为音响-结构耦合（Vibroacoustics）。
工作流是这样的。
1. 从火箭发动机尾喷气的音响功率用经验式推算（NASA SP-8072等）
2. 用BEM（边界元法）或FEM声学模型算音场
3. 音压转移到整流罩结构的FEA模型
4. 随机振动分析（Nastran SOL 111）计算搭载设备响应
音场也能用CFD直接算吗？
近年开始用LES算排气喷流近场声源，再用FWH面积分法算远场音压。但频率范围扩展到数kHz，网格分辨要求苛刻。实际上用经验音响功率谱更经济。
风荷载在设计规范中的应用
建筑领域CFD风荷载的应用现状如何？
建筑规范和各国设计规范（Eurocode 1、ASCE 7等）规定了风荷载系数，但特殊形状建筑或相邻建筑干扰效应用CFD评估的比例在增加。
单向FSI的流程是这样。
1. 用LES/DES算建筑周围非定常风压分布
2. 抽取壁面压力的PSD和空间相关性
3. 风荷载加到结构FEM模型
4. 评估层间位移和峰值加速度
CFD结果直接用于结构设计的荷载吗？
严格来说，把CFD替代风洞试验用于设计需要认证。很多国家仍然把风洞当基准，但CTBUH（国际高层建筑与城市居住委员会）指南在一定条件下（CFD与风洞对比验证）认可CFD设计应用的趋势在增长。
单向与双向的边界
再详细一点说单向FSI的适用边界。
根据Causin et al.（2005）的研究，给出这样的判断标准。
| 参数 | 单向足够 | 需要双向 |
|---|---|---|
| 位移/特征长度 $\delta/L$ | < 1% | > 5% |
| 密度比 $\rho_s/\rho_f$ | > 100 | < 10 |
| 换算速度 $U_r = U/(f_n D)$ | 涡激锁定域外 | 锁定域内 |
| 质量比 $m^* = m/(\rho_f D^2 L)$ | > 10 | < 2 |
水中结构基本都得用双向了。
完全同意。海洋工程的立管、系泊浮体、桥梁缆索等密度比小，涡激振动（VIV）的锁定现象是问题，一定要双向FSI。但航机外板和建筑物在空气中，一般单向就够。
声学FSI中的"看不见的振动"
管道内高速流体产生的压力变动作为声波激励管壁，这是声学FSI，特别是在原子能发电厂和航空航天燃料系配管中是严重风险。有趣的是，变形通常在"听不到声音"之前就已经开始疲劳损伤。用声学FSI（声→结构）做振动疲劳评估，在原子能的定期检查中实际应用了，是为数不多的"看不见的振动"定量化方法之一。
单向流体结构耦合的故障排除
映射错误
CFD转FEA的数据映射出错了...
常见原因和对策整理如下。
| 现象 | 原因 | 对策 |
|---|---|---|
| 有未映射节点 | CFD和FEA面在空间上偏移 | 统一CAD坐标系 |
| 某区域压力为零 | CFD壁面比FEA面小 | 适度扩大CFD壁面 |
| 力合计不一致 | 非保守映射 | 改用守恒映射 |
| 压力分布多噪声 | 网格尺度差异大 | 导出前粗化CFD壁面网格 |
坐标系偏移很容易忽视呢。
Fluent和Mechanical的Workbench连联的话坐标自动对齐，但异软件间（STAR-CCM+→Abaqus等）需要在CAD导入时对齐坐标系。单位系不同（m vs mm）也要注意。
压力滤波
非定常CFD的压力噪声很多，结构分析变得不稳定...
用以下滤波方法有效。
1. 时间滤波：结构固有振动频率的10倍以上的高频用低通滤波除去
2. 空间滤波：CFD壁面网格比结构网格细得多时，做空间平均
3. 模态滤波：对结构响应无贡献的高阶模态荷载分量除去
时间滤波的截止频率怎么定？
关注的结构固有频率 $f_n$ 基础上，$f_{cutoff} = 5\sim 10 \times f_n$ 左右。比如 $f_n = 100$ Hz 的结构就设 $f_{cutoff} = 500\sim 1000$ Hz。这个频率以上的压力变动基本激励不了结构，删了也没影响。
FEA边界条件不足
CFD压力加到结构后位移特别大...
检查结构的约束条件。CFD的压力是绝对压还是相对压要确认。管道内面加相对压时，外面忘了加大气压，差压就过大了。
检查清单：
1. CFD压力输出是相对压还是绝对压（Fluent：注意Reference Pressure）
2. 结构模型有没有适当的约束条件（支撑点、螺栓孔等）
3. FEA材料参数单位与FEA单位一致吗（Pa vs MPa）
4. FEA反力合计与CFD荷载合计对不对
Fluent的Reference Pressure是零的话，输出压力是相对压吧。
正确。Operating Pressure是0就相对压=绝对压。Operating Pressure是101325 Pa就壁面输出是相对压（大气基准）。转给结构时这个设置得弄清楚，很重要。
"荷载对、变形不对"——单向FSI的盲点
单向FSI分析中的常见困惑是"流体荷载合力与实验一致，但结构变形不符"。原因多是荷载空间分布太粗。合力对了，但压力分布峰值位置微微偏差，矩就改变了，变形模态也就跟着变。为了降成本而粗化CFD网格后，荷载局部分布丢失的例子很常见。"只验证荷载合力就安心"，这是单向FSI的典型陷阱。
相关主题
价值
更详细
错误