Calculation of hole spacing and surrounding rock damage analysis under the action of in situ stress and joints

Dec 19, 2023

Scientific Reports volume 12, Article number: 22331 (2022) Cite this article

864 Accesses

1 Altmetric

Metrics details

In situ stress and joints have a significant impact on the propagation and attenuation pattern of blast stress waves, and they are two important factors that must be considered for tunnel blasting hole network deployment. This paper proposes a blast stress wave attenuation equation and a peripheral hole distance calculation method under the combined action of in situ stress and joints. First, the static and dynamic parameters of the jointed slate are obtained by drilling core samples in the field and conducting indoor tests. Next, considering the geometric and physical attenuation of the blast stress wave, the attenuation formula of the blast stress wave under the combined action of in situ stress and joints is derived. Based on the theory of the combined action of stress waves and explosive gas, a formula for calculating the peripheral hole distance that integrates the effects of in situ stress, joints, and tensile strength of the rock body is proposed. Finally, LS-PREPOST software is used to analyze the damage to the surrounding rock, verified by an on-site blasting test. The results show that the blast stress wave attenuation formula proposed in this paper can accurately predict the stress wave peak value under the combined action of in situ stress and joints. Combining the geological conditions and blasting parameters of the Bayueshan Tunnel study section, the optimal peripheral hole spacing is calculated to be 45 cm. The average over-excavation value of the grade IV surrounding rock is controlled within 22 cm and the over-consumption of concrete per linear meter is controlled within 100% using the peripheral hole layout method and the hole network layout parameters proposed in this paper. The research results provide a reference for the control of over-excavation and under-excavation in large-section tunnel blasting.

The drilling and blasting method has the advantages of economic efficiency and is currently the main method for the geotechnical excavation of highway and railway tunnels. However, in the tunnel drilling and blasting method construction process, the problem of over-excavation and under-excavation is always a top-priority scientific problem that needs to be solved. Over-excavation and under-excavation not only affect the quality of tunnel blasting and increase tunnel support structure requirements1,2 but also raise operating costs and reduce construction progress. Therefore, it is of great importance for safe and efficient tunnel blasting to propose effective measures to control over-excavation and under-excavation.

A reasonable arrangement of peripheral holes is an effective means to control the over-excavation and under-excavation in tunnel blasting. Wang et al.3 established a numerical analysis model of tunnel smooth blasting using LS-DYNA software. By comparing the results of numerical simulation and field monitoring, it was shown that the spacing of peripheral holes should be mainly considered in the massive structure and the fractured structure. Paul and Peter4 analyzed the causes, hazards, and control measures of tunnel over-excavation and under-excavation caused by the drilling and blasting method through a model test system, and proposed a reasonable peripheral hole design scheme. Adel et al.5 established a tunnel over-excavation prediction model considering factors such as tunnel cross-sectional area and peripheral hole parameters. Yin et al.6 proposed a drilling method combining long and short holes in the peripheral holes, which effectively solved the problem of over-excavation and under-excavation. Liu et al.7 indicated that the peripheral hole parameters should be adjusted in time according to changes in the surrounding rock conditions. Huang et al.8 studied the application of a double-layer peripheral hole retaining wall blasting design in the blasting engineering of the soft and broken surrounding rock, and optimized the peripheral hole layout and charge structure. Based on the principle of smooth blasting, Man and Liu9 studied the blast hole arrangement scheme of peripheral holes in the case of empty hole spacing. Man et al.10 focused on the effect of peripheral hole spacing on the effect of smooth blasting.

However, a rock mass in nature exhibits anisotropy due to the existence of structural planes such as joints and fissures11 (1988), when the blast stress wave propagates to these structural surfaces, multiple reflections occur, resulting in the attenuation of the blast stress wave. Hyongdoo et al.12 proposed the over-excavation resistance coefficient (ORF) based on the relationship between the over-excavation measured values and geological parameters, indicating that the discontinuity factor has the greatest impact on blasting and over-excavation. Zhao and Cai13, Cai et al.14, and Zhao et al.15,16 studied the propagation and attenuation pattern of stress waves after vertically traversing groups of parallel joints. Li et al.17,18,19 derived the propagation equation of multiple parallel joints with the oblique incidence of stress waves. Perino et al.20 derived the propagation and attenuation pattern of stress waves in cross joints. Chai et al.21 derived the propagation equation of plane P-waves in a rock mass with two intersecting joints.

At the same time, the rock surrounding a tunnel is already in a certain initial stress state before it is subjected to a blast load. When the initial stress state changes, the micropores inside the rock mass evolve, affecting the propagation of blast stress waves. Mandal and Singh22 believed that the initial stress field of the surrounding rock had a great influence on the blast load. Mandal et al.23 indicated that the influence of initial in situ stress and stress redistribution should be considered during tunnel excavation, and based on this, an empirical formula for stress wave attenuation in media was proposed. Zheng 24 carried out a three-dimensional high in situ stress rock mass ultra-deep hole blasting test, and concluded that the initial in situ stress has an inhibitory effect on the expansion of blast cracks. Zhang et al.25 studied the damage to surrounding rock under blast load, showing that in situ stress has an inhibitory effect on the blast tension effect. Li et al.26 studied the propagation and attenuation pattern of elastic waves in a deeply fractured rock mass through model tests. Li et al.27 examined the propagation pattern of stress waves with the vertical incidence of multiple parallel jointed rock masses under different in situ stresses based on indoor model tests. Fan et al.28 analyzed the propagation pattern of a one-dimensional stress wave vertically incident on a single joint under the action of in situ stress. Liu and Xu29 established the differential equation of the motion of rock mass under the action of an explosion load and under the action of in situ stress, and completed the numerical simulation of a rock mass explosion under the action of in situ stress through the implicit–explicit continuous solution. Fan et al.30 studied the propagation pattern of stress waves in the rock medium under different initial stresses based on a numerical analysis method. Zhang et al.31 researched the propagation pattern of stress waves and the action mechanism of joint surfaces under the combined action of in situ stress and joints.

At present, most of propagation theories about stress waves in jointed rock masses treat the stress wave as a plane wave, and research on the attenuation pattern of cylindrical waves under the combined action of in situ stress and joints is rarely reported. In this paper, we took the Tongliang–Anyue Expressway Bayueshan Tunnel as the research background, drilled rock samples in the field, and conducted laboratory tests to obtain the static and dynamic parameters of the jointed slate. Considering the geometric and physical attenuation of the blast stress wave at the same time, the attenuation formula of the blast stress wave under the combined action of in situ stress and joints is proposed. Based on the combined action of stress wave and explosive gas theory, a calculation formula for peripheral hole spacing is proposed that comprehensively considers the effects of in situ stress, joints, and rock mass tensile strength. LS-PREPOST software is used to establish a single-hole blasting numerical analysis model of an in situ stress-jointed rock mass, and to verify the correctness of the explosion stress wave attenuation formula. We established numerical analysis models for tunnel blasting with peripheral hole spacings of 40 cm, 45 cm, 50 cm, and 55 cm, respectively, and compared the damage depth and peak vibration velocity (PPV) of the retained rock mass to verify the rationality of the calculation formula for the peripheral hole spacing. We conducted field blasting tests based on the theoretical derivation and numerical analysis results to provide a reference for the control of over-excavation and under-excavation in large-section tunnel blasting.

