# Experimental Study on Peak Pressure of Shock Waves in Quasi-Shallow Water

Authors: Zhenxiong Wang; Wenbin Gu; Jianqing Liu
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Based on the similarity laws of the explosion, this research develops similarity requirements of the small-scale experiments of underwater explosions and establishes a regression model for peak pressure of underwater shock waves under experimental condition. Small-scale experiments are carried out with two types of media at the bottom of the water and for different water depths. The peak pressure of underwater shock waves at different measuring points is acquired. A formula consistent with the similarity law of explosions is obtained and an analysis of the regression precision of the formula confirms its accuracy. Significance experiment indicates that the influence of distance between measuring points and charge on peak pressure of underwater shock wave is the greatest and that of water depth is the least within the range of geometric parameters. An analysis of data from experiments with different media at the bottom of the water reveals an influence on the peak pressure, as the peak pressure of a shock wave in a body of water with a bottom soft mud and rocks is about 1.33 times that of the case where the bottom material is only soft mud.


---

## Body

## 1. Introduction

Underwater shock wave generated by blasting charge in quasi-shallow water is an important reason causing damage to the target. And shock wave’s damage is not only influenced by shape, structure, and material property of target buildings, but also affected by intensity of shock wave which is usually judged by its maximum pressure.Development of experiment technology of underwater explosive load has brought much convenience to research on maximum pressure of underwater shock wave [1], and scholars at home and abroad have conducted numerous experiments and researched on the shock wave generated by the underwater explosion. With abundant experiments and predecessors’ research production, Cole has comprehensively elaborated physical phenomenon and basic law of underwater explosion and established calculating formula on maximum pressure of explosion flow field, specific impulse, and energy density in certain range [2]. Slifko has carried out profound experiment research on pressure characteristics of explosion shock wave in the limitless water mediums, revealing that functional relationships exist among pressure of underwater shock wave, depth, and other factors [3]. Hanasaki and Ito have calculated maximum pressure of underwater pressure wave by measuring the largest particle velocity of water surface when using 3# detonator in underwater explosion [4], but it may not be practical in the actual measure for being influenced by duration and transmission distance of pressure wave. Takayama and Onodera have studied transmission of underwater spherical shock wave through micro explosion induced by shock of ruby laser to plumbum board, quantitatively describing movement and isopycnic line of underwater spherical shock wave [5]; Swisdak Jr. has analyzed the experimental results and discovered that charging radius and density affect the pressure of underwater shock wave [6]; Adushkin et al. have also carried out researches on shock wave and seismic wave of underwater explosion [7]; Ostapenko’s analysis of stable, discontinuous, two-layer waves is of a purely theoretical nature, and the possibility of the existence of such waves requires a special experimental investigation [8]; through analyzing maximum pressure of underwater shock wave from four charging structures’ underwater explosion, Zhou et al. have selected appropriate explosion for underwater arms [9]; Kwang has reversed explosion of the Cheonan through analyzing its underwater shock wave [10]. In all, researches about underwater shock wave are mainly about experiments on maximum pressure of shock wave and pressure features of shock wave when it is under the influence of free surface as well as reversing some parameters of explosions through underwater shock wave.In quasi-shallow water, shock wave from the explosion is influenced by water and bottom interfaces. And the shock wave goes through multiple reflections and transmissions, generating interactions with the wave system and the problem is even more complicated [11]. Main factors influencing maximum pressure are water depth, explosion height, and measuring points’ position (horizontal distance from measuring point to charging center as well as distance from measuring point to water bottom); meanwhile, loading dose, performance parameter, and medium parameter of water bottom largely affect maximum pressure of shock wave. Due to complexity of water bottom factors, researches of pressure features of shock wave impacted by the water bottom are relatively few, and studies about pressure features of shock wave from the quasi-shallow water explosion impacted by both water surface and bottom are fewer [12]. Hence, on account of various factors’ influence on shock wave, a reasonable theoretical calculating model is still extremely difficult to establish.Building regression model of maximum pressure based on similarity theory, the present paper conducted the real explosion simulation experiments with 400 g press-fit TNT explosives. With different water bottom media, maximum pressure of shock wave from measuring points influenced by different factors is collected; then data fitting is adopted to regress calculation formula on maximum pressure of shock wave. Through significance experiment analysis, greater influential factors are displayed and calculating formula of the maximum pressure under different circumstances of the water bottom is demonstrated.


## 2. Determination of Similar Requirements and Similar Conditions of Simulation Experiment

To verify the similar law of underwater shock wave, Cole has conducted a large number of experiments to indicate that all parameters of underwater shock wave due to changes in dose and distance can be estimated accurately, and it proves that similarity law of the underwater explosion is the universal law of shock wave. But due to the special nature of near zone of charge, authentication area of similarity law is generally within the range of water pressure 180 MPa, equivalent to about seven times the area outside the radius of the spherical charge [13, 14]. According to explosion similarity theory, at least three basic similar requirements between prototype experiment and simulation experiment exist, namely, geometric similarity, motion similarity, and dynamic similarity. The three basic requirements are necessary, but not sufficient. Relationship between similar items and corresponding coefficient scale is shown in Table 1.Table 1
Relationship between similar items and coefficient scale.
Similar items
Coefficient scales
Geometric similarity
Length coefficient scaleλ
Motion similarity
Velocity coefficient scaleλ/τ; acceleration coefficient scale λ/τ2
Dynamic similarity
Pressure coefficient scaleπ=ρφ2; energy coefficient scale ε=πλ2
τ refers to time coefficient scale and ρ is density of medium.From the viewpoint of energy similarity, in prototype experiment and simulation experiment, there exist the following relationships between those corresponding parameters in explosive field:(1)Parameter:Pmtm,Xm,Ym,Zm=Pt,X,Y,Z,Velocity:umtm,Xm,Ym,Zm=ut,X,Y,Z,Time:tm=WmW1/3t.


### 2.1. Major Physical Parameters Restricting Explosion Parameters in Quasi-Shallow Water

Without considering water viscosity and gravity, parameters restricting explosion in quasi-shallow water include geometric parameters, water medium, parameters of water bottom medium, explosion parameters, and parameters of underwater shock wave. Restricting factors of shock wave in the quasi-shallow water explosion are in total 17, which are listed together with dimensions in Table2.Table 2
17 Major physical parameters restricting explosion.
Number
Symbol
Name
Dimension
1
E
0
Total energy of charge
ML
2
T
−2
2
r
0
Radius of charge
L
3
d
Depth of water
L
4
H
Distance from charge to water bottom (explosion height)
L
5
L
Horizontal distance from charging center to measuring point
L
6
h
Distance from measuring point to water bottom (measuring height)
L
7
ρ
0
Initial density of water bottom medium
ML
−3
8
σ
Compression strength of water bottom medium
ML
−1
T
−2
9
k
0
Cohesion of water bottom medium
ML
−1
T
−2
10
k
1
Modulus of plastic loading volume of water bottom medium
ML
−1
T
−2
11
p
0
Hydro-static pressure of water
ML
−1
T
−2
12
C
0
Acoustic velocity of water
LT
−1
13
t
Time parameter
T
14
p
Pressure of shock wave
LT
−1
15
D
Velocity of shock wave
ML
−1
T
−2
16
u
Particle velocity rear shock wave
LT
−1
17
ρ
Particle intensity rear shock wave
ML
−3


### 2.2. Similarity Criterion

For 17 main physical parameters necessary to study, the research selects the 3 physical parametersE(ML2T-2), r0(L), and C0(LT-1) to be the basic physical parameters. According to dimension’s harmonious condition as well as π theorem, 14 similarity criteria of independent dimensionless parameters can be listed in 3 categories:Similarity criteria of the 1st category:(2)π1=dr0;π2=Hr0;π3=Lr0;π4=hr0;π10=tr0c0.Similarity criteria of the 2nd category:(3)π5=ρ0E0r03c02;π6=σE0r03;π7=k0E0r03;π8=k1E0r03;π9=p0E0r03;π11=pE0r03;π14=ρE0r03c02.Similarity criteria of the 3rd category:(4)π12=Dc0;π13=uc0.In line with basic theorem of dimension theory, any dimensionless characteristic parameter π in and after explosion process can be above 14 independent dimensionless parameters:(5)π=fπ1,π2,…,π14.


### 2.3. Analysis on Similarity Criterion

The above 14 dimensionless number groups constitute the model law. According to explosion similarity theory, the 14 dimensionless parameters between the model and the real should remain the same. Here, analysis on similarity criterion is shown with subscript “1” standing for the real and subscript “2” for the model:Byπ1: d1/r01=d2/r02⇒d1/d2=r01/r02; the depth of water is similar to charging radius;Byπ2: H1/r01=H2/r02⇒H1/H2=r01/r02; distance from charging center to water bottom is similar to charging radius;Byπ3: L1/r01=L2/r02⇒L1/L2=r01/r02; horizontal distance from charging to measuring point is similar to charging radius;Byπ4: h1/r01=h2/r02⇒h1/h2=r01/r02; distance from charging to water bottom is similar to charging radius;Byπ5: ρ01r013C012/E01=ρ02r023C022/E02; substituting with equation E0=CQv, then ρ01r013C012/C1Qv1=ρ02r023C022/C2Qv2; supposing ρ01Qv2C012/ρ02Qv1C021=r023C1/r013C2=1, then water medium in experiment should be the same; ρ01=ρ02, C01=C02, and with the same explosive Qv1=Qv2;Byπ6: σ1r013/E01=σ2r023/E02⇒σ1r013/C1Qv1=σ2r023/C2Qv2; if C2r013/C1r023=1, the same explosive as well as medium of water bottom is required;Byπ7: k01r013/E01=k02r023/E02⇒k01r013/C1Qv1=k02r023/C2Qv2; provided C1r013/C2r023=1, thus the same explosive and medium of water bottom should be realized;Byπ8: k1r013/E01=k2r023/E02⇒k1r013/C1Qv1=k2r023/C2Qv2; provided C1r013/C2r023=1, then the same explosive and medium of water bottom are required;Byπ9: P01r013/E01=P02r023/E02⇒P01r013/C1Qv1=P02r023/C2Qv2; supposing C1r013/C2r023=1, the same explosive as well as hydrostatic pressure of water bottom is required;Byπ10: t1r013/E01=t2r023/E02⇒t1/t2=r01C02/r02C01; if water media are the same (C01=C02), time parameter is geometrically similar to charging radius;Byπ11: Pm1r013/E01=Pm2r023/E02⇒Pm1r013/C1Qv1=Pm2r023/C2Qv2; if C1r013/C2r023=1, with Pm1=Pm2 required, the same explosive should be adopted;Byπ12: D1/C01=D2/C02, with D1=D2 required; then water media in experiment should be the same;Byπ13: u1/C01=u2/C02, and, with u1=u2 required, water media in experiment should be the same;Byπ14: ρ1r013C012/E01=ρ2r023C022/E02⇒ρ1r013C012/C1Qv1=ρ2r023C022/C2Qv2, on condition that C1r013/C2r023=1, and with ρ1=ρ2 required; therefore, the same explosive and medium of water bottom should be realized.From the above analysis, similarity conditions of simulation experiment can be concluded as(1)
employing the same water medium in experiment (including equal hydrostatic pressure) and the same geological conditions,(2)
employing explosives of the same type and density,(3)
depth of water which is similar to distance from charging to water bottom; horizontal distance from charging center to measuring point is similar to charging radius and distance from measuring point to water bottom is geometrically similar to charging radius,(4)
charging shapes which are geometrically similar and charging dose is directly proportional to the cube of charging radius.If the above similarity conditions for simulation experiment are tenable, that is, the prototype experiment is physically similar to simulation experiment, the experiment results can be possible to achieve and meaningful parameters are the four similarity criteria of the 1st category; then(6)π=fdr0,Hr0,Lr0,tr0c0.


