一方向流体構造連成
理论与物理
单向FSI的概念
老师，单向（One-Way）FSI和双向（Two-Way）FSI有什么区别呢？
在单向FSI中，仅将流体分析得到的压力、剪切力作为结构分析的载荷条件进行传递，并假设结构变形不会影响流体。也就是说信息流是“流体→结构”的单向流动。
这个假设在结构变形不会显著改变流场的情况下是合理的。具体来说：
1. 结构位移相对于特征长度非常小（$\delta/L \ll 1$）
2. 结构固有频率与流体激励频率充分分离
3. 结构密度相对于流体密度足够大（$\rho_s / \rho_f \gg 1$）
具体可以用于哪些问题呢？
在以下场景中，单向FSI足以获得足够的精度。
| 应用 | 流体 | 结构 | 评估项目 |
|---|---|---|---|
| 汽车外装面板 | 行驶风 | 钢板面板 | 风压引起的变形与应力 |
| 管道的流体力 | 内部流动 | 管道结构 | 压力载荷引起的应力 |
| 建筑物的风荷载 | 风 | 钢结构/钢筋混凝土 | 风压分布、层间位移 |
| 涡轮机械叶片 | 稳态空气动力 | 叶片 | 离心力+气动力的静应力 |
| 电子设备的热变形 | 冷却气流 | 基板 | 温度分布→热应力 |
也就是说，在变形微小且无需反馈给流场的情况下使用对吧。计算成本也比双向FSI大幅降低吗？
成本要低得多。因为只需要分别求解一次CFD和FEA，与双向FSI（需要迭代求解数十次）相比，计算时间大约是其1/10到1/100的量级。对于设计初期的筛选或结构强度的概算来说非常高效。
桥梁设计中“省略”单向FSI的代价
单向FSI是“流体对结构影响较小时”有效的简化方法，但如果误判了这一点，后果会很严重。实际上，在某座吊桥的设计中，将流体力分析和结构分析完全独立进行，并用静力换算处理流体载荷的动态成分，结果在投入使用后，强风时的振动大幅超过了设计预测。这就是“分开求解答案应该相同”的固有观念导致的失败。在选择单向FSI之前，必须定量确认其假设的合理性。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水流会不稳定地喷溅，过一会儿才会变成稳定的水流，对吧？描述这个“变化过程中”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非稳态现象。那么稳态分析是什么？就是只看“经过足够时间流动稳定之后”——也就是将此项设为零。计算成本会大幅下降，因此先用稳态求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶扔进河里会怎样？会被水流带着往下游漂，对吧？这就是“对流”——流体运动搬运物体的效果。暖风的暖气能到达房间的另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含了“速度×速度”，因此是非线性的。也就是说，流速变快时此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水。因为蜂蜜的粘度（$\mu$）高，所以不易流动。粘度越大，扩散项越强，流体的运动就越“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：推动注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差就是推动流体的力。大坝放水也是同样的原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD中的“压力”通常指表压而非绝对压力。切换到可压缩分析时结果突然变得奇怪，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：受热的空气会上升——为什么呢？因为比周围空气轻（密度低），所以被浮力推上去。这个浮力作为源项添加到方程中。此外，燃气灶的火焰产生化学反应热、工厂的电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项来表示。忘记源项会怎样？在自然对流分析中忘记加入浮力，流体就完全不动——就像冬天在房间里开了暖气，暖空气却不上浮一样，会得到物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用的情形：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等方法）
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 从入口条件的体积流量换算时，注意截面面积的单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析中使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气：约1.225 kg/m³@20°C，水：约998 kg/m³@20°C |
| 粘度系数 $\mu$ | Pa·s | 注意与运动粘度系数 $\nu = \mu/\rho$ [m²/s] 的混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判断指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数据传递方法
将CFD结果传递给FEA时，网格应该不同吧，数据是怎么传递的呢？
需要从CFD网格（通常较细）到FEA网格（通常较粗）进行数据映射。主要方法如下。
| 方法 | 原理 | 精度 | 守恒性 |
|---|---|---|---|
| 最近邻 | 复制最近邻点的值 | 低 | 无 |
| 反距离加权 | 按距离倒数进行加权 | 中 | 无 |
| 守恒映射 | 面积/体积加权 | 高 | 有（力的总和守恒） |
| 分布保持 | 基于形函数的插值 | 高 | 无（分布形状保持） |
| RBF插值 | 径向基函数插值 | 非常高 | 取决于设置 |
对于力（压力、剪切力）的传递，守恒映射很重要。这保证了CFD面上的合力与FEA面上的合力严格一致。Ansys System Coupling或STAR-CCM+的Co-Simulation中都提供了这个选项。
没有守恒性会怎样？
例如，将升力100 N的CFD结果传递给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格式输出壁面压力，然后作为Imported Pressure导入Ansys Mechanical或Abaqus。映射精度取决于软件的导入功能，因此传递后确认力的总和很重要。
“数据映射”的精度决定单向FSI的命运
单向FSI中最容易被忽视的问题是接口的数据映射。CFD网格和结构网格通常形状和密度都不同，节点不匹配。如果对力或压力的插值映射不加考虑地使用最近邻插值，可能会产生非平凡的误差，导致载荷积分值偏差几个百分点。有报告称，这会使结构分析的峰值应力产生10%~20%的误差。“数据传递是简单作业”这种认识，是导致意想不到精度下降的原因。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必不可少。
中心差分（Central Differencing）
二阶精度，但Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
FVM：自然满足守恒定律。CFD的主流。FEM：对复杂形状、多物理场有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式方法：CFL ≤ 1为稳定条件。隐式方法：即使CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：一个时间步内信息传播不超过一个网格。
残差监控
连续性方程、动量、能量各项残差下降3~4个数量级可判断为收敛。质量守恒的残差尤其重要。
松弛因子
压力：0.2~0.3，速度：0.5~0.7为一般初始值。发散时降低松弛因子。收敛后可提高以加速。
非稳态计算的内部迭代
在每个时间步内迭代直到收敛到稳态解。内部迭代次数：5~20次为参考值。如果残差在时间步之间波动，则需要重新审视时间步长。
SIMPLE法的比喻
SIMPLE法是“交替调整”的方法。先假设求解速度（预测步），然后根据该速度修正压力以满足质量守恒（修正步），再用修正后的压力修正速度——重复这种“投接球”过程以逼近正确答案。类似于两人调整架子水平的作业：一人调整高度，另一人调整平衡，如此交替进行。
迎风格式的比喻
迎风格式是“站在河流中重视上游信息”的方法。站在河里的人看下游也看不到水的来源。
相关主题
なった
詳しく
報告