1. 课程导入与工程模态分析概览

当你站在一座大跨径斜拉桥上，脚下偶尔传来车辆通过时的轻微起伏；当你在风电场看到叶片在风中缓慢摆动；当你关上车门，听到车身发出低沉的“嗡”声——这些现象背后都指向同一个工程问题：结构在特定频率下会以特定形态振动。模态分析就是研究这些“特定频率”和“特定形态”的学科，而特征值计算则是求解它们的数学核心。

本课程围绕“特征值计算”与“工程模态分析”两条主线展开。前者是数值线性代数的经典问题，后者是结构动力学、机械振动、声学等领域的核心工具。二者通过一个简洁的数学关系紧密相连：Kφ = λMφ。其中 K 为刚度矩阵，M 为质量矩阵，λ 为特征值（对应固有频率的平方），φ 为特征向量（对应振型）。理解这个方程的含义、求解方法及其工程解释，是本课程的全部起点。

为什么工程结构需要“模态”视角

任何连续体在离散化之后都拥有有限个自由度，也就对应有限个固有频率与振型。以一座中等跨径斜拉桥为例，其有限元模型可能有数十万个自由度，但真正决定风致振动、车辆共振、地震响应的，往往只是前几十阶模态。工程师关心的是：哪些频率落在激励源的频带内？对应振型在哪些位置位移最大？ 例如，某斜拉桥一阶竖弯频率若为 0.28 Hz，而该地区常见风速对应的涡激振动频率恰好在 0.25～0.30 Hz 区间，就需要评估共振风险并考虑加装阻尼器。

风机叶片的情况更为典型。一台 2 MW 机组叶片长度约 50 m，其一阶挥舞频率通常在 0.5～1.0 Hz 量级。叶片旋转时产生的气动激励频率与转速相关，若叶片一阶频率与 3 倍转频接近，就可能出现共振。设计阶段通过调整铺层厚度与纤维方向改变刚度分布，本质上就是在“调”特征值——让 λ 避开激励频率。

汽车车身则涉及另一类问题：NVH（噪声、振动与声振粗糙度）。车身一阶扭转频率若偏低，在发动机怠速激励下容易产生方向盘抖动；一阶弯曲频率若与路面激励耦合，则影响乘坐舒适性。工程上常以“频率分离”为原则，要求车身低阶模态频率与主要激励频率保持足够裕度。

模态分析到底在做什么

从数学上看，模态分析就是求解广义特征值问题 Kφ = λMφ。对于无阻尼自由振动系统，假设解的形式为 x(t) = φ·sin(ωt)，代入运动方程 Mẍ + Kx = 0，即可得到上述特征值方程。λ = ω²，ω 为圆频率，φ 为对应的振型向量。每个特征对 (λᵢ, φᵢ) 描述结构的一个“固有振动模式”。

工程实践中，模态分析分为两条路径。计算模态分析（CAE）先建立有限元模型，组装 K、M 矩阵，再用 Lanczos、子空间迭代或 Arnoldi 等方法求解前若干阶特征对。实验模态分析（EMA）则通过锤击法或激振器对结构施加激励，测量响应，经频响函数估计与曲线拟合提取模态参数。两条路径互为验证：计算模型需要用实验模态修正，实验结果的物理意义需要计算模型来解释。

以锤击法测试为例，典型操作流程是：在车身或叶片表面布置加速度传感器，用力锤逐点敲击并记录力信号与响应信号，得到频响函数矩阵；通过峰值拾取或最小二乘复指数法识别固有频率、阻尼比与振型。若某阶频率在多次敲击中重复性差，往往提示传感器粘贴位置恰好靠近该阶振型节点，需要调整测点。

本课程的学习路径

课程将按“数学基础—数值方法—工程应用”三层推进。首先建立特征值问题的数学表述与矩阵性质（对称性、正定性、广义问题转化）；随后进入数值算法，重点讲清 Lanczos 迭代、子空间迭代与 QR 算法的适用场景和收敛判据；最后落到工程模态分析，包括有限元建模要点、边界条件处理、模态截断、实验模态与计算模态的相关性分析。

学习过程中建议始终带着一个具体对象：可以是一座桥、一片叶片或一个车身简化模型。每学完一种算法，就用手算小矩阵或开源工具验证一遍；每学完一个工程概念，就回到对象上问一句“这个参数变了，频率和振型会怎么变”。特征值计算不是孤立的数学练习，它是工程师判断结构动态行为的底层语言。