Abstract

Explosion overpressure and its duration are important parameters describing blast loading. However, their values obtained using the currently available methods differ significantly. To address this problem, we simulated the overpressure reflection coefficient (i.e., the ratio of reflection overpressure to incident overpressure) and fitted the analytical expressions to the test data. We then collected a large amount of theoretical and experimental blast loading data from literature and studied the probability distribution of blast loading for different scaled distances. The results show that the probability density function of explosion overpressure obeys an exponential distribution when the scaled distance is between 0.1  and 0.5 . On the other hand, it obeys a normal distribution when the scaled distance is greater than 0.5 . The coefficient of variation reaches the maximum value of 1 when the scaled distance is between 0.1  and 0.5 , while it is between 0.13 and 0.2 when the scaled distance is from 1.5  to 6 . The coefficient of variation of reflected overpressure is slightly bigger than that of incident overpressure. Nevertheless, the coefficient of variation of proportional duration (i.e., the ratio of actual duration to cubic root of charge mass) is generally smaller, between 0.408 and 1.017. It was also found that there is no clear probability distribution for the blast loading duration. Finally, formulas for the explosion overpressure and proportional duration were determined together with their 95% confidence level bounds.

1. Introduction

With the continuous occurrence of terrorist attacks and accidental explosions, the research on blast resistance of building structures has become an important topic in civil engineering. Especially since the 9/11 incident, a significant amount of research has been carried out worldwide on blast resistance of building structures, including blast loading, structural dynamic response, failure modes, simplified calculation methods, blast resistance analysis methods, and blast resistance reinforcement measures. This led to formulating important guidelines for the blast resistance of engineering structures.

The accuracy of blast loading estimation directly affects the results of blast resistant design and is, therefore, the core consideration. A significant amount of research on the blast loading itself has been conducted, and relevant specifications have been formulated to provide efficient guidance for blast resistant design of structures. These guidelines include Fundamentals of Protective Design for Conventional Weapons TM5-855-1 [1], Structures to Resist the Effects of Accidental Explosion TM5-1300 [2], Design of Buildings to Resist Progressive Collapse UFC3-340-02 [3], Design and Assessment of Buildings Subjected to Blast Loads CSA/S850-12 [4], China National Standards Blasting Safety Regulations GB6722-2014 [5], and Code for Design of Civil Air Defense Basement (GB50038-2005) [6]. These standards and specifications summarize theoretical research results and experimental data of blast loading in the form of formulas and graphs, which are widely used in blast resistant design and analysis.

Although calculation methods for blast loading are given in these specifications, the load values vary considerably, even by an order of magnitude, which makes selecting a design load calculation method challenging. Thus, probabilistic models of blast loading calibrated using experimental data were proposed to account for the uncertainty. Low et al. [7] performed a statistical analysis of eight commonly used blast loading determination methods (including theoretical formulas, empirical methods, and experimental data) and found that when the scaled distance was in the range of 0.24–40 , the average coefficients of variation of overpressure and its duration were 0.3227 and 0.130, respectively. Bogosian et al. [8] found that the coefficient of variation of overpressure was 0.24 when the scaled distance was in the range of 1.2–23.8 . Li et al. [911] established a relationship between the average value, standard value, and coefficient of variation of the reflected overpressure and the scaled distance from a statistical analysis of 77 samples of blast loading data and applied them to the reliability analysis of steel beams. The obtained formulas for the reflection overpressure and its duration used a 90% confidence level and can be used for the calculation of blast loading values in building structural design. Thus, these results provided a good reference for the determination of loads in blast resistant design of building structures.

In the previous studies, the mean or variance of blast loading was calculated directly from discrete random data, or the blast loading parameter variables were assumed to obey a normal distribution. The experimental values corresponding to different scaled distances were interpolated from small experimental samples. The reflection overpressure was also directly obtained from the incident overpressure using a theoretical formula without considering the actual changes in the reflection coefficient (the ratio of reflection overpressure to incident overpressure) [12]. Because of these shortcomings, in this paper, the relevant blast resistance specifications, theoretical formulas, and experimental datasets on the explosion overpressure and its duration are reviewed. The reflection overpressure is determined using the reflection coefficient obtained experimentally. The measured blast impulse is converted into an equivalent triangle impulse. From the overpressure and proportional duration (i.e., the ratio of actual duration to cubic root of charge mass) data and histograms for different scaled distances, probability distribution for the explosion overpressure and proportional duration for different scaled distances are determined and validated. The maximum likelihood estimators of the overpressure parameters are obtained and their variations are analyzed. Finally, formulas and confidence intervals for the incident overpressure, reflected overpressure, and proportional duration are given.

