4、共振的物理与数学基础：机械谐振、固有频率、阻尼比、传递函数
要理解伺服系统中的共振抑制，首先必须从物理和数学层面弄清楚“共振”到底是什么。伺服系统由电机、联轴器、丝杠、负载等机械部件构成，这些部件并非刚体，而是具有质量、刚度和阻尼的弹性系统。当电机输出力矩驱动负载时，机械传动链会发生扭转形变，形成一个典型的“质量-弹簧-阻尼”二阶振荡系统。本章将从机械谐振的物理本质出发，推导其数学模型，并重点讲解固有频率、阻尼比和传递函数这三个核心概念，为后续的共振识别与滤波器设计打下坚实的理论基础。
4.1 机械谐振的物理本质
在伺服系统中，机械谐振源于传动链的有限刚度。以典型的“电机-联轴器-丝杠-工作台”结构为例，联轴器和丝杠并非刚性，它们在传递力矩时会发生扭转变形。当电机以某一特定频率输出交变力矩时，如果该频率与机械系统的固有频率接近，变形量会被反复放大，导致工作台或负载出现剧烈的振荡，这就是机械谐振。
从能量角度看，机械谐振是动能与势能之间的周期性转换。电机输入的能量在惯量（质量）与刚度（弹簧）之间来回交换，而阻尼则不断消耗能量。如果阻尼很小，能量交换效率高，谐振幅度就会很大。在伺服系统中，常见的谐振频率范围从几十赫兹到几百赫兹不等，具体取决于机械结构的尺寸、材料和连接方式。
4.2 固有频率与阻尼比：决定谐振特性的两个核心参数
对于一个标准的二阶机械系统，其运动方程可以表示为：
J·θ''(t) + C·θ'(t) + K·θ(t) = T(t)其中，J为等效转动惯量（kg·m²），C为粘性阻尼系数（N·m·s/rad），K为扭转刚度（N·m/rad），θ为角位移，T为驱动力矩。将上式改写为标准二阶系统形式：
θ''(t) + 2ζωₙ·θ'(t) + ωₙ²·θ(t) = (1/J)·T(t)由此可以定义两个最重要的参数：
- 固有频率（无阻尼自然频率）：ωₙ = √(K/J)，单位为 rad/s。它反映了系统在没有阻尼时自由振荡的频率。对应的频率值 fₙ = ωₙ / (2π) Hz。固有频率越高，系统响应越快，但越容易激发高频谐振。
- 阻尼比：ζ = C / (2√(J·K))，无量纲。它反映了系统耗散能量的能力。ζ < 1 时系统欠阻尼，谐振峰明显；ζ = 1 时临界阻尼，响应最快且无超调；ζ > 1 时过阻尼，响应缓慢。
在实际伺服系统中，阻尼比通常较小（0.01~0.1），因此谐振峰非常尖锐。这意味着在固有频率附近，微小的激励力矩就会引起很大的速度或位置波动。
| 参数 | 物理意义 | 对伺服系统的影响 | 典型取值范围 |
|---|---|---|---|
| 固有频率 ωₙ | 系统自由振荡的角频率 | 决定谐振峰的位置，影响速度环带宽上限 | 50~500 Hz（视机械结构而定） |
| 阻尼比 ζ | 能量耗散与临界阻尼之比 | 决定谐振峰的尖锐程度，影响系统稳定性裕度 | 0.01~0.2（金属传动链） |
| 刚度 K | 抵抗扭转变形的能力 | K 越大，固有频率越高，谐振越难被激发 | 10³~10⁶ N·m/rad |
| 惯量 J | 旋转运动的惯性度量 | J 越大，固有频率越低，响应越慢 | 10⁻⁴~10⁻¹ kg·m² |
4.3 传递函数：从时域到频域的桥梁
传递函数是分析线性时不变系统最有力的数学工具。对上述运动方程进行拉普拉斯变换（初始条件为零），可以得到从电机力矩到负载角速度的传递函数：
G(s) = ω(s) / T(s) = 1 / (J·s² + C·s + K)将分母写成标准形式：
G(s) = (1/J) / (s² + 2ζωₙ·s + ωₙ²)在频域分析中，令 s = jω，可以得到系统的频率响应特性。当 ω 接近 ωₙ 时，分母的模值达到最小值，此时系统增益达到最大值，即产生谐振峰。谐振峰的幅值约为：
|G(jωₙ)| ≈ 1 / (2ζ·K)这意味着阻尼比越小，谐振峰越高。例如，当 ζ = 0.05 时，谐振峰增益是低频增益的 10 倍（20 dB），这足以使伺服系统在谐振频率附近产生严重的振荡甚至失稳。
4.4 案例拆解：双惯量系统的谐振分析
为了更直观地理解上述概念，我们分析一个典型的双惯量伺服系统。该系统由电机（J₁ = 0.002 kg·m²）通过一个弹性联轴器（K = 200 N·m/rad，C = 0.02 N·m·s/rad）驱动负载（J₂ = 0.008 kg·m²）。
第一步：计算等效惯量。对于双惯量系统，等效惯量 J_eq = J₁·J₂ / (J₁ + J₂) = 0.002 × 0.008 / 0.01 = 0.0016 kg·m²。
第二步：计算固有频率。ωₙ = √(K / J_eq) = √(200 / 0.0016) = √125000 ≈ 353.6 rad/s，对应 fₙ ≈ 56.3 Hz。
第三步：计算阻尼比。ζ = C / (2√(J_eq·K)) = 0.02 / (2√(0.0016 × 200)) = 0.02 / (2 × 0.566) ≈ 0.0177。
第四步：分析谐振峰。由于 ζ 很小，谐振峰增益约为 1/(2×0.0177×200) ≈ 0.141，相对于低频增益（1/K = 0.005）放大了约 28 倍（约 29 dB）。这是一个非常严重的谐振，如果不加抑制，速度环增益稍高就会导致系统啸叫或振荡。
第五步：设计陷波滤波器。针对 56.3 Hz 的谐振频率，在速度环输出端串联一个陷波滤波器，中心频率设为 56 Hz，陷波深度 20 dB，品质因数 Q = 2。滤波器传递函数为：
H(s) = (s² + ωₙ²) / (s² + (ωₙ/Q)·s + ωₙ²)加入滤波器后，谐振峰被有效压制，速度环增益可以提升约 10 dB，系统带宽从 8 Hz 提升到 15 Hz，同时保持足够的相位裕度（>45°）。
4.5 常见误区与总结复盘
在理解和处理机械谐振时，工程师常犯以下错误：
- 误区一：认为谐振频率是固定不变的。实际上，随着负载变化、温度变化或机械磨损，刚度 K 和惯量 J 都会改变，谐振频率会漂移。因此，固定参数的陷波滤波器可能失效，需要采用自适应陷波或在线辨识技术。
- 误区二：盲目降低速度环增益来避免谐振。这虽然能抑制振荡，但会牺牲系统的动态响应和跟踪精度。正确的做法是先识别谐振频率，再用滤波器针对性抑制，然后尽可能提高增益。
- 误区三：忽略阻尼比的影响。有些工程师只关注固有频率，而忽略了阻尼比。实际上，阻尼比决定了谐振峰的尖锐程度，也决定了系统对参数变化的敏感度。阻尼比越小，系统越“脆”，越容易因微小扰动而振荡。
- 误区四：将机械谐振与电气噪声混淆。机械谐振通常表现为特定频率的振荡，且与速度或位置反馈相关；而电气噪声往往是宽频带的随机信号。通过FFT频谱分析可以清晰区分两者。
总结复盘：本章的核心在于建立“物理现象→数学模型→频域工具”的分析链路。机械谐振是弹性传动系统的固有属性，由固有频率和阻尼比两个参数完全描述。传递函数将时域微分方程转化为频域代数表达式，使我们能够直观地看到谐振峰的位置和幅值。在实际伺服调试中，测量频率响应、识别谐振峰、设计陷波滤波器是标准的三步流程。理解这些基础，是后续学习共振抑制算法（如陷波滤波器、低通滤波器、加速度反馈等）的前提。