3. 材料本构与失效模型入门
在显式动力学冲击仿真中，材料模型决定了计算结果是否可信。冲击问题的特点是高应变率、大变形、局部温升和快速失效，因此不能简单套用静力分析中的线弹性模型。本节从最基础的弹性开始，逐层引入塑性、应变率效应、Johnson-Cook本构以及损伤起始与演化准则，帮助建立“从可逆变形到断裂”的完整认识。
一、弹性：一切本构的起点
弹性阶段描述材料在卸载后能完全恢复的变形行为，其核心参数是弹性模量 E 和泊松比 ν。在 Abaqus/Explicit 中，弹性通过 *ELASTIC 关键字定义，各向同性材料只需给出 E 和 ν。例如结构钢常取 E = 210 GPa、ν = 0.3，密度约 7850 kg/m³。弹性模型本身不产生永久变形，也不消耗能量，它只负责把应力与弹性应变线性关联起来。冲击仿真中弹性参数主要影响波速和初始刚度，波速 c = √(E/ρ)，这直接决定了稳定时间增量的大小。
二、塑性：屈服与硬化
当应力超过屈服极限，材料进入塑性阶段，产生不可恢复的变形。Abaqus 中塑性通过 *PLASTIC 定义，输入的是真实屈服应力与真实塑性应变的关系曲线。以常用的金属为例，可以给出如下硬化数据：
*PLASTIC
355.0, 0.000
420.0, 0.020
480.0, 0.050
520.0, 0.100
550.0, 0.200这里第一列是真实屈服应力（MPa），第二列是真实塑性应变。注意 Abaqus 要求塑性应变从 0 开始，且必须单调递增。塑性模型决定了材料在冲击下能承受多大的应力、发生多大的塑性流动。若只给弹性不给塑性，材料在屈服后会表现为理想弹性回弹，显然不符合金属冲击的真实情况。
三、应变率相关模型：冲击仿真的关键
冲击载荷下应变率可达 10²～10⁴ s⁻¹，远高于准静态实验的 10⁻³ s⁻¹。许多金属的屈服应力随应变率升高而显著提高，这就是应变率强化效应。Abaqus 中常用的应变率相关模型包括：
- Cowper-Symonds 模型：通过
*PLASTIC中的RATE参数定义，形式为σ_d/σ_s = 1 + (ε̇/D)^(1/p)，其中D、p为材料常数。 - Johnson-Cook 模型：将应变硬化、应变率强化和温度软化统一在一个表达式中，是冲击仿真中最常用的本构之一。
- 表格型应变率模型：直接给出不同应变率下的屈服曲线，适合实验数据充分的情况。
选择哪种模型取决于材料、应变率范围和可用数据。对于大多数金属冲击问题，Johnson-Cook 是首选。
四、Johnson-Cook 本构模型
Johnson-Cook 模型的流动应力表达式为：
σ = (A + Bεⁿ)(1 + C ln ε̇*)(1 - T*ᵐ)
其中 A 是准静态屈服应力，B 是硬化模量，n 是硬化指数，C 是应变率敏感系数，m 是温度软化指数。ε̇* = ε̇/ε̇₀ 为无量纲应变率，T* = (T - T_r)/(T_m - T_r) 为同系温度。在 Abaqus 中通过 *PLASTIC, HARDENING=JOHNSON COOK 定义，并配合 *RATE DEPENDENT 和 *INELASTIC HEAT FRACTION 使用。以某装甲钢为例，一组典型参数为：A = 792 MPa、B = 510 MPa、n = 0.26、C = 0.014、m = 1.03、T_m = 1793 K、T_r = 293 K。这些参数通常由霍普金森压杆实验拟合得到，属于示例性数据，实际使用应以材料实验为准。
五、损伤起始与演化准则
损伤起始描述材料从“完整”到“开始出现微裂纹”的临界条件。Johnson-Cook 损伤起始准则的等效塑性应变表达式为：
ε_pl^f = [d₁ + d₂ exp(d₃ η)] [1 + d₄ ln ε̇*] [1 + d₅ T*]
其中 η = -p/q 为应力三轴度，p 为静水压力，q 为 Mises 等效应力。当累积损伤参数 ω = Σ Δε_pl / ε_pl^f 达到 1 时，损伤起始。在 Abaqus 中通过 *DAMAGE INITIATION, CRITERION=JOHNSON COOK 定义，需要给出 d₁～d₅ 五个参数。例如某铝合金的示例取值为 d₁ = 0.071、d₂ = 1.248、d₃ = -1.142、d₄ = 0.147、d₅ = 0.0。
损伤起始之后，材料并未立即失去承载能力，而是进入损伤演化阶段。Abaqus 中通过 *DAMAGE EVOLUTION 定义，常用方式有两种：基于等效塑性位移的线性软化，或基于断裂能 G_f 的能量型演化。能量型演化能有效缓解网格依赖性，推荐优先使用。当损伤变量 D 达到 1 时，单元被删除，裂纹形成。
*SECTION CONTROLS 中的 ELEMENT DELETION=YES。若仿真中出现单元过度畸变导致计算中止，可适当调整损伤演化中的断裂能或引入黏性正则化。
六、参数获取与验证思路
Johnson-Cook 参数通常通过以下实验组合获取：准静态拉伸得到 A、B、n；霍普金森压杆（SHPB）在不同应变率下得到 C；高温拉伸或压缩得到 m；缺口拉伸与剪切实验结合得到损伤参数 d₁～d₅。若缺乏实验条件，可参考同类材料的文献值，但必须通过标准试样的仿真-实验对比进行校核。验证时重点关注：屈服应力-应变曲线是否吻合、颈缩位置是否一致、断裂应变是否在合理范围。只有本构与损伤参数同时可信，冲击仿真结果才具备工程参考价值。