双惯量伺服系统机械谐振抑制与陷波滤波器设计

1. 双惯量伺服系统机械谐振问题解析

在工业伺服控制领域，双惯量系统是最常见的机械传动模型之一。这种系统由电机侧惯量（通常指电机转子）和负载侧惯量通过弹性联轴器或传动轴连接构成。我曾在某半导体设备公司的运动控制项目中，亲眼目睹这种系统在高速启停时产生的剧烈振动——当时我们测量到的振幅超过了安全阈值的3倍，直接导致定位精度从±5μm恶化到±50μm。

机械谐振的本质在于弹性元件（如联轴器）的刚度k与系统惯量J形成的二阶振荡环节。当控制频率接近系统的固有频率ω_n=√(k/J)时，就会引发能量积累。以一个典型的滚珠丝杠平台为例，假设电机惯量J_m=0.001kg·m²，负载惯量J_l=0.01kg·m²，轴刚度k=1000Nm/rad，其谐振频率计算如下：

code复制等效惯量 J_eq = J_m + J_l = 0.011 kg·m²
固有频率 ω_n = √(k/J_eq) ≈ 301.5 rad/s ≈ 48Hz

这个频率范围内的控制信号会被机械系统放大，造成速度环震荡。我在调试某型号贴片机时，就遇到过48Hz附近的指令会引发机械臂剧烈抖动的情况，这正是谐振频率匹配导致的典型现象。

2. 陷波滤波器的原理与设计方法

陷波滤波器（Notch Filter）是解决谐振问题的利器，它的频率响应就像在频谱上"挖"出一个深坑。我更喜欢把它比作音频工程师使用的降噪耳机——能够精准屏蔽特定频段的干扰。其传递函数的标准形式为：

code复制G(s) = (s² + ω_n²) / (s² + 2ζω_ns + ω_n²)

其中ω_n是需要抑制的谐振频率，ζ控制滤波器的带宽。在去年参与的机器人关节控制项目中，我们通过实验测得系统在52Hz处存在谐振峰，于是设计了如下参数的陷波器：

matlab复制wn = 2*pi*52;  % 谐振频率(rad/s)
zeta = 0.1;    % 阻尼系数
num = [1 0 wn^2];
den = [1 2*zeta*wn wn^2];
notch_tf = tf(num, den);

实际调试中有个关键经验：ζ取值需要权衡。太小会导致抑制带宽不足，太大则会过度影响相邻频率。我们通过扫频测试发现，对于大多数工业伺服系统，ζ=0.1~0.3是最佳范围。某次在数控机床改造中，我们意外发现