## 2.1. Major Physical Parameters Restricting Explosion Parameters in Quasi-Shallow Water

Without considering water viscosity and gravity, parameters restricting explosion in quasi-shallow water include geometric parameters, water medium, parameters of water bottom medium, explosion parameters, and parameters of underwater shock wave. Restricting factors of shock wave in the quasi-shallow water explosion are in total 17, which are listed together with dimensions in Table2.Table 2
17 Major physical parameters restricting explosion.
Number
Symbol
Name
Dimension
1
E
0
Total energy of charge
ML
2
T
−2
2
r
0
Radius of charge
L
3
d
Depth of water
L
4
H
Distance from charge to water bottom (explosion height)
L
5
L
Horizontal distance from charging center to measuring point
L
6
h
Distance from measuring point to water bottom (measuring height)
L
7
ρ
0
Initial density of water bottom medium
ML
−3
8
σ
Compression strength of water bottom medium
ML
−1
T
−2
9
k
0
Cohesion of water bottom medium
ML
−1
T
−2
10
k
1
Modulus of plastic loading volume of water bottom medium
ML
−1
T
−2
11
p
0
Hydro-static pressure of water
ML
−1
T
−2
12
C
0
Acoustic velocity of water
LT
−1
13
t
Time parameter
T
14
p
Pressure of shock wave
LT
−1
15
D
Velocity of shock wave
ML
−1
T
−2
16
u
Particle velocity rear shock wave
LT
−1
17
ρ
Particle intensity rear shock wave
ML
−3


## 2.2. Similarity Criterion

For 17 main physical parameters necessary to study, the research selects the 3 physical parametersE(ML2T-2), r0(L), and C0(LT-1) to be the basic physical parameters. According to dimension’s harmonious condition as well as π theorem, 14 similarity criteria of independent dimensionless parameters can be listed in 3 categories:Similarity criteria of the 1st category:(2)π1=dr0;π2=Hr0;π3=Lr0;π4=hr0;π10=tr0c0.Similarity criteria of the 2nd category:(3)π5=ρ0E0r03c02;π6=σE0r03;π7=k0E0r03;π8=k1E0r03;π9=p0E0r03;π11=pE0r03;π14=ρE0r03c02.Similarity criteria of the 3rd category:(4)π12=Dc0;π13=uc0.In line with basic theorem of dimension theory, any dimensionless characteristic parameter π in and after explosion process can be above 14 independent dimensionless parameters:(5)π=fπ1,π2,…,π14.


## 2.3. Analysis on Similarity Criterion

The above 14 dimensionless number groups constitute the model law. According to explosion similarity theory, the 14 dimensionless parameters between the model and the real should remain the same. Here, analysis on similarity criterion is shown with subscript “1” standing for the real and subscript “2” for the model:Byπ1: d1/r01=d2/r02⇒d1/d2=r01/r02; the depth of water is similar to charging radius;Byπ2: H1/r01=H2/r02⇒H1/H2=r01/r02; distance from charging center to water bottom is similar to charging radius;Byπ3: L1/r01=L2/r02⇒L1/L2=r01/r02; horizontal distance from charging to measuring point is similar to charging radius;Byπ4: h1/r01=h2/r02⇒h1/h2=r01/r02; distance from charging to water bottom is similar to charging radius;Byπ5: ρ01r013C012/E01=ρ02r023C022/E02; substituting with equation E0=CQv, then ρ01r013C012/C1Qv1=ρ02r023C022/C2Qv2; supposing ρ01Qv2C012/ρ02Qv1C021=r023C1/r013C2=1, then water medium in experiment should be the same; ρ01=ρ02, C01=C02, and with the same explosive Qv1=Qv2;Byπ6: σ1r013/E01=σ2r023/E02⇒σ1r013/C1Qv1=σ2r023/C2Qv2; if C2r013/C1r023=1, the same explosive as well as medium of water bottom is required;Byπ7: k01r013/E01=k02r023/E02⇒k01r013/C1Qv1=k02r023/C2Qv2; provided C1r013/C2r023=1, thus the same explosive and medium of water bottom should be realized;Byπ8: k1r013/E01=k2r023/E02⇒k1r013/C1Qv1=k2r023/C2Qv2; provided C1r013/C2r023=1, then the same explosive and medium of water bottom are required;Byπ9: P01r013/E01=P02r023/E02⇒P01r013/C1Qv1=P02r023/C2Qv2; supposing C1r013/C2r023=1, the same explosive as well as hydrostatic pressure of water bottom is required;Byπ10: t1r013/E01=t2r023/E02⇒t1/t2=r01C02/r02C01; if water media are the same (C01=C02), time parameter is geometrically similar to charging radius;Byπ11: Pm1r013/E01=Pm2r023/E02⇒Pm1r013/C1Qv1=Pm2r023/C2Qv2; if C1r013/C2r023=1, with Pm1=Pm2 required, the same explosive should be adopted;Byπ12: D1/C01=D2/C02, with D1=D2 required; then water media in experiment should be the same;Byπ13: u1/C01=u2/C02, and, with u1=u2 required, water media in experiment should be the same;Byπ14: ρ1r013C012/E01=ρ2r023C022/E02⇒ρ1r013C012/C1Qv1=ρ2r023C022/C2Qv2, on condition that C1r013/C2r023=1, and with ρ1=ρ2 required; therefore, the same explosive and medium of water bottom should be realized.From the above analysis, similarity conditions of simulation experiment can be concluded as(1)
employing the same water medium in experiment (including equal hydrostatic pressure) and the same geological conditions,(2)
employing explosives of the same type and density,(3)
depth of water which is similar to distance from charging to water bottom; horizontal distance from charging center to measuring point is similar to charging radius and distance from measuring point to water bottom is geometrically similar to charging radius,(4)
charging shapes which are geometrically similar and charging dose is directly proportional to the cube of charging radius.If the above similarity conditions for simulation experiment are tenable, that is, the prototype experiment is physically similar to simulation experiment, the experiment results can be possible to achieve and meaningful parameters are the four similarity criteria of the 1st category; then(6)π=fdr0,Hr0,Lr0,tr0c0.


## 3. Establishment of Regression Model of Maximum Pressure and Regression Analysis Method

### 3.1. Fitted Model of Maximum Pressure

On basis of satisfying similar conditions, assuming that group charge explodes in quasi-shallow water, static pressure and density of water arep0 and ρ0. The overpressure of spatial point Δpm=pm-p0 (pm is maximum pressure of incident wave surface; as pm≫p0, Δpm≈pi) is determined by total energy of charging E0, distance from charging to water bottom H, distance from measuring point to water bottom h, depth of water d, and horizontal distance from measuring point to charging L as well as p0 and ρ0, with relative location diagram shown in Figure 1.Figure 1
Relative location diagram of measuring points.And its function expression is(7)Δpm=fH,h,d,L,E0,p0,ρ0.That is to say, (8)pm=fH,h,d,L,E0,p0,ρ0.If E, p0, and ρ0 are selected as independent variables, according to π theorem:(9)π1=pmp0,π2=Hp01/3E01/3,π3=hp01/3E01/3,π4=dp01/3E01/3,π5=Lp01/3E01/3.Hence, relationship among maximum pressure, charging, and measuring points can be demonstrated as(10)Fpmp0,Hp01/3E01/3,hp01/3E01/3dp01/3E01/3,Lp01/3E01/3=0.As E=CQv (C refers to mass of charging and Qv explosion heat), when Qv and p0 are determined, the above equation can be transformed into(11)pm=F1HC3,hC3,dC3,LC3.If H¯=H/C3, h¯=h/C3, d¯=d/C3, and L¯=L/C3 are, respectively, scale shocking height, scale measurement height, scale depth of water, and scale horizontal distance, then F1 can be(12)F1H¯,h¯,d¯,L¯=AH¯b1h¯b2d¯b3L¯b4.Then,(13)pm=AH¯b1h¯b2d¯b3L¯b4.In the equation, coefficient A and indexes b1, b2, b3, and b4 are determined by the experiment.


### 3.2. Regression Method of Experimental Coefficient

Take logarithm of both sides of (13) and nonlinear equation turns into linear equation:(14)ln⁡pm=ln⁡A+b1ln⁡H¯+b2ln⁡h¯+b3ln⁡d¯+b4ln⁡L¯.When X1=ln⁡H¯, X2=ln⁡h¯, X3=ln⁡d¯, X4=ln⁡L¯, X5=Y=ln⁡pm, and b0=ln⁡A, then (14) is(15)Y=b0+b1X1+b2X2+b3X3+b4X4.According to least square method, select b0, b1, b2, b3, and b4 to get residual sum of square:(16)Q=∑j=1nYj-b0+b1X1j+b2X2j+b3X3j+b4X4j2.Here, n is the number of measuring points in maximum pressure statistics; Yj is the jth maximum pressure; X1j, X2j, X3j, and X4j are, respectively, set parameters of Yj corresponding to variables X1, X2, X3, and X4 in (15). If Q has a minimum value, b0, b1, b2, b3, and b4 need to meet(17)∂Q∂bi=0i=0,1,…,4.Equation set can be induced from (17):(18)l11b1+l12b2+l13b3+l14b4=l15,l21b1+l22b2+l23b3+l24b4=l25,l31b1+l32b2+l33b3+l34b4=l35,l41b1+l42b2+l43b3+l44b4=l45.In the equation set,(19)lij=lji=∑XiXj-1n∑Xi∑Xj.And solving coefficient matrix of equation set (18), b0, b1, b2, b3, and b4 can be obtained.


### 3.3. Analysis on Regression Accuracy

