4、多自由度系统：多自由度系统运动方程、模态叠加法、固有频率与振型
好，咱们进入多自由度系统。说实话，单自由度系统你搞明白了，多自由度系统就是它的「升级版」。但这一升级，很多新东西就冒出来了——比如模态、振型这些概念。我在做汽车底盘振动分析时，第一次遇到多自由度系统，当时真有点懵。后来发现，只要抓住「矩阵」和「模态」这两把钥匙，一切就清晰了。
4.1 多自由度系统运动方程
先看运动方程怎么写。对于一个有n个自由度的系统，它的运动方程不再是单个方程，而是一个方程组。写成矩阵形式，就是：
[M]{ẍ} + [C]{ẋ} + [K]{x} = {F(t)}这里：
- [M] —— 质量矩阵，n×n 对称正定矩阵
- [C] —— 阻尼矩阵，n×n 对称矩阵（通常假设为比例阻尼）
- [K] —— 刚度矩阵，n×n 对称正定或半正定矩阵
- {x} —— 位移向量，n×1
- {F(t)} —— 外力向量，n×1
举个例子，一个两自由度弹簧-质量系统：
质量 m1 和 m2，弹簧 k1、k2、k3
运动方程：
m1·ẍ1 + (k1+k2)·x1 - k2·x2 = F1(t)
m2·ẍ2 - k2·x1 + (k2+k3)·x2 = F2(t)
写成矩阵：
[ m1 0 ]{ẍ1} [ k1+k2 -k2 ]{x1} {F1}
[ 0 m2 ]{ẍ2} + [ -k2 k2+k3 ]{x2} = {F2}4.2 无阻尼自由振动与特征值问题
先不考虑阻尼和外力，看看系统本身的特性。方程简化为：
[M]{ẍ} + [K]{x} = {0}假设解的形式为 {x} = {φ}·sin(ωt)，代入后得到：
([K] - ω²[M]){φ} = {0}这就是特征值问题。要使{φ}有非零解，系数矩阵的行列式必须为零：
det([K] - ω²[M]) = 0解这个方程，得到n个特征值 ω₁², ω₂², ..., ωₙ²，开方后就是系统的固有频率。每个特征值对应一个特征向量 {φᵢ}，就是振型。
4.3 模态叠加法
好，现在讲重点——模态叠加法。说白了，就是把复杂的多自由度振动，分解成若干个单自由度振动的叠加。
为什么能这么做？因为振型有一个重要性质：正交性。即：
{φᵢ}ᵀ[M]{φⱼ} = 0 (i ≠ j)
{φᵢ}ᵀ[K]{φⱼ} = 0 (i ≠ j)利用这个性质，我们可以做坐标变换。设：
{x} = [Φ]{q}其中 [Φ] 是振型矩阵，{q} 是模态坐标。代入运动方程，左乘 [Φ]ᵀ，利用正交性，方程就解耦了：
Mᵢ·q̈ᵢ + Cᵢ·q̇ᵢ + Kᵢ·qᵢ = Fᵢ(t) (i=1,2,...,n)每个方程都是独立的单自由度方程！这就是模态叠加法的精髓。
4.4 固有频率与振型的物理意义
固有频率告诉你：系统在哪些频率下容易「激动」。振型告诉你：激动起来是什么样子。
举个例子，一个两自由度系统：
- 第一阶振型：两个质量同向运动，频率最低
- 第二阶振型：两个质量反向运动，频率较高
我在做发动机悬置设计时，就利用这个原理——让发动机的激励频率避开悬置系统的固有频率，同时让振型与车身解耦。嗯，这里要注意，振型节点位置很关键，节点处振幅为零，是布置传感器的好地方。
4.5 模态分析步骤总结
实际工程中，做模态分析一般按这几步走：
- 建立有限元模型 —— 划分网格，定义材料属性
- 组装质量矩阵和刚度矩阵 —— 软件自动完成
- 求解特征值问题 —— 得到固有频率和振型
- 检查结果合理性 —— 看振型是否连续，频率是否合理
- 提取关键模态 —— 通常只关心前几阶低频模态
4.6 知识体系图
下面这张图，把多自由度系统的核心逻辑串起来了：
这张图你看懂了吗？从运动方程出发，分两条路走：一条走无阻尼自由振动，得到特征值问题，算出固有频率和振型；另一条走有阻尼受迫振动，用模态叠加法把方程解耦。最后两条路汇合，就得到了系统的完整动态特性。
4.7 工程应用小贴士
| 应用场景 | 关注点 | 我的建议 |
|---|---|---|
| 汽车车身模态 | 第一阶弯曲/扭转频率 | 避开发动机怠速频率（约25-30 Hz） |
| 桥梁风振 | 低阶模态阻尼比 | 注意涡激振动，必要时加TMD |
| 精密仪器隔振 | 高阶模态影响 | 模态截断时保留质量参与系数大的模态 |
| 转子动力学 | 临界转速 | 考虑陀螺效应，Campbell图是必备工具 |
公众号:蓝海资料掘金营，微信deep3321