📚 课程目录第2章金属塑性力学基础

2. 金属塑性力学基础

金属成形仿真的可信度，取决于材料本构模型是否真实反映金属在塑性阶段的力学行为。Abaqus 的塑性材料模块看似只是一组输入框，但每个输入框背后都对应着连续介质塑性力学中的一个明确概念。本节围绕应力应变张量、屈服准则、硬化规律与流动法则四条主线展开，并说明它们如何一一映射到 Abaqus 的材料定义中。

一、应力与应变张量：描述一点的真实受力状态

在三维连续体中，一点的应力状态由二阶对称张量 σ_ij 描述，共 6 个独立分量（3 个正应力、3 个剪应力）。工程计算中常将其分解为球张量与偏张量：

σ_ij = p·δ_ij + s_ij
p  = (σ_11 + σ_22 + σ_33)/3      （静水压力，只改变体积）
s_ij = σ_ij − p·δ_ij              （偏应力，驱动塑性变形）

金属塑性变形在常压下近似体积不可压，因此屈服与塑性流动主要由偏应力 s_ij 控制，这也是 Abaqus 中金属塑性默认采用不可压塑性流动的原因。应变同样分解为弹性应变 ε^e 与塑性应变 ε^p，总应变 ε = ε^e + ε^p。Abaqus 输出的 PEEQ（等效塑性应变）正是由塑性应变张量按等效应变定义累积得到的标量，用于衡量塑性变形程度。

需要特别强调参照约定：本文及后续章节统一采用拉应力为正、压应力为负的符号约定，Abaqus 中 S11、S22、S33 遵循同一约定。若从其他软件或实验数据导入，必须先确认符号方向一致，否则屈服判断会整体反向。

二、屈服准则：判断材料何时进入塑性

屈服准则给出应力状态进入塑性的临界条件。金属最常用的是 von Mises 准则，其核心是等效应力（Mises 应力）：

σ_eq = sqrt( 0.5·[(σ1−σ2)² + (σ2−σ3)² + (σ3−σ1)²] )

当 σ_eq 达到当前屈服应力 σ_y 时材料屈服。在 Abaqus 中，Plastic 材料模型默认即采用 Mises 屈服面，用户只需输入屈服应力与塑性应变的关系，无需手动指定屈服函数。对于各向异性明显的板材（如深冲用铝板），可改用 Hill 或 Barlat 屈服准则，此时需要在 Mechanical → Plasticity → Hill Plasticity 中额外输入各向异性参数 R₁₁、R₂₂、R₃₃、R₁₂ 等。

三、硬化规律：屈服面如何随塑性变形演化

金属屈服后继续加载，屈服应力会随塑性应变增大而提高，这就是硬化。Abaqus 中最常用的是各向同性硬化，通过一条“屈服应力—塑性应变”曲线定义：

*Plastic
 250.0,   0.000
 320.0,   0.020
 380.0,   0.050
 420.0,   0.100
 450.0,   0.200

表中第一列为真实屈服应力（MPa），第二列为对应的塑性应变（非总应变）。注意 Abaqus 要求输入的是塑性应变，若实验给出的是工程应力—应变曲线，需先转换为真实应力—真实应变，再扣除弹性部分。对于循环加载或包辛格效应明显的工况，应改用随动硬化（Kinematic Hardening）或混合硬化，并输入 Cyclic Hardening 参数。

四、流动法则：塑性应变增量的方向

流动法则规定塑性应变增量的方向。关联流动法则假设塑性势函数与屈服函数相同，塑性应变增量沿屈服面外法线方向：

dε^p_ij = dλ · ∂f/∂σ_ij

对 Mises 屈服面，这意味着塑性流动方向与偏应力方向一致，从而保证体积不可压。Abaqus 金属塑性默认采用关联流动法则，用户无需额外设置；只有在岩土、混凝土等摩擦类材料中才需要引入非关联流动。理解这一点有助于解释仿真中常见的“塑性变形体积不变”现象。

五、从实验曲线到 Abaqus 输入：一次完整的参数准备

假设某低碳钢单轴拉伸实验给出工程应力—应变数据：屈服点 250 MPa，工程应变 0.02 时应力 320 MPa，0.05 时 380 MPa。要填入 *Plastic，需先转换：

真实应力 σ_true = σ_eng · (1 + ε_eng)
真实应变 ε_true = ln(1 + ε_eng)
塑性应变 ε_p = ε_true − σ_true / E

取 E = 210000 MPa，对 320 MPa、ε_eng=0.02 这一点：σ_true ≈ 326.4 MPa，ε_true ≈ 0.0198，ε_p ≈ 0.0198 − 326.4/210000 ≈ 0.0182。按同样方法处理其余数据点，即可得到 Abaqus 可用的屈服应力—塑性应变表。这一步是初学者最容易出错的地方：直接把工程应力应变填入 *Plastic，会导致屈服过早、硬化偏软，仿真结果与实验明显偏离。

在 Abaqus/CAE 中，可通过 Material → Evaluate 快速检查输入曲线的斜率是否合理；若塑性段斜率出现负值或剧烈跳变，通常说明数据点顺序或单位有误。

六、常见误区与检查要点

单位不统一：屈服应力用 MPa、弹性模量用 Pa，会导致塑性应变计算完全错误。建议统一采用 mm-N-MPa 或 m-N-Pa 体系。
忽略弹性段扣除：未扣除弹性应变的“塑性应变”偏大，材料显得比实际更硬。
硬化曲线外推：当仿真塑性应变超过输入表最后一点时，Abaqus 默认保持最后一点的屈服应力不变（理想塑性）。若实际工况应变很大，应补充外推数据或改用硬化公式。
屈服准则选择：各向同性材料用 Mises 即可；板材深冲、管材成形等各向异性明显的场景，应通过拉伸试验测定 R 值并选用 Hill 或 Barlat。

把上述四条主线串起来看：应力应变张量提供描述语言，屈服准则给出进入塑性的判据，硬化规律描述屈服面演化，流动法则决定塑性应变方向。Abaqus 的 Elastic、Plastic、Hill Plasticity、Cyclic Hardening 等选项正是这些概念的工程化封装。理解映射关系，才能在材料定义阶段做出有依据的选择，而不是凭经验试参数。