简谐运动（SHM）— 固有频率・阻尼比与有限元模态分析
目录
1. 共振：最危险的振动现象
塔科马海峡大桥为什么会被风吹塌？那风速也没多快啊，风能有那么大的力吗？
不是风力有多大，而是共振。风绕过桥体产生卡门涡街（交替脱落的涡旋），涡旋脱落频率恰好和桥的扭转固有频率接近——大约0.2 Hz。每一次涡旋脱落都给桥一个小推力，因为频率吻合，这些小推力的效果会累积叠加，振幅越来越大，最终导致破坏。不是力大，是时机准。
简谐运动（Simple Harmonic Motion, SHM）是最基本的振动形式，也是工程振动分析的基础。理解SHM和共振现象，对于避免结构失效、设计隔振系统至关重要。
2. 简谐运动的方程
无阻尼弹簧-质量系统的运动方程：
解为正弦（余弦）形式：
其中：
- $A$：振幅（m），由初始条件决定
- $\omega_n$：固有圆频率（rad/s）
- $\phi$：初相位（rad）
速度和加速度：
加速度与位移成正比、方向相反——这是SHM的本质特征。
3. 固有频率与圆频率
将 $x = A\cos(\omega_n t)$ 代入运动方程，得到固有圆频率：
固有频率（Hz）和固有周期（s）：
| 结构/系统 | 典型固有频率 | 说明 |
|---|---|---|
| 高层建筑（第1阶） | 0.1 - 0.5 Hz | 层数越多，频率越低 |
| 大跨桥梁（第1阶弯曲） | 0.1 - 1 Hz | 跨径越大，频率越低 |
| 汽车车身（第1阶） | 20 - 50 Hz | 整车弯扭模态 |
| 发动机悬置 | 8 - 15 Hz | 需避开怠速激励（约25 Hz） |
| 旋转机械叶片 | 100 - 5000 Hz | 需避开转频及其倍频 |
| 人体站立自然频率 | 约 10 Hz | 全身振动标准的参考 |
4. 阻尼：振动的消耗机制
有阻尼弹簧-质量系统：
阻尼比（Damping Ratio）：
其中 $c_c = 2\sqrt{mk} = 2m\omega_n$ 是临界阻尼系数。
根据阻尼比大小，系统的响应分为三类：
| 阻尼状态 | 条件 | 响应特征 | 工程实例 |
|---|---|---|---|
| 欠阻尼 | $\xi < 1$ | 振荡衰减 | 大多数结构（金属：0.01~0.05） |
| 临界阻尼 | $\xi = 1$ | 最快无振荡回零 | 理想精密仪器设计目标 |
| 过阻尼 | $\xi > 1$ | 缓慢无振荡回零 | 缓冲器（门铰链） |
欠阻尼系统的阻尼固有频率：
常见材料/结构的阻尼比（$\xi$）：
| 材料/结构 | 典型阻尼比 ξ |
|---|---|
| 钢结构（焊接） | 0.01 ~ 0.02 |
| 混凝土结构 | 0.03 ~ 0.07 |
| 木结构 | 0.05 ~ 0.10 |
| 橡胶减振器 | 0.05 ~ 0.20 |
| 流体阻尼器（液压） | 0.20 ~ 0.40 |
5. 共振：激励频率与固有频率的相遇
有阻尼受迫振动方程：
稳态响应的幅值放大系数（Dynamic Amplification Factor, DAF）：
其中 $r = \omega/\omega_n$ 是频率比。
在共振点（$r = 1$）：
对于阻尼比 $\xi = 0.02$（典型钢结构）：$DAF_{\max} = 1/(2 \times 0.02) = 25$！同样的静力可以产生25倍的动力响应——这就是共振的危险所在。
那实际工程里如何避免共振？频率完全错开的话，设计起来不是很麻烦吗？
一般要求固有频率与激励频率保持1.3倍以上的分离比——要么低很多（隔振区），要么高很多（刚性区）。在激励频率附近的±30%范围内是"危险区域"，要通过改变结构质量、刚度或增加阻尼来避开。旋转机械里这叫"Campbell图分析"，在所有工作转速范围内检查固有频率和激励倍频有没有交叉点。
6. FEM模态分析：特征值问题
无阻尼自由振动方程（设解为 $\{u\} = \{\phi\}e^{i\omega t}$）：
这是广义特征值问题。非零解的条件：
解出 $n$ 个特征值 $\omega_1^2 \leq \omega_2^2 \leq \cdots \leq \omega_n^2$（固有角频率的平方），对应 $n$ 个特征向量 $\{\phi_i\}$（振型向量）。
模态分析在Nastran/Abaqus中的典型步骤：
- 建立有限元模型（网格、材料、边界条件）
- 指定提取的模态数（通常20~100阶）
- 选择求解器（Lanczos法最常用）
- 查看固有频率列表和振型动画
- 验证：对比试验数据（模态测试MAC值 > 0.8）
7. 模态振型与模态质量
振型（Mode Shape）描述结构在某一固有频率下的相对变形形状。振型之间满足正交性：
模态质量和模态刚度（质量归一化下）：
有效质量（Effective Mass）是判断某阶模态重要性的指标：
所有模态有效质量之和等于系统总质量。通常要求前若干阶模态的累积有效质量超过总质量的90%，以确保模态分析的完整性。
8. 历史案例：塔科马海峡大桥（1940年）
1940年11月7日，美国华盛顿州塔科马海峡大桥在42 mph（约19 m/s）的风速下发生颤振，最终坍塌。
| 参数 | 数值 |
|---|---|
| 主跨长度 | 853 m |
| 桥体扭转固有频率 | 约 0.2 Hz |
| 卡门涡街脱落频率（灾害风速下） | 约 0.2 Hz |
| 颤振风速 | 约 19 m/s（42 mph） |
| 最大扭转角（坍塌前） | 约 ±45° |
坍塌根本原因：桥梁设计采用实腹板桁架（Solid girder），风对桥面产生周期性升力，与桥梁扭转频率共振。现代悬索桥采用流线型箱形断面（Aerodynamic box girder），可有效避免卡门涡街共振。
此事故推动了气动弹性理论（Aeroelasticity）在桥梁工程中的系统应用，今天所有大跨桥梁设计都要做风洞试验和气动弹性数值仿真。
9. 工程应用：隔振与调谐质量阻尼器
隔振设计（Vibration Isolation）：
隔振系统的力传递率（Transmissibility）：
当 $r > \sqrt{2}$ 时，$TR < 1$，即隔振有效。要使 $TR < 0.1$，需要 $r > 3.2$——也就是说，隔振系统固有频率需低于激励频率的 1/3.2 以下。
调谐质量阻尼器（Tuned Mass Damper, TMD）：
在主体结构上附加小质量-弹簧-阻尼子系统，调谐到主结构固有频率，将振动能量转移到TMD中耗散。
典型工程案例：
- 台北101大楼：钢球型TMD，质量660吨，可使顶层加速度减小40%
- 上海中心大厦：摆型TMD，质量1000吨，减振25-30%
- 东京晴空塔：油压阻尼器系统，减少地震和风载振动
总结
简谐运动是一切振动分析的基础。从 $\omega_n = \sqrt{k/m}$ 到FEM模态特征值问题，从共振放大系数 $1/(2\xi)$ 到TMD减振设计——理解SHM就掌握了结构动力学的核心物理图像。所有复杂振动响应都可以分解为各阶模态SHM的线性叠加（模态叠加法）。