欧拉-伯努利梁理论
欧拉-伯努利梁理论的理论基础
最基本的梁理论
老师，"欧拉-伯努利梁"是什么？欧拉在座屈理论中也出现过呢。
莱昂哈德·欧拉和丹尼尔·伯努利在18世纪建立的梁弯曲理论。它是结构力学中最基本的理论，也是有限元梁单元的出发点。
基本假设
做了什么假设？
三个基本假设：
1. 平面保持假设 — 变形前平面的截面，变形后仍保持平面
2. 直交假设 — 变形后截面仍与梁轴（中性轴）正交
3. 微小变形 — 变形足够小
假设2看起来很重要。截面不倾斜意味着…
对。忽略剪切变形。梁仅由弯曲变形，剪切引起的截面倾斜为零。这是欧拉-伯努利梁最大的特点，同时也是最大的限制。
支配方程
请教我弯曲的微分方程。
挠度 $w(x)$ 的四阶常微分方程：
其中 $EI$ 是弯曲刚度，$q(x)$ 是分布荷载。
四阶微分方程！积分四次就能得到解？
对。四个积分常数由四个边界条件确定。每端需要2个条件（位移和转角，或剪切力和弯矩）。
每个微分的物理意义：
| 微分 | 物理量 | 公式 |
|---|---|---|
| $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''''$ |
挠度微分四次就回到荷载。很优雅。
这个关系如果记住，对验证有限元结果很有帮助。挠度→转角→曲率→弯矩→剪切力→荷载，形成链条。
忽略剪切变形的影响
忽略剪切变形会产生多大的误差？
梁的跨度/高度比（$L/h$）越大，误差越小：
| $L/h$ | 剪切变形贡献 | 欧拉-伯努利的精度 |
|---|---|---|
| > 20 | < 1% | 精度足够 |
| 10 〜 20 | 1 〜 5% | 工程上可接受 |
| 5 〜 10 | 5 〜 20% | 需要注意 |
| < 5 | > 20% | 不精确。应使用Timoshenko梁 |
从 $L/h < 10$ 开始就会有差异呀。
一般钢梁（H形钢）$L/h = 15 \sim 25$，用欧拉-伯努利就够了。但对于夹芯板或短连接梁（$L/h < 5$），需要用Timoshenko梁理论。
FEM中的梁单元
有限元中的欧拉-伯努利梁单元长什么样？
两节点梁单元，每节点3个自由度（2D情况）：挠度 $w$、转角 $\theta$、轴向位移 $u$。
重要特点是形状函数采用Hermite多项式（三次多项式）。与普通有限元单元（Lagrange多项式）不同，挠度和转角都是节点变量。因此仅用2个节点就能精确求解弯曲的四阶微分方程。
两个节点就能解四阶方程！一个单元就精确了吗？
等截面、均匀分布荷载的情况下，一个单元就能得到精确解。这是欧拉-伯努利梁单元的重大优势。集中荷载作用处需要放节点，除此之外网格可以很粗。
总结
请总结一下欧拉-伯努利梁理论。
要点：
- 忽略剪切变形的经典梁理论 — 截面始终与中性轴正交
- $EI w'''' = q$ — 四阶常微分方程
- $L/h > 10$ 时精度足够 — 适用于细长梁
- 有限元中用Hermite补间 — 两个节点给出等梁的精确解
- $L/h < 5$ 时用Timoshenko梁 — 剪切变形不可忽略
材料力学中求解的悬臂梁挠度 $\delta = PL^3/(3EI)$ 正是这个理论的解。有限元梁单元就是这个理论的离散化版本，原理完全相同。
欧拉梁理论的诞生
欧拉-伯努利梁理论源于欧拉1744年发表的《De Curvis Elasticis》。"变形后截面保持平面且与梁轴正交"的假设起初引起争议，但对于细长梁（L/h>10），精度至今仍在误差1%以内。
欧拉-伯努利梁理论的数值计算方法
梁单元刚度矩阵
请教我欧拉-伯努利梁单元的刚度矩阵。
两节点、每节点弯曲自由度2个（$w_i, \theta_i$）的4×4刚度矩阵：
这个 $12, 6L, 4L^2$ 的数字从哪来的？
用Hermite形状函数求导得到B矩阵，然后进行 $\int_0^L B^T EI B \, dx$ 积分得到的。Hermite形状函数是三次多项式，曲率（二阶导数）是一次多项式，所以积分可以精确执行。
加上轴力自由度（$u_i$）后，轴力独立的2×2刚度：
将轴向刚度[K_a]和弯曲刚度[K_b]组合，构成2D的6×6或3D的12×12梁单元刚度矩阵。
等效节点荷载
分布荷载怎样应用到梁单元上？
将分布荷载 $q$ 转换为等效节点荷载。均匀分布荷载的情况：
每个节点 $qL/2$ 的力和 $\pm qL^2/12$ 的力矩……这是两端固定梁的固定端反力吧！
完美理解！等效节点荷载等于固定端反力的符号相反。了解这个对应关系，就能直观验证等效节点荷载计算是否正确。
求解器的单元名称
| 单元 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 两节点梁（EB） | CBAR | B23（2D）, B33（3D） | BEAM3（2D）, BEAM4 |
| 两节点梁（Timoshenko） | CBEAM | B21（2D）, B31（3D） | BEAM188/189 |
Nastran中的CBAR和CBEAM有区别吗？
CBAR是欧拉-伯努利梁（不含剪切变形），CBEAM是Timoshenko梁（含剪切变形+翘曲处理）。CBEAM功能更强，但对细长梁，CBAR就足够了。
Abaqus中B33是欧拉-伯努利梁，B31是Timoshenko梁。三维梁单元推荐用B31（Timoshenko）。B33的使用场景很少。
截面定义
梁单元的截面怎样定义？
主要参数：
- $A$ — 截面积（用于轴力）
- $I_y, I_z$ — 截面二次矩（围绕主轴）
- $J$ — 扭转常数（Saint-Venant扭转）
- $C_w$ — 翘曲常数（仅CBEAM）
定义方式有两种：
1. 直接输入 — 手动输入 $A, I_y, I_z, J$ 的值
2. 从截面形状自动计算 — 给定H形钢、箱形、圆形等形状参数，自动计算截面特性。Abaqus中用 *BEAM SECTION，Ansys中用 SECTYPE/SECDATA
手动输入的风险是什么？
单位错误最常见。$I$ 是长度的四次方，mm系则为 mm⁴。mm和m混用会差 $10^{12}$ 倍。从截面形状自动计算更安全。
总结
梁单元的数值方法，请总结一下。
要点：
- Hermite补间给出精确弯曲解 — 等梁1个单元就够
- 等效节点荷载 = 固定端反力符号相反 — 可直观验证
- CBAR是EB梁，CBEAM是Timoshenko梁 — Nastran中的区别
- 截面诸元从形状自动计算更安全 — 避免单位错误
- 三维梁推荐用Timoshenko梁 — EB梁使用场景少
Hermite补间的威力
欧拉-伯努利梁单元因为要求挠度和转角都连续，采用三次Hermite多项式作为形状函数。这个四自由度定式在1965年被Turner等人应用于NASA宇宙结构，仅用6个单元对梁的解析解的误差就控制在0.3%。这在当时是革命性的成果。
欧拉-伯努利梁理论的实务应用
梁单元的实务应用
梁单元在工程中怎样应用？
梁单元是结构分析中使用最广泛的单元。建筑框架、桥梁、机械轴、管道……到处都用。
框架分析
建筑骨架分析怎样做？
钢结构和混凝土框架分析是梁单元最大的应用领域。柱和梁用梁单元表示，接合部用节点处理。
要点：
- 柱梁接合部 — 刚接、铰接、半刚接的区分
- 荷载流向 — 荷载→楼板→梁→柱→基础的传递
- 支撑的处理 — 桁架单元（仅轴力）或弹簧单元
楼板的面刚度怎样建模？
三种方案：
1. 刚楼板假设 — 各层节点用刚体链接连接。最简单
2. 用壳单元建立楼板模型 — 精确考虑面刚度
3. 用斜支撑表示支撑效果 — 用支撑代替楼板的面内剪切
梁单元结果的解释
梁单元的结果能像实体单元那样读吗？
完全不同。梁单元结果输出为截面力（轴力 $N$、剪切力 $V$、弯矩 $M$、扭矩 $T$）。应力需要从截面力计算：
应力是后处理才计算的呀。
对。很多求解器会自动输出梁的应力，但那是基于上述公式的名义应力，不包含局部应力集中（圆角、螺栓孔、焊接）。应力集中的评估需要用壳或实体的子模型。
梁单元的局限性
梁单元分析不了什么情况？
| 情况 | 梁单元可以分析？ | 替代方案 |
|---|---|---|
| 弯曲+轴力 | ○ | |
| 扭转 | △（CBAR无翘曲） | CBEAM或壳 |
| 局部座屈 | × | 壳单元 |
| 接合部的详细应力 | × | 实体单元 |
| 开孔腹板 | × | 壳单元 |
| 复合截面（组合梁） | △（用等截面） | 壳+实体 |
梁单元是"把握整体"，壳和实体做"详细"，这样分层次吗？
正是。结构分析的工作流通常是"梁框架→壳子模型→实体局部模型"的分层结构。
实务检查清单
使用梁单元的检查清单请给我。
截面主轴错误好像很严重。
H形钢的强轴和弱轴搞反，截面二次矩差10倍以上。挠度也变10倍，马上就能看出问题。必须用手计算核对一下位移量级。
桥梁设计中的实务应用
明石海峡大桥（1998年完成，主塔高度298m）的设计分析中，主塔初期刚性评估采用欧拉-伯努利梁模型，约500个单元，几小时内就能试算风荷载下的水平变位。与3D壳模型对比，位移误差仅约2%。
欧拉-伯努利梁理论的软件对比
梁单元工具
梁单元框架分析有哪些工具？
通用有限元的梁单元之外，还有结构设计专用的框架分析软件。
结构设计专用软件
| 软件 | 地区 | 特点 |
|---|---|---|
| SS3/SuperBuild | 日本 | 建筑一体化结构计算。包含保有水平耐力计算 |
| MIDAS Gen | 韩国/全球 | 建筑土木兼用。支持各国规范 |
| ETABS | 美国 | 建筑框架专用。CSI的旗舰产品 |
| SAP2000 | 美国 | 通用框架分析。也支持桥梁 |
| ROBOT Structural Analysis | Autodesk | 与Revit集成。BIM流程 |
| STAAD.Pro | Bentley | 通用一般结构分析 |
SS3在日本建筑设计中很常见。
SS3（超高层对应版本为SuperBuild/SS7）是完全符合日本建筑基准法的一体化结构计算程序。建设许可申请需要的保有水平耐力计算功能内置。日本建筑设计事务所的事实标准。
通用有限元和专用设计软件怎样分工？
专用设计软件内置规范合规性检验。梁的许可应力、座屈耐力、挠度限制……全自动检验。通用有限元只输出应力，规范检验需手工进行。
因此：
- 建筑设计实务 → SS3, ETABS, MIDAS Gen（包含规范检验）
- 非标结构的详细分析 → 通用有限元（Nastran, Abaqus, Ansys）
- 耐震研究 → OpenSees
- 教育 → SAP2000, MASTAN2（柔活且直观）
开源软件
有开源的框架分析软件吗？
选型指南
总结一下？
梁单元是最简单的单元，但工具选择最多的领域呢。
梁单元框架分析是建筑土木的日常业务，效率化需求大。工具的选择往往由"用什么规范"决定。
各求解器的梁单元实现
NASTRAN的CBEAM单元继承自1963年NASA原始版本，能明确指定截面剪切中心与中性轴的偏移。而Abaqus B23单元作为三维Euler-Bernoulli梁采用三节点（二次补间），曲率表达精度比CBEAM更高。
欧拉-伯努利梁理论的先端研究
梁理论的先端研究
欧拉-伯努利梁是18世纪的理论，现在还有研究题目吗？
梁理论本身已完善，但其扩展和新材料应用很活跃。
功能梯度材料（FGM）梁
材料特性沿截面方向连续变化的FGM梁的弯曲理论。中性轴偏离几何中心，通常的弯曲公式 $M = EI\kappa$ 不能直接用。FGM梁的有限元定式化在研究中。
纳米梁（非局部弹性理论）
纳米尺度也能用梁理论吗？
碳纳米管和石墨烯薄片的力学可以用非局部弹性理论（Eringen理论）基础的欧拉-伯努利梁描述。普通弹性论中加入"非局部参数"，考虑分子间力的长程效应。
其中 $e_0 a$ 是非局部参数。变成六阶微分方程。
18世纪的理论竟然用于纳米技术！
概率梁分析
将材料特性和荷载作为随机变量的概率有限元法（SFEM）的梁单元在研究中。用于结构可靠性评估，定量化 $E$ 和 $I$ 的波动如何传播到挠度和应力分布。
共旋转公式
梁的大变形怎样处理？
共旋转公式（co-rotational formulation）是大变形梁的标准方法。每个单元设有随动的局部坐标系，单元内仍用小变形梁理论。全局上能应对大旋转、大位移。
OpenSees的非线性梁单元和Abaqus的B31单元都采用共旋转定式。地震下钢框架的大变形崩溃分析必不可少。
总结
梁理论的先端研究，请总结。
18世纪的理论从纳米到宏观，从确定性到概率，至今仍在结构力学中心。
剪切锁定问题的回避
把欧拉-伯努利单元当Timoshenko梁用会产生"剪切锁定"，高细长比时刚性矩阵过刚。1981年Prathap和Bhashyam证明了Reduced Integration法可以消除二节点Timoshenko单元的锁定，即使L/h=1000的细长梁也能保证精度。
欧拉-伯努利梁理论的故障排除
梁单元故障
梁单元分析常见故障是什么？
梁单元设置简单，反而容易忽视错误。
挠度与理论值不符
简支梁的挠度与 $\delta = 5qL^4/(384EI)$ 不符。
检查项：
1. $I$ 是否正确 — 强轴/弱轴混淆、单位错误（mm⁴ vs m⁴）
2. 支持条件 — 本意是两端铰接，实际上变成单端固定了吗。检查旋转自由度的约束
3. 荷载 — 线荷载（N/m）和集中力（N）混淆。等效节点荷载是否正确
4. 网格 — 集中荷载作用点有节点吗（没有节点的话荷载传递错误）
5. 剪切变形 — 用了Timoshenko梁单元，剪切变形分会额外增加挠度
5号是盲点呀。Timoshenko梁用欧拉-伯努利理论值比对，当然不符…
正是。Abaqus的B31（Timoshenko梁）用欧拉-伯努利理论值比对，会差几个百分点。要么用B33（欧拉-伯努利），要么用Timoshenko的理论解。
截面力图不正常
剪切力图间断了。
集中荷载作用点处剪切力间断是正常的。梁单元的截面力是按单元输出，同一节点的左右单元值可能略有不同。分布荷载的情况，单元边界处也会有小波动。
整体的弯矩图和剪力图怎样绘制？
各单元端部值接起来画。有限元后处理工具通常自动做这个，但梁的BMD/SFD显示功能有没有取决于求解器的后处理器。Abaqus/CAE能绘制梁的截面力图。Nastran的f06文件需手工绘制。
刚体链接的处理
梁和其他单元接合时的注意点？
梁单元有挠度和旋转两个自由度，但实体单元只有位移。直接连接会旋转自由度不受约束。
对策：
- RBE2/RBE3（Nastran）、*COUPLING/MPC（Abaqus） 耦合梁和实体的自由度
- 用壳单元中介 — 壳有旋转自由度，自然连接
RBE2和RBE3的区别？
另篇详述，简言：
- RBE2 — 刚体结合。独立节点运动以刚体方式传到从属节点
- RBE3 — 荷载分配。独立节点的力按权重分配到从属节点。不增加刚度
座屈分析的座屈荷载不正确
用梁单元做座屈分析，侧向座屈没出现。
欧拉-伯努利梁单元（CBAR、B33等）没有翘曲自由度，不能正确表现侧向座屈（弯扭座屈）。要评估侧向座屈：
- Timoshenko梁单元（CBEAM、B31OS等） — 有翘曲自由度
- 用壳单元建模 — 侧向座屈自动呈现
总结
梁单元故障对应，请总结。
梁单元项目少，反而一个错误影响很大呢。
对。梁单元的截面诸元和支持条件2个要素把对，结果的可信度非常高。这两个要素确保好就是梁单元分析的全部。
线性梁在大变形时的界限
线性欧拉梁单元超过梁长度约5%时，几何非线性影响不可忽略。1980年代NASA航天飞机主翼桁初期分析忽视这一点，热变形评估产生约8%的低估，在NASA-TM-89865中被指出。