In multiple linear regressions, correlative coefficients between variables are complex for correlation may exist between any two variables. When independent variables are highly correlated, it may be really difficult for factor analysis. Correlation coefficient betweenXi and Xj can be obtained by the following equation:(20)rij=∑i=04∑j=1nXij-Xi¯Yj-Y¯bi∑j=1nYj-Y¯2.Here, Xi¯=1/n∑j=1nXiji=0,1,…,4.When the regression coefficient is obtained, regression accuracy should also be calculated and analyzed, including the residual sums of square:(21)U=∑Y^-Y¯2=∑i=0kbilikk=0,1,…,4.The bigger it is, the stronger the regression regularity is, and so are the effects. And the residual standard deviation is(22)s=Qn-k-1.Regression mean square to residual mean square ratio is(23)F=U/kQ/n-k-1=Uks2.Significance experiment can be conducted on the whole regression with value of F, and it would analyze whether linear relationship between Y and k independent variables X1,X2,X3,…,Xk is significant.


## 3.1. Fitted Model of Maximum Pressure

On basis of satisfying similar conditions, assuming that group charge explodes in quasi-shallow water, static pressure and density of water arep0 and ρ0. The overpressure of spatial point Δpm=pm-p0 (pm is maximum pressure of incident wave surface; as pm≫p0, Δpm≈pi) is determined by total energy of charging E0, distance from charging to water bottom H, distance from measuring point to water bottom h, depth of water d, and horizontal distance from measuring point to charging L as well as p0 and ρ0, with relative location diagram shown in Figure 1.Figure 1
Relative location diagram of measuring points.And its function expression is(7)Δpm=fH,h,d,L,E0,p0,ρ0.That is to say, (8)pm=fH,h,d,L,E0,p0,ρ0.If E, p0, and ρ0 are selected as independent variables, according to π theorem:(9)π1=pmp0,π2=Hp01/3E01/3,π3=hp01/3E01/3,π4=dp01/3E01/3,π5=Lp01/3E01/3.Hence, relationship among maximum pressure, charging, and measuring points can be demonstrated as(10)Fpmp0,Hp01/3E01/3,hp01/3E01/3dp01/3E01/3,Lp01/3E01/3=0.As E=CQv (C refers to mass of charging and Qv explosion heat), when Qv and p0 are determined, the above equation can be transformed into(11)pm=F1HC3,hC3,dC3,LC3.If H¯=H/C3, h¯=h/C3, d¯=d/C3, and L¯=L/C3 are, respectively, scale shocking height, scale measurement height, scale depth of water, and scale horizontal distance, then F1 can be(12)F1H¯,h¯,d¯,L¯=AH¯b1h¯b2d¯b3L¯b4.Then,(13)pm=AH¯b1h¯b2d¯b3L¯b4.In the equation, coefficient A and indexes b1, b2, b3, and b4 are determined by the experiment.


## 3.2. Regression Method of Experimental Coefficient

Take logarithm of both sides of (13) and nonlinear equation turns into linear equation:(14)ln⁡pm=ln⁡A+b1ln⁡H¯+b2ln⁡h¯+b3ln⁡d¯+b4ln⁡L¯.When X1=ln⁡H¯, X2=ln⁡h¯, X3=ln⁡d¯, X4=ln⁡L¯, X5=Y=ln⁡pm, and b0=ln⁡A, then (14) is(15)Y=b0+b1X1+b2X2+b3X3+b4X4.According to least square method, select b0, b1, b2, b3, and b4 to get residual sum of square:(16)Q=∑j=1nYj-b0+b1X1j+b2X2j+b3X3j+b4X4j2.Here, n is the number of measuring points in maximum pressure statistics; Yj is the jth maximum pressure; X1j, X2j, X3j, and X4j are, respectively, set parameters of Yj corresponding to variables X1, X2, X3, and X4 in (15). If Q has a minimum value, b0, b1, b2, b3, and b4 need to meet(17)∂Q∂bi=0i=0,1,…,4.Equation set can be induced from (17):(18)l11b1+l12b2+l13b3+l14b4=l15,l21b1+l22b2+l23b3+l24b4=l25,l31b1+l32b2+l33b3+l34b4=l35,l41b1+l42b2+l43b3+l44b4=l45.In the equation set,(19)lij=lji=∑XiXj-1n∑Xi∑Xj.And solving coefficient matrix of equation set (18), b0, b1, b2, b3, and b4 can be obtained.


## 3.3. Analysis on Regression Accuracy

In multiple linear regressions, correlative coefficients between variables are complex for correlation may exist between any two variables. When independent variables are highly correlated, it may be really difficult for factor analysis. Correlation coefficient betweenXi and Xj can be obtained by the following equation:(20)rij=∑i=04∑j=1nXij-Xi¯Yj-Y¯bi∑j=1nYj-Y¯2.Here, Xi¯=1/n∑j=1nXiji=0,1,…,4.When the regression coefficient is obtained, regression accuracy should also be calculated and analyzed, including the residual sums of square:(21)U=∑Y^-Y¯2=∑i=0kbilikk=0,1,…,4.The bigger it is, the stronger the regression regularity is, and so are the effects. And the residual standard deviation is(22)s=Qn-k-1.Regression mean square to residual mean square ratio is(23)F=U/kQ/n-k-1=Uks2.Significance experiment can be conducted on the whole regression with value of F, and it would analyze whether linear relationship between Y and k independent variables X1,X2,X3,…,Xk is significant.


## 4. Simulation and Analysis on Experiment Data

To study influential factors of peak pressure of underwater shock wave, three-dimensional press shot of 400 g is employed in simulation experiment, and settings of blast experiment area and sensor are shown in Figures2 and 3. Through measuring peak pressure of measuring points with different water depths, explosion height, measuring height, and horizontal distance and monitoring measuring points in different underwater medium, peak pressure under different circumstances can be obtained. Export the measured data and process it with MATLAB; then history curve of shock wave pressure is attained and history curve of typical shock wave pressure at measuring points is shown in Figure 4. The data are processed according to regression model of Section 3; afterwards main factors influencing peak pressure at measuring points of the explosion in quasi-shallow water are analyzed. Formulas on peak pressure in quasi-shallow water are obtained through fitting, and then regression accuracy from the equation is analyzed and verified.Figure 2
Experiment area.Figure 3
Layout area of charging blast and sensor.Figure 4
Curve of typical flush water shock wave collected by sensor.


### 4.1. Research on Peak Pressure of Shock Wave of Slime Water Bottom

#### 4.1.1. Measuring Data of Underwater Shock Wave Pressure of Slime Water Bottom

Figure4 shows that, in simulation experiment, MATLAB helps obtain the history curve of typical shock wave pressure at measuring points in water [15], which shows that the change of shock wave pressure is more complex than attenuation of shock wave pressure in the unlimited water explosion. Due to impact of water bottom and surface, there appears obvious attenuation of the shock wave at measuring points, influenced by change of measuring sites. There are many factors influencing shock wave [16], and the present paper mainly studies the influence of water depth, measuring depth, distance, and so on upon peak pressure of shock wave; hence, factors affecting wave change at measuring points would not be analyzed in detail, and peak pressure would be collected and processed.Under the circumstance of slime water bottom, experiment results of peak pressure of part of measuring points in explosion field of simulation experiment are obtained, as is shown in Table3, where p∞ refers to peak pressure of shock wave at underwater measuring points when charge explodes in unlimited water medium with rational distance r calculated by Cole formula [14]:(24)p∞=52.3r¯-1.13106Pa0.637<r¯<12.821.In Table 3, pm is measured pressure of incident wave front of explosion in quasi-shallow water; p∞ is pressure of incident wave front at a distance in the same proportion in infinite water medium with unit of MPa.Table 3
Statistics of incident pressure at part of the measuring points in slime water bottom.
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
0.05
0.01
0.81
1.45
18.0
34.4
0.14
0.26
0.96
1.78
16.6
27.2
0.05
0.01
0.84
1.94
13.3
24.7
0.14
0.65
0.81
1.45
21.5
32.2
0.05
0.26
0.84
1.45
18.5
34.0
0.14
0.65
0.84
1.45
20.3
32.2
0.05
0.26
0.96
1.62
16.4
30.0
0.14
0.65
0.96
1.62
18.7
28.7
0.05
0.52
0.81
1.94
14.5
23.9
0.26
0.01
0.81
1.62
18.0
29.9
0.05
0.65
0.84
1.45
20.7
31.4
0.26
0.01
0.84
1.45
25.0
33.8
0.05
0.65
0.96
1.62
18.8
28.2
0.26
0.26
0.84
1.45
20.3
34.4
0.14
0.01
0.84
1.45
20.7
34.2
0.26
0.26
0.96
1.62
22.3
30.3
0.14
0.01
0.96
2.10
14.1
22.6
0.26
0.65
0.81
1.45
24.6
33.0
0.14
0.26
0.84
1.45
21.9
34.2
0.26
0.65
0.96
1.94
19.5
24.2


#### 4.1.2. Data Fitting and Regression Analysis

