Ansys状态非线性分析一一接触

1. 基本概念

1）接触状态

        接触状态分为远离（Far)、接近（Near）、黏接（Sticking）和滑动(Sliding）。接触状态在数值分析时可能出现两种状态：穿透（Penetration）和间隙（Gap）。物体在接触过程中是不允许穿透的，有时过大的穿透也会导致接触脱离，如果发生接触脱离（状态突变）时，会出现阶跃函数，导致收敛困难：通过设置可以允许一些轻微的穿透，使接触状态没有发生突变，则较容易收敛。

2）接触面与目标面

        接触面(Contact Surface）在接触时不能穿透到目标面 (Target Surface），目标面则可以穿透到接触面之中（详见Ansys帮助手册）。与Abaqus一致的是，可以将节点组（命名选择）定义为接触面（生成CONTA177单元），但是目标面（生成TARGE170）只能是面单元类型。对于刚-柔接触分析，明显应该把刚体定义为目标面，柔性体定义为接触面。如果柔-柔接触分析，接触面与目标面的一般有利规则为：

①凸面对凹面/平面，凸面定义为接触面；

②精细网格对粗糙网格，精细网格定义为接触面；

③刚度不同的接触对，软的定义为接触面；

④高阶单元对低阶单元，高阶单元面定义为接触面；

⑤一个面大于另外一个面，小面定义为接触面。

3）接触类型（接触行为）

        接触行为分为线性接触和非线性接触。线性接触包括绑定（Bonded）、不分离 (No Separation）；非线性接触包括无摩擦 （Frictionless)、粗糙(Rough）、摩擦（Frictional）等。这几种用接触状态表现如下：

-绑定 （Bonded）为不穿透、不分离、不滑移；

-不分离 （No Separation）表现为不穿透、不分离、无摩擦的小滑移；

-无摩擦（Frictionless ）为不穿透，自由的分离和滑移；

-粗糙（Rough）为不穿透、可分离、不滑移；

-摩擦（Friotional） 为不穿透、可分离、带摩擦的小滑移。

2. 接触定义与高级设置

1）（接触行为）公式化

        ansys提供以下几种公式模式：

①Pure Penalty（纯罚函数）
Pure Penalty 计算原理：在一对接触单元上添加接触刚度矩阵Knormal（K是非线性刚度），法向和切向力应用线弹性F=K.x的方法进行计算。

②Normal Lagrange（法向拉格朗日）
这种方法通过增加一个额外的自由度（接触压力）来满足接触兼容性。因此，接触力（接触压力）是作为额外的自由度明确求解的。Normal Lagrange 计算原理：F=接触压力（DOF）。强制实施零/近乎零的穿透与压力自由度（DOF）；不需要接触刚度（零弹性滑移）；需要直接求解器（Direct Solver）。

③Augmented Lagrange(增广拉格朗日）——系统默认值
Auemented Lagrange 计算原理：F=K.x+C。在罚函数法的基础上加一个额外项C，相较纯罚函数法，并不需要很大的法向刚度就可以保证较小的穿透，所以易收敛，但是如果网格变形得过于扭曲，则计算需要迭代的步数较多。

④Multi-Point Contraint (MPC）
MPC计算原理：内部通过约束方程 (Constraint Equations）以连接(Tie）接触面之间的位置，这种方法直接有效地连接接触区域，而且可以适用于大变形，主要用于Bonded 和 No Separation。特别适用于处理 Solid 与 Shell、Shell 与 Shell 接触时易出现的接触面法向错误，例如，提示为 "The normal of contact element XXX is not consistent with the normal of contact element XXX. Please use the ENORM command to correct it.” 。
另外，MPC 算法是整个接触算法中求解速度最慢的，当 Solid 与 Shell 模型需要定义接触时，可以采用下面方法提高计算速度：在Solid 的表面手动新建一层非常薄的 Shell模型，Solid 模型与新建的 Shell 模型采用 Bonded 连接，新建的 Shell 模型与原来的 Shell 模型采用Mesh Connections 连接。

⑤Beam （梁）
接触物体之间采用无质量的梁绑定连接，由于没有生成单元，所以不能使用 Beam Tool 查看结果，但是可以用 Probe- Bcam工具得到相应结果。一般不建议使用，如有需要，请查看专栏另一篇文章《Ansys Workbench螺栓模拟方法》，进行更合理的梁单元连接。

2）接触行为

        Mechanical提供三种行为：对称、不对称、自动不对称。

①对称：接触对的两个面都可以既是接触面(Contact Surface）也是目标面 (Target Surface），两个面的接触参数（网格粗细、接触刚度、面积、曲率等）接近时，接触容易建立时使用，可以同时提取接触面和目标面的结果。计算时间较长，在复杂情况下不易收敛。不适用与Normal-Lagrange和MPC公式。

②不对称：设置框中手动定义的接触面和目标面生效，在后处理的接触状态结果中，仅接触面的状态结果会被给出。翻转接触-目标面可能会导致问题，需要用户对接触-目标对有清晰认知。

③自动不对称：系统根据算法自动定义接触面和目标面，系统默认。

3）修剪接触

        当“Trim Contact”设置为程序控制时，对于自动接触修剪接触默认启用，对于手动接触则默认不进行修剪。对于自动接触，修剪接触半径会是自动值且不可编辑；对于手动接触，需要用户输入大于零的修剪半径。

        修剪接触实际上是减少了接触过程中系统检索的积分点数量，在小滑移的假设前提下，应用修剪接触一般不会对计算过程和结果造成影响，反而会加快计算速度；然而在有限滑移的情况下，应用修剪接触很可能会导致计算不收敛。

4）检测方法

①程序控制（Program Controlled）：
        这是默认设置。当公式设置为纯罚函数（Pure Penalty）和增强拉格朗日（Augmented Lagrange）时，应用程序使用高斯积分点（On Gauss Point）。对于MPC（多点约束）和法向拉格朗日（Normal Lagrange）公式，它使用节点点（Nodal-Normal to Target）。
②高斯点（On Gauss Point）：
        接触检测位置位于高斯积分点上。这个选项不适用于使用MPC或法向拉格朗日公式的接触。
③节点 - 从接触面法线出发（Nodal - Normal From Contact）：
        接触检测位置位于节点上，其中接触法线与接触面垂直。这意味着接触检测是基于接触表面的法线方向进行的。
④节点 - 目标面法线（Nodal - Normal To Target）：
        接触检测位置位于节点上，其中接触法线与目标面垂直。
⑤节点 - 接触面法线投影（Nodal - Projected Normal From Contact）：
        接触检测位置位于接触和目标表面重叠区域（基于投影的方法）。这种方法使用投影来确定接触点，通常用于处理复杂的几何形状或接触条件。

5）渗透容差

        渗透容差控制了接触计算的收敛精度，影响了计算速度。容差值越大，计算收敛越快，但计算精度越差，法向接触应力值误差越大，反之收敛越慢甚至可能无法收敛。允许在公式化中选择“程序控制（Program Controlled）”、“纯罚函数（Pure Penalty）”或“增强拉格朗日（Augmented Lagrange）”时为接触对指定穿透容忍度值或穿透容忍度因子。穿透容忍度因子最大值1.0，基于用于计算穿透值的第一层单元厚度。

6）法向刚度

        接触刚度的因子范围一般是0.01~10，10表示绑定/无分离；数值越小表示接触刚度越软，可能造成的接触穿透量越大，如果刚度过小而没有相应地修改渗透容差，则会导致计算无法收敛。还没熟悉ansys底层接触刚度矩阵求解过程前，最好不要胡乱定义。

7）稳定阻尼因数

        APDL的关键词为FDMN，我们定义的接触可能由于接触和目标单元之间的节点或积分点之间的初始几何间隙导致接触对处于“开放”的状态。在分析过程中，“开放‘”状态的接触将不会被检测并计算，因此可能导致接触定义失败而产生刚体运动。稳定阻尼因子通过提供一定的阻力来阻止接触表面之间的相对运动，从而防止刚体运动。这个因子仅应用于接触法线方向，并且仅对接触类型为“无摩擦（Frictionless）”、“粗糙（Rough）”和“有摩擦（Frictional）”有效。阻尼应用于每个初始接触状态为开放的加载步。稳定化阻尼因子的值应足够大以防止刚体运动，但又足够小以确保解的存在。通常，值为1是合适的。