1. 接触非线性分析概述

在机械装配、密封连接、齿轮啮合、螺栓预紧、冲压成形、轮胎与路面作用等大量工程问题中，载荷并不是通过单一连续体传递的，而是通过两个或多个物体之间的表面相互作用传递。这类问题的共同特征是：接触区域的大小、位置和压力分布事先未知，且会随载荷、变形和材料状态不断变化。ANSYS 把这类问题归入接触非线性，它是结构非线性分析中最常见、也最容易出现收敛困难的一类。

一、接触问题的工程背景

以法兰螺栓连接为例：拧紧螺栓后，上下法兰的密封面被压紧，接触压力必须超过内部介质压力才能保证密封。此时接触面既不能互相穿透，又要能传递法向压力和切向摩擦力。若把两个法兰简化为一个整体网格，就无法描述“可能分开、可能滑移”的真实界面行为。接触分析的核心任务，就是让求解器自动判断：哪些节点对发生了接触、接触力有多大、是否发生相对滑移或分离。

接触问题的难点来自三个方面。第一，边界条件非线性：接触状态（分离、粘着、滑移）在求解前未知，只能在迭代中确定。第二，几何非线性：接触区往往伴随大位移、大转动，刚度矩阵需不断更新。第三，材料非线性：接触压力可能超过屈服极限，塑性变形反过来改变接触面积。三者耦合，使接触分析成为对建模经验和求解策略要求较高的任务。

二、非线性分类与接触的定位

结构非线性通常分为几何非线性、材料非线性和状态非线性（边界非线性）三类。接触属于状态非线性，其刚度矩阵随接触状态切换而突变，因此不能像线性静力分析那样一次求解，必须采用增量迭代。ANSYS 中常用的非线性求解控制包括载荷步、子步、平衡迭代和收敛准则，接触分析对这些设置尤其敏感。

三、ANSYS 接触分析能力

ANSYS 提供多种接触建模方式。面—面接触使用 TARGE170 定义目标面、CONTA174 定义接触面，适用于大多数三维装配问题；点—面接触适合螺栓头与法兰、铆钉等局部区域；点—点接触多用于梁、杆系或已知接触点对的问题。接触算法上，罚函数法通过引入法向刚度近似约束穿透，计算效率高但会产生少量穿透；拉格朗日乘子法严格满足不可穿透条件，但可能引入额外自由度并增加收敛难度；增广拉格朗日法在罚函数基础上迭代修正，兼顾精度与稳定性，是工程中较常用的选择。

摩擦方面，ANSYS 支持库仑摩擦模型，通过摩擦系数 μ 描述切向力与法向力的关系，并可设置等效剪应力极限，避免在高压区出现不合理的过大摩擦力。对于绑定或不分离的连接，可使用绑定接触或 MPC 算法，将接触面自由度耦合，从而把非线性接触退化为线性约束，显著降低求解成本。

四、典型应用场景与建模要点

接触分析在工程中覆盖面很广。密封与法兰连接关注接触压力是否连续、是否超过密封所需阈值；齿轮与轴承关注接触斑位置、接触应力和次表面剪应力；冲压成形关注板料与模具的贴合、回弹和起皱；轮胎接地关注接地印痕形状与摩擦分布；螺栓预紧关注夹紧力在接触面上的分配和螺纹附近的应力集中。

建模时需注意几个关键点。第一，接触面网格应适当加密，目标面可略粗，但两者网格密度不宜相差过大，否则容易引起接触压力振荡。第二，法向刚度因子 FKN 不宜盲目取大，一般先取 1.0 左右，再根据穿透量和收敛情况调整。第三，初始接触状态要合理，可通过调整几何间隙或使用“adjust to touch”让接触对在求解开始时处于刚好接触状态。第四，求解应打开自动时间步和线性搜索，必要时启用大变形效应。

五、结果判读与常见问题

接触分析的结果不能只看等效应力云图。应重点检查：接触状态图（NearContact、Sliding、Sticking、FarOpen）是否与预期一致；接触压力是否出现尖峰或负值；穿透量是否在可接受范围内（一般应小于单元尺寸的 5%～10%）；力收敛曲线是否在允许容差内平稳下降。若出现穿透过大，可提高 FKN 或加密接触面网格；若出现收敛困难，可减小初始时间步、降低 FKN、改用增广拉格朗日法，或检查是否存在刚体位移和初始穿透。