本页目录

动力 III · 多自由度与模态分析

层次：本科核心 + 业界日常工具｜ 🔗 前置：线性代数特征值问题（数学站）。真实结构有无穷多自由度。模态分析的思想极其漂亮：把复杂的耦合振动分解成一组互不干扰的单自由度振动——于是上一页的全部直觉都能直接复用。这是振动分析的核心工具，也是有限元动力学计算的标准输出。

一、多自由度系统与特征值问题

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

无阻尼自由振动 \(\mathbf{x} = \boldsymbol{\phi}e^{i\omega t}\) 代入，得广义特征值问题：

\[\left(\mathbf{K} - \omega^2\mathbf{M}\right)\boldsymbol{\phi} = 0\]

特征值 \(\omega_i^2\) → 固有频率；
特征向量 \(\boldsymbol{\phi}_i\) → 振型（模态）：系统按该频率振动时各点的相对位移形状。

\(n\) 个自由度给出 \(n\) 阶模态。

二、模态的正交性：为什么可以解耦

关键性质：不同模态关于质量矩阵与刚度矩阵正交：

\[\boldsymbol{\phi}_i^T\mathbf{M}\boldsymbol{\phi}_j = 0, \qquad \boldsymbol{\phi}_i^T\mathbf{K}\boldsymbol{\phi}_j = 0 \quad (i\neq j)\]

后果（这是模态分析的全部威力所在）：用模态坐标 \(\mathbf{x} = \sum q_i\boldsymbol{\phi}_i\) 展开，方程组解耦为 \(n\) 个独立的单自由度方程：

\[\ddot{q}_i + 2\zeta_i\omega_i\dot{q}_i + \omega_i^2 q_i = f_i(t)\]

于是上一页关于共振、阻尼、频段的一切直觉，逐模态适用。 复杂系统的振动被理解为"若干个独立单自由度系统的叠加"。

注意阻尼的前提：只有比例阻尼（\(\mathbf{C}=\alpha\mathbf{M}+\beta\mathbf{K}\)，即瑞利阻尼）才能被同一组模态解耦。真实结构的阻尼往往不是比例的（局部阻尼器、连接摩擦），此时需要复模态分析。工程上常直接假设比例阻尼并给各阶模态指定实测阻尼比——这是一个务实但需要意识到的近似。

三、模态的工程意义

① 只有少数几阶模态重要。 高阶模态频率高、且通常被激励不到。工程上常只保留前若干阶（模态截断），把上万自由度的模型降到几十个模态——这使得动力学仿真可行，也是 🔗 自动化站第七页"平衡截断降阶"的力学版本。

② 振型决定"哪里该加固"。 若某阶模态的振型在某处位移大，则：

在位移大处加质量 → 降低该阶频率；
在变形大（应变能高）处加刚度 → 提高该阶频率。

在振型的节点处加固几乎无效——这是很多"加了筋却没用"的原因。必须先看振型，再决定加哪里。

③ 振型决定"哪里该测量"。 传感器放在节点上会测不到该阶模态。这与 🔗 自动化站第七页的"能观性"是同一个概念——传感器位置决定了哪些模态可观测。

四、实验模态分析：业界的标准手段

理论计算需要模型，实验模态分析（EMA）直接从实测得到模态参数：

流程：

激励：力锤敲击（快速、适合现场）或激振器（可控、适合精细测量）；
测量：加速度计或激光测振仪采集响应；
计算频响函数（FRF） \(H(\omega) = X(\omega)/F(\omega)\)；
参数识别：从 FRF 的峰值提取固有频率、阻尼比、振型。

它的用途：

验证有限元模型（计算与实测对不上，说明模型的边界条件或连接刚度假设错了——这是最常见的原因）；
诊断现场振动问题：先测出结构固有频率，再与激励频率比对；
模型修正：用实测数据反调有限元模型参数。

一个实践要点：有限元算出的固有频率往往偏高，因为模型中的连接（螺栓、焊缝、导轨）通常被建得过于刚硬，而真实连接有柔性与间隙。"计算频率比实测高 10–30%"是常见现象，模型修正的主要工作就是调这些连接刚度。

五、连续体的振动

杆、梁、板、壳有无穷多自由度。经典解给出解析的频率与振型。

梁的横向振动（欧拉-伯努利）：

\[EI\frac{\partial^4 w}{\partial x^4} + \rho A\frac{\partial^2 w}{\partial t^2}=0 \quad\Longrightarrow\quad \omega_n = \beta_n^2\sqrt{\frac{EI}{\rho A}}\]

注意频率与刚度的平方根、与长度的平方成反比（\(\beta_n \propto 1/L\)，故 \(\omega \propto 1/L^2\)）。长度加倍，固有频率降为 1/4——细长件极易低频共振，这是长轴、悬臂刀具、机械臂的通病。

振动的"频率上限"直觉：结构越轻越刚，固有频率越高。这就是为什么高速机床追求"轻质高刚"——它要把结构的第一阶固有频率推到工作转速的谐波之上（第九页、第十八页）。

六、模态在设计中的用法

避开共振是首要目标。标准做法：

计算/实测结构的前几阶固有频率；
列出所有激励频率（转速及其谐波、齿轮啮合频率、叶片通过频率——第七页那张表）；
画坎贝尔图（Campbell diagram）：横轴转速、纵轴频率，画出各阶固有频率线与各激励谐波线（斜线），交点即为潜在共振转速；
通过改刚度/质量，把共振点移出工作转速范围；移不出去的，要么快速通过（不长时间停留），要么加阻尼。

坎贝尔图是旋转机械设计的核心工具（下一页的转子动力学会正式用到）。

七、要点

模态分析把耦合振动解耦成独立的单自由度问题——上一页的全部直觉可逐模态复用；
正交性是解耦的数学基础；真实阻尼常非比例，需注意近似；
振型决定加固位置与测点位置——节点处加固无效、测量无感；
实验模态分析是验证有限元模型的标准手段；计算频率常因连接过刚而偏高；
梁的固有频率 \(\propto 1/L^2\)，细长件极易低频共振；
坎贝尔图用于系统地避开共振转速。

下一页：把振动理论用到旋转机械上。转子动力学有一些反直觉的现象——比如转速超过临界转速后，转子反而会"自动对中"。