2. 金属塑性力学基础

在 Abaqus 中做金属成形仿真,真正决定结果可信度的并不是网格画得多细,而是本构模型是否与材料真实行为匹配。金属塑性力学的五个核心概念——应力应变度量、屈服准则、硬化规律、流动法则与塑性势——共同构成了 Abaqus 塑性材料定义的骨架。理解它们之间的逻辑关系,才能在材料卡片中填对每一个参数。

一、应力与应变的度量:为什么大变形要用真应力—真应变

在弹性小变形阶段,工程应力 σeng=F/A0 与工程应变 εeng=ΔL/L0 足够使用。但金属成形中试样截面持续收缩,若仍用初始面积 A0 计算,应力会被系统性低估。因此塑性阶段必须切换到真应力与真应变:

σ_true = F / A          (A 为当前截面积)
ε_true = ln(L / L0)     (对数应变,可累加)
ε_true = ln(1 + ε_eng)
σ_true = σ_eng · (1 + ε_eng)   (颈缩前近似)

举例:某低碳钢拉伸到工程应变 0.20 时工程应力为 420 MPa,则真应变约为 ln(1.20)=0.182,真应力约为 420×1.20=504 MPa。Abaqus 的 *PLASTIC 卡片要求输入的正是真应力—塑性应变对,其中塑性应变 εpl=εtrue−σtrue/E。若误填工程值,屈服后应力会被低估约 15%~20%,成形载荷预测随之偏低。

此外,Abaqus/Standard 默认使用小应变理论,当成形中累积等效塑性应变超过约 5% 或局部转动较大时,应开启 *STEP, NLGEOM=ON,让程序采用有限应变公式(基于 Jaumann 或 Green-Naghdi 客观应力率),否则会出现虚假的应力振荡。

二、屈服准则:从 Tresca 到 von Mises

屈服准则回答“什么应力状态下材料开始塑性流动”。对延性金属,最常用的是 von Mises 准则,其等效应力为:

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

当 σvm 达到单轴屈服强度 σy 时进入塑性。它在主应力空间中是一个圆柱面,与静水压力无关,因此适合描述大多数金属。Tresca 准则以最大剪应力为判据,结果偏保守(屈服面为内接六棱柱),在 Abaqus 中较少直接使用,但可用于校核。对于各向异性明显的轧制板材,则应改用 Hill48 或 Barlat 准则,通过 *POTENTIAL 与 *ANISOTROPIC YIELD 定义 L、M、N 等参数。

三、硬化规律:等向硬化与随动硬化

屈服面随塑性变形如何演化,由硬化规律决定。等向硬化(Isotropic Hardening)让屈服面均匀膨胀,适合单调加载的成形分析;随动硬化(Kinematic Hardening)让屈服面平移,用于描述包辛格效应,适合循环加载或回弹反复路径。Abaqus 中:

  • *PLASTIC, HARDENING=ISOTROPIC:输入真应力—塑性应变表,程序线性插值。
  • *PLASTIC, HARDENING=KINEMATIC:需配合 *CYCLIC HARDENING 给出背应力演化。
  • *PLASTIC, HARDENING=COMBINED:等向与随动按比例混合,参数由循环试验拟合。

硬化曲线的斜率 H=dσ/dεpl 直接决定颈缩后的承载能力。若曲线在较大应变处被截断,Abaqus 会保持最后一段斜率外推,可能造成载荷虚高,因此建议将试验数据外推到仿真中可能出现的最大等效塑性应变。

四、流动法则与塑性势:塑性应变往哪个方向走

流动法则规定塑性应变增量的方向。关联流动(Associated Flow)假设塑性势函数 g 等于屈服函数 f,于是:

dε_pl = dλ · ∂g/∂σ = dλ · ∂f/∂σ

对 von Mises 屈服面,∂f/∂σ 的方向正是偏应力方向,因此塑性变形体积不可压缩(dεvolpl=0),这与金属实验一致。Abaqus 的金属塑性默认采用关联流动,用户无需额外指定。若采用非关联流动(g≠f),则需通过 *POTENTIAL 单独定义,常见于岩土或粉末材料,金属成形中一般不用。

在 Abaqus/Explicit 中,塑性流动通过径向返回算法在每个增量步内迭代求解,增量步过大时可能不收敛或产生应力漂移,建议控制质量缩放并监控 ALLAE 与 ALLKE 的比值。

五、五要素的协同关系

这五个概念并非孤立:应力应变度量提供输入数据,屈服准则定义“何时开始塑性”,硬化规律定义“屈服面如何变”,流动法则与塑性势定义“塑性应变朝哪个方向、以多大规模发展”。它们共同决定 Abaqus 每个增量步中应力更新与切线刚度的计算路径。

金属塑性本构五要素在 Abaqus 增量步中的协同 应力应变度量 真应力 σ_true 塑性应变 ε_pl 屈服准则 σ_vm ≥ σ_y ? von Mises / Hill48 硬化规律 等向 / 随动 / 混合 σ_y = σ_y(ε_pl) 流动法则 dε_pl = dλ·∂g/∂σ 方向 + 大小 塑性势 g 关联:g = f 非关联:g ≠ f 应力更新 径向返回算法 切线刚度 D_ep ε_pl 累积反馈 实线:本构计算主路径 虚线:塑性应变累积反馈 橙色:屈服面演化

在 Abaqus 材料定义中,这五要素的落地方式可以概括为:*ELASTIC 给出 E 与 ν;*PLASTIC 给出硬化曲线;*POTENTIAL 与各向异性屈服参数给出屈服面形状;流动法则由程序按关联流动自动处理。只要这四类数据与试验曲线一致,仿真才能真实反映金属从屈服、硬化到颈缩、断裂前的力学响应。

常见错误:把工程应力—应变直接填入 *PLASTIC;硬化曲线未外推到足够大的塑性应变;大变形工况未开启 NLGEOM。这三类问题会分别导致屈服偏低、载荷虚高和应力振荡,是金属成形仿真中最典型的偏差来源。