大变形下的有限应变理论简介

在大变形条件下，传统的线性小应变理论已无法准确描述材料的几何非线性和物理非线性。有限应变理论（Finite Strain Theory）提供了一套完整的数学框架，用以刻画任意尺度的变形过程。本文从基本概念出发，系统阐述有限应变的几何描述、常用应变张量、应力测度以及典型本构模型，并给出简明的数值示例，帮助读者快速建立起大变形分析的理论基础。

• 参考构型与当前构型

  • 参考构型（Reference configuration）记为 (\mathcal{B}_0)，用坐标 (\mathbf{X}) 表示。
  • 当前构型（Current configuration）记为 (\mathcal{B})，用坐标 (\mathbf{x}) 表示。
  • 两者之间的映射为 位移映射 (\boldsymbol{\chi} : \mathcal{B}_0 \rightarrow \mathcal{B})，满足 (\mathbf{x}= \boldsymbol{\chi}(\mathbf{X}))。

• 变形梯度
  [
  \mathbf{F}= \frac{\partial \mathbf{x}}{\partial \mathbf{X}} = \nabla_{!0}\boldsymbol{\chi}
  ]
  (\mathbf{F}) 同时包含了旋转、拉伸和剪切信息，是有限应变理论的核心变量。

• 不可压缩条件
  [
  \det \mathbf{F}=J=1 \quad (\text{体积保持})
  ]

2. 常用应变张量

注：在实际计算中，往往选取 (\mathbf{E}) 或 (\boldsymbol{\varepsilon}) 作为应变测度，因为它们在 (\mathbf{F}) 接近单位张量时退化为线性小应变。

3. 应力的测度

• 第一 Piola–Kirchhoff 应力 (\mathbf{P})
  [
  \mathbf{P}= J,\boldsymbol{\sigma},\mathbf{F}^{-\mathrm{T}}
  ]
  与参考构型面积元素对应，常用于弱形式的有限元实现。

• 第二 Piola–Kirchhoff 应力 (\mathbf{S})
  [
  \mathbf{S}= \mathbf{F}^{-1}\mathbf{P}= \mathbf{F}^{-1}\boldsymbol{\sigma}\mathbf{F}^{-\mathrm{T}}J
  ]
  与 (\mathbf{E}) 配对使用，便于构造能量函数。

• Cauchy（真）应力 (\boldsymbol{\sigma})
  [
  \boldsymbol{\sigma}= \frac{1}{J}\mathbf{P}\mathbf{F}^{\mathrm{T}}
  ]
  与当前构型面积元素对应，直观反映实际受力状态。

4. 本构关系的表述

在有限应变框架下，本构模型通常以 应变能密度函数 (W(\mathbf{C})) 或 (W(\mathbf{F})) 为基础：

[
\mathbf{S}= 2\frac{\partial W}{\partial \mathbf{C}}, \qquad
\mathbf{P}= \frac{\partial W}{\partial \mathbf{F}}
]

4.1 超弹性模型（Neo‑Hookean）

[
W = \frac{\mu}{2}\left(\mathrm{tr},\mathbf{C}-3\right) - \mu\ln J + \frac{\lambda}{2}(\ln J)^2
]

• (\mu,\lambda) 为拉梅常数。
• 对应的第一 Piola–Kirchhoff 应力
  [
  \mathbf{P}= \mu\left(\mathbf{F}-\mathbf{F}^{-\mathrm{T}}\right) + \lambda \ln J,\mathbf{F}^{-\mathrm{T}}
  ]

4.2 Mooney–Rivlin 模型

[
W = C_{10}(\bar I_1-3) + C_{01}(\bar I_2-3) + \frac{1}{D_1}(J-1)^2
]
其中 (\bar I_1,\bar I_2) 为去体积化的不变量。

4.3 大变形塑性（J2 流动）

• 等效应变能
  [
  W = \int_0^{\varepsilon_e} \sigma_y(\varepsilon), d\varepsilon
  ]
• 流动法则（基于对数应变）
  [
  \dot{\boldsymbol{\varepsilon}} = \dot{\gamma},\frac{\partial f}{\partial \boldsymbol{\sigma}}
  ]
  其中 (f) 为屈服函数，(\dot{\gamma}) 为塑性乘子。

5. 数值示例：单轴拉伸

假设一根弹性杆在参考长度 (L_0=1) m 上施加 20% 拉伸，即终长度 (L=1.2) m。取材料为 Neo‑Hookean，(\mu=0.8) MPa，(\lambda=1.2) MPa。

1. 变形梯度（一维情形）
   [
   \mathbf{F}= \begin{bmatrix} \lambda_z & 0 & 0\0&1&0\0&0&1 \end{bmatrix},
   \quad \lambda_z = \frac{L}{L_0}=1.2
   ]
2. 体积比 (J = \det\mathbf{F}=1.2)。
3. 第一 Piola–Kirchhoff 应力
   [
   P_{zz}= \mu(\lambda_z-\lambda_z^{-1}) + \lambda\ln J,\lambda_z^{-1}
   ]
   代入数值
   [
   P_{zz}=0.8(1.2-0.8333)+1.2\ln(1.2)\times0.8333\approx0.8(0.3667)+1.2(0.1823)\times0.8333\approx0.293+0.182\approx0.475;\text{MPa}
   ]
4. Cauchy 应力
   [
   \sigma_{zz}= \frac{P_{zz}}{J}= \frac{0.475}{1.2}\approx0.396;\text{MPa}
   ]

该结果表明，在 20% 拉伸下，材料的真实应力约为 0.4 MPa，远低于线性弹性预测的 (\sigma=E\varepsilon)（若取 (E=2\mu(1+\nu))），体现了几何非线性的软化效应。

6. 实现要点与常见误区

• 保持张量对称性：在数值积分时，(\mathbf{C}) 与 (\mathbf{E}) 必须保持对称，否则会导致能量不守恒。
• 避免除零：在计算 (\mathbf{F}^{-1}) 或 (\mathbf{C}^{-1}) 时，若出现局部压缩导致 (J) 接近零，需要使用正则化技术（如添加小的正数 (\epsilon)）。
• 应变测度的选取：对数应变在大塑性流动中更具物理意义，但其数值实现稍显复杂；绿拉格朗日应变实现最为直接。
• 本构模型的可微性：在有限元求解非线性方程时，需要提供应力张量对 (\mathbf{F}) 的一致切线模量（consistent tangent），否则收敛速度会显著下降。

7. 小结

• 大变形下的有限应变理论以 变形梯度 (\mathbf{F}) 为基本变量，构建了右/左 Cauchy–Green 张量、Green‑Lagrange 应变以及对应的 Piola‑Kirchhoff 应力体系。
• 通过应变能密度函数 (W) 可以统一描述弹性、超弹性以及塑性材料的本构行为。
• 示例演示了在单轴拉伸中从位移到应力的完整计算链条，验证了理论的可操作性。
• 实际工程实现时需关注张量对称、数值稳定性以及一致切线的提供，以保证大变形数值模拟的准确性与收敛性。

掌握上述核心概念后，读者即可在有限元软件或自编代码中实现大变形分析，为后续的非线性结构、软体机器人以及生物组织力学等前沿研究奠定坚实的理论基础。