振动与动力学基础 — 单自由度到随机振动完整解析

作者：NovaSolver Contributors (Anonymous Engineers & AI)  |  分类：基础理论  |  日文版 →  |  英文版 →

🧑‍🎓

老师，我做汽车车身模态分析，报告要给出前10阶固有频率和振型。这些数据有什么工程意义？

🎓

模态分析给出的固有频率和振型，是结构振动的"DNA"。工程意义很直接：第一，避开激励频率——发动机振动频率、路面激励频率等不能和车身固有频率重合，否则共振会引起噪声和疲劳。第二，振型告诉你哪里动得最多——这就是加强筋和减振块应该放的位置。

1. 振动基本概念——固有频率与阻尼

1.1 固有角频率

一个弹簧-质量系统（弹簧刚度 $k$，质量 $m$）的无阻尼固有角频率：

$$\omega_n = \sqrt{\frac{k}{m}}$$

固有频率（Hz）：$f_n = \omega_n / (2\pi)$，固有周期：$T_n = 1/f_n$。

直观理解：质量越大、刚度越小 → 固有频率越低（"笨重又软"的系统振动慢）。汽车悬架固有频率约1～2 Hz，工程结构桥梁约0.5～3 Hz，电子芯片封装可达数百 Hz。

1.2 阻尼与阻尼比

实际系统总有阻尼，用阻尼比 $\zeta$ 衡量阻尼程度：

$$\zeta = \frac{c}{c_{cr}} = \frac{c}{2\sqrt{km}} = \frac{c}{2m\omega_n}$$

$c$ 是粘性阻尼系数，$c_{cr}$ 是临界阻尼系数。阻尼类型：

阻尼比 ζ运动状态物理描述工程举例
ζ = 0无阻尼永远振动不衰减理论情况
0 < ζ < 1欠阻尼振幅逐渐衰减金属结构ζ≈1%～5%，橡胶ζ≈10%
ζ = 1临界阻尼最快返回平衡不振荡精密仪器返回零位
ζ > 1过阻尼缓慢返回平衡液压减振器（某些条件）

2. 单自由度（SDOF）系统详细分析

🧑‍🎓

SDOF系统的运动方程是什么？为什么要从这个最简单的情况入手？

🎓

SDOF是理解所有振动问题的基石。它的运动方程是：$m\ddot{x} + c\dot{x} + kx = F(t)$。虽然简单，但所有重要概念——共振、阻尼、频率响应——都能在SDOF里清晰地看到。多自由度系统经过模态分解，本质上也是一堆SDOF的叠加。

2.1 运动方程

$$m\ddot{x} + c\dot{x} + kx = F(t)$$

无量纲化（除以 $m$，用 $\omega_n$ 和 $\zeta$ 表示）：

$$\ddot{x} + 2\zeta\omega_n\dot{x} + \omega_n^2 x = \frac{F(t)}{m}$$

2.2 自由振动响应（欠阻尼）

有阻尼固有角频率：$\omega_d = \omega_n\sqrt{1-\zeta^2}$

自由振动响应（初位移 $x_0$，初速度 $v_0$）：

$$x(t) = e^{-\zeta\omega_n t}\left(x_0\cos\omega_d t + \frac{v_0 + \zeta\omega_n x_0}{\omega_d}\sin\omega_d t\right)$$

从衰减振动的包络中，可以用对数衰减率实验测定阻尼比：

$$\delta = \ln\frac{x_n}{x_{n+1}} = \frac{2\pi\zeta}{\sqrt{1-\zeta^2}} \approx 2\pi\zeta \quad\text{（小阻尼时）}$$

2.3 谐波激励下的稳态响应——动力放大因子

激励 $F(t) = F_0\sin\Omega t$，稳态响应的幅值：

$$X = \frac{F_0/k}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}, \quad r = \frac{\Omega}{\omega_n}$$

动力放大因子（DAF）= $X/(F_0/k)$，在共振时（$r \approx 1$）：

$$\text{DAF}_{\max} \approx \frac{1}{2\zeta} \quad\text{（共振放大倍数）}$$

阻尼比越小，共振峰越高越尖锐。钢结构 $\zeta=0.02$，共振放大约25倍；橡胶 $\zeta=0.1$，放大仅5倍——这是减振设计中加橡胶垫的物理依据。

3. 多自由度（MDOF）系统与特征值分析

3.1 MDOF运动方程

$$\mathbf{M}\ddot{\mathbf{u}} + \mathbf{C}\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} = \mathbf{f}(t)$$

