ANSYS Mechanical非线性接触功能基础概述，了解一下下吗？

导读：ANSYS Mechanical具有强大的非线性计算能力，能够对几何非线性、材料非线性、接触非线性、混合非线性（塑性和接触等）、非线性屈曲等计算问题进行非常好的模拟仿真，是目前最强大的非线性问题计算软件之一。 

限于篇幅，本文仅对非线性求解与收敛、接触非线性问题中涉及的设置等内容进行简要整理和说明，以期望对于ANSYS Mechanical接触非线性设置的使用有一定的帮助作用。 
 

本文原稿基于ANSYS Documentation-ANSYS Mechanical User's Guide中部分章节内容以及其他相关技术文档进行整理编写，另外原稿整理时间较早，部分设置与ANSYS18/19版本不尽相同。 

本文共  6737 个字，预计用时超过 10分钟 ，您将学到以下内容： 

1、非线性求解与收敛 
 

-收敛与收敛判据 

-载荷步，时间步与平衡迭代 

-求解控制 

-重启动控制 

-非线性控制 
 

2、接触与接触设置 

-接触协调 

-探测方法 

-修剪接触 

-穿透和华移容差 

-法相接触刚度 

-Pinball 区域 

-对称/非对称行为 

-接触中的体类型 

-界面处理与接触几何修正 

-接触工具 

-自接触设置 

-单元死生与接触过程控制 

一、非线性分析背景 

1、结构非线性的定义 

江边渔者垂钓的鱼竿、承载运动器重量的轮胎、重型机械零件的锻造、薄板零件的冲压与折弯等，都涉及结构非线性问题。 

对这些结构绘制载荷与变形的曲线，均不再符合胡克定理 ，而表现出非线性结构的基本特性，结构刚度不再是常量，而成为函数变量。 

2、非线性行为类型 

① 几何非线性 

当结构经受大变形时，变化的几何形状可能会引起结构的非线性响应。例如鱼竿提线的过程、大型壳体结构非线性屈曲的过程，都属于几何非线性问题。 

② 材料非线性 

金属塑性、橡胶材料超弹体与粘弹性、混凝土问题、率相关性蠕变等问题都是典型的材料非线性问题。 

图1 混合非线性 

③ 接触非线性 

接触是种普遍的非线性行为，是状态变化非线性类型**殊而重要的子集。当零件彼此之间接触与分离，接触刚度会发生相应的改变。 

非线性问题通常不会单一出现，ANSYS Mechanical具有同时处理混合非线性问题的能力，如图1所示就是混合非线性的问题，同时考虑金属的塑性与接触的非线性问题。 
 

3、构建非线性模型 

适度的非线性行为求解采用的几何设置、网格通常无需修正，与线性行为一致。 

多数情况需要对若干特性进行修正，包括单元的特性、非线性材料数据（塑性、蠕变、超弹体等）、几何建模（奇异引起不收敛问题）等。 

如下设置与修改，可能有益于非线性分析问题的收敛。 

① 修改【Mesh】→【Shape Checking】选项为“Aggressive Mechanical”或“Nonlinear Mechanical”，如图2所示，这能够提供一种增强的单元质量处理大应变分析中过度扭曲。 

② 修改【Geometry】→【Element Control】选项为“Manual”，可以在完全积分与减缩积分策略中进行切换。该选项影响单元中积分点的数量，如图3所示。例如强制修改完全积分应用于高阶单元（默认采用一致减缩积分），这通常有助于厚度上仅有一层单元的结构，以提高计算求解精度。 

③ 默认情况下，结构单元将默认采用具有中间节点的高阶单元。通过图4中所示进行高阶单元（Kept）与低阶单元（Dropped）的转换。 

图2 Shape Checking 

图3 Element Control 
 

图4 Element Midside Nodes 

二、非线性求解与收敛 

（一）牛顿-辛普森方程 

 

非线性分析中，不能直接由线性方程求得响应，需要将载荷分解成许多增量求解，每一增量确定一平衡条件。 

 

ANSYS Mechanical 使用牛顿-辛普森方程（Newton-Raphson）平衡迭代法进行非线性问题的求解。每个载荷增量步结束时，平衡迭代驱使解回到平衡状态。 

 

其中： 

 

为切向刚度矩阵； 

为位移增量； 

是施加的载荷矢量； 

