3. 几何非线性基础

在 Abaqus/Standard 与 Abaqus/Explicit 中，几何非线性并不是一个“精度选项”，而是决定平衡方程建立在哪个构型上的根本假设。线性分析默认在初始构型上建立平衡，位移和转动都假定足够小，应变与位移近似为线性关系；一旦结构出现大位移、大转动或大应变，线性化就会引入不可忽略的误差，此时必须打开几何非线性，让平衡方程在变形后的构型上迭代求解。

大位移、大转动与小应变假设

“大位移”指结构整体位移量级与自身特征尺寸可比，例如悬臂梁端部挠度达到梁长的 20% 以上；“大转动”指局部截面或单元发生显著刚体转动，即使应变很小，转动也会改变载荷方向与刚度贡献；“小应变”则指材料本身的应变仍处于弹性或小塑性范围，本构关系无需引入大应变张量。工程中大量橡胶密封、柔性铰链、薄板屈曲问题属于“大位移 + 大转动 + 小应变”组合，这正是几何非线性最典型的适用场景。

需要区分三种几何非线性来源：一是应变-位移关系的非线性（大应变），二是转动对刚度矩阵的贡献（大转动），三是载荷方向随变形更新（跟随力）。Abaqus 的 NLGEOM 开关同时激活前两者，而压力、随动载荷等跟随力行为则与载荷类型和 NLGEOM 状态共同决定。

NLGEOM 开关与求解器行为

在 Abaqus/Standard 中，NLGEOM=ON 写在 *STEP 行，表示该分析步采用几何非线性。此时求解器使用 Newton-Raphson 迭代，每轮迭代根据当前变形构型更新刚度矩阵与内力向量，直到残差满足收敛容差。若 NLGEOM=OFF，则整个分析步保持线性小变形假设，即使后续分析步打开 NLGEOM，前面步骤的构型也不会被追溯修正。

*STEP, NLGEOM=ON, INC=200
*STATIC
0.01, 1.0, 1e-5, 0.05
*BOUNDARY
FIXED, 1, 3
*CLOAD
TIP, 2, -500.0
*END STEP

在 Abaqus/Explicit 中不存在 NLGEOM 开关，因为显式求解天然基于当前构型更新，几何非线性始终开启。因此从 Standard 迁移到 Explicit 时，不必也不能再设置 NLGEOM，但要注意质量缩放、稳定时间增量与几何非线性共同影响结果。

跟随力与压力载荷方向更新

跟随力（follower force）指载荷方向随结构变形而转动。Abaqus 中压力载荷 *DLOAD 的 P 类型默认按当前构型法向施加，属于典型的跟随力；集中力 *CLOAD 默认保持全局方向不变，除非使用 *FOLLOW FORCE 或指定非零方向更新。压力载荷在 NLGEOM=ON 时随表面转动，在 NLGEOM=OFF 时固定在初始法向，这一差异对橡胶密封圈、气囊、薄膜结构影响极大。

橡胶密封圈压缩案例：几何非线性效应观察

取一个轴对称橡胶 O 形密封圈，截面直径 5 mm，安装于宽 4.2 mm 的矩形沟槽中，上下刚性板压缩量 1.0 mm。材料采用 Mooney-Rivlin 超弹性模型，C10=0.8 MPa，C01=0.2 MPa，D1=0.02 MPa⁻¹。网格使用 CAX4H 杂交单元，接触采用硬接触与罚摩擦，摩擦系数 0.2。

分析分两个分析步：第一步 NLGEOM=OFF，施加 0.2 mm 压缩；第二步 NLGEOM=ON，继续压缩至 1.0 mm。观察三个量：密封圈与沟槽侧壁的接触压力、截面最大主应变、以及压缩反力。线性步中反力随压缩近似线性增长；打开 NLGEOM 后，由于截面大转动与接触面积快速增加，反力曲线明显上翘，接触压力峰值从约 1.6 MPa 升至约 2.4 MPa，最大主应变从 0.18 增至 0.27。若忽略几何非线性，会低估密封接触压力约 30%，导致密封能力判断偏乐观。

结果判读与常见误区

打开 NLGEOM 后，应重点检查三项：一是反力-位移曲线是否出现明显非线性上翘或软化；二是接触压力分布是否随变形重新分布；三是能量平衡中伪应变能（ALLAE）与总能量之比是否低于 5%。若比值过高，说明网格畸变严重或接触定义不当，需要细化网格或调整接触刚度。

对于密封圈这类“大位移 + 大转动 + 小应变”问题，NLGEOM=ON 是必选项，压力载荷方向更新是必选项，而材料本构仍可用小应变超弹性模型。理解这三者的分工，才能在 Abaqus 中正确设置分析步、载荷与输出，避免把几何非线性误当成材料非线性来处理。