几何非线性有限元：从理论到工程应用的关键认知

几何非线性（Geometric Nonlinearity）在有限元软件里通常是一个复选框：Abaqus 里的 NLGEOM，ANSYS 里的 NLGEOM,ON。很多人对于要不要勾选这个选项的决策，取决于一个简单的判断——变形大不大。变形大了就勾，变形不大就不勾。

但这个判断其实是把”大位移”和”大应变”混在了一起——它们是几何非线性的两个不同来源，对计算结果的影响方式和程度也不同。把这个概念上的混淆理清楚之后，才能避免在不该开非线性的时候浪费计算资源，或者在需要开的时候漏掉关键效应。

大位移 vs 大应变——你面对的到底是什么

大位移小应变是最常见的几何非线性场景。一个典型的例子是细长梁的大挠度弯曲：结构的位移量级可以达到梁长度的一半，但材料本身的应变仍然在弹性小应变范围内（< 0.1%）。这种情况下，线性分析把平衡方程建立在未变形的初始构型上，会低估结构的刚度——因为弯曲之后力臂变短了，力造成的弯矩变小了，结构实际上比线性预测更”硬”。

大位移大应变在高分子材料、橡胶件、金属成型模拟中更常见。不光是平衡需要建立在变形后的构型上，应变张量也不能再用线性化的柯西无穷小应变。这里 Green-Lagrange 应变或对数应变（真应变）才是正确的应变度量。

Abaqus 里的 NLGEOM 选项同时处理了大位移和大应变两种情况——但 ANSYS Mechanical 里需要额外注意：打开 Large Deflection 之后，SHELL181/BEAM188 单元默认是考虑有限转动的，但如果用了低阶实体单元 SOLID185，默认仍然是小应变公式，需要单独通过 KEYOPT 切换到大应变公式。

什么时候不开几何非线性是”对的”

不是说所有计算都应该开几何非线性。以下场景完全可以线性处理，开了反而浪费时间：

模态分析。 固有频率和振型对结构刚度矩阵的线性化处理高度依赖——在未变形构型下求解广义特征值问题是标准做法。除非你在做预应力模态分析（先做非线性静力加预载、冻结刚度、再做线性模态），否则开 NLGEOM 跑模态分析会得到一个非线性振动的初值问题，而不是线性模态分析的结果。

疲劳分析。 绝大多数疲劳分析基于线性应力-应变叠加原理。如果单次加载的塑性变形很小，几何非线性和材料的非线性效应都可以忽略——因为你关心的不是绝对应力的精确值，而是应力幅值的变化趋势。

初始方案对比。 设计初期在多个方案之间做筛选时，线性分析能快速给出相对优劣的判断。这个时候开非线性反而拖慢了迭代节奏——反正不同方案之间的线性对比结论和最终的精确非线性结果，在趋势上通常是同向的。

收敛困难时的几个经验性调整

几何非线性计算的收敛失败，最常见的原因是增量步过大导致切线刚度矩阵偏离真实路径太远。Abaqus 默认的自动增量策略在大多数情况下能处理，但对于以下场景几乎必死：

接触 + 大位移。 结构的变形改变了接触面的法向方向，接触的”打开-闭合”过程可能在两个增量步之间完成。如果增量步太大，上一个增量步是接触打开、下一个增量步已经变成接触关闭——算法在中间找不到连续的切线路径。解决方法是手动限制最大增量步，让接触状态的变化在 3-5 个增量步内平滑过渡。

屈曲或失稳。 薄壳结构在屈曲临界点附近，切线刚度矩阵接近奇异，Newton-Raphson 法直接发散。Abaqus 里有一个被低估的选项是 Stabilization（自动增加一个粘性阻尼项来抑制局部不稳定），设置一个很小的耗散能分数（如 1E-6）就足以让程序跳过屈曲点附近的刚度奇异，而耗散能对最终结果的污染几乎不超出数值误差的量级。

材料软化。 如果材料模型定义了损伤或软化（如混凝土的 CDP 模型、金属的延性损伤），切线刚度可能从正变负——Newton-Raphson 在负刚度区必然发散。弧长法（Riks Analysis）在这里是唯一能拿到后屈曲路径的方法。

说到底，几何非线性的收敛问题，根源在于 “沿着切线方向走一步”的思路在切线本身剧烈变化时失效了。减小步长、增加阻尼、换弧长法——这些补救手段都是为了让切线在每一步内的变化尽可能小，从而让 Newton-Raphson 至少”假装”这个步内是线性的。理解了这一点，收敛策略就不会变成无头苍蝇一样乱调参数。