By least square method, data in Table1 are analyzed by regression as b0=3.83, b1=0.1507, b2=0.0243, b3=0.0403, and b4=-1.054. And regression equation is(25)Y^=3.83+0.151X1+0.024X2+0.040X3-1.054X4.Regression variables X1, X2, X3, and X4 arrange importance of peak pressure’s expected value Y^, as is in Table 2. Then significance experiment (F test) is carried out. In the table, pi=bi2/cii, cii is the value on diagonal in inverse matrix of coefficient matrix in (18), and Fi=pi/s2.Table4 shows that X4 affects Y most meaningfully, X1 comes second, and X3 does not importantly impact Y. Therefore, the influence of horizontal distance from charging center to measuring point on peak pressure of incident front wave is the most consequential, explosion height comes second, and then comes measuring height; water depth does not have important influence.Table 4
F test.
X
1
X
2
X
3
X
4
p
i
0.20070
0.0862
0.0004
0.2484
F
i
30.5011
13.1025
0.0631
37.7578From the viewpoint of energy distribution, since the equivalent charging radius of 400 g, TNT Charging isr0=4 cm; water depth H=47 cm. Clearly, H>7r0, so energy generated by the explosion is mainly used for forming shock wave and underwater pits of explosion and the energy does not overflow water or form air shock wave. And it is consistent with regression conclusion “water depth does not significantly influence peak pressure of incident.” Hence, in explosion experiment of quasi-shallow water, water depth can be ignored in researching peak value of measuring points within a flow field [1]. After removing items including X3 in (25), experiment data are reanalyzed by regression as b0=3.8225, b1=0.1515, b2=0.0246, and b4=-1.0462. Thus, regression equation of the experiment results is (26)Y^=3.823+0.152X1+0.025X2-1.046X4.Calculating results of simple correlation coefficients among regression variables are in Table 5. As can be seen from the table, variables have little correlation between each other. Replace Y^(X5) in formula (26) with ln⁡pm and X1, X2, and X4 with ln⁡H¯, ln⁡h¯, and ln⁡L¯, respectively; then, in experiment of soft mud, when single charge explodes in quasi-shallow water, calculating formula of shock wave pressure of measuring points in flow field is (27)pm=45.6H¯0.152h¯0.025L¯-1.046.Fitting range of (27) is 0.0529≤H¯≤0.258; 0.0068≤h¯≤0.244; 0.814≤d¯≤0.964; 1.45≤L¯≤2.10.Table 5
Calculation of simple correlation coefficients among regression variables.
r
i
j
X
1
X
2
X
4
X
5
=
Y
X
1
1.0000
0.0013
-
0.0803
0.5430
X
2
0.0013
1.0000
-
0.0000
0.2276
X
4
-
0.0803
-
0.0000
1.0000
-
0.7617
X
5
=
Y
0.5430
0.2276
-
0.7617
1.0000Calculate peak pressure of some measuring points with (27) peak pressure shock wave equation (27) and then compare results with actually measured peak pressure and the details are shown in Table 6, where the unit of measured and calculated pm is MPa.Table 6
Comparisons of calculating results by regression and measured results.
H
¯
h
¯
L
¯
A-pm
C-pm
R
H
¯
h
¯
L
¯
A-pm
C-pm
R
0.05
0.0068
1.45
17.91
17.42
2.74
0.14
0.65
1.45
22.48
22.59
0.49
0.05
0.0068
1.94
12.34
12.95
4.94
0.14
0.65
1.62
19.48
20.27
4.06
0.05
0.26
1.45
18.51
19.11
3.24
0.14
0.65
1.94
17.90
16.8
6.15
0.05
0.26
1.62
16.45
17.14
4.20
0.14
0.65
2.10
16.40
15.47
5.67
0.05
0.26
1.78
16.55
15.54
6.10
0.26
0.0068
1.45
21.95
22.18
1.05
0.05
0.26
2.10
14.44
13.08
9.42
0.26
0.0068
1.62
18.07
19.90
10.13
0.05
0.52
1.94
14.63
14.46
1.16
0.26
0.0068
1.78
19.47
18.04
7.34
0.05
0.52
2.10
12.63
13.32
5.46
0.26
0.0068
2.10
15.2
15.19
0.66
0.05
0.65
1.45
20.20
19.57
3.12
0.26
0.26
1.45
23.88
24.34
1.93
0.05
0.65
1.62
18.05
17.55
2.77
0.26
0.26
1.62
22.33
21.83
2.24
0.05
0.65
2.10
13.09
13.40
2.37
0.26
0.26
1.78
21.02
19.79
5.85
0.14
0.0068
1.45
20.86
20.11
3.60
0.26
0.26
1.94
18.52
18.09
2.32
0.14
0.0068
1.94
12.93
14.95
15.62
0.26
0.65
1.45
25.63
24.92
2.77
0.14
0.0068
2.10
14.13
13.77
5.48
0.26
0.65
1.62
22.77
22.35
1.84
0.14
0.26
1.45
21.98
22.06
0.36
0.26
0.65
1.94
19.55
18.53
5.22
0.14
0.26
1.78
18.72
17.94
4.17
0.26
0.65
2.10
15.59
17.06
9.43
0.14
0.26
1.94
16.61
16.40
1.26
A
-
p
m: actually measured pm; C-pm: calculated pm; R: relative error (%).Analysis in Table6 shows that relative error between calculated peak pressure at 33 measuring points and corresponding actual data is between 0.49% and 15.62%, and there is only one point where the relative error is beyond 15%, and it should be regarded as singular experiment data. Therefore, singular results in measurement can be eliminated through comparing back-stepped peak pressure calculated by regression formula with actually measured pressure and analyzing experiment wave form. There are seven measuring points where the relative error is between 5% and 10%, accounting for 20%, and those of the remaining points are smaller than 5%. In a word, regression formula can be employed in calculating peak pressure of the shock wave of the single charge explosion in the quasi-shallow water flow field. And accuracy of the calculating results is relatively high in formula scope.


### 4.2. Studies on Peak Pressure of Shock Wave in Water Bottom of Soft Mud with Rocks

#### 4.2.1. Peak Pressure of Shock Wave in Water Bottom of Soft Mud with Rocks

Under the condition at the water bottom of soft mud with rocks, in model experiment of actual explosion, actual results of experimented peak pressure of some measuring points within the explosion flow field are listed in Table7.Table 7
Incident wave pressure of some measuring points under the condition of water bottom of soft mud with rocks.
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
0.05
0.01
0.81
1.45
26.6
34.4
0.14
0.65
0.84
1.45
25.8
32.2
0.05
0.01
0.81
1.62
21.1
30.3
0.14
0.65
0.81
1.62
25.0
28.7
0.05
0.01
0.81
1.78
18.8
27.3
0.14
0.65
0.81
1.78
22.3
26.1
0.05
0.01
0.96
2.10
15.2
22.6
0.14
0.65
0.81
1.94
22.7
23.8
0.05
0.26
0.84
1.45
24.8
34.0
0.14
0.65
0.81
2.10
20.3
21.9
0.05
0.26
0.96
1.62
22.7
30.0
0.26
0.01
0.81
1.45
28.9
33.8
0.05
0.26
0.96
1.78
21.7
27.0
0.26
0.01
0.81
1.94
20.3
24.5
0.05
0.26
0.96
1.94
20.6
24.6
0.26
0.01
0.96
1.62
28.9
29.9
0.05
0.26
0.81
2.10
20.3
22.5
0.26
0.26
0.84
1.45
29.7
34.4
0.05
0.65
0.81
1.94
23.8
23.9
0.26
0.26
0.96
1.62
26.6
30.3
0.05
0.65
0.96
1.78
23.1
25.7
0.26
0.65
0.81
1.94
21.5
24.2
0.14
0.01
0.81
1.45
24.4
34.2
0.26
0.65
0.96
1.62
27.3
29.4
0.14
0.01
0.81
1.62
21.4
30.2
0.26
0.65
0.96
1.78
24.6
26.5
0.14
0.26
0.81
2.10
19.5
22.6
0.26
0.65
0.96
2.10
19.9
21.3


#### 4.2.2. Data Fitting and Regression Analysis

The least square method can be adopted to calculate data in Table7 on regression. Through F Experiment and analysis, it is shown that water depth X3 does not significantly influence pressure regression parameter X5. Thus, under the experiment condition, water depth could be ignored in regression of peak pressure of the shock wave. Experiment data of regression are b0=3.83 MPa, b1=0.097, b2=0.017, and b4=-0.903. Regression equation is(28)Y^=3.83+0.097X1+0.017X2-0.903X4.Therefore, under the condition at the water bottom of soft mud with rocks, in the explosion of single charge in quasi-shallow water, calculating formula of peak pressure of the shock wave at measuring points within the flow field is(29)pm=48.4H¯0.097h¯0.017L¯-0.903.In the equation, fitting range is 0.0529≤H¯≤0.258; 0.0068≤h¯≤0.244; 0.814≤d¯≤0.964; 1.45≤L¯≤2.10.When (29) is employed in calculating explosion in the water bottom of mud with rocks or in quasi-shallow water, actually measured peak pressure and peak pressure of the shock wave at some measuring points are shown in Table 8.Table 8
Actually measured peak pressure and calculated peak pressure of shock wave in water bottom of mud with rocks.
H
¯
h
¯
L
¯
A-pm
C-pm
R
H
¯
h
¯
L
¯
A-pm
C-pm
R
0.05
0.0068
1.45
26.61
23.93
10.07
0.14
0.0068
2.10
17.08
18.24
6.79
0.05
0.0068
1.62
21.12
21.55
2.00
0.14
0.26
1.62
26.08
25.28
3.07
0.05
0.0068
1.78
19.04
19.61
2.99
0.14
0.26
1.62
20.85
21.10
1.20
0.05
0.0068
1.94
17.98
17.99
0.06
0.14
0.26
2.10
19.48
19.49
0.05
0.05
0.26
1.45
24.90
25.57
2.69
0.14
0.65
1.62
25.04
25.71
2.67
0.05
0.26
1.62
22.69
23.04
1.54
0.14
0.65
1.78
22.30
23.39
4.89
0.05
0.26
1.78
21.78
20.96
3.76
0.14
0.65
1.94
22.27
21.46
3.64
0.05
0.26
1.94
20.64
19.23
6.83
0.14
0.65
2.10
20.36
19.82
2.65
0.05
0.52
1.94
19.52
19.47
0.26
0.26
0.0068
1.62
25.96
25.19
2.97
0.05
0.52
2.10
17.66
17.99
1.87
0.26
0.0068
1.78
22.82
22.92
0.44
0.05
0.65
1.45
24.18
26.01
7.57
0.26
0.0068
1.94
20.65
21.03
1.84
0.05
0.65
1.62
22.59
23.43
3.72
0.26
0.0068
2.10
19.26
19.43
0.88
0.05
0.65
1.78
22.22
21.32
4.05
0.26
0.26
1.45
30.68
29.89
2.58
0.05
0.65
1.94
20.53
19.56
4.72
0.26
0.26
1.62
26.60
26.93
1.24
0.14
0.0068
1.45
27.13
26.25
3.24
0.26
0.26
2.10
21.31
20.76
2.58
0.14
0.0068
1.62
23.05
23.65
2.60
0.26
0.65
1.62
27.33
27.39
0.22
0.14
0.0068
1.78
22.05
21.52
2.40
0.26
0.65
1.78
25.41
24.92
1.93
0.14
0.0068
1.94
19.70
19.74
0.20
0.26
0.65
1.94
23.02
22.86
0.70
A
-
p
m: actually measured pm; C-pm: calculated pm; R: relative error (%).Comparison on measured value and calculated value from the regression formula at the same measuring points displays that comparative error is between 0.06% and 10.07%. There are 32 measuring points with proportional error smaller than 5% and only one point larger than 10%. Consequently, regression formula can be employed in calculating peak pressure of the shock wave of a single charge explosion in the quasi-shallow water flow field. And accuracy of the calculating results is relatively high in formula scope.


### 4.3. Analysis on Influence of Water Bottom Medium on Peak Pressure of Shock Wave

