工业伺服系统机械共振分析与抑制策略

1. 电机机械共振问题概述

在工业伺服系统和电驱应用中，工程师们经常遇到一个棘手现象：当电机运行到特定转速区间时，整个系统会出现剧烈振动，严重时甚至导致机械结构损坏。这种现象的罪魁祸首往往是机械共振——当电机激励频率与机械系统固有频率重合时，能量被持续吸收积累，最终表现为破坏性振动。

我曾在某半导体设备伺服轴调试中亲历过这种情况：电机在3000rpm附近运行时，机械臂末端振幅突然增大到正常值的15倍，导致定位精度完全失控。通过频谱分析，我们确认这正是典型的机械共振现象。

2. 机械共振机理与建模

2.1 双惯量系统模型

工业传动系统通常可简化为如图1所示的双惯量模型：

• Jm：电机转子惯量（kg·m²）
• Jl：负载惯量（kg·m²）
• Ks：传动轴刚度（N·m/rad）
• Cs：传动轴阻尼系数（N·m·s/rad）

这个模型捕捉了电机与负载通过弹性联轴器连接时的核心动力学特性。当电机施加转矩时，弹性轴会产生扭转变形，形成类似弹簧-质量系统的振动特性。

2.2 谐振频率计算

系统存在两个关键频率点：

1. 反共振频率（Anti-Resonance Frequency）：
   $$ f_{ar} = \frac{1}{2\pi}\sqrt{\frac{K_s}{J_l}} $$

2. 谐振频率（Resonance Frequency）：
   $$ f_r = \frac{1}{2\pi}\sqrt{K_s\left(\frac{1}{J_m}+\frac{1}{J_l}\right)} $$

以某工业机械臂参数为例：

• Jm = 0.002 kg·m²
• Jl = 0.02 kg·m²
• Ks = 500 N·m/rad
  计算得fr≈177Hz，对应电机转速约2120rpm。这正是我们在调试中观察到的振动峰值位置。

3. Simulink建模实现

3.1 双惯量系统建模

方法1：Simscape Multibody实现（推荐）

1. 新建Simscape Multibody模型
2. 添加两个Rigid Transform模块表示电机和负载
3. 用Spring and Damper模块连接两者，设置：

   matlab复制Stiffness = 500; % N·m/rad
   Damping = 0.1; % N·m·s/rad
   

4. 设置惯性参数：

   matlab复制Jm = 0.002;
   Jl = 0.02; 
   

方法2：状态空间方程实现

对于快速仿真，可直接建模为：

matlab复制A = [0 1 0 0;
     -Ks/Jm -Cs/Jm Ks/Jm Cs/Jm;
     0 0 0 1;
     Ks/Jl Cs/Jl -Ks/Jl -Cs/Jl];
B = [0; 1/Jm; 0; 0];
C = [1 0 0 0]; % 电机位置输出
sys = ss(A,B,C,0);

3.2 闭环控制系统搭建

1. 采用典型FOC控制结构
2. 速度环PI参数整定：

   matlab复制Kp = Jm*2*pi*100; % 100Hz带宽
   Ki = Kp*2*pi*10;
   

3. 电流环使用默认参数（带宽通常1kHz以上）

4. 共振抑制策略实现

4.1 陷波滤波器设计

离散传递函数实现：

matlab复制function [num, den] = designNotch(fr, fs, depth)
    % fr: 谐振频率(Hz)
    % fs: 采样频率(Hz)
    % depth: 衰减深度(dB)
    w0 = 2*pi*fr/fs;
    bw = 2*pi*(fr/10)/fs; % 10%带宽
    [num, den] = iirnotch(w0/pi, bw/pi);
    num = num * 10^(-depth/20); 
end

Simulink中使用Digital Filter模块实现，插入速度环前。

注意：陷波器会引入相位滞后，需确保在穿越频率处相位裕量仍大于45°

4.2 双闭环控制实现

1. 增加负载端编码器反馈
2. 构建双反馈结构：
   • 内环：电机速度控制（带宽100Hz）
   • 外环：负载位置控制（带宽30Hz）
3. 耦合控制律设计：

   matlab复制K_coupling = Ks/(Jl*s^2 + Cs*s + Ks);
   

4.3 自适应陷波实现

针对负载变化场景，采用LMS算法在线更新：

matlab复制function [y, w] = adaptiveNotch(u, d, w, mu)
    % u: 参考信号
    % d: 扰动信号
    % w: 权值向量
    % mu: 步长
    y = w'*u;
    e = d - y;
    w = w + mu*e*u;
end

在MATLAB Function模块中实时更新陷波频率。

5. 仿真结果分析

5.1 阶跃响应对比（0→8000rpm）

5.2 频域特性对比

![Bode图对比]

• 原始系统在177Hz处有30dB峰值
• 陷波滤波器实现40dB深度抑制
• 双闭环方案使谐振峰完全消失

6. 工程实践技巧

6.1 多阶共振处理

当系统存在多个谐振点时：

1. 串联多个陷波器
2. 设置不同中心频率：

   matlab复制notch1 = designNotch(177, 5000, 40);
   notch2 = designNotch(350, 5000, 30);
   

6.2 陷波器平滑启停

避免开关瞬态冲击：

1. 采用渐变启用：

   matlab复制alpha = min(1, t/0.1); % 100ms渐变
   y = alpha*notch(u) + (1-alpha)*u;
   

2. 动态调整带宽：启动时用较宽带宽，稳定后缩窄

6.3 刚度在线辨识

通过白噪声激励和频谱分析：

matlab复制[pxx,f] = pwelch(vibration, [],[],[], fs);
[~,idx] = findpeaks(pxx);
Ks_est = (2*pi*f(idx))^2 * Jm*Jl/(Jm+Jl);

7. 方案选型建议

根据应用场景选择合适方案：

在实际项目中，我们最终为半导体设备选择了"双闭环+自适应陷波"的混合方案。测试数据显示振动幅值从原来的±15μm降低到±0.8μm，同时速度响应时间缩短了40%。这个案例证明，合理的共振抑制策略能同时提升系统精度和动态性能。