气动弹性分析 — 故障排除
更充实的内容请查看 aeroelasticity.html。
耦合计算发散
CFD-CSD耦合计算发散了,原因是什么?
气动弹性耦合发散原因TOP 3。
| 原因 | 症状 | 对策 |
|---|---|---|
| 附加质量效应 | 网格变位发散 | 将Under-relaxation降低至0.3-0.5 |
| 网格变形破裂 | 负体积单元发生 | Laplacian smoothing使用diffusivity=inverseDist |
| 时间步长不匹配 | 解呈振荡发散 | 流体与结构的时间步长保持一致 |
附加质量效应具体是什么意思?
机翼表面运动时,周围流体也随之运动。这股流体的惯性作为"附加质量"作用于结构。密度比 $\rho_s / \rho_f$ 越小(结构越轻),附加质量的影响越大,显式耦合方案(一步内流体→结构→流体)就越容易不稳定。
对策如下。
1. 隐式耦合:每个时间步内运行FSI子迭代(System Coupling的coupling iterations)
2. Aitken加速:动态优化under-relaxation因子
3. Robin-Neumann法:界面条件引入Robin条件以增强稳定性
网格变形质量劣化
机翼变形量较大时,网格会被压扁...
网格变形方法的选择是关键。
| 方法 | 大变形耐受性 | 计算成本 | 实现软件 |
|---|---|---|---|
| Spring analogy | 低~中 | 低 | Fluent |
| Diffusion-based smoothing | 中~高 | 中 | Fluent, STAR-CCM+ |
| RBF(径向基函数) | 高 | 高 | OpenFOAM(RBFMotionSolver) |
| Overset mesh | 非常高 | 高 | Fluent, STAR-CCM+ |
当预期出现大变形时,Overset mesh(Chimera法)最安全。翼周围的网格作为独立体进行移动和旋转,背景网格在重叠区域进行插值。无需网格变形,品质不会下降。
颤振速度在模型间差异很大
DLM和CFD的颤振速度相差超过20%,应该相信哪个?
在跨音速范围内,相差20%并不少见。DLM基于势流假设,冲击波效应没有精确反映。另一方面,CFD的结果也依赖于湍流模型和网格密度。
判断标准如下。
1. 首先确认CFD的定常压力分布与实验数据(如有)一致
2. 确认CFD的颤振结果已进行网格收敛验证(两个以上网格水平)
3. 在跨音速范围内以CFD结果为基础进行判断,但理解DLM差异的原因
4. 认证时以DLM为基础,将CFD作为补正数据活用是行业标准
最后认证还是以DLM为基础啊。
是的。不过FAA/EASA也在推进CFD颤振分析认证指南的整备。将来CFD有可能成为主角。
颤振分析"收敛后答案反向"的陷阱
气动弹性分析的常见故障中,固有模式符号或位相设置错误导致计算继续进行的情况意外较多。颤振速度的计算中,流体与结构耦合矩阵的特征值虚部符号至关重要,符号错误会导致"稳定"与"不稳定"反向。"分析结果表明有余量,但实机试验发生颤振"的局面是任何担当者的噩梦。耦合分析以符号·位相确认为生命。
相关主题
更详细
错误