Figure5 demonstrates change curve where theoretical and experimental values change with proportional distance. With the water bottom of mud with rocks, peak pressure of measuring points distributes near the theoretical values. Comparison of two pictures shows that peak pressure at measuring points on the water bottom of mud with rocks is higher than that of the water bottom of soft mud, which is caused by reflective wave to the water bottom and surface cut off pressure of shock wave and different interactions. To study the influence of different media of the water bottom on peak pressure, experiments are carried out in two different media of water bottom and peak pressure at the same measuring points which are shown in Table 9.Table 9
Comparison of peak pressure in two different media of water bottom (unit: MPa).
H
¯
h
¯
L
¯
r
¯
p
m
Water bottom of soft mudp1m
Water bottom of soft mud with rocksp2m
0.053
0.0068
1.4522
1.4529
17.97
26.56
1.9408
1.9410
12.30
17.97
0.2579
1.4522
1.4666
18.51
24.80
1.6151
1.6280
16.41
22.66
1.7779
1.7897
16.41
21.70
2.1037
2.1137
14.45
20.31
0.5157
1.9408
1.9952
14.45
23.83
2.1037
2.1540
12.50
20.31
0.6515
1.4522
1.5707
20.70
23.83
1.6151
1.7224
18.36
22.36
0.136
0.0068
1.4522
1.4579
20.70
27.16
1.9408
1.9451
12.89
19.66
2.1037
2.1077
14.06
17.10
0.6515
1.4522
1.5410
20.69
25.78
1.9408
2.0081
17.97
22.27
2.1037
2.1659
16.41
20.31
0.258
0.0068
1.4522
1.4738
25.78
29.30
1.6151
1.6345
17.97
28.91
2.1037
2.1186
15.23
21.88
0.2579
1.4522
1.4522
20.31
30.66
1.6151
1.6151
22.31
26.56
1.6151
1.6151
22.79
27.34
1.9408
1.9408
19.53
23.05
2.1037
2.1037
15.63
20.51Figure 5
Change curve of peak pressure of shock wave’s change with proportional distance in different media of water bottom.
(a)
Water bottom of soft mud
(b)
Water bottom of mud with rocksStatistical results of Table9 display that peak pressure at the water bottom of soft mud with rocks is greater than that of soft mud, indicating that the harder the medium at the water bottom is, the greater the explosion load of charging in quasi-shallow water is, and the better the destruction of the explosion can be improved. Based on data in Table 9 and numerous experimental data, relationship betweenp1m and p2m can be obtained by linear regression:(30)p2m=1.33p1m.Under the same experiment conditions in quasi-shallow water explosion, peak pressure of the shock wave in the water bottom of mud with rocks is 1.33 times that of the water bottom of soft mud,(31)r¯=H¯-h¯2+L¯2.Fitting curve of experimentally measured ln⁡pm of different water bottom and ln⁡r¯ can be shown in Figure 6 and responding functional relations arewith water bottom of soft mud(32)ln⁡pm1=3.42-0.99ln⁡r¯,with water bottom of soft mud with rocks(33)ln⁡pm2=3.63-0.85ln⁡r¯.And exponential forms are (34)pm1=30.55r¯-0.99MPa,(35)pm2=37.84r¯-0.85MPa.Analysis on experimentally measured values and results of regression reveals the following:(1)
In (27) and (29) of two different water bottoms, peak pressure of shock wave increases as the increase of explosion height and measuring height, and the influence of explosion height is greater than that of measuring height; besides, peak pressure of measuring points decreases with the increase of horizontal distance, and the influence is greater than explosion height and measuring height; meanwhile, peak pressure at measuring points in water is smaller than that at the same position in the infinite water explosion.(2)
Comparison of (33) and (34) indicates that, in the explosion of quasi-shallow water at the water bottom of soft mud, peak pressure of the shock wave attenuates faster than that in the water bottom of mud with rocks, and amplitude of peak pressure at the water bottom of soft mud is smaller than that of the water bottom of soft mud with rocks.Figure 6
Fitting curve of peak pressure of shock wave with two media of water bottom.


## 4.1. Research on Peak Pressure of Shock Wave of Slime Water Bottom

### 4.1.1. Measuring Data of Underwater Shock Wave Pressure of Slime Water Bottom

Figure4 shows that, in simulation experiment, MATLAB helps obtain the history curve of typical shock wave pressure at measuring points in water [15], which shows that the change of shock wave pressure is more complex than attenuation of shock wave pressure in the unlimited water explosion. Due to impact of water bottom and surface, there appears obvious attenuation of the shock wave at measuring points, influenced by change of measuring sites. There are many factors influencing shock wave [16], and the present paper mainly studies the influence of water depth, measuring depth, distance, and so on upon peak pressure of shock wave; hence, factors affecting wave change at measuring points would not be analyzed in detail, and peak pressure would be collected and processed.Under the circumstance of slime water bottom, experiment results of peak pressure of part of measuring points in explosion field of simulation experiment are obtained, as is shown in Table3, where p∞ refers to peak pressure of shock wave at underwater measuring points when charge explodes in unlimited water medium with rational distance r calculated by Cole formula [14]:(24)p∞=52.3r¯-1.13106Pa0.637<r¯<12.821.In Table 3, pm is measured pressure of incident wave front of explosion in quasi-shallow water; p∞ is pressure of incident wave front at a distance in the same proportion in infinite water medium with unit of MPa.Table 3
Statistics of incident pressure at part of the measuring points in slime water bottom.
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
0.05
0.01
0.81
1.45
18.0
34.4
0.14
0.26
0.96
1.78
16.6
27.2
0.05
0.01
0.84
1.94
13.3
24.7
0.14
0.65
0.81
1.45
21.5
32.2
0.05
0.26
0.84
1.45
18.5
34.0
0.14
0.65
0.84
1.45
20.3
32.2
0.05
0.26
0.96
1.62
16.4
30.0
0.14
0.65
0.96
1.62
18.7
28.7
0.05
0.52
0.81
1.94
14.5
23.9
0.26
0.01
0.81
1.62
18.0
29.9
0.05
0.65
0.84
1.45
20.7
31.4
0.26
0.01
0.84
1.45
25.0
33.8
0.05
0.65
0.96
1.62
18.8
28.2
0.26
0.26
0.84
1.45
20.3
34.4
0.14
0.01
0.84
1.45
20.7
34.2
0.26
0.26
0.96
1.62
22.3
30.3
0.14
0.01
0.96
2.10
14.1
22.6
0.26
0.65
0.81
1.45
24.6
33.0
0.14
0.26
0.84
1.45
21.9
34.2
0.26
0.65
0.96
1.94
19.5
24.2


### 4.1.2. Data Fitting and Regression Analysis

By least square method, data in Table1 are analyzed by regression as b0=3.83, b1=0.1507, b2=0.0243, b3=0.0403, and b4=-1.054. And regression equation is(25)Y^=3.83+0.151X1+0.024X2+0.040X3-1.054X4.Regression variables X1, X2, X3, and X4 arrange importance of peak pressure’s expected value Y^, as is in Table 2. Then significance experiment (F test) is carried out. In the table, pi=bi2/cii, cii is the value on diagonal in inverse matrix of coefficient matrix in (18), and Fi=pi/s2.Table4 shows that X4 affects Y most meaningfully, X1 comes second, and X3 does not importantly impact Y. Therefore, the influence of horizontal distance from charging center to measuring point on peak pressure of incident front wave is the most consequential, explosion height comes second, and then comes measuring height; water depth does not have important influence.Table 4
F test.
X
1
X
2
X
3
X
4
p
i
0.20070
0.0862
0.0004
0.2484
F
i
30.5011
13.1025
0.0631
37.7578From the viewpoint of energy distribution, since the equivalent charging radius of 400 g, TNT Charging isr0=4 cm; water depth H=47 cm. Clearly, H>7r0, so energy generated by the explosion is mainly used for forming shock wave and underwater pits of explosion and the energy does not overflow water or form air shock wave. And it is consistent with regression conclusion “water depth does not significantly influence peak pressure of incident.” Hence, in explosion experiment of quasi-shallow water, water depth can be ignored in researching peak value of measuring points within a flow field [1]. After removing items including X3 in (25), experiment data are reanalyzed by regression as b0=3.8225, b1=0.1515, b2=0.0246, and b4=-1.0462. Thus, regression equation of the experiment results is (26)Y^=3.823+0.152X1+0.025X2-1.046X4.Calculating results of simple correlation coefficients among regression variables are in Table 5. As can be seen from the table, variables have little correlation between each other. Replace Y^(X5) in formula (26) with ln⁡pm and X1, X2, and X4 with ln⁡H¯, ln⁡h¯, and ln⁡L¯, respectively; then, in experiment of soft mud, when single charge explodes in quasi-shallow water, calculating formula of shock wave pressure of measuring points in flow field is (27)pm=45.6H¯0.152h¯0.025L¯-1.046.Fitting range of (27) is 0.0529≤H¯≤0.258; 0.0068≤h¯≤0.244; 0.814≤d¯≤0.964; 1.45≤L¯≤2.10.Table 5
Calculation of simple correlation coefficients among regression variables.
r
i
j
X
1
X
2
X
4
X
5
=
Y
X
1
1.0000
0.0013
-
0.0803
0.5430
X
2
0.0013
1.0000
-
0.0000
0.2276
X
4
-
0.0803
-
0.0000
1.0000
-
0.7617
X
5
=
Y
0.5430
0.2276
-
0.7617
1.0000Calculate peak pressure of some measuring points with (27) peak pressure shock wave equation (27) and then compare results with actually measured peak pressure and the details are shown in Table 6, where the unit of measured and calculated pm is MPa.Table 6
Comparisons of calculating results by regression and measured results.
H
¯
h
¯
L
¯
A-pm
C-pm
R
H
¯
h
¯
L
¯
A-pm
C-pm
R
0.05
0.0068
1.45
17.91
17.42
2.74
0.14
0.65
1.45
22.48
22.59
0.49
0.05
0.0068
1.94
12.34
12.95
4.94
0.14
0.65
1.62
19.48
20.27
4.06
0.05
0.26
1.45
18.51
19.11
3.24
0.14
0.65
1.94
17.90
16.8
6.15
0.05
0.26
1.62
16.45
17.14
4.20
0.14
0.65
2.10
16.40
15.47
5.67
0.05
0.26
1.78
16.55
15.54
6.10
0.26
0.0068
1.45
21.95
22.18
1.05
0.05
0.26
2.10
14.44
13.08
9.42
0.26
0.0068
1.62
18.07
19.90
10.13
0.05
0.52
1.94
14.63
14.46
1.16
0.26
0.0068
1.78
19.47
18.04
7.34
0.05
0.52
2.10
12.63
13.32
5.46
0.26
0.0068
2.10
15.2
15.19
0.66
0.05
0.65
1.45
20.20
19.57
3.12
0.26
0.26
1.45
23.88
24.34
1.93
0.05
0.65
1.62
18.05
17.55
2.77
0.26
0.26
1.62
22.33
21.83
2.24
0.05
0.65
2.10
13.09
13.40
2.37
0.26
0.26
1.78
21.02
19.79
5.85
0.14
0.0068
1.45
20.86
20.11
3.60
0.26
0.26
1.94
18.52
18.09
2.32
0.14
0.0068
1.94
12.93
14.95
15.62
0.26
0.65
1.45
25.63
24.92
2.77
0.14
0.0068
2.10
14.13
13.77
5.48
0.26
0.65
1.62
22.77
22.35
1.84
0.14
0.26
1.45
21.98
22.06
0.36
0.26
0.65
1.94
19.55
18.53
5.22
0.14
0.26
1.78
18.72
17.94
4.17
0.26
0.65
2.10
15.59
17.06
9.43
0.14
0.26
1.94
16.61
16.40
1.26
A
-
p
m: actually measured pm; C-pm: calculated pm; R: relative error (%).Analysis in Table6 shows that relative error between calculated peak pressure at 33 measuring points and corresponding actual data is between 0.49% and 15.62%, and there is only one point where the relative error is beyond 15%, and it should be regarded as singular experiment data. Therefore, singular results in measurement can be eliminated through comparing back-stepped peak pressure calculated by regression formula with actually measured pressure and analyzing experiment wave form. There are seven measuring points where the relative error is between 5% and 10%, accounting for 20%, and those of the remaining points are smaller than 5%. In a word, regression formula can be employed in calculating peak pressure of the shock wave of the single charge explosion in the quasi-shallow water flow field. And accuracy of the calculating results is relatively high in formula scope.


