Next Article in Journal
Opportunity to Use New Ways of Working: Do Sectors and Organizational Characteristics Shape Employee Perceptions?
Next Article in Special Issue
A Constitutive Model Study of Chemical Corrosion Sandstone Based on Support Vector Machine and Artificial Bee Colony Algorithm
Previous Article in Journal
The Impact of Wearable Devices on the Construction Safety of Building Workers: A Systematic Review
Previous Article in Special Issue
Study on the Prediction of the Uniaxial Compressive Strength of Rock Based on the SSA-XGBoost Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of an Explosion-Proof Wall on Shock Wave Parameters and Safe Area Prediction

1
Shaanxi Key Laboratory of Safety and Durability of Concrete Structures, Xijing University, Xi’an 710123, China
2
School of Environment and Resource, Southwest University of Science and Technology, Mianyang 621010, China
3
Chengdu Institute of Urban Safety and Emergency Management, Chengdu 610011, China
4
Mechanical Engineering Research Institute, Xi’an 710123, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(14), 11164; https://doi.org/10.3390/su151411164
Submission received: 6 June 2023 / Revised: 2 July 2023 / Accepted: 8 July 2023 / Published: 18 July 2023

Abstract

:
To study the influences of an explosion-proof wall on shock wave parameters, an air explosion protection experiment was performed, the time history of shock wave pressure at different positions before and after the explosion-proof wall was established, and the characteristics of shock wave impulse and dynamic pressure were analyzed. The explosion-proof working conditions of five different diffraction angles were simulated and analyzed using Autodyn software(2019R3). Results indicated the following findings. The explosion-proof wall exerted an evident attenuation effect on the explosion shock wave, but considerable pressure still existed at the top of the explosion-proof wall. Overpressure behind the wall initially increased and then decreased. The larger the diffraction angle, the faster the attenuation speed of the diffraction overpressure of the shock wave in the air behind the wall. The history curve of shock wave pressure exhibited an evident bimodal structure. The shock wave diffraction of the wall made the shock wave bimodal structure behind the wall more prominent. The characteristics of the bimodal structure behind the wall (the interval time of overpressure peak Δt was less than the normal phase time of the diffracted shock wave T+) caused the shock wave impulse to stack rapidly, significantly improving its damage capability. The peak value of dynamic pressure on the oncoming surface was approximately two times the peak value of overpressure, and the inertia of air molecules resulted in a longer positive duration of dynamic pressure than overpressure. The maximum overpressure on the ground behind the explosion-proof wall appeared at approximately two times the height of the explosion-proof wall, while the maximum overpressure in the air behind the explosion-proof wall appeared at approximately one times the height of the explosion-proof wall. The relatively safe areas on the ground and in the air behind the wall were approximately 4–4.5 times and 3.5–4 times the height of the explosion-proof wall, respectively.

1. Introduction

As a modern field protection facility, an explosion-proof wall can effectively isolate the explosion source and protect personnel, considerably reducing the damage and destructive effects of explosive shock waves and shrapnel on important military targets and personnel [1,2,3], thus improving the survivability of military personnel. As a major wartime damage element, the blast shock wave will be reflected near the ground, superimposed with the incident shock wave, and act on the blast-facing surface of the blast wall, causing the blast wall to be severely damaged or even overturned [4]. Simultaneously, the explosion shock wave diffracts through the upper and both sides of the wall after coming in contact with the explosion-proof wall, changing the shock wave pressure flow field behind the wall and even forming a higher overpressure in a specific area behind the wall [5,6], causing damage to the evading personnel. Therefore, studying the influences of an explosion-proof wall on the explosion shock wave flow field and the law of shock wave parameters is highly significant to evaluate the protective performance of an explosion-proof wall [7,8,9,10].
The anti-overturning and stability of explosion-proof walls have achieved certain progress in domestic and foreign research. Zhang [11] and Scherbatiuk K [3,4,6] established a finite element model of the explosion load and the response characteristics of an explosion-proof wall, proposed an analytical model of such a wall based on the rotation of a rigid body, and analyzed the response characteristics of an explosion-proof wall on the basis of the displacement time curve measured in the experiment. Explosion-proof walls of different materials exert varying blocking effects on the explosion shock wave. When the explosion shock wave passes through a concrete wall [12], it will produce an evident circulation phenomenon, which exerts a significant effect on the flow field behind the wall; the buffer effect of the soft medium explosion-proof wall of the water pack is evident [13,14]. The environment of an explosion-proof wall also affects shock wave overpressure: when explosives detonate in different environments, such as air, rigid ground, and sandy ground, they affect the distribution of overpressure in different directions [15]. The different geometric dimensions of an explosion-proof wall affect the diffraction angle and also exert a considerable influence on the propagation of shock waves [16,17,18]. Zhu and Mu [19,20] determined that the maximum value of the diffraction overpressure behind the wall is related to the value of the diffraction angle. However, the relationship between the diffraction angle and the shock wave diffracted flow field has not yet been studied in detail. In summary, the environmental and geometric parameters of an explosion-proof wall exhibit an effect on its anti-explosive performance, and thus, the shock wave overpressure [21,22], specific shock volume, and dynamic pressure of such a wall [23,24] should be comprehensively analyzed when evaluating its anti-explosive performance. Most studies have focused only on one of the three parameters and cannot completely investigate the impact protection performance of an explosion-proof wall. The choice of an appropriate material dynamic model has a significant impact on the results of numerical calculations. Ma [25] proposed an improved Johnson–Holmquist–Beissel model that can effectively predict the projectile penetration behavior of brittle materials. A new modeling strategy with inclusions and an FDEM method is also proposed to study the dynamic response and fracture behavior of geomaterials [26]. In this paper, the compaction nonlinear state equation is used to describe the effect of blast walls in blast loading.
In the current study, two types of sensors were arranged on the ground and in the air behind an explosion-proof wall, considerably covering the hazard area of the shock wave behind the explosion-proof wall for personnel and equipment. The calculation of shock wave overpressure, impulse, and dynamic pressure at each measuring point was analyzed in detail. The diffraction angle range was set to 22°–44°. The corresponding explosion distance can be determined in accordance with its geometric relationship. The pressure variation law of diffracted shock waves from varying angles at different positions of the explosion-proof wall was studied through an experiment and numerical simulation to identify the parameters related to explosion safety protection and to provide a reference for designing the size and structure of an explosion-proof wall, safely avoiding the damage caused by shock waves.

2. Experiment

2.1. Shock Wave Testing System

