ISO 16940 标准深度解读2 | 用阻抗峰值识别共振频率？你可能一直在用"反共振"

这一误识问题在振动测试领域由来已久。Ewins 在经典著作《Modal Testing》中明确警告：阻抗 FRF 中的反共振峰与共振峰同样显著，极易被误识为模态共振，这是模态参数提取中最常见且后果最严重的错误之一。本文将从实验现象、理论曲线和有限元仿真三个角度，系统澄清这个问题。

一、实验现象

1.1 实验设置

本文测试样件为夹层玻璃，尺寸 $L = 300\ \text{mm}$，$b = 25\ \text{mm}$，$h = 6.38\ \text{mm}$（含 PVB 层），和裸钢条，尺寸 $L = 300\ \text{mm}$，$b = 20\ \text{mm}$，$h = 1.0\ \text{mm}$，样件自由-自由支撑，中心点激励，中心点测振动加速度。

测试仪器为聿钛科技的CDP型阻尼测试系统，可以采集激振点原点导纳 $v/F$ 和机械阻抗 $F/v$，通过 dB 标度输出。

1.2 阻抗曲线与导纳曲线的镜像对称

根据定义，导纳 $Y(\omega)$ 和阻抗 $Z(\omega)$ 的表达式为：

$$Y(\omega) = \frac{v(\omega)}{F(\omega)}, \qquad Z(\omega) = \frac{F(\omega)}{v(\omega)} = \frac{1}{Y(\omega)}$$

取 dB 标度，当参考值匹配（$Y_{\text{ref}} = 1$ (mm/s)/N，$Z_{\text{ref}} = 1$ N/(mm/s)）时：

$$Z_{\text{dB}} = -Y_{\text{dB}}$$

即阻抗曲线与导纳曲线关于 0 dB 完全对称镜像——这是倒数关系的纯数学推论，不涉及任何近似。

1.3 实验观察与困惑

对阻尼较大的夹胶玻璃（laminated glass）进行测试，分别导出实测的导纳曲线和阻抗曲线，并且根据半功率带宽法在曲线上识别峰值频率和损耗因子。

以第一阶为例：导纳曲线上识别的峰值频率 Fn1 = 224.5 Hz（$\eta$ = 0.0600），而阻抗曲线上识别的峰值频率（也就是导纳曲线上的谷值频率）Anti-1 = 139.0 Hz（$\eta$ = 0.0671）——两者频率相差 85.5 Hz。而损耗因子 $\eta$ 的值相近（0.0600 vs 0.0671），这正是"将错就错"现象的根源。

其余各阶对比同样明显：

关键观察：峰值频率的偏差随阶次升高而增大（85→499 Hz），而损耗因子的偏差始终较小（0.01~0.02 以内）。这说明阻抗法"看似能用"仅限于损耗因子，共振频率的识别是完全错误的。

二、简要理论分析

2.1 理论导纳曲线与模态对应

基于 Euler-Bernoulli 梁理论，可以推导出自由-自由梁中心点速度导纳的解析表达式：

$$Y(\omega) = j\omega \left[ \frac{-1}{\omega^2 M} + \sum_{n=3,5,7,\ldots}^{N} \frac{\phi_n^2(L/2)}{M(\omega_n^2 - \omega^2 + j\eta_n\omega_n^2)} \right]$$

其中 $M$ 为梁的总质量，$\omega_n$ 为第 $n$ 阶固有频率，$\phi_n$ 为模态振型，$\eta_n$ 为损耗因子。

模态频率计算

自由-自由梁的固有频率由特征方程确定：

$$\cos(\beta L) \cdot \cosh(\beta L) = 1$$

数值求解得到前几阶特征根 $\beta_n L$，进而计算固有频率：

$$f_n = \frac{(\beta_n L)^2}{2\pi L^2} \sqrt{\frac{EI}{\rho A}}$$

模态的对称性——决定中心点 FRF 形态的关键：

• 刚体模态（Mode 1, 2）这里不作展开
• 对称模态（Mode 3, 5, 7, …）在中心点振型不为零，被激发产生共振峰，在导纳曲线上为真正的共振峰
• 反对称模态（Mode 4, 6, 8, …）在中心点振型为零，不会被激发，对导纳曲线或阻抗曲线都无贡献

因此，中心点 FRF 只能看到奇数阶柔性模态的共振峰。

理论导纳曲线清晰地展示了以下特征：

