第 8–9 章 · 连续体振动与近似

杆、轴、梁、膜、板——无限多质点的连续体有自己的精确振动理论（偏微分方程），而工程上常把它离散化回多自由度。这一节是振动学的收官：连续与离散在此合龙。

连续体的精确理论

连续体与多自由度的对应：每个振型 = 一个模态，正交性（$\int\rho\phi_i\phi_j dx = 0$）与模态叠加法原样成立——离散理论是连续理论的采样版。

• 杆的纵振：$\frac{\partial^2 u}{\partial t^2} = c^2\frac{\partial^2 u}{\partial x^2}$（$c = \sqrt{E/\rho}$）——波动方程 + 边界条件（固定端 $u=0$、自由端 $\partial u/\partial x = 0$）给出无穷阶固有频率：两端自由杆 $\omega_n = \frac{n\pi c}{l}$；
• 轴的扭转：同一方程，$c = \sqrt{G/\rho}$；
• 梁的弯曲（Euler-Bernoulli）：$EI\frac{\partial^4 y}{\partial x^4} + \rho A\frac{\partial^2 y}{\partial t^2} = 0$——四阶方程、四个边界条件组合（简支/固支/自由），悬臂梁 $\beta_1 = 1.875$（一阶频率 $\omega_1 = \beta_1^2\sqrt{EI/(\rho A l^4)}$）；细长比假设是它的适用边界（剪切与转动惯量被忽略，厚梁需 Timoshenko 修正：频率下调）；
• 膜与板：二维推广，板因泊松效应出现"无简并"的方形板频率分裂。

离散化三法：Rayleigh、Ritz、有限元

• 瑞利法（估基频）：假设振型 $\phi$，用能量比 $\omega^2 = \frac{\int EI\phi''^2 dx}{\int\rho A\phi^2 dx}$——基频上界（假设振型偏刚）。用静挠度曲线当试函数，误差常在百分之几内；
• 瑞利-里茨法（多模态）：$\phi = \sum a_i\psi_i$，变分得广义特征值问题——前几阶频率的上界与振型近似；
• 有限元：Ritz 的自动化（单元形函数即 $\psi_i$）——结构动力学软件的引擎。三法共同的思想：在有限维子空间里找最好的近似。

老师的提醒

• 记住两个"一阶频率"手算公式备用：简支梁 $\omega_1 = \pi^2\sqrt{EI/(\rho A l^4)}$、悬臂梁用 $\beta_1^2 = 1.875^2$——现场估算的第一反应。
• 瑞利法的魂：软的试函数估不出上界误差，刚的假设给上界——想要准，让试函数"能弯得出来"（满足几何边界条件即可）。

自检

• 欧拉-伯努利梁与 Timoshenko 梁的差别与适用范围？
• 瑞利法为什么给基频上界？
• 连续体振型的正交性怎么写、怎么用？

更深一层：连续体振动的"离散-连续合龙"

连续体（杆/轴/梁/膜/板）的振动方程是偏微分方程，但其振型与频率的结构和多自由度完全同构：振型正交性（$\int\rho A\phi_i\phi_j dx = 0$）、模态叠加、无穷阶频率（连续体的"模态密度"随频率增长）。欧拉-伯努利梁 $\omega_1 = \beta_1^2\sqrt{EI/(\rho A l^4)}$（简支 $\beta_1 = \pi$、悬臂 $\beta_1 = 1.875$）是手算第一反应。Timoshenko 修正（剪切变形 + 转动惯量）在厚梁/高频时下调频率——"细长比假设"的适用边界是梁理论的第一体检项。

连续体与离散的合龙点：有限元 = 连续振型的里茨近似——单元形函数就是试函数，网格加密 = 子空间扩张。连续体理论为有限元提供"真值参照"（简单形状的精确解用于验证单元）。

推导细读：瑞利法的"上界魔术"与试函数手艺

瑞利法 $\omega^2 = \frac{\int EI\phi''^2dx}{\int\rho A\phi^2dx}$ 给基频上界（假设振型比真实"刚"），误差极小的秘诀：试函数只需满足几何边界条件（位移/转角），不必满足力边界——"能弯得出来"比"弯得完全正确"重要。静挠度曲线作试函数是工程速算标配（误差常在百分之几）。瑞利商的极值性（基频 = 商在所有允许函数上的最小值）是这一切的理论根：商的最小化自动"纠偏"，假设得越好纠得越少。

在真实世界里：从手算到软件的"验证链"

振动工程师的验证链：手算（瑞利法/单自由度折算——量级与方向）→ 简化模型（梁/板解析解——趋势与参数敏感度）→ 有限元（全模型——细节与耦合）→ 实验模态（锤击/激振器——真实频率阻尼振型）。四层之间的每一次对账（手算 vs 有限元差 < 20%、有限元 vs 实验差 < 10%）才是"分析可信"的证据链——跳层的人不是快，是在裸奔。这条验证链也是从本书通向工程实践的最后一段路：理论给你算得快，验证链给你算得对。

专题深入：连续体振动的"离散-连续合龙"完整版

连续体（杆/轴/梁/膜/板）的振动方程是偏微分方程，但其振型与频率的结构和多自由度完全同构：振型正交性（$\int\rho A\phi_i\phi_j dx = 0$）、模态叠加、无穷阶频率（连续体的"模态密度"随频率增长）。欧拉-伯努利梁 $\omega_1 = \beta_1^2\sqrt{EI/(\rho A l^4)}$（简支 $\beta_1 = \pi$、悬臂 $\beta_1 = 1.875$）是手算第一反应。Timoshenko 修正（剪切变形 + 转动惯量）在厚梁/高频时下调频率——"细长比假设"的适用边界是梁理论的第一体检项。

连续体与离散的合龙点：有限元 = 连续振型的里茨近似——单元形函数就是试函数，网格加密 = 子空间扩张。连续体理论为有限元提供"真值参照"（简单形状的精确解用于验证单元）。

数学细节展开：瑞利法的"上界魔术"与试函数手艺完整版

瑞利法 $\omega^2 = \frac{\int EI\phi''^2dx}{\int\rho A\phi^2dx}$ 给基频上界（假设振型比真实"刚"），误差极小的秘诀：试函数只需满足几何边界条件（位移/转角），不必满足力边界——"能弯得出来"比"弯得完全正确"重要。静挠度曲线作试函数是工程速算标配（误差常在百分之几）。瑞利商的极值性（基频 = 商在所有允许函数上的最小值）是这一切的理论根：商的最小化自动"纠偏"，假设得越好纠得越少。

工程案例细节：从手算到软件的"验证链"工程版

振动工程师的验证链：手算（瑞利法/单自由度折算——量级与方向）→ 简化模型（梁/板解析解——趋势与参数敏感度）→ 有限元（全模型——细节与耦合）→ 实验模态（锤击/激振器——真实频率阻尼振型）。四层之间的每一次对账（手算 vs 有限元差 < 20%、有限元 vs 实验差 < 10%）才是"分析可信"的证据链——跳层的人不是快，是在裸奔。这条验证链也是从本书通向工程实践的最后一段路：理论给你算得快，验证链给你算得对。