📚 课程目录第2章结构动力学基础与NVH建模语言

2. 结构动力学基础与NVH建模语言

在NVH工程中，我们面对的结构几乎都是连续体，但真正用于建模、计算和判断的，往往是一组被离散化后的“质量—弹簧—阻尼”系统。理解单自由度（SDOF）与多自由度（MDOF）系统，是读懂频响函数、模态叠加、阻抗与导纳这些NVH核心语言的前提。本节从最基本的动力学方程出发，逐步建立从物理坐标到模态坐标、从时域到频域的完整表达链条。

一、单自由度系统：NVH语言的最小单元

单自由度系统由一个质量块 m、一个弹簧 k 和一个阻尼器 c 组成，其运动方程为：

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

其中 x 为位移，f 为外力。对两边做傅里叶变换，令 s = jω，得到频域代数方程：

(-mω² + jcω + k)·X(ω) = F(ω)

由此定义位移频响函数（FRF）：

H(ω) = X(ω)/F(ω) = 1 / (k - mω² + jcω)

这个复数函数同时携带幅值和相位信息。其幅值在 ω_n = √(k/m) 附近出现峰值，峰值高度由阻尼比 ζ = c/(2√(km)) 控制。在NVH中，我们常关注三个量：固有频率 f_n = ω_n/(2π)、阻尼比 ζ、以及共振放大倍数 Q ≈ 1/(2ζ)。例如，一个 m=1 kg、k=39478 N/m、c=2.5 N·s/m 的系统，其 f_n≈31.6 Hz，ζ≈0.0063，共振峰非常尖锐，意味着该频率附近很小的激励就会引起大幅振动——这正是NVH中需要避开的危险区。

二、多自由度系统与模态叠加

实际结构离散为 N 个自由度后，运动方程变为矩阵形式：

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

其中 [M]、[C]、[K] 分别为质量、阻尼和刚度矩阵。直接求解耦合方程组计算量大且物理意义不直观。模态叠加法的核心思想是：通过求解无阻尼特征值问题

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

得到 N 个固有频率 ω_r 和对应的模态振型 {φ_r}。这些振型具有关于质量矩阵和刚度矩阵的正交性。将物理坐标 {x} 用模态坐标 {q} 展开：

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

代入原方程并左乘 [Φ]^T，利用正交性，耦合方程组解耦为 N 个独立的单自由度方程：

m_r·q̈_r + c_r·q̇_r + k_r·q_r = φ_r^T·{F(t)}

其中 m_r = φ_r^T[M]φ_r 为模态质量，k_r = φ_r^T[K]φ_r 为模态刚度，右端为模态力。这意味着：任何复杂结构的响应，都可以看作各阶模态独立响应的线性叠加。在NVH仿真中，通常只保留前若干阶对问题频段有贡献的模态，这就是“模态截断”。例如，分析0–200 Hz的车内噪声，若结构前10阶模态覆盖到250 Hz，则取前10阶即可，高阶模态对低频段贡献很小。

三、频响函数、阻抗与导纳

对解耦后的单自由度方程做频域变换，可得第 r 阶模态的频响函数。将各阶叠加，得到原点或跨点频响函数：

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

其中 H_jk 表示在 k 点激励、在 j 点响应的频响函数。这个公式是实验模态分析和仿真后处理的核心：频响函数的每个峰值对应一阶模态，峰值的形状给出阻尼，峰值的符号和相对大小给出振型信息。

与频响函数紧密相关的两个概念是机械阻抗和机械导纳。阻抗定义为力与速度之比：

Z(ω) = F(ω)/V(ω) = jω·F(ω)/X(ω) = jω / H(ω)

导纳则是阻抗的倒数，即 Y(ω) = 1/Z(ω) = V(ω)/F(ω)。在NVH中，阻抗常用于描述基础或安装点对振动的“抵抗能力”：阻抗越大，同样力下速度响应越小，隔振效果越好。导纳则常用于描述结构对激励的“敏感程度”。两者互为倒数，选择哪个取决于分析视角——从激励源看用阻抗，从响应点看用导纳。

四、从方程到工程判断：一个简化的NVH分析示例

假设某支架结构离散为3自由度系统，质量矩阵为对角阵 diag(2, 1.5, 1) kg，刚度矩阵为：

[K] = [ 8000  -4000    0
        -4000   8000  -4000
            0  -4000   8000 ] N/m

求解特征值可得三阶固有频率约为 f₁≈7.2 Hz、f₂≈19.8 Hz、f₃≈29.5 Hz。若激励源在第二自由度上施加10 N的简谐力，频率扫描范围0–40 Hz，则通过模态叠加计算各点响应。在 f₂ 附近，第二自由度响应最大，因为该阶模态在第二自由度的振型分量最大。若在该点布置加速度传感器，实测频响曲线应在19.8 Hz附近出现明显峰值。若峰值过高，可通过增加阻尼（如粘贴阻尼片）使 ζ 从0.01提升到0.05，共振幅值约下降为原来的1/5。这就是从动力学方程到NVH改进的完整逻辑。

建模提示：在NVH仿真中，模态叠加法的精度取决于模态截断阶数和阻尼模型。对于比例阻尼，模态解耦严格成立；对于非比例阻尼，需采用复模态或状态空间法。工程中常先用比例阻尼快速评估，再对关键频段做精细验证。

五、阻抗与导纳的工程选择

在隔振设计中，我们关心的是“基础阻抗”与“设备导纳”的匹配。若设备导纳大（即阻抗小），则设备对基础的反作用力小，隔振器容易发挥作用；若基础阻抗大，则基础对振动的吸收能力强，传递到地面的振动小。实际测试中，常用力锤激励、加速度计测量，通过 H = a/F 得到加速度导纳，再积分得到速度导纳或位移导纳。注意：不同导纳之间相差 jω 因子，比较时必须统一类型。

总结来说，SDOF给出最基本的共振与阻尼概念；MDOF通过模态叠加将复杂问题解耦；FRF是连接仿真与实验的桥梁；阻抗与导纳则从能量传递角度描述结构特性。掌握这四者之间的数学关系和物理意义，就掌握了NVH建模语言的核心语法。