3. 金属塑性理论基础与Abaqus实现

金属材料在超过屈服点后，应力与应变不再保持线性关系，卸载后还会留下永久变形。要在Abaqus中把这种行为算准，必须同时回答三个问题：材料在什么应力组合下开始屈服（屈服准则）、屈服面之后如何变化（硬化法则）、塑性应变沿什么方向增长（流动法则）。这三者共同构成塑性本构模型的骨架，缺一不可。本节围绕这三块内容展开，并给出在Abaqus中定义弹塑性材料、验证单轴曲线的完整做法。

一、屈服准则：从单向应力到复杂应力状态

单轴拉伸时，屈服点就是应力达到屈服强度。但产品实际受力往往是多向的，需要用一个标量函数来判断是否屈服。对大多数延性金属（钢、铝、铜），静水压力对屈服影响很小，因此采用不受静水压力影响的准则。

von Mises准则是最常用的各向同性屈服准则，它认为当偏应力第二不变量J₂达到临界值时材料屈服，等价于等效应力σ_e达到屈服强度：

σ_e = sqrt( 0.5*[(σ1-σ2)^2 + (σ2-σ3)^2 + (σ3-σ1)^2] )
屈服条件： σ_e = σ_y

在Abaqus中，*PLASTIC配合各向同性硬化默认就使用von Mises。它适合各向同性、无明显织构的金属。

Hill准则用于各向异性材料，典型对象是轧制钢板、拉拔棒材、挤压型材。它把von Mises推广为带方向权重的二次型：

f = F(σ22-σ33)^2 + G(σ33-σ11)^2 + H(σ11-σ22)^2
    + 2L*σ23^2 + 2M*σ31^2 + 2N*σ12^2
屈服条件： f = 1

其中F、G、H、L、M、N由不同方向的屈服应力比值确定。若各方向屈服强度相等，Hill退化为von Mises。Abaqus中通过*POTENTIAL定义Hill系数，需先确定材料坐标系（通常与轧制方向对齐）。

二、硬化法则：屈服面如何随塑性变形演化

屈服之后继续加载，屈服面会变化，这就是硬化。三种典型法则的差别在于屈服面是“均匀膨胀”“整体平移”还是“既膨胀又平移”。

• 各向同性硬化：屈服面中心不动、半径增大。适合单调加载、比例加载，如一次性压溃、成形分析。Abaqus中用*PLASTIC给出真实应力—塑性应变数据点即可。
• 随动硬化：屈服面大小不变、中心平移（背应力α）。能描述Bauschinger效应，适合循环加载、低周疲劳、反向加载。Abaqus中用*KINEMATIC（线性）或*COMBINED配合非线性随动参数。
• 混合硬化：屈服面同时膨胀和平移，用参数控制各向同性部分与随动部分的比例，最贴近真实循环行为。Abaqus中*COMBINED的R参数即各向同性权重。

选择依据很直接：加载是否反向、是否循环。单调加载用各向同性硬化足够；出现反向屈服且需捕捉包辛格效应时，必须用随动或混合硬化。

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

流动法则规定塑性应变增量dε^p的方向。Abaqus金属塑性默认采用关联流动法则，即塑性势函数与屈服函数相同，塑性应变增量沿屈服面外法线方向：

dε^p = dλ * ∂f/∂σ

对von Mises，这意味着塑性变形不产生体积变化（不可压），且方向由偏应力决定。关联流动是金属塑性的标准假设，一般不需要用户修改；只有在岩土、粉末等摩擦材料中才考虑非关联流动。

四、在Abaqus中定义弹塑性材料

以Q345钢为例，弹性模量210000 MPa，泊松比0.3，屈服强度345 MPa。塑性段必须输入真实应力—塑性应变，而不是工程应力—应变。转换关系为：

真实应力 σ_true = σ_eng * (1 + ε_eng)
真实应变 ε_true = ln(1 + ε_eng)
塑性应变 ε_pl = ε_true - σ_true / E

假设单轴试验得到工程数据，转换后输入Abaqus的*PLASTIC如下（示例数据，仅用于演示格式）：

*MATERIAL, NAME=Q345
*ELASTIC
210000., 0.3
*PLASTIC
345., 0.0
420., 0.02
480., 0.05
520., 0.10
550., 0.20

第一列是屈服应力，第二列是塑性应变，第一行塑性应变必须为0。若使用随动硬化，则改用*KINEMATIC并给出背应力—塑性应变数据；混合硬化用*COMBINED并指定R值。

五、单轴曲线验证

建模时取一个单元或小试件，一端固定、一端施加位移，输出反力与位移，换算成应力—应变后与输入曲线对比。判断要点：弹性段斜率应等于E；屈服点应出现在345 MPa附近；塑性段应与输入数据点吻合；卸载后应保留塑性应变。若曲线在屈服点出现明显折角而非平滑过渡，说明输入数据点太少或未做真实应力转换。若反向加载时没有包辛格效应，说明硬化法则选成了各向同性。