The research takes the ZK13 + 760 ~ ZK13 + 940 section of the Tongliang–Anyue Expressway Bayueshan Tunnel as the engineering background. The total length of the research section is 180 m, the average burial depth is 56 m, the rock mass integrity coefficient Kv is 0.68, and the revised basic quality index of the surrounding rock (BQ) is 321, which is classified as grade IV. The upper layer of the tunnel site is covered with residual slope deposits (Qel+dl), containing gravelly powder clay, with 35–40% gravel content and 0.5–2 cm particle size. The underlying bedrock of the tunnel site is tuffaceous slate of the second section (Ptbnbf2) of the Fanzhao Formation of the upper Banxi Group in the Yuangu realm. According to the characteristics of rock joints and fissure development, hardness and integrity, it is divided into two layers: a strong weathering layer (4.5–24.5 m) and a medium weathering layer. Figure 1 shows the rock surrounding the tunnel face in the study section. The lithology of the surrounding rock is mainly slate, and there are many joints and fissures at 60° to the horizontal plane. To obtain the static and dynamic parameters of slate under the combined action of vertical load σV and horizontal load σH, we drilled core rock samples on site and made standard slate samples with 60° joint fissures, and we carried out uniaxial compression, triaxial compression, Brazilian splitting, and impact dynamics tests.

Engineering background. (a) surrounding rock condition of tunnel face; (b) standard specimen for static test; (c) standard specimen for dynamic test; (d) Tong’an Expressway Bayueshan Tunnel.

We measured the diameter and the height of the rock samples with vernier calipers, and weighed the rock samples with an electronic balance. The densities of the nine rock samples were 2746 kg/m3, 2747 kg/m3, 2744 kg/m3, 2749 kg/m3, 2754 kg/m3, 2778 kg/m3, 2735 kg/m3, 2758 kg/m3, and 2752 kg/m3, with an average density of 2752 kg/m3. The uniaxial compression test, triaxial compression test, and Brazilian splitting test were carried out on the microcomputer-controlled electro-hydraulic servo rock triaxial testing machine TAJW-2000, as shown in Fig. 2.

Static test. (a) Uniaxial compression test; (b) triaxial compression test; (c) Brazilian splitting test.

In the uniaxial compression test, a preload of 0.2 KN was first applied and then converted to displacement loading with a loading rate of 0.12 mm/min after the instrument was stabilized. The load was continuously applied until the specimen was damaged, and followed by immediate unloading. In the uniaxial compression test process, shown in Fig. 2a, the displacement change process of the rock specimen was collected by an axial displacement transducer and a circumferential displacement transducer, with the stress–strain curve shown in Fig. 3. The uniaxial compressive strengths of the three rock samples were 54.2 MPa, 50.8 MPa, and 37.2 MPa, respectively, and the average uniaxial compressive strength fc was 47.4 MPa. In Fig. 1, the slate is subjected to the combined action of the vertical load σV and the horizontal load σH, which is two-way loading. The literature Huang32 shows that the ratio β of bidirectional compressive strength fbc to uniaxial compressive strength fc is 1.493fc– 0.0634. Therefore, the bidirectional compressive strength fbc of slate is 55.4 MPa. The elastic moduli of the three rock samples were 40.2 GPa, 38.6 GPa, and 35.3 GPa, respectively, with an average elastic modulus of 38.03 GPa. The Poisson's ratios of the three rock samples were 0.23, 0.28, and 0.29, respectively, with an average Poisson's ratio of 0.27.

Uniaxial compressive stress–strain curve.

In the Brazilian splitting test, a linear load was applied in the diameter direction of the cylindrical rock specimen with a loading rate of 0.15 mm/min until the specimen was damaged. The Brazilian splitting test process is shown in Fig. 2c. The tensile strength \(\sigma_{t}\) of the rock sample can be calculated as follows:

where P is the load corresponding to the splitting failure of the rock sample, and D and L represent the diameter and height of the rock sample, respectively. The loads of the three rock samples at failure were 43.5 kN, 41.6 kN, and 56.5 kN, respectively. According to Eq. (1), the tensile strength of the three rock samples were 11.1 MPa, 10.6 MPa, and 14.4 MPa, respectively, and the average tensile strength of the rock samples was 12.03 MPa.

The slate impact dynamic test was carried out using the separate Hopkinson pressure bar test system ALT100 developed by Archimedes Industrial Technology Co., Ltd., as shown in Fig. 4. The diameters of the impact bar, incident bar, and transmission bar were all 50 mm, and their lengths were 400, 2000, and 2000 mm, respectively. The density of the pressure bar was 7.81 g/cm3, the elastic modulus was 210 GPa, the Poisson's ratio was 0.28, and the longitudinal wave velocity was 5410 m/s. The data acquisition system comprised strain gauges pasted on the pressure bar, a wheatstone bridge (strain gauge junction bridge box), super dynamic strain gauges, and a high-speed acquisition system.

SHPB test system.

Based on the assumption of one-dimensional stress waves and dynamic balance, the dynamic stress \(\sigma \left( t \right)\) and dynamic strain \(\varepsilon \left( t \right)\) can be accurately calculated using the "two-wave method", with the incident strain signal \(\varepsilon_{I} \left( t \right)\), the reflected strain signal \(\varepsilon_{R} \left( t \right)\), and the transmitted strain signal \(\varepsilon_{T} \left( t \right)\) measured on the pressure bar. The calculation method is as follows:

where A and L are the cross-sectional area and length of the sample, respectively; E and A0 represent the elastic modulus and cross-sectional area of the pressure bar, respectively; and C denotes the longitudinal wave velocity of the pressure bar.

Impact compression tests were conducted on the slate samples under the impact pressure of 0.15 MPa, 0.2 MPa, and 0.3 MPa, respectively, and the impact compression stress–strain curve was drawn through Eq. (2), as shown in Fig. 5. The average peak stress of slate under 0.15 MPa, 0.2 MPa, and 0.3 MPa impact pressures was 76.8 MPa, 95.7 MPa, and 112.2 MPa, respectively. The dynamic tensile strength of rock has little sensitivity to the loading strain rate, and it can be considered that the dynamic tensile strength of rock is the same as the uniaxial tensile strength. The mechanical parameters of the 60° jointed slate are shown in Table 1.

Impact compressive stress–strain curves: (a)–(c) are impact pressures of 0.15 MPa, 0.2 MPa, and 0.3 MPa, respectively.