2. Blast Loading Characteristics

Explosions produce extremely high-pressure and high-temperature gases. As their volume expands rapidly, air shock waves are formed in front of the compressed gas. The overpressure is a part of the shock wave pressure exceeding the ambient air pressure. The overpressure value decays rapidly with time. When the pressure drops below the ambient air pressure, a negative pressure zone is formed. The overpressure attenuation curve of a typical shock wave is shown in Figure 1. In structural blast resistant design, this negative pressure is often ignored; thus, the main relevant physical parameters are the positive overpressure, its duration, and impulse of blast loading. The common exponential attenuation function is as follows:

In the formula, is the peak value of positive overpressure, is the duration of positive overpressure, and is the attenuation coefficient.

In engineering applications, to simplify the calculations, the overpressure attenuation curve is often simplified as an equivalent descending triangular model [1, 2, 13, 14], as shown in Figure 2. The equivalent triangular loads have the same overpressure peak value, but different simplifications of their duration are adopted according to the expected time of the maximum response of the structure considered as follows: when the maximum response of the structure or component occurs at an early stage of the overpressure application, it is assumed that the slope of the triangular loads is tangent to the actual overpressure curve, and the equivalent action duration is denoted as . This method is often used for long explosion relative to the time of the maximum structural response. When the maximum response of the structure or component occurs after the overpressure drops to zero, it is necessary to ensure that the triangular impulse is consistent with the actual overpressure impulse. The equivalent action time is denoted as , and this method is often used for short explosion relative to the time of the maximum structural response. For the intermediate cases, it is assumed that the triangular impulse passes through of the actual overpressure time history plot to calculate the duration . For the commonly used advanced explosives, the second simplified method for obtaining the equivalent impulse is often used because the overpressure action time is often very short compared to the relative time of structural response.

3. Explosion Overpressure Study

3.1. Theoretical Overpressure Data

A large number of researchers, including Brode [1517], Henrych [18, 19], Mills [12], Kinney [20], and Newmark [21], studied the shock wave overpressure and derived formulas by introducing empirical coefficients under ideal air explosion conditions. However, there are still significant differences among the different formulas, especially close to the explosion center, as shown in Figure 3. Therefore, this study examined as many as 17 commonly used explosion overpressure formulas, as shown in (2)–(18). Of them, (2)–(17) can be used to directly calculate the incident overpressure, and (18) can be used to directly calculate the reflected overpressure.

Kinney et al. [20] gave the formula for calculating peak incident overpressure (Mpa) under free air explosion as follows:

Henrych et al. [18, 19] gave the expression of incident overpressure peak (bar) for ideal air explosion based on the experimental method as follows:

Baker et al. [22] gave the expression of incident overpressure peak (MPa) for ideal air explosion as follows:

Wu and How [23] gave the expression of incident overpressure peak (MPa) during the ground explosion as follows:

Naumyenko et al. [19, 24, 25] gave the expression of incident overpressure peak (bar) under ideal air explosion according to the model similarity theory and determined the coefficients by experimental method. Among them, the first equation is obtained by Naumyenko and Petrovsky, and the second equation is derived by Sadovskii. The specific calculation expression is as follows:

Brode [1517] gave the expression of incident overpressure peak (bar) for ideal air explosion by using similar model theory as follows:

Mills [12] combined the similarity theory and numerical simulation and gave the expression of incident overpressure peak (MPa) for ideal air explosion as follows:

Ye [26] gave the expression of incident overpressure peak (MPa) for ideal air explosion as follows:

The expression of incident overpressure peak (bar) for explosive ground explosion provided in Chinese National Standard Blasting Safety Regulations [5] is as follows:

The expression of incident overpressure peak (MPa) for explosive ground explosion provided in the National Standard Code for Design of Civil Air Defense Basement [6] is as follows:

The expression of incident overpressure peak (MPa) for ground explosion proposed by Newmark and Hansen [21] is as follows:

