4、摩擦模型基础：静摩擦、库仑摩擦、粘性摩擦、Stribeck曲线的数学描述与物理意义
速度环的低速平稳性优化，本质上是在与摩擦这一非线性扰动做斗争。在伺服系统的低速区间（通常指转速低于额定转速的5%~10%），摩擦力的影响远大于惯性力和电磁推力，成为限制轮廓精度、引发爬行（stick-slip）现象和换向过零畸变的主要根源。要设计出有效的摩擦补偿算法，首先必须对摩擦的静态与动态特性建立精确的数学模型。本章将系统梳理静摩擦、库仑摩擦、粘性摩擦以及Stribeck曲线的数学表达与物理内涵，为后续的摩擦前馈补偿、扰动观测器设计奠定理论基础。
上图展示了经典摩擦模型的完整构成。在零速附近，静摩擦（F_s）起主导作用；当外力克服静摩擦后，系统进入边界润滑区，摩擦力随速度增加而下降，形成Stribeck负斜率段；随后进入部分润滑区，摩擦力达到最小值（即库仑摩擦F_c）；随着速度继续升高，粘性摩擦项B·v逐渐占据主导，摩擦力近似线性上升。这张图是理解低速爬行、换向过零误差以及摩擦前馈补偿设计的核心参照。
4.1 静摩擦与库仑摩擦：低速域的两大“顽固分子”
4.1.1 核心原理：静摩擦的“死区”特性与库仑摩擦的“恒定阻力”
静摩擦（Static Friction）是指物体在接触面之间尚未发生相对运动时产生的阻力。其数学描述为：当速度v=0且外力|F_ext| < F_s时，摩擦力F_f = F_ext（即摩擦力完全抵消外力，保持物体静止）；当|F_ext| ≥ F_s时，物体开始运动，摩擦力瞬间降至库仑摩擦水平。静摩擦的物理本质是接触面微凸体之间的机械互锁和分子粘着效应。在伺服系统中，静摩擦表现为位置环的“死区”——当指令位置增量对应的驱动力矩小于静摩擦力矩时，电机轴无法转动，导致轮廓加工中的圆角误差和跟随误差。
库仑摩擦（Coulomb Friction）是物体在运动过程中，与速度方向相反、大小近似恒定的阻力。其数学表达式为：F_f = F_c · sign(v)，其中F_c为库仑摩擦力常数，sign(v)为符号函数。库仑摩擦来源于接触面粗糙度引起的机械啮合和犁沟效应，其大小与法向载荷成正比，与速度大小基本无关。在速度环中，库仑摩擦是换向过零时产生“平顶”或“削顶”现象的直接原因——当电机反向运动时，摩擦力方向突变，导致速度环输出力矩出现阶跃性跳变。
| 摩擦类型 | 数学表达式 | 物理来源 | 对速度环的影响 |
|---|---|---|---|
| 静摩擦 F_s | v=0, |F_ext|<F_s → F_f=F_ext | 微凸体互锁、分子粘着 | 死区效应、启动延迟、爬行 |
| 库仑摩擦 F_c | F_f = F_c·sign(v) | 粗糙度啮合、犁沟效应 | 换向过零畸变、稳态误差 |
| 粘性摩擦 B·v | F_f = B·v | 润滑油膜剪切阻力 | 高速段阻尼、增益裕度下降 |
| Stribeck效应 | F(v)=F_c+(F_s-F_c)·e^{-(v/v_s)^2} | 润滑状态转变 | 负阻尼、极限环振荡 |
4.1.2 判断方法与操作流程：如何通过实验辨识F_s和F_c
工程实践中，通常采用“恒速扫描法”或“力矩斜坡法”来辨识静摩擦和库仑摩擦参数。推荐以下操作流程：
- 准备阶段：将伺服系统置于速度环模式，设定速度指令为0，记录电机编码器位置反馈，确保系统处于静止状态。
- 静摩擦辨识：以极慢的速度（如0.1 rpm/s）线性增加力矩指令，同时监测电机轴是否开始转动。记录轴开始转动瞬间的力矩指令值，即为静摩擦力矩F_s。重复5次取平均值，以消除随机误差。
- 库仑摩擦辨识：设定多个不同的恒定速度（如5、10、20、50 rpm），分别记录稳态时的力矩指令平均值。将不同速度下的力矩值绘制成曲线，外推至速度为零时的截距即为库仑摩擦力矩F_c。
- 数据记录：使用伺服驱动器的内部示波器功能，同步采集速度反馈和力矩指令，采样率建议不低于1 kHz，记录时间不少于2秒。
4.1.3 案例拆解：数控转台低速爬行问题的摩擦参数辨识
背景与约束：某五轴加工中心的数控转台（A轴）在低速（<2 rpm）旋转时出现明显的“一顿一顿”爬行现象，导致工件表面产生振纹。转台采用力矩电机直驱，配备高分辨率圆光栅（分辨率0.036°），减速比为1:1。系统额定力矩为120 Nm，最大转速为30 rpm。
初始现象与数据：通过驱动器采集速度环数据，发现当速度指令为1 rpm时，实际速度在0.2~1.8 rpm之间周期性波动，波动频率约为3 Hz。力矩指令在±15 Nm之间振荡，明显高于正常切削力矩（约5 Nm）。
分析判断过程：根据Stribeck理论，爬行现象通常发生在边界润滑区，此时摩擦力随速度增加而下降，形成负阻尼。负阻尼会与系统弹性环节（如联轴器、轴承刚度）耦合，产生极限环振荡。初步判断静摩擦F_s过大且Stribeck速度v_s较低。
实际操作步骤：① 按照上述流程辨识参数，测得F_s=18.5 Nm，F_c=6.2 Nm，v_s=0.8 rpm；② 对比正常转台（F_s≈8 Nm），发现静摩擦异常偏高；③ 检查转台轴承预紧力，发现预紧力矩过大，导致接触面压力增加；④ 适当降低轴承预紧力至厂家推荐值的80%，并更换低粘度润滑脂。
结果对比：调整后重新辨识，F_s降至7.8 Nm，F_c降至5.5 Nm。在1 rpm速度指令下，实际速度波动减小至±0.05 rpm以内，力矩指令稳定在±6 Nm。爬行现象完全消除，工件表面粗糙度从Ra 1.6 μm改善至Ra 0.4 μm。
风险与代价：降低轴承预紧力可能影响转台刚性，在重切削工况下需重新验证动态刚度。该案例表明，摩擦参数辨识不仅是算法问题，更是机械系统健康状态的“体检表”。
4.1.4 常见误区与总结复盘
- 误区一：认为F_s和F_c是固定常数。实际上，它们随温度、载荷、磨损状态变化，需要定期重新辨识。
- 误区二：忽略静摩擦的“预位移”效应。在微观尺度下，静摩擦并非刚性约束，而是存在微米级的弹性变形，这会影响纳米级定位精度。
- 误区三：将库仑摩擦补偿直接加在速度环输出端而不考虑方向切换的平滑处理，容易引起换向冲击。
4.2 粘性摩擦与Stribeck曲线：从线性到非线性的跨越
4.2.1 核心原理：粘性摩擦的线性模型与Stribeck的指数衰减描述
粘性摩擦（Viscous Friction）源于润滑油膜的内摩擦，其大小与相对运动速度成正比，数学表达式为F_f = B·v，其中B为粘性摩擦系数（单位：N·m·s/rad）。在高速区间，粘性摩擦是主要的阻尼来源，有助于抑制机械谐振。然而，过大的粘性摩擦会降低系统增益裕度，限制速度环带宽。
Stribeck曲线完整描述了摩擦力从静摩擦到库仑摩擦再到粘性摩擦的过渡过程。经典Stribeck模型（Tustin模型）的数学表达式为：
F(v) = F_c + (F_s - F_c) · e^{-(v/v_s)^2} + B·v (v ≠ 0)
F(v) = F_s (v = 0)其中，v_s为Stribeck特征速度，决定了负斜率段的宽度。当v < v_s时，摩擦力随速度增加而下降，形成负阻尼区；当v > v_s时，摩擦力随速度增加而上升，进入正阻尼区。v_s的典型值在0.1~10 rpm之间，取决于润滑方式和接触面材料。
Stribeck效应的物理本质是润滑状态的转变：在零速附近，接触面处于边界润滑状态，微凸体直接接触，摩擦力大；随着速度增加，润滑油被卷入接触面，形成部分弹流润滑，摩擦力下降；当速度足够高时，形成完全弹流润滑，摩擦力由粘性剪切主导，随速度线性上升。
4.2.2 判断方法与操作流程：Stribeck曲线的实验拟合
获取准确的Stribeck曲线是设计摩擦补偿器的前提。推荐以下实验流程：
- 速度扫描：在速度环模式下，从0.01 rpm到额定速度，按对数间隔选取15~20个速度点（如0.01、0.02、0.05、0.1、0.2、0.5、1、2、5、10、20、50、100、200、500、1000 rpm）。
- 稳态力矩采集：在每个速度点保持恒定速度运行2秒，取后1秒的力矩指令平均值作为该速度下的摩擦力矩。
- 正反向分别测量：由于摩擦的方向依赖性，需分别测量正向和反向运动时的摩擦力矩，取平均值作为最终数据。
- 曲线拟合：使用最小二乘法将实验数据拟合到Tustin模型，求解F_c、F_s、v_s和B四个参数。可使用MATLAB的lsqcurvefit函数或Python的scipy.optimize.curve_fit。
4.2.3 案例拆解：高速主轴Stribeck参数辨识与摩擦前馈补偿
背景与约束：某精密磨床的电主轴（额定转速8000 rpm）在低速磨削（50~200 rpm）时出现圆度误差超差。磨削力较小（约2 Nm），但主轴摩擦力矩在低速时波动明显。主轴采用永磁同步电机，配备旋转变压器（分辨率12位），驱动器速度环带宽为50 Hz。
初始现象与数据：在100 rpm恒速磨削时，实测主轴力矩指令在3.5~5.5 Nm之间波动，波动频率约8 Hz。通过FFT分析，发现波动频率与主轴转频（1.67 Hz）无关，而与Stribeck负阻尼区对应的极限环频率相关。
分析判断过程：初步怀疑Stribeck效应导致负阻尼，与主轴-砂轮系统的扭转刚度耦合形成振荡。通过实验辨识，得到F_s=4.8 Nm，F_c=2.1 Nm，v_s=35 rpm，B=0.0008 Nm·s/rad。注意到v_s=35 rpm，正好落在磨削工作区间（50~200 rpm）的边缘，说明磨削点处于部分润滑区，负阻尼效应显著。
实际操作步骤：① 在速度环中增加基于Stribeck模型的前馈补偿项，补偿力矩F_comp = F_c + (F_s - F_c)·e^{-(v/v_s)^2} + B·v；② 对补偿项进行低通滤波（截止频率200 Hz），避免引入高频噪声；③ 在换向时增加力矩斜坡（斜坡时间5 ms），平滑摩擦力的方向切换；④ 重新整定速度环PI参数，将比例增益降低10%，积分时间常数增加15%，以增加系统阻尼。
结果对比：补偿后，100 rpm磨削时的力矩波动从±1 Nm降至±0.2 Nm，圆度误差从12 μm改善至3 μm。同时，主轴在低速段的加减速响应时间缩短了20%，因为摩擦前馈减少了速度环的调节负担。
风险与代价：摩擦前馈补偿依赖于精确的模型参数，若参数漂移（如温度变化导致B值变化），补偿效果会下降。因此，该方案需要配合在线参数辨识或自适应补偿策略。此外，前馈补偿无法消除摩擦的动态效应（如预滑动位移），对于纳米级定位精度要求，还需引入扰动观测器。
4.2.4 常见误区与总结复盘
- 误区一：将Stribeck曲线简化为线性模型，忽略负阻尼区。这会导致低速振荡问题无法从根源上解决。
- 误区二：v_s的辨识不准确。v_s决定了负阻尼区的宽度，若v_s辨识偏大，补偿器会在高速区过度补偿，引入额外扰动。
- 误区三：忽视温度对B值的影响。在长时间运行后，电机温度升高，润滑脂粘度下降，B值可能下降30%以上，导致补偿不足。
4.3 摩擦模型的工程应用与选择策略
4.3.1 不同工况下的模型选择
并非所有应用都需要完整的Stribeck模型。根据控制精度和速度范围，可以选择不同复杂度的摩擦模型：
| 应用场景 | 推荐模型 | 理由 |
|---|---|---|
| 高速主轴（>1000 rpm） | 库仑+粘性摩擦 | Stribeck效应仅在低速区显著，高速区可忽略 |
| 通用伺服（10~1000 rpm） | 库仑+粘性+简单Stribeck | 兼顾精度与计算复杂度 |
| 精密转台/光刻机（<10 rpm） | 完整Stribeck+预滑动模型 | 低速精度要求极高，需考虑微观摩擦行为 |
| 机器人关节（频繁换向） | Stribeck+方向切换平滑 | 换向冲击是主要问题，需平滑处理 |
4.3.2 摩擦参数的温度补偿策略
温度是影响摩擦参数的主要因素。实测数据显示，当电机温度从25℃升至80℃时，粘性摩擦系数B可下降40%，静摩擦F_s下降15%。工程上可采用以下策略：
- 查表法：在不同温度点（如25℃、40℃、60℃、80℃）分别辨识摩擦参数，建立温度-参数查找表，运行时根据温度传感器实时插值。
- 在线辨识法：利用递推最小二乘法（RLS）在线估计B值，但需注意激励条件和收敛速度。
- 模型参考自适应：将摩擦参数作为自适应参数，通过Lyapunov稳定性理论设计自适应律，实时调整补偿量。
4.3.3 摩擦补偿与速度环整定的协同
摩擦补偿不能替代速度环PI参数的整定，两者需要协同工作。基本原则是：摩擦前馈负责消除已知的确定性扰动，速度环PI负责处理模型误差和外部负载扰动。在整定顺序上，应先进行摩擦补偿，再整定PI参数。否则，PI参数会为了克服摩擦而设置过大的积分增益，导致系统在换向时产生超调或振荡。
具体操作时，可先启用摩擦前馈，观察速度环的阶跃响应和抗扰性能，然后逐步增加PI增益，直到出现轻微振荡再回调10%~15%。这样既能保证系统刚度，又能避免摩擦引起的极限环。
4.3.4 案例复盘：摩擦模型选择不当导致的补偿失效
背景：某包装机械的伺服摆臂机构，工作速度范围0.5~300 rpm，频繁正反转。工程师最初采用简单的库仑+粘性摩擦模型进行前馈补偿，但在低速（<5 rpm）时出现明显的“点头”现象，即摆臂在启动瞬间先向前冲一小段再停止。
分析：“点头”现象源于静摩擦补偿不足。库仑+粘性模型在零速附近提供的补偿力矩仅为F_c，而实际需要F_s（约为F_c的2倍）。当速度环输出力矩超过F_c但小于F_s时，摆臂无法启动，导致位置误差累积；当误差足够大时，速度环输出力矩超过F_s，摆臂突然启动并冲过目标位置。
解决：将摩擦模型升级为完整Stribeck模型，并增加静摩擦的“预补偿”逻辑——当速度指令接近零且位置误差小于阈值时，额外施加一个短暂的静摩擦补偿脉冲（持续50 ms），帮助系统克服静摩擦死区。
结果：升级后，“点头”现象完全消失，摆臂定位精度从±0.5 mm提升至±0.05 mm，且换向平稳性显著改善。
复盘：摩擦模型的选择必须与实际速度范围匹配。对于频繁穿越零速的应用，静摩擦和Stribeck效应不可忽略。同时，静摩擦补偿脉冲的幅值和持续时间需要根据负载惯量仔细调整，过大或过小都会引起新的问题。