In the area near the blasting, the stress wave generated by the blasting of the blasthole usually propagates outward in the form of a cylindrical stress wave (two-dimensional wave), and its wavefront series of cylindrical surfaces take the central axis of the blasthole as the axis. Cylindrical waves in intact rocks undergo both geometric and physical attenuation with increasing propagation distance33, where geometric attenuation is the attenuation due to the increase in the spatial distribution of stress wave energy. The stress on the wavefront of the cylindrical wave and the vibration velocity of the particle are both attenuated by \(1/\sqrt r\)34, where r is the distance between the wavefront and the wave source. The particle vibration velocity vr on the wavefront distanced from the wave source r is

where r0 and v0 are the radius of the initial wavefront of the cylindrical wave and the vibration velocity of the particle, respectively. The physical attenuation of cylindrical waves in intact rocks is caused by the friction effect on the surface of microcracks inside the propagating medium35, and the particle vibration velocity on the cylindrical wave wavefront decays negatively exponentially with the increase in propagation distance:

where \(\alpha\) is the physical attenuation coefficient of the cylindrical wave.

In summary, considering both physical attenuation and geometric attenuation, the propagation attenuation equation of cylindrical waves in intact rock is

According to the stress wave propagation theory, radial stress \(\sigma_{r}\) at a certain point caused by the propagation of the cylindrical stress waves in the rock mass is

It can be seen from the elastic mechanics that, without considering the tangential strain \(\varepsilon_{\theta }\), it can be considered that the radial stress \(\sigma_{r}\) and the radial strain \(\varepsilon_{r}\) at a certain point in the rock body approximately satisfy the following relationship:

Combining Eqs. (6) and (7), we can obtain

where Ed is the dynamic elastic modulus.

Assuming that the density, radial wave velocity, and dynamic elastic mode of the model material remain unchanged during the blasting process, the radial strain on the cylindrical wave front is approximately proportional to the particle vibration velocity. Therefore, it can be considered that the attenuation coefficient of the radial strain of the particle and the radial vibration velocity are approximately the same:

Dong produced test blocks with single and double joints, and carried out model tests36. The schematic diagram of the model is shown in Fig. 6.

Schematic diagram of the model36: (a) single-joint model; (b) double-joint model.

Through analysis of the Nos. 13–16 measuring point strain peaks under different confining pressures, the relationship between the physical attenuation coefficient \(\alpha\) of the cylindrical wave in the complete rock mass and the confining pressure \(\sigma\) is obtained:

Substituting Eq. (10) into Eq. (9), the attenuation pattern of the particle point radial strain when the cylindrical wave propagates in the complete rock mass under different confining pressures is obtained:

Substituting Eq. (11) into Eq. (7), the attenuation pattern of the particle radial stress when the cylindrical wave propagates in the complete rock mass under different confining pressures is obtained:

Equation (12) is the attenuation pattern of the particle radial stress when the cylindrical wave propagates in the complete rock mass, without considering the influence of the joint on the attenuation of the radial stress. However, the rock mass is often rich in joint planes, which have a significant impact on the attenuation of radial stress at the particle point. Therefore, it is necessary to establish the attenuation relationship of the radial stress of the particle when the cylindrical wave propagates in the jointed rock mass under different confining pressures:

where K is the correction coefficient of the joint, which is related to the confining pressure and the property of the joint, and its value is between 0 and 1.

Dong36 set the measurement points 17–24 as the stress measurement points in the model test. The ratio of the peak values of the time history curves recorded at the stress measurement points before and after the joints with different angles is defined as the transmission coefficient of the cylindrical wave, and the variation pattern of the cylindrical wave transmission coefficient of the single-joint and double-joint test blocks under different confining pressures is obtained, as shown in Fig. 7.

Cylindrical wave transmission coefficients under different confining pressures: (a) single-joint test block; (b) double-joint test block.

The joint correction coefficient K is defined as the ratio of the peak value of the stress time history curve of measuring points 11, 18, and 19 to the peak value of the stress calculated by Eq. (12). Equation (12) is used to calculate the peak stress at measuring points 17–24, and the joint correction coefficient K can be calculated by combining with the peak data of the stress time-history curve at measuring points 17–24, as shown in Table 2. The relationship between the joint correction coefficient K and the confining pressure of the single-joint test block and the double-joint test block is plotted in Fig. 8.

Variation of joint correction coefficient K with confining pressure: (a) single-joint test block; (b) double-joint test block.

It can be seen from Fig. 8 that the joint correction coefficient K has a good correlation with the confining pressure, the joint inclination angle, and the number of joints. When there is no confining pressure, the joint correction coefficient decreases with the increase in the joint inclination angle. When the confining pressure increases from 0 to 1.5 MPa, the joints are gradually closed under the action of the confining pressure, the stiffness of the joints increases rapidly, and the attenuation of the stress wave decreases, so the joint correction coefficient K increases. When the confining pressure continues to increase to 3.0 MPa, micro-cracks are formed in the closed joints, and the expansion of new cracks leads to an increase in the attenuation of stress waves, so that the joint correction coefficient K decreases. The joint correction coefficient K in the double-joint test block is lower than that in the single-joint test block. The greater the number of joints, the stronger the transmission and reflection effect of the stress wave between the joints, the attenuation amplitude of the stress wave after passing through the joint group decreases, and the attenuation amplitude of the joint correction coefficient K decreases37. Therefore, this paper only analyzes the single-joint test block and the double-joint test block.

When the uncoupled continuous charge is used, the impact stress of the blast hole wall after blasting is

where \(\rho_{0}\) is the density of the explosive (kg/m3); D represents the detonation velocity of the explosive (m/s); dc denotes the diameter of the cartridge (m); db refers to the diameter of the blast hole (m); lc is the charge length and lb is the length of the blast hole (m); and n depicts the pressure increase coefficient of the detonation gas impacting the hole wall, generally taking 10.

According to Eq. (13), considering the influence of in situ stress and joint properties, the radial stress and tangential stress of the particle at the distance r from the center of the blast hole are

where rb is the radius of the blast hole (m); b represents the lateral pressure coefficient of the rock mass (\(b = \mu_{d} /\left( {1 - \mu_{d} } \right)\)); μd stands for the dynamic Poisson's ratio; and the other symbols are the same as defined above.

The condition of tensile failure of the rock mass under the action of tangential stress \(\sigma_{\theta }\) is \(\sigma_{\theta } \ge \sigma_{td}\), and \(\sigma_{td}\) is the dynamic tensile strength of the rock mass. Therefore, the calculation expression of the fracture circle radius R is

After determining the relevant parameters, the value of the fracture circle radius R can be calculated using the MATLAB program.

According to the theory of the combined action of stress wave and blasting gas, the formation of through cracks between blasting holes is due to the static pressure of the blasting gas, and the formation conditions of through cracks are