The shock wave pressure sensors used in this work were 10 wall-type PCB shock wave pressure sensors with a minimum resolution of 0.10 kPa and 10 free-field PCB shock wave pressure sensors with a minimum resolution of 0.07 kPa. The shock wave data acquisition instrument adopted a Chengdu Tytest Blast-PRO shock wave tester with a sampling rate of 4 MHz and a Jiangsu Donghua DH5960 data acquisition instrument with a sampling rate of 20 MHz. The analysis system was primarily composed of the PC terminal and its supporting software. The shock wave test system is shown in Figure 1.

2.2. Experimental Design

Two sets of explosion-proof wall anti-explosion experiments were established. The trinitrotoluene (TNT) equivalent was 1.5 kg, and the explosion height was 1.2 m. The proportional distance between the explosion-proof wall and the explosion source, the distance, and the corresponding relationship of the diffraction angle α are provided in Table 1 and Figure 2. The experimental explosion-proof unit consisted of five individual sand-filled blast walls with a single blast wall geometry of 1 m × 1 m × 2.1 m.
The types of deployed sensors were mostly wall-type shock wave sensors (measurement points 1–10) and free-field shock wave sensors (measurement points 11–20). The height of the free-field pressure sensors from the ground was 1.2 m. The shock wave sensors were arranged at 1, 2, 3, 4, 5, 6, 7, and 8 m from the back of the explosion-proof wall. The literature [19] pointed out that the maximum overpressure behind the wall occurs within 1.5–2.5 times the height of the explosion-proof wall behind the wall. The distance between measuring points 19 and 20 and the charge is consistent with the horizontal distance from measuring points 5 (13) and 9 (17) to the explosion source, respectively. The outer layer of the explosion-proof wall used in the experiment was a reinforced mesh frame, an inner layer lined with a geotextile pocket protective sheet, and an inner layer filled with an energy-absorbing material (i.e., sand). The height of the wall was 2.1 m, the width was 5 m, and three soil samples were collected. The measured physical parameters are listed in Table 2. The experimental explosion-proof unit consisted of five individual sand-filled blast walls with a single blast wall geometry of 1 m × 1 m × 2.1 m.

2.3. Layout of Measuring Points

This explosion experiment was performed in an open field. The sensor should be level when fixing the air free-field sensor. In accordance with the actual situation, the shock wave will be diffracted on the top, left, and right sides of the explosion-proof wall. To ensure the integrity of the shock wave signal test, the plane part of the sensor was directed upward, and the sensor was pointed toward the grain tested and on the same horizontal plane as the center of the explosion source. Meanwhile, to ensure the stability of the free-field sensor in the air, a fixed sandbag was placed at the bottom of the fixed support and secured with a nylon rope. Finally, the on-site layout is depicted in Figure 3.

3. Results and Analysis

3.1. Shock Wave Overpressure

Considering the large number of measuring points in this experiment, the pressure time history curves of some typical measuring points in Experiment 1-1 were selected for explanation below in order to avoid this paper being too lengthy, as shown in Figure 4.
Figure 4 shows the shock wave pressure time history curves of several measuring points in Experiment 1-1, where P1 is the peak pressure of the diffraction wave; P2 is the peak pressure of the reflected wave; and t1 and t2 are the arrival time of the diffracted and reflected waves, respectively. The aforementioned measuring points were located on the ground (9) and in the air (14) of the back explosion surface, and in the air of the front explosion surface (19,20). The diffraction overpressure peak, reflection overpressure peak, and the corresponding time of the two overpressure peaks at each measuring point were provided. The shock wave pressure time history curve of the two measuring points in the air of the blast-facing surface under the action of a 1.5 kg TNT air explosion is similar to the ideal shock wave pressure time history curve. The peak incident overpressure of the shock wave at measuring points 19 and 20 was 31.57 kPa and 15.06 kPa, respectively. In accordance with the empirical formula of Hengrych J [27] (Formula (1)), the shock wave overpressure at measuring point 19 was 34.91 kPa and that at measuring point 20 was 14.56 kPa, which are basically the same as the measured shock wave overpressure, demonstrating the reliability of this test. Evident “double peak” waveforms can be observed in the shock wave pressure time history curve of each measuring point behind the wall, and this time history curve differs from the typical air free-field shock wave pressure time history curve. The reflected wave overpressure appeared when the diffracted wave attenuation process did not enter the negative phase. The shock wave parameters of the second experiment are provided in Table 3 and Table 4.
Δ P m = 1407.17 / R ¯ + 553.97 / R ¯ 2 35.72 / R ¯ 3 + 0.625 / R ¯ 4 , 0.05 R ¯ 0.3 619.38 / R ¯ 32.62 / R ¯ 2 213.24 / R ¯ 3 , 0.3 R ¯ 1 66.2 / R ¯ + 405 / R ¯ 2 + 328.8 / R ¯ 3 , 1 R ¯ 10
where ∆Pm is the peak overpressure, kPa; R ¯ = R / W 3 ; R ¯ is the proportional distance, m·kg−1/3; R is the blast distance, m; and W is the TNT equivalent, kg.
Table 3 and Table 4 show that the explosion-proof wall exerted a significant effect on reducing the shock wave of the explosion caused by 1.5 kg of TNT. The overpressure of the front explosion surface of the two actual explosion experiments was 3845.50 kPa and 3754.10 kPa. A large overpressure remained when the top of the explosion-proof wall was not protected, and the shock wave overpressure was lower on the back side of the explosion-proof wall. The overpressure distribution behind the explosion-proof wall is shown in Figure 5 and Figure 6. The time relationship between the reflected and diffracted waves is illustrated in Figure 7.
Figure 5 suggests that the diffracted and reflected waves can hardly reach the bottom behind the explosion-proof wall when the diffraction angle was 32°. Moreover, the diffracted overpressure was 16.47 kPa and the reflected overpressure was 1.51 kPa at measuring point 3. In accordance with the TNT air explosion conditions and the propagation path of the shock wave bypassing the explosion-proof wall, the diffraction overpressure was primarily caused by the explosion shock wave diffracting from the top of the explosion-proof wall, and the overpressure reflected was mostly formed by the shock wave being diffracted and reflected along the two sides of the explosion-proof wall. The overpressure reflected by the ground measure point behind the wall would be greater than the diffraction overpressure after approximately one times the height of the explosion-proof wall. In the subsequent measurement points, the reflected overpressure would be greater than the diffraction overpressure, and then the maximum value of 19.69 kPa would be reached at a distance of 5 m (approximately 2.4 times the height of the explosion-proof wall). Here, the diffraction overpressure also reached the maximum value of 18.51 kPa, exhibiting a trend of initially increasing and then decreasing. In the air of the back explosion surface, diffraction overpressure gradually decreased with an increase in the distance from the explosion-proof wall. The maximum diffraction overpressure was 22.49 kPa at 1 m behind the wall, and reflection overpressure initially increased and then decreased. The maximum reflection overpressure was 22.75 kPa at 5 m away from the explosion-proof wall. At 3 m behind the wall (approximately 1.5 times the height of the explosion-proof wall), air reflection overpressure began to become greater than diffraction overpressure. As shown in Figure 6, the shock wave overpressure on the ground behind the wall initially increased and then decreased when the diffraction angle was small, with a slight difference between diffraction overpressure and reflection overpressure. In the measuring points after approximately two times the height of the explosion-proof wall, reflection overpressure was always greater than diffraction overpressure, and the maximum diffraction overpressure was 13.97 kPa at 5–6 m away from the explosion-proof wall (approximately 2.4–2.6 times the height of the explosion-proof wall). In the air above the back explosion surface, diffraction overpressure also decreased as the distance from the explosion-proof wall increased. The maximum diffraction overpressure was 20.92 kPa at 1 m behind the wall. Figure 5b and Figure 6b show that as the diffraction angle decreased, the attenuation rate of diffraction overpressure ( n = P 11 P 18 / P 11 ) in the air behind the wall decreased by 6.9%, indicating that the greater the diffraction angle, the faster the attenuation rate of diffraction overpressure in the air behind the wall increased with the distance from the explosion-proof wall. Reflection overpressure initially increased and then decreased. The maximum reflection overpressure was 20.14 kPa at 3 m from the explosion-proof wall (approximately 1.5 times the height of the explosion-proof wall). Simultaneously, reflection overpressure in the air began to become greater than diffraction overpressure, and reflection overpressure was always greater than diffraction overpressure.
Figure 4, Figure 5 and Figure 6 suggest that the TNT explosion shock wave will exhibit two more evident overpressure peaks in the pressure time history curve. Figure 7 shows that the interval time between diffraction overpressure and reflection overpressure in the pressure time history curve gradually decreased with an increase in the distance from the explosion-proof wall. This interval time is also a key factor that affects the destructive effect of the explosion shock wave. If the interval time is less than the positive phase time of the diffraction shock wave, then the two overpressure peaks are in the stage of rapid impulse rise, and the impulse is superimposed rapidly, exerting a strong destructive effect.

