过度迭代（收敛延迟）

不只是变慢而已。例如在非线性解析中，收敛判定值（比如残差范数）要达到目标值1e-3，需要远超过默认最大迭代次数（Abaqus/Standard默认为16次），达到50次、100次等反复迭代的状态。在物理上讲，求解器正在「四处搜索」平衡点，计算成本变得巨大。

有三个主要原因。第一是「材料非线性的急剧变化」。例如，降伏应力为200MPa的钢材，某个单元从199MPa到201MPa时，刚度矩阵会急剧变化。第二是「接触条件的不连续变化」。第三是「不稳定物理现象本身」。在座屈或snap-through发生的区域，解可能不唯一或不稳定。

完全正确。Newton-Raphson法的收敛性在很大程度上取决于该切线刚度的精度。每当材料进入塑性区、接触开启/关闭时，

例如在Abaqus/Standard的Full Newton法中，每增分内重复以下循环：1) 形成和分解切线刚度矩阵（成本大），2) 计算残差力，3) 求解位移增量，4) 收敛判定。过度迭代是指该循环无法满足「残差力为平均力的0.5%以下」等收敛条件而无限循环的状态。特别是第1)步的矩阵分解（如LDL^T）在每次迭代中进行，计算时间会随迭代次数成比例爆炸。

有。那就是「修正Newton法」。在Abaqus中，可通过`STEP`选项指定`ITERATION=MODIFIED`。这种方法只在第一次迭代时计算和分解切线刚度，之后重用同一矩阵，仅更新残差力。收敛会变慢（迭代次数增加），但单次迭代的成本大幅下降，总计算时间可能减少。不过在强非线性问题中易发散，存在权衡。

影响很大。Ansys Mechanical的默认直接法求解器（Sparse）对刚度矩阵的病态性相对鲁棒，但大规模模型会大量消耗内存。迭代法求解器（如PCG）内存效率更高，但条件数恶化时收敛会受阻，与非线性迭代合并后形成「双重迭代地狱」。网格质量差（宽深比>10）时，迭代法求解器的收敛极端恶化，全体迭代数增加，这是典型案例。

首先打开日志文件（Abaqus的.msg、Ansys的.out），确定在哪个增分（或时间步）迭代次数急剧增加。接下来，绘制该增分的「力残差」和「位移修正量」历史。在Ansys中，在`SOLVE`命令前后设置`OUTPR,ITER`。观察残差是在特定模式（如局部座屈模式）上振荡，还是在某个值以下无法进一步下降。

可视化该增分开始时的模型状态。在Abaqus中使用`restart`，在Ansys中用`RESUME`读取该步，检查应力、塑性应变、接触状态。例如，某一接触对的接触压力是否在10MPa和0之间振荡？或某个单元的塑性应变超过0.2后问题才开始？这是「物理原因」确认工作。

基本方法是「减小增分大小」。在Abaqus中设置`INITIAL=1E-4, MINIMUM=1E-8`等极小值，让求解器缓慢求解问题。其次是「放松收敛判定」。将默认的残差基准0.5%改为1%或2%（精度和收敛的权衡）。在Ansys中通过`CNVTOL`命令调整。然后，「对不稳定现象导入弧长法」。Abaqus的`STATIC, RIKS`或Ansys的`ARCLEN`命令，这些对追踪解路径很有效，但设置复杂。

检查以下5点：1) 网格：宽深比是否≤10，特别是在塑性区和接触面？2) 材料数据：塑性曲线有无不连续折点（直接输入实验数据吗？）。3) 接触：初始接触状态是否恰当，罚函数刚度是否太大（默认100倍等）？4) 约束：是否存在剩余刚体模式？5) 单位制：力[N]和长度[mm]混用，刚度是否被实质乘以1e6？这些都会恶化刚度矩阵的条件数。

基本思路相同，但设置方法和默认行为不同。Abaqus/Standard默认收敛判定相对严格，最大迭代次数也少（16次），容易因「无法收敛」而停止。而Ansys Mechanical的默认最大迭代次数较多，会进入「过度迭代」状态继续计算，忽然发现已过了很长时间。

有几个。首先是`*CONTROLS`参数，可调整刚度矩阵的更新频率和收敛判定。例如，`ANALYSIS=DISCONTINUOUS`对接触或摩擦导致刚度不连续变化的情况有效。还可用`*VISCO`引入伪阻尼，抑制振荡式收敛延迟。切换到Explicit求解器（`DYNAMIC, EXPLICIT`）在接触极其复杂时是最后手段，但计算时间就成了另一个问题。

非常有效的功能之一。这个功能在迭代中判断收敛缓慢的区域（如应变集中部），在该增分内自动细化网格。R19及以后版本本格导入。如果是由局部变形集中引起的刚度恶化，仅在该部分提高分辨率可改善切线刚度精度，减少全体迭代数。但要权衡自动重网格的计算开销。

是这样。COMSOL需求解耦合的多物理场，不同物理量的尺度差异大。例如热应力分析，温度场变化（数量级10^2 K）与位移场（数量级10^-3 m）耦合。默认缩放不当时，Jacobian矩阵条件数会极端恶化，成为过度迭代的主因。在COMSOL中必须设置「位移缩放」或「手动缩放因子」，使所有自由度的大小相当，这是关键对策。

这是「收敛速率为线性」状态，Newton法本应有的二阶收敛性丧失。很可能是切线刚度更新功能失效。在Abaqus中尝试将`STEP`的`EXTRAPOLATION`设为`OFF`。初始预测可能偏离太远。接下来，在`CONTROLS`中开启`LINE SEARCH`或调整系数（将默认0.25增大到0.5等）。线搜索对位移增量

这表示不稳定的物理现象或接触「打颤」（chattering）在发生。先确定振荡的自由度。用Ansys的`PRNSOL`或Abaqus的字段输出，看每次迭代的位移。如果是局部座屈等模式振荡，认真考虑导入弧长法（`ARCLEN`）。接触打颤时，用Ansys的`CNCHECK`命令详查接触状态，降低罚函数刚度（默认值的0.1倍等），或改用「Boltzmann增分」等平滑接触模型。

很可能是跨越了材料降伏或损伤启动等「状态变化阈值」。例如，某增分模型1%的单元进入塑性，下一增分变成10%。这种急剧状态变化让求解器无法应对。对策分两个阶段。首先，从接近该阈值的步开始，手动将增分大小减小10倍或100倍。其次，重新审视材料模型。例如降伏后的塑性硬化曲线斜率（硬化系数）过小，数值上可能病态（ill-conditioned）。即使是JIS G 3131 SPHC这样的材料，降伏后的数据点也应至少定义到0.2% Proof Stress的数倍应变。

现实中有3个「逃脱方案」。1) **求解器切换**：从Abaqus/Standard改用Abaqus/Explicit（动力问题时）。或从Ansys Mechanical改用Ansys LS-DYNA。2) **物理模型简化**：先暂时移除收敛差的非线性性（如摩擦系数0.2的摩擦），去掉它后运行，从无摩擦情况开始理解现象机制，再重新挑战。3) **人工稳定化**：用Ansys的`STABILIZE`命令或Abaqus的`*VISCO`加入虚拟阻尼力，强制数值收敛。这仅为「让其收敛」的手法，结果的物理正当性需另外验证。作为工程师，此时需判断是放弃还是从基础重建模型。