ANSYS非线性求解过程中迭代曲线的各项含义

对于在ANSYS中进行过非线性求解的读者，下面这张图片肯定不会陌生。从开始求解到求解结束，这个界面会一直显示在求解界面中。下面就来介绍下这个图片中各个项目的含义。

1、坐标轴含义

（1）横轴：Cumulative Iteration Number

Cumulative Iteration Number 翻译为累积迭代数，即求解器在整个非线性求解计算过程中所进行迭代计算的次数，当一个问题的非线性越强，其数目就会越多，故而横坐标数量的大小直接与项目的非线性程度相关。

（2）纵轴： Absolute Convergence Norm

Absolute Convergence Norm 翻译为绝对收敛范数。我们在建模过程中输入的各个参数都不是范数的形式，而ANSYS在求解过程中，会将相应的变量进行类似规范化处理，也即是这儿所体现的范数，很多情况来讲，我们并不需要关心绝对收敛范数的具体值。

2、时间

上部的Time 相对来讲是一个比较重要、需要我们关心的数值。在非线性求解过程中，time只是一个时间标记，与真实时间无关。一般来讲，我们在非线性求解之前应设置好荷载步、荷载子步，具体求解时，time会按照我们设置的时间步进行相应的变化。如果在具体求解过程中，发生不收敛情况导致求解退出，我们可以根据时间数据定位到发生求解失败的地方从而方便我们调整模型。

3、曲线

曲线的走向是我们在进行非线性求解过程中很重视的一个点，这里我们比较常用的有三个收敛准则：力收敛准则F、位移收敛准则U、弯矩收敛准则M。对于收敛准则的选取和设置方法可参考本站文章《提高ANSYS非线性求解收敛性能的常用方法》。

在曲线的右侧有CRIT和L2标签，其中CRIT是 criteria的缩写，代表收敛判断准则。L2指的是L2级范数，也即是通常所称的计算残差，除了L2，还有L0、L1，具体可以参考命令Cnvtol的帮助。

求解器在每一步迭代计算过程中，会进行相应的收敛判别准则和计算残差，计算残差是所有单元内力的范数，只有当残差小于准则时（也即L2<CRIT），本荷载步计算收敛，并进行下一步的迭代计算。 故在计算过程中，我们会注意到L2有许多的峰值，并且L2会逐渐下降，向CRIT靠拢并最终交汇小于CRIT，本次迭代计算收敛，进行下一步。

那么如何解释CRIT曲线不断在缓慢上升呢？个人认为可以这么理解：假定力收敛准则为设置为2000，荷载为10，则对于第一个荷载步而言，力收敛的绝对值就是2000/10=200，是总体力收敛值的1/10，在随后计算累积的过程中，逐渐向最终值2000发展，当然，这只是个人理解。

其实迭代曲线为我们提供了可以监视求解器运算状况的窗口，如果在整个计算过程中发现曲线不正常，无法收敛，这时候就需要对我们的模型或者非线性求解设置进行调整了。