梁自由振动分析
理论与物理
梁的振动 — 动态分析的起源
老师，梁的自由振动是固有频率分析最基本的问题吗？
是的。梁的自由振动是存在理论解的少数问题之一，对FEM精度验证至关重要。
控制方程
欧拉-伯努利梁的自由振动方程：
分离变量法 $w(x,t) = W(x) e^{i\omega t}$ 得到空间部分：
是四阶常微分方程呢。解是怎样的？
通解：
其中 $\beta^4 = \rho A \omega^2 / (EI)$。
四个常数 $C_1 \sim C_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 | 两端为自由端 |
悬臂梁和两端固定梁的一阶 $(\beta L)^2$ 是3.516 vs. 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~2%），但十阶模态有时会有10~30%的差异。需要精确计算到高阶模态时，必须使用铁木辛柯梁。
FEM的验证
梁的自由振动是FEM基准测试的理想问题：
1. 理论解已知且精确
2. 可改变单元数验证网格收敛性
3. 可体验EB梁单元与铁木辛柯梁单元的差异
4. 可比较梁单元与壳/实体单元
对于学习FEM的人来说是最好的练习题呢。
用FEM计算悬臂梁的一阶固有频率，并与理论值 $f_1 = 3.516/(2\pi L^2) \sqrt{EI/\rho A}$ 比较。仅此一项就能验证单元精度、边界条件、质量设置是否正确。
总结
我来整理一下梁的自由振动。
要点：
- $EI W'''' = \rho A \omega^2 W$ — 四阶ODE特征值问题
- 存在理论解 — FEM的最佳基准
- 边界条件导致频率变化6倍以上 — 首先检查
- 高阶模态需要铁木辛柯梁 — EB梁在高阶时不准确
- 记住 $(\beta_n L)^2$ 的值 — 悬臂梁: 3.516, 简支梁: $\pi^2$
Euler-Bernoulli梁的历史
“欧拉-伯努利梁理论”是莱昂哈德·欧拉（1744年）和雅各布·伯努利合作的成果。最初伯努利研究平躺梁的弯曲，欧拉将其推广。直到铁木辛柯在1921年发表包含剪切变形和旋转惯性的改进版本之前，这是唯一的梁振动理论。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动而建筑物的质量“被落下”。静力分析中此项设为零，这是“缓慢施力故加速度可忽略”的假设。冲击载荷或振动问题中绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时能感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸长更多？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，决定了刚度。常见误解：“刚度高=强度高”是不对的。刚度是“不易变形性”，强度是“不易破坏性”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（如压力、接触力）。可以这样想——桥上卡车的重量是“作用在整个内容物上的力”（体积力），轮胎压路面的力是“仅作用在表面的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却成了“压缩”——听起来像笑话，但在3D空间坐标系旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他弦试试。声音会一直持续吗？不，会逐渐变小。因为振动能量通过空气阻力和弦的内部摩擦转化为热能。汽车的减震器也是同样原理——故意吸收振动能量以提高乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
基于FEM的梁振动分析
用梁单元进行振动分析时，需要多少单元？
EB梁单元（埃尔米特插值）能精确表示弯曲的四次多项式，所以静力分析中一个单元就足够。但振动分析中为了表现模态振型需要多个单元。
要表现 $n$ 阶模态的半波长需要：
- 1阶模态: 至少4个单元
- 5阶模态: 至少20个单元
- $n$ 阶模态: 至少 $4n$ 个单元（参考）
比静力分析需要更多单元呢。
振动的模态振型是正弦/余弦的组合，模态阶次越高波长越短。要解析短波长就需要更细的网格。板和实体单元也是如此。
质量矩阵的影响
一致质量矩阵和集中质量矩阵结果会不同吗？
EB梁单元的情况：
- 一致质量 — 倾向于高估 $f$（刚度的上界定理）
- 集中质量 — 倾向于低估 $f$
- 真值介于两者之间
哪个更准确？
通常一致质量更准确。不过如果单元数足够多，差异就很小。用4个EB梁单元求悬臂梁的一阶固有频率，一致质量能达到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$ 个
相关主题
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告