where DP is the peripheral hole spacing (m); and Pb represents the pressure when the blast hole is filled with explosive gas (Pa).

Based on the theory of entropy expansion, the pressure Pb when the blast hole is filled with explosive gas is

where Pa is the explosion pressure (Pa); Pk stands for the critical pressure during the expansion of the explosion gas (taking 100 MPa); Vc and Vb represent the volume of the cartridge and the blast hole, respectively (m3); k denotes the adiabatic coefficient of the explosive; and h refers to the isentropic coefficient of the explosive (k = h = 3.0).

From Eqs. (17)–(19), the calculation expression of the peripheral hole spacing under the action of in situ stress and joints can be obtained:

According to the blasting manual, when DP = 0.8 W, the blasting effect is better, and the value of the minimum resistance line W can be calculated.

Based on the original blast stress wave attenuation formula, this paper proposes a new blast stress wave attenuation formula by considering the combined effect of in situ stress and joints. To verify the correctness of Eq. (15), LS-PREPOST software was used to establish a three-dimensional numerical calculation model of the double joint under 1.5 MPa in situ stress, as shown in Fig. 9a–c. In Fig. 9a–c, the in situ stress is 1.5 MPa, and the physical attenuation coefficient α is 0.56224. The dip angles of the joints are 30°, 60°, and 90°, and the joint correction coefficients K are 0.681, 0.711, and 0.719, respectively. To save calculation time, the size of the model is 5 m × 5 m × 2 m. The diameter of the blast hole is 42 mm, the depth of the blast hole is 1.4 m, and the diameter of the cartridge is 32 mm. The first joint is 1 m away from the center of the blast hole, and the two joints are separated by 1 m. The keyword *DEFINE is used to define the load curve CURVE, 0–1.5 MPa is loaded, the keyword *INTERFACE is used to output the DYINA file with in situ stress, the original k file is replaced, and 1.5 MPa of stress is applied on the top and side surfaces of the model. The rock is defined as a solid and the Holmquist–Johnson–Cook (HJC) constitutive model is adopted, with the parameters shown in Table 3. The joint material parameters are shown in Table 4. Explosives and air are defined as fluids, and the parameters are shown in Tables 5 and 6. The fluid is divided into meshes using common nodes, and the solid and the fluid are connected using fluid–solid coupling. Except for the free surface, all other surfaces are set to no reflection boundary conditions. Four monitoring points (A, B, C, and D) were arranged at 0.6 m, 1.2 m, 1.8 m, and 2.4 m from the center of the blast hole to record the stress attenuation pattern during the blasting. The numerical simulation results and the calculation results of Eq. (15) are compared and analyzed, as shown in Fig. 9d –f.

Verification of the explosion stress wave attenuation formula: (a)–(c) the three-dimensional numerical calculation models under different joint inclination angles; (d)–(c) the comparison between the theoretical calculation results and the numerical simulation results.

Under the three joint inclination angles, the theoretical calculation results of the stress wave attenuation are consistent with the numerical simulation results, which verifies the correctness of the explosion stress wave attenuation formula proposed in this paper. With the increase in blast center distance, the radial stress wave peak value gradually decreases. During the propagation from point A to point D, the attenuation amplitude of the peak value of the radial stress wave is 60%. The attenuation amplitude of the radial stress wave passing through the first joint is 30%, and the attenuation amplitude when passing through the second joint is 18%. The explosion stress wave occurs with multiple reflections at the joint, resulting in the attenuation of the peak value of the radial stress wave. The first joint is close to the explosion source, and the attenuation amplitude is large. In the calculation results of Eq. (15), the peak values of radial stress waves at point A under the three joint inclination angles are different, because all the Eq. (15) calculation results have considered the influence of joints.

According to the derivation of the peripheral hole spacing calculation formula under the action of in situ stress and joints in section "Determination of the peripheral hole spacing", combined with the blasting parameters of the Bayueshan Tunnel, the values of the peripheral hole spacing DP and the minimum resistance line W are determined. The average burial depth of the Bayueshan Tunnel study section is 56 m, the slate density is 2752 kg/m3, and σV = σH = ρgh = 1.5 MPa. The physical attenuation coefficient \(\alpha = 0.56224\) is calculated by Eq. (10). The No. 2 rock emulsion explosive is used for the on-site blasting, the density of the explosive is 1.24 g/cm3, the detonation speed is D = 4200 m/s, the diameter of the blast hole is 42 mm, and the diameter of the cartridge is 32 mm. Peripheral holes are equipped with 1.5 cartridges, the charge length lc is 0.45 m, and the blast hole length lb is 1.4 m. Calculated using Eq. (14), the impact stress of the blast hole wall after the explosion is P0 = 177.6 MPa. It can be seen from Table 2 that the joint correction coefficient K of the double-joint model is 0.711 when the confining pressure is 1.5 MPa and the joint inclination angle is 60°. In engineering blasting, it can be considered that the relationship between the dynamic Poisson's ratio μd and the static Poisson's ratio μ is μd = 0.8μ38. The calculated dynamic Poisson's ratio μd is 0.22, and the lateral pressure coefficient of the rock mass is \(b = \mu_{d} /\left( {1 - \mu_{d} } \right) = 0.28\). It can be seen from Table 2 that the dynamic tensile strength of 60° slate is 12.03 MPa. Calculated using Eq. (17), the radius of the crack ring is R = 18.1 cm. Using Eq. (19), Pa is calculated as 3.7774 GPa, and Pb = 24.6 MPa. Equation (20) calculates that the peripheral hole distance DP = 45 cm under the action of in situ stress and joints.

To verify that 45 cm is the optimal value of the peripheral hole spacing around the research section of the Bayueshan Tunnel, the peripheral hole distances DP are set as 40 cm, 45 cm, 50 cm, and 55 cm, respectively, and the corresponding minimum resistance lines W are 50 cm, 56 cm, 62 cm, and 68 cm, respectively. By selecting the material constitutive model and determining the model parameters, the numerical analysis model of Bayueshan Tunnel blasting is established.

60° slate constitutive model and parameter determination.

The HJC (Holmquist–Johnson–Cook) constitutive model is a rate-dependent constitutive model proposed by Holmquist, Johnson, and Cook to solve the problem of the large deformation of concrete under a high strain rate and a high pressure load39, and is widely used in the dynamic analysis of rock materials. The HJC constitutive model totally includes 21 parameters: basic mechanical parameters—R0, fc, G, and T; strength parameters—A, B, C, N, SFMAX, and EPSO; damage parameters—D1, D2, EFMIN, and FS; and pressure parameters—Pc, \(\mu_{c}\), Pl, \(\mu_{l}\), K1, K2, and K3.

The basic mechanical parameter R0 is taken as 2752 kg/m3, fc is taken as the bidirectional compressive strength of 55.4 MPa, \(G = E/2\left( {1 + \mu } \right) = 15.212\;\;{\text{GPa}}\), \(K = E/3(1 - 2\mu ) = 25.35\;\;{\text{GPa}}\), and T takes 12.03 MPa.

