从频响函数中提取模态参数的理论基础

1. 多项式形式的频响函数

一般振动系统的频响函数的任意一个元素为：
$$H_{ij}(\omega )=\sum _{r=1}^N{Y}_r{\phi }_{ir}{\phi }_{jr}=\sum _{r=1}^N\frac{{\phi }_{ir}{\phi }_{jr}}{k_r-m_r{\omega }^2+j\omega c_r}$$
$H_{ij}$表达的物理含义是对j点施加冲击激励，i点的响应的傅里叶变换。因此，频响函数中的任意一行为：
[ H i 1 H i 2 ⋯ H i N ] = ∑ r = 1 N Y r [ φ i r { φ j r } T ] = ∑ r = 1 N [ r H i 1 r H i 2 ⋯ r H i N ] = ∑ r = 1 N Y r φ i r [ φ 1 r φ 2 r ⋯ φ N r ] = ∑ r = 1 N φ i r k r − ω 2 m r + j ω c r [ φ 1 r φ 2 r ⋯ φ N r ] \left[\begin{array}{llll} H_{i 1} & H_{i 2} & \cdots & H_{i N} \end{array}\right]=\sum_{r=1}^{N} Y_{r}\left[\varphi_{i r}\left\{\varphi_{j r}\right\}^{T}\right]=\sum_{r=1}^{N}\left[\begin{array}{llll} { }_{r} H_{i 1} & { }_{r} H_{i 2} & \cdots & { }_{r} H_{i N} \end{array}\right]\\=\sum_{r=1}^{N} Y_{r} \varphi_{i r}\left[\begin{array}{llll} \varphi_{1 r} & \varphi_{2 r} & \cdots & \varphi_{N r} \end{array}\right] \\=\sum_{r=1}^{N} \frac{\varphi_{i r}}{k_{r}-\omega^{2} m_{r}+j \omega c_{r}}\left[\begin{array}{llll} \varphi_{1 r} & \varphi_{2 r} & \cdots & \varphi_{N r} \end{array}\right] [Hi1​​Hi2​​⋯​HiN​​]=r=1∑N​Yr​[φir​{φjr​}T]=r=1∑N​[r​Hi1​​r​Hi2​​⋯​r​HiN​​]=r=1∑N​Yr​φir​[φ1r​​φ2r​​⋯​φNr​​]=r=1∑N​kr​−ω2mr​+jωcr​φir​​[φ1r​​φ2r​​⋯​φNr​​]
上式表明，频响函数为系统所有阶模态的叠加。并且，频响函数的任一行，都包含了该系统所有阶模态振型的信息。因此，我们可以通过激励各点，在i点拾振的方式，来得到频响函数的一行，从而通过曲线拟合的方式（例如，最小二乘原理）来得到各阶模态的振型。通过测量频响函数的一列来得到系统模态振型的方式，原理是一样的。

2. 部分分式形式的频响函数

除了采用多项式形式的频响函数，我们也可以采用部分分式形式表述的频响函数。对于一个振动系统，其传递函数为：
$$H(s)=\frac{1}{m{s}^2+cs+k}$$
该函数可以做如下形式的分解：
$$H(s)=\frac{1}{m(s-p)\left(s-p^{\ast }\right)}$$
其中，$p,p^{\ast }=-c/2m\pm \sqrt{(c/2m{)}^2-k/m}=-\xi {\omega }_n\pm j{\omega }_d$，上式可以进一步分解为部分分式的形式：
$$H(s)=\frac{r_1}{s-p}+\frac{r_1^{\ast }}{s-p^{\ast }}$$
其中，$r,r^{\ast }=\pm \frac{1}{2jm{\omega }_d}$。该数称为留数。令$s=j\omega$，则系统的频响函数为：
$$H(j\omega )=\frac{1}{2jm{\omega }_d(j\omega -p)}+\frac{-1}{2jm{\omega }_d\left(j\omega -p^{\ast }\right)}$$
对于多自由度系统，系统频响函数中的任意元素为：
$$H_{ij}(j\omega )=\sum _{r=1}^N\left(\frac{{}_r{A}_{ij}}{j\omega -p_r}+\frac{{}_r{A}_{ij}{}^{\ast }}{j\omega -p_r{}^{\ast }}\right)$$
完整的频响函数表达式为：
$$H_{ij}(j\omega )=\sum _{r=1}^N\left(\frac{\left[{}_{r}A\right]}{j\omega -p_r}+\frac{\left[{}_r{A}^{\ast }\right]}{j\omega -p_r{}^{\ast }}\right)$$
其中，${}_r{A}_{ij}=\frac{{\phi }_{ir}{\phi }_{jr}}{2j{m}_r{\omega }_{ch}},{}_r{A}_{ij}^{\ast }=-\frac{{\phi }_{ir}{\phi }_{jr}}{2j{m}_r{\omega }_{cb}}$。这虽然和多项式形式的频响函数表达式不同，但是本质上是一样的。
综上所述，我们只要测得了频响函数的任意一行或者一列，就能够从中得到系统的各阶模态，至于曲线拟合的方法，这就是后话了。