振动·动解析的基础 — 从固有振动到随机振动·冲击解析

分类：基础理论 | 2026-03-25 | 网站地图

CAE visualization for vibration dynamics - technical simulation diagram — Vibration Dynamics

振动的基本概念

振动解析是「预测结构物对动荷载（随时间变化的力）如何应答」的解析。与静力分析不同，需要加入惯性力、阻尼力以及时间的概念。汽车、飞机、桥梁、家电产品 — 任何东西都与振动问题息息相关。

🧑‍🎓

老师，振动解析是为了了解什么呢？和静力分析有什么区别呢...

🎓

最重要的区别是「共振（resonance）」的有无。静荷重100N和具有相同100N振幅的动荷重 — 如果这与结构物的固有振动数一致会怎样？应答在理论上会变为无穷大。也就是说，即使100N也会导致建筑物倒塌。塔科马悬索桥的崩塌（1940年）是风引起的自激振动所致。因此在振动解析中，「确认固有振动数在哪里、外部扰动的频率是否重合」是最重要的。

1.1 振动的分类

分类轴	种类	说明
外力有无	自由振动	仅给予初始变位/速度，无外力。以固有振动数振动
外力有无	强制振动	持续外力下振动。稳态应答和过渡应答重叠
阻尼有无	无阻尼振动	无能量耗散。永久振动（理想）
	欠阻尼	$\zeta < 1$：振幅逐渐衰减的振动。大多数结构
	过阻尼	$\zeta > 1$：不振动，指数衰减。防振隔振器
自由度	SDOF	单自由度系。理论分析的基础
自由度	MDOF	多自由度系。FEM离散化后可达数百万自由度

1.2 固有圆频率和固有频率

质量 $m$、弹簧常数 $k$ 的单自由度系的固有圆频率：

$$\omega_n = \sqrt{\frac{k}{m}} \quad [\text{rad/s}]$$

固有频率（Hz）和固有周期（s）：

$$f_n = \frac{\omega_n}{2\pi} \quad [\text{Hz}], \qquad T_n = \frac{1}{f_n} = \frac{2\pi}{\omega_n} \quad [\text{s}]$$

单自由度（SDOF）系的详细分析

2.1 运动方程和解的形态

SDOF系的运动方程：

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

引入阻尼比 $\zeta = c / (2\sqrt{km}) = c / c_c$（其中 $c_c = 2\sqrt{km}$ 是临界阻尼系数）：

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

自由振动（$f=0$）的解：

条件	ζ值	解的形式	行为
欠阻尼	0 < ζ < 1	$x(t) = e^{-\zeta\omega_n t}(A\cos\omega_d t + B\sin\omega_d t)$	衰减振动 ($\omega_d = \omega_n\sqrt{1-\zeta^2}$)
临界阻尼	ζ = 1	$x(t) = (A + Bt)e^{-\omega_n t}$	最快回到原点（无振动）
过阻尼	ζ > 1	$x(t) = Ae^{s_1 t} + Be^{s_2 t}$ ($s_{1,2}$ < 0)	缓慢回到原点（无振动）

🧑‍🎓

实际结构物的阻尼比 ζ 大约是多少呢？能测量吗？

🎓

焊接钢结构物 ζ ≈ 0.01〜0.03（1〜3%）、混凝土结构 ζ ≈ 0.05（5%）是经验值。测量通常在实务中采用「半值宽度法」或「对数衰减率法」。用冲击锤敲击，加速度传感器计测应答，从FRF进行同定的实验模态分析是标准手法。橡胶隔振器系统 ζ ≈ 0.1〜0.3 会更高。

2.2 正弦波加振的稳态应答和共振

$f(t) = F_0 \sin\omega t$ 的稳态应答：

$$x(t) = \frac{F_0/k}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}} \sin(\omega t - \phi), \quad r = \frac{\omega}{\omega_n}$$

位相延迟：

$$\phi = \arctan\left(\frac{2\zeta r}{1 - r^2}\right)$$

动态放大系数（DAF, Dynamic Amplification Factor）：

