4. 材料模型基础与金属塑性

在整车碰撞仿真中，金属薄壁件的吸能行为几乎决定了乘员舱的完整性与加速度波形。LS-DYNA 提供了大量金属材料模型，其中 *MAT_PIECEWISE_LINEAR_PLASTICITY（常写作 MAT_024）是车身钣金件最常用的弹塑性本构之一。它用一条分段线性硬化曲线描述屈服后的应力—塑性应变关系，并可选地叠加应变率缩放与失效应变判据，从而在保证计算效率的同时覆盖从准静态到高速冲击的宽应变率范围。

一、MAT_024 的本构结构

MAT_024 属于各向同性弹塑性模型，其核心由三部分组成：弹性段、塑性硬化段和失效段。弹性段由弹性模量 E 与泊松比 PR 定义；塑性段由屈服应力 SIGY 与一条“塑性应变—真实应力”曲线（通过 LCSS 指定曲线 ID）共同描述；失效段则由失效应变 FAIL 或等效塑性应变判据控制。当单元等效塑性应变达到 FAIL 时，该单元被删除，不再参与后续承载。

需要特别强调的是：MAT_024 的硬化曲线输入的是真实应力—真实塑性应变，而不是工程应力—工程应变。工程数据需先转换：真实应力 σ_t=σ_e(1+ε_e)，真实应变 ε_t=ln(1+ε_e)，塑性应变 ε_p=ε_t−σ_t/E。若直接输入工程曲线，屈服后应力会被系统性低估，导致吸能量偏小、侵入量偏大。

二、应变率效应与 Cowper-Symonds 模型

车身用钢在碰撞中应变率可达 10²～10³ s⁻¹，屈服强度相比准静态可提升 30%～80%。MAT_024 通过两个参数 C 和 P 启用 Cowper-Symonds 缩放因子：

σ_y(ε̇) = σ_y0 · [1 + (ε̇ / C)^(1/P)]

其中 σ_y0 为当前塑性应变对应的准静态屈服应力，ε̇ 为等效塑性应变率。典型低碳钢 C≈40 s⁻¹、P≈5；高强钢 C≈3200 s⁻¹、P≈5。若 C 与 P 均置零，则关闭应变率效应，模型退化为率无关塑性。工程上常通过分离式霍普金森压杆（SHPB）试验获得不同应变率下的应力—应变曲线，再拟合 C、P。

三、失效应变与单元删除

MAT_024 的失效判据基于等效塑性应变：当 ε_p ≥ FAIL 时单元失效并被删除。删除机制通过 *CONTROL_TERMINATION 或关键字中的 ERODE 选项控制。删除后，原本由该单元承担的载荷会重新分配到相邻单元，因此网格尺寸对结果影响显著——粗网格下单个单元删除会造成更大的“缺口”，使结构过早失去承载能力。

失效应变的取值通常来自单轴拉伸试验的断裂塑性应变，但需注意：试验测得的是工程断裂应变，需转换为真实塑性应变；同时薄壁件在碰撞中处于多轴应力状态，单轴失效应变往往偏保守。工程实践中常结合 *MAT_ADD_EROSION 引入应力三轴度相关的失效准则，或直接采用 GISSMO 模型进行更精细的损伤演化描述。

四、参数设置示例

以某 590DP 高强钢为例，典型 MAT_024 卡片参数如下（示例值，需按实际材料试验替换）：

硬化曲线 LCSS=100 以“塑性应变—真实应力”对的形式给出，例如：0→420、0.02→510、0.05→580、0.10→650、0.20→720 MPa。曲线应覆盖至失效应变附近，避免外插。

五、网格与删除的耦合影响

单元删除并非“越早越好”。若 FAIL 偏小，结构在吸能不足时即被删除，碰撞力波形会出现异常跌落；若 FAIL 偏大，单元过度拉伸导致负体积或时间步骤降。建议在标定阶段用单个帽形梁或三点弯曲模型做网格敏感性分析：分别取 2 mm、4 mm、6 mm 特征尺寸，观察力—位移曲线与失效位置是否收敛。通常 3～5 mm 的壳单元在精度与效率间较平衡。

六、常见问题与排查

• 力—位移曲线出现锯齿或突降：多为单元删除后载荷重分配导致，可适当加密网格或改用 *MAT_ADD_EROSION 平滑失效。
• 时间步骤降：删除后剩余单元畸变严重，检查 FAIL 是否过大，或启用 *CONTROL_TIMESTEP 的 DT2MS 质量缩放。
• 吸能量偏低：确认硬化曲线是否误用工程应力，或 C/P 是否遗漏导致率效应未激活。
• 单元负体积：通常伴随 FAIL 过大或接触穿透，可先降低 FAIL 并检查接触刚度。

掌握 MAT_024 的关键在于理解“弹性—硬化—失效”三段式的参数来源与相互制约：硬化曲线决定吸能水平，C/P 决定动态强化幅度，FAIL 决定失效时机。三者需与网格尺寸、接触设置协同标定，才能得到既稳定又贴近试验的碰撞响应。