冲击波超压峰值曲线绘制
不同距离超压峰值曲线
ω \omega ω是装药质量, r r r为爆炸距离, r ˉ = r / ω 3 \bar{r}=r/\sqrt[3]{\omega} rˉ=r/3ω称为对比距离
几个超压峰值经验公式如下:
我
国
规
范
:
Δ
p
=
0.084
r
ˉ
+
0.27
r
2
ˉ
+
0.7
r
3
ˉ
B
r
o
d
e
公
式
:
Δ
p
=
0.0975
r
ˉ
+
0.01455
r
2
ˉ
+
0.585
r
3
ˉ
−
0.0019
阿
连
绍
夫
公
式
:
Δ
p
=
0.079
r
ˉ
+
0.158
r
2
ˉ
+
0.65
r
3
ˉ
我国规范:\Delta_p=\frac{0.084}{\bar{r}}+\frac{0.27}{\bar{r^2}}+\frac{0.7}{\bar{r^3}}\\ Brode公式:\Delta_p=\frac{0.0975}{\bar{r}}+\frac{0.01455}{\bar{r^2}}+\frac{0.585}{\bar{r^3}}-0.0019\\ 阿连绍夫公式:\Delta_p=\frac{0.079}{\bar{r}}+\frac{0.158}{\bar{r^2}}+\frac{0.65}{\bar{r^3}}\\
我国规范:Δp=rˉ0.084+r2ˉ0.27+r3ˉ0.7Brode公式:Δp=rˉ0.0975+r2ˉ0.01455+r3ˉ0.585−0.0019阿连绍夫公式:Δp=rˉ0.079+r2ˉ0.158+r3ˉ0.65
r
ˉ
⊂
(
1
/
15
,
1
)
\bar{r}\subset(1/15, 1)
rˉ⊂(1/15,1)
bar_r = (1/15 : 0.01: 1);
delta_p1 = 0.084./bar_r+0.27./power(bar_r, 2)+0.7./power(bar_r, 3);
delta_p2 = 0.0975./bar_r+0.01455./power(bar_r, 2)+0.585./power(bar_r, 3)-0.0019;
delta_p3 = 0.079./bar_r+0.158./power(bar_r, 2)+0.65./power(bar_r, 3);
plot(bar_r, delta_p1,'linewidth',1),hold on;
plot(bar_r, delta_p2,'linewidth',1),hold on;
plot(bar_r, delta_p3,'linewidth',1),hold on;
xlabel('r'), ylabel('\Delta_p');
title('冲击波不同距离的超压峰值曲线');
legend('我国标准','Brode公式','阿连绍夫公式');
grid on;