$\mathbf{M}$、$\mathbf{C}$、$\mathbf{K}$ 分别是质量矩阵、阻尼矩阵和刚度矩阵，维数为 $n \times n$（$n$ = 自由度数，FEM模型可达数百万）。

3.2 无阻尼自由振动——广义特征值问题

自由振动 ($\mathbf{f}=0$，$\mathbf{C}=0$) 假设解为 $\mathbf{u} = \boldsymbol{\phi}e^{i\omega t}$：

$$(\mathbf{K} - \omega^2\mathbf{M})\boldsymbol{\phi} = \mathbf{0}$$

求解得 $n$ 对特征值-特征向量 $(\omega_i^2, \boldsymbol{\phi}_i)$，即固有频率和振型。振型关于质量矩阵和刚度矩阵正交：

$$\boldsymbol{\Phi}^T \mathbf{M} \boldsymbol{\Phi} = \mathbf{I}, \quad \boldsymbol{\Phi}^T \mathbf{K} \boldsymbol{\Phi} = \text{diag}(\omega_1^2, \omega_2^2, \ldots, \omega_n^2)$$

$\boldsymbol{\Phi} = [\boldsymbol{\phi}_1, \boldsymbol{\phi}_2, \ldots, \boldsymbol{\phi}_n]$ 是振型矩阵。

3.3 模态坐标变换（模态叠加法）

令 $\mathbf{u} = \boldsymbol{\Phi}\mathbf{q}$，代入运动方程并利用正交性，解耦为 $n$ 个独立的 SDOF 方程：

$$\ddot{q}_i + 2\zeta_i\omega_i\dot{q}_i + \omega_i^2 q_i = \boldsymbol{\phi}_i^T \mathbf{f}(t), \quad i = 1, 2, \ldots, n$$

物理响应：$\mathbf{u}(t) = \sum_{i=1}^{m} \boldsymbol{\phi}_i q_i(t)$（只需取前 $m$ 个主导模态，$m \ll n$，称为模态截断）。这是频率响应分析和随机振动分析的基础。

🧑‍🎓

实际仿真中，模态分析要取多少阶模态才够？

🎓

经验法则：提取的最高阶模态频率应达到关注频率范围上限的1.5～2倍。比如做振动测试到200 Hz，模态就取到300～400 Hz。另外看模态有效质量——各阶模态有效质量之和应达到总质量的80%以上，这样才能保证截断误差可控。实际工程中，汽车NVH分析通常取0～200 Hz范围内所有模态，可能有几百阶。

3.4 振型相关系数 MAC

模态振型一致性系数（Modal Assurance Criterion）用于比较两个振型（仿真vs测试，或细化前后）的相似性：

$$MAC(\boldsymbol{\phi}_A, \boldsymbol{\phi}_B) = \frac{|\boldsymbol{\phi}_A^T \boldsymbol{\phi}_B|^2}{(\boldsymbol{\phi}_A^T \boldsymbol{\phi}_A)(\boldsymbol{\phi}_B^T \boldsymbol{\phi}_B)} \in [0, 1]$$

MAC = 1.0 表示两振型完全一致，MAC < 0.8 说明振型差异较大，需要检查边界条件或模型精度。

4. 数值时间积分法

🧑‍🎓

Abaqus里的瞬态动力分析，有"隐式"和"显式"两种，怎么选？

🎓

简单判断：慢速/中速动力问题选隐式（振动、地震），快速冲击碰撞选显式。汽车碰撞试验（持续约100ms，含大变形和接触）是显式的典型应用；发动机振动、模态分析是隐式的典型应用。显式每步计算便宜但时间步很小，适合短时高速事件；隐式每步贵但步长可以大，适合低频长时分析。

4.1 Newmark-β法（隐式）

最广泛使用的隐式时间积分格式：

$$\dot{\mathbf{u}}_{n+1} = \dot{\mathbf{u}}_n + \Delta t\left[(1-\gamma)\ddot{\mathbf{u}}_n + \gamma\ddot{\mathbf{u}}_{n+1}\right]$$
$$\mathbf{u}_{n+1} = \mathbf{u}_n + \Delta t\dot{\mathbf{u}}_n + \Delta t^2\left[\left(\frac{1}{2}-\beta\right)\ddot{\mathbf{u}}_n + \beta\ddot{\mathbf{u}}_{n+1}\right]$$