为内力矢量。 

 图5 牛顿-辛普森方程 

 

如图5所示是一个载荷增量四个迭代步的迭代求解过程： 

 

第一次迭代施加总载荷 

 ，对应的位移结果为 ，根据位移 ，计算力，若是 ，系统不收敛，将进行刚度矩阵的修正，如图中点虚线所示。其中 的差值即外力与内力的偏差成为残差力，残差力需要足够小才能得以收敛，然后进行第二次迭代求解，第三次迭代……直至 收敛。 

 

（二）收敛与收敛判据 

 

牛顿-辛普森方程需要一个收敛度量以决定迭代的过程，定义“残差”是外部载荷与内部载荷的差值： 

 

 

通常认为当残差足够小的时候 

 

 

即当残差的范数小于指定容差因子乘以参考力值时，就认为得到了收敛的求解数值（实际上残差从未真正等于0）。 

 

其中： 

 

 称为力收敛准则。 

 为容差因子。 

是载荷与反力的范数。 

 

（三）载荷步，时间步与平衡迭代 

 

1、载荷步，时间步与平衡迭代 

 

非线性求解过程可以按照三个层次进行组织： 

 

① 载荷步是顶层，求解选项、载荷与边界条件都施加于某个载荷步内。 

② 时间步是载荷步中的载荷增量，时间步用于逐步施加载荷。 

③ 平衡迭代是为得到给定时间步（载荷增量）的收敛解而采用的方法。 

 

图6 两个载荷步的求解 

 

率相关分析（蠕变，粘塑性）与瞬态分析中，“时间”代表真实的时间。 

率无关的静态分析，“时间”表示加载次序。 

 

2、自动时间步 

 

“自动时间步”基于结构对于应用载荷的响应，在每个子步结束将会计算一个优化的时间步长。 

 

 “自动时间步”通过求解过程自动调节载荷增量，可以定义初始载荷步、最小载荷步和最大载荷步。 

 

难收敛的求解问题，将会施加更小的增量。当收敛困难时，自动时间载荷步算法将进行“二分”，求解将返回最后一个成功收敛的子载荷步，按照更小增量施加载荷(通常设置为上次的一半)。 

 

（四）求解控制 

 

 图7 求解控制 
 

结构静力学【Static Structural】和结构瞬态动力学【Transient Structural】分析模块均可以考虑非线性分析问题。 

1、【Solver Type】求解类型 

(1) Program Controlled程序控制能够根据求解问题自动选择求解器。 

(2) Direct (Sparse)稀疏求解器更稳健，推荐用于梁壳等结构中。 

(3) Iterative (PCG)更有效率，适合用于大型结构求解。 

2、【Weak Springs】弱弹簧 

用来防止结构刚**移导致的计算不收敛，引入Weak Springs能帮助求解收敛。由于刚度很小，不会明显影响结构的计算结果。 
 

3、【Large Deflection】大变形 
 

① 对于结构瞬态动力学，默认设置为Large Deflectio=On。

② 对于静力学分析模块，需要修改Large Deflection=On。 

Large Deflectio=On，使分析考虑大变形、大旋转和大应变引起的单元形状和方向的改变，计算结果更为准确。 

（五）重启动控制 

重启动技术用于已经停止计算的问题再次提交求解。能够对不收敛的计算进行修改，例如分析设置、载荷等，也能对已经完成求解的模型进行扩展分析，例如增加新载荷步、求解时间等。 

如图8所示，重启动设置的过程一般如下： 

1、修改Generate Restart Points=Manual。 

2、 修改Load Step，设置为Last或者All 

3、修改Substep选项用于指定“Load Step”内生成多少个重启动点。可选择“Last、All、Specified、Equally Spaced”等方法进行定义。 

4、默认情况下，重启动文件是自动删除的，需要修改Retain Files After Full Solve=Yes，设置Future Analysis=Prestressed analysis或Delete Unneeded Files=No，都可以保留重启动文件。
 

5、当求解计算不收敛时，通过图8步骤的方式，建立重启动点的设置，再次提交求解计算。 

 

图8 重启动控制 

（六）非线性控制 

 

点击完整阅读全文

ANSYS Mechanical非线性接触功能基础概述，了解一下下吗？

相关阅读推荐：​​​​​​​

ANSYS Workbench2022R1有限元理论与工程实践方法120讲（免费）