3. 装配界面与结合面物理力学
精密机构中，真正决定装配后运动品质的往往不是零件本体的尺寸精度，而是两个零件“贴合在一起”的那一层界面。无论导轨滑块与导轨、轴承内外圈与轴/座、法兰端面、还是螺纹副与定位销，载荷与运动都要穿过这些结合面传递。结合面在微观上并不是理想平面，而是由大量微凸体峰顶接触形成的离散接触点集合，因此它同时表现出接触刚度有限、摩擦耗能、微滑移、迟滞与非线性等特征。理解并建模这些特征，是预测机构定位精度、重复性、动态响应和热稳定性的基础。
一、表面形貌：结合面物理行为的起点
任何加工表面都存在多尺度起伏。工程上通常把表面形貌分解为三个尺度：形状误差（宏观，如平面度、圆柱度，波长毫米级）、波纹度（中观，波长0.1～1 mm）、粗糙度（微观，波长微米级）。粗糙度参数中，Ra（算术平均偏差）最常用，但它只反映高度统计，不反映峰顶曲率与分布；对接触刚度更敏感的是峰顶密度、峰顶曲率半径和高度分布。因此仅凭图纸标注 Ra 0.4 μm 并不能直接推算接触刚度，需要结合加工工艺（磨削、研磨、刮研、珩磨）判断形貌类型。
以典型精密磨削钢平面为例：Ra 约 0.2～0.4 μm，峰顶密度约 10⁴～10⁵ 个/mm²，微凸体峰顶曲率半径约 1～10 μm，高度分布接近高斯。这些量级决定了：在名义压应力 0.1～1 MPa 下，真实接触面积只占名义面积的 0.1%～5%，单个微凸体接触直径常在亚微米到几微米量级。
二、接触刚度：从微凸体到结合面
单个微凸体在弹性阶段可近似为球-平面接触，其法向力与变形满足 Hertz 关系 F ∝ δ^(3/2)。大量微凸体叠加后，结合面法向刚度 Kn 随法向载荷呈幂律增长，典型形式为 Kn = a·Fm，m 约 0.5～0.8。这意味着预紧力越大，结合面越“硬”，但增益递减。切向刚度 Kt 通常小于 Kn，比值 Kt/Kn 约 0.2～0.7，取决于表面状态与润滑。
工程估算中常用单位面积刚度：磨削钢-钢干结合面，法向刚度约 100～500 MPa/μm（即 10¹¹～5×10¹¹ N/m³），切向约 30～200 MPa/μm。若结合面名义面积 1000 mm²，则 Kn 约 1×10⁸～5×10⁸ N/m。这个量级与一根 100 mm 长、直径 20 mm 钢轴的轴向刚度（约 3×10⁸ N/m）相当，说明结合面刚度不可忽略，常是整机刚度链的薄弱环节。
三、摩擦、微滑移与迟滞
结合面在切向载荷下并非“要么不动、要么滑动”。当切向力小于极限摩擦力 μFn 时，接触边缘先发生局部微滑移，中心仍粘着，形成部分滑移。卸载时切向位移不沿原路返回，形成迟滞回线，其面积即一个循环的耗散能量。这一机制带来两个后果：一是微动磨损与颗粒生成，污染精密导轨与轴承；二是定位重复性下降，因为每次加载路径不同，残余切向位移不同。
精密机构中，微滑移幅值常在 0.1～5 μm。若定位重复性要求 1 μm，则必须控制切向载荷不超过极限摩擦力的 30%～50%，或采用高预紧、高摩擦系数（如刮研面、涂层）来扩大粘着区。润滑会降低 μ，扩大微滑移区，因此精密定位结合面常采用干接触或微量润滑。
四、结合面建模方法的选择
工程上不存在“万能模型”，选择取决于精度要求、可用参数和计算资源。常见四类方法及其适用边界如下：
| 方法 | 核心输入 | 输出 | 适用场景 | 主要局限 |
|---|---|---|---|---|
| GW 统计模型 | Ra、峰顶密度、曲率半径、弹性模量 | 真实接触面积、Kₙ(F) | 弹性干接触、快速估算 | 忽略塑性、粘着与尺度效应 |
| 分形模型 | 分形维数 D、尺度系数 G | 跨尺度 Kₙ、真实面积 | 磨削/研磨面、多尺度分析 | D、G 标定困难，物理意义争议 |
| 有限元接触 | 真实形貌或等效粗糙面、弹塑性参数、μ | 局部应力、位移、滑移区 | 关键结合面强度/刚度校核 | 计算量大，形貌测量与网格成本高 |
| 等效弹簧-阻尼 | Kₙ、Kₜ、cₙ、cₜ（实验或经验） | 整机模态、频响、定位误差 | 系统级动力学与精度仿真 | 参数不准则结果失真 |
五、实操：结合面参数辨识与验证
以一台精密直线导轨副为例，说明如何获得可用于整机模型的结合面参数。目标：辨识滑块-导轨结合面的法向刚度 Kn、切向刚度 Kt 与摩擦系数 μ。
- 法向刚度辨识：将滑块固定在刚性基座上，沿法向施加阶梯载荷（如 0→500 N→1000 N→1500 N→0），用分辨率 0.1 μm 的电感测微仪测量滑块相对导轨的法向位移。绘制 F-δ 曲线，对每一级载荷取割线刚度。示例数据（合理假设）：500 N 时 δ=2.1 μm，K=238 N/μm；1000 N 时 δ=3.6 μm，K=278 N/μm；1500 N 时 δ=4.9 μm，K=306 N/μm。可见刚度随载荷增大而提高，符合幂律趋势。
- 切向刚度与微滑移：保持法向预紧 1000 N，施加切向循环载荷 0→200 N→0→-200 N→0，记录切向位移。若回线面积小、残余位移小于 0.2 μm，说明处于粘着主导区；若回线明显张开，则进入部分滑移。示例：切向力 200 N 时位移 1.8 μm，卸载后残余 0.3 μm，回线面积约 0.12 mJ，提示微滑移已不可忽略。
- 摩擦系数估计：逐步增大切向力直至出现连续滑移，记录临界力 Fs。若 Fn=1000 N 时 Fs=180 N，则 μ≈0.18。该值偏低，说明导轨面可能残留润滑膜，需清洁或调整预紧。
- 模型验证：将辨识得到的 Kn=280 N/μm、Kt=90 N/μm、μ=0.18 代入整机有限元模型，计算一阶固有频率并与锤击法实测对比。若误差超过 10%，应检查结合面面积折算、螺栓预紧分布或阻尼参数。
综上，结合面建模的核心逻辑是：形貌决定接触，接触决定刚度与摩擦，摩擦与微滑移决定迟滞和重复性。工程实践中，应根据精度目标选择建模层级——快速估算用 GW 或经验公式，关键结合面用有限元校核，整机动力学用等效弹簧-阻尼，并通过实验辨识与验证闭环，才能让模型真正服务于装配质量与运动品质的预测。