汽车为什么会有自己的“振动形状”？——从固有频率到模态降阶理解模态分析
看到整车的两种典型振动形态以后，我觉得模态分析其实可以先从一个很直观的问题开始理解：一辆车如果发生自由振动，它会以什么频率振，又会以什么形状振？ 一种常见的低阶整体振动形式是 1st Global Torsion Mode（第一阶整体扭转模态），车身整体以扭转为主；另一种是 1st Global Bending Mode（第一阶整体弯曲模态），车身主要发生整体弯曲。
这里实际上已经包含了 Modal Analysis（模态分析）的两个核心结果：Natural Frequency（固有频率）和 Mode Shape（振型）。固有频率回答“结构天然倾向于以多快的频率振动”，振型则回答“在这个频率下结构以什么空间形状振动”。因此，一个完整模态最好理解成一对信息：
其中 是第 阶固有频率， 是对应振型。对于汽车 NVH（Noise, Vibration and Harshness，噪声、振动与声振粗糙度）来说，这两类信息非常重要，因为道路、动力总成等外部激励如果接近某个重要固有频率，就可能明显激发与之对应的振动形态。所以低阶整体弯曲、扭转等模态，往往会成为整车和白车身动力学设计最先关注的对象。
从实际结构动力学出发，一个线性系统通常可以写成
其中 是 Mass Matrix（质量矩阵）， 是 Damping Matrix（阻尼矩阵）， 是 Stiffness Matrix（刚度矩阵）， 是外部动态激励。如果研究的是“道路或发动机激励作用以后结构到底振多大”，就需要保留右侧的 ，这属于 Forced Vibration（受迫振动）。但模态分析首先想知道的是：把外部持续激励拿掉以后，这个结构自身具有哪些基本振动特性？ 因此先令 ，得到自由振动方程
如果再暂时忽略阻尼，即令 ，就得到常规模态分析中最基础的 Undamped Free Vibration（无阻尼自由振动）方程：
这里忽略阻尼，并不是认为真实汽车没有阻尼，而是先把问题集中到质量和刚度这两个最基本的动力学属性上。对于许多轻阻尼结构，无阻尼固有频率与有阻尼振动频率本身非常接近。以单自由度系统为例，有阻尼自由振动的特征根为
对应的有阻尼振动频率为
当阻尼比 很小时， 。因此在大量常规结构 NVH 问题中，先求无阻尼模态，可以比较清楚地提取结构自身由质量和刚度决定的基本动力学骨架；而当系统存在强非比例阻尼、摩擦或非对称刚度时，则需要进一步考虑 Damped Eigenvalue（有阻尼特征值）甚至 Complex Eigenvalue（复特征值）。
无阻尼固有频率可以从运动方程中直接推出来。假设结构在某一个模态上进行 Harmonic Motion（简谐运动），可以写成
其中 描述空间振动形状，而 描述随时间的周期变化。对时间求两次导数：
代回无阻尼自由振动方程，得到
定义 Eigenvalue（特征值）
就得到标准的 Generalized Eigenvalue Problem（广义特征值问题）：
要存在 的 Nontrivial Solution（非平凡解），必须满足
求出一系列特征值 后，可以进一步得到
以及以 Hz 表示的固有频率
再将每一个 代回
就可以求得与之对应的 Eigenvector（特征向量） ，也就是 Mode Shape（振型）。所以固有频率和振型并不是两类彼此独立的结果，而是同一个特征值问题的两个方面：特征值决定“振多快”，特征向量决定“怎么振”。
求出这些特征值和特征向量以后，它们的作用并不只是生成一张固有频率表和几幅振型图。对于经典线性系统，不同固有频率对应的特征向量具有非常重要的 Orthogonality（正交性）。对于 ，有质量正交关系
以及刚度正交关系
如果采用 Mass Normalization（质量归一化），进一步令
那么把各阶振型组成 Modal Matrix（模态矩阵）
以后，就得到
同时，由特征值关系
可以进一步得到
其中
是一个对角矩阵。也就是说，原来在物理坐标下复杂、耦合的质量矩阵和刚度矩阵，经过模态变换以后，可以被同时转换成非常简单的对角形式。这就是所谓的 Diagonalization（对角化），也是模态分析真正具有强大计算价值的关键一步。
有了这组特征向量以后，还可以把它们作为一组 Modal Basis（模态基底）。如果原有限元模型具有 个物理自由度，而在关注的频率范围内只需要保留前 个重要模态，就可以构造
并用
来近似描述结构运动。其中 不再是原来每一个节点对应的 Physical Coordinates（物理坐标），而是各个模态对应的 Modal Coordinates（模态坐标）。如果一个有限元模型原来有 100 万个自由度，而目标频率范围内只需要 1000 个主要模态，那么动力学问题就有可能从百万级物理自由度降到千级模态自由度，这就是 Modal Reduction（模态降阶）的基本思想。
更进一步，如果把
代入完整动力学方程
并左乘 ，就得到
利用前面的正交性和对角化关系，可以写成
其中
这时可以看到一个很关键的事实：用无阻尼振型作为模态基底，并不意味着后续动力学分析中把阻尼和激励也丢掉了。 阻尼和外部载荷仍然存在，只是从原来的 Physical Space（物理空间）被投影到了 Modal Space（模态空间）。如果阻尼还是 Proportional Damping（比例阻尼），那么 也可以成为对角矩阵，于是整个多自由度系统就进一步分解成一组彼此独立的单自由度振子。这个问题正好是后面值得继续展开的一步：为什么正交性能让矩阵对角化，为什么对角化又能把一个复杂耦合动力学问题拆成许多简单的单自由度问题。
这样把整个逻辑串起来以后，我现在更愿意把模态分析理解成一套新的动力学坐标系：
这样再回头看整车第一阶弯曲和扭转振型，它们就不只是两幅变形图了。每一个振型实际上都可以看成描述结构动力学的一根“坐标轴”。当许多重要振型组合在一起以后，就形成了一套特别适合结构动力学的坐标系，而后续的模态叠加、频率响应以及大规模 NVH 计算，本质上都在利用这套坐标系，把原来高维、耦合的问题转换成更低维、更容易求解的问题。
以上内容是个人学习汽车 NVH 分析过程中的整理，属于学习记录性质。部分理解可能还比较粗浅，欢迎相关领域老师和同行交流指正。
