1.2 现场勘测第一课：守恒定律与控制方程

1.2 现场勘测第一课：守恒定律与控制方程

质量、动量、能量三大守恒定律是爆炸力学全部计算的出发点：对任取的控制体，流入减流出加累积等于源项。勘测中的每一条计算——药量反演、破片初速、容器破裂能量——本质上都是在解这组方程的某种简化形式。本节把守恒定律翻译成勘测语言。

上一节我们判定"能量从哪来"，这一节回答"能量去了哪"。账本的两面合起来，才是一次完整勘测。

守恒定律为什么是勘测的地基

勘测面对的永远是不完整的信息：传感器只记下了半截波形，弹坑被雨水泡塌了一角。守恒定律的价值在于它不依赖过程细节——无论爆炸过程多么混乱，爆炸前后的总质量、总动量、总能量必须平衡。这意味着勘测人员可以用"边界上的测量"去约束"内部发生的过程"：测出破片打墙的冲量，就限定了爆炸气体带走动量的份额；量出弹坑抛土的分布，就限定了用于抛掷的能量比例。

具体写出来，对固定控制体 V（边界为 A），三大守恒的积分形式是：

• 质量守恒：控制体内质量的累积率，等于通过边界流入的质量流。紧凑装药爆炸后产物质量不变，这条看似平凡，却在核对"实际装药量是否与台账一致"时直接可用；
• 动量守恒：作用在控制体上的合外力（压力加应力），等于动量流出率加累积率。破片初速反推全靠它；
• 能量守恒：进入控制体的能量流加做功加化学反应放热，等于内能加动能的累积率。TNT 当量法的骨架就是它。

💡 关键直觉：守恒方程是"记账规则"，不是"剧情描述"。爆炸过程可以千变万化，账必须平——勘测就是找账本上对不上的那一栏。

从方程到勘测计算：破片初速演练

以 1.1 节的储罐爆破为例，演示动量与能量守恒如何落地。储罐破裂抛出 12 块大破片，其中一块 8 千克的法兰盖飞出 120 米。想知道它的初速，不需要模拟爆炸全程，只需要两步：先用抛射运动学反初速，再用能量账核对合理性。

# 破片初速的勘测反演（平抛 + 空气阻力修正的简化版）
import math
m = 8.0            # 破片质量 kg
L = 120.0          # 飞行距离 m
h = 4.0            # 破裂时离地高度 m
g = 9.8
# 不计阻力：飞行时间由高度决定
t = math.sqrt(2*h/g)
v0 = L/t
print(f"不计阻力初速 v0 = {v0:.0f} m/s")
# 空气阻力修正（破片阻力面积 A，阻力系数 Cd，空气密度 rho）
rho, Cd, A = 1.2, 1.0, 0.030
k = rho*Cd*A*L/(2*m)         # 阻力参数
v_corr = v0 * math.exp(k/2)  # 指数阻力模型的近似修正
print(f"计阻力修正后 v0 ≈ {v_corr:.0f} m/s")
# 能量核对：全部动能占膨胀功的比例
W = 1.02e8                     # 1.1 节算出的膨胀功
KE = 0.5*m*v_corr**2
print(f"该破片动能 {KE:.2e} J，占膨胀功 {12*KE/W*100:.1f}%（按12块同量级估）")
​

输出：不计阻力初速约 133 米/秒，修正后约 170 米/秒；该破片动能约 1.2×10^5 焦。若 12 块破片同量级，破片总动能约 1.4×10^6 焦，只占膨胀功的 1.4%——绝大多数能量留在了冲击波里。这个 1%~2% 的破片动能份额是压力容器爆破的经典经验值，能量账本在此完成了一次自洽核对：如果算出的份额离谱（比如 50%），说明初速反演或当量估计有错，回头查。

控制方程组与后续各章的接口

把守恒定律用到连续介质上、再补一条本构关系（状态方程），就得到爆炸力学的控制方程组——欧拉方程组加能量方程。本教程后续所有章节的模型都是它的化身：

⚠️ 常见坑：初学者常把"守恒"误当成"各部分均分"。能量账平了，不代表每一栏都如你直觉所料——破片只分走 1%~2%，冲击波拿走大头，这个分配比由过程决定，不能靠猜。勘测结论必须写明"该份额引自经验值"还是"由测量反演"。

控制方程的勘测边界

守恒框架强大，但要清楚它在哪里开始"失明"。第一条边界是本构关系：守恒定律本身不描述材料，混凝土在动态载荷下裂而未碎的阶段，其应力应变关系强烈依赖损伤历史，此时账本仍平，但把损伤演化算清要额外的损伤方程。第二条边界是多相耦合：含水饱和土壤爆炸后孔隙水压力骤升、有效应力丧失而液化，宏观行为像流体，账本里若仍按两相分别记账再简单相加，误差可达量级。

勘测上对应的纪律是：凡结论涉及这些边界区，报告必须注明所用补充模型及其标定来源。举例来说，某湿砂场地爆破振动预测与实测偏差三倍，复盘发现预测用的是干砂经验参数——不是守恒定律错了，是给定律喂的性格描述错了。把"方程可信"与"参数可信"分开表述，是勘测报告专业度的分水岭。

再补一个实用技巧：动量账可以直接用于快速核对破片分布对称性。把现场回收破片按方位角分组，统计各象限的总动量（质量乘反演初速），若某象限系统性偏大一倍，指向该侧存在附加能量源（如次生爆炸或反射增强），这是账本思维在完全不引入新公式时也能给出的现场判断。

最后补一组常用换算，勘测现场可直接默算：一焦耳等于一牛米；一兆帕等于一牛每平方毫米；重力加速度取九点八，估算破片落地时间用它足够。控制体的选取技巧也值得一句：选控制体时把"想知道的量"放在边界上（比如想知道界面载荷，就把控制体切在界面两侧各取一半），未知量自动出现在守恒方程里，省去中间推导。

本节要点回顾

• 守恒定律是记账规则：不依赖过程细节，用边界测量约束内部过程，特别适合信息残缺的勘测场景；
• 三大守恒的积分形式：质量、动量、能量各对应一条勘测计算主线（药量核对、破片初速、当量评估）；
• 破片动能份额 1%~2%：压力容器爆破的经典经验值，是能量账自洽核对的关键一栏；
• 控制方程组是后续所有模型的公共祖先：第 2 到第 8 章的每张公式脸，都是这张脸的变体。

下一节暂时放下公式，看三百年爆炸技术史——你会发现每一条守恒规则的背后，几乎都站着一个被炸醒的工程师。