1. 几何非线性大变形悬臂梁仿真

问题定义与建模思路

本案例研究一根矩形截面悬臂梁在自由端集中载荷作用下的几何非线性大变形响应。梁的几何与材料参数如下：长度 L = 500 mm，截面宽 b = 20 mm，高 h = 10 mm，弹性模量 E = 210000 MPa，泊松比 ν = 0.3。自由端施加竖向集中力 P = 50 N。

先用欧拉-伯努利梁理论估算小变形挠度：

I = b·h³/12 = 20 × 10³ / 12 = 1666.67 mm⁴
δ_small = P·L³ / (3·E·I)
        = 50 × 500³ / (3 × 210000 × 1666.67)
        ≈ 5.95 mm

挠度与梁长之比约为 5.95/500 ≈ 1.2%，看似很小，但本案例的目的是通过对比验证：当载荷逐步放大到使挠度超过梁长的 10% 时，几何非线性效应如何显著改变梁端挠度与根部应力分布。因此实际分析中会施加多级载荷，观察小变形假设何时失效。

有限元模型搭建

在 Abaqus/Standard 中建立二维平面应力模型（也可用三维实体，但二维足以展示几何非线性效应且计算高效）。关键设置如下：

• 部件：二维平面、可变形、壳/平面应力单元。绘制矩形 500 × 10 mm。
• 材料：线弹性，E = 210000 MPa，ν = 0.3。本案例只关注几何非线性，材料保持线弹性以隔离变量。
• 截面：平面应力，厚度 20 mm（对应梁宽 b）。
• 网格：CPS4R 四节点双线性平面应力单元，沿长度方向布置 100 个单元，沿高度方向 4 个单元，共 400 个单元。大变形分析中网格密度对弯曲精度影响较大，长度方向单元过少会导致自锁或精度不足。
• 边界条件：左端面完全固定（U1 = U2 = UR3 = 0）。
• 载荷：自由端右端面施加竖向集中力，通过参考点耦合或直接施加在端面节点上。使用 Amplitude 或直接分多个分析步逐级加载。

分析步与 NLGEOM 设置

创建两个静力通用分析步（Static, General），分别对应小变形与几何非线性：

1. Step-1（小变形）：Nlgeom = Off，施加载荷 P = 50 N。
2. Step-2（几何非线性）：Nlgeom = On，施加载荷从 50 N 逐步增加到 500 N，使用固定增量步长 0.1，共 10 个增量步。

在 Step 模块中，通过 Step → Edit → Basic → Nlgeom 切换开关。注意：NLGEOM 一旦在某个分析步开启，后续分析步默认继承，但可以在后续步中关闭。

结果对比与判读

提取自由端中面节点的竖向位移 U2，以及固定端根部上下表面的轴向应力 S11。下表给出预期趋势（示例数据，用于说明判读逻辑，非实测值）：

判读要点：随着载荷增大，几何非线性使梁端挠度小于线性预测值，这是因为大变形下梁的轴向刚度参与承载，且载荷方向随变形发生偏转，等效弯曲力臂减小。根部应力也略低于线性结果。若差异超过 5%，说明小变形假设已不适用。

关键操作与验证清单

• 确认 Step-2 的 Nlgeom 已设为 On，且增量步长足够小（建议 ≤ 0.1），避免大变形下收敛困难。
• 在 Field Output 中请求 U（位移）和 S（应力），并在 History Output 中记录自由端节点的 U2 和根部单元的 S11。
• 对比 Step-1 与 Step-2 在相同载荷 P = 50 N 下的结果，两者应几乎一致，验证 NLGEOM 在微小变形下不引入额外误差。
• 绘制载荷-位移曲线，观察 NLGEOM 曲线在线性曲线下方且随载荷增大偏离加剧，这是几何非线性的典型特征。
• 检查根部应力分布：NLGEOM 下中性轴可能发生偏移，上下表面应力不再严格对称。