3.2. Overpressure Impulse

(1)
Impulse
In addition to the shock wave overpressure describing the power and destructive capability of the ammunition when it explodes, the impulse is also one of the important indicators. Its positive impulse is expressed as
I + = t 1 T + P t P 0 d t
where I+ is the impulse of the positive shock wave, Pa·s; t1 is the arrival time of the diffraction wave, ms; T+ is the normal phase duration of the diffracted wave, ms; P t is the pressure at this location as a function of time, kPa; and P0 is the ambient pressure, kPa. The pressure time history curve obtained in the experiment was based on atmospheric pressure.
The impact of ground reflected waves on the explosion impulse cannot be ignored. If they are located below the three-wave point traces, then the ground reflected waves will obviously affect the pressure peak, the explosion directly generated by the incident wave, and will be superimposed and enhanced. If located above the three-wave point traces, the incident wave and the ground reflected waves will be significantly separated, weakening the pressure peak. The explosion-proof wall was installed in a field explosion; thus, the shock wave reflection situation was complicated and evidently different from the attenuation history in the free-field pressure time history curve of an unobstructed air explosion. The overpressure of the reflected wave at some positions was greater than that of the diffraction. Therefore, the explosive impulse studied in this work was not only the positive impulse phase a of the diffracted wave attenuation history but also included the negative impulse phase b between the diffracted and reflected waves and impulse phase c of the reflected wave attenuation history. If the reflected wave arrives before the pressure of the diffracted wave enters the negative phase, then phase b will not occur. If the overpressure of the reflected wave has a negative value, then only the positive impulse a of the diffracted wave is calculated, as shown in Figure 8.
(2)
Characteristics of the impulse curve at different positions
Figure 9 shows the impulse time history curve at several positions, illustrating the evident difference between the arrival of the reflected wave before and after the pressure of the diffracted wave decays into the negative phase. In Figure 9, a and b illustrates the impulse curves of the two measuring points behind the wall. The reflected wave reached the measuring point before entering the negative phase in the attenuation process of the diffracted wave. The impulse increased rapidly after the arrival of the reflected wave. For example, the impulse of the diffracted wave at measuring point 9 was 15.17 Pa·s. The impulse continued to increase with the arrival of the reflected wave, and the maximum impulse was 43.19 Pa·s. At measuring point 14, the impulse of the diffracted wave was 30.08 Pa·s and the maximum impulse was 57.13 Pa·s. By contrast, for c and d in Figure 9, the reflected wave arrived after the diffracted wave entered the negative phase, indicating that its impulse decreased. In this case, the impulse of the diffracted wave was generally considerably greater than that of the reflected wave. For example, the impulse of the diffracted wave at measuring point 19 was 49.41 Pa·s, and the maximum impulse was 54.65 Pa·s. At measuring point 20, the impulse of the diffracted wave was 35.66 Pa·s, which was also the maximum impulse. The peak overpressure of the reflected wave was small because the reflected wave arrived late, and a continuous impulse rising phase did not occur. Thus, the impulse at the measuring point dropped rapidly.
Table 5 and Table 6 provide the positive phase impulse of the diffracted wave, the cumulative impulse of the reflected wave before reaching the negative phase, and the maximum impulse of the reflected wave before entering the negative phase after the arrival of the reflected wave in Experiments 1-1 and 1-2. In the two experiments, the maximum impulse of most measuring points was equal to the accumulated impulse before the reflected wave reached the negative phase. A considerable difference in value was noted compared with the positive phase impulse of the diffracted wave because the reflected wave of the explosion shock wave arrived before the diffracted wave entered the negative phase. Thus, the impulse history curve continued to rise. At measuring points 2, 18, and 19, the maximum impulse of the reflected wave before entering the negative phase was close to the positive impulse of the diffracted wave, or even less than the diffracted impulse. This condition was attributed to the reflected wave reaching the measuring point after the diffracted wave entered the negative phase. The negative phase impulse offset a large amount of reflected impulse.
Figure 10 shows the impulse comparison chart of the air and ground measurement points in the two experiments at different distances behind the explosion-proof wall and the average value of the impulse of the air and ground measurement points. When the diffraction angle was 32°, the mean impulse values of the air and ground measuring points were basically equal, and the impulse values at the two measuring points at different heights were close at the same horizontal distance from the explosion source. When the diffraction angle was 28°, a certain difference was noted in the mean impulse value between the air and ground measuring points. The recorded data show that the impulse values recorded by the air and ground measuring points gradually approached each other with an increase in distance. The maximum impulse in the air behind the wall of the two experiments was 3–5 m behind the wall (approximately 1.5–2.4 times the height of the wall). These values are consistent with the diffraction and reflection overpressure in b of Figure 5 and Figure 6. At the two diffraction angles, the average air impulse was greater than the ground average impulse, and the change in diffraction angle exerted a minimal effect on the air diffraction impulse.

