3. 特征值问题的数学表述

在工程模态分析中，我们真正关心的往往不是某个瞬时响应，而是结构“天生”具有的振动特性：它倾向于以哪些频率振动、这些频率对应的振型是什么形状。要把这种直觉转化为可计算的数学问题，核心一步就是把振动微分方程整理成特征值问题。本节围绕标准与广义特征值问题的形式、振动微分方程的导出、质量矩阵与刚度矩阵的物理含义，以及边界条件如何改变问题规模与解，建立后续数值求解的数学基础。

一、从标准特征值问题到广义特征值问题

线性代数中的标准特征值问题写作 A x = λ x，其中 A 为 n×n 方阵，λ 为特征值，x 为非零特征向量。它的几何含义是：矩阵 A 作用在向量 x 上，只改变其长度（缩放 λ 倍），不改变方向。求解即令 (A − λI)x = 0 有非零解，等价于 det(A − λI) = 0，展开得到关于 λ 的 n 次特征多项式。

结构动力学中的问题通常不是这个形式，而是广义特征值问题：

K φ = λ M φ        （λ = ω²）

等价形式： (K − λ M) φ = 0
非零解条件： det(K − λ M) = 0

其中 K 为刚度矩阵，M 为质量矩阵，φ 为振型向量，λ 为特征值，其物理意义是圆频率的平方 ω²。与标准问题相比，多了一个矩阵 M 参与“加权”。当 M 可逆时，左乘 M⁻¹ 可得 M⁻¹K φ = λ φ，形式上化为标准问题，但 M⁻¹K 一般不再对称，数值上并不推荐这样直接处理；工程软件多采用 Cholesky 分解 M = LLᵀ 后转为对称标准问题，或直接使用 Lanczos、子空间迭代等广义求解算法。

二、振动微分方程的导出

以多自由度集中质量系统为例。设位移向量为 u(t)，对每个质量块列牛顿第二定律，弹性恢复力与相对位移成正比，阻尼力与速度成正比，外力为 F(t)，整理成矩阵形式即得：

M ü(t) + C u̇(t) + K u(t) = F(t)

这是结构动力学的一般运动方程。M、C、K 分别为质量、阻尼、刚度矩阵，ü、u̇、u 为加速度、速度、位移向量。模态分析关注系统自身的自由振动特性，因此令外力和阻尼为零（无阻尼自由振动）：

M ü(t) + K u(t) = 0

假设解具有简谐形式 u(t) = φ sin(ωt + θ)，则 ü(t) = −ω² φ sin(ωt + θ)。代入并约去公共因子 sin(ωt + θ)，得到 −ω² M φ + K φ = 0，即 K φ = ω² M φ。至此，一个动力学问题被转化为纯代数特征值问题——这正是模态分析的关键一步：时间被消去，剩下的是关于空间振型的代数方程。

三、质量矩阵与刚度矩阵的物理意义

质量矩阵描述系统对加速度的惯性抵抗。集中质量模型下它通常是对角阵，对角元为各自由度上的集中质量或转动惯量；一致质量矩阵则由形函数积分得到，含非对角耦合项，反映相邻自由度惯性上的相互影响。M 应为对称正定矩阵，其二次型 φᵀMφ 对应动能相关的广义质量。

刚度矩阵描述系统对位移的弹性抵抗。K 对称且一般半正定；其元素 K_ij 的物理含义是：使第 j 个自由度产生单位位移、其余自由度固定时，在第 i 个自由度上所需施加的力。二次型 φᵀKφ 对应应变能相关的广义刚度。特征值可写成 Rayleigh 商：

λ = ω² = (φᵀ K φ) / (φᵀ M φ)

该式直观说明：刚度越大、质量越小，频率越高；反之频率越低。这也解释了为什么轻质高刚度结构（如碳纤维构件）固有频率偏高，而重而软的结构频率偏低。

四、边界条件的影响

边界条件在数学上通过两种方式进入问题：一是减少自由度，二是改变刚度矩阵的构成。以悬臂梁为例，固定端三个方向位移与转角均为零，这些自由度在离散模型中直接划去，矩阵规模随之缩小；而自由端无约束，对应自由度保留。若把同一根梁改为两端简支，则两端的竖向位移被约束，但转角自由，得到的刚度矩阵与悬臂情形完全不同，固有频率也随之改变。

一个便于记忆的规律是：增加约束通常提高各阶固有频率。直观理解是约束限制了结构变形方式，使其“更硬”。例如一根钢梁，自由—自由状态第一阶弯曲频率可能只有几十赫兹，而两端固支后可提高数倍。此外，边界条件还会影响振型形状：悬臂梁一阶振型为根部不动、端部摆幅最大的弯曲形态，而简支梁一阶为半正弦波。工程中若发现实测频率明显低于仿真值，首先要排查的就是边界条件是否被过度理想化——实际螺栓连接往往既非完全固支也非完全自由，而是具有一定柔度的弹性约束。

五、一个可手算的示例

设两自由度系统 m₁ = m₂ = 1 kg，k₁ = k₂ = 1000 N/m，一端固定。其质量矩阵与刚度矩阵为：

M = [1  0]      K = [ 2000  -1000]
    [0  1]          [-1000   1000]

求解 det(K − λM) = 0：

(2000 − λ)(1000 − λ) − 1000² = 0
λ² − 3000λ + 1 000 000 = 0
λ = (3000 ± √(9 000 000 − 4 000 000)) / 2
  = (3000 ± 2236) / 2
λ₁ ≈ 382 /s²  → ω₁ ≈ 19.5 rad/s → f₁ ≈ 3.11 Hz
λ₂ ≈ 2618 /s² → ω₂ ≈ 51.2 rad/s → f₂ ≈ 8.14 Hz

对应振型可由 (K − λM)φ = 0 求得：一阶振型两质量同向运动、幅值比约 1:1.62；二阶振型反向运动。若把第二端也固定（增加约束），则系统变为两个独立弹簧—质量单元，频率分别约为 5.03 Hz 与 5.03 Hz（此处两质量相同），可见约束方式对结果影响显著。该示例为教学假设参数，实际工程需按真实材料与几何建模。