Take the 0.15 MPa impact air pressure test result in Fig. 5a as an example. To eliminate the influence of hydrostatic pressure on the dynamic strength of limestone, starting from the characteristic tensile strength \(T^{*} = T/f_{c}\), straight lines are drawn through the data points of equivalent strength under different strain rates, and a straight line is drawn perpendicular to the horizontal axis at the constant characteristic hydrostatic pressure P* = 1/3. The intersection point of the straight lines at different slopes represents the characteristic equivalent stress of 60° jointed slate under different strain rates, as shown in Fig. 10a. By fitting the data points of the characteristic equivalent stress under different strain rates using the straight-line equation, the strain rate influence coefficient C = 0.000624 was obtained, as shown in Fig. 10b.

Determination of the strain rate influence coefficient C. Relationships between: (a) characteristic equivalent stress and characteristic hydrostatic pressure under different strain rates; and (b) characteristic equivalent stress and strain rate.

According to the triaxial compression test results of 60° slate, the cohesion c = 18.4 MPa was calculated. The characterized cohesion strength is A = c/(1 + Cln10−4) fc, and by substituting the value of the strain rate influence coefficient C, the characteristic cohesion strength A = 0.33 can be obtained. According to \(\sigma^{*} = \left( {\sigma_{1} - \sigma_{3} } \right)/f_{c}\) and \(p^{*} = \left( {2\sigma_{1} + \sigma_{3} } \right)/3f_{c}\), \(\sigma^{*} - P^{*}\) curves are drawn, and the fitting can be obtained such that the values of parameters B and N are 0.5765 and 0.413, respectively. SFMAX and EPSO are, respectively, 20.0 and 1.0 according to the literature40.

For the damage parameters, \(D_{1} = 0.01/(1/6 + T^{*} ) = 0.026\) and D2 takes a constant value of 1.0. According to the literature40, EFMIN and FS are taken as 0.01 and 0.004, respectively.

The pressure parameter is Pc = fc/3 = 55.4/3 = 18.47 MPa. μc = Pc/K = 7.3 × 10–4, \(\mu_{l}\) is obtained from \(\mu_{l} = \rho_{g} /\rho_{0} - 1 = 0.054\), where \(\rho_{g}\) is the compacted density, which is 2900 kg/m3. The parameters Pl, K1, K2, and K3 are non-sensitive parameters, and the research results of this paper refer to the literature40 for their values. The parameters of the 60°slate constitutive model are shown in Table 3.

Joint constitutive model and parameter determination.

The No. 003 material model *MAT_PLASTIC_KINEMATI is selected as the material constitutive model of the joint, and its basic parameters are shown in Table 4.

Explosive constitutive model and parameter determination.

The material constitutive model of the explosive material selects the No. 008 material model *MAT_HIGH_EXPLOSIVE_BURN, and the pressure, volume, and energy characteristics of the explosive during the explosion process are simulated by defining the explosive detonation parameters in the explosive state equation *EOS_JWL. The JWL state equation is:

where P is the unit detonation pressure; V represents the initial relative volume; E0 stands for the initial specific internal energy; and A, B, R1, R2, and \(\omega\) are the basic parameters of the material.

The No. 2 rock emulsion explosive is used for blasting the Bayueshan Tunnel, and its basic parameters are shown in Table 5.

Air constitutive model and parameter determination.

The air material constitutive model selects the No. 009 material model *MAT_NULL, and the EOS_LINEAR_POLYNOMIAL state equation is:

where E is the internal energy of the material; and C0, C1, C2, C3, C4, C5, and C6 represent the coefficients of the state equation. The basic parameters are shown in Table 6.

Determination of the constitutive model and parameters of the clay.

The No. 005 material model *MAT_SOIL_AND_FOAM is selected as the material constitutive model of the clay material, and its basic parameters are shown in Table 7.

The research section of the Bayueshan Tunnel adopts the upper and lower bench blasting construction method. In this paper, numerical analysis models of upper bench blasting are established to discuss the influence of peripheral hole spacing on the damage depth and peak vibration velocity of the retained rock mass. LS-PREPOST software is used to establish a three-dimensional numerical multi-joint model, when the in situ stress is 1.5 MPa, the joint angle is 60°, and the physical attenuation coefficient α and the joint correction coefficient K are 0.56224 and 0.711, respectively. The model size is 20 m × 20 m × 2 m, and the mesh is divided by mapping to achieve common nodes between fluids and between solids. The rock mass, joints, explosives, air, and clay use the material models and parameters in Tables 3, 4, 5, 6 and 7. Stress of 1.5 MPa is applied to the top and side surfaces of the model. To avoid the influence of reflected waves caused by artificial boundaries on the calculation results, the tunnel surface is set as free boundary conditions, and the other planes are set as non-reflection boundary conditions. The three-dimensional numerical analysis model is shown in Fig. 11, and the parameters of each blast hole in the model are shown in Table 8.

Three-dimensional numerical analysis model. (a) E = 40 cm; (b) E = 45 cm; (c) E = 50 cm; (d) E = 55 cm.

Figure 12 shows the cumulative damage cloud map of the remaining rock mass after all blasts are completed, with the damage levels decreasing from red to green to blue. The cumulative damage cloud map shows that the damage depth gradually degenerates from the explosion center to the outer free surface. Wang et al.45 showed that the red damage depth line can be used to estimate the over-excavation and under-excavation values after tunnel blasting. Therefore, the red damage depth contour lines under blasting with different peripheral hole distances are drawn, as shown in Fig. 13. It can be seen from Fig. 13 that the remaining rock mass after blasting with different peripheral hole spacings is damaged to different degrees, resulting in different levels of over-excavation and under-excavation. When the peripheral hole spacing is 55 cm and 50 cm, although the damage depth of the remaining rock mass in each part is small, there are different degrees of under-excavation. Under-excavation will affect the steel arch erection process, and in severe cases, additional explosives are required, increasing the operating cost. The main reason for this phenomenon is that the distance between the peripheral holes is too far, the explosion stress wave cannot produce the superposition enhancement effect, and no through cracks are formed between the blast holes. When the peripheral hole spacing is 45 cm, the damage depth of the remaining rock mass at each part is small, and there is no under-excavation. When the peripheral hole spacing is 40 cm, the damage depth of the remaining rock mass in each part is slightly greater than 45 cm, and there is also no under-excavation. This is because the explosion stress wave produces a large superposition enhancement effect between the blast holes, and more residual energy causes more damage to the retained rock mass. When the peripheral hole distance is 40 cm, the drilling workload increases and the operating cost increases. Therefore, from the perspective of retaining the damage depth of the rock mass, the optimal value of peripheral hole spacing is 45 cm, which is consistent with the calculation result of Eq. (20).

