2. 大变形分析的理论基石

在 Abaqus 中做小变形分析时，我们习惯用工程应变 ε=(L−L₀)/L₀ 和 Cauchy 应力 σ 描述问题，因为位移足够小，变形前后的构型几乎重合，应变与应力的定义不会引起歧义。但一旦结构进入大变形、大转动或大应变状态——例如橡胶密封圈压缩 40%、金属板料冲压翻边、软体驱动器弯曲——变形前后构型差异显著，必须区分参考构型（未变形，坐标 X）与当前构型（变形后，坐标 x）。此时“应变”和“应力”不再是唯一量，而是一族随构型与转动行为不同而不同的度量。本节要建立的，正是这套度量体系及其相互换算关系。

变形梯度：一切应变度量的源头

描述变形的核心量是变形梯度张量 F，定义为当前坐标对参考坐标的偏导：

F_iJ = ∂x_i / ∂X_J      x = φ(X, t)

F = I + ∂u/∂X           （u 为位移场）

F 把参考构型中的微元线矢量 dX 映射为当前构型中的 dx：dx = F·dX。它同时包含了刚体转动与真实变形两部分信息，这正是大变形问题的麻烦所在：直接用 F 构造的应变会混入转动，导致“刚体转动也产生应变”的错误。因此需要从 F 中剥离转动，得到纯变形度量。

体积变化由 Jacobian 决定：J = det(F)。J=1 表示不可压缩（橡胶、塑性流动常近似如此），J>1 表示体积膨胀，J<1 表示压缩。Abaqus 的杂交单元（如 C3D8H）正是通过独立插值压力来约束 J，从而处理近不可压大变形。

Green-Lagrange 应变：参考构型下的客观度量

把 F 做极分解 F = R·U（R 为转动，U 为右伸长），可以证明右 Cauchy-Green 张量 C = Fᵀ·F = U² 完全消去了转动 R。基于 C 定义Green-Lagrange 应变：

E = ½ (C − I) = ½ (FᵀF − I)

分量形式：E_IJ = ½ (∂u_I/∂X_J + ∂u_J/∂X_I + ∂u_k/∂X_I · ∂u_k/∂X_J)

最后一项是位移梯度的二次项，正是它让 E 在大转动下仍保持客观：刚体转动时 E=0。E 定义在参考构型上，与第二类 Piola-Kirchhoff 应力 S 能量共轭，是 Abaqus/Standard 中几何非线性（NLGEOM=ON）默认采用的应变度量之一，适合金属塑性、总拉格朗日格式。

对数应变：当前构型下的“真应变”

Green-Lagrange 应变在单轴拉伸下与工程应变的关系为 E = ε + ε²/2，当 ε=0.5 时 E=0.625，与真实伸长偏离明显。工程实践中更常用对数应变（Hencky 应变）：

ε_ln = ln(λ)，  λ = L / L₀  （单轴拉伸）

张量形式：ε_ln = ½ ln(FᵀF) = ln(U)

对数应变具有可加性：分两步拉伸 λ₁、λ₂，总对数应变等于 ln(λ₁)+ln(λ₂)，而工程应变不满足。它定义在当前构型，与 Cauchy 应力 σ 能量共轭，是 Abaqus/Explicit 及橡胶超弹性、大应变塑性中更自然的度量。下表对比三种应变在单轴拉伸下的数值差异：

可见 λ=2 时三者相差超过一倍，选错度量会直接导致本构响应失真。

客观应力率：消除刚体转动的影响

应力同样面临转动污染问题。Cauchy 应力 σ 是当前构型的真实应力，但在大转动下其时间导数 σ̇ 不客观——刚体转动也会让 σ̇ 非零。若本构关系写成 σ̇ = C:ε̇，就会把转动误判为变形。解决办法是采用客观应力率（共旋率），把应力先“随材料一起转”，再求导：

Jaumann 率：  σ̊ = σ̇ − W·σ + σ·W
Green-Naghdi 率：σ̊ = σ̇ − Ω·σ + σ·Ω
Truesdell 率：  σ̊ = σ̇ − L·σ − σ·Lᵀ + (tr L)σ

其中 W 是旋率张量（速度梯度的反对称部分），Ω 由 R 的转动率给出。Abaqus 在金属大变形塑性中默认使用 Jaumann 率配合对数应变，在超弹性中则直接以应变能函数对 C 求导，回避了应力率问题。选择哪种率取决于材料模型：金属塑性常用 Jaumann，橡胶、生物软组织用超弹性格式更稳健。

在 Abaqus 中如何落地这些度量

理论量最终要落到关键字与输出变量上。开启几何非线性只需在 Step 中设置 NLGEOM=ON（Standard）或使用 Explicit 默认的大变形格式。此时 Abaqus 内部按材料类型自动选择度量：

• 金属塑性：采用 Green-Lagrange 应变与第二类 P-K 应力，配合 Jaumann 客观率，输出变量 NE、LE、S、PE 等。
• 超弹性（橡胶）：以应变能函数 W(C) 直接定义，输出 LE（对数应变）、S、U、V 等。
• 大应变塑性：可切换为对数应变格式，输出 LE 与 Cauchy 应力。

后处理时务必注意：E 是工程应变（小变形假设），LE 是对数应变，NE 是名义应变。大变形问题应读取 LE 而非 E，否则 λ=2 时读数会从 0.693 变成 1.000，误差 44%。同理，应力应读取 S（Cauchy）而非名义应力。

总结本节逻辑链：变形梯度 F 是源头，它同时携带转动与变形；C=FᵀF 剥离转动得到 Green-Lagrange 应变 E，适合参考构型描述；对数应变 ε_ln=ln(U) 适合当前构型描述且具可加性；应力侧则需用客观率消除转动污染，才能写出正确的本构增量关系。这四件工具共同构成 Abaqus 大变形分析的数学地基，后续几何非线性求解、超弹性本构与稳定性分析都建立在此之上。