3. 几何非线性与增量迭代

在 Abaqus/Standard 中，几何非线性并不是一个可以随意勾选的“开关”，而是一整套关于应变度量、平衡方程建立位置以及载荷步推进方式的假设体系。理解这套体系，才能判断一个分析到底该用 NLGEOM=ON 还是保持线性，以及为什么打开几何非线性后收敛会突然变难。

大位移、大转动、小应变：三个独立假设

工程中常说的“几何非线性”实际上包含三个可以彼此独立出现的现象：大位移（结构整体位置显著改变，如悬臂梁端部挠度达到梁长的 10% 以上）、大转动（截面或单元发生显著刚体转动，如薄板翻卷、壳单元折叠）、小应变（材料本身仍处于弹性或小塑性应变范围，应变分量远小于 1）。Abaqus 的几何非线性选项同时处理这三者，但它们的数学来源不同：大位移改变的是平衡方程中力的方向与力臂，大转动改变的是应变—位移关系的非线性项，而小应变假设决定了本构关系是否仍可用工程应变与工程应力直接对应。

关键点在于：大转动并不要求大应变。一块薄板可以发生 90° 的刚体转动而材料应变仍小于 0.2%，此时若使用线性分析，单元会因为忽略转动引起的几何刚度变化而给出完全错误的结果。反过来，橡胶类材料可能在小位移下就出现大应变，那属于材料非线性而非几何非线性。判断是否需要 NLGEOM=ON 的实用准则是：当最大位移超过结构特征尺寸的 5%，或任何单元发生超过 5° 的转动时，就应打开几何非线性并重新评估。

Newton-Raphson 迭代：增量步内的平衡搜索

Abaqus/Standard 采用增量—迭代格式求解非线性平衡方程。每个载荷增量步内，Newton-Raphson（N-R）迭代的基本过程是：以当前构型为参考，计算内力 I(u) 与外载 P 的不平衡力 R = P − I(u)；若 R 不满足收敛判据，则用切线刚度矩阵 K_t 求解修正量 Δu = K_t⁻¹ R，更新位移后重新计算内力，直到不平衡力足够小。对于几何非线性问题，K_t 必须包含材料刚度、几何刚度（初应力刚度）和载荷刚度三部分，这也是打开 NLGEOM 后计算量明显增加的原因。

实际使用中，增量步长由 Abaqus 自动控制。若某增量步在默认 16 次迭代内未收敛，求解器会将增量步长减半重试；若连续多次成功收敛，则逐步放大步长。用户可通过 *STEP, INC=200 限制最大增量数，通过 *CONTROLS 调整迭代上限与收敛容差。对于强烈非线性问题，手动设置初始增量（如 *STATIC, STABILIZE 或固定小步长）往往比完全依赖自动控制更可靠。

收敛判据：残差与位移的双重检查

Abaqus 默认同时检查两个量：残差力（不平衡力相对于时间平均外载的比值）和位移修正量（本次迭代位移增量相对于总位移增量的比值）。两者都必须小于容差才判定收敛。默认残差容差为 5×10⁻³，位移容差为 10⁻²。对于以力加载为主的问题，残差判据更关键；对于以位移加载或接近极限承载力的屈曲问题，位移判据往往先触发。

当出现“收敛但结果异常”时，常见原因是容差过松导致迭代停在错误的平衡点附近。此时可收紧容差（如 *CONTROLS, PARAMETERS=FIELD 中设置 R_n = 1e-4），或改用位移控制加载。另一个实用技巧是查看 .msg 文件中每次迭代的残差与位移修正量，判断是残差下降缓慢（可能刚度矩阵病态）还是位移修正量震荡（可能接近分叉点）。

弧长法：越过极限点追踪后屈曲路径

当结构存在极限点（如受压圆柱壳的失稳、浅拱的跳跃）时，载荷—位移曲线会出现下降段。此时常规的载荷控制或位移控制 N-R 迭代会在极限点附近失效，因为切线刚度矩阵趋于奇异，且下降段上同一个载荷对应多个位移解。弧长法（Riks 法）通过引入弧长约束，将载荷因子也作为未知量，沿平衡路径的弧长前进，从而自然越过极限点并追踪后屈曲路径。

在 Abaqus 中使用弧长法的典型设置是：*STEP, NLGEOM=ON 配合 *STATIC, RIKS，并指定弧长增量 *RIKS 参数。初始弧长建议取总弧长的 1%～5%，最大弧长不超过 10%。弧长法输出的是载荷因子—位移曲线，用户需要关注 .sta 文件中的弧长步与载荷比例因子，以及 .fil 或 .odb 中的 LPF（Load Proportionality Factor）。当 LPF 开始下降时，说明已越过极限点进入后屈曲段。

本节核心关系图解

实操中的参数设置与结果判读

以一个受压薄壁圆柱壳的线性屈曲与后屈曲分析为例。首先用 *BUCKLE 提取前 10 阶特征值，得到临界载荷估计值。随后进行后屈曲分析：*STEP, NLGEOM=ON, INC=200，*STATIC, RIKS，*RIKS, 0.02, , , 0.1 设置初始弧长 0.02、最大弧长 0.1。加载方式采用参考载荷 *CLOAD 施加单位压力，弧长法会自动缩放载荷因子。

提交后检查 .sta 文件：若 LPF 从 1.0 逐步上升至峰值后开始下降，说明成功追踪到后屈曲路径；若 LPF 在峰值附近反复震荡且弧长步不断减小，说明模型存在局部失稳或接触突变，应改用 *STATIC, STABILIZE 并设置阻尼因子。在 .odb 中提取 U 与 RF，绘制载荷—位移曲线，与理论临界载荷对比，偏差在 10% 以内通常可接受。

总结本节要点：几何非线性由大位移、大转动、小应变三个假设共同定义，判断准则是位移超过特征尺寸 5% 或转动超过 5°；N-R 迭代在每个增量步内通过切线刚度矩阵搜索平衡，收敛判据同时检查残差力与位移修正量；弧长法通过引入弧长约束越过极限点，是追踪后屈曲路径的标准工具，但不适用于接触突变或材料软化问题。掌握这三者的适用边界与参数设置，是进行可靠非线性结构分析的基础。