Siddiqui [27] Ahmad [28], and Iqbal [29] gave the expressions of incident overpressure peak (MPa) for ground explosion based on experimental data as follows:

The expression of incident overpressure peak (KPa) of air explosion given in the US military protection technology manual TM5-855-1 [1] is as follows:

Held et al. [30] gave the expression of incident overpressure peak (MPa) in air explosion based on theoretical analysis as follows:

The expression of reflection overpressure peak (KPa) for ground explosions given in the Canadian specification CSA/S 850-12 [4] is as follows:

It can be seen that there is a very limited number of theoretical or empirical formulas for the reflection overpressure; thus, more reflection overpressure data need to be obtained using the reflection coefficient and the incident overpressure. Mills proposed a theoretical expression for the reflection coefficient, as shown in (19). In the US Military Protection Technology Manual TM5-855-1 [1], the variation of reflection coefficient with the incident angle is considered (Figure 4). A comparison of the reflection coefficients obtained by (19) and Figure 4 is shown in Figure 5 for the vertical incident angle (0°). Figure 5 shows that when the incident overpressure is small, the theoretical and experimental reflection coefficients are equal. However, when the incident overpressure is large, the maximum experimental reflection coefficient is 13.2, which is much larger than the theoretical value of 8. Therefore, when the reflection overpressure is calculated using the incident overpressure and the reflection coefficient, the reflection coefficient should be determined according to the experimental value. To obtain the experimental reflection coefficient conveniently for the case of vertical incidence, a formula for the reflection coefficient, μ, for different overpressures obtained by fitting experimental data is shown in (20), and the regression coefficient is 0.9871, as shown in the imaginary line in Figure 5. Thus, the relationship formula equation (21) between the reflection overpressure and the incident overpressure is obtained.

3.2. Experimental Overpressure Data

Several researchers have carried out field blast resistant tests. The test data gathered will help to analyze the explosion overpressure more accurately. In this study, 125 reliable explosion overpressure experimental datasets were collected from literature. The data cover the charge mass from 30 g to 50 kg, and the scaled distance from 0.295 to 55 , as shown in Table 1. It can be seen from the table that the blast loading is sensitive to the environmental factors, and the overpressure values obtained by multiple tests under the same conditions are different. The explosion distances in the tests were relatively small, mostly within 5 m, which was mainly because of the small equivalent explosive masses and desire to reduce the influence of environmental factors. At the same time, considering the stability and accessibility of explosives, some tests used ANFO explosives (equivalent coefficient of 0.82), emulsion explosives (0.71), Composite C-4 (1.37), and Gelamon VF80 (0.8), rather than TNT explosives, which also increased the variability of the test data. In some tests, the charge was placed at a certain height from the ground, and the reflection of the shock wave off the ground also affected the overpressure. In some tests, the charge shape was nonspherical or irregular, especially when the equivalent explosive mass was small, which had some influence on the overpressure. Therefore, for the experimental explosion data listed in Table 1, the above factors cannot be completely avoided, and these differences are also important factors leading to explosion overpressure uncertainty.

In the collected experimental datasets, there are both vertical and nonvertical incidence cases. However, for the sake of consistency, only the vertical incidence cases were retained for analysis. Furthermore, some datasets had only the incident overpressure samples, some had only the reflected overpressure samples, and only 11 datasets had both the incident and the reflected overpressure. For the same batch of tests, 22 datasets exhibited significant differences and were excluded from the subsequent analysis. Eventually, a total of 103 datasets were used.

The existing formulas for the incident explosion overpressure (2)–(18) have the general forms as shown in (22) and (23). In each formula, there are three undetermined coefficients. Equation (22) is obtained from the mechanical explosion theory and its correlation coefficient is determined by experiments. This method has been more widely applied in the existing formula for the incident overpressure. Equation (23) is obtained by fitting the experimental data, and it also has certain applications. Therefore, the reflection coefficient formula (20) and the incident overpressure formula (22) can be used together to gain the overpressure data under certain scaled distance, which provides statistical samples for analyzing the probability distribution of the explosion overpressure.

3.3. Overpressure Probability Distribution