3.3. Shock Wave Dynamic Pressure

With the rapid development of new high-energy weapons and their battlefield applications, in contrast with shock wave overpressure, diffraction pressure and reflected pressure last only for a short time (typically much less than 1 s). However, dynamic pressure may last longer (up to 2–3 s), and the damaging effect of dynamic pressure becomes more prominent. At present, dynamic pressure measurement experiences many technical difficulties. The literature [23] asserted that air can be regarded as a hot complete gas when shock wave overpressure was less than 4 MPa and wave front temperature was less than 2000 K. The current study used the literature [24] as a reference to present a method for calculating typical air impact fluctuating pressure under an ideal explosion condition (i.e., the air is a complete hot gas, the ground has no thermal layer, and a large amount of dust is absent from the air).
q s = 2.45 × Δ p m 2 7.2 + Δ P m
where qs is the peak value of dynamic pressure, kPa.
Through calculation, Table 7 shows the dynamic pressure peaks at several measuring points on the blasting face of the two experiments.
The closer the distance between the dynamic pressure peak on the blasting face of the explosion-proof wall and the explosive, the greater the dynamic pressure peak. In accordance with the overpressure time history curve, each time corresponds to the overpressure of a wave front, and the dynamic pressure time history curve can be drawn as shown in Figure 11 and Figure 12.
Table 7 and Figure 11 and Figure 12 show that the peak value of dynamic pressure at measuring points 1, 19, and 20 on the explosion-facing surface was approximately two times the peak overpressure. From the overpressure and dynamic pressure time history curves, the dynamic pressure peak decreased with a decrease in the overpressure peak in the positive phase of overpressure. In the two experiments, when shock wave overpressure at measuring points 19 and 20 entered the negative phase, the air molecules behind the shock wave front possessed inertia and would not immediately drop into the negative phase. Measuring point 1 was located on the explosion-facing surface of the explosion-proof wall. It would not enter the negative phase immediately with large shock wave reflection overpressure, and impact fluctuating pressure would also have a long positive phase duration.

4. Numerical Simulation

4.1. Simulation Design

The diffraction angle range of this test was set to 22°–44°. The corresponding explosion distance can be obtained in accordance with the geometric relationship. The effects of the explosion-proof wall on shock wave diffraction under five proportional explosion distances and diffraction angles were calculated via numerical simulation. The specific layout is provided in Table 8.

4.2. Material Model

One of the simplest forms of an equation of state is that for an ideal polytropic gas which may be used in many applications involving the motion of gases. The state equation of air is approximated using the ideal gas state equation, and the specific expression is as follows:
P = γ 1 ρ ρ 0 E 0
where p is the air pressure, kPa; ρ is the air density after compression or expansion, g/cm−3; ρ0 is the initial air density, g/cm−3; γ is the adiabatic index; and E0 is the initial specific internal energy of air, kJ/m3. The main parameters of air are shown in Table 9.
TNT was described using the JWL state equation, and its specific form is as follows:
P = C 1 1 ω r 1 v e r 1 v + C 2 1 ω r 2 v e r 2 v + ω E 0 v
where P is the pressure, kPa; e is the internal energy of the explosive, kJ/m3; v is the relative volume of the explosive, cm3/g; E0 is the initial internal energy, kJ/m3; and C 1 , r 1 , C 2 , r 2 , and ω are the five constants of the state equation. The main parameters of TNT are shown in Table 10.
Sand was described using the compaction nonlinear state equation, the Drucker–Prager Strength Model, and the Tensile Stress Failure Model for Blast Resistance Simulation of Blast Walls. The material parameters are listed in Table 2. The specific form is as follows:
K ρ = ρ d P d ρ
where P is the current pressure, kPa; and ρ is the density of the material under zero pressure, g/cm3. The main parameters of sand are shown in Table 11 and the relationship between sand pressure change and density is shown in Figure 13.

4.3. Model Size and Meshing

In accordance with the design model size of the field experiment, the size of the air model was 22 m × 10 m × 5 m, the size of the sand model was 1 m × 5 m × 2.1 m, and the height of the medicine bag was 1.2 m above the ground. The TNT equivalent was consistent with the physical experiment. The air grid was 5 cm, and the explosion-proof wall grid was 10 cm. The 1/2 model was used for calculation because of the large model and its symmetry. The established model is shown in Figure 14.
To save calculation time and improve calculation accuracy, the unique remap technology of Autodyn software (2019R3).was used. This technology maps the calculation results obtained after a certain period in 1D and 2D onto another 3D model. This process not only considerably reduces calculation time, but also ensures the accuracy of the calculation results of the explosion problem. The specific process is illustrated in Figure 15. In the simulation process adopted in the current study, the remap mapping range was 90 cm, and the grid size was 0.1 cm.

4.4. Calculation Results and Analysis

