4. 模态参数与模态空间

在完成频响函数测量之后，我们得到的是一组以频率为自变量的复数曲线。这些曲线本身只是“现象”，真正描述结构内在动力学特性的，是一组模态参数。本节要建立的，就是“测量数据—模态参数—物理坐标—模态坐标”之间的完整映射关系。理解这套映射，是后续做模态参数识别、模型修正和响应预测的基础。

一、五个核心模态参数的定义

对任意一个具有 N 个自由度的线性时不变结构，其动力学方程可写为：

[M]{ẍ} + [C]{ẋ} + [K]{x} = {f(t)}

其中 [M]、[C]、[K] 分别为质量、阻尼、刚度矩阵，{x} 为物理坐标下的位移向量。对无阻尼系统求解广义特征值问题 [K]{φ} = ω²[M]{φ}，可得到 N 个特征值 ω_r² 和对应的特征向量 {φ}_r。这组解直接引出五个模态参数：

• 固有频率 ω_r（或 f_r = ω_r/2π）：第 r 阶模态自由振动时的角频率，由 ω_r² = k_r/m_r 决定，是模态的“时间尺度”。
• 模态振型 {φ}_r：第 r 阶模态下各测点位移的相对比例关系，是一个只保留形状、不含绝对幅值的向量。它描述“结构以什么形状振动”。
• 模态阻尼比 ζ_r：第 r 阶模态的能量耗散程度，决定共振峰的尖锐程度与自由振动衰减速率。工程金属结构常在 0.1%～2%，橡胶隔振件可达 5%～20%。
• 模态质量 m_r：第 r 阶模态的等效惯性量，满足 {φ}_r^T[M]{φ}_r = m_r。它并非结构总质量，而是该阶振型下的“参与质量”。
• 模态刚度 k_r：满足 {φ}_r^T[K]{φ}_r = k_r，与模态质量通过 k_r = ω_r²m_r 严格绑定。

需要特别强调：模态质量与模态刚度不是唯一确定的，它们随振型的归一化方式而变。若将振型按最大位移归一化为 1，则 m_r、k_r 取某一组值；若按模态质量归一化（使 m_r = 1），则 k_r = ω_r²。因此工程报告中必须注明归一化约定，否则模态质量数值没有可比性。但无论怎样归一化，ω_r、ζ_r 和振型的形状比例是不变的，这才是模态参数的“物理内核”。

二、模态坐标变换：从耦合到解耦

物理坐标 {x} 下，N 自由度系统的 [M]、[C]、[K] 通常是耦合的——某一个测点的运动同时受多个模态影响。模态分析的核心思想，是用振型矩阵 [Φ] = [{φ}₁, {φ}₂, …, {φ}_N] 作为变换基，引入模态坐标 {q}：

{x(t)} = [Φ]{q(t)} = Σ<sub>r=1</sub><sup>N</sup> {φ}<sub>r</sub> q<sub>r</sub>(t)

这里 q_r(t) 就是第 r 阶模态坐标，表示该阶模态在总响应中的“参与量”。把该式代入动力学方程，并左乘 [Φ]^T，利用振型关于 [M]、[K] 的正交性（{φ}_r^T[M]{φ}_s = 0，r ≠ s），可得一组完全解耦的单自由度方程：

m<sub>r</sub> q̈<sub>r</sub> + c<sub>r</sub> q̇<sub>r</sub> + k<sub>r</sub> q<sub>r</sub> = {φ}<sub>r</sub><sup>T</sup>{f} = F<sub>r</sub>(t)

其中 c_r 为模态阻尼，F_r 为模态力。这一步的意义在于：原本 N 个互相耦合的方程，被拆成 N 个独立的一阶/二阶方程，每个方程只含一个模态参数。这就是“模态空间”的实质——在模态坐标下，结构表现为一组互不干扰的单自由度振子。

三、物理坐标到模态坐标的映射关系

实际测量中，我们关心的是“在某个物理测点 i 施加激励、在测点 j 拾取响应”的频响函数 H_ij(ω)。利用模态叠加，它可以写成：

H<sub>ij</sub>(ω) = Σ<sub>r=1</sub><sup>N</sup> [ φ<sub>ir</sub> φ<sub>jr</sub> / ( m<sub>r</sub>(ω<sub>r</sub>² − ω² + j2ζ<sub>r</sub>ω<sub>r</sub>ω) ) ]

这个公式是整门课程最关键的桥梁之一。它说明：

• 物理坐标下的频响函数，是各阶模态贡献的线性叠加；
• 每一阶的贡献由“振型乘积 φ_irφ_jr”决定，即激励点与响应点在该阶振型上的幅值乘积；
• 分母只含该阶的 ω_r、ζ_r、m_r，与测点无关。

由此可以反推：如果某阶振型在测点 i 或 j 上恰好有节点（φ ≈ 0），那么该阶在这一对测点的频响中几乎不出现。这正是模态实验中“激励点不能选在振型节点上”的原因。反过来，若想突出某一阶，应把激励点和响应点都布置在该阶振型幅值较大的位置。

四、一个可复现的示例

设一悬臂梁被离散为 3 个自由度，已通过有限元或实验得到前两阶模态参数（示例数据，非实测）：

现在要预测在测点 1 激励、测点 3 拾取的频响 H₃₁(ω)。按公式逐阶计算：

第1阶贡献：φ<sub>31</sub>φ<sub>11</sub> = 1.00 × 0.42 = 0.42
第2阶贡献：φ<sub>32</sub>φ<sub>12</sub> = 1.00 × (−0.95) = −0.95

H<sub>31</sub>(ω) = 0.42 / (ω<sub>1</sub>² − ω² + j2ζ<sub>1</sub>ω<sub>1</sub>ω)
            + (−0.95) / (ω<sub>2</sub>² − ω² + j2ζ<sub>2</sub>ω<sub>2</sub>ω)

在 ω = ω₁ = 2π×12.5 ≈ 78.5 rad/s 处，第一项分母实部为零，只剩虚部 j2ζ₁ω₁²，幅值约为 0.42/(2×0.012×78.5²) ≈ 0.028；第二项因 ω 远离 ω₂ 而贡献很小。因此 H₃₁ 在 12.5 Hz 附近出现明显峰值，且相位在共振点附近快速翻转约 180°。若把测点 1 换成测点 2，则第一阶贡献变为 φ₃₂φ₂₂ = 1.00×0.78 = 0.78，峰值幅值随之增大——这正说明频响函数的幅值直接反映振型乘积。

五、模态空间的工程意义

把物理坐标映射到模态坐标，带来的实际好处有三点。第一，降阶：多数工程问题中，前 10～20 阶模态已覆盖 90% 以上的响应能量，只需保留这些模态即可高精度重建频响，计算量远小于完整有限元模型。第二，解耦：每个模态可独立修改。若某阶共振峰过高，只需针对该阶调整 ζ_r 或改变 k_r（例如加筋、加阻尼层），不必重新求解整个耦合系统。第三，可诊断：当实测频响与预测不符时，可逐阶比对 ω_r、ζ_r 和振型乘积，定位是频率偏移、阻尼异常还是某测点振型符号错误。