Detonation Of Explosives Research Articles

Study on shock wave load characteristics from condensed phase explosive detonations in deep-water explosions

This study investigates shock wave load characteristics from condensed phase explosive detonations in deep-water environments using a high-order compressible multiphase solver. Spatial terms of the solver are discretized by fifth-order weighted essentially non-oscillatory reconstruction in characteristic space, while a third-order total variation diminishing Runge–Kutta method is adopted to deal with the temporal terms. The level-set method captures multiphase interfaces, while a programed burn model describes detonation wave generation. Numerical and experimental validations focus on shock waves in explosives interacting with water. Validations across shallow and deep-water conditions align numerical results with theoretical and experimental values. The solver examines shock wave characteristics across varied charge masses and water depths, revealing peak pressure deviations under identical conditions. The numerical simulation results indicate that the similarity of shock wave loads in underwater explosions is evident not only in peak pressures but also in the pressure–time history curves. Through extensive comparative analysis of results, it has been found that existing formulas for calculating shock wave positive pressure durations are not applicable to deep-water explosions. The research findings and analytical methods presented in this paper can serve as crucial references for further studies on the characteristics of shock wave loads in deep-water explosions.

Development of Best Practice Guideline for the management of hot holes in surface coal mines

Surface coal operations in which mining activities are conducted above old underground workings experience hot holes after drilling. In a coal mine, hot holes are defined as shot holes, which after being drilled have an in-hole ambient temperature of 40°C or above or show a temperature increase of 3°C or more during monitoring. Hot holes and other cavities, such as cracks, pose health and safety risks to workers, such as exposure to hot air and high concentrations of noxious gases released from these holes. In addition, workers may be exposed to premature detonation of explosives caused by in-hole temperature increases and chemical reactions. To this end, Coaltech Research Association commissioned a project to develop a Best Practice Guideline for the management of hot holes. This paper is a compilation of the activities conducted in the development of the Best Practice Guideline between 2021 and 2022. The activities included a review of standard operating procedures and hot-hole temperature-measuring and -monitoring devices, and the assessment of current hot-hole procedures at selected mines. The results indicated that management of hot holes requires a focus on pre-emptive risk assessment of mining blocks, identification of hot holes using the correct temperature-measuring devices, and continuous monitoring of hot holes from the time of drilling until just before blasting. Hot-hole management accessories, such as polyvinyl chloride sleeves, were found to be effective in insulating the hot-hole emulsion from the rock mass temperature, thus preventing the potential for premature detonation.

Enhancing structural protection of reinforced concrete structures against high dynamic loads using post-installed UHPFRC/UHPFRSC

Recently there has been an observed change in security policy across the globe which implies that humanity must prepare for hybrid and asymmetric acts of war. This security situation will result in spontaneous attacks from non-military groups. The threats include ballistic attacks from commercially available weapons, detonation of explosives, and impact of vehicles. These attacks are focused on critical infrastructures. Reinforced concrete (RC) elements with ballistic and explosion protection are of huge relevance. This requires new concepts and innovative solutions for the protection of new buildings and/or for the strengthening of existing structures. In this paper a method for the increase of structural safety is shown using Ultra-High Performance Fiber Reinforced Concrete (UHPFRC) or sprayed UHPFRC, so called UHPFRSC, as a surface layer for reinforced concrete structures. The UHPFRC/UHPFRSC is applied in thin layers of 30 - 80 mm on the protection side. The interaction of the two different cementitious materials results in increased resistance to high dynamic loads. The traditional reinforced concrete on the impact side absorbs a major fraction of the energy of impact. On the back side (protection side), tensile forces are generated by the impact, which can excellently be handled by the UHPFRC/UHPFRSC. For proof of concept of the hybrid strengthening method, tests were carried out on UHPFRC - RC composite elements. In the test campaign, ballistic tests were performed using test specimens with dimensions of 500 × 500 mm. The layer thicknesses were 120 mm reinforced concrete and 40 mm UHPFRC. Additional steel rebars in the UHPFRC layer were applied for half of the specimens. The investigation demonstrated the effectiveness of the strengthening method, showing high resistance to high-dynamic loads. This behavior is attributed to the absorption and transfer of tensile stresses by the UHPFRC, that reduces penetration depth and damage, ensuring remaining load-bearing capacity after gunfire.