Cloud map of accumulated damage of retained rock mass. (a) E = 40 cm; (b) E = 45 cm; (c) E = 50 cm; (d) E = 55 cm.

Damage depth contour. (a) E = 40 cm; (b) E = 45 cm; (c) E = 50 cm; (d) E = 55 cm.

The PPV method is the most commonly used method for the actual measurement of surrounding rock damage in the field46,47,48. Four groups of PPV measuring points were arranged at 0.5 m, 1.0 m, 1.5 m, and 2.0 m away from the tunnel outline, with 41 measuring points in each group evenly distributed around the tunnel outline, as shown in Fig. 14. The PPV curves of different distances are shown in Fig. 15. The average PPV values of each group when the peripheral hole distances are 40 cm, 45 cm, 50 cm, and 55 cm are shown in Table 9. As the distance from the tunnel contour increases, the average PPV value gradually decreases, and the attenuation amplitude is related to the distribution of joints. With the increase in the peripheral hole spacing, the average PPV value of each group gradually decreases. Compared with the average PPV value when the peripheral hole spacing is 40 cm, the average PPV value attenuation amplitude is smaller when the peripheral hole spacing is 45 cm, and the average PPV value attenuation amplitude is larger when the peripheral hole spacing is 50 cm and 55 cm. According to the PPV damage criterion of jointed rock mass49, for hard rock, the threshold is 70 cm/s. When the peripheral hole spacing is 40 cm, 45 cm, 50 cm, and 55 cm, the damage range is 183 cm, 167 cm, 115 cm, and 62 cm, respectively. The greater the peak vibration velocity PPV, the greater the damage depth. Combined with the analysis results in section "Retained rock mass damage depth", it can be seen that when the peripheral hole spacing is 40 cm, the average PPV value is larger, the damage depth is larger, and it is easy to cause greater over-excavation. When the peripheral hole spacing is 50 cm or 55 cm, the average PPV value is small, the damage depth is shallow, and the phenomenon of under-excavation occurs easily. Therefore, considering the peak vibration speed, the optimal value for the peripheral hole distance is 45 cm.

Layout of measuring points.

PPV curve. (a) E = 40 cm; (b) E = 45 cm; (c) E = 50 cm; (d) E = 55 cm.

Based on the theoretical derivation results in section "Determination of the peripheral hole spacing" and the numerical analysis results in section "Numerical simulation", the peripheral hole spacing DP around the study section of the Bayueshan Tunnel is 45 cm, the resistance line W is 56 cm, the hole network layout is shown in Fig. 16, and the hole network layout parameters are shown in Table 7 (the peripheral hole distance is 45 cm).

Layout of hole mesh for field blasting test.

Three field blasting tests were carried out using the above methods and parameters. After the completion of the first blasting, the residual holes in the peripheral holes of each part are clear, as shown in Fig. 17a. The effect of over-excavation and under-excavation is shown in Fig. 17b. Measurements of the over-excavation and under-excavation for each part of the upper step are shown in Fig. 17c, d. It can be seen from the measurement results that, after the blasting of the upper steps is completed, all parts are over-excavated and there is no under-excavation. The average over-excavation is 21.8 cm, the over-excavation per linear meter is 5.82 m3, the design concrete dosage for the upper steps is 6.41 m3 per linear meter, and the concrete over-consumption per linear meter is 90.8%. After the completion of the second blasting, the average over-excavation is 22.1 cm, the over-excavation per linear meter is 5.78 m3, and the concrete over-consumption per linear meter is 90.2%. After the third blasting, the average over-excavation is 21.4 cm, the over-excavation per linear meter is 5.9 m3, and the concrete over-consumption per linear meter is 92%. Under the combined action of in situ stress and joints, the blasting construction is carried out using the peripheral hole layout method and hole network layout parameters proposed in this paper, the concrete excess consumption per linear meter is controlled within 100%, and the concrete excess consumption control effect is good.

Blasting effect. (a) Peripheral hole residual hole; (b) steel arch erection effect; (c) and (d) measurement of over-excavation and under-excavation.

Seven representative locations of the cumulative damage cloud map were selected for zoom-in analysis, as shown in Fig. 12. Take Fig. 12b as an example for analysis. Position 1: When the line connecting the peripheral holes is parallel to the joint surface, the damage degree of a single blast hole increases, because the explosion stress wave causes multiple reflections at the joint surface, forming a superposition of stress waves. Position 2, Position 3, and Position 4: When the peripheral holes are just located on the joint surface, the explosion stress wave first propagates along the joint surface, and further expands and penetrates based on the original joint surface, resulting in serious over-excavation. Position 5, Position 6, and Position 7: When the joint surface is located between two peripheral holes, the explosion stress wave has a large attenuation effect on the joint surface, so the over-excavation phenomenon near these positions is effectively controlled. In summary, the relative position of the joint surface and the peripheral holes will have a great impact on rock mass damage, and on over-excavation and under-excavation50,51, which can be overcome by increasing or decreasing the charge.

Uniaxial compression, triaxial compression, Brazilian splitting, and dynamic shock compression tests were carried out in this paper, and the static and dynamic parameters of 60° jointed slate were obtained. The HJC constitutive model was selected as the jointed slate material model, the method for determining the parameters of the HJC constitutive model was provided, and the parameters of the HJC constitutive model were determined.

Considering the geometric attenuation and physical attenuation of the blast stress wave, the attenuation formula of the blast stress wave under the combined action of in situ stress and joints was deduced. Using LS-PREPOST software, a three-dimensional numerical calculation model under the combined action of in situ stress and joints was established, and the radial stress wave peaks at 0.6 m, 1.2 m, 1.8 m, and 2.4 m from the center of the blast hole were monitored. These were consistent with the predicted value of the explosion stress wave attenuation formula proposed in this paper, verifying the correctness of the attenuation formula. Based on the theory of the interaction between stress waves and explosive gas, a peripheral hole spacing calculation formula was proposed that comprehensively considers the effects of in situ stress, joints, and rock mass tensile strength.

Combined with the geological conditions and blasting parameters of the Bayueshan Tunnel study section, it was calculated that the optimal peripheral hole spacing around the study section was 45 cm. A three-dimensional numerical analysis model was established using LS-PREPOST software. The damage depth of the retained rock mass and the peak vibration velocity (PPV) were compared and analyzed when the peripheral hole spacings were 40 cm, 45 cm, 50 cm, and 55 cm, respectively. The numerical simulation results verified the correctness of the peripheral hole spacing calculation formula proposed in this paper.

Taking the 45 cm peripheral hole spacing and the 56 cm resistance line as the standard, a field blasting test was carried out in the research section of the Bayueshan Tunnel. Under the combined action of in situ stress and joints, blasting construction was carried out using the peripheral hole layout method and the hole mesh layout parameters proposed in this paper. The test results show that the average over-excavation value of the grade IV surrounding rock was controlled at 22 cm, the concrete over-consumption per linear meter was controlled within 100%, and the concrete over-consumption control effect was good.

