4. 塑性本构模型基础选型

塑性本构模型不是一组可以随意拼装的“选项”,它由三个相互耦合的部分构成:屈服准则决定材料在什么应力状态下从弹性进入塑性;流动法则决定塑性应变增量的方向;硬化规律决定屈服面随塑性变形如何演化。选型的本质,是把这三者与具体金属成形问题的应力状态、应变路径和加载历史匹配起来。三者中任何一项选错,都会在成形载荷、回弹量或残余应力上产生系统性偏差。

一、屈服准则:先看应力三轴度与静水压力敏感性

金属塑性最常用的起点是 von Mises 准则,其屈服函数为

f = sqrt(3/2 · s:s) − σ_y(ε_p) = σ_eq − σ_y(ε_p)
s = σ − (1/3)·tr(σ)·I

它只依赖偏应力第二不变量 J₂,认为静水压力不影响屈服。对绝大多数致密金属(钢、铝、铜)在常规冷成形中,这一假设成立,因此 von Mises 是板材冲压、锻造、挤压的首选。Tresca 准则以最大剪应力为判据,形式更简单,结果偏保守(屈服面内接于 von Mises 六边形),适合解析估算和需要保守设计的场合。

当问题涉及高应力三轴度或孔隙演化时,必须换准则。GTN(Gurson–Tvergaard–Needleman)模型引入孔洞体积分数 f,屈服面随 f 增大而收缩,用于预测韧性断裂;Drucker–Prager 类准则引入静水压力项,适用于粉末冶金坯料、多孔材料或岩石类介质。判断依据可以概括为:若成形过程中孔洞、损伤或静水压力对屈服有可测影响,就不能用 von Mises。

二、流动法则:关联还是非关联

流动法则给出塑性应变增量方向:

dε_p = dλ · ∂g/∂σ

当塑性势 g 等于屈服函数 f 时,称为关联流动法则。对金属,关联流动法则与实验吻合良好,且保证材料稳定性(Drucker 公设),因此是默认选择。它带来一个重要推论:塑性应变增量垂直于屈服面,且塑性变形不产生体积变化(不可压)。在有限元中,这直接要求使用不可压或近似不可压单元,否则会出现体积锁死。

非关联流动(g ≠ f)主要用于摩擦材料、土体和某些聚合物,金属成形中很少需要。若在金属问题中强行使用非关联法则,会破坏本构的对称性,增加计算成本且缺乏实验支撑。因此本节的原则是:金属成形默认关联流动,除非有明确的实验证据表明塑性势与屈服面不一致。

三、硬化规律:等向、随动还是混合

硬化规律描述屈服面尺寸和位置的变化。等向硬化假设屈服面中心不动、半径扩大,适合单调加载或比例加载路径,如简单拉伸、单向轧制。随动硬化假设屈服面尺寸不变、中心平移,能捕捉包辛格效应,适合反向加载、循环加载,如拉弯矫直、循环弯曲。混合硬化同时允许尺寸扩大和中心平移,用参数 β 分配两者比例,是复杂应变路径下最通用的选择。

硬化曲线的输入形式通常为 σ_y(ε_p) 表格或 Swift、Voce 等解析式。例如 Swift 形式:

σ_y = K · (ε_0 + ε_p)^n

其中 K 为强度系数,n 为硬化指数,ε_0 为预应变。对深冲用低碳钢板,n 值通常在 0.2 左右;对高强钢,n 值偏低但 K 值高。选型时要确认硬化曲线是否覆盖实际成形中的应变范围,外推段不宜过长。

塑性本构三要素与金属成形问题匹配关系 成形问题特征 致密金属 · 单调加载 低应力三轴度 反向/循环加载 包辛格效应显著 高应力三轴度 孔洞/损伤演化 多孔/粉末坯料 静水压力敏感 非比例应变路径 屈服面旋转/畸变 本构三要素选型 ① 屈服准则 von Mises / Tresca(致密金属) GTN(孔洞损伤) Drucker–Prager(压力敏感) ② 流动法则 关联流动 g = f(金属默认) 塑性应变增量 ⊥ 屈服面 塑性变形不可压 ③ 硬化规律 等向硬化(单调加载) 随动硬化(反向加载) 混合硬化(复杂路径) 对结果的影响 成形载荷精度 屈服面形状 回弹预测 残余应力分布 断裂应变预测 损伤演化 致密度/压实 体积变化 循环软化/硬化 棘轮效应

四、匹配方法:从问题反推模型组合

实际选型可以按以下顺序推进。第一步,判断应力三轴度范围:若成形区三轴度长期低于 0.5 且无损伤累积,选 von Mises;若局部三轴度超过 1.0 且需预测断裂,选 GTN 或带损伤的 von Mises。第二步,判断加载路径是否反向:单调加载选等向硬化;存在明显反向加载选随动或混合硬化。第三步,确认流动法则:金属默认关联流动。第四步,标定硬化曲线:用单轴拉伸数据拟合 Swift 或 Voce 参数,并检查拟合区间是否覆盖目标应变。

以一个汽车结构件冲压为例:材料为 DP590 高强钢,板厚 1.5 mm,成形区最大等效塑性应变约 0.25,局部圆角处存在反向弯曲。选型结果为:von Mises 屈服准则 + 关联流动法则 + 混合硬化(β 取 0.3~0.5)。若忽略随动硬化部分,回弹预测的侧壁卷曲量会偏小,因为包辛格效应未被捕捉。若把硬化曲线只标定到 0.15 应变,则圆角处 0.25 应变的流动应力会被外推低估,导致成形载荷偏小。

选型检查口诀:先看三轴度定准则,再看反向定硬化,金属默认关联流,硬化曲线覆盖应变。

五、常见误配与判断依据

误配一:用等向硬化模拟循环弯曲。结果是反向屈服应力被高估,回弹量偏小。判断依据:若加载历史中存在应力反向且反向屈服明显低于正向,必须引入随动硬化。

误配二:用 von Mises 模拟粉末压制。结果是静水压力对屈服的影响被忽略,压坯密度分布预测失真。判断依据:若材料存在初始孔隙且压制中孔隙显著闭合,应选 Drucker–Prager 或 Cam-Clay 类模型。

误配三:硬化曲线外推过长。结果是高应变区流动应力失真。判断依据:外推段应变超过标定最大应变 30% 以上时,应补充更高应变的试验数据(如压缩试验或 bulge 试验)。

误配四:在金属中采用非关联流动。结果是本构切线不对称,收敛变差且无实验依据。判断依据:除非有明确的塑性势实验,否则金属一律关联流动。

选型的最终验证不在模型本身,而在结果:用同一组试验(如载荷–位移曲线、回弹轮廓、断裂位置)对比模拟输出,若偏差超出工程允许范围,回到三要素中逐项排查,而不是盲目加密网格或调整求解器参数。