应变不变量及其工程应用
在固体力学与连续介质力学中，应变张量 $\boldsymbol{\varepsilon}$ 是描述物体变形状态的核心物理量。然而，应变张量的分量值依赖于坐标系的选取，这给不同方向上的变形比较带来了困难。为了解决这一问题，引入了应变不变量（Strain Invariants）的概念。应变不变量是应变张量的多项式函数，其数值在任意正交坐标系变换下保持恒定。这意味着，无论观察者如何旋转坐标系，这些不变量的值都不发生改变，从而提供了描述材料变形状态的客观、本征特征。
对于三阶对称应变张量 $\boldsymbol{\varepsilon}$，存在三个主要的不变量，通常记为 $I_1, I_2, I_3$：
第一应变不变量 ($I_1$)：等于应变张量的迹（Trace），即主应变之和。
$$ I_1 = \text{tr}(\boldsymbol{\varepsilon}) = \varepsilon_{11} + \varepsilon_{22} + \varepsilon_{33} = \varepsilon_1 + \varepsilon_2 + \varepsilon_3 $$
其中 $\varepsilon_1, \varepsilon_2, \varepsilon_3$ 为主应变。$I_1$ 直接关联于体积变化率，在不可压缩材料中，$I_1 = 0$。第二应变不变量 ($I_2$)：等于应变张量所有二阶主子式的和，或者主应变两两乘积之和。
$$ I_2 = \sum_{i<j} \varepsilon_i \varepsilon_j = \varepsilon_1\varepsilon_2 + \varepsilon_2\varepsilon_3 + \varepsilon_3\varepsilon_1 $$
该不变量与形状变化（畸变）密切相关，常用于构建畸变能密度函数。第三应变不变量 ($I_3$)：等于应变张量的行列式，即主应变的乘积。
$$ I_3 = \det(\boldsymbol{\varepsilon}) = \varepsilon_1 \varepsilon_2 \varepsilon_3 $$
$I_3$ 反映了变形后体积相对于原始体积的变化比例。
基于不变量的本构关系构建
在建立材料的本构方程（Constitutive Equations）时，利用应变不变量可以确保方程满足物体客观性（Objectivity）原则，即本构关系不应依赖于观察者的空间取向。对于各向同性超弹性材料，其应变能密度函数 $W$ 通常被表述为应变不变量的函数：
$$ W = W(I_1, I_2, I_3) $$
通过这种形式，工程师可以简洁地描述复杂材料的非线性弹性行为。例如，在橡胶等超弹性材料中，常用的 Mooney-Rivlin 模型或 Ogden 模型均基于不变量构建。应力张量 $\boldsymbol{\sigma}$ 可以通过对应变能密度函数求导获得：
$$ \boldsymbol{\sigma} = \frac{\partial W}{\partial \boldsymbol{\varepsilon}} = \frac{\partial W}{\partial I_1} \frac{\partial I_1}{\partial \boldsymbol{\varepsilon}} + \frac{\partial W}{\partial I_2} \frac{\partial I_2}{\partial \boldsymbol{\varepsilon}} + \frac{\partial W}{\partial I_3} \frac{\partial I_3}{\partial \boldsymbol{\varepsilon}} $$
这种基于不变量的推导方法不仅简化了数学处理，还保证了物理模型的普适性。
工程应用实例：畸变能与屈服准则
在金属塑性力学中，应变不变量在定义屈服准则和能量准则方面发挥着关键作用。
畸变能密度（Distortion Energy Density）：
在弹性力学中，总应变能可以分解为体积改变能和形状改变能（畸变能）。畸变能密度 $U_d$ 与第二应变不变量 $I_2$ 以及第一不变量 $I_1$ 有关。对于各向同性线弹性材料，畸变能密度可表示为：
$$ U_d = \frac{1}{2} \mu (I_2 - \frac{1}{3}I_1^2) $$
其中 $\mu$ 为剪切模量。这一概念是 von Mises 屈服准则 的理论基础。von Mises 准则认为，当畸变能密度达到某一临界值时，材料开始发生塑性流动。该准则在有限元分析（FEA）软件中广泛应用于预测金属结构的屈服行为。最大剪应变理论：
虽然最大剪应变理论直接基于主应变差，但其本质也与应变张量的特征值（即主应变）相关，而主应变正是通过求解特征方程 $\det(\boldsymbol{\varepsilon} - \lambda \mathbf{I}) = 0$ 得到的，该方程的系数正是由应变不变量 $I_1, I_2, I_3$ 构成的。
计算示例
假设某点处的应变张量在笛卡尔坐标系下为：
$$ \boldsymbol{\varepsilon} = \begin{bmatrix} 0.01 & 0.005 & 0 \ 0.005 & 0.02 & 0 \ 0 & 0 & 0.015 \end{bmatrix} $$
我们可以直接计算其不变量：
- $I_1 = 0.01 + 0.02 + 0.015 = 0.045$
- $I_3 = \det(\boldsymbol{\varepsilon}) = 0.01 \times (0.02 \times 0.015 - 0) - 0.005 \times (0.005 \times 0.015 - 0) + 0 = 0.000003 - 0.000000375 = 0.000002625$
- $I_2$ 可通过二阶主子式计算：
$$ I_2 = \begin{vmatrix} 0.01 & 0.005 \ 0.005 & 0.02 \end{vmatrix} + \begin{vmatrix} 0.01 & 0 \ 0 & 0.015 \end{vmatrix} + \begin{vmatrix} 0.02 & 0 \ 0 & 0.015 \end{vmatrix} $$
$$ I_2 = (0.0002 - 0.000025) + (0.00015) + (0.0003) = 0.000175 + 0.00015 + 0.0003 = 0.000625 $$
这些不变量值可用于后续的本构方程计算或屈服判断，且若坐标系发生旋转，计算结果将保持一致。
总结
应变不变量是连接张量分析与工程应用的桥梁。它们消除了坐标系依赖性，使得本构关系的建立更加严谨和简洁。从超弹性材料的能量函数到金属塑性流动的屈服准则，应变不变量在固体力学理论及数值模拟中占据着核心地位。掌握不变量的计算与物理意义，是深入理解复杂材料变形行为的基础。