手算核心方法（二）：能量法（瑞利法）求解基频

上一章我们聊了刚度法，说白了就是硬碰硬地列平衡方程。但有些结构，比如变截面梁、带集中质量的复杂系统，列方程能把你写到手抽筋。这时候，能量法就派上用场了。

我个人习惯把能量法叫做「偷懒的艺术」。它不关心每个质点的瞬时平衡，只盯着系统的总能量。你想想看，一个结构振动时，动能和势能互相转换，在最大位移处势能最大、动能为零；在平衡位置动能最大、势能为零。抓住这个规律，基频就出来了。

瑞利法的核心思想

瑞利法基于一个很朴素的假设：我假设一个振型，然后让系统在这个假设振型下做简谐振动。根据能量守恒，最大动能等于最大势能，频率就解出来了。

公式长这样：

ω² = U_max / T*_max

其中：

• U_max —— 系统最大势能（应变能）
• T*_max —— 系统最大参考动能（不含 ω² 的部分）

为什么叫「参考动能」？因为真正的最大动能 T_max = ½ ω² T*_max，我们把 ω² 提出来，等式两边一约，就得到频率了。

操作步骤：四步走

我在项目中用过不下几十次瑞利法，总结下来就四个步骤：

1. 假设振型 —— 选一个满足几何边界条件的形状函数 Y(x)。比如简支梁用正弦波，悬臂梁用抛物线。
2. 写势能 —— 对于弯曲梁，U_max = ½ ∫ EI (Y'')² dx。对于弹簧质量系统，U_max = ½ k Y_max²。
3. 写动能 —— T*_max = ½ ∫ m Y² dx（分布质量）或 ½ Σ m_i Y_i²（集中质量）。
4. 代入公式 —— ω² = U_max / T*_max，开方得频率。

实战案例：简支梁基频

拿最简单的简支梁开刀。梁长 L，抗弯刚度 EI，线质量 m。假设振型为正弦波：

Y(x) = sin(πx / L)

这个形状满足两端位移为零、弯矩为零的边界条件。

第一步：算势能

Y''(x) = -(π²/L²) sin(πx/L)
U_max = ½ ∫₀ᴸ EI [Y''(x)]² dx
      = ½ EI (π⁴/L⁴) ∫₀ᴸ sin²(πx/L) dx
      = ½ EI (π⁴/L⁴) · (L/2)
      = (π⁴ EI) / (4 L³)

第二步：算参考动能

T*_max = ½ ∫₀ᴸ m [Y(x)]² dx
       = ½ m ∫₀ᴸ sin²(πx/L) dx
       = ½ m · (L/2)
       = m L / 4

第三步：求频率

ω² = U_max / T*_max
   = (π⁴ EI / 4L³) / (m L / 4)
   = π⁴ EI / (m L⁴)

ω = π² √(EI / m L⁴)

这个结果跟精确解一模一样！因为正弦波就是简支梁的真实基频振型。但如果你换个假设，比如抛物线，结果就会偏大一点。

集中质量系统的处理

实际工程中，梁上经常挂着设备、电机这些集中质量。瑞利法处理起来很顺手。

假设梁上在位置 x₁, x₂, …, xₙ 处有集中质量 m₁, m₂, …, mₙ，梁本身分布质量为 m(x)。那么参考动能写成：

T*_max = ½ ∫ m(x) [Y(x)]² dx + ½ Σ m_i [Y(x_i)]²

势能不变，还是梁的弯曲应变能。代入公式就完事了。

瑞利法的误差特性

我用表格总结一下不同假设振型的误差情况，方便你心里有数：

结构类型	假设振型	基频误差	说明
简支梁	正弦波	0%	恰好是精确振型
简支梁	抛物线	+0.5%	工程上可接受
悬臂梁	静力挠曲线	+0.1%	推荐使用
悬臂梁	抛物线	+1.5%	误差稍大
两端固支梁	静力挠曲线	+0.3%	效果不错

你看，只要振型选得靠谱，误差都在 2% 以内。对于手算来说，这个精度已经非常好了。

知识结构图

下面这张图把瑞利法的脉络理清楚了：

避坑指南

瑞利法用起来简单，但坑也不少。我把自己踩过的坑列出来：

• 振型不满足几何边界条件 —— 比如悬臂梁假设的振型在固定端斜率不为零，算出来的频率会严重偏大。我刚开始学的时候就犯过这个错。
• 忘了加集中质量的动能 —— 前面说过了，大质量设备必须考虑，否则误差大到离谱。
• 势能只算了一部分 —— 有些结构既有弯曲变形又有剪切变形，如果只算弯曲势能，频率会偏大。对于深梁，剪切变形不能忽略。
• 振型选得太粗糙 —— 比如用直线去近似悬臂梁的振型，误差能到 20% 以上。至少用抛物线或者静力挠曲线。

瑞利法的适用场景

不是所有问题都适合用瑞利法。我根据自己的经验画了个表：

场景	推荐程度	理由
变截面梁基频	★★★★★	手算几乎唯一选择
带集中质量的梁	★★★★★	处理集中质量很方便
复杂边界条件	★★★★	只要振型选对，精度不错
高阶频率	★★	误差大，不推荐
连续梁多跨	★★★	可以算，但振型难选

说白了，瑞利法就是用来快速拿基频的。你想想看，在方案阶段，你连有限元模型都没建好，手算一个基频就能判断结构会不会共振。这个价值，做工程的都懂。