第 6 章 · 多自由度与模态分析
把两自由度推广到 $n$ 自由度,方法论的顶点是模态分析:把耦合系统解耦成 $n$ 个独立单振子。这是实验模态测试、有限元动力学、随机振动的共同底座。
特征值问题与模态
无阻尼自由振动 $\mathbf{M}\ddot{\mathbf{x}} + \mathbf{K}\mathbf{x} = 0$,设 $\mathbf{x} = \boldsymbol\phi e^{i\omega t}$:
$n$ 组解:固有频率 $\omega_i$ 与振型 $\boldsymbol\phi_i$。模态的两大性质:关于 M 加权正交 $\boldsymbol\phi_i^T\mathbf{M}\boldsymbol\phi_j = 0\ (i\ne j)$、可按质量归一($\boldsymbol\phi_i^T\mathbf{M}\boldsymbol\phi_i = 1$)。振型的物理意义:系统"偏爱的变形模式",最低阶模态 = 最软的变形方式(结构动力学中它决定地震响应的主力)。
模态叠加:解耦流水线
坐标变换 $\mathbf{x} = \boldsymbol\Phi\mathbf{q}$($\boldsymbol\Phi$ = 模态矩阵),正交性让方程组对角化:
每个模态变成单振子被"模态激励"驱动——解耦完成。响应 = 模态响应叠加。工程惯例:只取低阶 $m$ 个模态($m \ll n$)——模态截断是有限元动力学的省钱法宝(低频激励用低阶模态足够)。
阻尼的模型化与实验模态
比例阻尼($C = \alpha M + \beta K$)才能用无阻尼模态精确解耦;一般阻尼用模态阻尼比近似(每个模态各配一个 $\zeta_i$,由实验测定)。实验模态分析的正向链路:敲击/激振器扫频 → 测频响函数(FRF)曲线 → 峰值读固有频率、半功率带宽读阻尼、峰值之比辨振型——理论与实验在这一节握手。
老师的提醒
- 模态分析四步模板:列 $M, K$ → 解特征值问题 → 归一模态 → 模态叠加求响应。任何 $n$ 自由度题都走这条流水线。
- 模态截断的检查习惯:算完看"参与因子/有效质量"——被截掉模态的有效质量占比大时结果会骗人。
自检
- 模态关于 M 加权正交的含义与用途?
- 模态截断为什么合法?什么时候危险?
- 实验 FRF 曲线上怎么读频率、阻尼、振型?
更深一层:模态分析是"耦合系统的语法书"
多自由度系统的模态分析之所以是"语法书",因为它把任何 $n$ 自由度系统的动力学翻译成 $n$ 个独立单振子的词典:特征值问题 $(K - \omega^2M)\phi = 0$ 解出固有频率与振型、M 加权正交($\phi_i^TM\phi_j = 0$)保证模态坐标下解耦、任意响应 = 模态叠加。语法书的通用性:结构有限元(模态叠加法降维)、实验模态(FRF 拟合)、控制工程(模态截断保留低阶)、随机振动(模态 PSD 叠加)——同一语法的四种方言。
模态截断的经济学值得单独标注:有限元百万自由度只取前几十阶模态(低频激励下高阶模态贡献 $\propto 1/\omega^2$ 急剧衰减),计算量从 $O(n^3)$ 降到小矩阵求解。截断的风险检查:模态有效质量占比(被截掉的模态参与质量 < 10% 才放心)、激励频段与截断频率的关系——"截断合法"是有验收标准的省略。
推导细读:阻尼模型的"诚实边界"
比例阻尼 $C = \alpha M + \beta K$ 是模态解耦的门票,但它的物理假设(各阶模态阻尼比 $\zeta_i = \frac{\alpha}{2\omega_i} + \frac{\beta\omega_i}{2}$ 随频率平滑变化)与真实结构(阻尼集中在接头/材料局部)常有出入。非比例阻尼(局部阻尼器、材料界面)需要复模态理论(状态空间特征值问题),工程折中是"模态应变能法"(按各部件储存的应变能比例分配阻尼比)。阻尼是振动建模中最不精确的参数——行业惯例是阻尼比实测 + 10% 保守余量,而不是迷信模型。
专题深入:模态分析是"耦合系统的语法书"完整版
多自由度系统的模态分析之所以是"语法书",因为它把任何 $n$ 自由度系统的动力学翻译成 $n$ 个独立单振子的词典:特征值问题 $(K - \omega^2M)\phi = 0$ 解出固有频率与振型、M 加权正交($\phi_i^TM\phi_j = 0$)保证模态坐标下解耦、任意响应 = 模态叠加。语法书的通用性:结构有限元(模态叠加法降维)、实验模态(FRF 拟合)、控制工程(模态截断保留低阶)、随机振动(模态 PSD 叠加)——同一语法的四种方言。
模态截断的经济学值得单独标注:有限元百万自由度只取前几十阶模态(低频激励下高阶模态贡献 $\propto 1/\omega^2$ 急剧衰减),计算量从 $O(n^3)$ 降到小矩阵求解。截断的风险检查:模态有效质量占比(被截掉的模态参与质量 < 10% 才放心)、激励频段与截断频率的关系——"截断合法"是有验收标准的省略。
数学细节展开:阻尼模型的"诚实边界"完整版
比例阻尼 $C = \alpha M + \beta K$ 是模态解耦的门票,但它的物理假设(各阶模态阻尼比 $\zeta_i = \frac{\alpha}{2\omega_i} + \frac{\beta\omega_i}{2}$ 随频率平滑变化)与真实结构(阻尼集中在接头/材料局部)常有出入。非比例阻尼(局部阻尼器、材料界面)需要复模态理论(状态空间特征值问题),工程折中是"模态应变能法"(按各部件储存的应变能比例分配阻尼比)。阻尼是振动建模中最不精确的参数——行业惯例是阻尼比实测 + 10% 保守余量,而不是迷信模型。
工程案例细节:实验模态分析的"正向链路"工程版
实验模态分析的正向链路:敲击/激振器扫频 → 测频响函数(FRF)曲线 → 峰值读固有频率、半功率带宽读阻尼、峰值之比辨振型——理论与实验在这一节握手。锤击法的细节纪律:力锤的刚度决定激励带宽、参考点的选择决定模态的可辨识性、窗函数(力窗与响应窗)的选择影响谱泄漏。FRF 曲线的"三读":峰定频率、带宽定阻尼、峰的形状定振型参与——实验模态分析的全部手艺都在这三读里。
例题
两自由度:$m_1 = m_2 = 1$,$k_1 = 100$(墙-$m_1$)、$k_2 = 100$($m_1$-$m_2$),$m_2$ 自由端。求固有频率。
解答
振型:低阶 $(1, 1.618)$(同向、黄金比!)、高阶 $(1, -0.618)$(反向)。
- 悬臂链的振型天然带 Fibonacci/黄金比味道——两自由度手算就能看到它。