## 4.1.1. Measuring Data of Underwater Shock Wave Pressure of Slime Water Bottom

Figure4 shows that, in simulation experiment, MATLAB helps obtain the history curve of typical shock wave pressure at measuring points in water [15], which shows that the change of shock wave pressure is more complex than attenuation of shock wave pressure in the unlimited water explosion. Due to impact of water bottom and surface, there appears obvious attenuation of the shock wave at measuring points, influenced by change of measuring sites. There are many factors influencing shock wave [16], and the present paper mainly studies the influence of water depth, measuring depth, distance, and so on upon peak pressure of shock wave; hence, factors affecting wave change at measuring points would not be analyzed in detail, and peak pressure would be collected and processed.Under the circumstance of slime water bottom, experiment results of peak pressure of part of measuring points in explosion field of simulation experiment are obtained, as is shown in Table3, where p∞ refers to peak pressure of shock wave at underwater measuring points when charge explodes in unlimited water medium with rational distance r calculated by Cole formula [14]:(24)p∞=52.3r¯-1.13106Pa0.637<r¯<12.821.In Table 3, pm is measured pressure of incident wave front of explosion in quasi-shallow water; p∞ is pressure of incident wave front at a distance in the same proportion in infinite water medium with unit of MPa.Table 3
Statistics of incident pressure at part of the measuring points in slime water bottom.
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
0.05
0.01
0.81
1.45
18.0
34.4
0.14
0.26
0.96
1.78
16.6
27.2
0.05
0.01
0.84
1.94
13.3
24.7
0.14
0.65
0.81
1.45
21.5
32.2
0.05
0.26
0.84
1.45
18.5
34.0
0.14
0.65
0.84
1.45
20.3
32.2
0.05
0.26
0.96
1.62
16.4
30.0
0.14
0.65
0.96
1.62
18.7
28.7
0.05
0.52
0.81
1.94
14.5
23.9
0.26
0.01
0.81
1.62
18.0
29.9
0.05
0.65
0.84
1.45
20.7
31.4
0.26
0.01
0.84
1.45
25.0
33.8
0.05
0.65
0.96
1.62
18.8
28.2
0.26
0.26
0.84
1.45
20.3
34.4
0.14
0.01
0.84
1.45
20.7
34.2
0.26
0.26
0.96
1.62
22.3
30.3
0.14
0.01
0.96
2.10
14.1
22.6
0.26
0.65
0.81
1.45
24.6
33.0
0.14
0.26
0.84
1.45
21.9
34.2
0.26
0.65
0.96
1.94
19.5
24.2


## 4.1.2. Data Fitting and Regression Analysis

By least square method, data in Table1 are analyzed by regression as b0=3.83, b1=0.1507, b2=0.0243, b3=0.0403, and b4=-1.054. And regression equation is(25)Y^=3.83+0.151X1+0.024X2+0.040X3-1.054X4.Regression variables X1, X2, X3, and X4 arrange importance of peak pressure’s expected value Y^, as is in Table 2. Then significance experiment (F test) is carried out. In the table, pi=bi2/cii, cii is the value on diagonal in inverse matrix of coefficient matrix in (18), and Fi=pi/s2.Table4 shows that X4 affects Y most meaningfully, X1 comes second, and X3 does not importantly impact Y. Therefore, the influence of horizontal distance from charging center to measuring point on peak pressure of incident front wave is the most consequential, explosion height comes second, and then comes measuring height; water depth does not have important influence.Table 4
F test.
X
1
X
2
X
3
X
4
p
i
0.20070
0.0862
0.0004
0.2484
F
i
30.5011
13.1025
0.0631
37.7578From the viewpoint of energy distribution, since the equivalent charging radius of 400 g, TNT Charging isr0=4 cm; water depth H=47 cm. Clearly, H>7r0, so energy generated by the explosion is mainly used for forming shock wave and underwater pits of explosion and the energy does not overflow water or form air shock wave. And it is consistent with regression conclusion “water depth does not significantly influence peak pressure of incident.” Hence, in explosion experiment of quasi-shallow water, water depth can be ignored in researching peak value of measuring points within a flow field [1]. After removing items including X3 in (25), experiment data are reanalyzed by regression as b0=3.8225, b1=0.1515, b2=0.0246, and b4=-1.0462. Thus, regression equation of the experiment results is (26)Y^=3.823+0.152X1+0.025X2-1.046X4.Calculating results of simple correlation coefficients among regression variables are in Table 5. As can be seen from the table, variables have little correlation between each other. Replace Y^(X5) in formula (26) with ln⁡pm and X1, X2, and X4 with ln⁡H¯, ln⁡h¯, and ln⁡L¯, respectively; then, in experiment of soft mud, when single charge explodes in quasi-shallow water, calculating formula of shock wave pressure of measuring points in flow field is (27)pm=45.6H¯0.152h¯0.025L¯-1.046.Fitting range of (27) is 0.0529≤H¯≤0.258; 0.0068≤h¯≤0.244; 0.814≤d¯≤0.964; 1.45≤L¯≤2.10.Table 5
Calculation of simple correlation coefficients among regression variables.
r
i
j
X
1
X
2
X
4
X
5
=
Y
X
1
1.0000
0.0013
-
0.0803
0.5430
X
2
0.0013
1.0000
-
0.0000
0.2276
X
4
-
0.0803
-
0.0000
1.0000
-
0.7617
X
5
=
Y
0.5430
0.2276
-
0.7617
1.0000Calculate peak pressure of some measuring points with (27) peak pressure shock wave equation (27) and then compare results with actually measured peak pressure and the details are shown in Table 6, where the unit of measured and calculated pm is MPa.Table 6
Comparisons of calculating results by regression and measured results.
H
¯
h
¯
L
¯
A-pm
C-pm
R
H
¯
h
¯
L
¯
A-pm
C-pm
R
0.05
0.0068
1.45
17.91
17.42
2.74
0.14
0.65
1.45
22.48
22.59
0.49
0.05
0.0068
1.94
12.34
12.95
4.94
0.14
0.65
1.62
19.48
20.27
4.06
0.05
0.26
1.45
18.51
19.11
3.24
0.14
0.65
1.94
17.90
16.8
6.15
0.05
0.26
1.62
16.45
17.14
4.20
0.14
0.65
2.10
16.40
15.47
5.67
0.05
0.26
1.78
16.55
15.54
6.10
0.26
0.0068
1.45
21.95
22.18
1.05
0.05
0.26
2.10
14.44
13.08
9.42
0.26
0.0068
1.62
18.07
19.90
10.13
0.05
0.52
1.94
14.63
14.46
1.16
0.26
0.0068
1.78
19.47
18.04
7.34
0.05
0.52
2.10
12.63
13.32
5.46
0.26
0.0068
2.10
15.2
15.19
0.66
0.05
0.65
1.45
20.20
19.57
3.12
0.26
0.26
1.45
23.88
24.34
1.93
0.05
0.65
1.62
18.05
17.55
2.77
0.26
0.26
1.62
22.33
21.83
2.24
0.05
0.65
2.10
13.09
13.40
2.37
0.26
0.26
1.78
21.02
19.79
5.85
0.14
0.0068
1.45
20.86
20.11
3.60
0.26
0.26
1.94
18.52
18.09
2.32
0.14
0.0068
1.94
12.93
14.95
15.62
0.26
0.65
1.45
25.63
24.92
2.77
0.14
0.0068
2.10
14.13
13.77
5.48
0.26
0.65
1.62
22.77
22.35
1.84
0.14
0.26
1.45
21.98
22.06
0.36
0.26
0.65
1.94
19.55
18.53
5.22
0.14
0.26
1.78
18.72
17.94
4.17
0.26
0.65
2.10
15.59
17.06
9.43
0.14
0.26
1.94
16.61
16.40
1.26
A
-
p
m: actually measured pm; C-pm: calculated pm; R: relative error (%).Analysis in Table6 shows that relative error between calculated peak pressure at 33 measuring points and corresponding actual data is between 0.49% and 15.62%, and there is only one point where the relative error is beyond 15%, and it should be regarded as singular experiment data. Therefore, singular results in measurement can be eliminated through comparing back-stepped peak pressure calculated by regression formula with actually measured pressure and analyzing experiment wave form. There are seven measuring points where the relative error is between 5% and 10%, accounting for 20%, and those of the remaining points are smaller than 5%. In a word, regression formula can be employed in calculating peak pressure of the shock wave of the single charge explosion in the quasi-shallow water flow field. And accuracy of the calculating results is relatively high in formula scope.


## 4.2. Studies on Peak Pressure of Shock Wave in Water Bottom of Soft Mud with Rocks

### 4.2.1. Peak Pressure of Shock Wave in Water Bottom of Soft Mud with Rocks

