Study on the Earth-Covered Magazine Models under the Internal Explosion

. Te use of earth-covered magazines (ECMs) is increasingly prevalent in protective engineering due to their concealment and cost-efectiveness. To explore the optimal thickness of earth covering for ECMs, scaled model tests were conducted under explosive charges equivalent to 30kilograms of TNT. Te resulting overpressure outside the model in the 180 ° direction was measured. Subsequently, computational analyses were conducted employing LS-DYNA software to examine these experimental fndings. Te fndings indicate that increasing the thickness of the rear soil can mitigate peak overpressure, delay the air shock wave’s arrival time, and reduce the impulse of the positive phase. Te numerical calculations closely align with experimental data, with peak overpressure deviation remaining under 10%. Te shock wave initially impacts the top of the model before reaching the rear, with soil scattering more pronounced in the 90 ° direction compared to the 180 ° direction. Furthermore, an analysis of soil energy absorption rate variation was conducted based on energy conservation principles. Tese results provide valuable insights for optimizing the design and construction of ECMs.


Introduction
Also referred to as a semiunderground storage, this structure features top, sides, and rear coverage with soil, while the front wall is equipped with entrances and exits.It is commonly utilized for storing explosive materials and ammunition [1].With urban areas expanding continuously, adjacent buildings often surround ammunition depots.Given the occasional internal explosion risks within ammunition depots and the insufcient safety distance between neighboring structures, it becomes imperative to devise methods for minimizing the shock wave impact radius.Such measures ofer crucial direction for designing weapons and ammunition storage facilities.
Experimental studies on earth-covered magazines trace back to 1945, notably with the Igloo Magazine tests conducted in Arco, Idaho.Tese endeavors aimed to establish safe distances while conserving land resources [2,3].Tese tests distinctly focused on examining the overall structural response, air explosion measurements, and the density and distribution patterns of debris.Determining the closest safe distance between the earth-covered magazines was guided by assessing debris dispersion and overpressure peaks resulting from internal explosion tests.
Existing research suggests that the strength and thickness of the earth covering are the crucial factors in mitigating external shock wave.Kingery et al. [4] conducted 1 : 50 scaled internal explosion tests on earth-covered magazine models, meticulously examining the impact of varying soil cover thickness above the model on the external air pressure in the 0 °, 90 °, and 180 °directions.Te fndings reveal that augmenting the soil layer thickness does not lead to increased air pressure in the 0 °direction, while the alteration pattern of air pressure in the 90 °remains inconspicuous, and the decay of air pressure in the 180 °direction intensifes.Nonetheless, Charles' experimental inquiry primarily focuses on the thickness of earth cover above the magazine, with scant exploration into the infuence of changes in earth cover thickness behind the magazine on external pressure and positive pressure impulse.Moreover, Charles' experimental results indicate that the pressure sensors were susceptible to high temperatures and light, resulting in numerous missing experimental data points.Consequently, in this study, to reduce the infuence of light and heat on the pressure sensors, protective measures include shielding the sensors with aluminum flm and applying silicone grease.
Te United States Department of Defense (DOD) standard [5] highlights a pressure reduction of approximately 66.7% within an earth-covered magazine equipped with protective shielding.Similarly, according to the Chinese Standard [6], installing protective barriers around the earthcovered magazine is recommended to achieve a maximum reduction in external pressure of 60%.Beppu et al. [7] conducted experiments on earth-encased magazines with protective barriers on their front walls, observing that increased soil thickness alters the distribution of external pressure.As the depth of the soil increases, the intensity of the shock wave decreases at the rear.Liu et al. [8] and Li et al. [9] analyzed the factors infuencing the safety distance of earth-covered magazine, emphasizing the mitigating efect of soil on shock wave overpressure and its efcacy in blocking structural debris pathways.Furthermore, Xu et al. [10,11] conducted scaled tests on earth-covered magazines to investigate overpressure and debris distribution in various directions.Tese fndings were integrated with numerical simulations to estimate the distance that concrete fragments would project.
Te German military [12] conducted an internal explosion tests within earth-covered ISO containers to ascertain safe distances and debris ranges between adjacent containers.Results revealed that pressure relief and debris distribution predominantly occurred in the container's frontal area, with a notable reduction in debris range facilitated by soil mounding at the container's front.Tian and Deng [13] explored the infuence of earth covering thickness on structural dynamic responses and shock wave overpressure propagation in underground structures subjected to internal explosions.Findings indicate that earth covering thickness impacts the structure's dynamic response to some extent, while shock wave propagation during internal explosions exhibits complexity, with overpressure diminishing as the proportional distance from the blast center increases.Wang [14] investigated the impact of ground blast shock waves on soil dynamic properties using a combination of experiments and numerical simulations.Destructive efects leading to the deformation and destruction of external protective soil layers were analyzed.
Numerous scholars [15][16][17][18][19] have extensively studied the dynamic mechanical characteristics and energy dissipation behaviors of soils under varying strain rates, utilizing a combination of experiments and numerical simulations.Liu and Chen [20] examined the specifc energy absorption of compacted cohesive soils under diferent strain rates and water contents, employing SHPB (split Hopkinson pressure bar) test equipment.Ma et al. [21] investigated the energy dissipation properties of frozen sandy soils across various deformation stages via split Hopkinson pressure bar tests, revealing close correlations between energy dissipation parameters and moisture content, temperature, and strain rate.Alekseenko and Rykov [22] analyzed the impact of detonation location on the distribution of blast energy in soil, employing physical parameters such as axial stress and impulse as indicators.Grujicic et al. [23] and Grujicic et al. [24] investigated the explosion process within soils under varying water content conditions, combining experimental approaches with numerical simulations.Teir fndings indicate that low porosity and high saturation soils exhibit efective energy absorption capabilities.
In summary, while numerous scholars have primarily concentrated their research on internal pressure distribution and structural response resulting from internal explosions, there remains a paucity of studies on the efects of internal explosions within earth-covered magazines on induced damage to external surroundings, particularly concerning the propagation of pressure externally and the impact of earth covering thickness.In this work, internal explosion tests were performed using four scaled models of earthcovered magazines.Pressure data in 180 °direction and soil scattering were acquired utilizing pressure sensors, highspeed camera, and unmanned aerial vehicles (UAVs).Additionally, the simulation of internal explosions within the models was conducted using the LS-DYNA software.Te study evaluated the infuence of soil on explosive shock wave energy absorption and proposed the most suitable design value for rear soil thickness.Te fndings ofer valuable insights for engineering applications, providing guidance for optimizing earth-covered magazine designs and enhancing their protective capabilities.

