Blasting excavation of rock masses under high in-situ stress often encounters difficulties in rock fragmentation and a high boulder rate. To gain a deeper understanding of this issue, the stress distribution of rock masses under dynamic and static loads was first studied through theoretical analysis. Then, the ANSYS/LS-DYNA software was employed to simulate the blasting crack propagation in rock masses under various in-situ stress conditions. The fractal dimension was introduced to quantitatively analyze the influence of in-situ stress on the distribution of blasting cracks. The results indicate that in-situ stress primarily affects crack propagation in the later stages of the explosion, while crack initiation and propagation in the early stages are mainly driven by the explosion load. In-situ stress significantly influences the damage area and fractal dimension of cut blasting. Under hydrostatic in-situ stress, as the in-situ stress increases, the damage area and fractal dimension of blasting cracks gradually decrease. Under non-hydrostatic in-situ stress, when the principal stress difference is small, in-situ stress promotes the damage area and fractal dimension of the surrounding rock, enhancing rock fragmentation. However, when the principal stress difference is large, in-situ stress inhibits the damage area and fractal dimension of the surrounding rock, hindering effective rock breaking.

The dynamic characteristics of the filling body are the key parameters for designing the filling ratio and evaluating the stability of an underground stope. The different environment (water-bearing state) of the filling body in the underground stope exerts a complex impact on the mechanical behavior of the filling body. Therefore, six groups of cemented filling body specimens with different states were formed and subjected to dynamic uniaxial impact tests. The effects of water content on the mechanical properties, fractal dimension, and deformation damage characteristics of the cemented backfill under dynamic load were analyzed in depth, and a dynamic damage constitutive model that considers water damage and the compaction stage was established. The results indicate the following: (1) Due to the change of the specimen from the dry state to the water saturation state, the dynamic compressive strength of the cemented filling body decreases from 5.03 Mpa to 1.79 Mpa; however, the ductility of the specimen generally increases, and the filling body specimens with different water contents mainly exhibit tensile failure. (2) There is a significant nonlinear relationship between the water content and the fractal dimension D_b of the cemented backfill specimen, and the growth rate of the fractal dimension D_b tends to slow down with the increase in the water content. (3) From the energy evolution perspective, the water content of the specimen exerts a significant effect on the elastic deformation and failure stage of the stress–strain curve, and the slope of the dissipated energy-strain curve decreases with the increase in water content. (4) Based on the Weibull distribution and damage theory, a statistical damage constitutive model of cemented backfill was established, and it was compared with the experimental curve to verify the rationality of the model. Therefore, the relationship between stress and damage and the strain curves is discussed, and it is inferred that the damage evolution curve of cemented backfill is a typical S-shaped curve that exhibits a stable development-rapid increase-tending to be gentle. This study can provide a theoretical reference for further understanding the dynamic behavior and stability of backfill under different water conditions.

To precisely control the tunnel smooth blasting effect, this study conducts both model experiments and numerical simulations to investigate the impact of shaped charge jet initiation on emulsion explosives and surrounding rock damage fractal characteristics under different ratios of the main-to-secondary charge lengths (L₁/L₂). The study also includes field validation. The results indicate the following: (1) The Arbitrary Lagrangian–Eulerian (ALE) method can accurately reproduce the formation, motion, impact, initiation, and dynamic damage evolution of a shaped charge jet inside a blast hole, with a deviation of less than 6.4% compared to high-speed photography observations. (2) Under the working conditions in this study, when an axial aluminum energetic liner and two-stage air-segmented charge in the peripheral holes are used, the fractal dimension (D_f) initially increases from 1.57 to 1.66 and then decreases to 1.41 as the L₁/L₂ ratio increases. (3) Field test results demonstrate that, when using a two-segment explosive charge with a 20 cm gap between segments and an L₁/L₂ ratio of 2, the average over- or under-excavation is controlled within 7 cm, with the maximum deviation not exceeding 12 cm. The corresponding average fragment size (d₅₀) is minimized, resulting in an excellent smooth blasting effect and effectively controlling the fragmentation of the smooth blasting layer. The conclusions of this study provide valuable insights for the development of advanced shaped charge blasting techniques.

The propagation of stress waves and the evolution of blasting-induced cracks in discontinuous rock mass are strongly influenced by joint characteristics. This study focuses on intersecting-joint rock mass and employs the explicit dynamic finite element method with the RHT constitutive model to simulate blasting. The effects of joint thickness and angle on stress wave behavior and crack propagation were systematically analyzed, and fractal theory was applied using MATLAB 2023b to quantitatively characterize crack distribution. The results indicate that joint thickness significantly affects stress wave transmission and reflection: thicker joints enhance reflected wave intensity while attenuating transmitted waves due to successive reflections within the infill. Joint angle controls the number of reflections between joints, with smaller angles producing stronger interactions. Thickness and angle jointly govern crack propagation, causing cracks to concentrate between the borehole and joints. Thickness primarily determines crack density on the blast-facing side of horizontal joints, whereas joint angle exerts deflecting and inhibiting effects. Fractal analysis shows that the fractal dimension of cracks varies nonlinearly with joint characteristics, with higher values in the borehole-to-joint region, increasing with thickness and angle. These findings provide quantitative insight into stress wave–joint interactions and guidance for optimizing blasting design in jointed rock mass, helping to minimize overbreak and underbreak and improve excavation safety.

The eccentric non-coupled explosive loading structure is widely present in engineering rock blasting. The mechanism of its influence on the blasting effect and energy utilization still needs in-depth research. This paper establishes a theoretical model of peak pressure on the borehole wall for the eccentric non-coupled explosive loading structure and uses the ANSYS/LS-DYNA numerical simulation method to study the influence law of different radial non-coupling coefficients K on the damage characteristics of the rock mass. It also quantitatively analyzes the damage spatial distribution characteristics by combining the fractal theory. Finally, the theoretical and numerical analysis results are verified through field tests. The research results show: (1) The pressure distribution on the borehole wall of the eccentric loading structure is significantly asymmetric, and the pressure difference between the coupled side and the non-coupled side increases with the increase in K value; the error between theoretical calculations and numerical simulation results is less than 10%. (2) As the K value increases from 1.6 to 4.0, the damage areas and crack numbers of both loading structures decrease, but the eccentric loading structure always maintains a larger damage range on the coupled side. (3) Based on the fractal theory analysis, when K = 3.4–4.0, the rock mass damage spatial distribution is the most uniform, the standard deviation is the smallest, and the energy utilization efficiency is the highest. (4) The field test verification shows that when K = 3.75, the pre-splitting effect is the best, the half-hole rate reaches 85.7%, the wall surface flatness is controlled within ±5 cm, and compared with traditional methods, the cost is saved by 7% and the efficiency is increased by 45%. The research results provide an important theoretical basis and reference for the engineering application of the eccentric non-coupled loading structure.