Newton-Raphson法失败 — 故障排除
Newton-Raphson法失败 — 故障排除
故障排除步骤
当Newton-Raphson法不收敛时,请整理应该检查的内容。
步骤1:收敛监控
确认以下信息:
| 确认项目 | Abaqus | Nastran | Ansys |
|---|---|---|---|
| 力的残差 | .sta FORCE | .f06 OLOAD | Force Convergence |
| 位移增分 | .sta DISP | .f06 DISP | DOF Increment |
| 接触变化 | .sta SDI | .f06 CONTACT | Contact Status |
| 能量 | .sta ENERGY | .f06 ENERGY | Energy Convergence |
步骤2:收敛基准的确认
默认收敛基准:
- Abaqus: 力的残差 < 0.5% × 时刻历最大力、位移修正量 < 1% × 总增分位移
- Nastran: 力的残差 < 10^-3 × 外力范数(CONV=PW)
- Ansys: 力的残差 < 0.5% × 参考值(L2范数)
如果放松基准会收敛,但精度会降低,对吧?
正是如此。如果将力的残差基准从0.5%放松到5%,收敛会更容易,但结果精度会降低。按经验,如果要求应力精度±1%,残差基准应为0.5%以下;如果为±5%,则2%左右就足够了。
步骤3:荷载步长策略
根据非线性程度的增分设置:
步骤4:替代算法
当标准Newton-Raphson法无法求解时:
1. 线搜索: 抑制过冲。基本上无风险可激活
2. 拟Newton法(BFGS): 减少刚性更新成本。但不适合接触
3. Arc-Length法: 对座屈、跳跃问题必需
4. 显式法: 不需要Newton迭代。对动态问题或剧烈接触最优
5. 接触阻尼: 用粘性稳定接触不稳定性
所以最后还是要根据具体问题选择算法。
正确。没有万能方法。但在实践中,首先启用线搜索并减小增分。如果不行,再用接触稳定化。如果仍不行,再用显式法。按这个顺序升级最有效率。
花在求解器错误原因判断上的时间,应该可以更少。 — Project NovaSolver将错误诊断体验改善作为研究课题之一。
与实务人员共同思考CAE的未来
Project NovaSolver直面Newton-Raphson法失败中的实务课题本质,致力于开发支撑工程现场的工具。
查看项目最新信息 →相关主题
本文的评价
感谢您的回答!
有参考
价值
价值
希望
更详细
更详细
报告
错误
错误