• 共振峰仅出现在奇数阶固有频率处，峰值的中心频率与奇数阶固有频率一一对应
• 相邻共振峰之间存在反共振峰（共振谷，导纳极小值），频率与奇数阶和偶数阶固有频率都对不上
• 导纳曲线的谷（阻抗曲线上的峰）也并非自由梁的模态，在此识别模态参数没有理论依据

2.2 FEA 仿真验证：偏心激振如何激发偶数阶的反对称模态

前面的分析基于一个理想的中心点激振：激励力和响应测量点严格位于梁的正中心（$x = L/2$）。在这种理想情况下，偶数阶反对称模态在中心点的振型为零（$\phi_n(L/2) = 0$），因此对 FRF 没有任何贡献，中心点导纳曲线上只出现奇数阶对称模态的共振峰。

然而在实际测试和 FEA 仿真中，激励/响应点很难做到严格居中——即使是毫米级的偏差，也会改变这一结论。为了可以看到偶数阶的共振，FEA 仿真模型中，激励和响应点均不在严格的正中心，均偏差了半个网格的距离。

从图 7 可见：

• 奇数阶的对称模态 在导纳曲线上产生大幅值共振峰，因为它们在激振点响应点处于振型的波腹处
• 偶数阶反对称模态在导纳曲线上也可会产生可辨识的小峰，因为 FEA 模型中的激励/响应点与几何中心存在微小偏差

对本文主题的启示：

• FEA 结果进一步证实了中心点激振法的核心逻辑——当激振和响应点理想居中时，导纳曲线的峰确实是真正意义上的共振峰，且都是由奇数阶的对称模态所贡献，偶数阶的反对称模态的贡献被完全抑制。
• 导纳曲线的谷，常被称作反共振峰（导纳曲线上看是谷，对应阻抗曲线的峰），其频率值并不是偶数阶固有频率。因此，导纳曲线的谷（阻抗曲线的峰）并非真正意义的共振，并非奇数阶或偶数阶的模态所贡献。

2.3 导纳曲线上的谷：模态相消干涉

谷底既不对应任何对称模态，也不对应任何反对称模态。 它不是某个"隐藏模态"的体现，而是模态相消干涉（modal cancellation）的结果。

模态相消干涉的物理机制

驱动点导纳的模态叠加形式为：

$$Y(\omega) = j\omega \sum_{n} \frac{\phi_n^2(x_p)}{M(\omega_n^2 - \omega^2 + j\eta_n\omega_n^2)}$$

每一项的贡献系数 $\phi_n^2(x_p)$ 始终为正（振型平方）。但分母中的 $(\omega_n^2 - \omega^2)$ 项在跨过低阶模态的共振频率后符号翻转：

• 当 $\omega < \omega_n$：该项为正（刚度控制区）
• 当 $\omega > \omega_n$：该项为负（质量控制区）

因此，在两个相邻共振频率之间的某个频率处，较低阶模态的贡献为负、较高阶模态的贡献为正，二者恰好等量抵消，导致总响应出现极小值——这就是反共振（导纳曲线的谷底）。

数学上，反共振对应 FRF 传递函数的零点（而非极点）。它不是系统本身的固有振动模式，而是特定激励/响应位置处的局部响应特征。

阻抗曲线的峰值不对应任何模态的固有频率。 将阻抗峰值用于模态参数识别，本质上是用系统的"零点"去拟合"极点"，这是概念上的根本错误。

三、讨论：应该识别共振峰还是反共振峰？

3.1 阻抗曲线并非不能用——关键是用它的峰还是谷

由于 $Z = 1/Y$，$Z_{\text{dB}} = -Y_{\text{dB}}$，阻抗与导纳的峰谷完全互换。

关键结论：阻抗曲线可以用，但必须用它的谷（而非峰）来识别模态参数。

阻抗曲线的谷 = 导纳曲线的峰 = 共振频率（系统固有频率）。所以如果你手头只有阻抗数据，又不想翻转至导纳曲线，识别阻抗曲线的谷底同样可以得到正确的共振频率和损耗因子，这与识别导纳曲线的峰值是等价的。

3.2 为什么阻抗峰值识别"看起来也能用"？

1. 半功率带宽的数学性质：由于 $Z = 1/Y$ 的倒数关系，弱阻尼下反共振"谷"的形状与共振"峰"互为镜像，两者的 3 dB 带宽近似：

   $$\frac{\Delta f_{\text{res}}}{f_{\text{res}}} \approx \frac{\Delta f_{\text{anti}}}{f_{\text{anti}}} \approx 2\zeta$$

   因此即使测量的是反共振带宽，得到的 $\eta$ 值与共振带宽的 $\eta$ 接近。但这是反共振的带宽，不是共振的带宽。