A novel method for investigating the underwater explosion loads and bubble evolution

This paper presents an innovative experimental method for studying the evolution and energy output characteristics of underwater explosion bubbles. We independently constructed an experimental testing system for underwater electrical wire explosions (UEWE), in which electrodes connected to a metal wire serve as the load, and underwater explosions are initiated through instantaneous high-voltage discharge. By varying the diameter of the metal wire and configuring parallel wire arrays, we analyzed and discussed the explosion characteristic parameters and the current–voltage (I–V) signals under different conditions. The maximum bubble radius of the underwater metal wire explosion was compared with the corresponding equivalent explosive simulation results, and a numerical model for underwater metal wire explosion equivalent to explosive detonation was established. Subsequently, we discussed the characteristics of bubble generation and evolution under various conditions, clarifying the similarities and differences between wire explosions and explosive detonations. On this basis, we explored the propagation laws of shock waves and secondary pulsation waves (SPW) under different conditions. We also calculated and analyzed energy output characteristic parameters, such as shock wave energy and bubble energy. The results indicate that there are significant differences between copper wire and aluminum wire loads in UEWE. For copper wires with a diameter greater than 0.4 mm, the shock wave overpressure peak value significantly decreases, while for aluminum wires with a diameter greater than 0.5 mm, it slightly decreases. Both metals exhibit similar trends in parallel wire arrays, with the shock wave overpressure peak value initially increasing and then decreasing as the number of wires increases. Unlike underwater explosive detonations, the SPW peak value in UEWE may exceed that of the shock wave. For single wires, the SPW peak value of copper wires is generally higher than that of aluminum wires, but in wire arrays, the trend is reversed. The multi-wire parallel connection can improve the energy conversion efficiency of the shock waves. However, for bubble energy, under all conditions, a single aluminum wire with a diameter of 0.5 mm produced the maximum bubble energy, reaching 1023.1 J. These findings provide new insights into the energy features of UEWE.

Impact of prototype scale effect on the dynamic responses and damage of steel plates under explosions

Scaling model was widely used in engineering to reveal the responses of full-scale prototype model. This study provided insight into the impact of the scale effect on the dynamic response and damage of stiffened and unstiffened steel plates under explosions. A numerical model consider fluid-structure-interaction (FSI) was established using LS-DYNA platform, and the numerical model was verified by the blast test results. Moreover, the scaling factors (S) equal to 1 (prototype), 1/2, and 1/4, and 1/8, and 1/12 were considered in the study, the scaling effects on the breach size, plastic area, and the deformation shape were analyzed. The results showed that the breach size and plastic area of the scaled model deviated from the ideal results based on the replica law. However, the deformation shapes of scaled-down plates accurately reflected the deformation characteristics of the full prototype plate. The scaling factor significantly influenced the dynamic responses of stiffened and unstiffened steel plates; as the scaling factor (S) increased, the dimensionless displacement decreased. The results indicated that a smaller S led to a larger deviation in response. The influencing factors of the scaling effect were investigated, and it was found that stiffened steel plates were more sensitive to the scaling factor than unstiffened steel plates. By comparing the detonation of cuboid, cylindrical, and spherical explosives blast-loaded onto stiffened steel plates, it was determined that the detonation of cubic explosives caused the most severe scaling effect on stiffened steel plates under blast loads. The study also found that explosive mass and standoff distance influenced the scaling effect of stiffened steel plates under blast loads. Additionally, steel plates in the plastic stage and the critical state between the elastic and plastic stages were particularly sensitive to the scaling effect. Due to the strain rate effect, which can lead to a deviation from the ideal scale law, different steel materials exhibit varying scale sensitivities. Comparing stiffened plates made from Q235, Q345, DH36, and L907A, the order of sensitivity to the scaling factor was: Q235 > Q345 ≈ DH36 > L907A.

Application of Explosive Equivalency Approach in Blast-Induced Seismic Effect Prediction Using EXPLO5 Thermochemical Code

