在看各种有爆炸内容的视频时,评论区经常有人对爆炸对建筑物的破坏争论不休,遂产生好奇心,特地翻阅了一些资料,并写出这篇短文,供需要的人学习和估算。


这里贴出mathematica代码,方便函数调用。
Ztnt[R_, W_] := R/W^(1/3); (*返回值为比例距离*)
Henrych[z_] := Piecewise[{
{1.40717/z + 0.55397/z^2 - 0.03572/z^3 + 0.000625/z^4,
0.05 <= z <= 0.3},
{0.61938/z - 0.03262/z^2 + 0.21324/z^3, 0.3 < z < 1},
{0.0662/z + 0.405/z^2 + 0.3288/z^3, 1 <= z <= 10}}];
(*Henrych经验公式*)
keben[z_] := Piecewise[{
{1.07/z^3 - 0.1, z <= 0.9379},
{1.008/z^2.1, 0.9379 < z < 5.3781},
{0.084/z + 0.27/z^2 + 0.7/z^3, 5.3781 <= z <= 10}}]
(*课本推荐的公式*)