应力不变量及其工程应用
在固体力学与材料科学领域，应力张量 $\boldsymbol{\sigma}$ 是描述物体内部某一点受力状态的核心数学工具。然而，应力张量的分量值依赖于坐标系的选取，具有方向依赖性。为了客观地描述材料的力学行为，特别是判断材料是否屈服或破坏，我们需要寻找一组与坐标系无关、仅由应力状态本身决定的标量量，即应力不变量（Stress Invariants）。
应力不变量共有三个，通常记为 $I_1$、$I_2$ 和 $I_3$。它们分别对应应力张量特征值（即主应力 $\sigma_1, \sigma_2, \sigma_3$）的初等对称多项式。理解这些不变量对于建立本构方程、分析复杂应力状态下的材料响应至关重要。
应力不变量的数学定义
对于任意一点处的应力张量，其三个不变量定义如下：
第一应力不变量 ($I_1$)：
$$ I_1 = \text{tr}(\boldsymbol{\sigma}) = \sigma_{xx} + \sigma_{yy} + \sigma_{zz} = \sigma_1 + \sigma_2 + \sigma_3 $$
它等于应力张量的迹，也等于三个主应力之和。在体积变化分析中，$I_1$ 与体积应变直接相关。第二应力不变量 ($I_2$)：
$$ I_2 = \frac{1}{2}[(\text{tr}\boldsymbol{\sigma})^2 - \text{tr}(\boldsymbol{\sigma}^2)] $$
或者用主应力表示为：
$$ I_2 = \sigma_1\sigma_2 + \sigma_2\sigma_3 + \sigma_3\sigma_1 $$
它反映了应力状态中剪切分量的贡献，是许多屈服准则（如 von Mises 准则）的核心参数。第三应力不变量 ($I_3$)：
$$ I_3 = \det(\boldsymbol{\sigma}) = \sigma_1\sigma_2\sigma_3 $$
它等于应力张量的行列式，即三个主应力的乘积。
工程应用与屈服准则
应力不变量在工程实践中最广泛的应用是构建屈服准则（Yield Criteria）。屈服准则用于判断材料从弹性阶段进入塑性阶段的临界状态。由于屈服通常与体积变化无关（对于大多数金属塑性变形），因此基于偏应力张量（Deviatoric Stress Tensor）的不变量更为常用。
1. von Mises 屈服准则
这是工程中最常用的屈服准则，适用于延性金属。它基于第二偏应力不变量 $J_2$（由 $I_1$ 和 $I_2$ 导出）：
$$ J_2 = \frac{1}{2} \boldsymbol{s} : \boldsymbol{s} $$
其中 $\boldsymbol{s}$ 是偏应力张量。当 $J_2$ 达到临界值时，材料发生屈服。该准则假设屈服主要取决于剪切能，与静水压力无关。
2. Tresca 屈服准则
Tresca 准则基于最大剪应力理论。虽然它直接涉及主应力差值，但其本质也可通过不变量来理解。它认为当最大剪应力达到单向拉伸试验中的屈服剪应力时，材料屈服。
3. 考虑静水压力影响的准则
对于脆性材料或岩土材料，静水压力（由 $I_1$ 表征）对强度有显著影响。例如，Mohr-Coulomb 准则就同时考虑了 $I_1$ 和 $J_2$，能够更准确地描述土体或岩石在高压下的破坏行为。
计算示例
假设某点处的应力状态由以下主应力给出（单位：MPa）：
$$ \sigma_1 = 100, \quad \sigma_2 = 50, \quad \sigma_3 = -20 $$
我们可以计算其应力不变量：
计算 $I_1$：
$$ I_1 = 100 + 50 + (-20) = 130 , \text{MPa} $$计算 $I_2$：
$$ I_2 = (100)(50) + (50)(-20) + (-20)(100) $$
$$ I_2 = 5000 - 1000 - 2000 = 2000 , \text{MPa}^2 $$计算 $I_3$：
$$ I_3 = 100 \times 50 \times (-20) = -100,000 , \text{MPa}^3 $$
若该材料为钢，屈服强度 $\sigma_y = 250 , \text{MPa}$，使用 von Mises 准则判断是否屈服：
等效应力（von Mises stress）$\sigma_{vm}$ 与 $J_2$ 的关系为 $\sigma_{vm} = \sqrt{3J_2}$。
首先计算偏应力不变量 $J_2$：
$$ J_2 = I_2 - \frac{1}{3}I_1^2 = 2000 - \frac{1}{3}(130)^2 = 2000 - 5633.33 < 0 $$
注意：此处直接代入公式需小心，通常 $J_2$ 定义为偏应力张量的第二不变量。更直观的方法是使用主应力公式计算 von Mises 应力：
$$ \sigma_{vm} = \frac{1}{\sqrt{2}} \sqrt{(\sigma_1-\sigma_2)^2 + (\sigma_2-\sigma_3)^2 + (\sigma_3-\sigma_1)^2} $$
$$ \sigma_{vm} = \frac{1}{\sqrt{2}} \sqrt{(50)^2 + (70)^2 + (-120)^2} = \frac{1}{\sqrt{2}} \sqrt{2500 + 4900 + 14400} $$
$$ \sigma_{vm} = \frac{1}{\sqrt{2}} \sqrt{21800} \approx 104.4 , \text{MPa} $$
由于 $\sigma_{vm} (104.4 , \text{MPa}) < \sigma_y (250 , \text{MPa})$，该点处于弹性状态，未发生屈服。
总结
应力不变量提供了独立于坐标系的应力状态描述方法，是连接微观力学机制与宏观工程应用的桥梁。通过 $I_1$、$I_2$ 和 $I_3$，工程师可以精确量化体积效应和剪切效应，从而选择合适的本构模型和屈服准则，确保结构设计与安全评估的准确性。在有限元分析软件中，这些不变量也是后处理输出和材料失效预测的基础数据。