$$\text{DAF}(r, \zeta) = \frac{1}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$

$r=1$（共振）时的应答：

$$\text{DAF}(1, \zeta) = \frac{1}{2\zeta}$$

当阻尼比 $\zeta = 0.01$ 的钢结构发生共振时，静荷载的50倍应答会产生。

🧑‍🎓

那样的话，把固有振动数从外力频率离开就安全了对吧？需要离开多少呢？

🎓

通用的设计目标是在固有振动数和加振频率之间设置±20〜30%的「回避带域」。但是像汽车发动机那样转速在很大范围内变化的情况，全周波数域回避是困难的。因此有时会反过来让共振点通过，同时用充分的阻尼来限制应答。这就是动态吸振器（动态减震器）的思想，即附加副系统，把主系的共振分裂成两个，来降低应答的手法。

2.3 过渡应答 — 阶跃输入和脉冲响应

阶跃荷载 $f(t) = F_0 \cdot H(t)$（$H$：Heaviside函数）的应答（欠阻尼）：

$$x(t) = \frac{F_0}{k}\left[1 - e^{-\zeta\omega_n t}\left(\cos\omega_d t + \frac{\zeta}{\sqrt{1-\zeta^2}}\sin\omega_d t\right)\right]$$

脉冲响应函数（格林函数）$h(t)$：

$$h(t) = \frac{1}{m\omega_d} e^{-\zeta\omega_n t} \sin\omega_d t, \quad t \geq 0$$

任意荷载的应答可通过Duhamel积分（卷积积分）求得：

$$x(t) = \int_0^t h(t-\tau) f(\tau) \, d\tau$$

多自由度（MDOF）系和固有值解析

3.1 运动方程的矩阵形式

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

FEM离散化的结构物具有 $n$ 个自由度（节点数×方向数）的矩阵。实用结构分析中 $n = 10^4 \sim 10^7$。

3.2 固有值问题

假设无阻尼自由振动（$\mathbf{C}=\mathbf{0}$, $\mathbf{f}=\mathbf{0}$），设 $\mathbf{u} = \boldsymbol{\phi} e^{i\omega t}$：

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

这是固有值问题。非平凡解的条件：

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

第 $i$ 个固有圆频率 $\omega_i$ 和特征向量（模态形状）$\boldsymbol{\phi}_i$ 的对：

$\omega_i$ 按从小到大排列，分别称为一阶模态、二阶模态等
$\boldsymbol{\phi}_i$ 是模态形状，表示振幅的相对分布

3.3 模态的正交性和质量正规化

模态关于质量矩阵和刚度矩阵正交：

$$\boldsymbol{\phi}_i^T \mathbf{M} \boldsymbol{\phi}_j = \delta_{ij}, \qquad \boldsymbol{\phi}_i^T \mathbf{K} \boldsymbol{\phi}_j = \omega_i^2 \delta_{ij}$$

（质量正规化情况下）

🧑‍🎓

在Ansys中执行「模态解析」时，显示「1阶固有振动数：23.5 Hz」，这是什么意思呢？

🎓

该结构物能够以最缓慢速度（能量最低状态）振动的振动数就是 23.5 Hz。换句话说，用发动机转速来讲是 23.5 × 60 = 1410 rpm 附近会发生共振。同时显示的模态形状（变形动画）可以让你看到「是如何变形的」。例如，底盘的弯曲模式作为一阶出现，那么就要确认该方向是否有加振源。

3.4 模态叠加法（模态叠加）

引入模态坐标 $\mathbf{q}$ 使 $\mathbf{u} = \boldsymbol{\Phi}\mathbf{q}$（$\boldsymbol{\Phi}$ 为模态矩阵），则耦联 $n$ 自由度系统转换为 $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$$

实务中，仅用下位 $m$ 个模态（$m \ll n$）就能精确近似，因此可大幅降低计算成本（模态截断）。

截断模态数的目标：包含加振周波数3倍以上的固有振动数。

实用的数值积分法

4.1 Newmark-β 法

振动问题时刻历响应计算中最常用的无条件稳定隐式方法。速度和变位的更新式：

