单体（一体型）求解法 — 故障排除指南
更充实的内容可在 monolithic-approach.html 中查看。
故障排除
原来如此。那么如果一体型做成了，基本上就没问题了，对吧？
常见错误及对策
先生也在单体（一体型）求解法中通宵调试过吗？（笑）
1. 收敛失败
收敛失败具体是怎样的情况？
症状: 求解器在指定迭代次数内未收敛，异常终止
可能的原因:
- 网格质量不足（过度扭曲的单元）
- 材料参数设置不当
- 初始条件不恰当
- 非线性性过强（荷载步数不足）
对策:
- 进行网格质量检查（宽高比、雅可比行列式）
- 确认材料参数的单位制
- 将荷载分成多个步长（增加子步数）
- 放松收敛判定准则（但要注意精度）
也就是说，在收敛失败的地方偷工减料的话，后来会吃亏啊。铭记在心！
2. 非物理的结果
接下来是非物理结果的话题。请讲讲内容吧。
症状: 应力/位移/温度等物理上非现实的数值
可能的原因:
- 边界条件设置错误
- 单位制混杂（SI单位与工程单位混淆）
- 不恰当的单元类型选择
- 应力奇点的存在
对策:
- 确认反力的总和（力平衡）
- 确认单位制的一致性
- 重新检查单元类型的恰当性
- 奇点消除或子模型化
我前辈说过"收敛失败一定要好好做"，现在明白意思了。
3. 计算时间超过
计算时间超过具体是怎样的情况？
症状: 计算耗时远超预期
对策:
- 优化网格的粗细分布
- 利用对称性（1/2、1/4模型）
- 优化求解器设置（迭代法、预处理的选择）
- 使用并行计算
4. 内存不足
请讲解一下"内存不足"！
症状: 内存溢出错误
我前辈说过"收敛失败一定要好好做"，现在明白意思了。
对策:
- 使用核外求解法
- 减少网格规模
- 确认使用64位版求解器
- 增加内存分配
哦~，收敛失败的话题真的很有意思！ 请多讲些。
求解器别错误信息
想更详细地了解计算背后发生了什么！
| 工具名 | 开发商/现在 | 主要文件格式 |
|---|---|---|
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
| Simcenter STAR-CCM+ | Siemens Digital Industries Software | .sim, .java, .csv |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Mechanical (原ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
Nastran代表性错误
代表性错误具体是怎样的情况？
Abaqus代表性错误
请讲解一下"代表性错误"！
原来如此。那么如果工具名做成了，基本上就没问题了，对吧？
调试流程图
先生也在单体（一体型）求解法中通宵调试过吗？（笑）
1. 确认和分类错误信息
3. 用简化模型重现测试
4. 通过逐步复杂化来确定问题位置
5. 修正和重新分析
6. 确认结果的合理性
也就是说，在错误信息的确认上偷工减料的话，后来会吃亏啊。铭记在心！
质量保证检查清单
有没有教科书中没有的"现场的智慧"之类的东西啊？
嗯，做得很好！ 实际上动手操作是最好的学习方法。有不明白的地方随时可以问。
物性值不连续引发的发散
在单体方法中使用统一网格的情况下，固-液界面处的物性值（k·ρ·cp）的急剧变化会使单元内梯度异常增大，导致发散。不锈钢（k=16 W/mK）和空气（k=0.026 W/mK）的热传导率比超过600。对策是在界面单元使用调和平均数的"界面通量修正"。在OpenFOAM中，通过指定"laplacianScheme"为"Gauss linear corrected"可大幅改善数值稳定性。
相关主题
价值
一些
错误