第4章连续体振动基础

4. 连续体振动基础：梁的横向振动、板的振动、有限元法初步

各位工程师朋友，欢迎来到第四章。前面几章我们聊了单自由度、多自由度系统，那些都是离散化的模型。但现实中的结构，比如桥梁的桥面、飞机的机翼、咱们手里的电路板，它们都是连续体。说白了，质量和刚度是连续分布的，不是集中在几个点上。

这一章，我们就来啃这块硬骨头。我会从最简单的梁开始，讲到板，最后引入有限元法这个工程利器。我个人觉得，理解连续体振动，是区分“会算题”和“真懂结构”的分水岭。

梁的横向振动，是连续体振动里最经典的问题。你想想看，一根细长的杆子，你敲它一下，它就会弯来弯去地振动。这就是横向振动。

描述这个现象，最常用的模型是欧拉-伯努利梁理论。它的核心假设是：变形前垂直于轴线的截面，变形后仍然保持为平面，并且垂直于变形后的轴线。说白了，就是忽略剪切变形和转动惯量的影响。

这个模型的控制方程长这样：

EI * ∂⁴w(x,t)/∂x⁴ + ρA * ∂²w(x,t)/∂t² = f(x,t)

其中：

这个方程是四阶偏微分方程，求解起来比单自由度系统复杂得多。但它的物理意义很清晰：弯曲刚度（EI）抵抗变形，惯性力（ρA）抵抗加速度。

关键点：梁的振动是无穷多个模态的叠加。每个模态对应一个固有频率和一个振型。比如简支梁的一阶振型是半个正弦波，二阶振型是一个完整正弦波。

我在项目中遇到过一个问题：一个精密仪器的支撑梁，在某个转速下振动剧烈。一开始大家以为是电机不平衡，后来我用锤击法做模态测试，发现是梁的某阶固有频率和电机转速重合了。嗯，这就是典型的共振。解决办法很简单，要么改变梁的截面尺寸（调刚度），要么加个阻尼器。

从梁到板，维度从一维变成了二维。板可以看作是梁在宽度方向上的扩展。它的振动行为更丰富，也更复杂。

描述薄板横向振动的经典理论是基尔霍夫-勒夫板理论。它的控制方程是：

D * ∇⁴w(x,y,t) + ρh * ∂²w(x,y,t)/∂t² = f(x,y,t)

其中：

你看，这个方程比梁的方程多了一个维度。求解起来，边界条件的影响非常大。同样的板，四边简支和四边固支，固有频率能差出好几倍。

我的经验：处理板的振动问题时，我建议你先用解析公式估算一下前几阶频率。比如四边简支矩形板的固有频率公式：

ω_mn = π² * √(D/ρh) * [(m/a)² + (n/b)²]

其中 m、n 是半波数，a、b 是板的长宽。这个公式能帮你快速判断问题的大致范围。

我曾经处理过一个汽车天窗的振动问题。天窗玻璃在高速行驶时会产生低频轰鸣声。我一开始以为是风噪，后来用加速度传感器一测，发现是玻璃板的一阶模态被气流激励起来了。解决方案是在玻璃边缘贴了一层约束阻尼层，效果立竿见影。

好了，梁和板的解析解虽然漂亮，但只适用于简单几何和边界条件。现实中的结构，比如一个带加强筋的机翼、一个开孔的电路板，解析解基本算不出来。这时候，有限元法（FEM）就登场了。

有限元法的核心思想，说白了就是“化整为零，积零为整”。把连续体离散成有限个小单元，每个单元用简单的形函数近似，然后组装成整体方程求解。

对于振动问题，有限元法最终会得到一个多自由度系统的方程：

[M]{ü} + [C]{u̇} + [K]{u} = {F}

你看，这和第二章的多自由度系统方程一模一样。只不过这里的质量矩阵[M]和刚度矩阵[K]是通过单元矩阵组装得到的。

重要概念：有限元法求解振动问题的步骤：

我建议你刚开始学有限元时，不要急着用商业软件。先用一个简单的梁单元，手算一遍。比如一个悬臂梁，用两个梁单元，算它的前两阶频率。这个过程能让你深刻理解“形函数”、“单元刚度矩阵”、“组装”这些概念到底是什么意思。

避坑指南：我曾经犯过一个错误——网格划分太粗。当时算一个板的模态，前几阶频率和实验值差了30%。我以为是理论错了，折腾了两天。后来把网格加密一倍，结果就对了。记住：网格密度直接影响高阶模态的精度。一般来说，每个波长至少需要4-6个单元。

为了让你更直观地理解本章的知识结构，我画了一张图。你可以看到，我们从最简单的梁开始，逐步过渡到板，最后用有限元法解决工程实际问题。这三者之间是层层递进的关系。

这张图清晰地展示了本章的脉络。从左到右，是从简单到复杂、从理论到工具的演进。你学完这一章，应该能回答三个问题：梁怎么振？板怎么振？用有限元怎么算？

好了，这一章的内容就到这里。记住，连续体振动是结构动力学里最考验直觉的部分。多想想物理图像，少死记硬背公式。下次你看到一座桥、一块电路板，试着在脑子里想象它的振型——这就是工程师的“内功”。