📚 课程目录第2章接触力学基础与Ansys实现映射

2. 接触力学基础与Ansys实现映射

在 Mechanical 中设置接触，本质上是在把接触力学的几个基本量翻译成软件可识别的参数。如果只记住“用 Bonded 还是 Frictional”这类操作，遇到收敛困难或结果异常时往往无从下手。本节把法向接触、切向摩擦、接触刚度和穿透容差这四个核心概念，与 Mechanical 中对应的设置项逐一对应起来，说明每个参数背后的物理含义，以及它如何影响求解行为。

法向接触：压力只能“推”，不能“拉”

两个物体接触时，法向方向上的相互作用有一个基本约束：接触面之间只能传递压力，不能传递拉力。也就是说，法向接触力 F_n 必须满足 F_n ≥ 0。当两个面分离时，接触力为零；当它们相互压紧时，接触力为正。这个“单边约束”是接触问题非线性的根源——接触状态会随载荷变化在“分离”和“接触”之间切换。

在 Mechanical 中，法向行为由 Contact > Definition > Behavior 控制。选择 Asymmetric 或 Symmetric 决定接触对是否双向检测；选择 Program Controlled 时，软件会根据接触类型自动决定法向算法。真正决定“能不能穿透、允许多少穿透”的是 Advanced > Normal Stiffness 和 Normal Lagrange 等选项。纯拉格朗日法（Normal Lagrange）通过引入额外自由度强制满足不可穿透条件，穿透量理论上为零，但会带来更多方程和收敛难度；罚函数法（Pure Penalty）允许微小穿透，用罚刚度把穿透“压回去”，收敛更平稳，是默认选择。

切向摩擦：从粘着到滑移的过渡

切向行为描述接触面之间的摩擦。最常用的是库仑摩擦模型：当切向力 F_t 小于摩擦系数 μ 与法向力 F_n 的乘积时，接触点处于粘着（Sticking）状态，切向位移被约束；一旦 F_t 达到 μF_n，接触点开始滑移（Sliding），切向力维持在极限值附近。

在 Mechanical 中，切向行为由 Contact > Definition > Friction Coefficient 设置。摩擦系数为 0 表示无摩擦，接触面可自由滑动；输入正值则启用库仑摩擦。对于需要考虑粘着刚度的场景，可在 Advanced > Tangential Stiffness 中选择 Program Controlled 或手动指定。摩擦系数不是“越大越好”，它直接影响切向刚度矩阵和收敛性：摩擦系数过大且法向刚度很高时，粘着区难以释放，容易出现切向振荡。

接触刚度：罚函数法的核心参数

接触刚度 k_n 是罚函数法中把穿透量换算成接触力的比例系数：F_n = k_n · x_penetration。刚度越大，同样的穿透量产生的接触力越大，接触面越“硬”，结果越接近真实不可穿透；但刚度过大会导致刚度矩阵病态，收敛困难。刚度过小则穿透明显，接触力被低估，结果偏软。

Mechanical 默认使用 Program Controlled，由软件根据材料弹性模量和单元尺寸估算一个合适的罚刚度。用户可在 Advanced > Normal Stiffness > Manual 中指定 Normal Stiffness Factor，该因子乘以程序估算值得到实际刚度。工程经验是：先让程序控制跑通，若发现穿透量明显（例如超过单元尺寸的 5%～10%），再逐步增大因子（如从 1.0 调到 2.0、5.0），每次观察收敛曲线和穿透量，而不是一次性调到很大。

穿透容差：判断接触是否“合格”的尺子

穿透容差（Penetration Tolerance）是判断接触状态和收敛的参考量。在 Advanced > Penetration Tolerance 中，可以设置 Factor 或 Value。它决定了求解器认为“多大的穿透是可以接受的”。如果实际穿透超过容差，求解器会继续迭代调整接触力；如果穿透在容差内，则认为该接触点已满足约束。

需要区分两个概念：允许穿透是罚函数法固有的数值现象，穿透容差是判断收敛的阈值。前者由刚度决定，后者由用户或程序设定。查看结果时，可通过 Solution > Contact Tool > Penetration 和 Gap 云图检查实际穿透分布，若最大穿透远大于容差，说明刚度偏软或网格偏粗。

参数之间的耦合关系

这四个概念不是孤立的。法向刚度决定穿透量，穿透量又影响接触面积和法向力分布；法向力决定摩擦极限 μF_n，从而影响粘着/滑移区的比例；穿透容差则决定求解器在什么精度上停止迭代。实际调试时，常见的问题是“接触面穿透过大导致结果偏软”，此时应先检查网格是否足够细，再考虑增大法向刚度因子，而不是直接改摩擦系数。

实用判断依据：在 Contact Tool 中查看 Penetration 云图，若最大穿透量超过接触面附近单元特征尺寸的 10%，通常说明罚刚度偏软或网格偏粗。此时可先将 Normal Stiffness Factor 从 1.0 提高到 2.0，重新求解并对比穿透云图和等效应力，确认结果趋于稳定后再决定是否继续增大。

一个可复现的示例设置

假设一个钢制滑块压在钢制导轨上，滑块顶面施加 500 N 向下压力，摩擦系数取 0.15。在 Mechanical 中可这样设置：接触对类型选 Frictional，Friction Coefficient 填 0.15；Behavior 保持 Program Controlled；Advanced 中 Normal Stiffness 先保持 Program Controlled，Penetration Tolerance 保持默认。求解后查看 Contact Tool 的 Penetration 和 Sliding Distance 云图。若穿透最大值约为 0.002 mm，而接触区单元尺寸约 0.5 mm，则穿透占比约 0.4%，属于可接受范围；若穿透达到 0.05 mm 以上，则应把 Normal Stiffness Factor 调到 2.0 或 5.0 再试。