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）的转换。

二、非线性求解与收敛

（一）牛顿-辛普森方程

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

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

其中：

为切向刚度矩阵；

为位移增量；

是施加的载荷矢量；

为内力矢量。

 图5 牛顿-辛普森方程

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

第一次迭代施加总载荷 

，对应的位移结果为

，根据位移

 ，计算力

，若是

，系统不收敛，将进行刚度矩阵的修正，如图中点虚线所示。其中

的差值即外力与内力的偏差成为残差力，残差力需要足够小才能得以收敛，然后进行第二次迭代求解，第三次迭代……直至

收敛。

（二）收敛与收敛判据

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

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

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

其中：

称为力收敛准则。

为容差因子。

是载荷与反力的范数。

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

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

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

① 载荷步是顶层，求解选项、载荷与边界条件都施加于某个载荷步内。
② 时间步是载荷步中的载荷增量，时间步用于逐步施加载荷。
③ 平衡迭代是为得到给定时间步（载荷增量）的收敛解而采用的方法。

图6 两个载荷步的求解

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

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

2、自动时间步

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

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

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

（四）求解控制

结构静力学【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 重启动控制

（六）非线性控制

1、【Newton-Raphson Option】：略。

2、 默认的收敛准则包括：

①【Force Convergence】

② 【Moment Convergence】

③ 【Displacement Convergence】
④ 【Rotation Convergence】

默认的收敛准则适应于大多数的工程应用，可以重新修改默认收敛准则的容差，使其更加容易收敛，但这也会影响计算结果的精确程度。对于力 / 力矩缺省的容限是0.5%。对于位移 / 旋转增量的容限是 5%。

3、【Line Search】线性搜索：线性搜索有助于收敛，通过一个0-1的比例因子去影响位移增量帮助收敛，适合于施加力载荷、薄壳、细长杆结构或求解收敛震荡的情况。

4、【Stabilization】稳定性：用于非线性屈曲分析，限于本文宗旨和篇幅，不进行展开说明。

三、接触与接触设置

​ 点击完整阅读全文
ANSYS Mechanical非线性接触功能基础概述，了解一下下吗？
​