2. Abaqus材料非线性分析基础

在Abaqus中做非线性分析，最先要厘清的不是“怎么点按钮”，而是“非线性到底来自哪里”。结构仿真中的非线性通常来自三个源头：几何非线性、材料非线性和边界条件非线性（如接触）。本节聚焦前两者，并顺带把求解器选择、增量步与收敛、本构输入方式串成一条完整的认知链。

几何非线性与材料非线性：两种完全不同的“非线性”

几何非线性指的是结构发生大位移、大转动或大应变后，刚度矩阵必须随构型更新。典型场景是薄壁件失稳后屈曲、橡胶大变形、柔性机构翻转。此时平衡方程必须在变形后的构型上建立，Abaqus通过打开 NLGEOM=ON 来激活。它的本质是“几何描述变了”，材料本身可能仍是线弹性的。

材料非线性则指应力与应变不再满足线性关系，即使位移很小也可能发生。金属屈服后的塑性流动、橡胶的超弹性、混凝土开裂、蠕变松弛都属于此类。它的本质是“本构关系变了”，几何可以仍按小变形处理。

两者可以独立存在，也可以同时出现。例如一个金属薄板冲压成形：板料发生大位移（几何非线性），同时进入塑性（材料非线性），此时必须同时打开几何非线性和定义塑性本构。判断是否需要几何非线性，可以看最大位移与特征尺寸之比是否超过约5%，或是否出现大转动；判断是否需要材料非线性，则看工作应力是否超过屈服强度、材料是否为橡胶/泡沫等。

隐式与显式求解器：不是“哪个更好”，而是“哪个更合适”

Abaqus/Standard采用隐式积分，每个增量步内通过Newton-Raphson迭代求解平衡方程，适合静力、低频、接触规模适中的问题。它的优势是精度高、可处理强非线性，但每个增量步都要组装并求解整体刚度矩阵，内存和计算量随模型规模增长快，且可能因不收敛而中断。

Abaqus/Explicit采用显式中心差分，不需要迭代求解整体方程，而是由前一时刻的状态直接推进到下一时刻。它天然适合高速冲击、爆炸、复杂接触、大变形破坏等高度非线性问题，但要求时间增量步小于稳定极限，通常需要大量小步推进。对于准静态问题，可以用质量缩放或加载速率控制来提高效率，但必须监控动能与内能之比，一般建议动能不超过内能的5%～10%。

选择逻辑可以概括为：静力、低频、需要高精度且模型规模可控 → Standard；高速、强接触、大变形、破坏、规模巨大 → Explicit。如果一个问题既可以用Standard也可以用Explicit，优先用Standard，因为它的结果更“干净”，不需要额外控制惯性效应。

增量步与收敛：非线性求解的节奏

非线性分析不能一步到位，必须把总载荷或总位移分成若干增量步，逐步加载。每个增量步内，Standard通过迭代寻找满足平衡的解。收敛判据通常包括残差力与时间平均力之比、位移修正量与增量位移之比等，默认容差约为0.5%。如果迭代次数超过上限仍不收敛，Abaqus会尝试减小增量步重试；若多次减小仍失败，分析终止。

实际调试中，常见的不收敛原因包括：材料数据在软化段出现负刚度、接触定义过刚、网格畸变严重、增量步初始值过大。应对策略是：先检查材料曲线是否合理，再适当减小初始增量步（如从0.1降到0.01），打开自动增量控制，必要时引入阻尼或改用Explicit。

材料本构的数学表达与输入方式

以最常见的金属塑性为例，Abaqus中常用的本构描述包括：弹性段用杨氏模量E和泊松比ν；塑性段用真实应力-真实塑性应变数据对。注意，Abaqus塑性数据要求输入的是真实应力和塑性应变，而不是工程应力和总应变。转换关系为：真实应力 = 工程应力 × (1 + 工程应变)，真实应变 = ln(1 + 工程应变)，塑性应变 = 真实应变 − 真实应力/E。

例如，某低碳钢工程应力-应变数据为：屈服点工程应力250 MPa、工程应变0.0012；颈缩前某点工程应力400 MPa、工程应变0.15。转换后：真实应力约460 MPa，真实应变约0.1398，塑性应变约0.1375。这些数据对需要按从小到大的塑性应变顺序输入，且第一行塑性应变必须为0，对应屈服应力。

*Material, name=Steel
*Elastic
210000., 0.3
*Plastic
250., 0.
300., 0.02
400., 0.08
460., 0.1375

对于超弹性橡胶，常用Mooney-Rivlin或Ogden模型，输入的是材料常数而非应力-应变曲线；对于蠕变，则用Norton-Bailey等幂律模型，输入A、n、m等参数。无论哪种本构，关键都是：先确认材料试验数据的类型，再选择对应的Abaqus材料模型，最后按正确的变量含义输入。输入错误的本构数据，比不收敛更隐蔽，也更危险。