欧拉-伯努利梁理论
理论与物理
最基本的梁理论
老师，“欧拉-伯努利梁”是什么？欧拉在屈曲问题中也出现过呢。
这是莱昂哈德·欧拉和丹尼尔·伯努利在18世纪建立的梁弯曲理论。是结构力学中最基本的理论，也是FEM梁单元的出发点。
基本假定
它做了哪些假定呢？
3个基本假定：
1. 平截面假定 — 变形前为平面的横截面，变形后仍保持平面
2. 正交假定 — 变形后横截面仍与梁轴线（中性轴）正交
3. 小变形 — 变形足够小
假定2看起来很重要。截面不倾斜意味着…
是的。这意味着忽略剪切变形。梁仅因弯曲而变形，剪切引起的截面倾斜为零。这是欧拉-伯努利梁的最大特点，同时也是最大的限制。
控制方程
请告诉我弯曲的微分方程。
关于挠度 $w(x)$ 的四阶常微分方程：
其中 $EI$ 是弯曲刚度，$q(x)$ 是分布荷载。
四阶微分方程！积分四次就能得到解吗？
是的。四个积分常数由四个边界条件决定。每个端部需要两个条件（位移和转角，或剪力和弯矩）。
各阶导数的物理意义：
| 导数 | 物理量 | 公式 |
|---|---|---|
| $w$ | 挠度 | |
| $w' = dw/dx$ | 转角 $\theta$ | |
| $w'' = d^2w/dx^2$ | 曲率 $\kappa = M/(EI)$ | $M = EI w''$ |
| $w''' = d^3w/dx^3$ | 剪力 | $V = -EI w'''$ |
| $w'''' = d^4w/dx^4$ | 分布荷载 | $q = EI w''''$ |
挠度微分四次就回到荷载。很优雅呢。
记住这个关系对FEM结果的验证非常有用。挠度→转角→曲率→弯矩→剪力→荷载，形成连锁。
忽略剪切变形的影响
忽略剪切变形会产生多大程度的误差呢？
梁的跨高比（$L/h$）越大，误差越小：
| $L/h$ | 剪切变形的贡献 | 欧拉-伯努利梁的精度 |
|---|---|---|
| > 20 | < 1% | 足够精确 |
| 10 〜 20 | 1 〜 5% | 实用上没问题 |
| 5 〜 10 | 5 〜 20% | 需要注意 |
| < 5 | > 20% | 不精确。应使用铁木辛柯梁 |
从 $L/h < 10$ 左右开始出现差异了呢。
一般的钢梁（H型钢）$L/h = 15 \sim 25$ 左右，所以用欧拉-伯努利梁就足够了。但对于夹芯板或短的连接梁（$L/h < 5$），则需要铁木辛柯梁理论。
FEM中的梁单元
FEM的欧拉-伯努利梁单元是什么形式？
2节点梁单元，每个节点有3个自由度（2D情况）：挠度 $w$、转角 $\theta$、轴向位移 $u$。
重要特征是形函数为埃尔米特多项式（三次多项式）。与通常的FEM单元（拉格朗日多项式）不同，它不仅将位移，还将转角也作为节点变量。这使得仅用2个节点就能精确求解弯曲的四阶微分方程。
2个节点就能解四阶方程！一个单元就精确吗？
对于等截面、均布荷载的情况，1个单元就能得到精确解。这是欧拉-伯努利梁单元的一大优点。集中荷载作用的位置需要放置节点，但除此之外，粗网格就足够了。
总结
我来整理一下欧拉-伯努利梁理论。
要点：
- 忽略剪切变形的经典梁理论 — 截面始终与中性轴正交
- $EI w'''' = q$ — 四阶常微分方程
- $L/h > 10$ 则足够精确 — 适用于细长梁
- FEM中采用埃尔米特插值 — 2节点即可得到等截面梁的精确解
- $L/h < 5$ 时使用铁木辛柯梁 — 剪切变形不可忽略
材料力学课上解过的悬臂梁挠度 $\delta = PL^3/(3EI)$，正是这个理论的解。FEM的梁单元是将这个理论离散化后的产物，原理是相同的。
欧拉梁理论的诞生
欧拉-伯努利梁理论起源于1744年莱昂哈德·欧拉所著的《De Curvis Elasticis》。“变形后截面保持平面且与梁轴线正交”这一假定最初曾引发争议，但对于细长梁（长细比L/h>10），至今仍能发挥误差在1%以内的精度。
各项的物理意义
- 惯性项（质量项）：$\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 |
数值解法与实现
梁单元的刚度矩阵
请告诉我欧拉-伯努利梁单元的刚度矩阵。
2节点，每个节点有2个弯曲自由度（$w_i, \theta_i$）时，4×4的刚度矩阵：
这个 $12, 6L, 4L^2$ 的数字是从哪里来的？
对埃尔米特形函数微分得到 $B$ 矩阵，然后积分 $\int_0^L B^T EI B \, dx$ 即可得到。因为埃尔米特形函数是三次多项式，所以曲率（二阶导）是一次式，积分可以精确执行。
加上轴向力自由度（$u_i$）后，得到与弯曲独立的2×2刚度：
将此与弯曲的 $[K_b]$ 组合，构成6×6（2D）或12×12（3D）的梁单元刚度矩阵。
等效节点荷载
分布荷载如何施加到梁单元上？
将分布荷载 $q$ 转换为等效节点荷载。均布荷载的情况：
每个节点有 $qL/2$ 的力和 $\pm qL^2/12$ 的弯矩…这是两端固定梁的固端弯矩吧！
理解得很完美。等效节点荷载等于固定端反力变号后的值。知道这个对应关系，就能直观地验证等效节点荷载是否计算正确。
求解器对应的单元名称
| 单元 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 2节点梁（EB） | CBAR | B23（2D）, B33（3D） | BEAM3（2D）, BEAM4 |
| 2节点梁（Timoshenko） | CBEAM | B21（2D）, B31（3D） | BEAM188/189 |
Nastran的CBAR和CBEAM不一样吗？
CBAR是欧拉
なった
詳しく
報告