Using the explosion overpressure data, histograms corresponding to different blast scaled distances were drawn. In addition to the distribution characteristics observed in the histograms, the K-S and Lilliefors goodness-of-fit tests of overpressure probability distribution for different scaled distances were carried out. The results showed that when the scaled distance was less than 0.5 , the incident and reflected overpressure both followed exponential distributions. However, when the scaled distance was greater than 0.5 , the overpressure distributions were normal. Using these probability distributions, the maximum likelihood estimates of the expectation and standard deviation of the overpressure for different scaled distances were calculated. Furthermore, the 95% confidence interval of the overpressure was obtained, as shown in Tables 2 and 3. Then, the probability density functions (exponential and normal) of overpressure for different scaled distances were drawn as shown in Figures 6 and 7. It can be seen from Tables 2 and 3 that when the scaled distance is small, the standard deviation is the largest, and thus the overpressure distribution is wider.

To compare the dispersion degree of the incident and reflected overpressure for different scaled distances, the coefficients of variation of overpressure distribution were obtained from the expectations and standard deviations in Tables 2 and 3, as shown in Figure 8. It can be seen from Figure 8 that when the scaled distance was less than 0.5 , the coefficient of variation of overpressure reached the maximum value of 1, which is constant of the exponential distribution. When the scaled distance was about 1.5–6 , the coefficient of variation was small due to more experimental data. In the range 0.13–0.2  the coefficient of variation of the reflection overpressure was slightly larger than that of the incident overpressure.

3.4. Determination of Overpressure Formula

Using the maximum likelihood estimations of the overpressure expectations shown in Tables 2 and 3, the relationships between the upper and lower 95% confidence limits of the incident and reflected overpressure and the scaled distance can be obtained, as shown in Figures 9 and 10, respectively. Since the overpressure values for different scaled distances varied considerably, the scaled distances are plotted in three sections. It can be seen from the figure that the smaller the scaled distance, the more dispersed the overpressure distribution. A polynomial was fitted to the natural logarithm of the scaled distance and overpressure data using the least square method. As shown in Figures 11 and 12, the regression coefficients of the fitted curves of the incident and the reflected overpressure are greater than 0.99.

Finally, the formula for calculating incidence overpressure is as follows:

The upper and lower limits of 95% confidence interval for incidence overpressure are calculated as follows:

Similarly, the formula for calculating reflection overpressure is as follows:

The upper and lower limits of 95% confidence interval for reflection overpressure are calculated as follows:

In practical engineering applications, the recommended values of incident and reflected overpressure are calculated using (24) and (27), respectively. When there is sufficient experience supported by data, (25), (26), (28), and (29) can be used as the 95% confidence upper and lower limits of the overpressures under such conditions. It should be noted that when the scaled distance is very small (<0.1 ), the explosion overpressures should be determined from special studies on the proximity or contact explosion.

4. Explosion Duration Study

4.1. Theoretical Duration Data

Duration is another important parameter of blast loading, and different formulas were proposed using theoretical derivations and experimental data, as shown in (30)–(39). The relationships between the blast loading and time use either the actual exponential form (Figure 1) or simplified triangular form (Figure 2). The Code for Design of Civil Air Defense Basement [6], Held [30], and Canadian specification CSA/S850-12 [4] all use triangular form, while some other sources use the exponential form. In the statistics of theoretical data, the most used Conwep program calculation results are included. Then, taking 25 kg and 50 kg explosive equivalents as examples, the relationships between the duration and scaled distance for different explosive equivalents are shown in Figure 13. It can be seen that the overpressure durations for the different conditions are also quite different; thus, it is necessary to carry out statistical analysis on the duration of blast loading. It can be seen from the formulas that the duration, t, is mainly related to the explosive equivalent and scaled distance. Using the similarity rate, the functional relationship between parameter and the scaled distance can be obtained, and the statistical characteristics of the blast loading duration will be subsequently studied using this parameter.

For high-explosive TNT, Henrych [18, 19] obtained the calculation formula of overpressure duration for air explosion according to the experimental research work as follows:

Sadovskii [25] gave the calculation formula of explosion overpressure duration under the condition of air explosion as follows:

For coefficient B, Naumyenko [24] recommended 1.0, Pokrovskii [48] considered 1.3, and Baum [49] gave 1.5.

According to Yang et al. [50], through the numerical simulation results of parameter analysis and fitting, the relationship of positive overpressure duration is as follows:

In the process of data calculation, the formula has serious deviation when the scaled distance is large, so only the scaled distance less than 3  is considered.

