3.1 Rankine-Hugoniot 关系：波阵面上的守恒账本

3.1 Rankine-Hugoniot 关系：波阵面上的守恒账本

兰金-雨贡纽关系是把质量、动量、能量守恒直接跨过冲击波间断面写出的代数关系组：只需知道波速和波前状态，就能解出波后压力、密度、粒子速度与温度。它是全部冲击波理论的骨架，也是勘测者从"测到波速"推"波后载荷"的唯一桥梁。

本章开篇即本节：波阵面是冲击波身上唯一"可计算"的地方，读懂这张账本，后面两节的衰减与反射都是它的推论。

间断面上的三笔账

冲击波阵面是一个近乎理想的间断：厚度只有几个分子自由程，两侧状态发生跳变。正因为薄，守恒定律在阵面两侧取值即可列账，不涉及阵面内部细节——这与第 1 章"守恒不依赖过程"的精神一脉相承。

设波前状态为下标 0、波后为 1，波速为 D，跨阵面的三笔账写成人话：

• 质量账：穿过阵面的质量流连续，压缩前后"单位面积每秒通过的分子数"必须相等；
• 动量账：压力跳变等于质量流乘以速度跳变，压差就是把介质从静止推到粒子速度 u 的原因；
• 能量账：比内能的增加等于压力做的压缩功的平均值。

三笔账联立消元，得到最常用的两条：动量关系 p1 减 p0 等于波前密度乘波速乘粒子速度；以及雨贡纽曲线——波后可能状态在压力-密度平面上的轨迹。强波极限下还有一个极简结论：波后密度最多压缩到波前的六倍（空气，比热比 1.4 时），无论波多强都压不过这道"六倍墙"。

💡 关键直觉：雨贡纽曲线和等熵线是两条不同的线。声波沿等熵线走（弱到没有耗散），冲击波沿雨贡纽线走（内部有不可逆耗散）——两者在起点相切，波越强分岔越大。理解了这个分岔，就理解了冲击波为什么会衰减：耗散在持续付"熵税"。

勘测演练：从波速反推波后状态

现场布设了两组探针，测得某冲击波到达时差换算的波速为 1800 米每秒，环境为常压空气。用 Python 把波后状态完整算出来：

# 由波速 D 反推空气冲击波后状态（理想气体，gamma=1.4）
p0, rho0, gamma = 1.013e5, 1.225, 1.4
D = 1800.0                       # 波速 m/s
c0 = (gamma*p0/rho0)**0.5        # 波前声速
print(f"波前声速 c0 = {c0:.0f} m/s，马赫数 {D/c0:.2f}")
# 粒子速度（动量关系）：u = (p1-p0)/(rho0*D)，先由雨贡纽解 p1
# 理想气体雨贡纽解出密度比的表达式：
mu = (gamma+1)/(gamma-1)
r = (D/c0)**2                    # 马赫数平方
rho1_ratio = ((mu+1)*r - 2) / ((mu-1)*r + 2)   # 波后/波前密度比
u = 2*D/(mu+1)*(1 - 1/r)          # 粒子速度
p1 = p0 + rho0*D*u
T1_ratio = p1/p0/rho1_ratio       # 温度比（理想气体）
print(f"波后/波前密度比 = {rho1_ratio:.2f}（理论上限 {mu:.2f}）")
print(f"波后粒子速度 u = {u:.0f} m/s")
print(f"波后压力 p1 = {p1/1e5:.1f} atm")
print(f"波后温度约 {288*T1_ratio:.0f} K")
​

输出：波前声速 340 米每秒，马赫数 5.3；密度比 4.6（逼近 6 倍墙）；粒子速度约 1365 米每秒；波后压力约 15.9 个大气压；波后温度约 1650 开。一组探针数据就这样变成了完整的载荷谱——这就是阵面账本的威力。

波形的演化：从陡峭到平缓

阵面账本算的是"某一瞬间"，勘测还需要知道波形随距离怎么变。空气中爆炸波的典型压力历程是"陡升缓降"：阵面瞬间跳到峰值超压，随后近似指数衰减，尾巴上还有一段负压相（压力低于大气压）和第二次弱压缩。描述它用三个特征量：峰值超压、正压持续时间、正压冲量（压力对时间的积分）。

随距离增大，三条规律同时发生：峰值超压按距离的立方根比例律衰减（精确说法见 3.3 节）；正压持续时间拉长；波形整体钝化，最终退化为声波。用一张分层透视图总结全流程：

图：三条曲线对应三个距离档。注意负压相只在近中场显著，远场波形几乎是对称的衰减振荡。

⚠️ 常见坑：把波速当成常数去算远场参数。波速随传播衰减不断下降（趋于声速），用近场波速套远场距离会高估载荷。正确做法是对每一档距离分别用当地马赫数进阵面账本。

强波与弱波的两个极限

阵面账本在两个极限下退化成极简形式，勘测快速估算时极其好用。弱波极限（马赫数接近 1）：波后扰动退化为声波，超压与粒子速度之比趋于波阻抗（密度乘声速），一切线性声学关系成立——远场爆炸波的破坏评估就常按线性声学处理。强波极限（马赫数远大于 1）：波后密度锁死在六倍墙，粒子速度逼近波速的六分之五（空气），温度与速度平方成正比地飙升——近场火球的估算用这组关系。

介电离是强波极限的延伸话题：当波后温度高到几千开以上，空气分子开始电离，导电性出现，理想气体状态方程失效，需改用真实气体物态方程。核爆早期火球的雷达信号异常与电离层扰动，就是这一效应的宏观显现。勘测者遇到"爆炸后短波通信中断"的证词时，应意识到这是强波存在的旁证而非幻觉。

补充一条现场速算：马赫数小于一点二时，波后超压可用线性声学近似（超压约等于波阻抗乘粒子速度），误差小于百分之三；马赫数二时这个近似低估四成，必须回完整账本。判断"能不能偷懒"的分界线记在一处：波速接近声速就线性，明显超音速就完整算。

负压相虽强度小，却对轻质物体有独特破坏——屋面被掀、窗玻璃向外飞，都是负压相吸力或随后的稀疏回摆所致。勘测玻璃散落方向可以区分正压段破坏（玻璃向内）与负压段破坏（向外），两者对照能标定波形相位，为反演多一枚时间维度的指纹。

本节要点回顾

• 阵面账本：质量、动量、能量三笔账跨间断列写，波速加波前状态即可解出全部波后参数；
• 六倍密度墙：空气冲击波后密度最多压缩六倍，强波极限的硬边界；
• 雨贡纽线与等熵线分岔：耗散付"熵税"是波衰减的热力学根源；
• 波形三特征量：峰值超压、正压冲量、正压持续时间，勘测反演的三大输入。

下一节看波撞上结构：反射如何把载荷放大，以及勘测者最该警惕的马赫杆。