📚 课程目录第4章金属材料本构模型在Abaqus中的实现

4. 金属材料本构模型在Abaqus中的实现

在 Abaqus 中，金属材料的塑性行为并不是通过一个“材料类型”开关来选择的，而是由一组相互配合的选项共同定义：弹性段、塑性硬化段、率相关性、硬化方向（各向同性或随动）以及必要的损伤或失效判据。理解这些选项之间的层级关系，是正确输入材料参数的前提。本节以单轴拉伸试验曲线为起点，说明如何把一条工程应力—应变曲线转换为 Abaqus 可用的材料参数，并给出可直接套用的输入结构。

一、本构模型的层级结构

Abaqus/Standard 与 Abaqus/Explicit 中，金属塑性最常用的组合是：*ELASTIC 定义弹性模量与泊松比，*PLASTIC 定义屈服面与硬化数据，*RATE DEPENDENT 定义应变率敏感性，*KINEMATIC HARDENING 或 *COMBINED HARDENING 定义随动/混合硬化。这些关键字在输入文件中按顺序出现，但在物理上它们描述的是同一套本构关系的不同侧面。下图给出各选项与物理含义之间的对应关系，以及单轴拉伸数据进入模型的路径。

二、从工程曲线到真应力—真应变

试验机给出的是工程应力 σ_eng = F/A₀ 与工程应变 ε_eng = ΔL/L₀。Abaqus 的 *PLASTIC 需要的是真应力与塑性应变。转换关系为：

σ_true = σ_eng * (1 + ε_eng)
ε_true = ln(1 + ε_eng)
ε_pl   = ε_true - σ_true / E

以某低碳钢为例（示例数据，非实测）：E = 210000 MPa，ν = 0.3，屈服点 σ_eng = 250 MPa，对应 ε_eng = 0.00119。则 σ_true ≈ 250.3 MPa，ε_true ≈ 0.00119，ε_pl ≈ 0.00119 − 250.3/210000 ≈ 0.0000。继续取颈缩前若干点，例如 σ_eng = 400 MPa、ε_eng = 0.10，则 σ_true = 440 MPa，ε_true = 0.0953，ε_pl = 0.0953 − 440/210000 ≈ 0.0932。把这些点按“屈服应力—塑性应变”成对写入 *PLASTIC，第一行塑性应变必须为 0，对应初始屈服。

颈缩之后工程应力下降，但真应力仍在上升。Abaqus 的 *PLASTIC 数据要求塑性应变单调递增，因此通常只取到颈缩点，颈缩后的行为由损伤或失效模型处理，而不是继续外推。

三、硬化方向的选择

各向同性硬化假设屈服面中心固定、半径扩大，适合单调加载；随动硬化假设屈服面尺寸不变、中心随塑性应变平移，能描述包辛格效应，适合循环或反向加载。Abaqus 中：

*PLASTIC, HARDENING=ISOTROPIC：默认，输入一组 σ_y–ε_pl 数据。
*PLASTIC, HARDENING=KINEMATIC：配合 *KINEMATIC HARDENING，可用双线性或 Armstrong–Frederick 参数。
*PLASTIC, HARDENING=COMBINED：同时给出各向同性分量与随动分量，适合既有强化又有循环软化的场景。

选择依据是加载路径：若结构只经历单次成形，各向同性硬化足够；若经历拉—压反复，必须用随动或混合硬化，否则反向屈服会被高估。

四、率相关性的输入

当应变率对屈服应力有明显影响时，在 *PLASTIC 之后加 *RATE DEPENDENT。Abaqus 提供多种形式，最常用的是 Cowper–Symonds 幂律：

*RATE DEPENDENT, TYPE=COWPER-SYMONDS
  D, n

其中屈服应力放大系数为 1 + (ε̇/D)^(1/n)。例如 D = 40 s⁻¹、n = 5，表示应变率 40 s⁻¹ 时屈服应力约提高 1 倍。也可用 TYPE=YIELD RATIO 直接给出不同应变率下的屈服应力比表，适合有动态试验数据的场合。

五、可直接套用的输入模板

下面给出一个完整的材料定义片段，包含弹性、各向同性塑性、率相关与随动硬化切换注释，可直接粘贴到 .inp 中修改：

*MATERIAL, NAME=STEEL_Q235
*ELASTIC
 210000., 0.3
*PLASTIC, HARDENING=ISOTROPIC
 250.0, 0.0000
 290.0, 0.0100
 330.0, 0.0300
 370.0, 0.0600
 400.0, 0.1000
*RATE DEPENDENT, TYPE=COWPER-SYMONDS
 40.0, 5.0
** 若需随动硬化，将上一行 *PLASTIC 改为：
** *PLASTIC, HARDENING=KINEMATIC
** *KINEMATIC HARDENING, TYPE=ARMSTRONG-FREDERICK
**  C1, C2

输入 *PLASTIC 数据时，塑性应变列必须从 0 开始并严格递增；若出现负值或回退，Abaqus 会在数据检查阶段报错。真应力必须为正，且第一行屈服应力应等于弹性段结束时的真应力。

六、结果判断与常见问题

提交计算后，可在 .dat 文件中检查材料数据是否被正确读取；在 .msg 中确认没有“plastic strain not increasing”类警告。后处理时提取 PEEQ（等效塑性应变）与 MISES，若 PEEQ 在屈服前就非零，说明弹性模量或屈服应力输入有误；若反向加载时屈服应力与正向相同，说明误用了各向同性硬化。对于率相关模型，可对比不同加载速度下的力—位移曲线，确认屈服平台随速度提高而抬升，与 *RATE DEPENDENT 参数预期一致。

把试验曲线转成材料参数的过程，本质上是把“工程量”翻译成“真应力量”，再按 Abaqus 的关键字结构组织成弹性、塑性、率相关、硬化方向四个层次。只要转换公式正确、数据单调、硬化方向与加载路径匹配，单轴拉伸曲线就能可靠地驱动后续的成形仿真。