3. 材料本构与失效模型

冲击分析中，材料模型往往比网格和接触更能决定结果的可信度。准静态下一条弹塑性曲线就够用的金属，在高速撞击下会同时出现屈服提高、硬化减弱、温度上升、损伤累积直至断裂。Abaqus/Explicit 把这些行为拆成三块：塑性本构、损伤起始、损伤演化与单元删除。理解它们各自管什么、参数从哪来、在关键字里怎么落地，是本节的核心任务。

一、金属塑性：从屈服面到硬化规律

金属塑性在 Abaqus 中通过 *PLASTIC 定义，输入的是真实应力—真实塑性应变对。第一行必须给出初始屈服应力与零塑性应变，后续每行给出硬化段。对冲击问题，常用的是等向硬化（Isotropic），即屈服面随塑性应变均匀扩张；若需要考虑包辛格效应（如循环加载或回弹），才切换到随动硬化或组合硬化。

工程应力—应变必须先转换：真实应力 σ_t=σ_eng(1+ε_eng)，真实应变 ε_t=ln(1+ε_eng)，塑性应变 ε_pl=ε_t−σ_t/E。这一步做错，屈服应力和硬化斜率都会系统性偏移。例如某低碳钢 E=210 GPa，工程屈服 350 MPa、工程应变 0.0017，则真实屈服约 350.6 MPa，塑性应变约 0.00003，量级虽小但不可省略。

二、应变率效应与 Johnson-Cook 模型

冲击下应变率可达 10²～10⁴ s⁻¹，屈服应力随率升高而升高。Johnson-Cook 本构把这一现象写成三项乘积：

σ = (A + B·ε_pl^n) · (1 + C·ln(ε̇/ε̇₀)) · (1 − T*^m)
T* = (T − T_room) / (T_melt − T_room)

其中 A 为参考应变率下的屈服应力，B、n 控制硬化，C 控制率敏感度，m 控制热软化，ε̇₀ 通常取 1 s⁻¹。在 Abaqus 中通过 *PLASTIC, HARDENING=JOHNSON COOK 输入，并配合 *RATE DEPENDENT, TYPE=JOHNSON COOK 给出 C 与 ε̇₀，热软化项则需 *INELASTIC HEAT FRACTION 与材料热参数配合。典型示例参数（教学假设值）：A=350 MPa，B=275 MPa，n=0.36，C=0.022，m=1.0，T_melt=1800 K。这些数值仅用于演示输入格式与量级，不代表任何具体牌号的实测结果。

三、损伤起始：何时开始“坏”

损伤起始描述材料从均匀变形转向局部化的临界点。Abaqus 提供多种判据，冲击分析最常用的是 Johnson-Cook 损伤起始（JC Damage）：

ε_f = [D1 + D2·exp(D3·σ*) ] · [1 + D4·ln(ε̇/ε̇₀)] · [1 + D5·T*]
σ* = p / q   （应力三轴度，p 为静水压，q 为 Mises 等效应力）

当累积损伤 ω=Σ(Δε_pl/ε_f) 达到 1 时，损伤起始。D1～D5 需由不同应力三轴度下的断裂试验拟合。关键字为 *DAMAGE INITIATION, CRITERION=JOHNSON COOK。若缺少完整数据，也可用等效塑性应变判据 CRITERION=DUCTILE 作为简化，但应清楚它忽略了应力状态影响。

四、损伤演化与单元删除策略

损伤起始只是“开始”，演化才决定承载能力如何衰减。常用 *DAMAGE EVOLUTION, TYPE=DISPLACEMENT，用等效塑性位移 u_pl 控制软化段，避免网格尺寸敏感。软化曲线可选线性、指数或表格。当损伤变量 D 达到 1，单元刚度归零，若开启 *SECTION CONTROLS, ELEMENT DELETION=YES，该单元被删除，裂纹得以自然扩展。

单元删除不是越早越好。删除过早会丢失质量与能量，删除过晚则单元畸变导致时间步骤降。工程做法是：特征单元尺寸控制在 0.1～0.5 mm，u_pl 取 0.05～0.2 mm 量级，并监控 ALLKE/ALLIE 与伪应变能 ALLAE，确保沙漏与删除不主导响应。

五、参数落地与结果判读

把上述内容串成一次可执行的设置流程：先在材料卡片中定义密度、弹性模量与泊松比；再写 *PLASTIC, HARDENING=JOHNSON COOK 给出 A、B、n、m、T_melt、T_room；接着 *RATE DEPENDENT, TYPE=JOHNSON COOK 给 C 与 ε̇₀；然后 *DAMAGE INITIATION, CRITERION=JOHNSON COOK 给 D1～D5；最后 *DAMAGE EVOLUTION, TYPE=DISPLACEMENT 给 u_pl 并开启单元删除。每一步都对应一个物理机制，缺一项就会在结果中留下可识别的偏差。

结果判读时重点看三处：一是冲击波前沿的应力平台是否与率相关屈服一致；二是损伤云图是否从应力三轴度高的区域（如拉伸侧、剪切带）起始；三是删除区域是否形成连贯裂纹而非零散空洞。若删除过早，检查 u_pl 是否偏小、网格是否过粗；若裂纹不扩展，检查 D1～D5 是否过大或单元删除未开启。能量曲线中 ALLAE 持续上升通常意味着沙漏控制不足，需要调整 *SECTION CONTROLS 的沙漏刚度或细化网格。