Figure 16 shows the pressure cloud diagram of the shock wave diffracted jet flow field at different moments for five diffraction angles. At 0.5 ms, the two groups with a closer explosion distance had come into contact with the explosion-proof wall and formed a reflected wave. The explosion shock wave at 1–2 ms was reflected on the ground to form a Mach wave and advanced toward the explosion-proof wall. The reflection of the explosion-proof wall formed a downward and upward Mach wave. The downward Mach wave, the Mach wave formed by ground reflection, and the incident wave formed a superimposed wave at the wall foot, indicating that the wall foot was subjected to considerable shock wave pressure and continuous dynamic pressure. At 2.5–4.5 ms, the incident wave of the explosion shock wave was less than the reflected wave pressure formed on the wall. The incident wave continued to propagate forward along the top of the explosion-proof wall, and the sparse wave flowed to the high-pressure area, forming a cyclone on the top of the wall. The diffraction wave continued to propagate forward due to its interaction with the incident wave. At 12.5 ms, the diffracted wave was also reflected when it finally came into contact with the ground behind the wall and formed a Mach wave, which continued to propagate until it dissipated. The sequence of the main explosion shock wave on the front face of the explosion-proof wall is as follows: (1) incident wave of the explosion shock wave, (2) Mach wave reflected from the ground, and (3) superimposed wave of the Mach wave at the foot of the wall formed by the downward Mach wave reflected by the wall and the ground reflected wave. An inference can be made that if the explosion-proof wall is sufficiently far from the explosion source, then only the Mach wave reflected by the ground can reach the explosion-proof wall.
Figure 17 presents a comparison between the experimental and numerical simulation results when the diffraction angle was 32°. The simulation results are in good agreement with the experimental results, as shown in Figure 17. The experiment was performed on grassland and not on rigid ground, while the numerical calculation was conducted under ideal conditions without external interference and on rigid ground. An excessively large overpressure peak calculated during the simulation is normal. The pressure time history curve obtained by the experiment and the simulation indicated that the number and trend of each pressure peak were relatively consistent and could be used as the result of the trend analysis. The incident and reflected shock waves were superimposed significantly on the back side of the blast wall, and the peak pressure of the shock wave pressure behind the blast wall was tested at about 91% to 92% of the numerically calculated value, as shown in Figure 17.
When the explosive shock wave impacted the explosion-proof wall with different diffraction angles, the shock wave would generate varying pressures at different positions of the explosion-proof wall. The overpressure peak values calculated at different locations are provided in Table 12. A steel wire skeleton geotextile pocket was not added to the sand wall. Therefore, compactness was reduced and part of the pressure penetrated the sand wall, resulting in a small overpressure value at measuring point 1, but it minimally affected the study of the law of reflection overpressure at other measuring points behind the wall. The reduction of the shock wave by the explosion-proof wall is shown in Figure 18 and Table 13 (shock wave reduction Q 1 = P 19 P 13 , shock wave reduction Q 2 = P 20 P 17 , shock wave reduction rate η 1 = P 19 P 13 / P 19 , shock wave reduction rate η 2 = P 20 P 17 / P 20 , where the distances from P 19 and P 13 and from P 20 and P 17 to the explosion source were equal). This type of explosion-proof wall exhibits a good protection effect on the explosion shock wave. The reduction rate was more than 70% at the proportional distance position of the set measuring point.
Figure 19 shows that shock wave pressure changed at different positions behind the wall at the five diffraction angles. The maximum overpressure of the ground measuring point behind the explosion-proof wall appeared at measuring point 6 (approximately two times the height of the explosion-proof wall), which initially increased and then decreased, with two times the height of the explosion-proof wall as the apex. It was similar to the linear decrease of the quadratic function. The maximum overpressure of the aerial measuring point behind the explosion-proof wall appeared at measuring point 12 (approximately one times the height of the explosion-proof wall). After this distance, the quasi-linearity of the overpressure value continued to decrease. As the diffraction angle decreased, the pressure of the shock wave diffracted behind the wall also decreased. From Figure 19, a prediction can be made that as the measuring point moves farther away from the explosion-proof wall, the overpressure peaks will overlap regardless of whether the measuring point is in the air or on the ground. Therefore, shock wave pressure at different positions behind the wall was fitted.
By fitting the shock wave pressure at different positions on the ground behind the explosion-proof wall via a quadratic polynomial as shown in Figure 20a, the prediction formula (7) of overpressure position formed by the explosion shock wave with different diffraction angles on the ground behind the explosion-proof wall was obtained. By plotting the obtained overpressure prediction formula as shown in Figure 20b, most curves will intersect at 8.29–9.40 m behind the wall (approximately 4–4.5 times the height of the explosion-proof wall). Regardless of how the diffraction angle of the shock wave and the explosion proportion distance change, the pressure here will be roughly the same. Thereafter, shock wave pressure decreased rapidly. It can be considered a relatively safe position, which we call a relatively safe area, and the place that continues to remain far from the explosion source can acquire higher security.
y 1 = u x 2 + w x + z
where y 1 is the peak value of overpressure at x between the ground measurement point and the back surface of the wall under the current shock wave diffraction angle and proportional distance, kPa; x is the distance from the measuring point to the back surface of the wall, m; and u , w , and z are the constants related to the shock wave diffraction angle and the position of the measuring point, 0.20387 u 0.09137 , 0.72935 w 1.53256 , and 8.60411 z 10.47339 .
Linear fitting was performed for the points in the air of the explosion-proof wall after one times the height of the explosion-proof wall, as shown in Figure 21a. The prediction formula (8) of the overpressure position formed by the explosion shock wave with different diffraction angles on the ground behind the explosion-proof wall was obtained. By plotting the obtained overpressure prediction formula as shown in Figure 21b, most of the curves intersected at 7.42–8.60 m behind the wall (approximately 3.5–4 times the height of the explosion-proof wall). Similarly, it is called a relatively safe area.
y 2 = k x + h
where y 2 is the peak value of overpressure at x between the air measurement point and the back surface of the wall under the current shock wave diffraction angle and proportional distance, kPa; x is the distance from the measuring point to the back surface of the wall, m; and k and h are the constants related to the shock wave diffraction angle and the position of the measuring point, respectively, with 1.22714 k 0.72143 and 10.49000 k 14.73714 .

5. Conclusions

In the current study, an explosion-proof wall was set up to test a 1.5 kg TNT aerial explosion and perform a numerical simulation of the diffraction angle change. The conclusions drawn are as follows.
(1)
When the diffraction angle was 32°, the maximum overpressure on the ground and in the air was approximately 2.4 times the height of the explosion-proof wall. When the diffraction angle was 28°, the maximum diffraction overpressure was approximately 2.4–2.6 times the height of the explosion-proof wall from the explosion-proof wall, and the maximum reflection overpressure was approximately 1.5 times the height of the explosion-proof wall from the explosion-proof wall. Overpressure behind the wall initially increased and then decreased. The greater the diffraction angle, the faster the attenuation speed of the diffraction overpressure in the air behind the wall increased with distance from the explosion-proof wall.
(2)
Two evident overpressure peaks occurred in the pressure time history curve of the explosion shock wave. The time between the two pressure peaks Δt is a key factor that affects the destructive effect of the explosion shock wave. If Δt is less than the normal phase time T+ in the pressure time history curve of the diffracted shock wave, then the damage capability of the shock wave will be significantly improved. The impulses of the two overpressure peaks rise rapidly, and the impulses are superimposed rapidly, exhibiting a strong destructive effect.
(3)
By using the method for calculating the typical shock wave pressure in the ideal state, the peak value of dynamic pressure at three measuring points on the blast face was approximately two times that of overpressure. The air molecules behind the shock wave front possessed inertia, and dynamic pressure would have a longer positive duration than overpressure.
(4)
By using Autodyn software to simulate different diffraction angles of the shock wave, the maximum overpressure of the ground measuring points behind the explosion-proof wall appeared at approximately two times the height of the explosion-proof wall, and the maximum overpressure of the air measuring points behind the explosion-proof wall appeared at approximately one times the height of the explosion-proof wall. The overpressure prediction formula of the air and ground behind the wall was obtained by fitting shock wave pressure at different positions behind the wall. The relative safety area of the ground behind the wall was approximately 4–4.5 times the height of the explosion-proof wall. The relative safety area of the air behind the wall was approximately 3.5–4 times the height of the explosion-proof wall.
(5)
The relative safety area on the ground behind the wall is about 4 to 4.5 times the height of the blast wall, and the relative safety area in the air behind the wall is about 3.5 to 4 times the height of the explosion-proof wall.

