📚 课程目录第4章接触力学基本概念

4. 接触力学基本概念

在光机结构非线性接触仿真中，接触界面是决定整体刚度、热变形传递和微振动响应的关键区域。要建立可复现、可收敛的接触模型，必须先建立一套清晰的物理量体系。本节围绕六个核心概念展开：法向接触压力、切向摩擦剪应力、接触刚度、穿透容差、接触状态与接触斑。它们共同构成后续接触定义、求解设置与结果判读的基础。

一、法向接触压力：界面载荷的传递量

法向接触压力（Contact Pressure）是接触面之间沿公法线方向传递的单位面积力，量纲为压强（Pa 或 MPa）。在有限元中，它通常以节点或单元积分点的法向应力形式输出，记为 PRESSURE 或 CONTACT_PRESSURE。对于光学镜座与镜筒的螺纹压圈接触、镜面与柔性支撑的球头接触，法向压力直接决定界面是否发生局部屈服、是否产生压印或应力双折射。

判断法向压力是否合理，需要结合材料许用值。例如某铝合金镜座与玻璃镜体接触，铝合金屈服强度约 270 MPa，玻璃抗压强度约 800 MPa，但玻璃表面微裂纹对拉应力敏感。若仿真得到峰值接触压力 120 MPa，虽低于两者强度，仍需检查压力分布是否集中在边缘。工程上常将峰值压力控制在材料屈服强度的 30%～50% 以内，并尽量使压力分布均匀。

二、切向摩擦剪应力：抵抗相对滑移的界面力

切向摩擦剪应力（Friction Shear Stress）是接触面内抵抗相对滑动的剪应力，记为 TAU 或 SHEAR_STRESS。它由法向压力和摩擦模型共同决定。常用库仑摩擦模型为：

τ_max = μ · p_n
|τ| < τ_max  → 黏着（sticking）
|τ| = τ_max  → 滑移（sliding）

其中 μ 为摩擦系数，p_n 为法向接触压力。光机结构中，镜筒与镜座若采用无润滑的铝-铝接触，μ 可取 0.2～0.4；若表面有润滑脂或二硫化钼涂层，μ 可降至 0.05～0.15。摩擦剪应力过大可能导致镜片在温度循环中无法自由释放热变形，从而引入面形误差；过小则可能在振动下产生微滑移，引起指向漂移。

三、接触刚度：数值稳定与物理真实的平衡

接触刚度（Contact Stiffness）是罚函数法中用于约束穿透的惩罚参数，分为法向刚度 FKN 和切向刚度 FKT。在 ANSYS 等软件中，法向接触刚度常以系数形式给出，实际刚度由下式估算：

K_n = FKN · E · √A
K_t = FKT · K_n

其中 E 为接触对中较软材料的弹性模量，A 为接触面积。FKN 越大，穿透越小，但收敛越困难；FKN 越小，收敛容易，但穿透偏大，压力分布失真。经验上，以刚度为主的接触问题 FKN 取 0.1～1.0，以柔度为主或大变形问题取 0.01～0.1。切向刚度 FKT 通常取 FKN 的 0.1～1.0 倍。若出现“接触刚度太大导致不收敛”或“穿透过大导致结果不可信”，应优先调整 FKN 并观察穿透量与压力峰值的变化趋势。

四、穿透容差：判断接触是否成立的数值尺度

穿透容差（Penetration Tolerance）是允许的节点侵入量，记为 TOLN，默认值通常为 0.1×底层单元厚度。它决定了求解器何时认为“接触已建立”。若实际穿透小于容差，接触压力可能为零或很小；若大于容差，则产生惩罚力。对于光学镜面，穿透容差应远小于面形精度要求。例如面形 PV 要求 λ/10（约 63 nm），则穿透容差应控制在 10 nm 量级以下，否则数值穿透会直接污染面形结果。

穿透容差不是物理间隙，而是数值判据。设置过大会掩盖真实接触状态，设置过小会引发收敛震荡。建议先以默认值试算，再根据穿透云图与压力云图逐步收紧。

五、接触状态：从张开到滑移的连续谱

接触状态描述接触对在某一载荷步下的相对运动与受力状态，通常分为四类：

张开（Open / Far）：间隙大于零，法向压力为零，无摩擦剪应力。
黏着（Sticking）：间隙为零，法向压力大于零，切向剪应力小于极限摩擦应力，无相对滑移。
滑移（Sliding）：法向压力大于零，切向剪应力达到极限值，存在相对滑动。
穿透（Penetrating）：节点侵入目标面，通常由刚度过小或容差过大引起，需检查。

在光机装配中，温度升高时镜座膨胀可能使镜片与压圈由黏着转为滑移；温度降低时可能由黏着转为张开。接触状态的变化会直接改变结构刚度，因此非线性接触仿真必须输出每个载荷步的接触状态云图，而不能只看最终结果。

六、接触斑：真实承载区域的几何表达

接触斑（Contact Patch）是实际发生接触的区域，其面积、形状和位置决定了载荷传递路径。理想球面-平面接触在赫兹理论下为圆形接触斑，但光机结构中的接触斑常因加工误差、装配偏心和热变形而呈椭圆、月牙或离散多点。接触斑面积 A_c 与平均压力 p_avg 满足：

F_n = p_avg · A_c

若仿真得到的接触斑远小于设计预期，说明接触面存在局部高点或间隙，需要检查网格密度、几何公差和接触探测方式。对于镜面支撑，接触斑应尽量对称、连续，避免出现孤立小斑导致局部应力集中。

七、参数设置与结果判读的实操要点

将上述概念落到仿真操作中，可按以下顺序建立接触模型。首先，根据接触对材料确定摩擦系数 μ，并在接触定义中设置法向刚度系数 FKN 和切向刚度系数 FKT。其次，设置穿透容差 TOLN，建议初始取底层单元厚度的 0.1 倍。然后，在求解器中开启大变形和接触状态输出。求解完成后，依次检查：穿透云图是否超过 TOLN；接触压力峰值是否低于材料许用值；接触状态云图中滑移区是否位于预期位置；接触斑面积是否与设计承载区域一致。

以一个典型示例说明：某镜座与镜体接触，镜座材料为铝合金（E=70 GPa），镜体为微晶玻璃（E=90 GPa），接触面积约 200 mm²，摩擦系数 μ=0.2，法向载荷 500 N。若取 FKN=0.5，则法向接触刚度估算为 K_n = 0.5 × 70 GPa × √200 mm² ≈ 4.95×10⁵ N/mm。平均接触压力 p_avg = 500 N / 200 mm² = 2.5 MPa。极限摩擦剪应力 τ_max = 0.2 × 2.5 = 0.5 MPa。若仿真得到峰值压力 8 MPa，说明接触斑实际面积小于设计值，需检查几何公差或网格密度。若穿透量达到 0.05 mm，而 TOLN 设为 0.01 mm，则说明 FKN 偏小，应提高至 1.0 或更大后重新试算。

接触刚度、穿透容差和摩擦系数三者相互耦合。不要一次性大幅调整多个参数，否则难以判断收敛失败或结果异常的原因。建议每次只改一个参数，并记录穿透量、压力峰值和接触斑面积的变化。

掌握这六个物理量之后，后续的接触建模、求解设置和结果判读就有了统一的参照系。法向压力回答“力有多大”，摩擦剪应力回答“会不会滑”，接触刚度回答“数值上如何约束”，穿透容差回答“多小算接触”，接触状态回答“当前处于什么阶段”，接触斑回答“实际承载在哪里”。六者共同构成光机结构非线性接触仿真的物理量基础。