4. 金属薄壁材料本构与失效模型
薄壁结构在碰撞中能否被正确“压皱”而不是“撕裂”，取决于材料模型能否同时描述三件事：屈服面如何随塑性变形移动、流动应力如何随应变率升高、以及损伤在何种应力状态下累积到断裂。本节围绕弹塑性随动硬化、Cowper-Symonds 应变率模型、GISSMO 与 Johnson-Cook 失效准则、单元删除和应力三轴度影响展开，目标是让读者能为一根典型薄壁吸能管写出可运行的材料卡片，并理解每个参数在结果中的表现。
一、弹塑性随动硬化：让包辛格效应可被描述
金属薄壁件在压溃中经历反复的弯曲—反向弯曲，若只用各向同性硬化，反向屈服应力会被高估，回弹和残余变形失真。随动硬化把屈服面中心 α 随塑性应变平移，屈服条件写为
f = J₂(σ − α) − σ_y(ε_p) = 0在 LS-DYNA 中，*MAT_PLASTIC_KINEMATIC（MAT_003）用 β 参数在 0～1 之间插值：β=0 为纯随动硬化，β=1 为纯各向同性硬化。对低碳钢薄壁管，常用 β=0.2～0.5，兼顾包辛格效应与硬化记忆。典型 Q235 参数示例（教学假设值）：ρ=7.85e-9 t/mm³，E=210 GPa，ν=0.3，σ_y=235 MPa，E_t=1.2 GPa，β=0.3，失效应变 0.35。E_t 由单轴拉伸曲线塑性段斜率估算，若只有屈服强度，可按 E/100 量级初估再标定。
二、应变率效应：Cowper-Symonds 的用法
碰撞应变率可达 10²～10³ s⁻¹，流动应力相对准静态明显提高。Cowper-Symonds 模型给出
σ_d / σ_s = 1 + (ε̇ / C)^(1/P)其中 C 单位为 s⁻¹，P 为无量纲指数。对低碳钢常用 C=40 s⁻¹、P=5；对高强钢 C 可到 800、P≈7。在 MAT_003 中通过 *DEFINE_CURVE 定义应变率缩放曲线，或直接填写 C、P 由求解器计算。需要注意：应变率取自单元当前等效塑性应变率，网格尺寸过大会平滑掉局部高应变率，导致吸能偏低。建议压溃区单元特征尺寸 2～5 mm，并做一次网格敏感性检查。
三、GISSMO 与 Johnson-Cook 失效准则
Johnson-Cook 失效准则把失效应变写成三轴度、应变率和温度的函数：
ε_f = [D₁ + D₂ exp(D₃ η)] [1 + D₄ ln(ε̇*)] [1 + D₅ T*]对薄壁钢，D₁ 常取 0.1～0.3，D₂ 取 0.5～1.5，D₃ 取负值（如 −1.5），使高三轴度下失效应变迅速下降。GISSMO 则采用损伤变量 D 的增量累积：
ΔD = (n / ε_f(η)) · (ε_p)^(n−1) · Δε_p当 D 达到 1 时触发失效，并可通过 *MAT_ADD_EROSION 或 *MAT_GISSMO 控制单元删除。GISSMO 的优势在于能描述非线性损伤累积和“损伤—软化”耦合，适合薄壁撕裂路径预测。实操中建议先用 JC 准则快速标定，再用 GISSMO 修正三轴度曲线。
四、单元删除与应力三轴度影响
单元删除不是“删掉就完事”，它必须与失效应变、网格尺寸和沙漏控制配合。删除过早会让压溃力偏低、吸能不足；删除过晚则单元畸变、时间步骤降甚至负体积终止。应力三轴度 η=σ_m/σ_eq 是关键变量：薄壁弯曲外侧受拉 η≈0.33～0.58，失效应变低；内侧受压 η 为负，几乎不失效。因此同一根管在弯曲处先裂、在平直段后裂，正是 η 分布不均的结果。
五、可执行的材料卡片示例
*MAT_PLASTIC_KINEMATIC
$# mid ro e pr sigy etan beta
1 7.85e-09 210000 0.3 235 1200 0.3
$# src srp fs vp
40.0 5.0 0.35 0.0
*MAT_ADD_EROSION
$# mid exel mxpres
1 0 0
$# mneps effeps voleps numfip
0.0 0.35 0.0 1.0
该卡片中，src=40、srp=5 即 Cowper-Symonds 参数；fs=0.35 为失效应变；*MAT_ADD_EROSION 用 effeps=0.35 触发删除。若改用 GISSMO，则替换为 *MAT_GISSMO 并给出 η-ε_f 曲线和损伤指数 n。运行后重点看三个量：压溃力峰值、平台力均值、总吸能。若平台力波动超过 ±15%，优先检查网格尺寸和删除阈值，而不是先调材料强度。
把上述参数写成卡片只是起点。真正决定仿真可信度的是：随动硬化是否匹配反向加载、应变率曲线是否覆盖实际碰撞速率、失效准则是否在三轴度范围内标定、删除阈值是否与网格尺寸匹配。建议每改一个参数，只观察一个输出量，逐步逼近试验曲线，而不是一次性调完所有参数。