3. 冲压材料模型与参数获取

冲压仿真的可信度，首先取决于材料模型是否真实。几何、摩擦、压边力都可以调，但材料本构一旦失真，回弹预测、开裂判定和减薄分布都会系统性偏离。本节围绕“从实验数据到可求解材料卡片”这条主线，讲清弹塑性本构、各向异性屈服、硬化曲线拟合以及材料卡片建立的关键处理方式。

一、弹塑性本构：仿真里材料如何“先弹后塑”

金属薄板在冲压中经历弹性变形、屈服、塑性流动和卸载回弹。仿真中通常采用小弹性、大塑性的弹塑性本构：弹性阶段由弹性模量 E 和泊松比 ν 描述，塑性阶段由屈服函数与硬化规律共同决定。对多数汽车覆盖件用钢板，E 取 200～210 GPa、ν 取 0.3 是常见示例值；铝合金 E 约 70 GPa、ν 约 0.33。这些数值必须与实验或材料手册一致，不能随意套用。

塑性部分遵循三个基本假设：屈服面、流动法则和硬化律。屈服面定义“多大应力开始塑性变形”，流动法则决定塑性应变增量方向，硬化律描述屈服面随塑性应变如何扩大或移动。冲压仿真中最常用的是关联流动法则配合等向硬化，即塑性应变增量垂直于屈服面，屈服面随等效塑性应变单调扩大。对回弹敏感的高强钢，有时需要引入随动硬化或混合硬化来捕捉包辛格效应，否则反向加载时的屈服偏低会被忽略。

二、各向异性屈服：Hill 与 Barlat 的适用边界

轧制薄板存在明显的各向异性：0°、45°、90°方向的屈服强度和 r 值不同。若直接使用各向同性 Mises 屈服，仿真出的制耳、拉深高度和回弹方向都会失真。工程上常用两类各向异性屈服准则。

Hill 1948 是经典二次屈服准则，形式简单、参数少，适合各向异性不太强、以拉深为主的低碳钢。它通过 L、M、N 等各向异性系数把不同方向的屈服应力联系起来，参数可由单向拉伸的屈服应力和 r 值反算。缺点是二次形式对某些铝合金、高强钢的屈服面描述偏保守，尤其在平面应变和双拉状态下误差较大。

Barlat（如 Yld2000-2d） 属于非二次屈服准则，能更灵活地描述面内各向异性，适合铝合金、高强钢和复杂拉深件。它通过两个线性变换张量构造等效应力，需要 0°、45°、90°方向的屈服应力和 r 值，必要时补充双拉试验数据。参数越多，拟合越灵活，但对实验数据质量要求也越高。若只有单向拉伸数据，强行使用 Barlat 反而可能过拟合。

三、硬化曲线拟合：把拉伸曲线变成仿真可用的数学表达

实验得到的是真实应力-真实应变曲线，但仿真求解器需要的是等效塑性应变与流动应力的关系。第一步是数据转换：把工程应力-应变转换为真实应力-应变。真实应力 σ_t = σ_e(1+ε_e)，真实应变 ε_t = ln(1+ε_e)。第二步是扣除弹性段，得到塑性应变 ε_p = ε_t − σ_t/E。第三步才是拟合硬化模型。

常用模型包括 Swift、Ludwik、Hockett-Sherby 和分段线性。Swift 形式 σ = K(ε₀+ε_p)ⁿ 适合大多数低碳钢，参数 K、n、ε₀ 可通过最小二乘拟合得到。例如某 DC04 示例数据拟合后 K≈520 MPa、n≈0.22、ε₀≈0.005，R² 可达 0.998。但要注意：拉伸试验的均匀延伸率通常只有 20%～25%，而冲压局部应变可能超过 40%，外推段必须谨慎。若外推过高，仿真会低估开裂风险；外推过低，则过早判定开裂。工程上可结合液压胀形或剪切试验补充大应变段数据。

四、材料卡片建立与实验数据转换

材料卡片是求解器能直接读取的材料定义文件。以常用显式求解器为例，卡片通常包含：弹性参数（E、ν、密度）、塑性硬化曲线（分段点表或解析参数）、各向异性参数（r 值或屈服应力比）、单位制声明。建立时最容易出错的是单位制不统一：冲压仿真常用 mm-ton-s-MPa，密度应写 7.85e-9 t/mm³，而不是 7850 kg/m³。另一个常见问题是 r 值方向定义与轧制方向不一致，导致各向异性“转错方向”。

实验数据转换建议按固定流程执行：先确认试验方向与轧制方向对应关系，再统一转换为真实应力-真实应变，然后扣除弹性段得到塑性段，最后按所选屈服准则反算各向异性系数。若使用 Hill 1948，可由 r₀、r₄₅、r₉₀ 计算 L、M、N；若使用 Barlat，则需同时输入屈服应力和 r 值，并检查拟合后的屈服面是否与实验点吻合。

五、参数获取的实践判断

材料模型没有“越复杂越好”。若零件以拉深为主、各向异性中等，Hill 1948 配合 Swift 硬化往往足够；若零件涉及铝合金覆盖件、复杂翻边或回弹敏感区，Barlat 配合 Hockett-Sherby 更合适。关键是让模型复杂度与实验数据质量匹配：只有单向拉伸数据时，强行标定 Barlat 的多个参数容易过拟合，反而降低预测稳定性。

最终交付的材料卡片应附带一份参数来源说明：每个参数的实验来源、拟合方法、拟合优度、外推范围以及适用应变区间。这样在后续仿真出现异常时，才能快速判断是材料模型问题还是工艺参数问题，而不是盲目试错。