In the China Code for the Design of Civil Air Defense Basement [6], the suggested overpressure duration formula for the ground explosion is as follows:where is reflected overpressure, the unit is MPa, the unit of is ms, and the unit of is kg.

The calculation formula (triangular load) of incident overpressure and incident impulse of spherical explosive for air explosion given by Held [30] is as follows:

The duration after conversion is as follows:

In the equation, the unit of is ms, the unit of is kg, the unit of is MPa, and the unit of is Pa-s.

Wu and Hao [23] suggested using Henrych' s research results and giving the formula for calculating the overpressure duration of air explosions as follows:

Sui [51] obtained the calculation formula of duration according to explosion similarity law as follows:

Canadian specification CSA/S 850-12 [4] gives the calculation formula of reflection overpressure and reflection impulse when the ground explosion occurs. Thus, the duration of overpressure can be inversely calculated. Among them, the calculation formula of reflection impulse is as follows:

The unit of is KPa-ms.

4.2. Experimental Duration Data

Similar to the experimental overpressure data, 125 datasets were collected. The experimental data related to the blast loading duration are shown in Table 4. Only 87 datasets were useable. In the experiments, the measured load duration is the total duration of exponentially declining overpressure. To facilitate statistical comparisons, all durations were converted into triangular durations form. The conversion principle was that the positive impulse and overpressure peak value of the exponential blast loading and triangular blast loading were consistent. It can be seen from Table 4 that, under the same conditions, the duration of blast loading has high discreteness.

By analyzing the theoretical duration formulas (30)–(39), it was found that, unlike the formulas of overpressure, the ratio of duration to the cube root of charge mass (i.e., the proportional time) was only related to the single parameter of scaled distance. Therefore, the data from theoretical calculations and experimentally measured data were changed into the proportional time data, and the relationship between the proportional time and the scaled distance was studied to determine the probability distribution of duration.

4.3. Duration Probability Distribution and Formula Derivation

Using the explosion overpressure duration data, the histograms of the proportional time () for different scaled distances were drawn, as shown in Figure 14. Combined with the distribution characteristics of the duration histograms, the K-S and Lilliefors goodness-of-fit tests of proportional time probability distributions for different scaled distances were carried out. The results showed that the proportional times for different scaled distances differed considerably and no specific probability distribution was acceptable. Therefore, the statistical data of proportional time for different scaled distances were used directly. The relationship between the proportional time and the scaled distance can be obtained by the least square method, as shown in Figure 15.

The function formula of proportional time with scaled distance is as follows:

The curve formulas corresponding to the upper and lower limits of 95% confidence interval are as follows:

To compare the dispersion degree of the proportional time for different scaled distances, the coefficient of variation of the proportional time was calculated, as shown in Figure 16. As can be seen from Figure 16, for different scaled distances, the coefficients of variation of the proportional time were overall small, between 0.408 and 1.017. The variability becomes smaller with the increase in scaled distance. At the same time, the formula for the proportional time as a function of the scaled distance was highly reliable.

5. Conclusions

By analyzing a large amount of blast loading experimental and theoretical data, the probability distribution of blast loading was studied. The main conclusions are as follows:(1)Blast loading test data were affected by factors such as explosive type, equivalent mass and shape, and test environment. Both the overpressure and load duration varied. The experiments mostly used small equivalent explosives. The scaled distance was mostly concentrated in the range of 0.4–2.0 .(2)A formula for the reflection coefficient was proposed by curve fitting to the vertical incidence data, and then a formula for the reflection overpressure using the incident overpressure was obtained.(3)When the scaled distance was less than 0.5 , the probability density of explosion overpressure followed an exponential distribution, and the coefficient of variation of overpressure reached the maximum value of 1. When the scaled distance was greater than 0.5 , the probability density of explosion overpressure followed a normal distribution. When the scaled distance was about 1.5–6 , the coefficient of variation was small, between 0.13 and 0.2.(4)For different scaled distances, no specific probability distribution function for the proportional duration could be selected, but the coefficient of variation of the proportional duration was relatively small, between 0.408 and 1.017.(5)Using the data of explosion overpressure and proportional duration for different scaled distances, the formulas for the explosion overpressure and proportional duration and the 95% confidence upper and lower limit formulas were derived.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 51608229). And the authors would like to express their gratitude to EditSprings (https://www.editsprings.com/) for the expert linguistic services provided.