3. 模态分析理论核心

模态分析在数学上求解的是一个特征值问题。理解这一点，才能明白为什么求解器输出的“频率”和“振型”具有明确的物理意义，也才能正确设置求解参数、判断结果是否可信。本节从无阻尼系统出发，逐步引入阻尼、模态叠加、参与系数与有效质量，形成一条从方程到工程判据的完整链条。

一、无阻尼模态：特征值问题的原型

对无阻尼多自由度系统，运动方程为：

[M]{ẍ} + [K]{x} = {0}

其中 [M] 为质量矩阵，[K] 为刚度矩阵。假设系统做简谐振动 {x} = {φ}e^{jωt}，代入后得到：

([K] − ω²[M]){φ} = {0}

这是一个广义特征值问题。令 λ = ω²，则 [K]{φ} = λ[M]{φ}。求解得到 n 个特征值 λ_i 和对应的特征向量 {φ_i}。λ_i 的平方根即第 i 阶固有圆频率 ω_i，除以 2π 得到固有频率 f_i（Hz）；{φ_i} 即第 i 阶模态振型，它只表示各自由度相对振幅的比例，绝对大小由归一化方式决定。

特征向量具有正交性：{φ_i}ᵀ[M]{φ_j} = 0（i≠j），{φ_i}ᵀ[K]{φ_j} = 0（i≠j）。当 i=j 时，{φ_i}ᵀ[M]{φ_i} = m_i 称为模态质量，{φ_i}ᵀ[K]{φ_i} = k_i 称为模态刚度，且满足 ω_i² = k_i/m_i。正交性是模态叠加法能够成立的根本原因。

二、有阻尼模态：从实模态到复模态

实际结构总存在阻尼。引入阻尼矩阵 [C] 后方程为：

[M]{ẍ} + [C]{ẋ} + [K]{x} = {F(t)}

若采用比例阻尼假设，即 [C] = α[M] + β[K]（Rayleigh 阻尼），则阻尼矩阵也能被无阻尼振型对角化，系统仍具有实模态，第 i 阶模态阻尼比为 ζ_i = (α/ω_i + βω_i)/2。这是工程中最常用的处理方式，Ansys Mechanical 中可通过材料阻尼、常数阻尼比或 Rayleigh 系数输入。

若阻尼为一般非比例阻尼，振型将变为复数，各点之间存在相位差，此时需要求解复特征值问题，得到复模态。复模态的实部对应衰减率，虚部对应阻尼固有频率。对于大多数金属结构，阻尼比较小（常在 0.5%～3%），比例阻尼假设足够；但对于含大量橡胶、粘弹性材料或流体耦合的结构，必须考虑复模态。

三、模态叠加原理

模态叠加法的核心思想是：物理坐标下的响应可以表示为各阶模态振型的线性组合。令 {x} = [Φ]{q}，其中 [Φ] 为振型矩阵，{q} 为模态坐标。代入运动方程并左乘 [Φ]ᵀ，利用正交性可得一组解耦的单自由度方程：

q̈_i + 2ζ_iω_i q̇_i + ω_i² q_i = F_i/m_i

其中 F_i = {φ_i}ᵀ{F} 为第 i 阶模态力。这样，一个 n 自由度耦合系统被转化为 n 个独立单自由度系统，逐个求解后再叠加回物理坐标。Ansys Mechanical 的 Mode Superposition 求解器正是基于此原理，适用于线性系统、频响分析和瞬态分析。

模态叠加的截断是关键：只保留前 m 阶模态（m < n）。截断准则通常要求保留的模态在分析频段上限的 1.5～2 倍以内，或使各方向有效质量参与率达到 80%～90% 以上。

四、参与系数与有效质量

参与系数衡量某一阶模态在特定方向激励下被激发的程度。对方向向量 {r}（如 X 向平动），第 i 阶模态参与系数定义为：

Γ_i = {φ_i}ᵀ[M]{r} / m_i

有效模态质量则定义为：

m_eff,i = Γ_i² · m_i = ({φ_i}ᵀ[M]{r})² / m_i

有效质量的意义在于：它表示第 i 阶模态在 {r} 方向上“贡献”了多少可被激发的质量。将所有模态的有效质量累加，理论上等于结构在该方向的总质量。因此，工程上常用“累计有效质量比”判断模态是否取够：若某方向累计有效质量达到总质量的 90%，说明该方向的主要动态行为已被捕获。

下表给出一个简支梁示例（假设参数，仅用于说明判据）的模态参与情况：

可以看到，前两阶已贡献近 80% 的 X 向有效质量，但若要达到 90% 判据，需保留到第 5 阶。若只取前 3 阶，虽然频率覆盖看似足够，但有效质量不足，频响分析在共振区可能出现明显偏差。

五、工程实践中的关键判断

在 Ansys Mechanical 中执行模态分析时，求解器默认输出频率和振型。要判断结果是否可用，应重点检查三项：第一，累计有效质量比是否达到 90%；第二，是否存在局部模态（如某个小零件单独振动），局部模态会拉低有效质量贡献，必要时需通过网格或连接方式处理；第三，模态频率是否与分析频段匹配，若激励频率接近某阶模态，应确认该阶模态的阻尼比设置是否合理。

对于预应力模态（如旋转机械、张紧结构），需先进行静力分析得到应力刚度矩阵 [K_σ]，再求解 ([K] + [K_σ] − ω²[M]){φ} = {0}。此时频率会随载荷变化，不能直接用无预应力结果替代。

总结来看，无阻尼特征值问题给出频率和振型的数学基础；阻尼使振型可能变为复数，比例阻尼是工程简化；模态叠加将耦合系统解耦，使响应计算可行；参与系数和有效质量则提供了“取多少阶模态才够”的量化判据。四者环环相扣，构成模态分析从理论到工程应用的完整逻辑。