📚 课程目录第2章结构动力学基础回顾

2. 结构动力学基础回顾

NVH仿真的本质，是把连续的车身、悬架、动力总成等结构离散成有限自由度系统，再求解其在激励下的动态响应。因此，单自由度（SDOF）与多自由度（MDOF）系统、模态叠加、频响函数和阻尼模型，构成了整车NVH分析的数学与物理基础。本节不重复推导，而是把这些概念还原为仿真中真正会用到的量：质量、刚度、阻尼、固有频率、振型、频响峰值与相位。

一、单自由度系统：一切动态特性的最小单元

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

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

无阻尼时固有频率为 ω_n = √(k/m)，单位rad/s；工程中更常用 f_n = ω_n/(2π)，单位Hz。阻尼比定义为 ζ = c/(2√(km))。这两个参数决定了系统的全部动态行为：f_n 决定共振位置，ζ 决定共振峰的尖锐程度。

以车身悬置上的一个动力总成简化模型为例：质量 m = 200 kg，悬置垂向刚度 k = 320 N/mm = 3.2×10⁵ N/m，阻尼系数 c = 1200 N·s/m。则：

ωn = √(3.2e5 / 200) = √1600 = 40 rad/s
fn = 40 / (2π) ≈ 6.37 Hz
ζ  = 1200 / (2·√(200·3.2e5)) = 1200 / (2·8000) = 0.075

这意味着该悬置系统在约6.4 Hz附近会出现明显共振，阻尼比7.5%属于典型的橡胶悬置水平。若激励频率接近6.4 Hz，位移放大因子约为 1/(2ζ) ≈ 6.7，即激励幅值被放大约6.7倍。这个估算在仿真前就能帮助判断风险频段。

二、多自由度系统与模态

实际车身、车架、悬架都是多自由度系统，运动方程写成矩阵形式：

[M]·Ẍ + [C]·Ẋ + [K]·X = F(t)

其中 [M]、[C]、[K] 分别为质量、阻尼、刚度矩阵，维度等于自由度数量。有限元模型中，一个白车身可能有数百万自由度，直接求解代价极高。模态分析通过求解广义特征值问题：

([K] − ω²[M])·φ = 0

得到特征值 ω_r² 和特征向量 φ_r，即第 r 阶固有频率和振型。振型具有正交性：φ_r^T[M]φ_s = 0（r≠s），这是模态叠加法的数学基础。

工程上通常只提取前几十到几百阶模态。以某轿车白车身为示例（合理假设）：一阶扭转约28 Hz，一阶弯曲约42 Hz，前围板局部模态约55 Hz。若发动机怠速二阶激励为25 Hz，则一阶扭转模态可能被激发，需要在仿真中重点关注。

三、模态叠加与频响函数

模态叠加法将物理坐标响应表示为各阶模态的线性组合：

X(ω) = Σ φr · qr(ω)

其中 q_r 为模态坐标。代入运动方程并利用正交性，可将耦合的MDOF系统解耦为一系列独立的SDOF方程。每个模态坐标满足：

q̈r + 2ζrωr·q̇r + ωr²·qr = φrᵀ·F / mr

由此得到频响函数（FRF），即单位激励下的响应：

H(ω) = Σ [φr·φrᵀ] / [mr·(ωr² − ω² + j·2ζrωrω)]

FRF是NVH仿真的核心输出。它同时包含幅值和相位信息：幅值峰值出现在 ω ≈ ω_r 附近，峰值高度与 1/(2ζ_r) 成正比；相位在共振区发生约180°翻转。仿真中常用FRF判断：某测点在20–200 Hz内是否有超限峰值、峰值对应哪阶模态、阻尼是否足够。

判断技巧：若FRF峰值尖锐且相位快速翻转，说明该阶模态阻尼偏小；若峰值宽缓，则阻尼较大或模态密集耦合。结合模态振型云图，可定位是哪个部件在主导该峰值。

四、阻尼模型：从比例阻尼到结构阻尼

阻尼是NVH仿真中最难准确建模的部分。常用模型包括：

比例阻尼（Rayleigh阻尼）：[C] = α[M] + β[K]。α控制低频阻尼，β控制高频阻尼。适用于模态阻尼比随频率单调变化的场景，是整车仿真最常用的工程近似。
模态阻尼：直接给每阶模态指定阻尼比 ζ_r，如车身焊点结构取2%–5%，橡胶悬置取5%–10%。模态叠加法中直接使用，简单高效。
结构阻尼（滞回阻尼）：用复刚度 K(1+jη) 表示，η为损耗因子。适合橡胶、粘弹性材料，频域分析中与模态阻尼等价。
粘性阻尼：与速度成正比，物理意义明确，但实际结构阻尼往往更接近结构阻尼。

在仿真中，阻尼参数通常通过试验模态分析（EMA）或工作模态分析（OMA）反推。例如，对白车身进行锤击试验，测得某阶模态半功率带宽 Δf = 1.2 Hz，中心频率 f_r = 42 Hz，则阻尼比 ζ = Δf/(2f_r) = 1.2/84 ≈ 1.4%。这个值可直接输入仿真模型。

五、仿真中的实际应用要点

在整车NVH仿真中，上述理论落地为几个关键操作。第一，模态提取范围要覆盖关注频段。若分析目标为20–200 Hz，通常提取到300 Hz以上，确保模态叠加时高频截断误差可控。第二，模态叠加法适用于线性系统，若存在非线性接触、大变形或摩擦，需改用直接积分法或非线性模态方法。第三，FRF计算时激励点与响应点要对应实际工况：方向盘振动对应方向盘上12点方向，座椅导轨振动对应导轨安装点，动力总成悬置对应悬置主被动端。

阻尼设置直接影响FRF峰值预测精度。工程经验是：白车身焊点结构模态阻尼比取2%–4%，带涂胶和内饰后取4%–8%，橡胶悬置取5%–12%，液压悬置低频可达15%以上。若仿真峰值明显高于试验，优先检查阻尼是否偏小；若峰值频率偏移，则检查刚度或质量分布。

注意：比例阻尼中α和β不能随意取值。若已知两阶模态的阻尼比ζ₁、ζ₂和对应频率ω₁、ω₂，可反解：α = 2ω₁ω₂(ζ₁ω₂−ζ₂ω₁)/(ω₂²−ω₁²)，β = 2(ζ₂ω₁−ζ₁ω₂)/(ω₂²−ω₁²)。若只给一个目标阻尼比，通常取β=2ζ/ω_n，α取小值或零。

掌握SDOF与MDOF的对应关系、模态叠加的数学逻辑、FRF的幅相含义以及阻尼模型的适用边界，才能在后续的整车模态、悬置匹配、路噪和怠速振动仿真中，正确设置参数并合理解读结果。这些基础概念不是孤立的公式，而是贯穿NVH仿真全流程的判断依据。