3. 机械结构模态理论基础

模态试验与仿真验证的所有讨论，最终都要落到一组数学关系上：结构在什么频率下容易振动、振动时各点如何相对运动、能量以多快速度耗散。这一节把这三个问题用振动方程、频响函数和模态参数串起来，为后续锤击法测试、模态参数识别和仿真对标提供统一的数学语言。

单自由度系统：最小可用的振动模型

把一个质量块通过弹簧和阻尼器连接到固定基础上，就得到最基础的单自由度（SDOF）系统。设质量 m、刚度 k、粘性阻尼系数 c，外力 f(t) 作用在质量块上，位移 x(t) 以静平衡位置为原点、向下为正，则运动方程为：

m·ẍ(t) + c·ẋ(t) + k·x(t) = f(t)

两边除以 m，引入无阻尼固有频率 ω_n = √(k/m) 和阻尼比 ζ = c/(2√(km))，方程改写为标准形式：

ẍ(t) + 2ζω_n·ẋ(t) + ω_n²·x(t) = f(t)/m

这里 ω_n 的单位是 rad/s，工程上更常用赫兹表示固有频率 f_n = ω_n/(2π)。阻尼比 ζ 是无量纲量，金属结构通常在 0.001～0.02，橡胶隔振件可达 0.05～0.2。这两个参数加上振型，就是模态参数的完整内容。

对简谐激励 f(t) = F·e^jωt，稳态位移响应为 x(t) = X·e^jωt，代入方程消去时间项，得到位移频响函数：

H(ω) = X/F = 1 / (k - mω² + jωc)
     = 1 / [m(ω_n² - ω² + 2jζω_nω)]

频响函数是复函数，其幅值在 ω ≈ ω_n 附近出现峰值，峰值高度约为 1/(2ζk)；相位在共振区从 0° 快速过渡到 -180°。半功率点法正是利用幅值降到峰值的 1/√2 处对应的两个频率 ω₁、ω₂ 来估计阻尼比：ζ ≈ (ω₂ - ω₁)/(2ω_n)。

多自由度系统：从矩阵方程到模态叠加

真实机械结构是连续体，有限元离散后得到 N 个自由度的矩阵方程：

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

其中 [M]、[C]、[K] 分别为质量、阻尼、刚度矩阵，{x} 是节点位移向量。忽略阻尼时，令 {f} = 0 并假设 {x} = {φ}·e^jωt，得到广义特征值问题：

([K] - ω²[M])·{φ} = {0}

求解得到 N 个特征值 ω_r² 和对应的特征向量 {φ_r}，即第 r 阶固有频率和振型。振型是相对量，只描述各自由度振幅的比例和相位关系，需要归一化（常用质量归一化 {φ}^T[M]{φ} = 1）才能用于定量计算。

利用振型的正交性，可以把耦合的 N 自由度方程解耦成 N 个独立的单自由度方程。物理坐标 {x} 通过模态坐标 {q} 展开：

{x} = [Φ]·{q} = Σ φ_r · q_r

代入并左乘 [Φ]^T，得到解耦后的第 r 阶方程：

m_r·q̈_r + c_r·q̇_r + k_r·q_r = φ_r^T·{f}

其中 m_r、k_r、c_r 为模态质量、模态刚度、模态阻尼。由此可写出多自由度系统的频响函数矩阵：

H_jk(ω) = Σ_r [ φ_jr · φ_kr / (m_r(ω_r² - ω² + 2jζ_rω_rω)) ]

下标 j 表示响应点，k 表示激励点。当 j = k 时为原点频响，j ≠ k 时为跨点频响。这一表达式说明：任意一点的频响都是各阶模态贡献的叠加，靠近某阶固有频率时该阶模态主导响应。

模态参数与频响函数的对应关系

从上面的频响函数表达式可以看出，每一阶模态对频响的贡献由三个参数决定：ω_r 决定峰值在频率轴上的位置，ζ_r 决定峰值的尖锐程度和半功率带宽，{φ_r} 决定各测点之间的幅值比例和相位关系。反过来，在模态试验中，我们正是从实测频响曲线中提取这三个量：峰值定位固有频率，半功率带宽估算阻尼比，多测点拟合得到振型。

需要特别注意振型的归一化方式。质量归一化振型满足 {φ}_r^T[M]{φ}_r = 1，此时模态质量 m_r = 1，频响函数的分子直接是 φ_jr·φ_kr。如果采用最大位移归一化，则模态质量不再是 1，频响函数的幅值需要相应缩放。仿真与试验对标时，必须确认双方使用相同的归一化约定，否则振型相关系数（MAC）的计算会失真。

从方程到试验：一个可操作的对应关系

假设对一块铝制悬臂梁进行锤击试验，梁长 300 mm、截面 20 mm × 5 mm，材料弹性模量 70 GPa、密度 2700 kg/m³。用有限元或解析法可估算前三阶弯曲固有频率约为 82 Hz、515 Hz、1440 Hz（示例估算值，实际以试验为准）。在梁上布置 5 个测点，依次锤击每个测点、在固定参考点测量响应，得到 5 条跨点频响和 1 条原点频响。

对每条频响曲线，在 82 Hz 附近用半功率法读取 f₁ 和 f₂，计算 ζ = (f₂ - f₁)/(2f_n)。若测得 f₁ = 80.5 Hz、f₂ = 83.5 Hz，则 ζ ≈ 3/(2×82) ≈ 0.018，属于铝结构典型范围。将 5 个测点在该阶的幅值和相位组合起来，就得到第一阶弯曲振型——悬臂梁应呈现根部节点、端部最大幅值的形态。若某测点相位与相邻测点相差接近 180°，说明该点位于节点附近或测量方向接反，需要检查。

把试验得到的 f_r、ζ_r、{φ_r} 与仿真结果逐阶对比：频率误差通常要求控制在 5% 以内，阻尼比因仿真输入的不确定性可放宽，振型用 MAC 值判断，MAC 大于 0.9 认为相关性良好。这些判断标准的具体执行，将在后续模态试验与仿真验证章节中展开。