求解器收敛性确认与诊断
讨论收敛之前的理论基础
"收敛"一词的三种含义
教授，我汇报"分析收敛了"，结果前辈反问我："你说的是哪种收敛？"……
这个反问问得好。在CAE里，"收敛"至少指三件不同的事。①线性求解器的迭代收敛——用迭代法解 \( A\mathbf{x}=\mathbf{b} \) 时残差的下降；②非线性迭代的收敛——Newton法到达平衡状态；③离散化收敛——网格或时间步细化时解趋于极限。①②都收敛了，若没确认③，你可能只是"把错误的方程解得很准"。本文讲①②的确认与诊断；③请看网格收敛性验证。
残差的定义与收敛判据
非线性静力分析中，残差是外力与内力的不平衡量：
$$ \mathbf{r}^{(k)} = \mathbf{F}^{ext} - \mathbf{F}^{int}(\mathbf{u}^{(k)}) $$
求解器在残差范数降到阈值以下时判定"达到平衡"。实际会组合多种尺度：
| 判据 | 定义（概念） | 特点与局限 |
|---|---|---|
| 力残差判据 | \( \|\mathbf{r}\| / \|\mathbf{F}^{ext}\| \le \epsilon_F \) | 直接度量平衡，最重要。荷载接近零时相对化容易失效 |
| 位移修正判据 | \( \|\Delta\mathbf{u}^{(k)}\| / \|\Delta\mathbf{u}\| \le \epsilon_u \) | 检测解趋于稳定。无法区分停滞与收敛 |
| 能量判据 | \( |\mathbf{r}\cdot\Delta\mathbf{u}| \) 的下降 | 兼顾力与位移。绝不可单独使用 |
关键规则：绝不能只凭位移修正判据宣布收敛。迭代停滞时位移变化同样很小，力残差还很高的"表面收敛"就会蒙混过关。
什么决定收敛速度
Newton-Raphson法在解的邻域内二次收敛：残差应像 \( 10^{-2} \to 10^{-4} \to 10^{-8} \) 那样每次迭代位数翻倍。偏离这种"加速下降"模式就是诊断线索。而迭代线性求解器（CG、GMRES）的收敛速度由条件数 \( \kappa(A) \) 支配：细长单元、极端刚度比、薄壁结构都会恶化条件数、造成迭代停滞——线性求解器的表现本身就是关于网格、材料与约束的诊断信息。
控制收敛的数值方法
Newton法家族与荷载增量
非线性问题把荷载分成增量、每个增量内做Newton迭代。每次都更新切线刚度的完全Newton法二次收敛但单次迭代昂贵；省略更新的修正Newton法单次便宜、收敛降为线性。收敛困难时，务实的升级顺序是：
- 减小荷载增量——让每个增量的初始猜测落入Newton法的收敛半径
- 启用线搜索——缩放修正量，抑制向发散方向的过冲
- 加入稳定化（微小阻尼/黏性）——穿越屈曲、突跳前的不稳定；事后必须确认耗散能足够小
- 切换到弧长法（Riks法）——荷载-位移曲线折返（snap-through/snap-back）时荷载控制本身失效，用弧长参数追踪路径
线性迭代法与预条件
大规模模型中直接法（LU、Cholesky）受内存制约，预条件迭代法成为标准。预条件本质上是降低有效条件数的变换，代表是ILU系与代数多重网格（AMG）。迭代数百次仍停滞时，正道是更换预条件或改善网格质量，而不是提高迭代上限。求解器与矩阵的匹配也要确认：对称正定用CG，非对称（流体、部分接触格式）用GMRES。
收敛历史的读法
残差曲线只要在往下走就行了吧？
"怎么往下走"里才藏着信息。学会按模式解读，对策就能反查出来。
| 残差历史模式 | 解读 | 典型对策 |
|---|---|---|
| 每次迭代位数翻倍地下降 | 健康的二次收敛 | 保持现状即可 |
| 恒定速率的直线下降 | 修正Newton/拟Newton的线性收敛 | 可接受则无妨；太慢就提高切线更新频率 |
| 中途转平（停滞） | 切线与真实刚度偏离、接触抖动、病态矩阵 | 缩小增量、接触稳定化、检查预条件 |
| 振荡不下降 | 接触开合反复、材料加卸载往复 | 消除初始穿透、加阻尼、缩小增量 |
| 单调增大（发散） | 超出收敛半径、结构失稳（屈曲）、输入错误 | 初始增量取1/10、稳定化、重查模型 |
CFD收敛判定的特殊性
稳态CFD中，"残差降到规定值"与"答案定下来了"之间的偏差尤其大，必须残差与监测量并用。残差停在 \( 10^{-3} \) 但升力系数、压损完全稳定，实用上可能已收敛；反之残差很小、监测量却缓慢漂移（伪稳态）也常见。判定用三件套：①标度残差的下降，②积分量（力、流量、平均温度）达到定值，③边界守恒（质量·能量）残差小于输入的1%。分离流中监测量周期振荡，说明稳态解本就不存在——该转非定常分析了。
实务中的确认流程
不要盲信默认收敛判据
各求解器的默认值（如力残差0.5%）是"对多数问题都还行"的折中，不保证对你的问题足够。验证的定式是判据敏感性检查：把判据收紧一个量级重算，确认评估量（最大应力、位移、反力）的变化小于精度要求。变化大，说明常规判据太松；同样的方法也能发现过度收紧、白白浪费机时的情况。
把收敛历史当作交付物保存
收敛历史（每个增量的迭代数、cutback发生位置、最终残差）是分析的"体检报告"。报告里只贴最大应力云图、把历史扔掉的做法，事后无法审计结果的可信度。至少记录：①所有增量是否都在判据内结束（有无被截断的增量），②cutback集中在哪个荷载区间（那里有物理事件），③最终增量的残差值。只要有一行"达到最大迭代数，进入下一增量"，那个解就没有处于平衡。
切分收敛失败原因的标准流程
调参之前，先用二分法锁定原因：
- 退化为线性弹性、小变形——这样都解不了，问题在约束、连接或网格质量，与非线性无关
- 逐项恢复非线性——按材料、几何、接触的顺序启用，找出加入哪一项时崩坏
- 荷载减到1/10——小荷载能收敛⇒增量控制问题；仍不能⇒模型定义问题
- 可视化崩坏前一增量的变形与应力——局部屈曲、单元翻转、接触咬入都会现形
按此顺序，"不收敛"的九成能在几个增量的排查内定位。
各求解器的收敛控制与诊断输出
判据设置与诊断文件对照
| 求解器 | 收敛判据设置 | 诊断信息位置 | 实务备忘 |
|---|---|---|---|
| MSC Nastran / NX Nastran | NLPARM（CONV=UPW选择位移/荷载/功判据） | .f06的非线性迭代汇总（EUI·EPI·EWI） | 判据组合影响结果，别丢掉荷载判据 |
| Abaqus | 默认力残差0.5%（时间平均力之比），*CONTROLS修改 | .msg（迭代明细）·.sta（增量汇总） | 在.msg里确认cutback位置与奇异性警告 |
| Ansys Mechanical | CNVTOL（力·力矩·位移） | Solution Information·收敛曲线 | bisection发生点有物理事件；注意自动弱弹簧 |
| Ansys Fluent | 标度残差目标＋监测器 | 残差图·报告定义 | 不要只看残差；用Report Definition常态监测积分量 |
| OpenFOAM | fvSolution的solver tolerance／SIMPLE的residualControl | 日志＋solverInfo/residuals functionObject | 把松弛系数(URF)与收敛速度的取舍记录在案 |
Cutback（增量自动缩减）的读法
隐式求解器遇到收敛失败会自动缩小增量重试（Abaqus的cutback、Ansys的bisection）。这不是故障而是诊断信息：cutback集中的荷载水平上存在物理事件——屈曲、接触状态切换、塑性区扩展。把那个荷载区间的变形与接触状态可视化，模型想告诉你的行为就能读出来。把cutback当"浪费时间的麻烦"、粗暴放宽最小增量，等于扔掉诊断信息。
注意：显式求解没有"收敛判定"——盯能量
显式动力学（LS-DYNA、Abaqus/Explicit）不解联立方程，因此根本不存在收敛判据。替代的健全性指标是能量平衡：总能量守恒（内能＋动能＋耗散−外力功≈常数）、人工能（沙漏控制）低于内能的百分之几、质量缩放的质量增加率受控。"显式所以不用担心收敛"是误解——只是要审计的量换了而已。
前沿研究动态
大规模化与可扩展的预条件
千万自由度级分析的研究核心，是收敛性与可扩展性兼得的预条件：代数多重网格（AMG）在理想条件下实现与网格尺寸无关的迭代数，区域分解法（FETI-DP/BDDC）正在成为超大规模并行的标准。对实务的启示：许多"大规模收敛慢"的问题，通过求解器与预条件选型可获得数量级改善——在按默认设置扩充硬件之前，值得先检讨这一层。
机器学习加速收敛
用神经网络学习良好初值或预条件、削减迭代次数的研究很活跃。从相似分析结果生成初始猜测的方法，在"大量求解相似问题"的工作流（形状优化、参数研究）中已有可观的削减率报告。但分布外问题可能适得其反，因此当前健全的用法是：加速可以交给学习，收敛与否仍由传统残差判据裁决。
非线性鲁棒性的自动化
自动调节增量控制、稳定化与求解器切换的自适应策略在各求解器中不断完善：拟Newton（BFGS系）＋自动线搜索、有全局收敛保证的信赖域方法、接触的内点法格式等，都在压缩"靠工匠手感调参"的空间。而自动化程度越高，核查日志、确认自动做了什么的验证者角色就越重要——像耗散能事后审计这类对自动稳定化副作用的定量确认，今后只会更关键。
故障排查
按症状的诊断表
| 症状 | 可能原因 | 对策 |
|---|---|---|
| 第一个增量就发散 | 约束不足·单位制错误·初始穿透·荷载数量级错误 | 退化为线性分析排查；接触先做初始调整消除穿透 |
| 每次都在同一荷载水平停滞 | 屈曲·突跳·接触状态大规模切换 | 可视化其前一刻的变形；考虑弧长法、稳定化或转动力学 |
| cutback连发、推进不了 | 最小增量对该非线性来说仍太大 | 降低最小增量；同时排查局部网格与材料曲线拐点 |
| 线性求解器迭代停滞 | 病态（极端单元形状·刚度比·薄板） | 换AMG等预条件、改善网格，用直接法交叉验证 |
| 残差在降但结果随判据变化 | 判据太松（表面收敛） | 收紧一个量级做敏感性确认；以力残差为主判据 |
| CFD残差在10⁻³横盘 | 伪稳态·非定常性·网格引起的振荡 | 用监测量与守恒判定；振荡则转非定常 |
怎么都不收敛时的最终确认
全都试过了还是不收敛的话，该怎么办……？
最后要怀疑的是"这个问题究竟存不存在静力平衡解"。超过崩塌荷载的结构、动态突跳的屈曲、自激振荡的流动——让求解器去找一个不存在的稳态解，它当然永远收敛不了。这时"不收敛"本身就是物理给出的答案。唯一正确的处理，是换成追踪时间或路径的格式：动力学分析、弧长法、非定常分析。分清"求解器设置问题"还是"物理问题"，是收敛诊断的最后一关。
收敛确认不是"计算跑完了没有"，而是"平衡与守恒满足到了什么程度"的定量确认。把判据敏感性检查和历史存档变成习惯，收敛引发的质量事故基本可以杜绝。另请参阅边界条件的验证方法、网格收敛性验证。
帮助
更多
错误