Ansys Workbench工程应用之——结构非线性（终）：收敛方法总结

本文作为非线性系列文章的终章，也作为静力学系列文章的终章，本文内容主要摘抄自周老师的《Ansys Workbench有限元分析实例详解》。

非线性求解不收敛主要表现为4种形式：刚体运动，系统力不平衡，材料不稳定，单元匹配误差。

在求解信息中提示内部幅值超限（Internal solution magnitude limit was exceeded），或者DOF limit exceede...，Small/Negative Pivot error...，Negative main diagonal...，Max DOF INC=A very large number...等错误或警告，即表示系统有刚体平移，导致求解不收敛或无意义。在结果中查看最终（或不收敛结果的前一步）变形，会发现存在非常夸张的变形。

导致刚体运动的主要原因有约束不足、连接接触关系不足、耦合不正确（比如不同自由度的单元使用了Form New Parts组合）等。

刚体运动在不收敛问题中属于比较好检查的，可以通过模特分析找出0HZ的振形查看不足的约束，或者在结果变形中查看夸张的变形，补齐相应的约束与接触即可。常见处理方法如下。

此外，应尽量避免使用弱弹簧，因为有时弱弹簧会使刚体平移难以被快速检查出来，反而计算出不正确的结果。

在求解信息——力收敛中可查看力收敛图。紫色线表示计算的力残差，蓝色线表示力收敛容差，只有当力残差小于收敛容差时，才表示系统的力达到平衡要求，子步才完成收敛。

导致系统力不平衡的因素有：接触刚度太大导致震荡，收敛容差太小，加载太迅速，材料不稳定，结构屈曲失稳等，解决方法如下表。

材料不稳定往往伴随着其他收敛问题，所以关键需要区分导致不收敛的原因。例如，不收敛时查看力收敛结果，如果残余力非常大，一般是系统力不平衡所致；又如，模型发生高度翘曲、发生负体积进而导致不收敛，一般也不是材料的问题。

主要原因：计算塑性、超弹性、蠕变等本构时应变增量太大；选用不合适的本构模型。

处理方法：选择合适的本构模型,选择合适的单元参数，增加子步。

单元匹配误差相比前三种情况，问题出现在更深层，一般不会有明显的提示。主要表现为：模型高度翘曲、极大应变、体积自锁、沙漏现象、屈曲等。

处理方法：计算前使用Check命令检查，包括遗失的弹性特性、无约束的模型和单元形状检查，使用Mcheck命令检查网格中的孔洞和裂纹等缺陷；一般采用高阶单元计算更快更准确；对于Shell模型尽量全部采用四边形网格，可避免三角形网格计算出现的噪点误差；优先选用Direct (Sparse) Solver；尽可能增加子步；针对密封或材料成形等大应变分析使用Rezoning;对超弹性、大蠕变分析可以尝试将模型仅仅划分为一个单元。

注意，收敛完成并非万事大吉，还必须校核计算结果的准确性。可通过以下方法检查。

我看了下第一篇发的笔记，是2021年5月，现在完结的时间是2022年11月，静力学部分总共学习了1年半时间，学习版本也从WB18升级到2021R2。在这期间，认识了很多同道中人，也学习了很多知识，感谢大家的关注与陪伴。

以下文章均在仿真秀平台可查找。

静力学部分的一些坑还没有填，比如如何划分结构化网格，因为经过图惜的研究，发现结构相关的工程问题可能用不到结构化网格，普通质量的网格已经能完全应对，正如教主和一些老师所说，网格自动化划分是最有可能被软件取代的，何必花大量时间去学习呢，还有一个重要原因是WB的mesh模块功能越来越强大，操作越来越简单，真正让结构分析重点回归到工程问题，而不是数字网格上。（流体等问题另当别论，我也不懂。）

当然，图惜的总结都是一些基础性的知识，主要参考了周老师的《Ansys Workbench有限元分析实例详解》和一些其他书籍、视频、文章。很多问题并没有更深入的研究。由于图惜的知识和实践经验实在有限，所写文章也难免纰漏百出，敬请批评指正。

将来图惜还会继续和大家一起学习动力学相关内容，但是笔记风格可能改变下，毕竟长篇大论看的让人头痛。

很多朋友在问图惜的学习经历，下一期就写一篇有限元学习之路的水文吧，没有一点点理论和工程知识，也不植入广告，也不推荐书籍教程，只当朋友一样一边喝酒一边扯把子（四川话：吹牛）侃侃学习经历，敬请关注。