1. 项目背景与核心价值
电机机械共振是工业自动化领域长期存在的棘手问题。当电机驱动机械系统运行时，如果激励频率接近系统的固有频率，就会引发共振现象。这种现象轻则导致定位精度下降、表面加工质量恶化，重则直接造成机械结构损坏。我在某数控机床改造项目中就遇到过主轴电机在1800rpm附近剧烈振动的案例，当时导致加工件表面出现明显振纹，废品率飙升30%。
传统解决方案往往采用机械方式（如增加阻尼器、改变结构刚度）或简单降低转速避开共振点。但这些方法要么成本高昂，要么牺牲设备性能。而基于模型的控制算法设计能够在不改动硬件的前提下，通过软件手段有效抑制共振。Simulink作为控制系统仿真的事实标准工具，其可视化建模方式和丰富的电机库特别适合这类问题的快速验证。
这个仿真项目的核心价值在于：
- 提前在虚拟环境中验证控制算法有效性，避免实机调试风险
- 可快速对比PID控制、陷波滤波、自适应控制等不同方案的抑制效果
- 通过参数化建模实现不同型号电机的通用分析框架
- 为后续实际控制器代码生成提供可靠参考模型
2. 仿真模型构建要点
2.1 电机-负载系统建模
机械共振的本质是二阶振动系统，需要建立包含惯量、刚度、阻尼等参数的精确模型。在Simulink中推荐使用Simscape Multibody模块库构建物理模型，这比传统传递函数方法更能反映实际物理特性。关键参数包括：
- 电机转子惯量J_m（通常可在规格书找到）
- 负载惯量J_l（需实际测量计算）
- 轴系扭转刚度K（可通过频率响应测试获得）
- 粘性阻尼系数C（最难准确获取，建议先估算后修正）
经验提示：阻尼系数C的初始值可按0.01~0.05 N·m·s/rad估算，后期通过频响曲线拟合调整
2.2 共振特性分析方法
在施加控制算法前，必须先明确系统的共振特性。推荐采用扫频法进行频域分析：
- 构建开环测试模型，使用Chirp Signal模块生成0.1-100Hz线性扫频信号
- 通过PSD Estimator计算转速反馈信号的功率谱密度
- 使用Peak Finder定位共振峰对应的频率点
某400W伺服电机的实测频谱显示，在28.5Hz处出现明显共振峰，对应转速1710rpm。这个频率点将作为后续控制算法设计的重点抑制对象。
3. 主流抑制方案实现与对比
3.1 陷波滤波器设计
陷波滤波是最直接的解决方案，其传递函数为：
code复制 s² + ωₙ²
G(s) = ———————————
s² + 2ζωₙs + ωₙ²
在Simulink中可通过以下步骤实现：
- 使用Digital Notch Filter模块
- 设置中心频率ωₙ=2π×28.5（rad/s）
- 品质因数Q建议初始值设为5-10
- 通过Bode图验证滤波器在28.5Hz处有足够深度衰减
实测效果：共振峰值降低约12dB，但带来约5°的相位滞后，可能影响动态响应。
3.2 自适应陷波滤波进阶方案
固定参数的陷波滤波器在转速变化时效果会下降。采用LMS自适应算法可实现动态跟踪：
matlab复制function [y, e, w] = lms_notch(u, d, mu, order)
persistent w_hat;
if isempty(w_hat)
w_hat = zeros(order,1);
end
x = u:-1:u-order+1;
y = w_hat' * x';
e = d - y;
w_hat = w_hat + mu * e * x';
w = w_hat;
end
在Simulink中封装为MATLAB Function模块，配合Hilbert变换提取瞬时频率。实测表明该方法在±10%转速波动范围内保持良好抑制效果。
3.3 基于状态观测器的主动阻尼控制
更先进的方案是构建状态观测器估计谐振扭矩，然后前馈补偿：
- 设计降阶观测器估计负载端状态
- 计算谐振扭矩T_res = K(θ_m - θ_l)
- 在电流环给定中叠加补偿量ΔIq = T_res/Kt
关键参数关系：
code复制J_m·s²θ_m = K(θ_l - θ_m) + T_m
J_l·s²θ_l = K(θ_m - θ_l) - T_res
该方案在保持相同抑制效果的情况下，相位滞后减少到不足1°，特别适合高精度定位场景。
4. 仿真验证与结果分析
4.1 时域响应对比测试
设置转速指令从0阶跃到1800rpm（共振区），对比三种方案：
| 指标 | 无抑制 | 固定陷波 | 自适应陷波 | 主动阻尼 |
|---|---|---|---|---|
| 超调量(%) | 25.3 | 8.7 | 6.2 | 4.1 |
| 稳定时间(ms) | 420 | 380 | 350 | 310 |
| 残余振动(μm) | 15.2 | 5.8 | 3.3 | 1.2 |
4.2 频域特性验证
通过扫频测试获取闭环Bode图，重点关注：
- 共振峰衰减程度（至少-20dB）
- 相位裕度（保持>45°）
- 幅值裕度（保持>6dB）
实测数据表明主动阻尼方案在28.5Hz处实现-32dB抑制，同时相位裕度保持58°，综合性能最优。
5. 工程实施注意事项
-
参数敏感性分析：
- 刚度K误差±15%时，主动阻尼方案效果下降明显
- 建议配合在线参数辨识算法使用
-
离散化处理要点：
- 控制周期≤250μs（对应4kHz采样）
- 使用Tustin变换保持频率响应特性
- 注意量化误差积累问题
-
实机调试技巧：
- 先以1/10额定电流进行扫频测试
- 逐步提高滤波器强度观察效果
- 用听诊器辅助判断机械振动情况
某立式加工中心应用案例显示，采用该方案后：
- 表面粗糙度Ra从1.6μm提升到0.8μm
- 刀具寿命延长40%
- 极限进给速度提高15%