2. 共振频率不被重视：部分应用只关心损耗因子，不关心频率绝对值，导致问题长期被忽视。

3.3 用阻抗峰值识别的实际危害

1. 有限元模型修正：用反共振频率作为目标值修正 FEA 模型，会导致刚度参数严重偏离真实值。
2. 高频模态密集区：共振峰与反共振峰交错出现，阻抗峰值导致模态识别完全混乱。
3. 标准化和可重复性：不同实验室采用不同方法，结果无法横向对比。
4. 物理概念混淆：将反共振误认为共振，是对系统动力学行为的根本误解。

四、结论与建议

1. 识别模态参数时，关键是识别有意义的峰或谷。导纳的峰和阻抗的谷都对应共振（系统固有频率），是正确的识别对象；阻抗的峰和导纳的谷都对应反共振（模态相消干涉），不是固有频率。

2. 用阻抗峰值识别损耗因子"看起来差不多"是弱阻尼下驱动点 FRF 的物理推论，而非偶然巧合。但共振频率的识别是完全错误的，物理含义也是错误的。

3. 推荐做法：使用导纳曲线的峰值来识别共振频率和损耗因子。这是最标准最合理的方法，也是最直观的方法。如果只能使用阻抗数据，则应识别阻抗曲线的谷底（而非峰值）。

附录 A：频率响应函数的理论推导

A.1 Euler-Bernoulli 梁振动方程

等截面均匀梁横向自由振动的控制方程：

$$\rho A \frac{\partial^2 w}{\partial t^2} + EI \frac{\partial^4 w}{\partial x^4} = f(x,t)$$

分离变量 $w(x,t) = \phi(x) e^{j\omega t}$，得空间部分四阶 ODE：

$$\frac{d^4\phi}{dx^4} - \beta^4 \phi = 0, \qquad \beta^4 = \frac{\rho A \omega^2}{EI}$$

频率与波数关系：

$$\omega = \beta^2 c, \qquad c = \sqrt{\frac{EI}{\rho A}}$$

A.2 自由-自由边界条件与特征方程

自由-自由梁两端弯矩和剪力为零：

$$\phi''(0) = \phi'''(0) = \phi''(L) = \phi'''(L) = 0$$

代入通解，导出特征方程：

$$\boxed{\cos(\beta L) \cdot \cosh(\beta L) = 1}$$

A.3 振型函数与对称性

柔性模态振型函数：

$$\phi_n(x) = [\cosh(\beta_n x) + \cos(\beta_n x)] - \gamma_n [\sinh(\beta_n x) + \sin(\beta_n x)]$$

A.4 模态叠加求位移导纳

在中心点 $x = L/2$ 施加单位力 $F = 1$，位移导纳（位移/力）由模态叠加给出：

$$H_{\text{disp}}(\omega) = \frac{-1}{\omega^2 M} + \sum_{n=3,5,7,\ldots}^{N} \frac{\phi_n^2(L/2)}{M(\omega_n^2 - \omega^2 + j\eta_n\omega_n^2)}$$

A.5 速度导纳与阻抗

速度导纳：

$$Y(\omega) = \frac{v(L/2)}{F} = j\omega \cdot H_{\text{disp}}(\omega)$$

机械阻抗为速度导纳的倒数：

$$Z(\omega) = \frac{1}{Y(\omega)}$$

dB 标度下：

$$Z_{\text{dB}} = 20\log_{10}|Z| = -20\log_{10}|Y| = -Y_{\text{dB}}$$

A.6 交错性质（Interlacing Property）

对于驱动点导纳（原点 FRF），共振与反共振严格交替出现——每两个相邻共振峰之间必有一个反共振谷。这是因为驱动点导纳的分子多项式阶数比分母低一阶，且所有零点均为实数，必然交错分布在极点之间。

这一性质可作为实验数据质量的检验标准：若实测驱动点 FRF 不呈现峰谷交替，则说明测量有误。

A.7 半功率带宽法的理论依据

对于弱阻尼系统，速度导纳在共振频率 $\omega_n$ 附近近似单自由度行为。半功率点定义为幅值下降到峰值的 $1/\sqrt{2}$ 时的频率 $\omega_1$ 和 $\omega_2$：

$$|Y(\omega_{1,2})| = \frac{|Y|_{\text{max}}}{\sqrt{2}}$$

损耗因子由半功率带宽给出：

$$\eta = 2\zeta = \frac{\omega_2 - \omega_1}{\omega_n} = \frac{f_2 - f_1}{f_n}$$