$$\dot{\mathbf{u}}_{n+1} = \dot{\mathbf{u}}_n + [(1-\gamma)\ddot{\mathbf{u}}_n + \gamma\ddot{\mathbf{u}}_{n+1}]\Delta t$$

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

加速度的确定通过在 $t_{n+1}$ 处求解运动方程得到：

$$\mathbf{K}_\text{eff} \mathbf{u}_{n+1} = \mathbf{f}_\text{eff}^{n+1}$$

$\mathbf{K}_\text{eff} = \mathbf{K} + \frac{\gamma}{\beta\Delta t}\mathbf{C} + \frac{1}{\beta\Delta t^2}\mathbf{M}$ 是有效刚度矩阵。

4.2 显式法（中心差分法）的对比

项目	显式法（中心差分）	隐式法（Newmark-β）
稳定条件	$\Delta t \leq \Delta t_{cr} = \frac{2}{\omega_{max}}$	无（β≥1/4, γ≥1/2 时）
时间步长	最小单元尺寸/音速（非常小）	激励周波数的 1/10〜1/20 左右
单步成本	非常便宜（无需矩阵分解）	昂贵（每步求解联立方程）
非线性·接触	易处理	需要反复迭代
适用问题	冲击·爆炸（短时·高频）	地震·振动（长时·低频）
代表求解器	Abaqus/Explicit, LS-DYNA	Abaqus/Standard, Ansys Mechanical

4.3 HHT-α 法（高精度隐式法）

Hilber-Hughes-Taylor的 $\alpha$ 法是Newmark-β法的扩展，在保持低频精度的同时，对高周波数成分进行数值衰减：

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

$-1/3 \leq \alpha \leq 0$：当 $\alpha = 0$ 时退化为Newmark平均加速度法。$\alpha = -0.1$ 左右是实务标准设置。

频率响应分析（FRA）

5.1 频率响应函数（FRF）

定常正弦波加振下的应答的传递函数。复数表示：

$$H(\omega) = \frac{X(\omega)}{F(\omega)} = \frac{1}{k - m\omega^2 + ic\omega} = \frac{1/k}{(1 - r^2) + i(2\zeta r)}$$

FRF的形式：

导纳（导纳FRF）： $H = X/F$ [m/N]
流动度（流动度）： $Y = \dot{X}/F$ [m/(N·s)]
加速度度（加速度度）： $A = \ddot{X}/F$ [m/(N·s²)]

🧑‍🎓

FRF是在实验中测量的吗？和FEM分析有什么不同呢？

🎓

FRF既可以通过实验（脉冲锤试验或激振器加振）进行测量，也可以通过FEM进行计算。实务中两者的对比是很重要的。如果FEM与实验FRF大幅不符 — 比如固有振动数偏差10%以上 — 那么网格质量、质量分布或边界条件某处有问题的信号。这种「实验模态与分析模态的相关确认（MAC值）」的验证工作，在汽车·飞机设计中是必需的流程。

5.2 阻尼模型的种类

模型	数学表达式	特点	应用
粘性阻尼	$c\dot{x}$	与速度成比例。时间域中自然	流体阻力·阻尼器
结构阻尼（滞后）	$ig k x$（复数刚度）	与振幅成比例。频率域易处理	金属内部摩擦
Rayleigh阻尼	$\mathbf{C} = \alpha\mathbf{M} + \beta\mathbf{K}$	保持模态正交性	FEM通用
库仑摩擦阻尼	$\mu N \text{sign}(\dot{x})$	非线性。与振幅无关的常数	接头·滑动部位

5.3 MAC值（模态保证准则）

分析模态 $\boldsymbol{\phi}_i^\text{FEM}$ 与实验模态 $\boldsymbol{\phi}_j^\text{EXP}$ 的相关系数：

$$\text{MAC}_{ij} = \frac{|\boldsymbol{\phi}_i^{\text{FEM},T}\boldsymbol{\phi}_j^\text{EXP}|^2}{(\boldsymbol{\phi}_i^{\text{FEM},T}\boldsymbol{\phi}_i^\text{FEM})(\boldsymbol{\phi}_j^{\text{EXP},T}\boldsymbol{\phi}_j^\text{EXP})}$$