The datasets used during the current study available from the corresponding author on reasonable request.

Skrzypkowski, K., Korzeniowski, W., Zagórski, K. & Dudek, P. Application of long expansion rock bolt support in the underground mines of Legnica-Głogów copper district. Stud. Geotech. Mech. 39, 47–57. https://doi.org/10.1515/sgem-2017-0029 (2017).

Article ADS Google Scholar

Skrzypkowski, K., Korzeniowski, W., Zagórski, K. & Zagórska, A. Modified rock bolt support for mining method with controlled roof bending. Energies 2020, 13. https://doi.org/10.3390/en13081868 (1868).

Article CAS Google Scholar

Wang, J. X. et al. Numerical simulation of rock mass structure effect on tunnel smooth blasting quality: A case study. Appl. Sci. 2021, 11. https://doi.org/10.3390/app112210761 (2021).

Article CAS Google Scholar

Paul, S. & Peter, X. Causes, impact and control of overbreak in underground excavations. Tunn. Undergr. Space Technol. 20(5), 63–71. https://doi.org/10.1016/j.tust.2004.05.004 (2005).

Article Google Scholar

Adel, M., Farhang, S. & Mohammad, A. Overbreak prediction in underground excavations using hybrid ANFIS-PSO model. Tunn. Undergr. Space Technol. 8(6), 1–9. https://doi.org/10.1016/j.tust.2018.05.023 (2018).

Article Google Scholar

Yin, W. G., Wang, H. L., Hu, H. X., Zhang, F. Q. & Chen, J. H. Research on tunnel forming blasting control technology and surrounding rock damage range. Tunnel Constr. (Chin. Engl.) 38(5), 851–856 (2018).

Google Scholar

Liu, D., Gao, W. X. & Liu, M. G. Causes and control technology of tunnel over and under excavation. J. Undergr. Space Eng. 3(s2), 1468–1471 (2007).

Google Scholar

Huang, Z. Q., Cui, Z. Y., Liu, D. & Yang, X. F. Application of double-layer peripheral hole retaining wall blasting in weak and broken surrounding rock. Blasting 2022, 1–12 (2022).

CAS Google Scholar

Man, K. & Liu, X. L. Influence of surrounding void effect on blasting vibration velocity and surrounding rock damage. Eng. Mech. 37(11), 127–134 (2020).

Google Scholar

Man, K., Liu, X. L., Wang, X. Y., Chen, Z. M. & Wang, J. Influence of blasthole distance of peripheral holes on smooth blasting effect. Sci. Technol. Eng. 16(29), 47–53 (2016).

Google Scholar

Sun, G. Z. Mechanics of Rock Mass Structure (Science Press, 1988).

Google Scholar

Hyongdoo, J., Youhei, K. & Utsuki, S. An empirical approach of overbreak resistance factor for tunnel blasting. Tunn. Undergr. Space Technol. 13(3), 1–10. https://doi.org/10.1016/j.tust.2019.103060 (2019).

Article Google Scholar

Zhao, J. & Cai, J. G. Transmission of elastic P-waves across single fractures with a nonlinear normal deformational behavior. Rock Mech. Rock Eng. 34(1), 3–22. https://doi.org/10.1007/s006030170023 (2001).

Article MathSciNet Google Scholar

Cai, J. G. & Zhao, J. Effects of multiple parallel fractures on apparent attenuation of stress waves in rock masses. Int. J. Rock Mech. Min. Sci. 37(4), 661–682. https://doi.org/10.1016/s1365-1609(00)00013-7 (2000).

Article Google Scholar

Zhao, J., Zhao, X. B. & Cai, J. G. A further study of P-wave attenuation across parallel fractureswith linear deformational behaviour. Int. J. Rock Mech. Min. Sci. 43, 776–788. https://doi.org/10.1016/j.ijrmms.2005.12.007 (2006).

Article Google Scholar

Zhao, X. B., Zhao, J. & Cai, J. G. P-wave transmission across fractures with nonlineardeformational behaviour. Int. J. Numer. Anal. Meth. Geomech. 30(11), 1097–1112. https://doi.org/10.1002/nag.515 (2006).

Article MATH Google Scholar

Li, J. C. & Ma, G. W. Analysis of blast wave interaction with a rock joint. Rock Mech. Rock Eng. 43(6), 777–787. https://doi.org/10.1007/s00603-009-0062-0 (2010).

Article ADS Google Scholar

Li, J. C., Li, H. B., Ma, G. W. & Zhao, J. A time-domain recursive method to analyse transient wave propagation across rock joints. Geophys. J. Int. 188(2), 631–644. https://doi.org/10.1111/j.1365-246x.2011.05286.x (2012).

Article ADS Google Scholar

Li, J. C. Wave propagation across non-linear rock joints based on time-domain recursive method. Geophys. J. Int. 193(2), 970–985. https://doi.org/10.1093/gji/ggt020 (2013).

Article ADS Google Scholar

Perino, A., Orta, R. & Barla, G. Wave propagation in discontinuous media by the scatteringmatrix method. Rock Mech. Rock Eng. 45(5), 901–918. https://doi.org/10.1007/s00603-012-0286-2 (2012).

Article ADS Google Scholar

Chai, S. B., Li, J. C., Zhang, Q. B., Li, H. B. & Li, N. N. Stress wave propagation across a rock mass with two non-parallel joints. Rock Mech. Rock Eng. 2016, 1–10. https://doi.org/10.1007/s00603-016-1068-z (2016).

Article Google Scholar

Mandal, S. & Singh, M. Evaluating extent and causes of overbreak in tunnels. Tunn. Undergr. Space Technol. 24(1), 22–36. https://doi.org/10.1016/j.tust.2008.01.007 (2009).

Article Google Scholar

Mandal, S., Singh, M. & Dasgupta, S. Theoretical concept to understand plan and design smooth blasting pattern. Geotech. Geol. Eng. 26(4), 399–416. https://doi.org/10.1007/s10706-008-9177-4 (2008).

Article Google Scholar

Zheng, Z. T. Research on 3D model test of ultra-deep hole columnar charge blasting in rock mass with high geostress. In Anhui University of Science and Technology (2017).

Zhang, C., Yang, J. H., Yao, C. & Chen, Y. Y. Stress evolution pattern of surrounding rock during continuous blasting excavation of deep tunnels. J. Changjiang Acad. Sci. 35(12), 108–112 (2018).

Google Scholar

Li, X. P. et al. Experimental study on the pattern of elastic wave propagation and attenuation in deeply fractured rock mass. Chin. J. Rock Mech. Eng. 34(11), 2319–2326. https://doi.org/10.13722/j.cnki.jrme.2015.0868 (2015).

Article Google Scholar

