塑性流动法则与硬化模型

在固体力学研究中，描述材料在超过弹性极限后发生永久变形的过程，即塑性变形，是构建本构模型的核心任务。一个完整的塑性本构模型通常由三个基本要素组成：屈服准则（Yield Criterion）、塑性流动法则（Plastic Flow Rule）

在固体力学研究中，描述材料在超过弹性极限后发生永久变形的过程，即塑性变形，是构建本构模型的核心任务。一个完整的塑性本构模型通常由三个基本要素组成：屈服准则（Yield Criterion）、塑性流动法则（Plastic Flow Rule）以及硬化模型（Hardening Model）。

其中，屈服准则定义了材料进入塑性状态的边界，而流动法则决定了塑性应变增量的方向，硬化模型则描述了屈服面随塑性变形演化的规律。本文将重点深入探讨塑性流动法则与硬化模型的核心理论。
当材料的应力状态达到或超过屈服面时，材料开始发生塑性变形。塑性流动法则的任务是确定塑性应变率 $\dot{\varepsilon}{ij}^p$ 与当前应力状态 $\sigma{ij}$ 之间的数学关系。

1. 正交性原理与关联流则

在经典的塑性力学理论中，最常用的假设是正交性原理（Normality Rule）。该原理认为，塑性应变增量向量 $\Delta \varepsilon_{ij}^p$ 的方向必须垂直于当前的屈服面。

如果我们将屈服函数表示为 $f(\sigma_{ij}, \kappa) = 0$（其中 $\kappa$ 为硬化变量），那么根据正交性原理，塑性应变增量可以表示为：
$$\dot{\varepsilon}{ij}^p = \dot{\lambda} \frac{\partial f}{\partial \sigma{ij}}$$
其中，$\dot{\lambda} \ge 0$ 是塑性乘子，其大小取决于加载速率和材料的硬化特性。

在这种情况下，屈服函数 $f$ 本身也可以作为塑性势函数（Plastic Potential Function） $g$。当 $f = g$ 时，这种流动法则被称为关联流则（Associated Flow Rule）。对于许多金属材料（如遵循 Von Mises 准则的材料），关联流则能够很好地预测其塑性变形行为。

2. 非关联流则

然而，对于某些特定材料，如砂土、岩石等颗粒材料，关联流则会导致预测结果与实验观察严重不符（例如，关联流则往往会过度预测体积膨胀效应）。为了修正这一问题，引入了非关联流则（Non-associated Flow Rule）。

在非关联流则中，塑性势函数 $g$ 与屈服函数 $f$ 是不同的：
$$\dot{\varepsilon}{ij}^p = \dot{\lambda} \frac{\partial g}{\partial \sigma{ij}}, \quad f \neq g$$
通过独立定义 $g$，研究人员可以根据实验数据调整塑性应变的方向，从而更准确地模拟材料的剪胀（Dilatancy）或剪缩行为。

硬化模型

随着塑性变形的进行，材料的屈服边界会发生改变，这种现象称为硬化。硬化模型描述了屈服面在应力空间中是如何演化的。根据屈服面演化的几何特征，硬化模型主要分为以下三类：

1. 各向同性硬化 (Isotropic Hardening)

各向同性硬化假设屈服面在应力空间中是均匀地扩张或收缩的，而其中心位置保持不变。

  • 特征：屈服面的大小随等效塑性应变 $\bar{\varepsilon}^p$ 的增加而增大。
  • 数学表达：若初始屈服面为 $f(\sigma_{ij}) = 0$，则硬化后的屈服面为 $f(\sigma_{ij}, \bar{\varepsilon}^p) = 0$。
  • 适用场景：适用于单向拉伸或压缩加载过程，能够较好地描述材料在持续单向加载下的强度提升。

2. 运动学硬化 (Kinematic Hardening)

在循环加载（如拉-压交替加载）过程中，材料表现出一种特殊的现象——包辛格效应（Bauschinger Effect）：即在单向拉伸产生塑性变形后，材料在反向压缩时的屈服强度会显著降低。各向同性硬化模型无法解释这一现象，因此需要引入运动学硬化。

运动学硬化假设屈服面在应力空间中发生整体平移，而其形状和大小保持不变。

  • 特征：引入了一个称为“背应力”（Back Stress）的变量 $\alpha_{ij}$。新的屈服函数表示为：
    $$f(\sigma_{ij} - \alpha_{ij}) = 0$$
  • 物理意义：$\alpha_{ij}$ 代表了由于位错堆积等微观机制引起的内部应力场，它改变了应力空间中屈服面的中心位置。

3. 混合硬化 (Combined Hardening)

在实际工程应用中，金属材料往往同时表现出各向同性硬化和运动学硬化的特征。为了获得更精确的本构描述，通常采用混合硬化模型：

  • 各向同性部分：负责描述屈服面半径的扩张（强度提升）。
  • 运动学部分：负责描述屈服面的平移（捕捉包辛格效应）。

通过将两者结合，模型可以同时模拟单向加载下的强度增加以及循环加载下的滞回环（Hysteresis Loop）行为。

总结与应用示例

为了直观理解，我们可以通过一个简单的例子来总结：

示例：金属棒的循环拉压实验

  1. 初始阶段：金属棒在应力空间中有一个初始的 Von Mises 屈服圆。
  2. 单向拉伸：当应力超过屈服点，材料发生塑性变形。此时，各向同性硬化使圆圈变大，流动法则决定了应变沿径向方向增加。
  3. 反向压缩：当应力转为压缩时，如果仅使用各向同性硬化，预测的压缩屈服点会非常高；但由于实际存在运动学硬化，屈服圆向负应力方向平移，导致在较低的压缩应力下就发生了塑性变形，这完美契合了实验观察到的包辛格效应。

在有限元分析（FEA）软件（如 Abaqus 或 ANSYS）中，用户可以通过选择不同的硬化模型（如 Isotropic, Kinematic, 或 Combined）来模拟复杂的金属成形或疲劳问题。掌握这些理论对于建立准确的数值模型、预测结构失效具有至关重要的意义。