2. 接触力学基础

接触问题之所以在 ANSYS 中比普通线性静力分析更难处理，根本原因在于：两个物体之间的“接触”并不是一个始终存在的边界条件，而是一种随载荷、位移和间隙变化而不断切换的状态。理解接触界面的力学行为，是后续设置接触对、选择算法、判断收敛问题的前提。

一、接触界面的力学行为

当两个表面靠近时，界面上的力学行为可以分为三个层次：法向行为、切向行为和状态切换。法向行为决定两个面是否互相穿透、能否传递压力；切向行为决定接触面之间是否发生相对滑移、能否传递摩擦力；状态切换则决定接触点当前处于“分离”“粘着”还是“滑移”状态。

在 ANSYS 中，接触对通常由接触面（contact surface）和目标面（target surface）组成。接触面一般选择较软、网格较密的一侧，目标面选择较硬、网格较粗的一侧。接触单元（如 CONTA174）会监测接触面节点与目标面之间的相对位置，并根据法向间隙和切向滑移量判断当前接触状态。

接触界面的一个重要特征是：接触压力只能为压，不能为拉。当法向间隙大于零时，接触面分离，法向压力为零；当法向间隙等于零或试图穿透时，接触面闭合，法向压力大于零。这种“单边约束”是接触非线性的核心来源之一。

二、法向接触条件

法向接触条件可以用三个量描述：法向间隙 g_n、法向接触压力 p_n 和法向接触刚度 k_n。它们之间的关系为：

若 g_n > 0：p_n = 0（分离状态）
若 g_n ≤ 0：p_n = k_n · (-g_n)（闭合状态，压力与穿透量成正比）

这里 g_n 的正负号约定为：间隙为正表示分离，间隙为负表示穿透。在 ANSYS 中，程序通过罚函数法或拉格朗日乘子法来施加这个条件。罚函数法允许微小穿透，用接触刚度 k_n 将穿透量转化为接触压力；拉格朗日乘子法则引入额外自由度，强制穿透为零，但会增加方程规模和收敛难度。

实际设置中，法向接触刚度 k_n 通常用法向刚度因子 FKN 控制。FKN 的默认值一般为 1.0，表示程序根据材料弹性模量和接触面尺寸自动估算刚度。FKN 越大，穿透越小，但收敛越困难；FKN 越小，收敛容易，但穿透可能过大，影响结果精度。对于以弯曲为主的接触问题，FKN 可适当取小（如 0.1～0.5）；对于以压缩为主、需要精确接触压力的场合，FKN 可适当取大（如 1.0～10）。

三、切向接触条件与摩擦模型

切向接触条件描述接触面之间的摩擦行为。ANSYS 中常用的摩擦模型是库仑摩擦模型：当切向力 F_t 小于极限摩擦力 μ·p_n 时，接触面处于粘着状态，无相对滑移；当切向力达到极限摩擦力时，接触面开始滑移，切向力保持为 μ·p_n。

其中 μ 为摩擦系数，可通过材料参数 MP, MU 或接触实常数设置。摩擦系数越大，切向约束越强，但收敛也越困难。对于无摩擦接触，可将 μ 设为 0，此时接触面只能传递法向压力，不能传递切向力。

ANSYS 还提供了粘着刚度（FKT）来控制粘着状态的切向刚度。FKT 默认值为 1.0，表示程序自动估算。FKT 过小会导致粘着状态判断不准，FKT 过大则可能引起切向振荡。对于摩擦主导的接触问题，建议先使用默认值，若收敛困难再逐步调整。

四、接触非线性来源

接触非线性主要来自三个方面：

• 状态非线性：接触点可以在分离、粘着、滑移之间切换，每次切换都会改变系统的刚度矩阵。
• 几何非线性：大变形下接触面的法向和切向方向会随变形改变，接触搜索和压力方向需要不断更新。
• 材料非线性：接触区附近常伴随塑性、超弹性等材料非线性，与接触非线性耦合后收敛难度显著增加。

在求解过程中，ANSYS 通过 Newton-Raphson 迭代逐步施加荷载。每次迭代中，程序会检查接触状态，更新接触刚度，计算不平衡力，直到残差满足收敛准则。若接触状态频繁切换，或接触刚度突变，迭代可能发散。此时可通过调整 FKN、FKT、摩擦系数、时间步长或使用增广拉格朗日法来改善收敛。

五、接触刚度与摩擦参数的设置建议

在实际分析中，接触刚度和摩擦参数需要根据问题类型调整。以下是一份可直接参考的初始设置表：

理解接触界面的法向与切向条件、接触刚度的物理意义以及摩擦模型的数学表达，是后续在 ANSYS 中正确设置接触对、选择求解器选项和诊断收敛问题的理论基础。接触非线性不是“软件设置问题”，而是力学行为在数值计算中的真实反映。