Author Contributions

D.X. (Conceptualization, methodology, Writing—original draft, Writing—review & editing, Funding acquisition, Investigation, Supervision); W.Y. (Conceptualization, Formal analysis, methodology, Writing—original draft, Investigation); M.L. (methodology, software, Resources, Project administration); X.L. (Xiaoming Lü) (software, Project administration); K.L. (data curation, Visualization); J.Z. (software); X.L. (Xiaoshuang Li) (data curation); Y.L. (Visualization). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Open Foundation of the Shaanxi Key Laboratory of Safety and Durability of Concrete Structures (No. SZ02201).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge the Shock and Vibration of Engineering Materials and Structures Key Laboratory of Sichuan Province and the Shaanxi Key Laboratory of Safety and Durability of Concrete Structures.

Conflicts of Interest

The manuscript has not been published before and is not being considered for publication elsewhere. We declare there is no conflict of interest.

References

  1. Xu, R.; Chen, L.; Fang, Q.; Zheng, Y.; Li, Z.; Cao, M. Protective effects of gabion wall against blast waves from large tnt-equivalent explosions. Eng. Struct. 2021, 249, 113389. [Google Scholar] [CrossRef]
  2. Xu, W.; Wang, C.; Yuan, J.; Goh, W.; Li, T. Investigation on energy output structure of explosives near-ground explosion. Def. Technol. 2020, 16, 290–298. [Google Scholar] [CrossRef]
  3. Scherbatiuk, K.; Rattanawangcharoen, N. A hybrid rigid-body rotation model with sliding for calculating the response of a temporary soil-filled wall subjected to blast loading. Int. J. Impact Eng. 2011, 38, 637–652. [Google Scholar] [CrossRef]
  4. Scherbatiuk, K.; Rattanawangcharoen, N.; Pope, D.J.; Fowler, J. Generation of a pressure–impulse diagram for a temporary soil wall using an analytical rigid-body rotation model. Int. J. Impact Eng. 2008, 35, 530–539. [Google Scholar] [CrossRef]
  5. Remennikov, A.M.; Rose, T.A. Predicting the effectiveness of blast wall barriers using neural networks. Int. J. Impact Eng. 2007, 34, 1907–1923. [Google Scholar] [CrossRef]
  6. Scherbatiuk, K.; Rattanawangcharoen, N. Experimental testing and numerical modeling of soil-filled concertainer walls. Eng. Struct. 2008, 30, 3545–3554. [Google Scholar] [CrossRef]
  7. Li, J.; Ma, T.; Huang, C.; Huang, X.; Kang, Y.; Long, Z.; Liu, M. Protective mechanism of helmet under far-field shock wave. Int. J. Impact Eng. 2020, 143, 103617. [Google Scholar] [CrossRef]
  8. Mostert, F.J. Challenges in blast protection research. Def. Technol. 2018, 14, 426–432. [Google Scholar] [CrossRef]
  9. Türker, L. Thermobaric and enhanced blast explosives (tbx and ebx). Def. Technol. 2016, 12, 423–445. [Google Scholar] [CrossRef] [Green Version]
  10. Zhou, X.Q.; Hao, H. Prediction of airblast loads on structures behind a protective barrier. Int. J. Impact Eng. 2008, 35, 363–375. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, X.; Ding, Y.; Shi, Y. Numerical simulation of far-field blast loads arising from large tnt equivalent explosives. J. Loss Prev. Process Ind. 2021, 70, 104432. [Google Scholar] [CrossRef]
  12. Zhang, Q.L.; Zhang, Y.; Nian, X.Z. Effect of concrete protective wall on explosion shock wave. J. Vib. Shock. 2013, 32, 192–197. [Google Scholar] [CrossRef]
  13. Chen, L.; Zhang, L.; Fang, Q.; Mao, Y. Performance based investigation on the construction of anti-blast water wall. Int. J. Impact Eng. 2015, 81, 17–33. [Google Scholar] [CrossRef]
  14. Chen, L.; Fang, Q.; Zhang, L.; Zhang, Y.; Chen, W. Numerical investigation of a water barrier against blast loadings. Eng. Struct. 2016, 111, 199–216. [Google Scholar] [CrossRef]
  15. Li, X.; Wu, G.; Jia, H. The numerical simulation of explosion shock wave propagation subject to the retaining wall. Eng. Mech. 2012, 29, 245–250. [Google Scholar] [CrossRef]
  16. Shi, J.; Zhu, Y.; Chen, G.; Fu, J.; Liu, S. Assessment of blast resistance capacities of corrugated blast walls based on the P-I Model. J. Vib. Shock. 2017, 36, 188–195. [Google Scholar] [CrossRef]
  17. Li, Z.; Zhang, X.; Shi, Y.; Wu, C.; Li, J. Finite element modeling of frp retrofitted rc column against blast loading. Compos. Struct. 2021, 263, 113727. [Google Scholar] [CrossRef]
  18. Shi, Y.; Jiang, R.; Li, Z.; Ding, Y. A new macromodel of beam-column joints for progressive collapse analysis of reinforced concrete structures under blast loading. Int. J. Struct. Stab. Dyn. 2022, 22, 2250008. [Google Scholar] [CrossRef]
  19. Chaomin, M.; Huiqi, R.; Yongchi, L.; Kai, X. Experimental Study of Blast Wave Reflection and Diffraction on a Shelter Wall. J. Exp. Mech. 2008, 23, 169–174. [Google Scholar]
  20. Zhu, F. Research on Anti-Blast Performance of Portable Explosion-Proof Walls in the Field Environment; North University of China: Taiyuan, China, 2017. [Google Scholar]
  21. Mu, C.; Zhou, H.; Ma, H. Prediction method for ground shock parameters of explosion in concrete. Constr. Build. Mater. 2021, 291, 123372. [Google Scholar] [CrossRef]
  22. Ge, T.; Cao, H.; Zhang, Z. Wave suppression performance tests of rapidly assembled anti-blast walls. J. Vib. Shock. 2019, 38, 239–244. [Google Scholar] [CrossRef]
  23. Ren, H.Q.; Huang, K.; Wu, X.Y. Research progress of dynamic pressure damage to ground targets by air shock wave. Prot. Eng. 2021, 43, 1–9. [Google Scholar] [CrossRef]
  24. Zhou, H.; Ren, H.; Yu, S.; Chen, J.; Cao, Y. Review on measurement technology of dynamic pressure of air shock wave. Prot. Eng. 2019, 41, 61–69. [Google Scholar]
  25. Ma, J.; Chen, J.; Guan, J.; Lin, Y.; Chen, W.; Huang, L. Implementation of Johnson-Holmquist-Beissel model in four-dimensional lattice spring model and its application in projectile penetration. Int. J. Impact Eng. 2022, 170, 104340. [Google Scholar] [CrossRef]
  26. Lin, Y.; Ma, J.; Lai, Z.; Huang, L.; Lei, M. A FDEM approach to study mechanical and fracturing responses of geo-materials with high inclusion contents using a novel reconstruction strategy. Eng. Fract. Mech. 2023, 282, 109171. [Google Scholar] [CrossRef]
  27. Moszynski, J.R. The dynamics of explosion and its use. Nucl. Technol. 1983, 60, 167. [Google Scholar] [CrossRef]
