梁的自由振动解析
梁自由振动的理论基础
梁的振动 — 动力学分析的基础
先生，梁的自由振动是固有振动数解析中最基本的问题吗？
是的。梁的自由振动是存在理论解的少数几个问题之一，对FEM的精度验证至关重要。
控制方程
欧拉-伯努利梁的自由振动方程：
用分离变量法 $w(x,t) = W(x) e^{i\omega t}$ 求解空间部分：
这是4阶常微分方程。解是什么形式？
通解为：
其中 $\beta^4 = \rho A \omega^2 / (EI)$。
4个常数 $C_1 \sim C_4$ 由4个边界条件确定。非平凡解的存在条件（振动数方程）可得固有振动数。
各边界条件下的振动数
第 $n$ 阶固有振动数 $f_n = (\beta_n L)^2 / (2\pi L^2) \sqrt{EI/(\rho A)}$：
| 边界条件 | $(\beta_1 L)^2$ | $(\beta_2 L)^2$ | $(\beta_3 L)^2$ | 模式形状 |
|---|---|---|---|---|
| 悬臂梁 | 3.516 | 22.03 | 61.70 | 末端最大位移 |
| 单跨支撑 | $\pi^2 = 9.870$ | $4\pi^2 = 39.48$ | $9\pi^2 = 88.83$ | $\sin(n\pi x/L)$ |
| 两端固定 | 22.37 | 61.67 | 120.9 | 两端为零 |
| 自由-自由 | 22.37 | 61.67 | 120.9 | 两端自由 |
悬臂梁和两端固定梁的1次 $(\beta L)^2$ 分别为3.516和22.37…相差6倍以上。振动数是 $\sqrt{6} \approx 2.5$ 倍吗？
因为 $f \propto (\beta L)^2$，所以22.37/3.516 = 6.36倍。边界条件可使振动数变化6倍以上。当FEM结果与理论值不符时，首先检查边界条件。
蒂莫申科梁的振动
蒂莫申科梁的振动如何？
由于包含剪切变形和转动惯性，在高阶模式中与EB梁会产生差异：
蒂莫申科梁的振动数始终较低。当 $f \cdot h / c_s$（$c_s$ = 剪切波速）较大时差异明显。
高阶模式的差异越来越大吗？
是的。第1阶模式的差异很小（1～2%），但第10阶模式可能产生10～30%的差异。当需要精确求解高阶模式时需要采用蒂莫申科梁。
FEM的验证
梁的自由振动是FEM的最优基准问题：
1. 理论解精确已知
2. 可通过改变单元数验证网格收敛
3. 可体验EB梁单元与蒂莫申科梁单元的差异
4. 可比较梁单元与壳/实体单元
对于学习FEM的人来说最佳的练习问题。
用FEM计算悬臂梁的1次固有振动数，与理论值 $f_1 = 3.516/(2\pi L^2) \sqrt{EI/\rho A}$ 比较。仅通过这一步就能全面检查单元精度、边界条件和质量设置是否正确。
总结
整理梁的自由振动。
要点：
- $EI W'''' = \rho A \omega^2 W$ — 4阶ODE的特征值问题
- 存在理论解 — FEM的最佳基准
- 边界条件可使振动数相差6倍以上 — 首先检查
- 高阶模式需用蒂莫申科梁 — EB梁在高阶不准确
- 记住 $(\beta_n L)^2$ 值 — 悬臂梁：3.516，单跨梁：$\pi^2$
欧拉-伯努利梁理论的历史
"欧拉-伯努利梁理论"源于莱昂哈德·欧拉（1744年）与雅各布·伯努利的合作。伯努利最初从侧向挠度研究梁的弯曲，欧拉对其进行了普遍化。直到蒂莫申科在1921年加入剪切变形和转动惯性的改进版本前，这一直是唯一的梁振动理论。
梁自由振动的数值计算方法
FEM梁的振动解析
用梁单元进行振动解析时，单元数需要多少？
EB梁单元（埃尔米特插值）能精确表示弯曲的4次多项式，所以静力学中1个单元就够了。但在振动解析中需要多个单元来表示模式形状。
要表示第 $n$ 阶模式的半波长，需要：
- 1阶模式：最少4个单元
- 5阶模式：最少20个单元
- 第 $n$ 阶模式：最少 $4n$ 单元（经验值）
比静力学需要更多单元。
振动的模式形状由正弦/余弦组合而成，高阶模式的波长较短。要解析短波长需要细网格。这对板和实体也是同样的规律。
质量矩阵的影响
一致质量和集中质量的结果会不同吗？
在EB梁单元中：
- 一致质量 — 倾向于高估 $f$（刚性上界定理）
- 集中质量 — 倾向于低估 $f$
- 真值在两者之间
哪个更准确？
通常一致质量更精确。但单元数足够多时两者的差异很小。用EB梁单元4个的悬臂梁1次固有振动数，一致质量可达到0.1%以内的精度。
求解器的设置
```
SOL 103
CEND
METHOD = 10
BEGIN BULK
CBAR, ...
EIGRL, 10, , , 10
```
```
*BEAM SECTION, SECTION=RECT, ELSET=beam
0.01, 0.001
*STEP
*FREQUENCY, EIGENSOLVER=LANCZOS
10, ,
*END STEP
```
梁单元的固有振动数解析很简洁。
梁的振动"设置简单且能直接与理论解比较"，是理想的FEM练习。推荐新接触的求解器首先进行梁振动解析来验证设置。
总结
整理梁振动的数值方法。
要点：
- 第 $n$ 阶模式至少需 $4n$ 单元 — 比静力学需更多单元
- 一致质量更精确 — $f$ 倾向高估但精度高
- EB梁单元4个就达到1阶模式0.1%精度 — 效率极高
- FEM初期验证最优 — 用新求解器首先做梁振动验证
固有振动数的解析解与边界条件
两端固定梁的固有振动数由 $f_n=(\beta_n L)^2/(2\pi L^2)\sqrt{EI/\rho A)$ 给出，其中 $\beta_n L$ 是超越方程 $\cos(\beta_n L)\cosh(\beta_n L)=1$ 的解（第1阶=4.730）。两端简支梁的 $\beta_n L=n\pi$ 则可显式求得。与FEM结果比较时，单元数10个以上的1阶模式误差可控制在0.1%以内。
梁自由振动的实务应用
梁振动的实务应用
梁的自由振动解析在实务中如何应用？
配管振动的例子
与泵连接的配管跨度 $L = 3$ m、外径 $D = 100$ mm、壁厚 $t = 5$ mm、钢管。
如果泵转速为1750 rpm，则振动数为29.2 Hz…共振！
正好与泵的振动数一致！
对策：缩短跨度（增加支撑点）或增加配管的质量降低振动数。手计算就能做出这样的判断，这是梁振动理论的强大之处。
实务检查清单
梁振动解析的检查清单是什么？
"流体的附加质量"很容易被忽视。
装满水的配管的有效密度可能是钢管本身的近2倍。若在配管振动中忽略管内流体的质量，固有振动数会相差 $\sqrt{2} \approx 1.4$ 倍。
桥梁的固有振动数与行人共振
行人过桥时的步行频率为1.6～2.4 Hz（行走），2.4～4.0 Hz（跑步）。桥梁的1阶固有振动数落在此范围时会产生共振，2000年伦敦千禧桥开通第2天因行人横摇无法站立的事件著名。通过安装质量阻尼器（TMD）将频率从0.8 Hz外移解决了问题。
梁自由振动的软件对比
梁振动的工具
梁振动分析可以用什么工具？
从手算到FEM有多种选择。
| 工具 | 特点 |
|---|---|
| 手工计算（$f_n$ 公式） | 最快。最适合概算 |
| Excel/Python | 公式的参数化计算 |
| CAESAR II | 配管跨度的振动评估 |
| 通用FEM | 支持复杂边界条件、集中质量、变截面 |
很多情况下手算就够了。
匀截面梁用手算就够了。有变截面、集中质量、非匀匀边界条件时才需要FEM。
选型指南
梁振动"手算是第一工具"。FEM仅在复杂情况使用。
手算的理论解成为FEM结果的"验尺"。不知理论解就只用FEM做振动分析，就像无地图旅行一样。
COMSOL中的梁固有值分析
COMSOL的Beam物理界面可选择Euler-Bernoulli和Timoshenko，自动计算截面特性（I、A、J等）然后求固有值。在MEMS传感器设计中需要预测共振频率精确到0.1%，Timoshenko修正和满足 L/h>20的网格设置是精度保证的关键。
梁自由振动的先端研究
梁振动的先端主题
梁振动有先端研究吗？
梁本身是古典的，但新应用不断。
MEMS共振器
MEMS（微电机系统）的悬臂是梁自由振动的典型。用于质量传感器（生物传感器）和频率标准（共振器）。通过皮克克级质量变化检测固有振动数的改变。
纳米级振动
碳纳米管和石墨烯的振动用非局部弹性理论（Eringen理论）描述为EB/蒂莫申科梁。分子间长程力使经典梁理论给出不同的振动数。
能量收获
压电悬臂梁振动可从环境振动发电。把梁的固有振动数调匹配环境振动周期实现共振，最大化能量转换效率。
总结
梁振动的先端研究，总结。
18世纪的梁理论支撑21世纪的纳米技术和能源技术。
梁大变形非线性振动的频率转移
发生大变形时梁的弯曲刚性随变形而变，振动频率随振幅而变的"非线性振动"出现。振幅达梁长的10%以上时固有振动数最多上升20%的硬化型非线性。在MEMS纳米梁中此效应为主导，振动频率可作为位移传感器。
梁自由振动的故障排查
梁振动的故障
梁振动解析常见的故障有什么？
一次固有振动数与理论值不符
检查项目（优先级顺序）：
1. 密度 — $f \propto 1/\sqrt{\rho}$。确认单位
2. 边界条件 — 铰链（旋转自由）vs 固定（旋转约束）导致振动数大幅改变
3. 截面特性（$I$） — 强轴/弱轴错误
4. 梁理论类型 — EB梁 vs 蒂莫申科梁有差异（特别对短粗梁）
5. 单元数 — 高阶模式的单元数不足导致精度下降
蒂莫申科梁的结果与EB理论不一致
这是正常结果。蒂莫申科梁包含剪切变形，给出的振动数比EB梁低。验证时应让 $L/h$ 增大，看蒂莫申科梁结果是否收敛到EB梁理论值。
高阶模式不准确
高阶模式波长短，网格易不足。第 $n$ 阶模式至少需 $4n$ 单元。增加单元数检查收敛。
总结
梁振动故障对应，总结。
FEM梁解析高阶模式不准确时
FEM中高阶固有模式精度下降的主要原因是网格不足和数值积分次数选择不当。经验规则是目标最高模式次数为N时，单元数需最少N×4个。另外缩减积分单元（Abaqus B21R等）对低阶模式好但高阶精度差，推荐使用完全积分单元B22。
更详细
错误