Under the condition at the water bottom of soft mud with rocks, in model experiment of actual explosion, actual results of experimented peak pressure of some measuring points within the explosion flow field are listed in Table7.Table 7
Incident wave pressure of some measuring points under the condition of water bottom of soft mud with rocks.
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
0.05
0.01
0.81
1.45
26.6
34.4
0.14
0.65
0.84
1.45
25.8
32.2
0.05
0.01
0.81
1.62
21.1
30.3
0.14
0.65
0.81
1.62
25.0
28.7
0.05
0.01
0.81
1.78
18.8
27.3
0.14
0.65
0.81
1.78
22.3
26.1
0.05
0.01
0.96
2.10
15.2
22.6
0.14
0.65
0.81
1.94
22.7
23.8
0.05
0.26
0.84
1.45
24.8
34.0
0.14
0.65
0.81
2.10
20.3
21.9
0.05
0.26
0.96
1.62
22.7
30.0
0.26
0.01
0.81
1.45
28.9
33.8
0.05
0.26
0.96
1.78
21.7
27.0
0.26
0.01
0.81
1.94
20.3
24.5
0.05
0.26
0.96
1.94
20.6
24.6
0.26
0.01
0.96
1.62
28.9
29.9
0.05
0.26
0.81
2.10
20.3
22.5
0.26
0.26
0.84
1.45
29.7
34.4
0.05
0.65
0.81
1.94
23.8
23.9
0.26
0.26
0.96
1.62
26.6
30.3
0.05
0.65
0.96
1.78
23.1
25.7
0.26
0.65
0.81
1.94
21.5
24.2
0.14
0.01
0.81
1.45
24.4
34.2
0.26
0.65
0.96
1.62
27.3
29.4
0.14
0.01
0.81
1.62
21.4
30.2
0.26
0.65
0.96
1.78
24.6
26.5
0.14
0.26
0.81
2.10
19.5
22.6
0.26
0.65
0.96
2.10
19.9
21.3


### 4.2.2. Data Fitting and Regression Analysis

The least square method can be adopted to calculate data in Table7 on regression. Through F Experiment and analysis, it is shown that water depth X3 does not significantly influence pressure regression parameter X5. Thus, under the experiment condition, water depth could be ignored in regression of peak pressure of the shock wave. Experiment data of regression are b0=3.83 MPa, b1=0.097, b2=0.017, and b4=-0.903. Regression equation is(28)Y^=3.83+0.097X1+0.017X2-0.903X4.Therefore, under the condition at the water bottom of soft mud with rocks, in the explosion of single charge in quasi-shallow water, calculating formula of peak pressure of the shock wave at measuring points within the flow field is(29)pm=48.4H¯0.097h¯0.017L¯-0.903.In the equation, fitting range is 0.0529≤H¯≤0.258; 0.0068≤h¯≤0.244; 0.814≤d¯≤0.964; 1.45≤L¯≤2.10.When (29) is employed in calculating explosion in the water bottom of mud with rocks or in quasi-shallow water, actually measured peak pressure and peak pressure of the shock wave at some measuring points are shown in Table 8.Table 8
Actually measured peak pressure and calculated peak pressure of shock wave in water bottom of mud with rocks.
H
¯
h
¯
L
¯
A-pm
C-pm
R
H
¯
h
¯
L
¯
A-pm
C-pm
R
0.05
0.0068
1.45
26.61
23.93
10.07
0.14
0.0068
2.10
17.08
18.24
6.79
0.05
0.0068
1.62
21.12
21.55
2.00
0.14
0.26
1.62
26.08
25.28
3.07
0.05
0.0068
1.78
19.04
19.61
2.99
0.14
0.26
1.62
20.85
21.10
1.20
0.05
0.0068
1.94
17.98
17.99
0.06
0.14
0.26
2.10
19.48
19.49
0.05
0.05
0.26
1.45
24.90
25.57
2.69
0.14
0.65
1.62
25.04
25.71
2.67
0.05
0.26
1.62
22.69
23.04
1.54
0.14
0.65
1.78
22.30
23.39
4.89
0.05
0.26
1.78
21.78
20.96
3.76
0.14
0.65
1.94
22.27
21.46
3.64
0.05
0.26
1.94
20.64
19.23
6.83
0.14
0.65
2.10
20.36
19.82
2.65
0.05
0.52
1.94
19.52
19.47
0.26
0.26
0.0068
1.62
25.96
25.19
2.97
0.05
0.52
2.10
17.66
17.99
1.87
0.26
0.0068
1.78
22.82
22.92
0.44
0.05
0.65
1.45
24.18
26.01
7.57
0.26
0.0068
1.94
20.65
21.03
1.84
0.05
0.65
1.62
22.59
23.43
3.72
0.26
0.0068
2.10
19.26
19.43
0.88
0.05
0.65
1.78
22.22
21.32
4.05
0.26
0.26
1.45
30.68
29.89
2.58
0.05
0.65
1.94
20.53
19.56
4.72
0.26
0.26
1.62
26.60
26.93
1.24
0.14
0.0068
1.45
27.13
26.25
3.24
0.26
0.26
2.10
21.31
20.76
2.58
0.14
0.0068
1.62
23.05
23.65
2.60
0.26
0.65
1.62
27.33
27.39
0.22
0.14
0.0068
1.78
22.05
21.52
2.40
0.26
0.65
1.78
25.41
24.92
1.93
0.14
0.0068
1.94
19.70
19.74
0.20
0.26
0.65
1.94
23.02
22.86
0.70
A
-
p
m: actually measured pm; C-pm: calculated pm; R: relative error (%).Comparison on measured value and calculated value from the regression formula at the same measuring points displays that comparative error is between 0.06% and 10.07%. There are 32 measuring points with proportional error smaller than 5% and only one point larger than 10%. Consequently, regression formula can be employed in calculating peak pressure of the shock wave of a single charge explosion in the quasi-shallow water flow field. And accuracy of the calculating results is relatively high in formula scope.


## 4.2.1. Peak Pressure of Shock Wave in Water Bottom of Soft Mud with Rocks

Under the condition at the water bottom of soft mud with rocks, in model experiment of actual explosion, actual results of experimented peak pressure of some measuring points within the explosion flow field are listed in Table7.Table 7
Incident wave pressure of some measuring points under the condition of water bottom of soft mud with rocks.
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
H
¯
h
¯
d
¯
L
¯
p
m
p
∞
0.05
0.01
0.81
1.45
26.6
34.4
0.14
0.65
0.84
1.45
25.8
32.2
0.05
0.01
0.81
1.62
21.1
30.3
0.14
0.65
0.81
1.62
25.0
28.7
0.05
0.01
0.81
1.78
18.8
27.3
0.14
0.65
0.81
1.78
22.3
26.1
0.05
0.01
0.96
2.10
15.2
22.6
0.14
0.65
0.81
1.94
22.7
23.8
0.05
0.26
0.84
1.45
24.8
34.0
0.14
0.65
0.81
2.10
20.3
21.9
0.05
0.26
0.96
1.62
22.7
30.0
0.26
0.01
0.81
1.45
28.9
33.8
0.05
0.26
0.96
1.78
21.7
27.0
0.26
0.01
0.81
1.94
20.3
24.5
0.05
0.26
0.96
1.94
20.6
24.6
0.26
0.01
0.96
1.62
28.9
29.9
0.05
0.26
0.81
2.10
20.3
22.5
0.26
0.26
0.84
1.45
29.7
34.4
0.05
0.65
0.81
1.94
23.8
23.9
0.26
0.26
0.96
1.62
26.6
30.3
0.05
0.65
0.96
1.78
23.1
25.7
0.26
0.65
0.81
1.94
21.5
24.2
0.14
0.01
0.81
1.45
24.4
34.2
0.26
0.65
0.96
1.62
27.3
29.4
0.14
0.01
0.81
1.62
21.4
30.2
0.26
0.65
0.96
1.78
24.6
26.5
0.14
0.26
0.81
2.10
19.5
22.6
0.26
0.65
0.96
2.10
19.9
21.3


## 4.2.2. Data Fitting and Regression Analysis

The least square method can be adopted to calculate data in Table7 on regression. Through F Experiment and analysis, it is shown that water depth X3 does not significantly influence pressure regression parameter X5. Thus, under the experiment condition, water depth could be ignored in regression of peak pressure of the shock wave. Experiment data of regression are b0=3.83 MPa, b1=0.097, b2=0.017, and b4=-0.903. Regression equation is(28)Y^=3.83+0.097X1+0.017X2-0.903X4.Therefore, under the condition at the water bottom of soft mud with rocks, in the explosion of single charge in quasi-shallow water, calculating formula of peak pressure of the shock wave at measuring points within the flow field is(29)pm=48.4H¯0.097h¯0.017L¯-0.903.In the equation, fitting range is 0.0529≤H¯≤0.258; 0.0068≤h¯≤0.244; 0.814≤d¯≤0.964; 1.45≤L¯≤2.10.When (29) is employed in calculating explosion in the water bottom of mud with rocks or in quasi-shallow water, actually measured peak pressure and peak pressure of the shock wave at some measuring points are shown in Table 8.Table 8
Actually measured peak pressure and calculated peak pressure of shock wave in water bottom of mud with rocks.
H
¯
h
¯
L
¯
A-pm
C-pm
R
H
¯
h
¯
L
¯
A-pm
C-pm
R
0.05
0.0068
1.45
26.61
23.93
10.07
0.14
0.0068
2.10
17.08
18.24
6.79
0.05
0.0068
1.62
21.12
21.55
2.00
0.14
0.26
1.62
26.08
25.28
3.07
0.05
0.0068
1.78
19.04
19.61
2.99
0.14
0.26
1.62
20.85
21.10
1.20
0.05
0.0068
1.94
17.98
17.99
0.06
0.14
0.26
2.10
19.48
19.49
0.05
0.05
0.26
1.45
24.90
25.57
2.69
0.14
0.65
1.62
25.04
25.71
2.67
0.05
0.26
1.62
22.69
23.04
1.54
0.14
0.65
1.78
22.30
23.39
4.89
0.05
0.26
1.78
21.78
20.96
3.76
0.14
0.65
1.94
22.27
21.46
3.64
0.05
0.26
1.94
20.64
19.23
6.83
0.14
0.65
2.10
20.36
19.82
2.65
0.05
0.52
1.94
19.52
19.47
0.26
0.26
0.0068
1.62
25.96
25.19
2.97
0.05
0.52
2.10
17.66
17.99
1.87
0.26
0.0068
1.78
22.82
22.92
0.44
0.05
0.65
1.45
24.18
26.01
7.57
0.26
0.0068
1.94
20.65
21.03
1.84
0.05
0.65
1.62
22.59
23.43
3.72
0.26
0.0068
2.10
19.26
19.43
0.88
0.05
0.65
1.78
22.22
21.32
4.05
0.26
0.26
1.45
30.68
29.89
2.58
0.05
0.65
1.94
20.53
19.56
4.72
0.26
0.26
1.62
26.60
26.93
1.24
0.14
0.0068
1.45
27.13
26.25
3.24
0.26
0.26
2.10
21.31
20.76
2.58
0.14
0.0068
1.62
23.05
23.65
2.60
0.26
0.65
1.62
27.33
27.39
0.22
0.14
0.0068
1.78
22.05
21.52
2.40
0.26
0.65
1.78
25.41
24.92
1.93
0.14
0.0068
1.94
19.70
19.74
0.20
0.26
0.65
1.94
23.02
22.86
0.70
A
-
p
m: actually measured pm; C-pm: calculated pm; R: relative error (%).Comparison on measured value and calculated value from the regression formula at the same measuring points displays that comparative error is between 0.06% and 10.07%. There are 32 measuring points with proportional error smaller than 5% and only one point larger than 10%. Consequently, regression formula can be employed in calculating peak pressure of the shock wave of a single charge explosion in the quasi-shallow water flow field. And accuracy of the calculating results is relatively high in formula scope.