$\text{MAC} \geq 0.9$：协相好。$\text{MAC} < 0.8$：需要改进模型。

随机振动

火箭发射振动、行驶路面凹凸、飞机大气乱流等，难以用时间域描述的「概率性振动输入」对结构应答的预测是随机振动分析。

🧑‍🎓

随机振动怎样评价呢？时刻历随意变化，好像没法分析...

🎓

时间域虽然确实随意变化，但转换到周波数域就能看出统计特性。那就是功率谱密度（PSD）。PSD表示各个周波数有多少能量。给定输入PSD后，用系统的传递函数（FRF）计算输出PSD，进一步求RMS值（均方根值）就是基本流程。

6.1 功率谱密度（PSD）

时间信号 $x(t)$ 的PSD：

$$S_x(f) = \lim_{T\to\infty} \frac{1}{T} |X_T(f)|^2$$

$X_T(f)$ 是 $[-T/2, T/2]$ 区间的傅里叶变换。单位为 [g²/Hz] 或 [m²/s³] 等。

RMS值（均方根值）：

$$x_\text{rms} = \sqrt{\int_0^\infty S_x(f) \, df}$$

6.2 输入-输出关系（Wiener-Khinchin定理）

线性系统的输入PSD $S_F$ 与输出PSD $S_x$ 的关系：

$$S_x(f) = |H(f)|^2 S_F(f)$$

SDOF系白色噪声输入 $S_0 = \text{const}$ 的应答RMS：

$$x_\text{rms} = \sqrt{\frac{\pi S_0 \omega_n}{2k^2}\left(\frac{1}{2\zeta}\right)} = \sqrt{\frac{\pi f_n S_0}{2k^2\zeta}}$$

6.3 3σ设计基准

假设高斯分布，峰值的超越概率：

峰值	超越概率	应用基准
$1\sigma = x_\text{rms}$	31.7%	通常设计不足够
$2\sigma$	4.55%	一般设计
$3\sigma$	0.27%	航空航天·MIL规格的标准
$4\sigma$	0.0063%	高可靠性要求部件

随机振动分析结果通常与材料的屈服应力进行比较：$3\sigma_\text{stress} \leq \sigma_Y / \text{SF}$（SF：安全系数）。

冲击应答谱（SRS）

7.1 什么是SRS

对于冲击输入（落下冲击·爆炸等），「具有各种固有振动数的SDOF系分别在最大处的应答」用频率函数表示的就是冲击应答谱（Shock Response Spectrum, SRS）。

$$\text{SRS}(\omega_n) = \max_t |x(t)|, \quad x(t) = \frac{1}{\omega_d}\int_0^t \ddot{u}_g(\tau) e^{-\zeta\omega_n(t-\tau)}\sin\omega_d(t-\tau)\,d\tau$$

7.2 代表的冲击脉冲形状

脉冲形状	定义	SRS 的特点
半正弦（Half-sine）	$a(t) = A\sin(\pi t/D)$, $0 \leq t \leq D$	低周波域近似为 $A$
矩形（Rectangular）	$a(t) = A$, $0 \leq t \leq D$	最大值 = $2A$（对称）
锯齿波（Sawtooth）	线性增加或减少	高周波域容易出现大值

🧑‍🎓

SRS什么时候使用呢？只有火箭等航空航天产品吗？

🎓

航空航天产品最容易想象，但物流·运输行业也经常使用。商品运输中的落下冲击（30cm落下 ≈ 100G冲击）对电子基板的耐受性就用SRS评价。IEC 60068-2-27（冲击试验）和MIL-STD-810中规定了SRS的试验条件。军用设备·计测仪器·医疗设备的设计中是必需的评价手法。

实务应用例

8.1 防振隔振器设计

发动机隔振器的设计目标是「使发动机振动不传递到车体」。传递率（TR: transmissibility）：