常用参数设置：

  • $\gamma=0.5$，$\beta=0.25$（平均加速度法）：无条件稳定，二阶精度，无数值阻尼
  • $\gamma=0.5$，$\beta=0$（线性加速度法）：条件稳定（$\Delta t < T_n/\pi$）
  • $\gamma > 0.5$（如 $\gamma=0.6$，$\beta=0.3025$）：引入高频数值阻尼，滤除噪声

4.2 HHT-α法（隐式，改进型）

Hilber-Hughes-Taylor（HHT-α）法在Newmark基础上引入参数 $\alpha \in [-1/3, 0]$，对动力平衡方程加权：

$$\mathbf{M}\ddot{\mathbf{u}}_{n+1} + (1+\alpha)\mathbf{C}\dot{\mathbf{u}}_{n+1} - \alpha\mathbf{C}\dot{\mathbf{u}}_n + (1+\alpha)\mathbf{K}\mathbf{u}_{n+1} - \alpha\mathbf{K}\mathbf{u}_n = (1+\alpha)\mathbf{f}_{n+1} - \alpha\mathbf{f}_n$$

$\alpha < 0$ 时引入高频数值阻尼，同时保持低频精度。Abaqus隐式动力分析默认使用HHT-α（$\alpha = -0.05$）。

4.3 中心差分法（显式）

加速度用中心差分近似：

$$\ddot{\mathbf{u}}_n = \frac{\mathbf{u}_{n+1} - 2\mathbf{u}_n + \mathbf{u}_{n-1}}{\Delta t^2}$$

每步不需要求解线性方程组（用集中质量矩阵使 $\mathbf{M}$ 对角化），每步计算量极小。但有条件稳定，时间步必须小于临界值：

$$\Delta t_{\text{cr}} = \frac{2}{\omega_{\max}} \approx \frac{L_e}{c_d}$$

$L_e$ 是最小单元尺寸，$c_d$ 是波速。碰撞仿真中，需要1000步以上，每步约0.01 ms，总计约10 ms时间历程。

5. 频率响应分析（FRF）

🧑‍🎓

我看振动测试报告里有"FRF"（频率响应函数），它和模态分析有什么关系？

🎓

FRF是系统对单位谐波激励的稳态响应，是模态参数的"外在表现"——FRF的峰值对应固有频率，峰的宽窄对应阻尼比，峰的形状包含振型信息。模态分析就是从FRF曲线里"识别"出模态参数的过程。实验模态分析（EMA）正是通过锤击试验或激振器激励，测量FRF，再用曲线拟合算法提取模态。

5.1 频响函数的定义

在频域内，位移频响函数（柔度FRF）定义为：

$$H_{jk}(\omega) = \frac{U_j(\omega)}{F_k(\omega)} = \sum_{r=1}^{n}\frac{\phi_{jr}\phi_{kr}}{\omega_r^2 - \omega^2 + 2i\zeta_r\omega_r\omega}$$

$j$ 是响应自由度，$k$ 是激励自由度，$r$ 是模态阶次。FRF的虚部（Imaginary part）在共振频率处有峰值，符号变化对应振型的相位关系。

5.2 FRF类型

FRF类型物理量分子分母斜率（远离共振）
柔度（Receptance）位移/力$U$$F$-40 dB/decade
导纳（Mobility）速度/力$V$$F$-20 dB/decade
加速度导纳（Inertance）加速度/力$A$$F$0 dB/decade（平坦）

6. 随机振动——PSD与RMS

🧑‍🎓

做航空航天产品的振动环境试验，规范里给的是PSD曲线，这是什么？和我们平时说的振幅有什么关系？

🎓

功率谱密度（PSD）描述振动能量如何分布在各频率上，单位是 g²/Hz 或 (m/s²)²/Hz。简单说：PSD曲线下的面积 = 总均方值（均方加速度），开方就是RMS值（均方根）。比如你说振动量级是"1 g RMS"，就是指均方根加速度值。这比峰值更能代表随机振动的"持续强度"。

6.1 功率谱密度（PSD）

随机过程 $x(t)$ 的自功率谱密度定义为其自相关函数的傅里叶变换：

$$S_{xx}(\omega) = \int_{-\infty}^{\infty} R_{xx}(\tau) e^{-i\omega\tau} d\tau$$

总均方值（Parseval定理）：

$$\sigma_x^2 = \overline{x^2} = \int_{-\infty}^{\infty} S_{xx}(\omega) \frac{d\omega}{2\pi} = \int_0^{\infty} G_{xx}(f) df$$