Li, X. Q. et al. Model test study on propagation pattern of explosion stress wave in jointed rock mass under different in-situ stresses. Chin. J. Rock Mech. Eng. 35(11), 2188–2196. https://doi.org/10.13722/j.cnki.jrme.2016.1004 (2016).

Article Google Scholar

Fan, L. F. & Sun, H. Y. Seismic wave propagation through an in-situ stressed rock mass. J. Appl. Geophys. 121, 13–20. https://doi.org/10.1016/j.jappgeo.2015.07.002 (2015).

Article ADS Google Scholar

Liu, Y. & Xu, J. Y. Numerical simulation of rock mass explosion under in-situ stress field. Geotechn. Mech. 28(11), 2485–2488. https://doi.org/10.16285/j.rsm.2007.11.042 (2007).

Article Google Scholar

Fan, X., Wang, M. Y. & Shi, C. C. Influence of initial stress on stress wave propagation and block motion pattern. Chin. J. Rock Mech. Eng. 28(2), 3442–3446 (2009).

Google Scholar

Zhang, F. P., Peng, J. Y., Fan, G. H., Li, S. J. & Li, Y. H. Research on rock blasting mechanism of rock mass under different static stress and joint conditions. Rock Soil Mech. 37(7), 1839–1846. https://doi.org/10.16285/j.rsm.2016.07.003 (2016).

Article Google Scholar

Huang, H. Experimental Study on Double Uniaxial Strength Ratio of Concrete at Early Age (Dalian University of Technology, 2013).

Cui, X. Z., Chen, S. H. & Liu, D. C. Attenuation mechanism of stress wave propagating in fractured rock mass. Eng. Blast. 5(1), 18–21 (1999).

Google Scholar

Wang, L. L. Fundamentals of Stress Waves-Second Edition (National Defense Industry Press, 2005).

Google Scholar

Liu, Z. Y., Hu, Y. L. & Chen, Y. Attenuation of ultrasonic P-wave in dry and water-saturated rocks under uniaxial compression. Acta Geophys. 27(4), 349–359 (1984).

Google Scholar

Dong, Q. Study on Propagation and Attenuation Pattern of Explosion Stress Wave in Jointed Rock Mass Under Different In-Situ Stress (Wuhan University of Technology, 2018).

Dong, Q., Li, X., Jia, Y. & Sun, J. A numerical simulation of blasting stress wave propagation in a jointed rock mass under initial stresses. Appl. Sci. 11(17), 7873. https://doi.org/10.3390/app11177873 (2021).

Article CAS Google Scholar

Li, X. B. & Gu, D. S. Rock Impact Dynamics (Central South University of Technology Press, 1994).

Google Scholar

Holmquist, T., Johnson, G. & Cook, W. A computational constitutive model for concrete subjective to large strain, high strain rates, and high pressure. In 1th International Symposium on Ballistic, Quebec City, Canada 591–600 (1993).

Fang, Q., Kong, X. Z., Wu, H. & Gong, Z. M. Method for determining parameters of rock Holmquist Johnson cook model. Eng. Mech. 31(03), 197–204 (2014).

Google Scholar

Ding, H. P. Research on Blasting Effect of Jointed and Fractured Rock Mass Tunnel (Jilin University, 2009).

Liu, J. C. et al. Numerical analysis and application of structural optimization of water-sealed blasting charge. Vibr. Shock 39(09), 57–62. https://doi.org/10.13465/j.cnki.jvs.2020.09.008 (2020).

Article Google Scholar

Wang, P. Optimization of Blasting Parameters for Wedge-Shaped Undercuts in Large-Section Tunnels (Guangxi University, 2020).

Li, F. Q. Numerical Simulation of Wedge-Shaped Undercut Blasting in Hollow Holes (Beijing Institute of Technology, 2018).

Wang, J., Yin, Y. & Luo, C. Rock materials: Parameter determination and application in tunnel smooth blasting. Appl. Sci. 8, 9 (2018).

ADS Google Scholar

Murthy, V., Dey, K. & Raitani, R. Prediction of overbreak in underground tunnel blasting: A case study. J. Can. Tunn. Canadien 2003, 109–115 (2003).

Skrzypkowski, K. A new design of support for burst-prone rock mass in underground ore mining. In E3S Web Conferences; EDP Sciences: Les Ulis, France, volume 71 00006. https://doi.org/10.1051/e3sconf/20187100006 (2018).

Skrzypkowski, K., Korzeniowski, W., Krzysztof, Z. & Anna, Z. Adjustment of the yielding system of mechanical rock bolts for room and pillar mining method in stratified rock mass. Energies 13(8), 2082. https://doi.org/10.3390/en13082082 (2020).

Article CAS Google Scholar

Persson, P. The relationship between strain energy, rock damage, fragmentation, and throw in rock blasting. Fragblast 1(1), 99–110 (1997).

Article Google Scholar

Wang, J. X. et al. Numerical simulation of rock mass structure effect on tunnel smooth blasting quality: A case study. Appl. Sci.-Basel 11, 22. https://doi.org/10.3390/app112210761 (2021).

Article ADS CAS Google Scholar

Wang, J. X., Yin, Y. & Luo, C. W. Johnson–Holmquist-II (JH-2) constitutive model for rock materials: Parameter determination and application in tunnel smooth blasting. Appl. Sci.-Basel 8, 9. https://doi.org/10.3390/app8091675 (2018).

Article ADS CAS Google Scholar

Download references

Thanks to Poly Xinlian Blasting Engineering Group Civil Explosion Engineering Laboratory and Poly Union Group Corporation Civil Explosion Engineering Laboratory for providing experimental venues for this paper. This work was supported by the National Natural Science Foundation of China (52064008).

College of Civil Engineering, Guizhou University, Guiyang, 550025, China

Xingchao Tian, Xia Liu, Jian Jia, Caijin Xie, Qianxing Lou & Qingzhi Chen

College of Mining, Guizhou University, Guiyang, 550025, China

Tiejun Tao

Northwest Engineering Co., Ltd. of CCCC First Highway Engineering Co., Ltd., Xi’an, 710000, China

Zhenhua Zhao

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

X.-C.T.: theoretical derivation, method validation, writing original draft. T.T.: editing and checking of manuscript. X.L.: model parameter assignment. J.J.: uniaxial compression test, triaxial compression test, Brazilian splitting test, and dynamic shock compression test. C.-J.X.: formula validation. Q.-X.L., and Q.-Z.C.: the establishment of numerical analysis model. Z.-H.Z.: field blasting Test.

Correspondence to Tiejun Tao.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Tian, X., Tao, T., Liu, X. et al. Calculation of hole spacing and surrounding rock damage analysis under the action of in situ stress and joints. Sci Rep 12, 22331 (2022). https://doi.org/10.1038/s41598-022-27017-w

Download citation

Received: 19 August 2022

Accepted: 23 December 2022

Published: 25 December 2022

DOI: https://doi.org/10.1038/s41598-022-27017-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.