$$\text{TR} = \frac{|f_\text{transmitted}|}{|f_\text{input}|} = \sqrt{\frac{1 + (2\zeta r)^2}{(1-r^2)^2 + (2\zeta r)^2}}$$

当 $r > \sqrt{2}$（加振频率是固有振动数的$\sqrt{2}$倍以上）时，TR < 1，防振效果出现。汽车发动机（怠速：700rpm = 11.7Hz）的隔振器固有振动数设定为5〜8Hz，可确保怠速时 $r \approx 1.5$ 以上。

8.2 地震应答解析（免震结构）

建筑地震应答分析中地动加速度 $\ddot{u}_g(t)$ 与相对运动的关系：

$$m\ddot{x}_\text{rel} + c\dot{x}_\text{rel} + kx_\text{rel} = -m\ddot{u}_g(t)$$

免震结构通过层积胶垫支承把固有周期延长到3〜5秒（$f_n = 0.2\sim0.3$ Hz），远离地震波主要能量集中的1〜2Hz带，从而降低应答。

8.3 NVH分析（噪音·振动·乐音感）

汽车NVH分析是最大规模的振动分析之一：

车体刚性（BIW固有振动数）： 通常30〜50Hz — 发动机2阶激励（怠速700rpm → 23Hz）的干涉需要管理
面板的声辐射效率： 仪表板·地板的振动 → 车室内音压的结构-声学耦联分析
子架的传递路径分析（TPA）： 路面 → 轮胎 → 悬架 → 车体的振动传递路径的确认

振动·动解析的基础 — 从固有振动到随机振动·冲击解析

振动的基本概念

1.1 振动的分类

1.2 固有圆频率和固有频率

单自由度（SDOF）系的详细分析

2.1 运动方程和解的形态

2.2 正弦波加振的稳态应答和共振

2.3 过渡应答 — 阶跃输入和脉冲响应

多自由度（MDOF）系和固有值解析

3.1 运动方程的矩阵形式

3.2 固有值问题

3.3 模态的正交性和质量正规化

3.4 模态叠加法（模态叠加）

实用的数值积分法

4.1 Newmark-β 法

4.2 显式法（中心差分法）的对比

4.3 HHT-α 法（高精度隐式法）

频率响应分析（FRA）

5.1 频率响应函数（FRF）

5.2 阻尼模型的种类

5.3 MAC值（模态保证准则）

随机振动

6.1 功率谱密度（PSD）

6.2 输入-输出关系（Wiener-Khinchin定理）

6.3 3σ设计基准

冲击应答谱（SRS）

7.1 什么是SRS

7.2 代表的冲击脉冲形状

实务应用例

8.1 防振隔振器设计

8.2 地震应答解析（免震结构）

8.3 NVH分析（噪音·振动·乐音感）

相关交互式工具

动手验证理论

相关主题

振动·动解析的基础 — 从固有振动到随机振动·冲击解析

振动的基本概念

1.1 振动的分类

1.2 固有圆频率和固有频率

单自由度（SDOF）系的详细分析

2.1 运动方程和解的形态

2.2 正弦波加振的稳态应答和共振

2.3 过渡应答 — 阶跃输入和脉冲响应

多自由度（MDOF）系和固有值解析

3.1 运动方程的矩阵形式

3.2 固有值问题

3.3 模态的正交性和质量正规化

3.4 模态叠加法（模态叠加）

实用的数值积分法

4.1 Newmark-β 法

4.2 显式法（中心差分法）的对比

4.3 HHT-α 法（高精度隐式法）

频率响应分析（FRA）

5.1 频率响应函数（FRF）

5.2 阻尼模型的种类

5.3 MAC值（模态保证准则）

随机振动

6.1 功率谱密度（PSD）

6.2 输入-输出关系（Wiener-Khinchin定理）

6.3 3σ设计基准

冲击应答谱（SRS）

7.1 什么是SRS

7.2 代表的冲击脉冲形状

实务应用例

8.1 防振隔振器设计

8.2 地震应答解析（免震结构）

8.3 NVH分析（噪音·振动·乐音感）

相关交互式工具

动手验证理论

相关文章

相关主题