RMS值：$x_{RMS} = \sqrt{\sigma_x^2}$，峰值设计通常用 $3\sigma$（对高斯过程，约99.73%的时间内不超过这个值）。

6.2 SDOF系统的随机响应

基础激励 PSD 为 $S_{gg}(f)$（g²/Hz），响应加速度 PSD：

$$S_{aa}(f) = |H(f)|^2 \cdot S_{gg}(f)$$

$|H(f)|^2$ 是功率传递函数。对于白噪声激励（$S_{gg} = \text{const}$）的SDOF系统，响应RMS加速度：

$$\sigma_a = \omega_n \sqrt{\frac{\pi}{2\zeta} \cdot f_n \cdot S_{gg}}$$

6.3 随机振动疲劳（Miner法则）

随机振动下疲劳损伤用Miner累积损伤准则评估：$D = \sum n_i / N_i$，其中 $n_i$ 是应力幅 $\sigma_i$ 的循环次数，$N_i$ 是材料在该应力幅下的疲劳寿命。$D = 1$ 时判定失效（工程中常取 $D = 0.5$～0.7 作为安全限）。

7. 冲击响应谱（SRS）

🧑‍🎓

军标里有SRS曲线，这和PSD有什么区别？什么时候用SRS？

🎓

PSD描述持续随机振动（飞机颠簸、火箭发射时的持续振动），SRS描述单次冲击事件（炸弹爆炸冲击、跌落、导弹级间分离冲击）。SRS给出在各固有频率的SDOF系统上产生的最大响应加速度——相当于"把冲击对各种刚度结构的破坏潜力"可视化。军标（MIL-STD-810）和航天标准（GEVS）都用SRS定义冲击环境。

7.1 SRS的定义

冲击响应谱 $SRS(f_n)$ 定义为：对固有频率 $f_n$、阻尼比 $\zeta$ 的SDOF系统施加冲击激励 $a(t)$，系统产生的最大绝对响应加速度：

$$SRS(f_n) = \max_t \left|a_{abs}(t; f_n, \zeta)\right|$$

通常阻尼比取 $\zeta = 0.05$（Q=10）。通过扫描 $f_n$（如1 Hz ～ 10000 Hz），即得到SRS曲线。

7.2 半正弦冲击脉冲

最简单的冲击规范，冲击加速度为：

$$a(t) = A_0 \sin\left(\frac{\pi t}{T}\right), \quad 0 \leq t \leq T$$

$A_0$ 是峰值加速度（g），$T$ 是脉冲持续时间。常见规范：11 ms，10 g（电子产品耐久性测试）；6 ms，40 g（汽车电控模块）；3 ms，500 g（弹药引信级别冲击）。

8. 工程应用与实践

🧑‍🎓

振动仿真和实验结果经常对不上，主要是哪些地方容易出问题？

🎓

最常见的三个原因：① 边界条件不符——仿真里是固定支撑，实验里实际支撑有刚度（比如柔性夹具），导致频率偏低；② 阻尼估计不准——FEM里的阻尼通常是猜的，实验才能准确测量；③ 连接件简化过度——螺栓连接、焊缝的刚度在仿真里往往做成刚性连接，但实际有柔性。

分析类型软件功能输出量工程应用
模态分析（Modal）特征值求解固有频率、振型、有效质量共振点识别、结构改进
谐响应分析（Harmonic）频域稳态求解FRF、共振峰、相位旋转机械、发动机激振
随机振动（Random）PSD→响应PSDRMS应力/位移、3σ值航空航天、车辆运输
冲击分析（Transient）时间积分（隐式/显式）时历响应、最大应力跌落、爆炸冲击、碰撞
响应谱（Spectrum）SRS叠加各模态最大响应地震分析、冲击认证

振动分析常见错误

  • 提取模态阶数不够：没有覆盖关注频率范围的1.5～2倍，导致响应谱分析误差大
  • 阻尼模型不当：瑞利阻尼在频率范围两端可能过阻尼，导致高/低频响应失真
  • 质量矩阵选择错误：显式分析需要集中质量矩阵（对角），用了一致质量矩阵会大幅增加计算时间
  • 随机振动中忘记静态预载：有预载（如螺栓预紧力）的结构，刚度因接触状态不同而改变，需要先做非线性静力分析
  • 时间步长太大（瞬态）：时间步长应为最高关注频率周期的1/20以上
跨主题标签（Cross-topics）
结构力学 材料力学（疲劳） 数值方法 结构分析 SDOF响应计算器 随机振动RMS计算器

相关文章