简谐运动（SHM）— 固有频率·阻尼比到FEM特征值分析

分类: 物理基础 | 2026-03-25 | 网站地图

1. 简谐运动是什么 — 最重要的振动模型

🧑‍🎓

塔科玛海峡大桥倒塌是什么现象？我听说与振动有关系。

🎓

1940年美国发生的桥梁倒塌事故。桥梁的固有振动频率与风产生的涡激振动（Kármán涡）周期恰好一致，导致振动急剧增大。准确来说是颤振不稳定这样的气弹性现象，但共振是根本原因。理解简谐运动的共振机制，就能看出"为什么很小的风也能导致巨大的桥梁崩塌"。设计者没有考虑振动的悲剧性案例。

简谐运动（Simple Harmonic Motion, SHM）是指"复原力与位移成正比，指向位移中心"的运动。遵循胡克定律 $F=-kx$ 的弹簧-质量系统是最基本的例子。

简谐运动之所以重要：任何线性振动系统都可以表示为多个简谐运动（固有模式）的叠加（模式叠加法）。换句话说，理解简谐运动是理解所有复杂振动现象的基础。

2. 简谐运动方程和固有频率

将胡克定律 $F = -kx$ 代入牛顿第二定律 $F = ma = m\ddot{x}$：

$$m\ddot{x} + kx = 0$$

这个二阶线性常微分方程的通解为：

$$x(t) = A\cos(\omega_n t + \phi)$$

固有角频率（Natural Angular Frequency）：

$$\omega_n = \sqrt{\frac{k}{m}} \quad \text{(rad/s)}$$

固有频率（Natural Frequency）和固有周期（Natural Period）：

$$f_n = \frac{\omega_n}{2\pi} = \frac{1}{2\pi}\sqrt{\frac{k}{m}} \quad \text{(Hz)}, \qquad T = \frac{1}{f_n} = 2\pi\sqrt{\frac{m}{k}} \quad \text{(s)}$$

🧑‍🎓

想要提高固有振动频率时应该怎么做？

🎓

由$\omega_n=\sqrt{k/m}$可知，需要提高刚性$k$或降低质量$m$。对汽车来说，既要轻量化车体（降低$m$）又要提高车身刚性（提高$k$），这是"跑、弯、停"高性能化的关键。但还要考虑乘坐舒适性的平衡。CAE中通过模态分析找出问题振动模式，然后反复进行刚性提升与轻量化同步追求的设计迭代。

3. 单摆和近似

长度为 $L$ 的单摆（摆角较小时）：

$$\omega_n = \sqrt{\frac{g}{L}}, \qquad T = 2\pi\sqrt{\frac{L}{g}}$$

当 $L=1\,\text{m}$ 时：$T = 2\pi\sqrt{1/9.81} = 2.006\,\text{s}$。这就是钟摆（1来回2秒 → 计时1秒）的原理。

这里使用了小角近似（$\theta \ll 1\,\text{rad}$），摆角较大时 $T$ 会增长（非线性效应）。在FEM中，当梁的挠度很大时会出现几何非线性，固有频率会改变（预应力效应、旋转加速度效应）现象与此原理相同。

4. 阻尼振动（1自由度阻尼系统）

现实的振动系统总有阻尼存在（内部摩擦·空气阻力·结构接合部摩擦）：

$$m\ddot{x} + c\dot{x} + kx = 0$$

阻尼比（Damping Ratio）：

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

三个阻尼区域

欠阻尼（Underdamped, $\zeta < 1$）：振动同时衰减。最常见的情况。

$$x(t) = Ae^{-\zeta\omega_n t}\cos(\omega_d t + \phi), \quad \omega_d = \omega_n\sqrt{1-\zeta^2}$$

临界阻尼（Critically Damped, $\zeta = 1$）：不振动，最快趋向稳定。门自关闭器的设计目标。

过阻尼（Overdamped, $\zeta > 1$）：不振动，缓慢趋向稳定。粘度高的情况。

结构物类型	代表性阻尼比 $\zeta$	备注
焊接钢结构	0.01～0.02	阻尼微小，易发生共振
螺栓连接钢结构	0.02～0.04	接合部摩擦增大
RC（钢筋混凝土）	0.03～0.08	裂缝增加时增大
木造建筑	0.03～0.08	接合部多，相对较高
隔震装置（层积橡胶）	0.10～0.20	有意设计高阻尼
轮胎（橡胶）	0.10～0.15	内部损失高

5. 强迫振动和共振

外部施加周期力 $F_0\cos\Omega t$ 的情况：

$$m\ddot{x} + c\dot{x} + kx = F_0\cos\Omega t$$

稳态状态下的振幅放大因子（动态因子）：

$$M(\beta) = \frac{X}{X_{st}} = \frac{1}{\sqrt{(1-\beta^2)^2 + (2\zeta\beta)^2}}, \quad \beta = \frac{\Omega}{\omega_n}$$

