4. 摩擦接触与非线性求解控制

摩擦接触是装配体仿真中最常见的非线性来源之一。与绑定接触不同，摩擦接触的切向力不是由几何约束直接决定，而是取决于法向压力、摩擦系数以及接触点当前所处的粘着或滑移状态。Ansys Mechanical 在求解时需要在每个迭代步中反复判断接触状态，因此摩擦接触天然带有路径依赖和收敛困难。理解摩擦系数、粘滑行为、大变形开关、载荷步与子步以及收敛诊断之间的关系，是保证装配可靠性仿真结果可信的关键。

摩擦系数设置：从物理输入到求解器行为

在 Mechanical 中，摩擦接触通过 Contact 区域的 Definition 设置。将 Type 设为 Frictional 后，需要输入 Friction Coefficient。该系数是库仑摩擦模型中的 μ，切向滑移开始的条件为 τ = μ·p，其中 p 为接触法向压力，τ 为切向剪应力。当切向应力小于 μ·p 时，接触点处于粘着状态；达到或超过该阈值后进入滑移状态。

摩擦系数不是单一常数。工程中常区分静摩擦系数与动摩擦系数，Ansys 允许通过 Tabular 数据或 Command 对象定义随压力、温度或滑移速度变化的摩擦系数。对于钢-钢干摩擦，示例可取 μ = 0.15～0.20；对于有润滑的螺栓连接面，示例可取 μ = 0.08～0.12；对于橡胶与钢，示例可取 μ = 0.6～0.8。这些数值属于教学示例，实际项目应依据试验或标准手册确定。

还需注意两个与摩擦系数配合的选项：一是 Stiffness Update，摩擦接触的切向刚度会随粘滑状态变化，建议设为 Each Iteration 或 Aggressive，使求解器及时更新；二是 Normal Stiffness，法向刚度因子过大会导致接触压力振荡，过小则产生穿透，通常从 1.0 开始，收敛困难时降至 0.1～0.5 试探。

粘滑行为：路径依赖与能量耗散

粘滑行为是摩擦接触最典型的非线性特征。粘着阶段接触点无相对滑移，切向刚度承担载荷；滑移阶段切向力被限制在 μ·p 以内，多余能量通过摩擦耗散。装配体在螺栓预紧、过盈配合或循环载荷下，往往经历“粘着—局部滑移—整体滑移”的演化过程。

粘滑判断对网格和载荷步非常敏感。若载荷步过大，求解器可能在一个子步内跨越粘着到滑移的突变，导致残差剧烈振荡甚至发散。因此摩擦问题通常需要比线性问题更细的子步划分，并开启自动时间步。粘滑还会造成接触压力重分布，使装配体刚度在加载过程中发生变化，这也是摩擦接触必须使用非线性求解器的原因。

大变形开关：几何非线性与摩擦的耦合

Large Deflection 开关控制是否考虑几何非线性。当装配体在载荷下产生较大位移或转动时，例如悬臂支架、柔性卡扣或大行程机构，必须开启大变形。开启后，求解器在每个子步更新几何构型，接触法向和切向方向随之变化，摩擦力的方向也随接触面转动而更新。

若关闭大变形而实际变形较大，接触压力分布会失真，摩擦滑移方向也可能错误，导致装配体刚度被高估。反之，若变形确实很小却开启大变形，会增加计算量但通常不影响结果。判断依据是最大位移与特征尺寸之比：当该比值超过约 5% 时，建议开启大变形。对于螺栓预紧这类局部接触问题，即使整体位移小，接触面局部变形也可能显著，因此常建议开启。

载荷步与子步：控制非线性路径

载荷步用于区分不同加载阶段，例如“预紧—外载—卸载”。每个载荷步内可设置初始子步数、最小子步数和最大子步数。摩擦接触问题建议初始子步不少于 10，最小子步 5～10，最大子步 100～200，并开启 Auto Time Stepping。自动时间步会根据收敛情况自动二分或放大步长：收敛顺利时放大，收敛困难时缩小。

子步过少会使求解器在一个步内跨越粘滑突变，残差无法下降；子步过多则计算时间增加。一个实用的起点是：预紧步用 10～20 个子步，外载步用 20～50 个子步，并允许自动时间步在 5～200 之间调整。若出现“接触状态变化数”在某子步突然增大，说明该处发生粘滑切换，应进一步细分该载荷区间。

收敛诊断与调试策略

Mechanical 的 Solution Information 提供力残差、位移残差、接触状态变化数、穿透量等曲线。判断收敛是否健康，不能只看“是否收敛”，还要看残差是否单调下降、接触状态变化数是否在后期趋于零。若残差在 1e-2 附近振荡，常见原因是法向刚度因子过大或子步过粗。

调试可按以下顺序进行：第一，将法向刚度因子从 1.0 降至 0.1～0.5；第二，增加子步并开启自动时间步；第三，将摩擦接触的 Stiffness Update 设为 Each Iteration；第四，若仍不收敛，将求解器切换为 Unsymmetric，以处理摩擦引起的非对称刚度矩阵；第五，检查 Pinball 区域是否过大或过小，过大导致远距离误接触，过小导致接触丢失。每一步调整后都应观察残差曲线和接触状态变化数，而不是盲目叠加设置。

一个可复用的求解控制配置示例

以下为某钢制支架螺栓连接装配体的教学示例配置，参数为合理假设，实际项目需按试验或规范调整：

接触设置：
  Type              = Frictional
  Friction Coeff    = 0.15（钢-钢干摩擦示例）
  Normal Stiffness  = 0.5（初始试探值）
  Stiffness Update  = Each Iteration
  Pinball Region    = 2 mm

分析设置：
  Large Deflection  = On
  Auto Time Step    = On
  Initial Substeps  = 20
  Min Substeps      = 10
  Max Substeps      = 200

载荷步：
  Step 1  螺栓预紧  10 子步
  Step 2  外载      20 子步
  Step 3  卸载      10 子步

收敛准则：
  力残差  1e-3
  位移残差 1e-3
  接触容差 默认

求解后应重点查看接触压力分布是否连续、滑移距离是否集中在预期区域、装配体整体刚度是否与手算或试验趋势一致。若接触压力在边缘出现尖峰且穿透量超过允许值，说明法向刚度仍偏大或网格过粗；若滑移距离远大于预期，应复核摩擦系数和载荷方向。通过“设置—求解—诊断—调整”的闭环，才能让摩擦接触结果既收敛又可信。