Blasting is a key process that plays a significant role in various industries, including mining and construction. To measure the effectiveness and potential impact of a blast generated by different explosives, industry professionals use a widely accepted parameter known as TNT (trinitrotoluene) equivalent. This manuscript provides an overview of the approach based on the application of the explosive equivalency principle in the prediction of the seismic effects caused by the detonation of different explosives. The explosive equivalents of studied explosives are derived from the results of thermochemical calculations using the EXPLO5 code and compared to field tests. The results have demonstrated that the equivalency approach can potentially be a useful tool in the assessment of blast-induced seismic effects.

Shock-induced nanoscale pore collapse and hotspot in cyclotetramethylene tetranitramine (HMX)

Pore collapse induced hotspot formation is a key mechanism for initiating detonation of shocked high explosives. Accurate continuum models that faithfully capture the pore collapse dynamics and resulting hotspot temperatures for multiscale initiation of high explosive are still lacking. Here, an atomistically informed dislocation plasticity model is developed for cyclotetramethylene tetranitramine (HMX), which includes nonlinear thermoelasticity, pressure-dependent and temperature-dependent dislocation motion and homogeneous dislocation nucleation, and new melting criteria. Nanoscale pore collapse simulations based on the proposed model could well reproduce the pore collapse transition from a strength-dominated regime to a hydrodynamic regime observed in molecular dynamics results. During the pore collapse transition, dislocation motion and dislocation generation induced plasticity is decoupled and evaluated. It indicates that the homogeneous dislocation nucleation suppresses the dislocation motion. Through linkage analysis of dislocation physical quantity histories of tracer particles, this work clarifies the interplay between physical dependences and dislocation mechanisms on pore collapse behaviors for the first time. The results reveal that the transition of pore collapse regime emerges from the strong pressure-dependent homogeneous dislocation nucleation. In addition, the pore collapse responses and energy localizations of the proposed model are further compared with different strength models to prove the importance of pressure dependence and homogeneous dislocation nucleation on pore collapse behaviors. Based on the simulations of pore collapse, the potential heat sources are quantified to provide a physical understanding of hotspot formation mechanisms. The present work demonstrates a key step in multiscale modeling of high explosive initiation in that atomistic simulation results are directly used to guide and refine the mesoscale model of energetic crystal used in the continuum.

Thermal Safety Study of Emulsion Explosive Matrix under the Coupled Effects of Environmental Pressure and Bubble Content with Internal Heat Source

Emulsion explosives have become a hot topic in various studies due to their explosive combustion characteristics and detonation performance under different environmental pressures. The thermal safety of an emulsified matrix was studied with ignition energy as the characterization. A minimum ignition energy test experimental system for emulsion matrices was established in this research. The system simulated the occurrence of hot spots inside emulsion matrices using an electric heating wire. The effect of bubbles on the thermal safety of the emulsified matrix was studied by adding expanded perlite additive to the emulsified matrix. This study investigated the variation trend in the minimum ignition energy of the emulsion matrix under the coupled effect of bubbles and ambient pressure using the orthogonal experimental method. The impacts of two factors on the thermal safety of the emulsion matrix were studied at different hot-spot temperatures. Coupled analysis experiments were conducted on emulsion matrices containing 0%, 1.5%, and 3% expanded perlite under pressure environments of 1 atm, 2 atm, and 3 atm. The critical hot-spot temperature of the emulsion matrix significantly decreases with increasing bubble content at 1 atm and 2 atm pressures, as revealed by intuitive analysis and analysis of variance. However, at 3 atm of pressure, the bubble content in the emulsion matrix has no significant effect on its critical hot-spot temperature. The results show that the thermal safety of the emulsified matrix decreases with the increase in the content of expanded perlite and environmental pressure, and the influence of environmental pressure is more significant than that of the bubble content. This paper’s research content serves as a reference for a safe emulsified matrix and as an experimental basis for establishing a production line for developing new equipment.

Effect of explosive detonation velocity on the degree of rock pre-fracturing during blasting

Effect of explosive detonation velocity on the degree of rock pre-fracturing during blasting

