4. 金属材料本构模型

在LS-DYNA显式动力学仿真中，材料本构模型决定了应力如何随应变、应变率和加载历史演化，是决定回弹预测精度、成形载荷量级乃至计算稳定性的核心输入。金属成形问题中，材料通常经历弹性加载、屈服、塑性流动、卸载与再加载的完整路径，因此不能只给一个弹性模量和密度了事。本节围绕弹塑性框架、硬化规律、应变率效应以及常用材料卡的选用展开，重点说明“什么工况该选哪张卡、卡里每个参数在物理上代表什么”。

一、弹塑性框架：从屈服面到流动法则

金属在低应力阶段表现为线弹性，应力应变关系由弹性模量 E 和泊松比 ν 描述。当等效应力达到屈服应力后进入塑性阶段，此时需要三个要素共同定义：屈服准则（判断何时进入塑性）、流动法则（塑性应变增量的方向）、硬化规律（屈服面如何随塑性变形演化）。LS-DYNA中金属成形最常用的是von Mises屈服准则，即用等效应力 σ_eq 与当前屈服面半径比较，超过则产生塑性流动。塑性变形假定不可压，因此体积模量控制体积响应，剪切模量控制偏斜响应，这也是为什么材料卡中通常输入 E、ν 而非直接输入 K、G。

需要特别强调的是，显式动力学中塑性本构采用“径向返回”算法：先按弹性预测应力，若超出屈服面，则沿屈服面法向把应力拉回，并更新等效塑性应变。这一过程在每个时间步内完成，因此材料卡的硬化曲线必须覆盖仿真中可能出现的最大塑性应变范围，否则LS-DYNA会外插甚至报错。

二、硬化规律：各向同性硬化与随动硬化

各向同性硬化（Isotropic Hardening）假定屈服面中心不动、半径随塑性应变增大而扩大。其数学表达为 σ_y = σ_y0 + K·(ε_p)^n 或分段线性形式。它适合单调加载、大塑性变形主导的工况，例如深冲、胀形、挤压，因为这类问题中包辛格效应不显著。缺点是反向加载时无法体现“反向屈服提前”的现象，回弹预测会偏硬。

随动硬化（Kinematic Hardening）假定屈服面半径不变、中心随塑性应变移动，从而自然描述包辛格效应。它适合存在明显反向加载的工况，例如拉弯矫直、反复弯曲、回弹后反向校形。LS-DYNA中 MAT_PLASTIC_KINEMATIC 通过一个硬化参数 β 在两者之间插值：β=0 为纯随动硬化，β=1 为纯各向同性硬化，中间值表示混合硬化。实际板料成形中常取 β=0.1~0.3，以兼顾单调段强度与反向段包辛格效应。

三、应变率相关模型

金属在高速成形（冲压速度高、爆炸成形、高速撞击）下屈服应力会随应变率提高而上升，这一现象称为应变率强化。LS-DYNA中常用Cowper-Symonds模型描述：σ_y = σ_y0·[1 + (ε̇/C)^(1/P)]，其中 C 和 P 为材料常数。例如低碳钢典型值 C≈40 s⁻¹、P≈5；铝合金 C≈6500 s⁻¹、P≈4。当仿真中最大应变率低于约 10² s⁻¹ 时，率效应通常可忽略；当冲压速度达到 5~10 m/s 以上，或涉及冲击、爆炸时，必须启用率相关项，否则会低估成形载荷和回弹。

更精细的率相关模型包括 MAT_JOHNSON_COOK（同时考虑应变硬化、应变率强化和温度软化）和 MAT_PIECEWISE_LINEAR_PLASTICITY（可输入多条不同应变率下的应力应变曲线）。Johnson-Cook模型形式为 σ = (A + B·ε_p^n)·(1 + C·ln ε̇*)·(1 - T*^m)，适合高温、高应变率耦合的成形与冲击问题，但参数标定成本高，需要SHPB实验和高温拉伸数据。

四、常用材料卡选用对照

五、参数获取与选用实操

以某低碳钢板料冲压为例，假设通过单轴拉伸获得：E=210 GPa、ν=0.3、ρ=7850 kg/m³、σ_y0=180 MPa、E_t=1.2 GPa（塑性段斜率）。若采用 MAT_PLASTIC_KINEMATIC，则输入 σ_y=180、E_t=1200（单位MPa），β=0.2，C=40、P=5。这里 E_t 是双线性硬化斜率，与真实曲线在最大塑性应变处的割线斜率对应；若仿真中最大等效塑性应变约0.25，则 E_t 应取该点割线值，而不是初始塑性段切线值。

选用流程可归纳为：先判断是否存在反向加载——有则选随动或混合硬化，无则各向同性即可；再判断应变率是否超过 10² s⁻¹——超过则启用率相关项；最后判断是否需要温度耦合——需要则升级到Johnson-Cook。对于大多数常温板料冲压，MAT_PLASTIC_KINEMATIC 配合 β=0.1~0.3 和Cowper-Symonds参数已足够；对于回弹精度要求极高的高强钢，建议改用 MAT_PIECEWISE_LINEAR_PLASTICITY 并输入实测多应变率曲线。