4.2 高应变率本构：Johnson-Cook 模型实战

4.2 高应变率本构：Johnson-Cook 模型实战

Johnson-Cook 模型用一个五参数乘积式描述金属的动态流动应力：应变硬化、应变率强化、温度软化三项相乘。它是爆炸与侵彻数值模拟中使用率最高的本构模型（第 5 章的直接输入），也是勘测者从断口反推载荷时评估"材料当时有多硬"的随身工具。

上一节说明了塑性波如何"施工"，本节回答施工阻力怎么算：同一种钢，静态拉伸的强度约几百兆帕，爆炸载荷下瞬时流动应力可能翻倍——差值来自应变率与温度的竞争。

三因子相乘的物理含义

流动应力写成三项乘积：第一项是应变硬化（塑性越大力发展越强，与静态一致）；第二项是应变率强化（变形速度每提高一个数量级，应力上升一截，对数关系）；第三项是温度软化（材料温度接近熔点时软化到趋于零）。乘积形式的好处是各项可独立标定：准静态试验定应变硬化项，霍普金森杆系列试验定率强化项，高温快速拉伸定软化项。

对勘测者的直觉意义：应变率让材料更硬，绝热升温让材料更软，两者在同一毫秒内赛跑。爆炸载荷初期率强化占上风（材料显得格外硬），随后塑性功转化为热量来不及散出（绝热假设），温度追上来把材料软化——这个此消彼长正是 4.3 节绝热剪切带的成因伏笔。

演算：一块装甲钢的动态强度

以典型装甲钢参数为例，用 Python 算三个工况对比：

# Johnson-Cook 模型：装甲钢动态流动应力估算
# sigma = (A + B*eps^n) * (1 + C*ln(eps_dot/eps0)) * (1 - Tm**m)
import math
A, B, n = 792e6, 510e6, 0.26      # 应变硬化参数 Pa
C, eps0 = 0.014, 1.0              # 率强化参数，参考应变率 1/s
Tmelt, Tr, m = 1793.0, 300.0, 1.03  # 熔点/参考温度/软化指数
eps = 0.2                          # 等效塑性应变
def jc(eps_dot, T):
    hard = A + B*eps**n
    rate = 1 + C*math.log(eps_dot/eps0)
    Tm = (T - Tr)/(Tmelt - Tr)
    soft = max(1 - Tm**m, 0.0)
    return hard*rate*soft
for label, ed, T in [("准静态 1/s 300K", 1, 300),
                     ("SHPB 2000/s 300K", 2000, 300),
                     ("爆炸级 1e5/s 600K", 1e5, 600)]:
    s = jc(ed, T)
    print(f"{label:>20}: 流动应力 {s/1e6:.0f} MPa")
print(f"准静态屈服参照 A = {A/1e6:.0f} MPa")
​

输出：准静态约 984 兆帕；两千每秒应变率下约 1077 兆帕；十万每秒应变率、六百开温度下约 1138 兆帕（率强化再加一成，但温度已开始抵消）。三个数字放在一起看出两个规律：率强化的增幅是对数型的，十个数量级的应变率只让应力翻倍左右；温度软化的下拉随塑性功累积越来越重，高应变率下温差竞争激烈。

参数敏感性与勘测应用

J-C 参数从公开手册查到后，勘测者最该做的是敏感性检查——哪一项参数错 20%，结论差多少：

⚠️ 常见坑：手册参数都标定于特定热处理状态。同一牌号钢，淬火与退火状态的 A 值可以差三倍——勘测复算前必须核对材料状态，不能只看牌号。

💡 关键直觉：J-C 是"工程够用"模型，不是"物理精确"模型。它的三项乘积假设互相独立，真实材料率强化与温度相关。但对勘测目的——估算载荷量级、核对断口反推——它给出的精度绰绰有余；需要更高精度时换用 Zerilli-Armstrong 等基于位错机制的模型，代价是参数更难获得。

岩土与混凝土：另一套性格

金属之外的介质（岩石、混凝土、土壤）的高压响应还多两个特征：孔隙压实（低压段应力-应变曲线呈凹形，大量能量耗在压扁孔隙上，冲击波在没有孔隙的密实段才真正形成）与压力相关性（波速与刚度随围压显著上升，地层深处的岩石"更硬更快"）。这就是为什么第 3 章的相似律在岩土近场要谨慎使用——耗散机制多了一个"压实吸能"通道。第 8 章爆破块度评估将直接引用本节的压实概念。

失效判据与本构的配对

本构回答"材料多硬"，失效判据回答"何时算坏"，两者必须配对使用。常用失效判据把累积塑性应变与应力三轴度（静水应力与等效应力之比）挂钩：三轴度高（缺口、裂纹尖端）材料延性骤降，失效应变小；纯剪切下延性最好。爆炸载荷的特殊性在于三轴度随位置剧烈变化——弹坑壁附近高压约束延性高，自由表面附近三轴度低易裂。勘测复算模拟结果时，除了看应力云图，务必看三轴度云图：失效集中区几乎总与高三轴度区重合。

模型族谱的选型经验：Johnson-Cook 适合中低应变率到中等温度区间的大量工程合金；Zerilli-Armstrong 基于位错热激活机制，外推性更好但参数稀缺；高应变率与高压并存的极端区间（超速撞击）直接进入流体动力学区，强度项退居次要，状态方程主导。选错区间比选错模型族更致命。

再补一个参数获取的现实路径：公开手册只有常见合金钢与铝合金的 J-C 参数，特种材料与具体热处理批次要靠专项标定——至少一组准静态拉伸、一档霍普金森杆、一次升温试验，才能把五参数锁住。预算紧张时的省法：A 与 B 由静态试验定，C 借用同类材料，m 用熔点近似估算，误差预算相应放宽并在报告中声明。诚实标注的近似，好过假装精确的借用。

本节要点回顾

• 三因子乘积：应变硬化、率强化、温度软化，各因子独立标定、乘积合成；
• 率强化对数弱、温度软化强：高应变率下两者赛跑，决定材料瞬时"性格"；
• 敏感性排序：A 线性最敏感，C 对数弱敏感，m 在高温场景强敏感；
• 岩土多两个通道：孔隙压实吸能与围压相关性，近场相似律慎用。

下一节到断口现场去：层裂与绝热剪切——材料把自己的受载历史写在了断面上，勘测者的工作是把它读出来。