Probabilistic analysis of near-field blast loads considering fireball surface instabilities and stochastic detonator location

High speed video analysis of near-field explosive detonations displays distinct stages of emergent hydrodynamic instabilities in the fireball/shock-air interface. Typically, beyond 10 charge radii, the instabilities experienced large growths giving rise to more chaotic behaviour of the interface and thus an increasing uncertainty in surface velocity. These surface instabilities are suggested as the primary cause of blast parameter variability in the near-field. However, as a deterministic tool, numerical simulation of the detonation process and subsequent blast wave propagation is not able to replicate the stochastic nature of fireball surface instabilities and hence near-field blast parameter variability. Therefore, it is necessary to develop new methods to simulate and characterise the stochastic features of the fireball/shock-air interface. This paper proposes an algorithm to generate an explosive charge element with random shape in finite element model in order to simulate irregularities in the fireball/shock-air interface, and therefore produce variabilities comparable to those from direct observation. The effect of chaotic fireball/shock-air interface on near-field loading is explored through a large number of numerical simulations in order to investigate the statistical distribution of parameters including peak overpressure and impulse. Subsequently, the effect of stochastic detonator location is explored in a similar manner. A computational procedure based on the Monte Carlo Method is proposed to establish a probabilistic model of near-field blast loads, termed PSL-Blast. The reliability of design blast loads calculated using the UFC 3-340-02 design manual is then estimated using PSL-Blast, which suggests that reliability decreases with decreasing scaled distance. Finally, reliability-based safety factors of blast loads are calculated based on different blast settings.

Design, preparation and firing characterization of PCB-EBW

Exploding Bridge Wire (EBW) is a special actuator that drives high-temperature and high-pressure plasma. To address the issues of low efficiency and high cost in the preparation process of EBW detonators, an EBW chip based on Printed Circuit Board (PCB) technology was designed and studied. By utilizing a specially designed "S"-shaped transition region and using copper (Cu) as the bridge wire material, the PCB method has enabled batch production of the PCB-EBW device. Four different types of PCB-EBW chips were constructed by designing bridge wires of varying lengths. Subsequently, electro-thermal simulations, electrical explosion performance studies, and shock wave simulations were conducted on the PCB-EBW. It is inferred from the research findings that the main factor triggering the explosive detonation of the PCB-EBW chip is the coupling effect between the high-temperature plasma generated by the electrical explosion and the electrical explosion shock wave. The speculation on the working mechanism of the PCB-EBW was validated by igniting PETN, with the minimum firing voltage determined to be 550 V (for a bridge wire length of 200 μm). This novel EBW design expands the application prospects and provides new insights into the mechanism of EBW for future research.

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

Detonation Soot: A New Class of Ice Nucleating Particle

AbstractTemperatures and pressures from high explosive detonations far exceed atmospheric conditions in typical combustion reactions, and consequently, detonation soot forms with physiochemical properties distinct from soot formed by combustion. In this study, samples of detonation soot from two high explosives, PBX 9502 and Composition B‐3, were analyzed. Ice nucleation experiments on soot collected after controlled detonations were conducted in the laboratory to probe immersion and contact mode freezing. Samples nucleated ice at temperatures warmer than commercially available nanodiamonds, which has a mean nucleation temperature of −20.7°C. Ice nucleation rate coefficients increase rapidly by two to three orders of magnitude below −20°C for every sample. Size‐selected 137 μm diameter particles produced during detonation in an ambient air atmosphere yield bimodal distributions of freezing temperature with primary and secondary nucleation modes centered at −20°C and −13°C, respectively. The presence of a secondary mode allows for enhanced ice nucleation rate coefficients (one to two orders of magnitude greater than samples without a secondary mode) at temperatures outside the influence of the primary mode (>−17°C). Given the observed onset nucleation temperatures of −9.2°C, our results imply that detonation soot of the type studied here would only need to reach an altitude of approximately 4 km to facilitate ice formation.

Simulation-based personal fatality risk assessment due to the fragmentation hazard