其中 $X_{st} = F_0/k$ 是静位移。

共振时（$\beta = 1$, $\Omega = \omega_n$）的振幅放大因子：

$$M\big|_{\beta=1} = \frac{1}{2\zeta}$$

钢结构的 $\zeta=0.01$ 时 $M=50$ 倍！这就是塔科玛海峡大桥等共振倒塌的原因。

共振回避设计指南

保持励振频率（发动机转速·外部振动源）与固有频率之比 $\beta < 0.7$ 或 $\beta > 1.3$
增加动力吸振器（消振器）
应用制振材料（粘弹性材料）

6. 多自由度系统（MDOF）的拓展

$n$ 个质量-弹簧系统串联的多自由度系统：

$$[M]\{\ddot{u}\} + [C]\{\dot{u}\} + [K]\{u\} = \{F(t)\}$$

自由振动、无阻尼情况下的特征值问题：

$$\left([K] - \omega_i^2[M]\right)\{\phi_i\} = \{0\}$$

求解 $i = 1, 2, \ldots, n$ 得到 $n$ 个固有角频率 $\omega_i$ 和固有模式 $\{\phi_i\}$。

模态叠加法（Modal Superposition）

强迫振动响应可表示为各固有模式独立响应的和：

$$\{u(t)\} = \sum_{i=1}^n q_i(t)\{\phi_i\}$$

每个 $q_i(t)$ 只需求解简谐运动（SDOF）方程，大幅提高计算效率。这是FEM模态分析的核心思想。

7. FEM特征值分析（模态分析）

大规模FEM模型（数百万自由度）的特征值问题通常用迭代法求解。最广泛使用的是Lanczos法。

Lanczos法概述

以随机初始向量开始，反复作用 $[K]^{-1}[M]$ 变换
每步正交化基向量（Gram-Schmidt过程）
得到的小三对角矩阵的特征值近似原大规模问题的特征值
能高效求出低阶固有模式

🧑‍🎓

汽车设计中要做多少次模态分析？应该很耗时吧。

🎓

从设计初期到最终确认要做几百次。车身刚性变化1Hz就会影响乘坐舒适性，引擎盖开闭的"嗙"声（闭合声）也可能导致需要调整固有值。行驶中的"嗡嗡声"（40～100Hz低频噪音）也通过模态分析追根溯源。CAE工程师每天都要多次运行模态分析来确认设计变更的效果。

特征值分析的FEM输入和输出

输入：质量矩阵 $[M]$（由密度和几何形状自动计算）、刚性矩阵 $[K]$（由杨氏模量、泊松比和形状自动计算）、边界条件

输出：固有频率列表（Hz）、固有模式形状（动画显示）、有效质量（各模式的质量参与率）

8. 实验模态分析（EMA）和V&V

FEM模态分析的验证需要实物振动试验数据。

试验步骤

加振：用冲击锤（锤击法）或电磁振动台对结构加振
测量：用压电加速度传感器在多点测量振动加速度
获得FRF：计算力传感器的输入力与加速度响应的传递函数（频率响应函数, FRF）
提取模态：从FRF识别固有频率、模式形状、阻尼比（Polymax、EFDD等算法）

FEM-EMA的相关性确认（MAC）

模式形状的相关性用模态保证准则（Modal Assurance Criterion, MAC）定量化：

$$MAC_{ij} = \frac{|\{\phi_i\}^T\{\psi_j\}|^2}{(\{\phi_i\}^T\{\phi_i\})(\{\psi_j\}^T\{\psi_j\})}$$

目标是 $MAC_{ii} > 0.9$ 表示良好相关性，$MAC_{ij} < 0.1$（$i\neq j$）确保模式间的正交性。

9. 实际案例：建筑物抗震设计

建筑物的第1阶固有周期随高度增长（经验公式）：

$$T_1 \approx 0.1 \times N \text{（楼层数）} \quad \text{或} \quad T_1 \approx H/40 \text{（高度 m）}$$

地震的优势周期（随地基种类变化）与建筑物固有周期一致时发生共振，地震灾害加剧。据说1995年阪神大地震中，固有周期为0.6～0.8 s的建筑物（相当于5～8层建筑）受灾最重。

隔震结构的原理

在建筑基础上安装隔震装置（层积橡胶支座），有意延长建筑固有周期（$T \approx 3～5\,\text{s}$）。使其偏离一般地震动的优势周期（0.5～1.5 s），从而大幅降低地震能量向建筑的输入。这是积极"回避"共振的设计。隔震建筑的设计必须进行包含层积橡胶非线性刚性和阻尼特性的详细FEM动态分析。

跨领域标签（Cross-topics）

#振动分析 #模态分析 #固有频率 #特征值分析 #波动