## 4.3. Analysis on Influence of Water Bottom Medium on Peak Pressure of Shock Wave

Figure5 demonstrates change curve where theoretical and experimental values change with proportional distance. With the water bottom of mud with rocks, peak pressure of measuring points distributes near the theoretical values. Comparison of two pictures shows that peak pressure at measuring points on the water bottom of mud with rocks is higher than that of the water bottom of soft mud, which is caused by reflective wave to the water bottom and surface cut off pressure of shock wave and different interactions. To study the influence of different media of the water bottom on peak pressure, experiments are carried out in two different media of water bottom and peak pressure at the same measuring points which are shown in Table 9.Table 9
Comparison of peak pressure in two different media of water bottom (unit: MPa).
H
¯
h
¯
L
¯
r
¯
p
m
Water bottom of soft mudp1m
Water bottom of soft mud with rocksp2m
0.053
0.0068
1.4522
1.4529
17.97
26.56
1.9408
1.9410
12.30
17.97
0.2579
1.4522
1.4666
18.51
24.80
1.6151
1.6280
16.41
22.66
1.7779
1.7897
16.41
21.70
2.1037
2.1137
14.45
20.31
0.5157
1.9408
1.9952
14.45
23.83
2.1037
2.1540
12.50
20.31
0.6515
1.4522
1.5707
20.70
23.83
1.6151
1.7224
18.36
22.36
0.136
0.0068
1.4522
1.4579
20.70
27.16
1.9408
1.9451
12.89
19.66
2.1037
2.1077
14.06
17.10
0.6515
1.4522
1.5410
20.69
25.78
1.9408
2.0081
17.97
22.27
2.1037
2.1659
16.41
20.31
0.258
0.0068
1.4522
1.4738
25.78
29.30
1.6151
1.6345
17.97
28.91
2.1037
2.1186
15.23
21.88
0.2579
1.4522
1.4522
20.31
30.66
1.6151
1.6151
22.31
26.56
1.6151
1.6151
22.79
27.34
1.9408
1.9408
19.53
23.05
2.1037
2.1037
15.63
20.51Figure 5
Change curve of peak pressure of shock wave’s change with proportional distance in different media of water bottom.
(a)
Water bottom of soft mud
(b)
Water bottom of mud with rocksStatistical results of Table9 display that peak pressure at the water bottom of soft mud with rocks is greater than that of soft mud, indicating that the harder the medium at the water bottom is, the greater the explosion load of charging in quasi-shallow water is, and the better the destruction of the explosion can be improved. Based on data in Table 9 and numerous experimental data, relationship betweenp1m and p2m can be obtained by linear regression:(30)p2m=1.33p1m.Under the same experiment conditions in quasi-shallow water explosion, peak pressure of the shock wave in the water bottom of mud with rocks is 1.33 times that of the water bottom of soft mud,(31)r¯=H¯-h¯2+L¯2.Fitting curve of experimentally measured ln⁡pm of different water bottom and ln⁡r¯ can be shown in Figure 6 and responding functional relations arewith water bottom of soft mud(32)ln⁡pm1=3.42-0.99ln⁡r¯,with water bottom of soft mud with rocks(33)ln⁡pm2=3.63-0.85ln⁡r¯.And exponential forms are (34)pm1=30.55r¯-0.99MPa,(35)pm2=37.84r¯-0.85MPa.Analysis on experimentally measured values and results of regression reveals the following:(1)
In (27) and (29) of two different water bottoms, peak pressure of shock wave increases as the increase of explosion height and measuring height, and the influence of explosion height is greater than that of measuring height; besides, peak pressure of measuring points decreases with the increase of horizontal distance, and the influence is greater than explosion height and measuring height; meanwhile, peak pressure at measuring points in water is smaller than that at the same position in the infinite water explosion.(2)
Comparison of (33) and (34) indicates that, in the explosion of quasi-shallow water at the water bottom of soft mud, peak pressure of the shock wave attenuates faster than that in the water bottom of mud with rocks, and amplitude of peak pressure at the water bottom of soft mud is smaller than that of the water bottom of soft mud with rocks.Figure 6
Fitting curve of peak pressure of shock wave with two media of water bottom.


## 5. Discussions

In quasi-shallow water explosion, part to the energy is consumed on crushing or throwing of a medium in the water bottom. Some energy transforms into vibration energy and transmits in the form of seismic wave propagation at the water bottom, and part of energy transmits in the form of the underwater shock wave. And distribution of the above three categories is difficult to study both theoretically and experimentally. Disturbance generated in water spreads between free surface and bottom surface to the water, reflects at the interface, interacts, and causes complex wave system. Therefore, the combined effect of free surface and bottom surface should be taken into consideration when studying explosion in quasi-shallow water [17]. Underwater shock wave generated by the charge explosion in water spreads outward spherically and directly to measuring points, as is shown in Figure 7(1). When the shock wave spreads to the free water surface, it reflects and, similar to the shock wave of the mirrored explosion in air, reflected wave continues to spread in water, as is shown in Figure 7(2). As influence of free surface, shock wave of the underwater explosion acts as the water bottom to make it transform quickly and compacted, generating the reflected wave to spread as Figure 7(3). With different relative position of charging and measuring points, pressure at the points may be unlikely. Peak pressure of measuring points is the result of joint action of Figures 7(1), (2), and (3). However, different relative position may cause unlikely influences of the three.Figure 7
Diagram of shock wave pressure in water at measuring points.Numerous documents reveal that peak pressure of the underwater shock wave exponentially attenuates with an increase of its distance with charging distance, thus, in influencing factors of peak pressure of underwater shock wave, distance between measuring points and charging has the greatest influence. Moreover, explosion height of charging also directly impacts peak pressure of the shock wave. Measuring height of measuring points and explosion height of charging as well as distance between measuring points and the charging has certain relationships, and measuring height can be determined by distance between measuring pints and charging, horizontal distance, and explosion height, so measuring heights also significantly affect peak pressure of the shock wave at measuring points.By analyzing the experimental data, it is obtained that water depth does not significantly affect peak pressure of shock wave, but it cannot be concluded that water depth has relatively smaller impact on underwater shock wave. Water depth directly impacts distance between charging water surface and water bottom. According to boundary effect of free water surface and water bottom on shock wave, water depth and medium at the water bottom affect pressure time history curve of the shock wave at measuring points. At different measuring points, influence of water surface and underwater reflected shock wave also varies. There is also the document showing similar conclusions [18], proposing that water depth has relatively less impact on peak pressure. And within a certain proportional distance, when measuring height is certain, increase of explosion height will increase pressure of measuring points with the same proportional distance.According to the literature [13], if charging explosion energy cannot be absorbed by the water bottom, charging explosion in the water bottom is equivalent to that of a double cartridge of infinite water. Assuming that energy is concentrated in water at the bottom of the water, as an approximation, peak pressure of the double cartridges at any distance increases to 21/3=1.26 times, duration and impulse will also increase, and energy flux density doubles. However, the water bottom cannot be rigid interface and will inevitably absorb part of explosion energy; different media of the water bottom may consume different explosion energy. Moreover, water bottom and reflected wave of water surface have impact; thus, intensity of the shock wave in water may be less than approximation of the double cartridges. Underwater reflected wave has the complex effect on transmission of the underwater shock wave, directly influencing peak pressure of water impact wave. And it is relevant to not only the angle of incidence and angle of reflection, but also properties of underwater media. According to the experimental studies during recent years, when the charging explodes near the water bottom of soft mud, what is, firstly, reflected is not the shock wave (or compression wave); it is reflected sparse wave. With compaction of a medium in the water bottom, reflected compression wave will be reflected later [19]. Effect of the water bottom on underwater shock wave is contrary to that of free surface. When the reflection is generated in underwater sediment layer, rock layer, or surface of other huge obstacles, interaction of reflected wave and waves spreading in water will affect peak pressure of the shock wave. Experiment data of the explosion in different media of the water bottom indicate that peak pressure of the shock wave in the water bottom of mud with rocks is about 1.33 times that of the water bottom with soft mud. Given that water bottom is the absolutely rigid boundary, it may be obtained that underwater shock wave reflects the compression wave, with the same intensity and incident. But it is not in real conditions; as media of water bottom and reflection coefficients are different, peak pressure at measuring points will also vary. Hence, influence of a medium at the water bottom on peak pressure of the underwater shock wave cannot be ignored. Peak pressure of the shock wave in the water bottom of soft mud decreases more quickly than that of the water bottom of soft mud with rocks, demonstrating that water bottom of soft mud absorbs more explosion energy. The higher the rigidity of a medium of the water bottom is, the more the peak pressure of the underwater shock wave will increase.


## 6. Conclusions

By establishing regression model and analysis method, factors impacting peak pressure of the shock wave in a quasi-shallow water explosion are discussed in detail. Regression analysis and corresponding experiment have been carried out on the model of experiment data, and, on the model, following conclusions can be drawn:(1)
Explosion height, measuring height, and water depth as well as horizontal distance all influence peak pressure of the shock wave at measuring points. Horizontal distance from charging center to measuring points has the most significant impact on peak pressure of incident wave front, followed by explosion height, and then comes measuring height, and water depth’s influence is not significant within the range of geometric parameters of these experiments.(2)
Data calculated by the formula of calculating models are more accurate and calculating results have relatively high accuracy within the formula scope.(3)
Explosion experiment of two different media of the water bottom reveals that the media significantly impact peak values of the shock wave. And, under the same experiment conditions, peak pressure through the medium of soft mud with rocks is 1.33 times that of the soft mud water bottom.(4)
Numerous factors influence peak pressure of the underwater shock wave, and plenty of experiment analysis is required for thoroughly exploring transmission law of the underwater shock wave. Due to financial constraints and experimental conditions, certain limitations exist in resulting experiment data in the present paper.



---

*Source: 702178-2015-12-28.xml*