In the military and chemical industry, modeling the fragmentation hazard field is of great significance in conducting the fatality risk assessment and calculating the safety distance for person. Although the ballistic methodology achieves the simulation of fragmentation flight trajectory, the acquisition of accurate initial projection data of fragmentation is a challenge. By coupling continuum-discontinuum element method and particle discrete element method, the fragmentation power algorithm is established, which achieves the integrated simulation of the initial projection data of fragmentation and the flight process of fragmentation. First, continuum-discontinuum element method is adopted to simulate the detonation products-driven fragmentation acceleration process by introducing the explosive detonation model, which achieves the acquisition of initial projection data of fragmentation. Then, the particle discrete element method is adopted to simulate the flight trajectory and power data of fragmentation. Based on the numerical and experimental result, the accuracy of fragmentation power algorithm is verified. In conjunction with the personal vulnerability model, a systematic numerical simulation framework is established to conduct the simulation-based personal fatality risk assessment when the metal-cased munition detonates accidentally, and the results indicate that the personal fatality risk due to the fragmentation becomes severer with the increase of munition curvature.

Experimental studies of confined detonations of plasticized high explosives in inert and reactive atmospheres

When explosives detonate in a confined space, repeated boundary reflections result in complex shock interactions and the formation of a uniform quasi-static pressure (QSP). For fuel-rich explosives, mixing of partially oxidized detonation products with an oxygen-rich atmosphere results in a further energy release through rapid secondary combustion or ‘afterburn’. While empirical formulae and thermochemical modelling approaches have been developed to predict QSP, a lack of high-fidelity experimental data means questions remain around the deterministic quality of confined explosions, and the magnitude and mechanisms of afterburn reactions. This article presents experimental data for RDX- and PETN-based plastic explosives, demonstrating the high repeatability of the QSP generated in a sealed chamber using pressure transducers and high-speed infrared thermometry. Detonations in air, nitrogen and argon atmospheres are used to identify the contribution of afterburn to total QSP, to estimate the duration of afterburn reactions and to speculate on the flame temperature associated with this mechanism. Computational fluid dynamic modelling of the experiments was also able to accurately predict these effects. Understanding and quantifying explosions in complex environments are critical for the design of effective protective structures: the mechanisms described here provide a significant step towards the development of fast-running engineering models for internal blast events.

Blast wave interaction during explosive detonation in a variable cross-sectional charge

Purpose. The research aims to assess the efficiency and performance of solid media destruction in directed blasting of a charge with variable cross-sectional shape. Methods. Numerical modelling of the blast wave interaction process is performed using the finite element method based on the Euler-Lagrange algorithm. The Johns-Wilkins-Lee equation of state is used to determine the pressure-volume dependences of medium destruction. Assessment of solid medium destruction mechanism during a directed blasting of a charge with variable cross-sectional shape is carried out based on polarization-optical method on models made of optically active material. Findings. Experimental studies of solid medium destruction by the action of directed blasting with a variable cross-section charge made it possible to determine the direction of blast wave propagation and its amplitude in stress wave, influencing the intensity of radial crack network formation in superposition areas, and directed perpendicularly to explosive cavity. At the same time, the average peak pressure in collision zone of two shock waves in centre of spherical cavity is approximately 1.48 and 1.84 times higher than that in weakly blast-loaded areas. It has been found that when two shock waves collide and superimpose on each other, the intensity of their impact increases. Moreover, the shock wave velocity in collision zone is higher than that of the radial shock wave. Originality. It has been determined that the maximum pressure values on the explosive cavity wall at the initiation points sharply increase and then gradually stabilize as the blast stress waves propagate and have an arbitrary distribution pattern. Three areas should be considered: not superimposed, weakly superimposed, and strongly superimposed. At each point of detonation, the pressure on explosive cavity wall will be minimal, while in the charge centre in the spherical insert zone, on the contrary, it will be maximal. In this case, the pressure in the central superposition area is about 2.84 times greater than at the initiation ends, and the nature of distribution changes according to a linear dependence. Practical implications. The performed research findings can serve as a basis for development of effective parameters of resource-saving methods for stripping hard rocks of complex structure in the conditions of ore mines.

Investigation of thermal damage in explosive bridgewire detonators via discrete element method simulations

AbstractExploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free‐standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. We present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

Reflection behavior of insensitive explosive detonation propagating around a cylinder

Reflection behavior of insensitive explosive detonation propagating around a cylinder

Thermal Hazard Analysis of Two Non-Ideal Explosives Based on Ammonium Perchlorate/Ammonium Nitrate and Aluminium Powder.

In recent years, various kinds of civil explosive detonation accidents have occurred frequently around the world, resulting in substantial human casualties and significant property losses. It is generally believed that thermal stimulation plays a critical role in triggering the detonation of explosives; consequently, the study of the thermal hazards of explosives is of great significance to many aspects of safety emergency management practices in the production, transportation, storage, and use of explosives. It is known that the thermal stability of the ammonium perchlorate-aluminium system and the ammonium nitrate-aluminium system has been extensively investigated previously in the literature. However, there is a paucity of research on the thermal hazard characteristics of non-ideal explosives under varying oxygen balance conditions within the academic sphere. Therefore, this research focused on the study of the thermal hazards of non-ideal explosives based on thermokinetic analysis. The thermal hazards of non-ideal explosive mixtures of ammonium perchlorate and aluminium and of ammonium nitrate and aluminium were studied by thermal analysis kinetics. The thermokinetic parameters were meticulously studied through differential scanning calorimetry (DSC) analysis. The results showed that the peak reaction temperature and activation energy of the ammonium perchlorate-aluminium system were significantly higher than those of the ammonium nitrate-aluminium system. Under the condition of zero oxygen balance, the peak reaction temperature of the ammonium nitrate-aluminium system was 259 °C (heating rate 5 °C/min), and the activation energy was 84.7 kJ/mol. Under the same conditions, the peak reaction temperature and activation energy of the ammonium perchlorate-aluminium system were 292 °C (heating rate 5 °C/min) and 94.9 kJ/mol, respectively. These results indicate that the ammonium perchlorate-aluminium system has higher safety under the same thermal stimulation conditions. Furthermore, research on both non-ideal explosive systems reveals that the activation energy is at its peak under negative oxygen balance conditions, recorded at 104.2 kJ/mol (ammonium perchlorate-aluminium) and 86.2 kJ/mol (ammonium nitrate-aluminium), which indicates a higher degree of safety. Therefore, the investigation into the thermal hazards of non-ideal explosive systems under different oxygen balance conditions is of utmost importance for the enhancement and improvement of safety emergency management practices.

Ensuring the protection of critical infrastructure energy facilities from explosions

Purpose. Analysis of the protection of critical energy infrastructure facilities, which are often destroyed under the influence of a shock wave, which affects the stability of Ukraine's energy systems. Determination of the state of security of critical infrastructure, which ensures the functionality, continuity of work, integrity and stability, which is reflected in the life support of the population. Methodology. To achieve this goal, a simulation of a detonation explosion and the distribution of shock waves within prismatic bodies simulating urban development using the ANSYS AUTODYN software product were carried out. In conditions of dense urban development, measures to ensure the protection of critical infrastructure facilities are selected based on the interaction of liquids, gases and solids, phase transitions, propagation of shock waves, etc. The results of the study were obtained. The calculated parameters of the shock wave distribution limit in 15 ms after detonation were determined, the pressure data corresponding to the measurement sensors T1 and T2 were recorded, it was determined that the values of the pressure of the blast wave at the manometric points T1 and T2 are the highest about 5 ms, аt the same time, the pressure on the underground structure is almost 50% less than on the ground. Protection of critical energy infrastructure facilities in the form of construction of underground transformer substations and their modernization in conditions of dense urban development is proposed. This is the optimal solution to eliminate problems, which arise with an increase in loads in city power grids and protect them from the action of a blast wave PCT. Originality. The dependencies of the pressure distribution on the shock wave in time for two variants of shading of critical energy infrastructure objects, namely underground location and ground location in the urban environment, are obtained. At the same time, the parameters of the materials of the structure and the characteristics and amount of explosives were taken into account. Practical value. The construction of underground transformer substations and their modernization in dense urban areas is the optimal solution for eliminating the problems that arise with an increase in loads in city power grids and protecting them from the effects of a blast wave.