Online First
Articles in press have been peer-reviewed and accepted, which are not yet assigned to volumes/issues, but are citable by Digital Object Identifier (DOI).
Display Method:
, Available online , doi: 10.11883/bzycj-2026-0057
Read OL
PDF 
Cited By
Abstract:
To address the issues of insufficient cycle advance and low blast-hole utilization resulting from inappropriate cut blasting layouts during roadway excavation, the empty-hole effect was employed to establish a three-dimensional cut blasting model using FEM-SPH（finite element method-smooth particle hydrodynamics）coupled numerical simulation. In the numerical modeling, the rock mass was described by the RHT constitutive model, which accounts for strain-rate sensitivity and confinement effects under high-strain-rate loading. The explosive was modeled using the Jones-Wilkins-Lee equation of state. A central charged hole was surrounded by four, five, or six empty holes with a spacing of 210 mm determined theoretically. SPH particles were applied to the central zone to capture large deformation, fragmentation, and rock ejection, while the FEM was used for the outer region to reduce computational cost. Non-reflecting boundaries were set to avoid stress wave reflection. The model was then used to simulate and analyze the influence of different cut blasting arrangements on rock mass damage, rock throw efficiency, and fragmentation patterns. Furthermore, field engineering tests were conducted in a copper mine roadway to validate the simulation results by examining the residual hole depth after blasting and the achieved cycle advance. The research findings indicate that a cut blasting layout with one central charged hole and five empty holes provides a moderate number of free surfaces, allowing for more sufficient superposition of stress waves and resulting in better blasting performance. Compared to layouts with four or six empty holes, the five-empty-hole configuration demonstrates optimal rock throw effects, with the slot cavity opening size increasing by 42.5% and 20.3%, respectively, and the slot cavity cross-sectional size improving by 52.4% and 34.7%, respectively. Field test results show that the cycle advance achieved with the five-empty-hole layout is 10.2% and 3.2% higher than that of the four-empty-hole and six-empty-hole layouts, respectively. The research reveals a nonlinear relationship between the number of empty holes and blasting effects, demonstrating that the optimal balance between stress wave superposition and energy concentration can be achieved using a moderate number of empty holes, providing a basis for the design and optimization of cut blasting schemes in roadway excavation.
To address the issues of insufficient cycle advance and low blast-hole utilization resulting from inappropriate cut blasting layouts during roadway excavation, the empty-hole effect was employed to establish a three-dimensional cut blasting model using FEM-SPH（finite element method-smooth particle hydrodynamics）coupled numerical simulation. In the numerical modeling, the rock mass was described by the RHT constitutive model, which accounts for strain-rate sensitivity and confinement effects under high-strain-rate loading. The explosive was modeled using the Jones-Wilkins-Lee equation of state. A central charged hole was surrounded by four, five, or six empty holes with a spacing of 210 mm determined theoretically. SPH particles were applied to the central zone to capture large deformation, fragmentation, and rock ejection, while the FEM was used for the outer region to reduce computational cost. Non-reflecting boundaries were set to avoid stress wave reflection. The model was then used to simulate and analyze the influence of different cut blasting arrangements on rock mass damage, rock throw efficiency, and fragmentation patterns. Furthermore, field engineering tests were conducted in a copper mine roadway to validate the simulation results by examining the residual hole depth after blasting and the achieved cycle advance. The research findings indicate that a cut blasting layout with one central charged hole and five empty holes provides a moderate number of free surfaces, allowing for more sufficient superposition of stress waves and resulting in better blasting performance. Compared to layouts with four or six empty holes, the five-empty-hole configuration demonstrates optimal rock throw effects, with the slot cavity opening size increasing by 42.5% and 20.3%, respectively, and the slot cavity cross-sectional size improving by 52.4% and 34.7%, respectively. Field test results show that the cycle advance achieved with the five-empty-hole layout is 10.2% and 3.2% higher than that of the four-empty-hole and six-empty-hole layouts, respectively. The research reveals a nonlinear relationship between the number of empty holes and blasting effects, demonstrating that the optimal balance between stress wave superposition and energy concentration can be achieved using a moderate number of empty holes, providing a basis for the design and optimization of cut blasting schemes in roadway excavation.
, Available online , doi: 10.11883/bzycj-2025-0177
Abstract:
The detonation-driven expansion and fracture behavior of single-layered and double-layered metallic cylindrical shells (with equal total thickness) were investigated via the DPS (Doppler pins system) array and high-speed photography. The velocity curves and high-speed photography images of the outer surface of the single- and double-layered cylindrical shells were obtained. The velocity curves show that spallation occurs at approximately 40% of the thickness from the outer surface in the single-layered cylindrical shell. In addition, the spallation part was caught up after ~1.89 μs by the fragments of the single-layered cylindrical shell, leading to the secondary impact loading. In the case of double-layered cylindrical shell, the inner and outer shells of the double-layered cylindrical shell separated quickly since they cannot withstand the tensile stress. The outer cylindrical shell expands freely with high velocity at initial stage. It is then caught up by the inner cylindrical shell after 26.13μs, leading to the secondary impact loading. There is a significant difference in secondary impact loading times the single- and double-layered shell expansion experiments. The premature separation of the double-layered cylindrical shell hinders the crack penetration from the inside to the outside, and delays the crack penetration until the complete fracture of the whole cylindrical shell. In contrast to the widely investigated fracture behavior in a single-layered cylindrical shell, the detonation product leakage time is delayed when the outer layer of the double-layered cylindrical shell expands. The fragment size of a double-layered shell is obviously smaller than that of a single-layered shell. The distribution of the overflowing detonation product is relatively scattered and rare, and the overall crack of the outer cylindrical shell is clearly visible. The results of SPH (smoothed particle hydrodynamics) numerical simulations demonstrate that under certain conditions, the lower the spallation strength in single-layer cylindrical shells, the later the residual fragments catch up with the spalled layer during secondary loading, the less likely the occurrence of penetration fracture. Meanwhile, the interfacial gaps in double-layered cylindrical shells significantly inhibit the penetration and propagation of cracks.
The detonation-driven expansion and fracture behavior of single-layered and double-layered metallic cylindrical shells (with equal total thickness) were investigated via the DPS (Doppler pins system) array and high-speed photography. The velocity curves and high-speed photography images of the outer surface of the single- and double-layered cylindrical shells were obtained. The velocity curves show that spallation occurs at approximately 40% of the thickness from the outer surface in the single-layered cylindrical shell. In addition, the spallation part was caught up after ~1.89 μs by the fragments of the single-layered cylindrical shell, leading to the secondary impact loading. In the case of double-layered cylindrical shell, the inner and outer shells of the double-layered cylindrical shell separated quickly since they cannot withstand the tensile stress. The outer cylindrical shell expands freely with high velocity at initial stage. It is then caught up by the inner cylindrical shell after 26.13μs, leading to the secondary impact loading. There is a significant difference in secondary impact loading times the single- and double-layered shell expansion experiments. The premature separation of the double-layered cylindrical shell hinders the crack penetration from the inside to the outside, and delays the crack penetration until the complete fracture of the whole cylindrical shell. In contrast to the widely investigated fracture behavior in a single-layered cylindrical shell, the detonation product leakage time is delayed when the outer layer of the double-layered cylindrical shell expands. The fragment size of a double-layered shell is obviously smaller than that of a single-layered shell. The distribution of the overflowing detonation product is relatively scattered and rare, and the overall crack of the outer cylindrical shell is clearly visible. The results of SPH (smoothed particle hydrodynamics) numerical simulations demonstrate that under certain conditions, the lower the spallation strength in single-layer cylindrical shells, the later the residual fragments catch up with the spalled layer during secondary loading, the less likely the occurrence of penetration fracture. Meanwhile, the interfacial gaps in double-layered cylindrical shells significantly inhibit the penetration and propagation of cracks.
, Available online , doi: 10.11883/bzycj-2026-0002
Abstract:
To further analyze the damage effects of projectile penetration into concrete targets and their influence on the penetration process, a systematic theoretical investigation was carried out. A spalling effect model was established based on the spalling phenomenon observed during normal perforation of concrete targets by projectiles. By integrating this model with the cratering effect model and the target resistance function, a theoretical model for penetration that comprehensively accounts for both cratering and spalling effects was developed. By means of the presented model, the influence of target thickness and impact velocity on the penetration process and target damage parameters were analyzed. The results show that as the concrete target thickness decreases or the impact velocity increases, the influence of the spalling effect on penetration resistance and residual velocity diminishes, while the proportion of the splash zone during the spalling stage increases, making the splash effect the dominant factor affecting resistance in the spalling stage. When the target thickness exceeds 4.5 times the projectile diameter, the influence of the spalling effect on the residual velocity exceeds 10%, rendering the spalling effect non-negligible. Cratering depth increases linearly with thickness when the thickness is not more than 4 times the projectile diameter, and then tends to stabilize. Spalling depth increases linearly with thickness when the thickness is not more than 6.9 times the projectile diameter, and then decreases linearly. During the linear increase stage, both cratering depth and spalling depth are approximately half of the concrete target thickness. Compared with perforation of thin concrete targets, variations in impact velocity exhibit a more pronounced influence on the crater depth and spalling depth during the penetration of thick concrete targets. The proposed model effectively characterizes the damage effects of concrete targets and can provide a theoretical basis for penetration analysis and protective design of concrete targets.
To further analyze the damage effects of projectile penetration into concrete targets and their influence on the penetration process, a systematic theoretical investigation was carried out. A spalling effect model was established based on the spalling phenomenon observed during normal perforation of concrete targets by projectiles. By integrating this model with the cratering effect model and the target resistance function, a theoretical model for penetration that comprehensively accounts for both cratering and spalling effects was developed. By means of the presented model, the influence of target thickness and impact velocity on the penetration process and target damage parameters were analyzed. The results show that as the concrete target thickness decreases or the impact velocity increases, the influence of the spalling effect on penetration resistance and residual velocity diminishes, while the proportion of the splash zone during the spalling stage increases, making the splash effect the dominant factor affecting resistance in the spalling stage. When the target thickness exceeds 4.5 times the projectile diameter, the influence of the spalling effect on the residual velocity exceeds 10%, rendering the spalling effect non-negligible. Cratering depth increases linearly with thickness when the thickness is not more than 4 times the projectile diameter, and then tends to stabilize. Spalling depth increases linearly with thickness when the thickness is not more than 6.9 times the projectile diameter, and then decreases linearly. During the linear increase stage, both cratering depth and spalling depth are approximately half of the concrete target thickness. Compared with perforation of thin concrete targets, variations in impact velocity exhibit a more pronounced influence on the crater depth and spalling depth during the penetration of thick concrete targets. The proposed model effectively characterizes the damage effects of concrete targets and can provide a theoretical basis for penetration analysis and protective design of concrete targets.
, Available online , doi: 10.11883/bzycj-2025-0052
Abstract:
To effectively prevent explosion hazards during propylene production, storage, and utilization, the explosion limits of propylene in air under varying initial temperatures (20 ℃–180 ℃) and initial pressures (0.1 MPa–0.9 MPa) were measured using 12L explosion limit tester. The study revealed that as the initial temperature and pressure increase, the upper explosion limit (UEL) of propylene rises significantly, while the lower explosion limit (LEL) decreases slightly, resulting in a marked broadening of the explosion limit range. At an initial temperature of 180 ℃, with the pressure increases, the carbon powder content in the explosive products significantly increases during the UEL test, and the LEL decline transitions from a linear to a sliding-curve pattern. Analysis using CHETAH 11.0 software revealed that the increase in carbon powder content is closely related to thermodynamic properties under high-temperature and high-pressure conditions. When the pressure rises from 0.1 MPa to 0.9 MPa, the carbon powder content in the explosion products surges from 3.82% to 25.88%. This significant growth primarily stems from two factors: first, high-pressure conditions promote the Boudouard reaction toward carbon formation; second, the increased total amount of reactants under high pressure leads to a multiplicative rise in product quantities.Using CHEMKIN software, the combustion characteristics of propylene under lower explosion limit conditions were investigated. In the fuel-lean region, propylene undergoes free radical chain reactions, ultimately generating CO2, which maintains the calculated adiabatic flame temperature (CAFT) above1400 K. As pressure increases, CAFT gradually decreases. Meanwhile, with rising temperature, CAFT exhibits a pressure-dependent transition: at low pressures (<0.5 MPa), CAFT increases, whereas at high pressures, it decreases, with 0.5 MPa serving as the critical transition threshold. The coupling effects of initial temperature and pressure on explosion limits are far more pronounced than those of individual factors, with a stronger impact on the UEL than the LEL. The coupling effects caused a 108% increase in the UEL and an 18.05% decrease in the LEL. For individual factors, initial temperature alone led to a 3.8% UEL increase and a 3.41% LEL decrease, while initial pressure alone resulted in a 51.3% UEL increase and 2.44% LEL reduction.
To effectively prevent explosion hazards during propylene production, storage, and utilization, the explosion limits of propylene in air under varying initial temperatures (20 ℃–180 ℃) and initial pressures (0.1 MPa–0.9 MPa) were measured using 12L explosion limit tester. The study revealed that as the initial temperature and pressure increase, the upper explosion limit (UEL) of propylene rises significantly, while the lower explosion limit (LEL) decreases slightly, resulting in a marked broadening of the explosion limit range. At an initial temperature of 180 ℃, with the pressure increases, the carbon powder content in the explosive products significantly increases during the UEL test, and the LEL decline transitions from a linear to a sliding-curve pattern. Analysis using CHETAH 11.0 software revealed that the increase in carbon powder content is closely related to thermodynamic properties under high-temperature and high-pressure conditions. When the pressure rises from 0.1 MPa to 0.9 MPa, the carbon powder content in the explosion products surges from 3.82% to 25.88%. This significant growth primarily stems from two factors: first, high-pressure conditions promote the Boudouard reaction toward carbon formation; second, the increased total amount of reactants under high pressure leads to a multiplicative rise in product quantities.Using CHEMKIN software, the combustion characteristics of propylene under lower explosion limit conditions were investigated. In the fuel-lean region, propylene undergoes free radical chain reactions, ultimately generating CO2, which maintains the calculated adiabatic flame temperature (CAFT) above
, Available online , doi: 10.11883/bzycj-2025-0351
Abstract:
Ti/steel clad plates are considered as ideal candidate materials for applications such as marine engineering and low-temperature pressure vessels operating at −70 ℃. However, their large-scale production and application are severely restricted by the manufacturing process and the quality of interfacial bonding. To investigate the influence mechanism of vacuum pressure on interfacial morphology and mechanical properties, TC4/09MnNiDR(Ti/steel) clad plates were fabricated using explosive welding under ambient pressures of 20 , 60, and 100 kPa, and the fabricated clad plates were named as CP-20, CP-60 and CP-100, respectively. A specially designed vacuum chamber was employed to precisely control the environmental pressure during the explosive process. After welding, all specimens were subjected to a uniform stress-relief annealing treatment at 550 ℃ for 2 h under a vacuum atmosphere to eliminate residual stresses and improve the reliability of subsequent characterization. The interfacial microstructure and chemical characteristics were systematically analyzed using multiple characterization techniques. Scanning electron microscopy (SEM) was used to observe interfacial morphology, including wave formation and defect distribution. Energy dispersive spectroscopy (EDS) was applied to examine elemental distribution across the bonding interface. Electron backscatter diffraction (EBSD) was employed to evaluate grain structure, grain boundary characteristics, recrystallization fraction, and texture evolution near the interface. Electron probe microanalysis (EPMA) was conducted to quantitatively determine the chemical composition of the interfacial melting zone and vortex regions, with both mapping and point analyses performed to identify the constituent phases. Mechanical properties at −70 ℃ were evaluated through tensile tests, Charpy impact tests, and three-point bending tests, all carried out in accordance with relevant national standards, with the results reported as mean values ± standard deviations. Results show that the vacuum environment significantly improves interfacial quality. With increasing vacuum degree, the interfacial wave becomes finer, more continuous, and more uniform, while the thickness of the molten layer decreases, and the presence of defects and brittle intermetallic compounds is reduced. EBSD analysis reveals grain refinement and an increased recrystallization fraction at the interface under vacuum conditions, demonstrating that the interfacial microstructure can be effectively controlled by adjusting the ambient pressure. EPMA further reveals that the vortex regions are mainly composed of Fe and Ti elements, while the weld seam region exhibits a stable chemical composition with Ti-to-Fe atom number ratios close to 1∶1 or 2∶1, indicating that the dominant intermetallic phases formed at the interface are TiFe and TiFe2. Benefiting from the optimized interfacial structure, the clad plates exhibit excellent mechanical properties at −70 ℃. The tensile strengths of the CP-20, CP-60 and CP-100 are 880, 911, and 867 MPa, respectively. The impact absorbed energies are 17.5, 10.2, and 6.4 J, respectively. And the flexural strengths are1469 , 1350 , and 1167 MPa, respectively. This study demonstrates that vacuum explosive welding is a reliable technique for producing high-performance low-temperature metal composites, with 60 kPa identified as the optimal processing window.
Ti/steel clad plates are considered as ideal candidate materials for applications such as marine engineering and low-temperature pressure vessels operating at −70 ℃. However, their large-scale production and application are severely restricted by the manufacturing process and the quality of interfacial bonding. To investigate the influence mechanism of vacuum pressure on interfacial morphology and mechanical properties, TC4/09MnNiDR(Ti/steel) clad plates were fabricated using explosive welding under ambient pressures of 20 , 60, and 100 kPa, and the fabricated clad plates were named as CP-20, CP-60 and CP-100, respectively. A specially designed vacuum chamber was employed to precisely control the environmental pressure during the explosive process. After welding, all specimens were subjected to a uniform stress-relief annealing treatment at 550 ℃ for 2 h under a vacuum atmosphere to eliminate residual stresses and improve the reliability of subsequent characterization. The interfacial microstructure and chemical characteristics were systematically analyzed using multiple characterization techniques. Scanning electron microscopy (SEM) was used to observe interfacial morphology, including wave formation and defect distribution. Energy dispersive spectroscopy (EDS) was applied to examine elemental distribution across the bonding interface. Electron backscatter diffraction (EBSD) was employed to evaluate grain structure, grain boundary characteristics, recrystallization fraction, and texture evolution near the interface. Electron probe microanalysis (EPMA) was conducted to quantitatively determine the chemical composition of the interfacial melting zone and vortex regions, with both mapping and point analyses performed to identify the constituent phases. Mechanical properties at −70 ℃ were evaluated through tensile tests, Charpy impact tests, and three-point bending tests, all carried out in accordance with relevant national standards, with the results reported as mean values ± standard deviations. Results show that the vacuum environment significantly improves interfacial quality. With increasing vacuum degree, the interfacial wave becomes finer, more continuous, and more uniform, while the thickness of the molten layer decreases, and the presence of defects and brittle intermetallic compounds is reduced. EBSD analysis reveals grain refinement and an increased recrystallization fraction at the interface under vacuum conditions, demonstrating that the interfacial microstructure can be effectively controlled by adjusting the ambient pressure. EPMA further reveals that the vortex regions are mainly composed of Fe and Ti elements, while the weld seam region exhibits a stable chemical composition with Ti-to-Fe atom number ratios close to 1∶1 or 2∶1, indicating that the dominant intermetallic phases formed at the interface are TiFe and TiFe2. Benefiting from the optimized interfacial structure, the clad plates exhibit excellent mechanical properties at −70 ℃. The tensile strengths of the CP-20, CP-60 and CP-100 are 880, 911, and 867 MPa, respectively. The impact absorbed energies are 17.5, 10.2, and 6.4 J, respectively. And the flexural strengths are
, Available online , doi: 10.11883/bzycj-2026-0022
Abstract:
Current transformers are core components of power systems used for current measurement and relay protection. During long-term service under coupled multiphysical-field conditions, the internal insulation structure is susceptible to breakdown under locally intensified electric fields, triggering arc discharge in the insulating oil. The resulting oil cracking and rapid gas expansion produce a sharp rise in internal pressure. If this pressure is not relieved promptly, the resulting pressure surge may lead to combustion and explosion accidents. Therefore, the rapid-response capability of pressure relief devices under explosive loading is critical to the operational safety of current transformers. This study focuses on the expander–rupture disc pressure relief assembly of an LVB-220 current transformer and systematically investigates its dynamic mechanical behavior and failure mechanism under explosive impact loading. An equivalent hydrogen–air premixed-gas explosion test platform was constructed. The test system consisted of two stainless-steel flame-acceleration tubes with an inner diameter of 168.3 mm and a total length of4250 mm, a gas-filling and mixing unit, an ignition device, and a high-speed data-acquisition system. A hydrogen-air mixture with a volume ratio of 27:75 was used as the combustible gas to reproduce the most severe explosion conditions caused by arc-generated cracking gas in transformer oil. Pressure and flame signals were recorded synchronously using a high-temperature pressure transducer and two flame detectors installed along the tube. The tested assembly consisted of a positive-arch, cross-scored rupture disc made of 316L stainless steel, with a static burst pressure of 0.2 MPa, and a nine-convolution expander. In addition, an explicit dynamic finite element model was developed using ANSYS/LS-DYNA. The Johnson-Cook constitutive model was adopted to describe the strain-rate-dependent behavior of 316L stainless steel. The cross-scored region of the rupture disc was finely meshed with a minimum element size of approximately 0.02 mm to accurately capture crack initiation and propagation. The measured pressure–time history was applied as the loading boundary condition. The model was validated by comparing the simulated opening pressure and expander deformation with the experimental results. The results show that, during the early stage of the explosion, the pressure wave reaches the relief port ahead of the flame front and triggers the opening of the rupture disc. The measured peak pressure before disc opening is 0.72 MPa, which is significantly higher than the static calibration value of 0.2 MPa because of strain-rate hardening, structural inertia, and the spatial nonuniformity of the pressure field under rapid loading. The expander exhibits limited plastic deformation, with the axial deformation of the nine convolutions ranging from 0.60 to 2.30 mm. The deformation distribution is characterized by larger values at both ends and smaller values in the middle, reflecting an energy-absorption mechanism dominated by low-order bending modes. The asymmetric curling of the rupture-disc petals is closely associated with stress-wave reflections and fluid–structure interaction induced by high-speed venting. The numerical simulation reproduces the crack-growth sequence, opening morphology, and stress distribution in good agreement with the experimental observations. The proposed integrated experimental-numerical approach provides a reliable mechanical basis and engineering guidance for optimizing the blast-resistant design of pressure relief devices for current transformers.
Current transformers are core components of power systems used for current measurement and relay protection. During long-term service under coupled multiphysical-field conditions, the internal insulation structure is susceptible to breakdown under locally intensified electric fields, triggering arc discharge in the insulating oil. The resulting oil cracking and rapid gas expansion produce a sharp rise in internal pressure. If this pressure is not relieved promptly, the resulting pressure surge may lead to combustion and explosion accidents. Therefore, the rapid-response capability of pressure relief devices under explosive loading is critical to the operational safety of current transformers. This study focuses on the expander–rupture disc pressure relief assembly of an LVB-220 current transformer and systematically investigates its dynamic mechanical behavior and failure mechanism under explosive impact loading. An equivalent hydrogen–air premixed-gas explosion test platform was constructed. The test system consisted of two stainless-steel flame-acceleration tubes with an inner diameter of 168.3 mm and a total length of
, Available online , doi: 10.11883/bzycj-2025-0042
Abstract:
To address the challenges posed by insufficient statistical sampling density and inherent irreducible uncertainties in multi-physical property detonation experiments, probability learning on manifold (PLoM) involving diffusion map and Itô projection sampling are used to generate sufficient dataset satisfying the detonation physical mechanism and therefore, to fulfill the uncertainty quantification of experiment. Firstly, scale transformation is implemented on the experimental data with multi-physical asset of insensitive high explosive PBX 9502. The training set is obtained through the normalization of the scale matrix by means of principal component analysis. Secondly, an improved high-dimensional Gaussian kernel density estimation is utilized to calibrate the probability measure of the random matrix associated with the training dataset. Diffusion map is used to deduce the nonlinear manifold based on the training dataset. Sampling on the manifold is fulfilled through Itô-MCMC generator defined by a dissipative Hamilton system driven by the Wiener process. Finally, the learning set is obtained via inverse transformation. The result shows that the Gaussian statistics obtained from random numbers generated by PLoM coincide with the statistical information of density of PBX9502 calibrated in the literature. Furthermore, the double logarithm model related to the distance to detonation and initial impact stress is constructed through the data generated. It also holds for the relationship between the time of detonation and initial shock stress. Fitting precision of the curve is almost equivalent to the accuracy of result in literature, however the cost is negligible. More accurate digital test result is obtained through the learning and processing of existing experimental data via PLoM. The PLoM method demonstrates strong generalization capability, enabling its extension to detonation experiments with various types of explosive.
To address the challenges posed by insufficient statistical sampling density and inherent irreducible uncertainties in multi-physical property detonation experiments, probability learning on manifold (PLoM) involving diffusion map and Itô projection sampling are used to generate sufficient dataset satisfying the detonation physical mechanism and therefore, to fulfill the uncertainty quantification of experiment. Firstly, scale transformation is implemented on the experimental data with multi-physical asset of insensitive high explosive PBX 9502. The training set is obtained through the normalization of the scale matrix by means of principal component analysis. Secondly, an improved high-dimensional Gaussian kernel density estimation is utilized to calibrate the probability measure of the random matrix associated with the training dataset. Diffusion map is used to deduce the nonlinear manifold based on the training dataset. Sampling on the manifold is fulfilled through Itô-MCMC generator defined by a dissipative Hamilton system driven by the Wiener process. Finally, the learning set is obtained via inverse transformation. The result shows that the Gaussian statistics obtained from random numbers generated by PLoM coincide with the statistical information of density of PBX
, Available online , doi: 10.11883/bzycj-2025-0139
Abstract:
In order to compare the similarities and differences among the ‘self-sharpening’ behavior as well as the corresponding ballistic performance of tungsten fiber / metallic glass matrix (WF/MG) composite long rod projectile with different fiber diameters, the deformation and failure characteristics of various WF/MG composite long rod projectile were investigated systematically integrated with related penetration test results and mesoscopic finite element (FE) simulation analysis. The two-dimensional (2D) and three-dimensional (3D) meso-scale FE geometric models of WF/MG composite long rod projectile were established according to the actual distribution characteristics of metallic glass matrix and tungsten fibers. A coupled ‘thermo-mechanical-free volume’ constitutive model was applied to characterize the high shear sensitivity as well as high strength of metallic glass matrix. Moreover, the transformation of deformation and failure modes of projectiles under different impact velocities was analyzed in detail. Related analysis demonstrated that due to the differences in the mechanical properties among the tungsten fibers with different diameters, i.e., bending resistance and shear resistance, etc., the corresponding composite long rod projectile exhibit different failure characteristics during the penetration process, which further affects their penetration/perforation performance. The impact velocity also plays an important role on the deformation and failure modes of composite long rod projectile and the corresponding penetration / perforation performance. When the impact velocity is relatively low, the fibers in the nose of composite rod projectile exhibits unstable buckling phenomenon and gradually disperse during the penetration process, which makes the projectile nose blunt at a certain extent. This induces the increase in the penetration resistance and decrease in the penetration performance. As the impact velocity increases, small-diameter fibers reflux along the nose periphery while large-diameter fibers fail in shear. The large-diameter fiber further strengthens the ‘self-sharpening’ behavior as well as the penetration performance of projectile. When the impact velocity and fiber diameter reach a certain upper threshold values, the thickness of ‘edge layer’ in the rod nose is sharply reduced, and thus the ‘self-sharpening’ performance is weakened, and the penetration capability of projectile is decayed again. Related research is beneficial to predicting the penetration / perforation performance of WF/MG composite long rod projectiles with different fiber diameters under different impact velocities, and optimizing the structural design of projectile as well as the impact velocity, etc.
In order to compare the similarities and differences among the ‘self-sharpening’ behavior as well as the corresponding ballistic performance of tungsten fiber / metallic glass matrix (WF/MG) composite long rod projectile with different fiber diameters, the deformation and failure characteristics of various WF/MG composite long rod projectile were investigated systematically integrated with related penetration test results and mesoscopic finite element (FE) simulation analysis. The two-dimensional (2D) and three-dimensional (3D) meso-scale FE geometric models of WF/MG composite long rod projectile were established according to the actual distribution characteristics of metallic glass matrix and tungsten fibers. A coupled ‘thermo-mechanical-free volume’ constitutive model was applied to characterize the high shear sensitivity as well as high strength of metallic glass matrix. Moreover, the transformation of deformation and failure modes of projectiles under different impact velocities was analyzed in detail. Related analysis demonstrated that due to the differences in the mechanical properties among the tungsten fibers with different diameters, i.e., bending resistance and shear resistance, etc., the corresponding composite long rod projectile exhibit different failure characteristics during the penetration process, which further affects their penetration/perforation performance. The impact velocity also plays an important role on the deformation and failure modes of composite long rod projectile and the corresponding penetration / perforation performance. When the impact velocity is relatively low, the fibers in the nose of composite rod projectile exhibits unstable buckling phenomenon and gradually disperse during the penetration process, which makes the projectile nose blunt at a certain extent. This induces the increase in the penetration resistance and decrease in the penetration performance. As the impact velocity increases, small-diameter fibers reflux along the nose periphery while large-diameter fibers fail in shear. The large-diameter fiber further strengthens the ‘self-sharpening’ behavior as well as the penetration performance of projectile. When the impact velocity and fiber diameter reach a certain upper threshold values, the thickness of ‘edge layer’ in the rod nose is sharply reduced, and thus the ‘self-sharpening’ performance is weakened, and the penetration capability of projectile is decayed again. Related research is beneficial to predicting the penetration / perforation performance of WF/MG composite long rod projectiles with different fiber diameters under different impact velocities, and optimizing the structural design of projectile as well as the impact velocity, etc.
, Available online , doi: 10.11883/bzycj-2025-0294
Abstract:
Magnesium powder, a commonly utilized metallic material, frequently gives rise to combustion and explosion incidents during its production. To identify efficient suppressants for magnesium dust explosions, six bimetallic supramolecular compounds were synthesized via co-precipitation, and their flame-suppression performance on magnesium dust flames was compared by means of flame morphology and flame-propagation velocity in a Hartmann tube; sodium bicarbonate, a traditional suppressant, served as the reference. In the tests, magnesium powder was premixed with each of the seven suppressants at a mass ratio of 1:6 and then ignited in the Hartmann tube. The average flame-propagation velocity was derived from the tube length and the time required for the flame to reach the tube top, whereas the instantaneous flame-propagation velocity was obtained by recording flame height at successive instants. A higher average velocity and a faster flame-front propagation velocity signify a weaker suppressant efficacy. When the seven suppressants were ranked according to these two indicators, however, the orders were not fully consistent. For instance, MgAl-Cl yielded a lower average velocity than NaHCO3, yet its maximum flame-propagation velocity exhibited the opposite trend. This discrepancy arises because certain suppressants only weakly inhibit the magnesium flame, allowing rapid flame acceleration near the tube outlet where oxygen is abundant. Integrating both the average and the maximum flame-propagation velocities, the suppressants can be arranged from the least to the most effective as follows: NaHCO3, MgAl-Cl, ZnCr-CO3, CaFe-Cl, MgAl-CO3, CuAl-CO3, CaFe-CO3. Moreover, the percentage reductions in average velocity imparted by the four bimetallic supramolecular compounds (MgAl-Cl versus MgAl-CO3, and CaFe-Cl versus CaFe-CO3) were 17.95 %, 27.30 %, 23.82 %, and 50.76 %, respectively, evidencing that carbonate as the interlayer anion have a superior suppression effect on magnesium powder deflagration compared to those with chloride as the interlayer anion. SEM and XRD analyses of the reaction products, together with TG-DSC traces of the bimetallic supramolecular compounds, revealed that during the decomposition process, bimetallic supramolecular compounds lower the flame temperature via the desorption of interlayer water molecules and the heat absorption associated with the decomposition of the layered structure. Moreover, the inert gases and metal oxides generated during decomposition can block heat transfer and inhibit the volatilization of combustible gases from the surface of magnesium powder particles. Meanwhile, the metal ions and interlayer anions participate in the combustion reaction, consuming free radicals and interrupting the chain reaction, thereby achieving the explosion suppression effect.
Magnesium powder, a commonly utilized metallic material, frequently gives rise to combustion and explosion incidents during its production. To identify efficient suppressants for magnesium dust explosions, six bimetallic supramolecular compounds were synthesized via co-precipitation, and their flame-suppression performance on magnesium dust flames was compared by means of flame morphology and flame-propagation velocity in a Hartmann tube; sodium bicarbonate, a traditional suppressant, served as the reference. In the tests, magnesium powder was premixed with each of the seven suppressants at a mass ratio of 1:6 and then ignited in the Hartmann tube. The average flame-propagation velocity was derived from the tube length and the time required for the flame to reach the tube top, whereas the instantaneous flame-propagation velocity was obtained by recording flame height at successive instants. A higher average velocity and a faster flame-front propagation velocity signify a weaker suppressant efficacy. When the seven suppressants were ranked according to these two indicators, however, the orders were not fully consistent. For instance, MgAl-Cl yielded a lower average velocity than NaHCO3, yet its maximum flame-propagation velocity exhibited the opposite trend. This discrepancy arises because certain suppressants only weakly inhibit the magnesium flame, allowing rapid flame acceleration near the tube outlet where oxygen is abundant. Integrating both the average and the maximum flame-propagation velocities, the suppressants can be arranged from the least to the most effective as follows: NaHCO3, MgAl-Cl, ZnCr-CO3, CaFe-Cl, MgAl-CO3, CuAl-CO3, CaFe-CO3. Moreover, the percentage reductions in average velocity imparted by the four bimetallic supramolecular compounds (MgAl-Cl versus MgAl-CO3, and CaFe-Cl versus CaFe-CO3) were 17.95 %, 27.30 %, 23.82 %, and 50.76 %, respectively, evidencing that carbonate as the interlayer anion have a superior suppression effect on magnesium powder deflagration compared to those with chloride as the interlayer anion. SEM and XRD analyses of the reaction products, together with TG-DSC traces of the bimetallic supramolecular compounds, revealed that during the decomposition process, bimetallic supramolecular compounds lower the flame temperature via the desorption of interlayer water molecules and the heat absorption associated with the decomposition of the layered structure. Moreover, the inert gases and metal oxides generated during decomposition can block heat transfer and inhibit the volatilization of combustible gases from the surface of magnesium powder particles. Meanwhile, the metal ions and interlayer anions participate in the combustion reaction, consuming free radicals and interrupting the chain reaction, thereby achieving the explosion suppression effect.
, Available online , doi: 10.11883/bzycj-2026-0132
Abstract:
Accurate sound velocity measurement under high pressure requires high-impedance window materials. However, no ideal candidates have been reported to date. In this study, the Hugoniot data and sound velocity of single-crystal lead fluoride (PbF2) in the pressure range from 64 GPa to 198 GPa were investigated by applying shock compression techniques and quantum molecular dynamics simulations. The theoretical simulation results of Hugoniot data and sound velocities are in excellent agreement with the experimental data. Along the Hugoniot curve, the sound velocities of liquid FbF2 increase linear with particle velocity. The signal from a displacement interferometer system for any reflector (DISAR) and simulation results confirm that PbF2 transforms into a liquid metal at 89.6 GPa, which enables efficient laser reflection at the shock front. Using a DISAR, both the shock wave velocity and the catch-up time of rarefaction waves can be accurately measured, providing critical technical support for the precise determination of sound velocity in materials under extremely high pressure. These superior properties, including high impedance, structural stability under high pressure, and pressure-induced metallization, enable PbF2 fully meet the core requirements for high-impedance standard windows.
Accurate sound velocity measurement under high pressure requires high-impedance window materials. However, no ideal candidates have been reported to date. In this study, the Hugoniot data and sound velocity of single-crystal lead fluoride (PbF2) in the pressure range from 64 GPa to 198 GPa were investigated by applying shock compression techniques and quantum molecular dynamics simulations. The theoretical simulation results of Hugoniot data and sound velocities are in excellent agreement with the experimental data. Along the Hugoniot curve, the sound velocities of liquid FbF2 increase linear with particle velocity. The signal from a displacement interferometer system for any reflector (DISAR) and simulation results confirm that PbF2 transforms into a liquid metal at 89.6 GPa, which enables efficient laser reflection at the shock front. Using a DISAR, both the shock wave velocity and the catch-up time of rarefaction waves can be accurately measured, providing critical technical support for the precise determination of sound velocity in materials under extremely high pressure. These superior properties, including high impedance, structural stability under high pressure, and pressure-induced metallization, enable PbF2 fully meet the core requirements for high-impedance standard windows.
, Available online , doi: 10.11883/bzycj-2025-0342
Abstract:
Efficiently achieving plane wave loading on the test specimen is a key technical issue in the design of blast wave simulators. Based on the existing device, testing data of large-scale shock tubes, and LS-DYNA software, numerical models of blast wave propagation in simulators were established. Numerical simulations were conducted to develop a safe, economical, and reusable blast wave simulator. A quantitative method for assessing the uniformity of overpressure load on the loading area was proposed. Numerical analyses were performed to investigate the influence of the shape and length of the expansion section and the length of the conditioning section on the uniformity of overpressure load on the loading area. Key structural parameters, including the wall thickness of the expansion and conditioning sections, as well as the spacing and height of the stiffeners, were numerically optimized. Validation tests were conducted on the blast wave simulator. It is found that the established numerical model can accurately reproduce the blast wave propagation, and the prediction results show good agreement with the testing data. Taking the errors of overpressure peak value and arrival time as indicators, the quantitative evaluation of overpressure load uniformity on the loading area is achieved. Considering the balance between technical and economic factors, the developed blast wave simulator is designed with a symmetrically configured expansion section 3 m in length. The length of the conditioning section can be extended as much as the actual investment allows. Based on the numerical results, the wall thickness of both the expansion and conditioning sections is determined to be 30 mm. The height and spacing of the stiffeners are both recommended to be 150 mm. The results of validation tests confirm that the geometric and structural design meets the requirements of blast loads and structural resistance, and the developed simulator is suitable for component-level tests.
Efficiently achieving plane wave loading on the test specimen is a key technical issue in the design of blast wave simulators. Based on the existing device, testing data of large-scale shock tubes, and LS-DYNA software, numerical models of blast wave propagation in simulators were established. Numerical simulations were conducted to develop a safe, economical, and reusable blast wave simulator. A quantitative method for assessing the uniformity of overpressure load on the loading area was proposed. Numerical analyses were performed to investigate the influence of the shape and length of the expansion section and the length of the conditioning section on the uniformity of overpressure load on the loading area. Key structural parameters, including the wall thickness of the expansion and conditioning sections, as well as the spacing and height of the stiffeners, were numerically optimized. Validation tests were conducted on the blast wave simulator. It is found that the established numerical model can accurately reproduce the blast wave propagation, and the prediction results show good agreement with the testing data. Taking the errors of overpressure peak value and arrival time as indicators, the quantitative evaluation of overpressure load uniformity on the loading area is achieved. Considering the balance between technical and economic factors, the developed blast wave simulator is designed with a symmetrically configured expansion section 3 m in length. The length of the conditioning section can be extended as much as the actual investment allows. Based on the numerical results, the wall thickness of both the expansion and conditioning sections is determined to be 30 mm. The height and spacing of the stiffeners are both recommended to be 150 mm. The results of validation tests confirm that the geometric and structural design meets the requirements of blast loads and structural resistance, and the developed simulator is suitable for component-level tests.
, Available online , doi: 10.11883/bzycj-2025-0355
Abstract:
The design of blast wave simulators for prototype or large-scale engineering structures and components must carefully balance safety, cost, and performance. To address this challenge, a blast wave simulator driven by methane–air deflagration was developed in this study. Numerical models of methane–air deflagration in large-scale tubes were established using the CFD software OpenFOAM and validated against experimental data. On the basis of the validated model, numerical simulations were performed to investigate the effects of gas cloud length, venting conditions, and obstacle configurations on overpressure loads, thereby revealing the mechanisms of blast wave and flame propagation in the methane–air deflagration-driven simulator. A testing scheme was subsequently proposed to generate different levels of blast loading. The results show that the developed numerical model can reasonably predict the overpressure time histories and spatial distributions of gas explosions in large-scale tubes. In the simulator, increasing the gas cloud volume and placing obstacles near the ignition zone of the driven section can significantly intensify the deflagration of methane–air mixtures. Moreover, side venting effectively separates the pressure wave from the flame, thereby preventing high-temperature effects in the loading area. An optimized loading configuration was proposed, in which four obstacles (1×50% blockage ratio and 3×25% blockage ratio, spaced at 3 m intervals) were arranged, and different blast load levels were achieved by varying the gas cloud length from 1.5 m to 6 m in the driven section. Trial tests of the developed simulator showed that the measured peak overpressures agree well with the numerical predictions, with deviations less than 12%, and the results exhibit good uniformity and high repeatability. These findings demonstrate that the proposed methane-air deflagration-driven blast wave simulator is suitable for blast testing of structural components such as reinforced concrete slabs.
The design of blast wave simulators for prototype or large-scale engineering structures and components must carefully balance safety, cost, and performance. To address this challenge, a blast wave simulator driven by methane–air deflagration was developed in this study. Numerical models of methane–air deflagration in large-scale tubes were established using the CFD software OpenFOAM and validated against experimental data. On the basis of the validated model, numerical simulations were performed to investigate the effects of gas cloud length, venting conditions, and obstacle configurations on overpressure loads, thereby revealing the mechanisms of blast wave and flame propagation in the methane–air deflagration-driven simulator. A testing scheme was subsequently proposed to generate different levels of blast loading. The results show that the developed numerical model can reasonably predict the overpressure time histories and spatial distributions of gas explosions in large-scale tubes. In the simulator, increasing the gas cloud volume and placing obstacles near the ignition zone of the driven section can significantly intensify the deflagration of methane–air mixtures. Moreover, side venting effectively separates the pressure wave from the flame, thereby preventing high-temperature effects in the loading area. An optimized loading configuration was proposed, in which four obstacles (1×50% blockage ratio and 3×25% blockage ratio, spaced at 3 m intervals) were arranged, and different blast load levels were achieved by varying the gas cloud length from 1.5 m to 6 m in the driven section. Trial tests of the developed simulator showed that the measured peak overpressures agree well with the numerical predictions, with deviations less than 12%, and the results exhibit good uniformity and high repeatability. These findings demonstrate that the proposed methane-air deflagration-driven blast wave simulator is suitable for blast testing of structural components such as reinforced concrete slabs.
, Available online , doi: 10.11883/bzycj-2025-0352
Abstract:
During service, numerous structural metals in nuclear reactors are subjected to high-energy particle irradiation under prestressed conditions and may also experience unpredictable impact loads. These conditions can lead to significant degradation of the mechanical properties of structural metals, including irradiation hardening, irradiation embrittlement, irradiation swelling, and irradiation creep, thereby reducing reactor safety and reliability. However, existing studies have predominantly focused on the effects of either irradiation or impact loading on the mechanical properties and deformation mechanisms of metals. A limited number of studies have investigated the degradation behavior of irradiated metals under dynamic loading, whereas the coupled effect of prestress, which is a typical reactor service condition, has been considered even less frequently. The current status of experimental research on the hardening and embrittlement behavior and underlying micromechanisms of metals subjected to the coupled effects of irradiation, prestress, and impact loading is reviewed and summarized. Shortcomings and major challenges in current research are identified, and perspectives and recommendations are proposed regarding the key scientific issues requiring further attention and the critical technological bottlenecks that must be overcome in future studies. Specifically, the following issues should be addressed: (1) the necessity and sufficiency of the causal relationship between dislocation channels and irradiation-induced hardening and embrittlement mechanisms require thorough investigation; (2) equivalent dose-rate assessment methods for metal hardening and embrittlement induced by different irradiation particle sources need to be strengthened; (3) the hardening and embrittlement mechanisms of structural metals in nuclear reactors under prestress and neutron irradiation conditions require further investigation; and (4) research on the impact behavior of prestressed and neutron-irradiated metals urgently needs to be initiated, while safe and efficient experimental techniques for characterizing the dynamic properties of neutron-irradiated metals must be developed immediately. This review is expected to provide a scientific basis and research methodology for future investigations into the hardening and embrittlement behavior and mechanisms of metals under irradiation-prestress conditions, as required in major national projects such as nuclear reactor life extension and the development of advanced reactors.
During service, numerous structural metals in nuclear reactors are subjected to high-energy particle irradiation under prestressed conditions and may also experience unpredictable impact loads. These conditions can lead to significant degradation of the mechanical properties of structural metals, including irradiation hardening, irradiation embrittlement, irradiation swelling, and irradiation creep, thereby reducing reactor safety and reliability. However, existing studies have predominantly focused on the effects of either irradiation or impact loading on the mechanical properties and deformation mechanisms of metals. A limited number of studies have investigated the degradation behavior of irradiated metals under dynamic loading, whereas the coupled effect of prestress, which is a typical reactor service condition, has been considered even less frequently. The current status of experimental research on the hardening and embrittlement behavior and underlying micromechanisms of metals subjected to the coupled effects of irradiation, prestress, and impact loading is reviewed and summarized. Shortcomings and major challenges in current research are identified, and perspectives and recommendations are proposed regarding the key scientific issues requiring further attention and the critical technological bottlenecks that must be overcome in future studies. Specifically, the following issues should be addressed: (1) the necessity and sufficiency of the causal relationship between dislocation channels and irradiation-induced hardening and embrittlement mechanisms require thorough investigation; (2) equivalent dose-rate assessment methods for metal hardening and embrittlement induced by different irradiation particle sources need to be strengthened; (3) the hardening and embrittlement mechanisms of structural metals in nuclear reactors under prestress and neutron irradiation conditions require further investigation; and (4) research on the impact behavior of prestressed and neutron-irradiated metals urgently needs to be initiated, while safe and efficient experimental techniques for characterizing the dynamic properties of neutron-irradiated metals must be developed immediately. This review is expected to provide a scientific basis and research methodology for future investigations into the hardening and embrittlement behavior and mechanisms of metals under irradiation-prestress conditions, as required in major national projects such as nuclear reactor life extension and the development of advanced reactors.
, Available online , doi: 10.11883/bzycj-2025-0419
Abstract:
Addressing the high probability of injury and fatality associated with traditionalanti-riot kinetic energy projectiles represented by rubber bullets, a novel low-kinetic-energy frangible projectile with the characteristics of reliable shell fragmentation and the energy release of powder dispersion was designed based on the existing launch platforms of non-lethal weapons, aimed to improve the application safety of anti-riot kinetic energy projectiles. Through shooting tests on anthropomorphic dummy targets and ballistic gelatin targets, the dynamic impact responses of the two targets under the blunt impacts of frangible projectiles were obtained. In the dummy target tests, the frangible projectile achieved reliable fragmentation and exhibited a significant kinetic energy slow-release effect; the calculated injury assessment indicators for different parts of the dummy were all within the safety threshold. Comparative analyses with the blunt impact test results of rubber bullets show that this type of frangible projectile has significantly higher safety than rubber bullets while maintaining effective pain-inducing performance. In the ballistic gelatin target shooting tests, under the same projectile mass and impact velocity, the peak pressure generated in the gelatin measured by the internal sensor differed by approximately one order of magnitude between the frangible projectile and the rubber bullet. The test results show that the novel low-kinetic-energy frangible projectile can effectively attenuate the peak impact pressure of the projectile and greatly reduce its injury potential through the fragmentation of the projectile shell and the dispersion of the powder filled inside the projectile, which provides technical support for improving the application safety of anti-riot kinetic energy ammunition.
Addressing the high probability of injury and fatality associated with traditionalanti-riot kinetic energy projectiles represented by rubber bullets, a novel low-kinetic-energy frangible projectile with the characteristics of reliable shell fragmentation and the energy release of powder dispersion was designed based on the existing launch platforms of non-lethal weapons, aimed to improve the application safety of anti-riot kinetic energy projectiles. Through shooting tests on anthropomorphic dummy targets and ballistic gelatin targets, the dynamic impact responses of the two targets under the blunt impacts of frangible projectiles were obtained. In the dummy target tests, the frangible projectile achieved reliable fragmentation and exhibited a significant kinetic energy slow-release effect; the calculated injury assessment indicators for different parts of the dummy were all within the safety threshold. Comparative analyses with the blunt impact test results of rubber bullets show that this type of frangible projectile has significantly higher safety than rubber bullets while maintaining effective pain-inducing performance. In the ballistic gelatin target shooting tests, under the same projectile mass and impact velocity, the peak pressure generated in the gelatin measured by the internal sensor differed by approximately one order of magnitude between the frangible projectile and the rubber bullet. The test results show that the novel low-kinetic-energy frangible projectile can effectively attenuate the peak impact pressure of the projectile and greatly reduce its injury potential through the fragmentation of the projectile shell and the dispersion of the powder filled inside the projectile, which provides technical support for improving the application safety of anti-riot kinetic energy ammunition.
, Available online , doi: 10.11883/bzycj-2026-0035
Abstract:
Aiming at the problems of complex characteristics of blast load induced by fuel-rich explosives in confined spaces, as well as the high cost of high-precision solution and engineering application, a numerical calculation method for confined blast loading considering the afterburning effect was established. The method was verified for its reliability by comparing the calculated quasi-static pressure, impulse within the saturation response time and residual deformation of the target with the experimental results, with all relative errors controlled within 10%. The spatiotemporal distribution law of confined blast loading in confined spaces was systematically analyzed, and a simplified equivalent loading method considering both the saturation response time and quasi-static pressure was proposed. The reliability of the simplification method was validated by comparing the calculated first peak deformation and residual deformation at the central point with the results from fully coupled calculation, with all errors within 10%. Through an investigation into the spatial distribution form of the equivalent load and the influence of quasi-static pressure on structural response, the results show that within the scope of the current research, the spatial distribution of the equivalent load has a relatively minor effect on structural response, while the contribution of quasi-static pressure cannot be neglected. On the basis of the above findings, a two-stage load simplification model based on the load characteristics at the central point of the target plate was proposed. The reliability of the simplified model was confirmed by comparing the residual deformation values obtained from 10 groups of simplified model calculations with the experimental data, with all errors within 15%. The research results indicate that the proposed model exhibits good applicability under different working conditions; it can significantly improve the calculation efficiency while ensuring the calculation accuracy, and thus provides a technical approach for the simplified analysis of engineering problems related to confined space explosions.
Aiming at the problems of complex characteristics of blast load induced by fuel-rich explosives in confined spaces, as well as the high cost of high-precision solution and engineering application, a numerical calculation method for confined blast loading considering the afterburning effect was established. The method was verified for its reliability by comparing the calculated quasi-static pressure, impulse within the saturation response time and residual deformation of the target with the experimental results, with all relative errors controlled within 10%. The spatiotemporal distribution law of confined blast loading in confined spaces was systematically analyzed, and a simplified equivalent loading method considering both the saturation response time and quasi-static pressure was proposed. The reliability of the simplification method was validated by comparing the calculated first peak deformation and residual deformation at the central point with the results from fully coupled calculation, with all errors within 10%. Through an investigation into the spatial distribution form of the equivalent load and the influence of quasi-static pressure on structural response, the results show that within the scope of the current research, the spatial distribution of the equivalent load has a relatively minor effect on structural response, while the contribution of quasi-static pressure cannot be neglected. On the basis of the above findings, a two-stage load simplification model based on the load characteristics at the central point of the target plate was proposed. The reliability of the simplified model was confirmed by comparing the residual deformation values obtained from 10 groups of simplified model calculations with the experimental data, with all errors within 15%. The research results indicate that the proposed model exhibits good applicability under different working conditions; it can significantly improve the calculation efficiency while ensuring the calculation accuracy, and thus provides a technical approach for the simplified analysis of engineering problems related to confined space explosions.
, Available online , doi: 10.11883/bzycj-2025-0386
Abstract:
To investigate the damage characteristics of concrete arch slabs subjected to underwater contact explosions, a series of underwater contact explosion tests was conducted on three concrete arch slabs. During the experiments, the spalled fragments generated by detonation were carefully collected and quantitatively analyzed using a mass-based partitioning statistical method, thereby enabling systematic characterization of the slabs’ fragmentation and spalling behavior. On this basis, a three-dimensional fluid–structure interaction numerical model was developed using the Coupled Eulerian-Lagrangian (CEL) method, which incorporates water, a concrete arch slab, and TNT. The reliability and effectiveness of the numerical model were validated by comparing the simulated macroscopic failure patterns and crack propagation paths with the experimental observations. Further analyses were conducted using the validated model to investigate the damage evolution process and the characteristics of stress wave propagation within the arch slabs under explosive loading, thereby systematically elucidating the damage mechanisms of concrete arch slabs subjected to underwater contact explosions. In addition, parametric studies were performed to examine the effects of explosive charge and span length on the damage characteristics of the arch slabs. The results indicate that, under underwater contact explosion, an annular damage zone composed of radial and circumferential cracks forms on the blast-facing surface of the arch slab. With increasing explosive charge, the number of spalled fragments increases, and their size becomes progressively finer. As the span length increases, the structure tends to generate larger fragments. When the TNT charge increases from 0.4 g to 1.6 g, the annular damage zone on the blast-facing surface expands from a localized region to the full width of the slab, and the failure mode transitions from localized damage to global failure. With increasing span length, the degree of damage at the arch ends decreases, while the damaged region near the mid-span expands. Compared with the slab with a span of 400 mm, the total energy dissipation of the slab with a span of 600 mm increases by 4.3%.
To investigate the damage characteristics of concrete arch slabs subjected to underwater contact explosions, a series of underwater contact explosion tests was conducted on three concrete arch slabs. During the experiments, the spalled fragments generated by detonation were carefully collected and quantitatively analyzed using a mass-based partitioning statistical method, thereby enabling systematic characterization of the slabs’ fragmentation and spalling behavior. On this basis, a three-dimensional fluid–structure interaction numerical model was developed using the Coupled Eulerian-Lagrangian (CEL) method, which incorporates water, a concrete arch slab, and TNT. The reliability and effectiveness of the numerical model were validated by comparing the simulated macroscopic failure patterns and crack propagation paths with the experimental observations. Further analyses were conducted using the validated model to investigate the damage evolution process and the characteristics of stress wave propagation within the arch slabs under explosive loading, thereby systematically elucidating the damage mechanisms of concrete arch slabs subjected to underwater contact explosions. In addition, parametric studies were performed to examine the effects of explosive charge and span length on the damage characteristics of the arch slabs. The results indicate that, under underwater contact explosion, an annular damage zone composed of radial and circumferential cracks forms on the blast-facing surface of the arch slab. With increasing explosive charge, the number of spalled fragments increases, and their size becomes progressively finer. As the span length increases, the structure tends to generate larger fragments. When the TNT charge increases from 0.4 g to 1.6 g, the annular damage zone on the blast-facing surface expands from a localized region to the full width of the slab, and the failure mode transitions from localized damage to global failure. With increasing span length, the degree of damage at the arch ends decreases, while the damaged region near the mid-span expands. Compared with the slab with a span of 400 mm, the total energy dissipation of the slab with a span of 600 mm increases by 4.3%.
, Available online , doi: 10.11883/bzycj-2026-0056
Abstract:
To clarify the ultimate load-bearing capacity characteristics of prestressed concrete T-girder bridges under typical blast loads and address the technical difficulties in the residual load-bearing capacity evaluation of blast-damaged bridge structures, a two-stage test method combining blast damage test and post-blast static loading test is adopted in this study. A full-scale (1:1) prestressed concrete T-girder bridge specimen consistent with practical engineering structural parameters is taken as the test object to carry out large-equivalent field blast tests. Multi-stage graded static loading tests are further implemented on the blast-damaged specimen to obtain structural damage evolution features and deflection response data under post-blast service conditions. Based on the measured damage morphologies, deformation data and mechanical response laws acquired from physical tests, a refined three-dimensional finite element (FE) numerical simulation model for prestressed concrete T-girder bridges is established. Key structural characteristics including concrete and reinforcement strain-rate effects, material nonlinearity, prestress loss and structural contact interaction are fully considered in model establishment. The accuracy and reliability of the numerical model are validated through comparative calibration with experimental results. The validated model is subsequently employed to quantitatively investigate the ultimate load-bearing performance of bridges under three typical blast working conditions, including deck central contact blast, combined blast of deck central explosion and inter-girder internal explosion, and bilateral inter-girder internal explosion. The research results demonstrate distinct structural damage and bearing degradation characteristics under different blast scenarios. Contact blast acting on the bridge deck mainly induces local slab perforation failure and concentrated plastic deformation of deck concrete. For internal explosion occurring between main girders, structural damage presents a more severe and extensive pattern. In addition to concrete spalling and fragmentation on the bottom surface of the bridge deck, obvious outward bulging deformation, concrete crushing and local collapse failure occur on the webs and transverse diaphragms of T-girders. Among the three typical blast damage states, the combined blast condition with deck central explosion followed by secondary inter-girder explosion causes the most significant degradation of structural ultimate load-bearing capacity due to the superimposed synergistic damage effect. Reasonable optimization of live load layout effectively avoids severe damage areas of blast-affected bridges, eliminates the adverse influence of blast-induced structural defects, and achieves a remarkable improvement in the ultimate load-bearing capacity of damaged bridges. An integrated evaluation system covering the whole process from blast damage characterization to ultimate load-bearing capacity quantification is established in this research. The proposed evaluation method provides a systematic theoretical basis and feasible technical support for residual load-bearing capacity assessment, damage prediction and blast-resistant reinforcement design of prestressed concrete T-girder bridges suffering blast loads.
To clarify the ultimate load-bearing capacity characteristics of prestressed concrete T-girder bridges under typical blast loads and address the technical difficulties in the residual load-bearing capacity evaluation of blast-damaged bridge structures, a two-stage test method combining blast damage test and post-blast static loading test is adopted in this study. A full-scale (1:1) prestressed concrete T-girder bridge specimen consistent with practical engineering structural parameters is taken as the test object to carry out large-equivalent field blast tests. Multi-stage graded static loading tests are further implemented on the blast-damaged specimen to obtain structural damage evolution features and deflection response data under post-blast service conditions. Based on the measured damage morphologies, deformation data and mechanical response laws acquired from physical tests, a refined three-dimensional finite element (FE) numerical simulation model for prestressed concrete T-girder bridges is established. Key structural characteristics including concrete and reinforcement strain-rate effects, material nonlinearity, prestress loss and structural contact interaction are fully considered in model establishment. The accuracy and reliability of the numerical model are validated through comparative calibration with experimental results. The validated model is subsequently employed to quantitatively investigate the ultimate load-bearing performance of bridges under three typical blast working conditions, including deck central contact blast, combined blast of deck central explosion and inter-girder internal explosion, and bilateral inter-girder internal explosion. The research results demonstrate distinct structural damage and bearing degradation characteristics under different blast scenarios. Contact blast acting on the bridge deck mainly induces local slab perforation failure and concentrated plastic deformation of deck concrete. For internal explosion occurring between main girders, structural damage presents a more severe and extensive pattern. In addition to concrete spalling and fragmentation on the bottom surface of the bridge deck, obvious outward bulging deformation, concrete crushing and local collapse failure occur on the webs and transverse diaphragms of T-girders. Among the three typical blast damage states, the combined blast condition with deck central explosion followed by secondary inter-girder explosion causes the most significant degradation of structural ultimate load-bearing capacity due to the superimposed synergistic damage effect. Reasonable optimization of live load layout effectively avoids severe damage areas of blast-affected bridges, eliminates the adverse influence of blast-induced structural defects, and achieves a remarkable improvement in the ultimate load-bearing capacity of damaged bridges. An integrated evaluation system covering the whole process from blast damage characterization to ultimate load-bearing capacity quantification is established in this research. The proposed evaluation method provides a systematic theoretical basis and feasible technical support for residual load-bearing capacity assessment, damage prediction and blast-resistant reinforcement design of prestressed concrete T-girder bridges suffering blast loads.
, Available online , doi: 10.11883/bzycj-2025-0313
Abstract:
Traditional homogenization models have difficulty accurately capturing the mesoscale effects of aggregate distribution, aggregate particle size, and reinforcement configuration on stress-wave propagation paths and energy-dissipation mechanisms in concrete, thereby limiting an in-depth understanding of the blast-induced failure mechanisms of reinforced concrete slabs. To address this issue, a three-dimensional mesoscale finite element model of reinforced concrete slabs incorporating reinforcement, aggregates, and mortar matrix was established through the combined use of MATLAB and LS-DYNA. The aggregates were modeled according to actual aggregate gradation characteristics, the reinforcing bars were accurately arranged based on practical engineering layout parameters, and appropriate contact algorithms were adopted among the matrix, aggregates, and reinforcing bars to simulate interfacial effects. The model was validated against contact explosion tests, and the results show that it can predict the failure modes and crater dimensions of reinforced concrete slabs subjected to contact explosion loads with reasonable accuracy. On this basis, the effects of aggregate characteristics, including distribution pattern and particle size, and reinforcement arrangement on the blast resistance and stress-wave propagation behavior of reinforced concrete were investigated through parametric mesoscale numerical simulations. Regarding the aggregate parameters, the particle-size distribution pattern and aggregate size govern the evolution of stress waves and energy-dissipation characteristics, thereby affecting the geometric dimensions of the craters on the top surface and the spalling craters on the bottom surface of the concrete slab. When the aggregate particle size decreases from the top surface to the bottom surface, the expansion of the top-surface crater and the development of bottom-surface spalling can be effectively suppressed; in contrast, an increasing particle-size distribution aggravates surface cratering and internal spalling damage. In terms of aggregate size, the bottom-surface spalling craters of slabs containing small aggregates exhibit shallow and wide characteristics, whereas those containing large aggregates exhibit deep and narrow morphologies. Compared with aggregates, reinforcing bars exert a weaker influence on the overall failure mode and stress-wave propagation behavior of the slab. At a low reinforcement ratio, the reinforcing bars have little effect on the dynamic transmission of compressive stress peaks, whereas under high explosion loads, they suppress slab fragmentation, mitigate flexural damage, and improve the structural integrity and damage resistance of the slab.
Traditional homogenization models have difficulty accurately capturing the mesoscale effects of aggregate distribution, aggregate particle size, and reinforcement configuration on stress-wave propagation paths and energy-dissipation mechanisms in concrete, thereby limiting an in-depth understanding of the blast-induced failure mechanisms of reinforced concrete slabs. To address this issue, a three-dimensional mesoscale finite element model of reinforced concrete slabs incorporating reinforcement, aggregates, and mortar matrix was established through the combined use of MATLAB and LS-DYNA. The aggregates were modeled according to actual aggregate gradation characteristics, the reinforcing bars were accurately arranged based on practical engineering layout parameters, and appropriate contact algorithms were adopted among the matrix, aggregates, and reinforcing bars to simulate interfacial effects. The model was validated against contact explosion tests, and the results show that it can predict the failure modes and crater dimensions of reinforced concrete slabs subjected to contact explosion loads with reasonable accuracy. On this basis, the effects of aggregate characteristics, including distribution pattern and particle size, and reinforcement arrangement on the blast resistance and stress-wave propagation behavior of reinforced concrete were investigated through parametric mesoscale numerical simulations. Regarding the aggregate parameters, the particle-size distribution pattern and aggregate size govern the evolution of stress waves and energy-dissipation characteristics, thereby affecting the geometric dimensions of the craters on the top surface and the spalling craters on the bottom surface of the concrete slab. When the aggregate particle size decreases from the top surface to the bottom surface, the expansion of the top-surface crater and the development of bottom-surface spalling can be effectively suppressed; in contrast, an increasing particle-size distribution aggravates surface cratering and internal spalling damage. In terms of aggregate size, the bottom-surface spalling craters of slabs containing small aggregates exhibit shallow and wide characteristics, whereas those containing large aggregates exhibit deep and narrow morphologies. Compared with aggregates, reinforcing bars exert a weaker influence on the overall failure mode and stress-wave propagation behavior of the slab. At a low reinforcement ratio, the reinforcing bars have little effect on the dynamic transmission of compressive stress peaks, whereas under high explosion loads, they suppress slab fragmentation, mitigate flexural damage, and improve the structural integrity and damage resistance of the slab.
, Available online , doi: 10.11883/bzycj-2026-0054
Abstract:
Impact fatigue refers to the phenomenon in which materials or structures, subjected to repeated impact loading, experience localized stress concentrations and rapid strain accumulation, leading to the initiation of internal micro-damage and ultimately the fracture failure. Impact fatigue loads are characterized by their short duration, rapid loading rates and significantly elevated strain-rates, which has greater threat than conventional fatigue. The wheel-rail dynamic contact forces of high-speed trains exhibit typical characteristics of impact fatigue loading, which induces the accumulation of impact fatigue damage, accelerates the deterioration of material’s mechanical properties; and consequently, compromises the operational safety of high-speed trains. The present study integrates a material-based impact fatigue damage-coupled constitutive model to develop a comprehensive three-dimensional wheel-rail rolling contact finite element model. The stress-strain states and adhesion-slip characteristics of wheel-rail rolling/sliding contact are clarified, and the distribution features and accumulation evolution of wheel-rail impact fatigue damage are analyzed. Meanwhile, the effects of train speed, friction coefficient, and traction coefficient on impact fatigue damage are studied, and the influence of material constitutive model on typical wheel-rail contact mechanical behavior is examined. The results clearly indicate that the proposed impact fatigue model is able to well represent the wheel-rail contact responses, adhesion-slip distribution characteristics and damage accumulation. Under repeated rolling contact, the impact fatigue damage of the rail exhibits a nonlinear cumulative increasing trend with the increase of rolling cycles; however, the growth rate gradually decreases and eventually tends to stabilize approximately. Compared with the elastoplastic constitutive model, the wheel-rail contact mechanical responses predicted by the impact fatigue constitutive model are more severe and dangerous. Moreover, the coupling effect of impact fatigue damage gradually intensifies with the increase of wheel passages. These findings provide valuable theoretical insights and technical support for fatigue damage assessment and life prediction of high-speed wheel-rail systems.
Impact fatigue refers to the phenomenon in which materials or structures, subjected to repeated impact loading, experience localized stress concentrations and rapid strain accumulation, leading to the initiation of internal micro-damage and ultimately the fracture failure. Impact fatigue loads are characterized by their short duration, rapid loading rates and significantly elevated strain-rates, which has greater threat than conventional fatigue. The wheel-rail dynamic contact forces of high-speed trains exhibit typical characteristics of impact fatigue loading, which induces the accumulation of impact fatigue damage, accelerates the deterioration of material’s mechanical properties; and consequently, compromises the operational safety of high-speed trains. The present study integrates a material-based impact fatigue damage-coupled constitutive model to develop a comprehensive three-dimensional wheel-rail rolling contact finite element model. The stress-strain states and adhesion-slip characteristics of wheel-rail rolling/sliding contact are clarified, and the distribution features and accumulation evolution of wheel-rail impact fatigue damage are analyzed. Meanwhile, the effects of train speed, friction coefficient, and traction coefficient on impact fatigue damage are studied, and the influence of material constitutive model on typical wheel-rail contact mechanical behavior is examined. The results clearly indicate that the proposed impact fatigue model is able to well represent the wheel-rail contact responses, adhesion-slip distribution characteristics and damage accumulation. Under repeated rolling contact, the impact fatigue damage of the rail exhibits a nonlinear cumulative increasing trend with the increase of rolling cycles; however, the growth rate gradually decreases and eventually tends to stabilize approximately. Compared with the elastoplastic constitutive model, the wheel-rail contact mechanical responses predicted by the impact fatigue constitutive model are more severe and dangerous. Moreover, the coupling effect of impact fatigue damage gradually intensifies with the increase of wheel passages. These findings provide valuable theoretical insights and technical support for fatigue damage assessment and life prediction of high-speed wheel-rail systems.
, Available online , doi: 10.11883/bzycj-2025-0357
Abstract:
Because of the difficulty in obtaining granite materials and the high cost of conducting penetration tests on granite targets, an equivalence study between reinforced concrete and granite targets was carried out. To establish the equivalence relationship between the two target types, dimensional analysis and a modified compensation method were adopted, with the projectile residual velocity taken as the equivalence criterion, and a computational method for determining equivalent thickness was derived. Based on existing experimental data, numerical models for medium-velocity projectile penetration into reinforced concrete and granite targets were developed and validated using the LS-DYNA software. By varying the projectile impact velocity and target thickness in the numerical simulations, the similarities in damage characteristics between reinforced concrete and granite targets were systematically investigated, and the corresponding failure regions were classified. On the basis of the simulation results, specific equivalent design formulas for granite and reinforced concrete targets were obtained through data fitting. The results show that the established numerical models can accurately predict the projectile residual velocity and reproduce the failure characteristics of both target types during penetration. Compared with reinforced concrete, granite exhibits a smaller compaction zone and tunnel diameter, finer and longer cracks with higher propagation velocities, larger surface crack areas, and a greater tendency to form large spallation craters. Under identical penetration conditions, granite targets and reinforced concrete targets of equivalent thickness display similar failure characteristics, and both can be divided into five distinct failure regions. Through dimensional analysis and compensation correction, a dimensionless residual-velocity function for projectile penetration into reinforced concrete and granite targets was derived, together with an equivalent thickness formula relating the two materials. The fitted equivalent thickness coefficient between granite and reinforced concrete was determined to be 1.69966. Validation of the proposed equivalence formula indicates that the residual-velocity error between the prototype and equivalent model targets is less than 5%. These results provide a useful reference for the equivalent design of rock targets subjected to medium-velocity projectile penetration and offer a systematic methodology for substituting reinforced concrete for granite in related experimental and engineering applications.
Because of the difficulty in obtaining granite materials and the high cost of conducting penetration tests on granite targets, an equivalence study between reinforced concrete and granite targets was carried out. To establish the equivalence relationship between the two target types, dimensional analysis and a modified compensation method were adopted, with the projectile residual velocity taken as the equivalence criterion, and a computational method for determining equivalent thickness was derived. Based on existing experimental data, numerical models for medium-velocity projectile penetration into reinforced concrete and granite targets were developed and validated using the LS-DYNA software. By varying the projectile impact velocity and target thickness in the numerical simulations, the similarities in damage characteristics between reinforced concrete and granite targets were systematically investigated, and the corresponding failure regions were classified. On the basis of the simulation results, specific equivalent design formulas for granite and reinforced concrete targets were obtained through data fitting. The results show that the established numerical models can accurately predict the projectile residual velocity and reproduce the failure characteristics of both target types during penetration. Compared with reinforced concrete, granite exhibits a smaller compaction zone and tunnel diameter, finer and longer cracks with higher propagation velocities, larger surface crack areas, and a greater tendency to form large spallation craters. Under identical penetration conditions, granite targets and reinforced concrete targets of equivalent thickness display similar failure characteristics, and both can be divided into five distinct failure regions. Through dimensional analysis and compensation correction, a dimensionless residual-velocity function for projectile penetration into reinforced concrete and granite targets was derived, together with an equivalent thickness formula relating the two materials. The fitted equivalent thickness coefficient between granite and reinforced concrete was determined to be 1.69966. Validation of the proposed equivalence formula indicates that the residual-velocity error between the prototype and equivalent model targets is less than 5%. These results provide a useful reference for the equivalent design of rock targets subjected to medium-velocity projectile penetration and offer a systematic methodology for substituting reinforced concrete for granite in related experimental and engineering applications.
, Available online , doi: 10.11883/bzycj-2025-0348
Abstract:
The present investigation examines the effectiveness of a composite cushion charging structure, placed at the bottom of underwater boreholes, to mitigate bedrock damage induced by blasting excavation during cross-sea engineering. With reference to the embedded open caisson blasting project for a major cross-sea bridge, this research adopted a combined methodology of field sampling and controlled underwater explosion model tests. A piezoelectric ceramic detection system was employed to quantitatively analyze how the composite cushion, composed of iron-sand concrete layers with varying wave impedances, affects rock sample damage. Piezoelectric signals were processed via wavelet packet analysis to calculate the axial damage factor (DI) of the blast hole. Furthermore, by integrating fractal dimension theory with damage mechanics, a quantitative assessment was performed on the propagation behavior of cracks formed on the top surfaces of rock specimens under different test conditions. Results demonstrate that the implemented composite bottom cushion effectively mitigates blast-induced damage by significantly inhibiting macroscopic crack propagation, reducing surface crack density, and decreasing the axial penetration depth of cracks into the rock mass. Detailed analysis of axial damage reveals a clear trend, showing that the reduction in DI becomes more pronounced as the wave impedance of the cushion material increases. Specifically, within the borehole-affected zone (0–12 cm), the maximum achievable reduction in DI reaches 10.70%. In the adjacent bedrock zone, damage mitigation is even more evident, with DI decreasing by 95.7% to 95.8%, depending on the specific wave impedance. In conclusion, this work confirms that employing a composite cushion charging structure in underwater drilling and blasting operations can significantly alleviate blast-induced damage to foundational bedrock. The findings also establish that axial damage to the rock mass can be effectively controlled by adjusting the wave impedance of the iron-sand concrete in the cushion layer. These outcomes offer valuable insights for optimizing blasting design and enhancing damage control in marine engineering projects.
The present investigation examines the effectiveness of a composite cushion charging structure, placed at the bottom of underwater boreholes, to mitigate bedrock damage induced by blasting excavation during cross-sea engineering. With reference to the embedded open caisson blasting project for a major cross-sea bridge, this research adopted a combined methodology of field sampling and controlled underwater explosion model tests. A piezoelectric ceramic detection system was employed to quantitatively analyze how the composite cushion, composed of iron-sand concrete layers with varying wave impedances, affects rock sample damage. Piezoelectric signals were processed via wavelet packet analysis to calculate the axial damage factor (DI) of the blast hole. Furthermore, by integrating fractal dimension theory with damage mechanics, a quantitative assessment was performed on the propagation behavior of cracks formed on the top surfaces of rock specimens under different test conditions. Results demonstrate that the implemented composite bottom cushion effectively mitigates blast-induced damage by significantly inhibiting macroscopic crack propagation, reducing surface crack density, and decreasing the axial penetration depth of cracks into the rock mass. Detailed analysis of axial damage reveals a clear trend, showing that the reduction in DI becomes more pronounced as the wave impedance of the cushion material increases. Specifically, within the borehole-affected zone (0–12 cm), the maximum achievable reduction in DI reaches 10.70%. In the adjacent bedrock zone, damage mitigation is even more evident, with DI decreasing by 95.7% to 95.8%, depending on the specific wave impedance. In conclusion, this work confirms that employing a composite cushion charging structure in underwater drilling and blasting operations can significantly alleviate blast-induced damage to foundational bedrock. The findings also establish that axial damage to the rock mass can be effectively controlled by adjusting the wave impedance of the iron-sand concrete in the cushion layer. These outcomes offer valuable insights for optimizing blasting design and enhancing damage control in marine engineering projects.
, Available online , doi: 10.11883/bzycj-2025-0263
Abstract:
To evaluate the crater damage effect of cylinder charge contact explosion on steel fiber reinforced concrete (SFRC) structures, a numerical model of an SFRC target was developed by the using the smooth particle Galerkin method and a structured arbitrary Lagrange-Euler (SPG-S-ALE) fluid-structure interaction algorithm. The failure modes and damage of the SFRC target under different charge mass Q and length-to-diameter ratios l/d were investigated. Based on contact explosion theory and dimensional analysis, cratering coefficients K1 and K2 were introduced to develop predictive models for the crater diameter D and depth H as functions of the effective charge mass Qₑ. The results indicate that under the combined effects of charge mass and length-to-to-diameter ratio, the primary failure mode of the SFRC targets is cratering damage. Under constant charge mass conditions, as the l/d ratio increases from 1 to 5, both the crater diameter D and depth H decrease by approximately 50%. Within the range of effective charge Qₑ less than 16 kg, the K1 and\begin{document}$ \sqrt{{K}_{2}} $\end{document}
To evaluate the crater damage effect of cylinder charge contact explosion on steel fiber reinforced concrete (SFRC) structures, a numerical model of an SFRC target was developed by the using the smooth particle Galerkin method and a structured arbitrary Lagrange-Euler (SPG-S-ALE) fluid-structure interaction algorithm. The failure modes and damage of the SFRC target under different charge mass Q and length-to-diameter ratios l/d were investigated. Based on contact explosion theory and dimensional analysis, cratering coefficients K1 and K2 were introduced to develop predictive models for the crater diameter D and depth H as functions of the effective charge mass Qₑ. The results indicate that under the combined effects of charge mass and length-to-to-diameter ratio, the primary failure mode of the SFRC targets is cratering damage. Under constant charge mass conditions, as the l/d ratio increases from 1 to 5, both the crater diameter D and depth H decrease by approximately 50%. Within the range of effective charge Qₑ less than 16 kg, the K1 and
, Available online , doi: 10.11883/bzycj-2025-0336
Abstract:
To explore the feasibility of nano aluminum thermite as a detonator substitute for achieving stable detonation in mixed biomass blasting agent systems, and to clarify the differences in detonation effects relative to detonators, digital electronic detonators (S0), Al/CuO thermite (S1), and Al/Bi2O3 thermite (S2) were selected to detonate a mixed biomass blasting agent composed of wood powder and peanut shell powder (mass ratio 1∶1). Based on theoretical analysis, qualitative analysis and quantitative calculation were conducted for the total energy release and energy release power of different detonation methods to estimate detonation performance. Industrial explosive performance testing methods-including aluminum thermite detonation test, orthogonal tests of detonation velocity and intensity, underwater explosion tests, and blasting funnel tests-were used to systematically test and compare the detonation response characteristics of the mixed biomass blasting agent under different detonation conditions. The evolutionary laws of the explosive performance were explored from perspectives of impact effect, energy release intensity, and spatial damage effect. Results indicate that detonators and thermites belong to two distinct energy release power level systems: detonators enable instantaneous power output in the MW range, while thermites only reach the kW range. Aluminum thermite detonation is a typical high-temperature explosive energy release process with high energy density, which can achieve effective energy coupling under limited constraints and possesses reliable detonation capability. Oxygen pressure is the dominant factor affecting detonation velocity and intensity, followed by steel pipe wall thickness. The impact of detonation methods is relatively weak, and detonation methods are substitutable. Synergistic enhancement of detonation velocity and intensity can be achieved by increasing oxygen pressure and optimizing constraint conditions. The excitation efficiency of the three detonation methods follows a consistent ranking: S0 is the strongest, S2 ranks second, and S1 is slightly weaker. This ranking has been verified by shock wave parameters and crater volumes (0.33 m3, 0.24 m3, 0.21 m3). This research provides experimental support for the optimization and application of biomass blasting technology.
To explore the feasibility of nano aluminum thermite as a detonator substitute for achieving stable detonation in mixed biomass blasting agent systems, and to clarify the differences in detonation effects relative to detonators, digital electronic detonators (S0), Al/CuO thermite (S1), and Al/Bi2O3 thermite (S2) were selected to detonate a mixed biomass blasting agent composed of wood powder and peanut shell powder (mass ratio 1∶1). Based on theoretical analysis, qualitative analysis and quantitative calculation were conducted for the total energy release and energy release power of different detonation methods to estimate detonation performance. Industrial explosive performance testing methods-including aluminum thermite detonation test, orthogonal tests of detonation velocity and intensity, underwater explosion tests, and blasting funnel tests-were used to systematically test and compare the detonation response characteristics of the mixed biomass blasting agent under different detonation conditions. The evolutionary laws of the explosive performance were explored from perspectives of impact effect, energy release intensity, and spatial damage effect. Results indicate that detonators and thermites belong to two distinct energy release power level systems: detonators enable instantaneous power output in the MW range, while thermites only reach the kW range. Aluminum thermite detonation is a typical high-temperature explosive energy release process with high energy density, which can achieve effective energy coupling under limited constraints and possesses reliable detonation capability. Oxygen pressure is the dominant factor affecting detonation velocity and intensity, followed by steel pipe wall thickness. The impact of detonation methods is relatively weak, and detonation methods are substitutable. Synergistic enhancement of detonation velocity and intensity can be achieved by increasing oxygen pressure and optimizing constraint conditions. The excitation efficiency of the three detonation methods follows a consistent ranking: S0 is the strongest, S2 ranks second, and S1 is slightly weaker. This ranking has been verified by shock wave parameters and crater volumes (0.33 m3, 0.24 m3, 0.21 m3). This research provides experimental support for the optimization and application of biomass blasting technology.
, Available online , doi: 10.11883/bzycj-2025-0308
Abstract:
To investigate the complex wake flow characteristics of multi-motor parallel rocket sleds, this study focuses on the mechanisms by which the nozzle horizontal spacing and the impingement height influence the flow structure and ground effect. A three-dimensional physical model was constructed for a dual-rail rocket sled system featuring three nozzles arranged in a triangular pyramid configuration. Four operating conditions were established, including large spacing (l = 7d), small spacing (l = 1d), low impingement height (h = 2d), and high impingement height (h=5.5d). The effects of nozzle center distance and impingement height on the flow field structure and ground effect were comparatively analyzed. Numerical simulations were performed using a computational fluid dynamics (CFD) method based on the Reynolds-Averaged Navier-Stokes equations, coupled with the Realizable k-ε turbulence model for transient solutions. The combustion chamber pressure-time curve derived from interior ballistic theory was applied to the nozzle inlet via a user-defined function (UDF). The sled velocity-time curve, determined from the exterior ballistic particle trajectory equation, was assigned as the far-field pressure boundary condition, enabling a coupled simulation framework of interior ballistics, exterior ballistics and flow field. The computational domain utilized a structured grid with refinement in the jet interaction region and near the ground to ensure calculation accuracy. The velocity and pressure fields obtained from numerical simulations were compared and validated against jet morphology, impingement height, and vortex core positions recorded by high-speed photography (2000 Hz). The flow field structure, pressure distribution, and thermal erosion behavior on the ground under different configurations are systematically revealed. The results indicate that the small-spacing nozzle arrangement triggers intense jet interference without ground effect participation, leading to a multi-peak and slow-recovery pressure evolution feature and substantially delays the flow field relaxation process. The coupling relationship between ground effect and jet interference is dominated by impingement height. At low impingement height conditions, the jet impinging on the ground induces intense reorganization and fragmentation of vortex structures, generating wall jets with velocities up to 960 m/s. Consequently, the peak ground surface temperatures reaches 1286.6 K with sustained high temperatures, which significantly elevates the risk of rail ablation. Conversely, a high impingement height effectively suppresses the ground effect, resulting in a more homogeneous and stable flow field structure. In this case, the peak ground temperature reduced by approximately 65% and maximum velocity reduced by 58%, significantly mitigating ablation risk. The initial phase (0－8 m) of the rocket sled is identified as the critical region subjected to the most severe thermomechanical loads. During this stage, the average acceleration reaches 832.7 m/s2, and the specific action time per unit distance is prolonged to 1.84 ms/m. Coupled with the transient complex flow field, this constitutes an extremely high risk for rail ablation. The numerical simulation results show excellent agreement with high-speed photographic experimental data regarding flow field morphology, impingement height, and vortex core positions, thereby validating the reliability of the established coupled model. This study elucidates the complex flow mechanisms of multi-nozzle parallel systems under strongly constrained conditions, and provids important theoretical foundations and design parameters for structural layout optimization and thermal protection design in high-acceleration, heavy-load rocket sled test systems.
To investigate the complex wake flow characteristics of multi-motor parallel rocket sleds, this study focuses on the mechanisms by which the nozzle horizontal spacing and the impingement height influence the flow structure and ground effect. A three-dimensional physical model was constructed for a dual-rail rocket sled system featuring three nozzles arranged in a triangular pyramid configuration. Four operating conditions were established, including large spacing (l = 7d), small spacing (l = 1d), low impingement height (h = 2d), and high impingement height (h=5.5d). The effects of nozzle center distance and impingement height on the flow field structure and ground effect were comparatively analyzed. Numerical simulations were performed using a computational fluid dynamics (CFD) method based on the Reynolds-Averaged Navier-Stokes equations, coupled with the Realizable k-ε turbulence model for transient solutions. The combustion chamber pressure-time curve derived from interior ballistic theory was applied to the nozzle inlet via a user-defined function (UDF). The sled velocity-time curve, determined from the exterior ballistic particle trajectory equation, was assigned as the far-field pressure boundary condition, enabling a coupled simulation framework of interior ballistics, exterior ballistics and flow field. The computational domain utilized a structured grid with refinement in the jet interaction region and near the ground to ensure calculation accuracy. The velocity and pressure fields obtained from numerical simulations were compared and validated against jet morphology, impingement height, and vortex core positions recorded by high-speed photography (
, Available online , doi: 10.11883/bzycj-2025-0401
Abstract:
Coupling coefficients and stemming coefficients are essential for predicting ground shock magnitude and damage zones from underground explosions, yet their variation with burial depth and dependence on media and explosion types remain insufficiently compared. Experimental and theoretical methods for determining these coefficients in concrete, rock, and soil under both chemical and nuclear explosions were systematically reviewed. Based on collected data from published tests, definitions of different coupling coefficients (energy coupling, ground shock parameter coupling, equivalent yield, and stemming coefficients) were clarified, and conversion relationships among them were derived using energy conservation and wave attenuation principles. The critical burial depth for full coupling was analyzed separately for peak quantities (e.g., peak stress, particle velocity) and integral quantities (e.g., impulse, surface vibration). Experimental results show that peak quantities reach full coupling at a depth approximately equal to the contained explosion cavity radius, which is much shallower than the critical depth for contained explosions. In contrast, impulse and ground motion require depths close to the contained explosion threshold. The coupling energy is proportional to the volume of the damaged zone, and the evolution of coupling coefficients with scaled depth of burial follows a Boltzmann function, requiring only two parameters to define the entire curve. For chemical explosions, the fully coupled equivalent yield coefficient (relative to contained explosions) in concrete, limestone, and soils was quantified, with stemming coefficients ranging from 1.4~1.8 in rocks and 2.4~2.6 in soils. It is concluded that different ground shock parameters exhibit distinct coupling behaviors, challenging the assumption of a universal coupling coefficient. Proposed empirical formulas provide conservative estimates for engineering design, while the underlying mechanisms of free-surface unloading require further quantitative investigation.
Coupling coefficients and stemming coefficients are essential for predicting ground shock magnitude and damage zones from underground explosions, yet their variation with burial depth and dependence on media and explosion types remain insufficiently compared. Experimental and theoretical methods for determining these coefficients in concrete, rock, and soil under both chemical and nuclear explosions were systematically reviewed. Based on collected data from published tests, definitions of different coupling coefficients (energy coupling, ground shock parameter coupling, equivalent yield, and stemming coefficients) were clarified, and conversion relationships among them were derived using energy conservation and wave attenuation principles. The critical burial depth for full coupling was analyzed separately for peak quantities (e.g., peak stress, particle velocity) and integral quantities (e.g., impulse, surface vibration). Experimental results show that peak quantities reach full coupling at a depth approximately equal to the contained explosion cavity radius, which is much shallower than the critical depth for contained explosions. In contrast, impulse and ground motion require depths close to the contained explosion threshold. The coupling energy is proportional to the volume of the damaged zone, and the evolution of coupling coefficients with scaled depth of burial follows a Boltzmann function, requiring only two parameters to define the entire curve. For chemical explosions, the fully coupled equivalent yield coefficient (relative to contained explosions) in concrete, limestone, and soils was quantified, with stemming coefficients ranging from 1.4~1.8 in rocks and 2.4~2.6 in soils. It is concluded that different ground shock parameters exhibit distinct coupling behaviors, challenging the assumption of a universal coupling coefficient. Proposed empirical formulas provide conservative estimates for engineering design, while the underlying mechanisms of free-surface unloading require further quantitative investigation.
, Available online , doi: 10.11883/bzycj-2025-0140
Abstract:
Hydrogen energy, as a zero-carbon energy source, holds broad application prospects in critical defense systems because of its high energy density and zero carbon emissions. To enhance energy utilization efficiency and ensure operational safety, an integrated approach combining experimental and numerical simulations was adopted to systematically examine the effects of hydrogen concentration on explosion dynamics in a confined space. Experiments were carried out in a cylindrical chamber equipped with high-frequency pressure sensors and a high-speed camera to record transient overpressure and track flame propagation behavior. Complementing the experimental setup, computational fluid dynamics (CFD) simulations were implemented using a detailed 19-step hydrogen/air chemical reaction mechanism to accurately reproduce the spatiotemporal evolution of flow field velocity during the premixed gas explosion process. Results indicate that the maximum explosion pressure occurred at a hydrogen volume fraction of 30%, peaking at 0.623 94 MPa. The peak flame area was largest at both 30% and 45%, exceeding results at 15% and 60% by 14.6% and 6.3%, respectively. The 30 % condition also achieved the peak flame area in the shortest time, at 8.2 ms. Furthermore, geometric constraints at the junction of the cylindrical sidewall and the endwall led to accumulation of unburned hydrogen, causing localized increases in density and pressure and resulting in four clearly discernible high-velocity regions within the flow field. At 9 ms, the flow velocity profile along the centerline exhibited symmetry with a dual-peak structure appearing unilaterally. While the 45% condition showed an early transient velocity advantage due to intensified local heat release, the 30% condition demonstrated superior late-stage velocity recovery owing to more stable and sustained combustion near the stoichiometric ratio. These findings underscore the high combustion efficiency and stability achievable near stoichiometric conditions, providing a scientific foundation for the design and optimization of high-efficiency hydrogen combustion systems..
Hydrogen energy, as a zero-carbon energy source, holds broad application prospects in critical defense systems because of its high energy density and zero carbon emissions. To enhance energy utilization efficiency and ensure operational safety, an integrated approach combining experimental and numerical simulations was adopted to systematically examine the effects of hydrogen concentration on explosion dynamics in a confined space. Experiments were carried out in a cylindrical chamber equipped with high-frequency pressure sensors and a high-speed camera to record transient overpressure and track flame propagation behavior. Complementing the experimental setup, computational fluid dynamics (CFD) simulations were implemented using a detailed 19-step hydrogen/air chemical reaction mechanism to accurately reproduce the spatiotemporal evolution of flow field velocity during the premixed gas explosion process. Results indicate that the maximum explosion pressure occurred at a hydrogen volume fraction of 30%, peaking at 0.623 94 MPa. The peak flame area was largest at both 30% and 45%, exceeding results at 15% and 60% by 14.6% and 6.3%, respectively. The 30 % condition also achieved the peak flame area in the shortest time, at 8.2 ms. Furthermore, geometric constraints at the junction of the cylindrical sidewall and the endwall led to accumulation of unburned hydrogen, causing localized increases in density and pressure and resulting in four clearly discernible high-velocity regions within the flow field. At 9 ms, the flow velocity profile along the centerline exhibited symmetry with a dual-peak structure appearing unilaterally. While the 45% condition showed an early transient velocity advantage due to intensified local heat release, the 30% condition demonstrated superior late-stage velocity recovery owing to more stable and sustained combustion near the stoichiometric ratio. These findings underscore the high combustion efficiency and stability achievable near stoichiometric conditions, providing a scientific foundation for the design and optimization of high-efficiency hydrogen combustion systems..