Experimental Set-Up
2.1.Models.Swisdak et al. [25] conducted explosion tests in earth-covered magazines, revealing that both reinforced concrete and masonry structures are destroyed when subjected to explosions with large charges.Additionally, Xiong's research [26] demonstrated that under high loading densities of 30 kg/m 3 , all reinforced concrete (RC) structures disintegrate, with minimal impact on structural integrity from the soil cover under such conditions.Terefore, in this study, structural strength considerations are disregarded, and high-density polyethylene (HDPE) pipes are utilized to represent the most unfavorable scenario.Half-round HDPE pipes were employed to construct circular arches, serving as the support system for the earth-covered magazine.Tese pipes are covered with sandy soil above and behind them to replicate the soil's attenuation efect on shock waves in the absence of concrete, representing the most unfavorable condition.It is important to note that HDPE pipe material possesses low strength and is directly damaged by explosive forces, thus exerting negligible infuence on shock wave propagation in the soil.Te sand utilized ranges 4.55 × 10 −5 ∼2 × 10 −5 m, with a density of 1.78 g/cm 3  Yang's research [28] indicates that the primary infuence on overpressure from internal explosions is determined from the following parameters: (1) Air properties: initial pressure (p 0 ), air density (ρ a ), air adiabatic index (c a ); (2) Explosive characteristics: mass (W), density (ρ e ), energy released per unit mass (e), and explosion product expansion index (c e ); (3) Structure attributes: length (l), width (b), height (h), strength (σ), Young's modulus (E), and material density (ρ s ); (4) Te distance (R) from the measurement point to the source of the explosion.
Citing the principle of similarity [29], the fundamental dimensions of mass (M), length (L), and time (T) are employed as the primary dimensions, as detailed in Table 1, to encapsulate the dimensional characteristics of each parameter.Te relationship between explosion overpressure inside the structure and its infuencing factors is described as follows: P � f W, p 0 , ρ a , c a , ρ e , e, c e , l, b, h, σ, E, ρ s , R . ( Given that E is not an independent factor for the same structural material, it can be omitted, along with c a and c e , which are dimensionless quantities.Hence, equation (1) simplifes to P � f W, p 0 , ρ a , ρ e , e, l, b, h, σ, ρ s , R . ( Revising equation ( 2) to incorporate W, ρ a , and σ as basic quantities, we obtain Under identical charge type and air conditions, with the structural material remaining consistent, if (W (1/3) Here, p refers to the prototype and m denotes the model, while λ represents the scaling ratio.
Hence, the dimensions of the test model structure were scaled down based on a prototype structure, illustrated in Figure 1.For security reasons, the specifcs of the prototype structure are undisclosed.In accordance with the stipulations outlined in the standard [6], the test incorporated four variations of rear soil thickness, detailed in Table 2.

Experimental Instruments and Equipment.
Te experimental instrumentation comprised pressure sensors, electricity provision setups, oscilloscopes, signal conditioning apparatus, and data collection systems.Figure 2 illustrates the detailed sensor layout diagram, utilizing Swiss KIS-TLER6233 pressure sensors.Tese sensors, with the resonance frequency of 500 kHz, featured a range of 0.17-6.9MPa.Signal acquisition from the pressure sensors occurred via signal acquisition occurred through the pressure sensors via the data acquisition tool, operating at a sampling rate of 1 MHz.Eleven sensors were strategically positioned along the 180 °direction at specifed distances from the blast center, as outlined in Table 3.Each sensor was securely fxed to a bracket via a casing pipe and horizontally adjusted for optimal positioning, as depicted in Figure 3.To minimize interference during testing, the bracket was buried 50 cm below the ground surface and secured with nails.An unmanned aerial vehicle (UAV) and a high-speed camera were used to document the test procedures, capturing vital data on soil scattering and shock wave propagation.
A box bulk TNT charge, with a density of 0.9 g/cm 3 and a total equivalent TNT mass of 30 kg, was selected as the explosive.To expedite the detonation process, an electrical initiator was employed to activate the RDX explosive column, subsequently causing the detonation of the TNT explosive.Positioned centrally, the explosive was depicted as shown in Figure 4.

Experiment Results
3.1.Pressure Parameters. Figure 5 illustrates the pressuretime histories of case number 1 recorded by the sensor, with data processing conducted to derive shock wave-related parameters as presented in Table 4 (blank sections indicate missing or abnormal data).Additionally, Figure 6 depicts a typical pressure-time history.Comparing the rear soil thickness of 30 mm to increased thickness of 60 mm, 150 mm, and 300 mm, when the scaled distance is small, pressure decreased by 19.07%, 25.04%, and 22.58%, respectively, while the shock wave arrival time was prolonged by 6.35%, 7.93%, and 5.86% and the positive phase impulse reduced by 7.62%, 6.61%, and 5.26%, respectively.For larger scaled distances, the pressure decreases by 36.34%,44.20%, and 37.71%, respectively, with essentially the same shock wave arrival time.Furthermore, the positive phase impulse decreases by 23.50%, 32.29%, and 17.43%, respectively.

Shock and Vibration
Using the peak overpressure equation for near-surface blast air summarized by Baker [30], the pressure-scaled distance and positive phase impulse-scaled distance were obtained from experimental data and plotted in double logarithmic coordinates (see Figures 7 and 8).As depicted in equations ( 1) and ( 2), each dataset was represented by a ftting function y � (A/x) + (B/x 2 ) + (C/x 3 ) with a confdence level exceeding 99%.Notably, for the thicknesses of 60 mm, 150 mm, and 300 mm, all exhibited a consistent trend in overpressure.However, upon comparing the impulse, it was found that in larger scaled distances, the minimum impulse occurred with a thickness of 60 mm.Hence, further analysis is warranted.
Shock and Vibration

Numerical Simulation
To investigate the impact of rear soil thickness on external pressure propagation, numerical simulation proves to be the prevalent and efcacious method.In this investigation, we utilized LS-DYNA fnite element simulation software, employing a fuid-structure interaction algorithm.Tis allowed us to calculate the temporal profles of air pressure outside the model and track the process of soil scattering.Shock and Vibration

Finite Element Model.
To explore the efects of explosions, numerical analysis tools must ensure coupling between hydrodynamic states and solids.Built upon the principles of continuum mechanics, LS-DYNA, a nonlinear fnite element analysis program, utilizes explicit time integration methods to capture dynamic mechanical behaviors with precision.Tis approach enables the study of instantaneous explosion dynamics, albeit requiring controlled time step (Δt) for stable calculations.
Efective employment of explicit time integration relies heavily on selecting an appropriate time step.Specifcally, ensuring a time step smaller than the critical time interval, which is determined by the time step control factor (TSSFAC), is essential for stability.LS-DYNA automatically partitions time (T) into T/Δt intervals.Te program manual recommends setting the time step control factor within the range of 0.1∼0.9, with a value lower than 0.67 suggested for explosion analysis [31].
where l mesh represents the minimum mesh dimension and C l denotes the wavelength.
For computational convenience, a symmetric structure was simulated, employing a half-model approach.Te software was used to construct models for air, explosives, pipes, and soil.Te numerical simulation involved a fuidstructure coupling algorithm, with the Euler mesh utilized for air and explosives, and the Lagrangian mesh employed for pipes.Tis method combines the benefts of the Lagrangian method and the Eulerian method, addressing material deformation and moving boundary complexities [32][33][34][35].During explosion, signifcant soil deformation occurs.Establishing a solid element mesh may result in negative volume generation, leading to substantial errors or halting the calculation.
Te soil characterization is established utilizing the keyword * INITIAL_VOLUME_FRACTION_GEOMETRY, facilitating modeling of complex structures and meshing.Te grid size of the soil is determined by * NTRACE, set to 17, resulting in subdivision of the ALE element into 353 regions.It is crucial to ensure the air region is sufciently large to capture a complete pressure-time curve when soil reaches the air boundary, causing a sudden decrease in step size.Nonrefection boundaries are established at air boundaries, except for the symmetry plane and ground.Te To ensure grid convergence, we chose six element sizes for the model, varying from 5 mm to 20 mm.Te soil element size matched that of the air element.Figure 10 illustrates the relative error between numerical simulation and experimental air pressure values at diferent earth covering thicknesses and air element sizes.As the element size decreases, the relative error diminishes, albeit with increased computational time.Balancing computational efciency and relative error, an element size of 12.5 mm is recommended, also applied to the circular pipe unit size.Te  Te characterization of air properties employs the * MAT_NULL material model combined with the EOS_LINEAR_POLYNOMIAL polynomial equation of state [31], which is formulated as In this equation, P represents the gas pressure and C 0 to C 6 denote the state parameters.For this model, Denoted as E, the internal energy of the gas is set at 2.5 × 10 5 J/m 3 , and V 1 represents the relative volume, which is assumed to be 1.0.
Te * MAT_HIGH_EXPLOSIVE_BURN model characterizes the explosive material, coupled with the EOS_JWL polynomial equation of state [31], which is formulated as In this equation, P represents the pressure, while the material parameters A, B, R 1 , and R 2 hold the values of 3.712 × 10 11 , 3.23 × 10 9 , 4.15, and 0.95, correspondingly.Te internal energy per unit volume, denoted by E, is prescribed as 7 × 10 9 J/m 3 , and V represents the relative volume, assumed to be 1.0.Shock and Vibration Te mechanical properties of HDPE material are viscoelastic, afected by stress, deformation, temperature, and time.Temperature efects on HDPE materials are disregarded in this study, indicating their mechanical properties are strain rate-dependent.Terefore, the * MAT_ PLASTICITY_ POLYMER model employed herein characterizes the dynamic response of polymers under high strain rates, with material parameters outlined in Table 5.
To accurately simulate the crushing or compacting behavior of the pore material and its dynamic response postexplosion, the soil constitutive law must be precise.Tis study utilizes the * MAT_SOIL_AND_FOAM_FAILURE model, tailored for soil and foam calculations, aligning with material calculations [31].Te material is constitutively modeled as follows: In this constitutive model, ϕ denotes the partial ideal plastic yield function, J 2 represents the second stress invariant, J 2 � (1/2)S ij S ij , S ij denotes the deviatoric stress tensor, and J 2 � (1/3)σ y 2 , where σ y is the uniaxial yield stress on the yield surface. a 0 , a 1 , and a 2 are constants; p denotes the pressure.In cases where the material's yield strength is disregarded, ϕ � 0. A summary of material parameters can be found in Table 6.

Numerical Results and Discussion
5.1.Overpressure.Te numerical methodologies, material parameters, and numerical models were validated by integrating them with experimental data.Ground explosion experiments were conducted to measure air overpressure, verifying the numerical analysis results and providing a reliable numerical model for ECM analysis.Due to the limited size of the air grid, data comparison is restricted to the frst measurement point.Figure 12 illustrates the comparison of pressure between experimental and simulation results at a scaled distance of 1.6 m/kg 1/3 .Te relative error (%) was calculated using the following formula: (numerical analysis value − experimental measurement value)/experimental measurement value × 100%.
Table 7 demonstrates that with the exception of specifc cases with signifcant positive pressure duration error, the remaining data errors are within 10%.Despite some

Shock and Vibration
discrepancies, these error values fall within acceptable limits.Te divergence between experimental fndings and numerical analysis results primarily stems from the inherent nonuniformity, nonlinearity, and dynamic fuctuations in materials, geometries, environmental conditions, and contact properties that are integral to the analysis of shock efects.Additionally, in fnite element analysis of instantaneous dynamics, it is imperative to ensure that material parameters are consistent with the evolving behavior observed in the material, which encompasses variations in internal energy, temperature, pressure, and volume as dictated by the EOS.Te meshed model also played a signifcant role in the numerical calculations.Overall, these fndings suggest good agreement between experimental and numerical models, afrming the suitability of the numerical model for further analysis.Moreover, the results validate the efectiveness of using the ALE fuid-structure coupling algorithm to examine the dynamic attributes of the explosion phenomenon [38][39][40][41].

Damage
Characteristics of the Pipe. Figure 13 illustrates the destruction process of the pipe with a soil cover thickness of 30 mm.Upon detonation of the explosive, a spherical blast wave rapidly forms, exerting pressure on the top and sides of the pipe at approximately 0.12 ms.As the soil's wave impedance is lower than that of the pipe, a segment of the compressive pulse reaching the pipe's surface is transmitted into the soil, while the remainder refects as a tensile wave.By 0.33 milliseconds, the tensile stress induced by this wave surpasses the pipe's dynamic tensile strength, resulting in cracks forming on both the upper and lateral surfaces.Over time, the blast wave's area of infuence expands, and deformation propagation shifts gradually from the pipe's center towards both ends.Te constrained elements at the pipe-soil interface experience signifcant deformation, culminating in failure as depicted in Figure 13    Shock and Vibration pipe exhibit the most severe damage, consistent with experimental fndings.

Soil Scattering.
Te soil is indirectly defned using the * INITIAL_VOLUME_FRACTION_GEOMETRY keyword, which notably streamlines modeling endeavors and alleviates limitations in mesh initialization by indicating the proportions of various materials within Euler background meshes [42].Te soil modeling process is depicted in Figure 14.Initially, the entire ALE part is flled with air, followed by replacing the material within the designated container with soil.Subsequently, the soil within the designated area is replaced with air.Notably, the CNTTYP value in the keyword is sequentially set to 5, 4, and 3 throughout the modeling process, defning the geometric shape of the container as cuboid, cylinder, and plane, respectively.Furthermore, the coordinates of the container's minimum and maximum points are input to determine its volume.For further details, consult the manual [31].
Figure 15 depicts the soil scattering process following the explosion.At 0.299 ms after detonation, deformation initiates at the top, with the top soil separating and dispersing at 0.899 ms, followed by deformation at the corner behind the top.Approximately 1.56 ms later, the soil at the corner is also dispersed and ejected, leading to deformation of the entire covering soil in all directions, reaching the air domain   consistent with high-speed photos (Figure 9) and validating the numerical simulation.
5.4.Shockwave Propagation.Figure 16 illustrates the process of shock wave transmission.Initially, emanating from the aperture, the shock wave encounters an unobstructed trajectory.Te shock waves reach the top and the rear at 0.09 ms and 0.18 ms, respectively.Subsequently, interaction with the soil causes refected waves to become apparent from 0.6 ms to 0.81 ms, as evidenced by the increased area of the red pressure zone in the fgure.At 1.50 ms, two distinct shock wave fronts at the top and corner, converging into a single wave front at 3.99 ms and propagating uniformly through the air.Finally, by 5.98 milliseconds, the shockwave extends to the boundary of the air domain.

Energy Absorption of the Soil.
Energy data in LS-DYNA are output in the d3hsp and glstat fles, and the simulation's energy change follows the energy conservation equation [43]: where E kin is the current kinetic energy, E int is the current internal energy, E si is the current sliding interface energy (including friction), E rw is the current rigid wall energy, E damp is the current damping energy, E hg is the current hourglass energy, E 0 kin is the initial kinetic energy, E 0 int is the initial internal energy, and W ext is the external work.
It is important to note that damping energy, sliding energy, and external work are not considered in this numerical simulation.Te dissipation of energy upon the shock wave's arrival at the air boundary is incorporated into the erosion energy calculation, which is output in the bndout fle.Te system starts with a total kinetic energy of 0. Te initial internal energy comprises the combined internal energies of both the air and the explosive.Tese energies are obtained from the given numerical simulation parameter E 0 , respectively.Various energy change curves are plotted as shown in Figure 17.

Shock and Vibration
Te analysis reveals that the initial internal energies of air and explosives are 1.07 × 10 7 J and 6.30 × 10 7 J, respectively.By 0.39 ms, the total energy decreases to 7.07 × 10 7 J.Both the pipe's kinetic and internal energies are minimal, with an hourglass energy of only 6.12 × 10 3 J.Importantly, the energy loss due to reduced integration in the display analysis is negligible, as it difers by four orders of magnitude from the total energy.Despite inevitable iterative errors and boundary condition infuences in the numerical simulation's iterative calculation process, the energy loss remains unavoidable.Te energy ratio in this numerical simulation stands at approximately 95.9%, ensuring energy conservation within an acceptable range.
Building upon this, energy-time curves for soil with varying thicknesses (ranging from 30 to 4800 mm) were generated, as depicted in Figure 18.Analysis reveals that the energy absorption value of the soil layer is approximately 8.0 × 10 6 J.Moreover, as the soil thickness increases, the spacing between adjacent horizontal data sections decreases.For further analysis, a soil thickness of 30 mm is established as the reference point, with the equation defned as follows: where η is the proportional multiple of soil thickness and d is the rear soil thickness.
To depict the energy absorption rate of the soil, a functional relationship between η and the absorbed energy by the soil is established.Te data points are ftted, resulting in ftted curves and relations with 96% confdence, as illustrated in Figure 19.Te ftted equation for the data takes the form of a logarithmic function with a second-order derivative less than 0, indicating a decreasing energy absorption rate of the soil.Tis suggests the existence of a suitable design value for the thickness of the rear cover soil that can satisfy both safety distance requirements and  engineering construction costs simultaneously.Further research is warranted to determine the optimal design value of the soil cover thickness.

Conclusion
In this paper, both explosion experiments and simulations were undertaken to explore the propagation of external pressure and assess the infuence of rear soil thickness on an earth-covered magazine during internal explosions.Te primary fndings are summarized as follows: (1) With an increase in the thickness of the earth covering behind the model compared to a rear soil thickness of 30 mm, signifcant reductions are noted in both air pressure surges and positive phase impulses.Specifcally, in the near zone of the explosion, air overpressure decreases by approximately 22%, accompanied by a 6% decrease in positive phase impulse.Similarly, in the far zone, there is a significant decrease of about 39% in air overpressure, coupled with a 24% decrease in positive phase impulse.Tese fndings underscore the efectiveness of augmenting rear soil thickness in mitigating blast shock wave damage.(2) Te disposal of soil from the model exhibits signifcant nonuniform characteristics, with soil accumulation primarily concentrated in a T-shaped area spanning ±90 °and 180 °.Te soil absorbs approximately 8.0 × 10 6 J of energy, representing roughly 10.85% of the cumulative internal energy within the system.Moreover, as the rear soil thickness increases, the energy absorption rate of the soil decreases.
(3) Te outcomes of the numerical simulations closely match the measured test results, afrming the credibility of the model.Tis alignment serves to afrm the precision and dependability of the model.Additionally, the ALE fuid-structure coupling algorithm proves highly efective in accurately simulating the experimental procedure for conducting internal explosion tests on the model of the earthcovered magazine.

4. 2 .
Material Properties.Te material's internal constitutive law establishes how the stress tensor correlates with the strain tensor.Additionally, following external force impact, pressure, temperature, density, internal energy, and volume changes are interrelated, as elucidated by the equation of state (EOS).For accurate analysis of explosion phenomena, matching the appropriate equations of state with the material ensures dynamic response and volume change alignment following exposure to explosion forces.

Figure 9 :
Figure 9: Te scattering and accumulation of soil.

Figure 11 :
Figure 11: Te numerical simulation model.(a) Geometry of the model.(b) Finite element mesh of pipe and explosive.(c) Finite element mesh of soil.

Figure 12 :
Figure 12: Comparison of test and numerical simulation.

Figure 14 :Figure 15 :Figure 16 Figure 16 :
Figure 14: Te process of soil modeling.(a) Set initial background mesh.(b) Fill soil in designated area.(c, d) Replace with air.(e) Modeling completed.

Table 1 :
Parameters and dimensions for internal explosion.

Table 3 :
Te distance from the measuring points of the blast center.Point number Scaled distance (m/kg 1/3 ) Actual distance (m)

Table 4 :
Shock wave-related parameters.were 5250 × 1600 × 5200 mm, while the TNT model dimensions were 450 × 100 × 200 mm. Figure 11 depicts the model diagram with a rear soil thickness of 30 mm.As the element size decreases, the relative error diminishes, albeit with increased computational time.

Table 7 :
Comparison of test and numerical simulation.