Figure 1. Shock wave testing system.
Figure 1. Shock wave testing system.
Sustainability 15 11164 g001
Figure 2. Diagram showing the layout of shock wave measuring points.
Figure 2. Diagram showing the layout of shock wave measuring points.
Sustainability 15 11164 g002
Figure 3. Physical layout of shock wave measure points.
Figure 3. Physical layout of shock wave measure points.
Sustainability 15 11164 g003
Figure 4. Diffraction angle 32°: time history curve of typical shock wave pressure: (a) Point 9; (b) Point 14; (c) Point 19; (d) Point 20.
Figure 4. Diffraction angle 32°: time history curve of typical shock wave pressure: (a) Point 9; (b) Point 14; (c) Point 19; (d) Point 20.
Sustainability 15 11164 g004
Figure 5. Diffraction angle 32°: overpressure distribution behind the wall: (a) ground; (b) air.
Figure 5. Diffraction angle 32°: overpressure distribution behind the wall: (a) ground; (b) air.
Sustainability 15 11164 g005
Figure 6. Diffraction angle 28°: overpressure distribution behind the wall: (a) ground; (b) air.
Figure 6. Diffraction angle 28°: overpressure distribution behind the wall: (a) ground; (b) air.
Sustainability 15 11164 g006
Figure 7. Interval time between the incident and reflected waves behind the wall: (a) diffraction angle 32°; (b) diffraction angle 28°.
Figure 7. Interval time between the incident and reflected waves behind the wall: (a) diffraction angle 32°; (b) diffraction angle 28°.
Sustainability 15 11164 g007
Figure 8. Impulse calculation rules.
Figure 8. Impulse calculation rules.
Sustainability 15 11164 g008
Figure 9. Diffraction angle 32°: typical impulse time history curve: (a) Point 9; (b) Point 14; (c) Point 19; (d) Point 20.
Figure 9. Diffraction angle 32°: typical impulse time history curve: (a) Point 9; (b) Point 14; (c) Point 19; (d) Point 20.
Sustainability 15 11164 g009
Figure 10. Comparison chart of air and ground impulse; (a) diffraction angle 32°; (b) diffraction angle 28°.
Figure 10. Comparison chart of air and ground impulse; (a) diffraction angle 32°; (b) diffraction angle 28°.
Sustainability 15 11164 g010
Figure 11. Diffraction angle 32°: overpressure and dynamic pressure: (a) Point 1; (b) Point 19; (c) Point 20.
Figure 11. Diffraction angle 32°: overpressure and dynamic pressure: (a) Point 1; (b) Point 19; (c) Point 20.
Sustainability 15 11164 g011
Figure 12. Diffraction angle 28°: overpressure and dynamic pressure: (a) Point 1; (b) Point 19; (c) Point 20.
Figure 12. Diffraction angle 28°: overpressure and dynamic pressure: (a) Point 1; (b) Point 19; (c) Point 20.
Sustainability 15 11164 g012
Figure 13. The relationship between sand density and pressure change.
Figure 13. The relationship between sand density and pressure change.
Sustainability 15 11164 g013
Figure 14. Establishment of the computational model.
Figure 14. Establishment of the computational model.
Sustainability 15 11164 g014
Figure 15. Remap mapping technology: (a) 1D explosion wave propagation; (b) 2D explosion wave propagation; (c) remap onto a 3D model.
Figure 15. Remap mapping technology: (a) 1D explosion wave propagation; (b) 2D explosion wave propagation; (c) remap onto a 3D model.
Sustainability 15 11164 g015
Figure 16. Pressure cloud of different diffraction angles: (a) 44°; (b) 38°; (c) 32°; (d) 28°; (e) 22°.
Figure 16. Pressure cloud of different diffraction angles: (a) 44°; (b) 38°; (c) 32°; (d) 28°; (e) 22°.
Sustainability 15 11164 g016aSustainability 15 11164 g016b
Figure 17. Comparison between the experiment and numerical simulation: (a) Point 4 and (b) Point 5.
Figure 17. Comparison between the experiment and numerical simulation: (a) Point 4 and (b) Point 5.
Sustainability 15 11164 g017
Figure 18. Shock wave attenuation of the explosion-proof wall with a diffraction angle of 44°: (a) shock wave attenuation 1; (b) shock wave attenuation 2.
Figure 18. Shock wave attenuation of the explosion-proof wall with a diffraction angle of 44°: (a) shock wave attenuation 1; (b) shock wave attenuation 2.
Sustainability 15 11164 g018
Figure 19. Shock wave pressure at different positions behind the wall: (a) ground behind the wall; (b) air behind the wall.
Figure 19. Shock wave pressure at different positions behind the wall: (a) ground behind the wall; (b) air behind the wall.
Sustainability 15 11164 g019
Figure 20. Safety zone prediction: (a) polynomial fitting; (b) safety point.
Figure 20. Safety zone prediction: (a) polynomial fitting; (b) safety point.
Sustainability 15 11164 g020
Figure 21. Safety zone prediction: (a) polynomial fitting; (b) safety point.
Figure 21. Safety zone prediction: (a) polynomial fitting; (b) safety point.
Sustainability 15 11164 g021
Table 1. Relationship between proportional distance and explosion distance.
Table 1. Relationship between proportional distance and explosion distance.
Test GroupProportional Distance/m·kg−1/3Burst Distance/mDiffraction Angle α/°
1-11.261.4432
1-21.481.6928
Table 2. Density and moisture content of sand.
Table 2. Density and moisture content of sand.
SandDensity/g·cm−1/3Water Content/%
12.01735.16
22.03331.06
32.03333.10
Table 3. Diffraction angle 32°: shock wave pressure parameters.
Table 3. Diffraction angle 32°: shock wave pressure parameters.
Measuring PointsPeak Value of Diffraction Wave/kPaPeak Value of Reflected Wave/kPaΔt/ms
13845.502628.310.18
2254.3828.731.23
316.471.516.05
411.2313.303.92
511.2513.202.91
610.4614.303.28
718.5119.693.20
813.8919.002.92
912.1517.922.52
108.4919.532.89
1122.4910.656.46
1220.6210.516.18
1320.6719.394.84
1412.9017.004.68
1512.8722.754.39
1610.0313.233.60
178.8113.343.36
186.8713.183.13
1931.576.984.81
2015.061.3614.96
Note: Interval between the reflected and diffracted waves, Δ t = t 2 t 1 .
Table 4. Diffraction angle 28°: shock wave pressure parameters.
Table 4. Diffraction angle 28°: shock wave pressure parameters.
Measuring PointsPeak Value of Diffraction Wave/kPaPeak Value of Reflected Wave/kPaΔt/ms
13754.102421.070.24
2271.20169.212.97
3---
49.43--
58.5011.324.46
610.6610.453.19
713.8614.293.62
813.9714.893.43
910.6715.602.93
109.2714.122.85
1120.9210.096.68
1218.5311.055.76
1319.7920.144.86
1414.3715.354.74
1513.0618.994.79
168.5212.643.98
178.3111.373.72
187.8410.263.53
1925.9710.465.67
2015.836.077.67
Table 5. Diffraction angle 32°: impulse.
Table 5. Diffraction angle 32°: impulse.
Measuring PointsPositive Impulse of
Diffraction Wave/Pa·s
Cumulative Impulse of Negative Reflection Wave Front/Pa·sMaximum Impulse/Pa·s
197.21217.52217.52
244.115.8647.16
316.8357.7957.79
415.2749.7449.74
515.0849.2149.21
615.1745.8745.87
718.8654.1954.19
813.9446.3946.39
915.1743.1943.19
1019.6760.6960.69
1122.9747.2347.23
1220.7149.1249.12
1323.7264.7764.77
1430.0857.1357.13
1512.2162.7162.71
1611.3147.9947.99
1710.8144.9844.98
189.4939.8439.84
1949.4154.6554.65
2035.662.4235.66
Table 6. Diffraction angle 28°: impulse.
Table 6. Diffraction angle 28°: impulse.
Measuring PointsPositive Impulse of
Diffraction Wave/Pa·s
Cumulative Impulse of Negative Reflection Wave Front/Pa·sMaximum Impulse/Pa·s
1189.27338.91338.91
235.14–39.1635.14
3---
411.57--
514.3441.4741.47
612.9333.0333.03
720.2350.9750.97
817.6945.9445.94
916.8342.1542.15
1013.6234.4334.43
1123.3250.3550.35
1214.2449.9049.90
1316.5568.2268.22
1416.4460.2860.28
1513.0055.7855.78
1610.1639.1639.16
179.0337.4537.45
187.7634.7934.79
1946.0047.8647.86
2030.5933.6233.62
Table 7. Dynamic pressure peak.
Table 7. Dynamic pressure peak.
Measuring
Points
Dynamic Pressure Peak/kPa
32°28°
19403.879179.94
1962.9849.81
2024.9626.67
Table 8. Relationship between diffraction angle and explosion distance.
Table 8. Relationship between diffraction angle and explosion distance.
Test GroupProportional Distance/m·g−1/3Burst Distance/mDiffraction Angle/°
1-10.810.9344
1-21.011.1538
1-31.261.4432
1-41.481.6928
1-51.952.2322
Table 9. Parameters for air.
Table 9. Parameters for air.
ParameterDensity/g·cm−1/3Adiabatic ExponentSpecific Heat/J·(kg·K)−1Specific Internal Energy/J·kg−1
Value0.0012251.4717.62.068 × 105
Table 10. Parameters for TNT.
Table 10. Parameters for TNT.
ParameterDensity/g·cm−1/3 C 1 /GPa C 2 /GPa r 1 r 2 ω VOD/m·s−1 E 0 /J·m−3 P C J /GPa
Value1.633.73 × 1023.744.150.90.3569306 × 10921
Table 11. Parameters for sand.
Table 11. Parameters for sand.
ParameterDensity/g·cm−1/3Shear Modulus/MPaHydro Tensile Limit/KPa
Value2.02876.9
Table 12. Peak pressure of test points in the numerical simulation.
Table 12. Peak pressure of test points in the numerical simulation.
Measuring PointsPeak Pressure/kPa
44°38°32°28°22°
12540.041680.39987.25604.49379.98
2110.49110.58103.7594.6474.99
311.7711.4710.8510.349.32
412.7212.2111.5810.959.60
513.2812.6811.9911.339.92
613.4112.8012.1111.4310.08
713.0412.4611.9011.2610.03
812.2011.6911.3110.769.70
911.2210.7710.4610.089.23
109.739.459.389.208.57
1112.0611.5110.539.658.13
1212.5712.0611.2010.488.99
1311.1110.6610.079.558.45
149.639.228.838.457.61
158.297.947.707.456.81
167.176.876.736.566.09
176.296.065.985.865.50
185.154.994.984.954.73
1965.3760.0254.7351.1444.06
2024.6223.8422.5321.8420.06
Table 13. Shock wave attenuation.
Table 13. Shock wave attenuation.
Diffraction Angle Q 1 /kPa η 1 /% Q 2 /kPa η 2 /%
44°54.2683.018.3374.5
38°49.3682.217.7874.6
32°44.6681.616.5573.5
28°41.5981.315.9873.2
22°35.6180.814.5672.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xiao, D.; Yang, W.; Lin, M.; Lü, X.; Liu, K.; Zhang, J.; Li, X.; Long, Y. Effects of an Explosion-Proof Wall on Shock Wave Parameters and Safe Area Prediction. Sustainability 2023, 15, 11164. https://doi.org/10.3390/su151411164

AMA Style

Xiao D, Yang W, Lin M, Lü X, Liu K, Zhang J, Li X, Long Y. Effects of an Explosion-Proof Wall on Shock Wave Parameters and Safe Area Prediction. Sustainability. 2023; 15(14):11164. https://doi.org/10.3390/su151411164

Chicago/Turabian Style

Xiao, Dingjun, Wentao Yang, Moujin Lin, Xiaoming Lü, Kaide Liu, Jin Zhang, Xiaoshuang Li, and Yu Long. 2023. "Effects of an Explosion-Proof Wall on Shock Wave Parameters and Safe Area Prediction" Sustainability 15, no. 14: 11164. https://doi.org/10.3390/su151411164

APA Style

Xiao, D., Yang, W., Lin, M., Lü, X., Liu, K., Zhang, J., Li, X., & Long, Y. (2023). Effects of an Explosion-Proof Wall on Shock Wave Parameters and Safe Area Prediction. Sustainability, 15(14), 11164. https://doi.org/10.3390/su151411164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop