3. 高速碰撞材料本构选择

在高速碰撞仿真中，材料本构模型决定了应力如何随应变、应变率和温度演化，也决定了单元何时失效、何时被删除。同一套网格、同一套接触参数，仅更换本构模型，峰值冲击力、侵彻深度和失效形貌就可能出现显著差异。因此，本构选择不是“选一个看起来高级的模型”，而是根据材料类型、应变率范围、温度变化幅度和失效模式，选择能覆盖这些物理机制的最小充分模型。

一、四类常用本构的适用边界

Johnson-Cook（JC）是高速碰撞中最常用的金属本构，其流动应力写成三个独立因子的乘积形式：

σ = (A + B·ε_p^n) · (1 + C·ln ε̇*) · (1 − T*^m)

其中：
ε̇* = ε̇_p / ε̇_0        （无量纲塑性应变率，ε̇_0 通常取 1 s⁻¹）
T*  = (T − T_r) / (T_m − T_r)  （同系温度）
A   ：准静态屈服应力
B, n：应变硬化系数与指数
C   ：应变率敏感系数
m   ：温度软化指数

JC 的优势在于参数物理意义清晰、参数标定路径成熟，且能同时表达应变硬化、应变率强化和温度软化。它适合应变率跨越 10²～10⁵ s⁻¹、温升可达数百开尔文的场景，例如弹体侵彻钢板、碎片撞击铝合金蒙皮。

Cowper-Symonds（CS）则把应变率效应写成动态放大因子形式：

σ_d / σ_s = 1 + (ε̇_p / D)^(1/p)

D：参考应变率相关常数
p：应变率敏感指数

CS 通常与塑性随动硬化模型组合使用，适合应变率效应明显但温度软化不占主导的工况，例如船体碰撞、汽车碰撞中的低碳钢薄壁结构。它的参数少、标定快，但对高温软化的表达能力弱于 JC。

塑性随动硬化（*MAT_PLASTIC_KINEMATIC）本身描述的是屈服面随塑性应变平移的硬化行为，可叠加 CS 应变率因子。它适合循环加载或反向加载明显的场景，例如多次撞击、回弹与再压缩。若只关心单次高速侵彻，各向同性硬化往往已足够；若存在明显包辛格效应，则应启用随动硬化并给出硬化参数 β。

应变率效应与温度软化并非独立模型，而是上述本构中的两个因子。高速碰撞中，塑性功大部分转化为热，温升 ΔT 可近似估算为：

ΔT = β_heat · ∫σ dε_p / (ρ · C_p)

β_heat：塑性功转热系数，常用 0.9
ρ     ：密度
C_p   ：比热容

当温升接近材料熔点的一定比例时，温度软化项不可忽略，否则会高估侵彻阻力。

二、失效参数标定方法

JC 失效模型常用等效塑性应变形式：

ε_f = (D1 + D2·exp(D3·σ*)) · (1 + D4·ln ε̇*) · (1 + D5·T*)

σ* = σ_m / σ_eq  （应力三轴度）

标定 D1～D3 需要至少三组不同应力三轴度的试样：光滑圆棒拉伸（σ*≈0.33）、缺口圆棒拉伸（σ*≈0.6～0.8）、纯剪切或扭转（σ*≈0）。通过试验获得各试样的断裂等效塑性应变，再拟合得到 D1、D2、D3。D4 和 D5 分别通过不同应变率下的拉伸试验和不同温度下的拉伸试验确定。若数据有限，可先固定 D4=D5=0，仅标定 D1～D3，并在仿真中通过单元删除阈值进行敏感性验证。

三、参数选择与验证流程

实际建模时，可按以下判断路径选择本构：若材料为金属且应变率超过 10³ s⁻¹、温升显著，优先选 JC；若为低碳钢薄壁结构、应变率在 10²～10³ s⁻¹、温升有限，可选 CS + 塑性随动硬化；若存在明显反向加载，启用随动硬化；若为复合材料或混凝土，则需换用相应本构，不能套用上述金属模型。

参数标定完成后，必须做三类验证：一是用单单元模型复现试验应力-应变曲线，检查各因子是否在目标应变率与温度下正确生效；二是用平板撞击或泰勒杆试验对比剩余长度与变形形貌；三是检查失效单元删除是否造成能量异常，若删除单元过多导致能量不守恒，应调整失效应变或网格尺寸。

图中左侧是工况判断入口，中间列出三类本构及其核心因子，右侧是分步标定与验证流程，底部虚线表示参数迭代修正。实际项目中，往往需要两到三轮“标定—验证—修正”才能得到可用于工程预测的参数集。

总结来说，本构选择的核心是“机制覆盖”：应变硬化、应变率强化、温度软化、随动硬化和失效，哪些机制在本工况中不可忽略，就选能表达这些机制的模型。参数标定则遵循“分步解耦、逐项验证”的原则，先静态后动态、先常温后高温、先硬化后失效，最终用独立试验验证参数集的预测能力。