悬臂梁弯曲（集中载荷）
理论与物理
概述
老师，我听说在悬臂梁自由端施加集中载荷的问题是V&V验证中的经典案例，实际上它是怎么被使用的呢？
悬臂梁弯曲是业界广泛用作FEA验证起点的基准问题。因为存在自由端挠度 $\delta = PL^3/(3EI)$、固定端最大应力 $\sigma_{max} = PLc/I$ 的精确解，可以定量检查数值解法的实现精度。NAFEMS的入门基准问题集也收录了这个问题。
因为有精确解所以可以进行“答案核对”，这就是代码验证（Code Verification）的核心吗？
正是如此。ASME V&V 10-2006将代码验证定位为确认分析代码数学正确性的阶段。使用具有精确解的问题来展示数值解与精确解的一致性是第一步。悬臂梁作为入门点非常合适，在欧拉-伯努利梁理论假设成立的范围内，使用一维梁单元可以得到精确一致的结果。
控制方程
那么请告诉我具体的方程。
基于欧拉-伯努利梁理论的挠曲线方程如下。
其中 $P$ 是自由端载荷，$L$ 是梁长度，$E$ 是杨氏模量，$I$ 是截面惯性矩。由固定端 $x=0$ 处 $w=0$, $w'=0$，自由端 $x=L$ 处弯矩和剪力为零的边界条件唯一确定。
应力在哪里达到最大？
弯矩在固定端最大 $M_{max} = PL$，因此最大弯曲应力出现在固定端的最外缘。
$c$ 是中性轴到最外缘的距离，$Z$ 是截面模量。对于矩形截面，$I = bh^3/12$，$c = h/2$。
如果使用铁木辛柯梁理论会有什么变化？
铁木辛柯梁考虑了剪切变形。会增加由剪切变形引起的挠度增量 $\delta_s = \kappa PL/(GA)$。$\kappa$ 是剪切修正系数（矩形截面为 $5/6$）。当跨高比 $L/h$ 大于10时，剪切变形的贡献小于1%，因此欧拉-伯努利理论就足够了。对于 $L/h < 5$ 的深梁，必须使用铁木辛柯梁或3D实体单元。
基准验证数据
我想用具体数值进行比较，请告诉我参数设置。
标准设置为 $L = 1$ m，$b = 0.1$ m，$h = 0.05$ m，$P = 1000$ N，$E = 200$ GPa，$\nu = 0.3$。此时理论值为 $\delta_{tip} = 0.160$ mm，$\sigma_{max} = 240$ MPa。
| 单元类型 | 网格 | DOF | δ_tip [mm] | σ_max [MPa] | 位移误差 [%] |
|---|---|---|---|---|---|
| BEAM2（线性梁） | 10单元 | 66 | 0.160 | 240.0 | 0.00 |
| QUAD8（二次壳） | 10×2 | 1,260 | 0.160 | 239.5 | 0.00 |
| HEX8（线性实体） | 40×8×4 | 15,120 | 0.155 | 228.1 | 3.13 |
| HEX20（二次实体） | 20×4×2 | 15,120 | 0.160 | 239.2 | 0.00 |
| TET10（二次四面体） | 自动 | ~25,000 | 0.159 | 237.5 | 0.63 |
为什么只有HEX8的误差比较大？
HEX8如果使用完全积分会发生剪切自锁，梁会表现得比实际更刚硬。在弯曲主导的问题中，低阶六面体单元本质上不利，除非使用减缩积分或B-bar法，否则收敛缓慢。二次单元因为中间节点能准确表现弯曲变形模式，即使在粗糙网格下也能保持高精度。
收敛性的理论依据
网格细化时收敛速度有理论依据吗？
有的。根据FEM的误差估计定理（由Céa定理导出的先验误差估计），$p$ 阶单元的能量范数误差以 $O(h^p)$ 的速度减小。也就是说，线性单元将单元尺寸减半，误差大约减半；二次单元则大约减至1/4。在实际的网格收敛中能否再现这个理论收敛率，是代码验证的本质。
收敛率偏离理论值的情况发生在什么时候？
典型情况是存在应力奇点时。悬臂梁的固定端实际上不是应力奇点，但根据固定约束的实现方式，有时会产生局部应力集中，导致收敛率下降。应对方法是意识到圣维南原理，在距离约束端足够远的位置进行评估，或者使用分布约束。
各项的物理含义
- 守恒量的时间变化项：表示所考虑物理量的时间变化率。在稳态问题中为零。【形象比喻】给浴缸放水时，水位随时间上升——这个“单位时间的变化速度”就是时间变化项。关闭阀门水位保持恒定的状态就是“稳态”，此时时间变化项为零。
- 通量项（流束项）：描述物理量的空间输运和扩散。大致分为对流和扩散两种。【形象比喻】对流是“河流水流运送小船”那样，物体随流动被运走。扩散是“墨水在静止水中自然扩散”那样，物体因浓度差而移动。这两种输运机制的竞争支配着许多物理现象。
- 源项（生成/消失项）：表示物理量局部生成或消失的外力/反应项。【形象比喻】在房间里打开暖气，该处就“生成”了热能。化学反应消耗燃料，质量就“消失”。表示从外部注入系统的物理量的项。
假设条件与适用范围
- 连续介质假设在空间尺度上成立
- 材料/流体的本构关系（应力-应变关系、牛顿流体定律等）在适用范围内
- 边界条件在物理上合理且在数学上正确定义
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 特征长度 $L$ | m | 需与CAD模型的单位制保持一致 |
| 特征时间 $t$ | s | 瞬态分析的时间步长需考虑CFL条件和物理时间常数 |
数值解法与实现
有限元公式化
用FEM求解悬臂梁时，通常使用哪种单元，怎么用？
根据目的区分使用。梁单元直接离散化欧拉-伯努利理论，因此能以最少的自由度达到精确解。为了验证目的而使用壳或实体单元时，基本方法是基于Galerkin法的弱形式离散化。
单元刚度矩阵通过数值积分计算。
$B$ 是应变-位移矩阵，$D$ 是材料刚度矩阵，$J$ 是雅可比矩阵。
整体刚度方程是 $[K]\{u\} = \{F\}$ 对吧？线性静力分析的话直接用直接法一次求解吗？
没错。对于悬臂梁这种规模的问题，直接法（Cholesky分解）没有问题。当DOF超过数万时，预处理共轭梯度法（PCG）的内存效率更高。这个问题的本质在于确认单元公式化的精度而非求解器选择，因此应更关注单元类型和网格密度。
单元选择的实现指南
请告诉我实际求解器中的设置。
积分方案的选择有什么影响？
完全积分的HEX8（2×2×2 Gauss点）在弯曲问题中会发生自锁。减缩积分（1×1×1）可以消除自锁，但有零能模式（沙漏模式）的风险。B-bar法或EAS（增强假定应变）法可以避免自锁同时抑制沙漏模式。Abaqus的 C3D8I（非协调模式）也
相关主题
なった
詳しく
報告