Effects of Unidirectional In Situ Stress on Crack Propagation of a Jointed Rock Mass Subjected to Stress Wave

Fan, Zhanfeng; Cai, Jianhua

doi:https://doi.org/10.1155/2021/5529540

Shock and Vibration

On this page

Abstract Introduction Results and Discussion Conclusions Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5529540 | https://doi.org/10.1155/2021/5529540

Effects of Unidirectional In Situ Stress on Crack Propagation of a Jointed Rock Mass Subjected to Stress Wave

Zhanfeng Fan^1,2and Jianhua Cai³

Academic Editor: Ahmed Ibrahim

Received21 Feb 2021

Revised27 Jul 2021

Accepted13 Aug 2021

Published30 Aug 2021

Abstract

This paper proposes a large-scale experiment combined with multiple cement mortar blocks to simulate stress wave propagation across a jointed rock mass under unidirectional in situ stress. Two identical mortar block models with smooth, dry, and unfilled joints were poured. The stress waves in Model 1 and Model 2 were generated by an electric spark source and a blast-induced source, respectively. The effects of these two excitation sources on stress wave propagation were compared through crack propagation experiments. The experimental results show that the peak value of the transmitted strain wave decreases as the in situ stress increases. The unidirectional in situ stress has a certain inhibitory effect on the stress wave propagation. It also indicates that for Model 1 with the electric spark source, no cracks on the upper surface, but a Livingstone blasting crater at the bottom is generated. For Model 2 with the blast-induced source, cracks on the upper surface and a blasting crater at the bottom are produced. The results verify the similarity between the electric spark source and the explosive source. The two-dimensional finite element program (ANSYS/LS-DYNA) was applied to further simulate the crack propagation of a jointed rock mass under different in situ stresses. The results of numerical simulation verify that the in situ stress has a clear guiding effect on the crack propagation.

1. Introduction

In recent years, the rapid development of tunneling and hydropower engineering has greatly promoted tunnel excavations to deeper depth [1]. Through this work, high in situ stress has been identified as a characteristic property of deep jointed rock masses [2–4]. Therefore, the stress wave propagation in a jointed rock mass is significantly influenced not only by the joints but also by the in situ stress. In addition, the influence of in situ stress on the stability and failure of the surrounding rock mass at depth has become increasingly clear, improving the prediction of rockburst, large deformation, and complex tunnel geology. The behavior of rock masses is also governed by the in situ geological structure and stress conditions [5]. Consequently, the study of the stress wave propagation across jointed rock masses under in situ stress conditions has continued to attract considerable attention.

A considerable amount of literature has focused on the following two aspects of stress wave propagation through jointed rock masses: the properties of stress waves and the nature of joints. Research on the properties of stress waves includes the incident angles [6, 7], the amplitude and frequency of incident waves [8], and different types of stress waves [6, 9]. Researchers have also examined the nature of joints, including the joint closure constitutive relationship [10], the joint opening-closing behavior [9, 11], the number of joints [7, 12, 13], the joint stiffness [14], and the joint filling properties [15, 16], etc. Among them, most of these studies adopted small-scale laboratory experiments, i.e., SHPB (split Hopkinson pressure bar) experiments. In general, the diameter and length of the cylindrical sample of SHPB are 50 mm and 25 mm, respectively. These experiments have significant side effects. Few studies have considered large-scale experiments under the coupling of in situ stress and dynamic stress loading.

The existing literature has indicated that in situ stress has strong effects on the failure strength and dynamic crack propagation of rock masses [17–19]. For instance, according to a biaxial compression test of a jointed rock mass with an opening, the propagation direction of the tensile cracks was found to be roughly normal to the joint surfaces [20]. High in situ stress was observed to resist the radially oriented pressure and damage extension around cut holes during tunneling excavation [21]. The levels and states of the in situ stress in the rock greatly affected the damage evolution, the crack initiation, and development and hence affected the fragmentation [22]. He et al. [23] have experimented with granite under several static and dynamic stresses to examine the crack propagation characteristics of blast-induced shock waves, and the results showed that the direction of the radial crack was related to the lateral pressure coefficient in the biaxial direction. The in situ stress could enhance the compression effect and weaken the tension effect in the radial direction of the borehole, and the jointed plane failure zone in coal masses decreased with increasing in situ stress [24]. Blasting of deeply buried structures is a complex dynamic process that can cause blast loading and in situ stress redistribution [25]. These studies mainly investigated the crack propagation of intact rock masses or rock masses containing microcracks or microvoids. The influence of the joints, large in extent and small in thickness compared with the wavelength of an incident wave, on crack propagation is rarely found in the available literature.

In addition to laboratory test methods, numerical simulation is an important method to investigate the blasting fracture mechanism. Currently, there are many numerical simulation methods, such as the finite element method (FEM), discrete element method (DEM), and boundary element method (BEM). Compared with the other two numerical methods mentioned above, FEM software ANSYS/LS-DYNA is advantageous in that the explicit dynamic analysis method can be used conveniently to simulate the large deformation of the material, and the explosives and other parameters can be input directly with a user subroutine interface. Since the FEM with a contact face can simulate crack growth in both the intact mass and the jointed plane, it may be a good method to study the fracture pattern in a jointed rock mass during the blasting process [26]. For example, the prediction of blast craters in brittle rock due to an underground explosion has been simulated using the LS-DYNA based on the Taylor–Chen–Kuszmaul continuum damage model [27]. The Johnson–Holmquist material model was implemented in LS-DYNA through user subroutines to simulate the blasting-induced rock fractures [28]. A constitutive model considering dynamic compressive and tensile failure was applied to analyze crack initiation and propagation in coal masses using LS-DYNA [26]. The propagation characteristics of shock waves subjected to blasting load in a parallel jointed rock mass were investigated using LS-DYNA, and the results showed that the joint properties had a great effect on the propagation of blast-induced vibration [29]. In the present study, ANSYS/LS-DYNA will be employed for crack propagation of a jointed rock mass under in situ stress.

As mentioned above, both the joints and the in situ stress influence the stress wave propagation across a jointed rock mass and consequently the crack propagation. Researchers have also realized the importance of these two factors in rock masses, but few have considered laboratory studies of wave propagation under unidirectional in situ stress. The theme of this study is to investigate crack propagation of a jointed rock mass caused by stress waves using two experimental models under unidirectional in situ stress. First, to simulate the joint distribution and apply different in situ stresses, a model consisting of multiple cement mortar blocks is presented. Then, we adopt a ZDF-3 (ZhiliuDian Fire-3) electric spark source to excite stress waves. This excitation source is compared with the blast-induced source through the law of crack propagation and the shape and size of the blasting crater in the model. Finally, the dynamic relaxation method of ANSYS/LS-DYNA is applied to study the law of crack propagation in jointed rock masses under unidirectional in situ stress. Our procedure provides a systematic methodology to quantify the crack propagation in a jointed rock mass under high in situ stress.

2. Laboratory Experiment

2.1. Experimental Model

A large-scale experimental model with multiple joints under unidirectional in situ stress was designed, as shown in Figure 1. The model materials mainly selected cement with a compressive strength of 30 MPa and medium sand. The ratio of water, cement, and medium sand was 1 : 4 : 0.5. The experiment did not consider the influence of the self-weight stress on stress wave propagation because the self-weight stress was very small compared to the in situ stress. Accordingly, it was not necessary to consider the weight similarity ratio of the model material.

(a)

(b)

(c)

The model was composed of a parent body and four jointed blocks. Hence, four parallel joints were formed into the model. A single joint was located on the left side of the parent body, and the remaining parallel joints were located on the right side. The experimental model mainly considered smooth, dry, and unfilled joints which mean no thickness joints. It is noted that the joints here are large in extent and small in thickness compared with the wavelength of an incident wave. The total size of the experimental model was 180 cm × 80 cm × 25 cm (length × width × height) and much larger than the size of the sample of SHPB. Among them, the length, width, and height of the parent body were 80 cm, 60 cm, and 25 cm, respectively, and the dimensions of the jointed block were the same as those of the parent body except that the width was 30 cm, as shown in Figure 1(a). Two loads were applied along the length direction to simulate in situ stress. The other two directions had a free boundary. The size of the model is large enough to avoid the influence of reflected stress wave propagation [23, 30, 31]. Eight monitoring points were arranged sequentially from left to right, numbered T1 to T8, as shown in Figure 1(b). As the distance of stress wave propagation increased, it was assumed that the stress wave measured at each pressure measurement point was a plane wave. To avoid the overflow of explosion gas from the bottom of the borehole during the experiment, a nonpenetrative blast hole was precut in the middle of the parent body and was surrounded by a cylindrical steel pipe with a base (5 cm × 12.5 cm), as shown in Figure 1(c).

The poured and assembled model consisted of the unidirectional in situ stress loading device, the dynamic stress loading device, and the dynamic strain acquisition system, as shown in Figure 2. The physical and mechanical parameters of the cement mortar block are listed in Table 1, where denotes the density of the cement mortar block; and denote Young’s modulus and Poisson’s ratio, respectively; , , and denote the uniaxial compressive strength, the tangent modulus, and the wave velocity, respectively.

2.2. ZDF-3 Electric Spark Source

Generally, the explosive is used as an excitation source in blasting tests. Once the explosive is detonated, the stress wave and explosion gas will coexist at the same time and are difficult to separate. Although explosive produces a wide frequency spectrum and a large amount of energy, it has poor controllability and is thus dangerous to work with. To study only the attenuation law of the stress wave, a new type of ZDF-3 (ZhiliuDian Fire-3) electric spark was selected as the excitation source, as shown in Figure 2(a).

The ZDF-3 electric spark source produces large energy and rich high-frequency components. The working principle of the ZDF-3 electric spark is to insert a special launch head into a blast hole filled with water, as shown in Figure 2(b), and then quickly charge the capacitor. Later, it uses ultrahigh voltage to instantly release energy into the water and thus generate stress waves. Compared with explosive detonation, its advantage is that only stress waves are produced without any explosion gas. The frequency range of the electric spark is 40 Hz to 300 Hz. The spark energy can be expressed as follows [32]:where , (=420 F), and are the emission energy, the capacitor capacitance, the voltage of the ZDF-3 electric spark, respectively. Due to the repeated application of dynamic loading to the model, a thick seamless steel pipe was embedded in the parent body of the model, as shown in Figure 2(c). The density of the seamless stainless steel pipe is 7900 kg/, Young’s modulus is 194.0 GPa, and Poisson’s ratio is 0.25, respectively. Meanwhile, eight steel bars were welded outside the stainless steel pipe to prevent it from separating from the cement mortar block.

2.3. Loading Equipment

An eight-channel ultradynamic data acquisition machine (DH5930) was used to collect dynamic strain measurements, as shown in Figure 2(d). The loading equipment was composed of two reaction frames, six steel strands, a ZB-500 electric fuel pump, and two center hole jacks, as shown in Figure 2(e)–2(g). The reaction frames were installed on the left and right sides of the model, and three steel strands were placed on the other sides. Uniform loads can be applied to the boundary of the model, and the maximum load was 5.0 MPa. To load the in situ stress, a ZB-500 electric fuel pump was used to control two center hole jacks to tension the steel strand. Two strain gauges were pasted on each small strain brick embedded in the model shown in Figure 2(h).

2.4. Experimental Process

Uniform loads of the same magnitude were applied in one direction through the loading equipment including the static and dynamic stages. The static loading onto the cement mortar block was approximately symmetrically distributed, and the stress wave propagated evenly from the center blast hole of the parent body. Two identical mortar block models were poured. The experimental process of Model 1 was first implemented when using the ZDF-3 electric spark as the excitation source, as shown in Figure 3. The experimental process of Model 2 is the same as that of Model 1 except that blasting in the blast hole is used as the excitation source.

During the static stage, the in situ stress was loaded in one direction with the loading equipment. Six different magnitudes of the unidirectional in situ stress, including 0.5 MPa, 1.0 MPa, 2.0 MPa, 3.0 MPa, 4.0 MPa, and 5.0 MPa, were considered. Each in situ stress loading was held for more than 10 minutes or until the cement mortar block was fully deformed.

During the dynamic stage, the blast hole was filled with water, and an electric spark cable with a copper core was inserted into the blast hole. The blast hole was blocked with plasticine and covered with an iron plate to seal the blast hole. Then, the ZDF-3 electric spark source was turned on and charged to the required voltage. Once the voltage was completely stable, it was suddenly released. At this point, the dynamic-static loading was complete. The dynamic strain curve of the measuring point was recorded through the DH5930 apparatus.

In the present study, three charging voltages were considered: 2500 V, 3000 V, and 3500 V. Each voltage was excited three times under each level of in situ stress, and thus 18 sets of dynamic strain data were obtained. Therefore, a total of 54 experiments were completed. When the electric spark was charged to 2500 V, the emission energy was obtained in Table 2. The maximum and minimum average emission energies were 746.3 J and 665.5 J, respectively. Overall, the electric spark excitation energy is adjustable and controllable. It is feasible to use the ZDF-3 electric spark source to replace the blast-induced source with exciting stress waves.

3. Results and Discussion

In this section, the dynamic strains of monitoring points T1 and T8 in Model 1 using the ZDF-3 electric spark source are presented. After that, the crack propagation patterns of the parent bodies of Model 1 using the ZDF-3 electric spark source and Model 2 using the blast-induced source are analyzed.

3.1. Dynamic Strain

Figure 4 shows the waveforms of the transmitted strain wave at T1 under unidirectional in situ stress when the stress wave passes through a single joint, where the symbol “-” indicates compression and “+” indicates tension. The charging voltage is 2500 V. Regardless of the in situ stress, the transmitted strain wave first exhibits a compressive wave and then a reflected tensile wave. When the in situ stress varies from 0.5 MPa to 5.0 MPa, the peak value of the transmitted strain wave decreases from 9.4E4 to 5.28E4 , and the time interval required to reach the peak of the transmitted strain wave is 217.1 s, 221.3 s, 216.4 s, 213.5 s, 211.6 s, 193.1 s, respectively. The greater the in situ stress is, the smaller the dynamic deformation of the rock mass caused by the stress wave. Therefore, the in situ stress has a certain inhibitory effect on stress wave propagation when the initial in situ stress and the stress waves are superimposed. Although the increasing in situ stress will lead to the closure of the joint to a certain extent, the jointed rock mass may produce new microcracks and plastic deformation under the combined action of dynamic and static stress.

Figure 5 shows the waveforms of the transmitted strain wave at T8 under different in situ stresses when the stress wave passes through multiple parallel joints. The attenuation law of the transmitted strain wave of the three joints is similar to that of the single joint; that is, the compressive wave appears first, and then the tensile wave appears. There are subsequent peaks after the first peak. The peak value of the transmitted strain wave decreases from 2.46E4 to 1.2E4 when the in situ stress gradually increases from 0.5 MPa to 5.0 MPa. The results show that the unidirectional in situ stress has a great influence on the stress wave propagation and attenuation. When the in situ stress is the same, the dynamic strain of T8 is significantly less than that of T1 because of the increase in the number of joints.

3.2. Crack Propagation of Model 1

The crack distribution patterns on the upper and bottom surfaces of the parent body of Model 1 are shown in Figure 6. No cracks are observed on the upper surface of Model 1, as shown in Figure 6(a). However, radial cracks and circumferential cracks appear on the lower surface of the matrix along the center of the model, as shown in Figure 6(b). Additionally, a Livingstone blasting crater formed during the experiment, and the funnel body is completely separated from the parent body. The diameter of the crater is approximately 52 cm, and the depth is 10 cm. The results of Model 1 show that a single stress wave loading may not be enough to cause damage to the jointed rock mass, but the repeated stress waves loading can cause the jointed rock mass instability and failure.

(a)

(b)

3.3. Crack Propagation of Model 2

Model 2 was blasted nine times, and the detonator was buried 12.5 cm deep. Most blast-induced sources used a delay detonator and water as the coupling medium, except for the second blast, which used a detonator and 4.9 g of explosives and blocked the blast hole with mud. The in situ stress loading increased from 0.5 MPa to 5.0 MPa. The specific cases and the parameters of filling delay detonator and explosive in Model 2 are listed in Table 3. The evolution of the crack distribution in the parent body of Model 2 from the first to ninth blast is shown in Figure 7. The following test results were observed:(1)No cracks were observed around the blast hole after blasting when the in situ stress loading was 0.5 MPa and 1.0 MPa, respectively. An initially very small crack around the blast hole had appeared since the 3rd blasting when the in situ stress was 2.0 MPa. After the 5th blasting, the crack tended to propagate along the loading direction when the in situ stress reached 3.0 MPa. This phenomenon indicates that the initial in situ stress can affect crack propagation. A superposition effect occurs if the stress wave encounters the in situ stress. Other researchers have reported similar conclusions. For example, Nicholls and Duvall [33] pointed out that the presplitting of test specimens along with the maximum in situ stress under confining pressure much more easily passed field tests than presplitting in other directions. Simha and Fourney [34] proposed the static stress field to blast fracture mode of testing and concluded that the preferred fracturing direction of polymethyl methacrylate (PMMA) was the direction of maximum principal stress.(2)The propagation direction of the cracks around the blast hole was not completely parallel to the loading direction but was at an angle to the loading direction. This phenomenon occurred after the 3rd blast, and the cracks became longer and wider with increasing in situ stress, as shown in Figure 7. For the so-called inclined deflection cracks, the angle between the crack and the in situ stress was between 15° and 45°. In other words, the crack propagation direction intersected the in situ stress loading direction.(3)From the first blast to the end of the experiment, no cracks were observed perpendicular to the loading direction. Due to the steel pipe surrounding the blast hole, the annular crack around the blast hole was very small, and the steel pipe also blocked the explosion gas from wedging into the cement mortar block. The extent of the blasting damage zone around the blast hole decreased as the in situ stress increased.(4)After nine blasts, an irregular circular area of deformation appeared at the bottom of Model 2, as shown in Figure 8. The diameter of the profile of the blasting crater was approximately 30 cm. This crack connected with the upper surface crack on the side of Model 2 (the yellow dotted line). However, the funnel body did not fall off.

The above results show that the in situ stress state of the rock mass greatly affects the damage evolution and crack initiation and development and therefore affects fragmentation [22]. The damage development tends towards the direction at which the unidirectional in situ stress is applied. In addition, the crack propagation patterns of Model 1 and Model 2 verify the similarity between the ZDF-3 electric spark source and the blast-induced source. Although the energy of the ZDF-3 electric spark is not as large as that of the explosive, it can play an important role in stress wave propagation and can be used as an excitation source in stress wave propagation.

4. Numerical Simulation of Crack Propagation

Since it is difficult to obtain an analytical solution for nonelastic materials under blasting, there has been a considerable amount of researches on the numerical simulation of blasting funnels [21, 27, 35, 36]. To further confirm that the blasting cracks in the jointed rock mass propagate along the loading direction under dynamic and static loading, a numerical simulation to take into account the nonlinear response of jointed rock masses is conducted in this research. A plane strain numerical calculation model is established using the finite element program ANSYS/LS-DYNA. The calculation model is the same size as the parent body of the models described in Section 2. Only unidirectional in situ stress (denoted as ) is considered in the calculation, as shown in Figure 9.

When the explosive is detonated, the strain of the rock mass near the charge is very large, and the strain rate effect is obvious. Therefore, the Cowper-Symonds model is considered. The elastic-plastic-kinematic model considers the effect of strain rate and strain hardening on the rock mass and can accurately describe the fracture pattern during blasting [37–39]. The expression is written as follows [40]:where , and are the same as before, is the initial stress, is the failure stress, is the strain rate, is the plastic hardening modulus, is the hardening parameter (), and and are Cowper-Symonds strain rate parameter, respectively. The effective plastic strain can be calculated as follows:where is the effective plastic strain, is the cumulative time of plastic strain, and is the partial component of plastic deviation strain of the rock mass. In this simulation, = 2.63; = 3.69; = 0.15; = 1.25.

In LS-DYNA, the Jones-Wilkens-Lee (JWL) equation of state is used to describe the chemical reactions after explosive detonation [40]. The JWL equation of state defines the pressure as follows:and is usually used for detonation products of high explosives, where is the pressure of the detonation products. is the relative specific volume of the detonation products. is the specific internal energy of the detonation products. and are parameters of JWL equation of state, respectively. In this study, the parameters of the JWL equation of state are listed in Table 4, where is the explosive density, is the Chapman-Jouguet pressure, and VoD is the velocity of detonation.

Four cases are considered: (i) Case I: = 0 MPa, and a steel pipe is either included or excluded in the model; (ii) Case II: = 0.5 MPa, and a steel pipe is included; (iii) Case III: = 1.0 MPa, and a steel pipe is included; (iv) Case IV: = 2.0 MPa, and a steel pipe is included. The parameters of the numerical calculation model are listed in Table 1. The SOLID185 with eight nodes of hexahedral element is adopted. Only two elements are set in the thickness direction due to the plane strain assumption. The numerical simulations have a total of 64 250 elements. In the finite element method, the joint is usually treated as a boundary. The boundary conditions were specified as follows: the upper and lower sides are nonreflecting boundaries, while the other two sides are free faces. The numerical simulation adopts the cm-g-s unit system.

4.1. Case I

Case I mainly compares the damage zones around the blast hole with and without a steel pipe. In this case, only the dynamic calculation steps of ANSYS/LS-DYNA are required due to the lack of in situ stress. The blast hole without a steel pipe is completely destroyed, and a crushing area is formed, as shown in Figure 10(a). The symmetrical cracks are distributed at the periphery of the crushing zone. The distribution of the damage zone around the blast hole with a steel pipe is shown in Figure 10(b). Compared with that in Figure 10(a), the range of the damage zone in Figure 10(b) is significantly reduced. Only small cracks around the steel pipe are observed, and the steel pipe separates from the parent body. The maximum dynamic strain is 1.87E4 . These findings explain why eight steel bars should be welded outside the stainless steel pipe in the model. The explosion gas does not wedge into the cement mortar block after the explosive is detonated, further confirming that the stress wave can damage the rock mass.

(a)

(b)

4.2. Case II

In Case II, the unidirectional in situ stress is loaded first, and then a dynamic loading is applied on the blast hole. The dynamic relaxation method of the ANSYS/LS-DYNA is used to obtain the nodal stress and displacement of each element under the initial in situ stress before and after blasting [24, 41, 42]. The distribution of the stress and the damage zone around the blast hole with unidirectional in situ stress ( = 0.5 MPa) is shown in Figure 11. At the static loading stage, the distribution of the maximum tensile and compressive stress is shown in Figure 11(a). The upper and lower sides of the rock mass are compressed, and the left and right sides are subjected to tension only when unidirectional in situ stress is applied. At the dynamic loading stage, the stress of each element obtained from the static analysis is input into the initial stress of the dynamic analysis, and then the explosion analysis is performed. Very small cracks form around the blast hole, as shown in Figure 11(b). The maximum dynamic strain is 3.34E4 . The range of crack propagation is small because the applied unidirectional in situ stress is low.

(a)

(b)

4.3. Case III

At the static stage, the maximum tensile and compressive stress under the action of the static stress field in the rock mass tends to increase when the unidirectional in situ stress increases to 1.0 MPa, as shown in Figure 12(a). After approximately 234.96 s of dynamic loading, the distribution of the damage zone around the blast hole is shown in Figure 12(b). The maximum dynamic strain is 6.68E4 . Two distinctly continuous cracks appear along the loading direction and then extend to the boundary. It is shown that the in situ stress has a clear guiding effect on the stress wave propagation. The crack propagation perpendicular to the direction of the unidirectional in situ stress is not obvious. These results are consistent with the crack propagation of Model 2 in Section 3. The difference between the two is that the crack propagation direction observed in Model 2 is at an angle to the loading direction, but the two cracks in the numerical calculation are parallel to the loading direction. The reason for this discrepancy may be related to the presence of microcracks and the plane strain assumption.

(a)

(b)

4.4. Case IV

The maximum tensile and compressive stress increase, and the cracks propagate farther when the unidirectional in situ stress is continuously increased to 2.0 MPa, as shown in Figure 13. The maximum dynamic strain is 9.38E4 . In addition to the crack observed in Case III, two reflective tensile cracks form along the direction of the unidirectional in situ stress loading. The reason for this phenomenon may be that the nonreflecting boundaries cannot completely absorb the reflected waves. The crack is still not observed in the vertical direction of in situ stress loading, further verifying that the in situ stress has a guiding effect on crack propagation. Xie et al. [21] indicate that the initial in situ stress has the effect of increasing the tensile stress and reducing the compressive stress; that is, the initial in situ stress has the effect of inhibiting the development of radial cracks.

(a)

(b)

5. Conclusions

This paper examines the stress wave propagation through a single joint and multiple parallel joints under unidirectional in situ stress by laboratory experiments and numerical simulations. The crack propagation of the ZDF-3 electric spark source and blast-induced source in the jointed rock mass are compared. Key findings are summarized below:(1)We investigate the propagation and attenuation characteristics of stress wave, which are generated by a ZDF-3 electric spark on the larger-scale model with multiple cement mortar blocks. The attenuation laws of the transmitted strain wave are similar for a single joint and multiple parallel joints. The peak value of the dynamic strain decreases as the in situ stress increases. Due to the increase of joint number, the dynamic strain of rock mass with multiple parallel joints is significantly less than that with a single joint. The in situ stress has a certain inhibitory effect on the stress wave propagation.(2)The crack propagation experiments in jointed rock masses under unidirectional in situ stress demonstrate the similarity between the ZDF-3 electric spark source and the blast-induced source. Whether the ZDF-3 electric spark excitation or the explosive excitation is used, a blasting crater appears at the bottom of the model. The blasting funnel body separates from the parent body in Model 1 but does not separate in Model 2. The ZDF-3 electric spark can be used as an excitation source in stress wave propagation. The in situ stress has a guiding effect on crack propagation, and an inclined deflection crack is formed on the upper surface of Model 2.(3)The crack propagation in the parent body of Model 2 is analyzed using the dynamic relaxation method of ANSYS/LS-DYNA. The numerical simulation results show that two parallel cracks are generated along the direction of the in situ stress loading, which is consistent with the phenomenon observed in the laboratory experiment. The numerical simulation further verifies that the in situ stress has a guiding effect on the crack propagation.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

X. Li, C. Li, W. Cao, and M. Tao, “Dynamic stress concentration and energy evolution of deep-buried tunnels under blasting loads,” International Journal of Rock Mechanics and Mining Sciences, vol. 104, pp. 131–146, 2018.
View at: Publisher Site | Google Scholar
M. He, “Latest progress of soft rock mechanics and engineering in China,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 6, no. 3, pp. 165–179, 2014.
View at: Publisher Site | Google Scholar
T. Yin, L. Bai, X. Li, X. Li, and S. Zhang, “Effect of thermal treatment on the mode I fracture toughness of granite under dynamic and static coupling load,” Engineering Fracture Mechanics, vol. 199, pp. 143–158, 2018.
View at: Publisher Site | Google Scholar
W. Liang, G. Zhao, H. Wu, and B. Dai, “Risk assessment of rockburst via an extended MABAC method under fuzzy environment,” Tunnelling and Underground Space Technology, vol. 83, pp. 533–544, 2019.
View at: Publisher Site | Google Scholar
D. Huang, J. Wang, and S. Liu, “A comprehensive study on the smooth joint model in DEM simulation of jointed rock masses,” Granular Matter, vol. 17, no. 6, pp. 775–791, 2015.
View at: Publisher Site | Google Scholar
J. C. Li, H. B. Li, G. W. Ma, and J. Zhao, “A time-domain recursive method to analyse transient wave propagation across rock joints,” Geophysical Journal International, vol. 188, no. 2, pp. 631–644, 2012.
View at: Publisher Site | Google Scholar
J. B. Zhu and J. Zhao, “Obliquely incident wave propagation across rock joints with virtual wave source method,” Journal of Applied Geophysics, vol. 88, pp. 23–30, 2013.
View at: Publisher Site | Google Scholar
L. F. Fan and H. Y. Sun, “Seismic wave propagation through an in-situ stressed rock mass,” Journal of Applied Geophysics, vol. 121, pp. 13–20, 2015.
View at: Publisher Site | Google Scholar
W. Wang, H. Hao, X. Li, Z. Yan, and F. Gong, “Effects of a single open joint on energy transmission coefficients of stress waves with different waveforms,” Rock Mechanics and Rock Engineering, vol. 48, no. 5, pp. 2157–2166, 2015.
View at: Publisher Site | Google Scholar
J. Zhao and J. G. Cai, “Transmission of elastic p-waves across single fractures with a nonlinear normal deformational behavior,” Rock Mechanics and Rock Engineering, vol. 34, no. 1, pp. 3–22, 2001.
View at: Publisher Site | Google Scholar
H. Li, Y. Liu, J. Li et al., “Numerical study on oblique incidence across rock masses with linear and nonlinear joints,” Arabian Journal of Geosciences, vol. 9, no. 1, p. 20, 2016.
View at: Publisher Site | Google Scholar
X. L. Huang, S. W. Qi, Y. S. Liu, and Z. F. Zhan, “Stress wave propagation through viscous-elastic jointed rock masses using propagator matrix method (PMM), geophys,” Journal of Intelligence, vol. 200, pp. 452–470, 2015.
View at: Google Scholar
W. Wu, J. B. Zhu, and J. Zhao, “A further study on seismic response of a set of parallel rock fractures filled with viscoelastic materials,” Geophysical Journal International, vol. 192, no. 2, pp. 671–675, 2013.
View at: Publisher Site | Google Scholar
Z.-H. Zhang, J.-H. Deng, and J.-B. Zhu, “A rapid and nondestructive method to determine normal and shear stiffness of a single rock joint based on 1d wave-propagation theory,” Geophysics, vol. 83, no. 1, pp. WA89–WA100, 2018.
View at: Publisher Site | Google Scholar
J. C. Li, W. Wu, H. B. Li, J. B. Zhu, and J. Zhao, “A thin-layer interface model for wave propagation through filled rock joints,” Journal of Applied Geophysics, vol. 91, pp. 31–38, 2013.
View at: Publisher Site | Google Scholar
L. F. Fan and L. N. Y. Wong, “Stress wave transmission across a filled joint with different loading/unloading behavior,” International Journal of Rock Mechanics and Mining Sciences, vol. 60, pp. 227–234, 2013.
View at: Publisher Site | Google Scholar
E. T. Brown and E. Hoek, “Trends in relationships between measured in-situ stresses and depth,” International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts, vol. 15, no. 4, pp. 211–215, 1978.
View at: Publisher Site | Google Scholar
F. V. Donzé, J. Bouchez, and S. A. Magnier, “Modeling fractures in rock blasting,” International Journal of Rock Mechanics and Mining Sciences, vol. 34, no. 8, pp. 1153–1163, 1997.
View at: Publisher Site | Google Scholar
J. B. Zhu, Z. Y. Liao, and C. A. Tang, “Numerical SHPB tests of rocks under combined static and dynamic loading conditions with application to dynamic behavior of rocks under in situ stresses,” Rock Mechanics and Rock Engineering, vol. 49, no. 10, pp. 3935–3946, 2016.
View at: Publisher Site | Google Scholar
M. Sagong, D. Park, J. Yoo, and J. S. Lee, “Experimental and numerical analyses of an opening in a jointed rock mass under biaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 48, no. 7, pp. 1055–1067, 2011.
View at: Publisher Site | Google Scholar
L. X. Xie, W. B. Lu, Q. B. Zhang, Q. H. Jiang, G. H. Wang, and J. Zhao, “Damage evolution mechanisms of rock in deep tunnels induced by cut blasting,” Tunnelling and Underground Space Technology, vol. 58, pp. 257–270, 2016.
View at: Publisher Site | Google Scholar
C. Yi, D. Johansson, and J. Greberg, “Effects of in-situ stresses on the fracturing of rock by blasting,” Computers and Geotechnics, vol. 104, pp. 321–330, 2018.
View at: Publisher Site | Google Scholar
C. He, J. Yang, and Q. Yu, “Laboratory study on the dynamic response of rock under blast loading with active confining pressure,” International Journal of Rock Mechanics and Mining Sciences, vol. 102, pp. 101–108, 2018.
View at: Publisher Site | Google Scholar
J. J. Zhao, Y. Zhang, and P. G. Ranjith, “Numerical modelling of blast-induced fractures in coal masses under high in-situ stresses,” Engineering Fracture Mechanics, vol. 225, Article ID 106749, 2019.
View at: Publisher Site | Google Scholar
X.-Q. Fang, T.-F. Zhang, B.-L. Li, and R.-J. Yuan, “Elastic-slip interface effect on dynamic stress around twin tunnels in soil medium subjected to blast waves,” Computers and Geotechnics, vol. 119, Article ID 103301, 2020.
View at: Publisher Site | Google Scholar
J.-J. Zhao, Y. Zhang, and P. G. Ranjith, “Numerical simulation of blasting-induced fracture expansion in coal masses,” International Journal of Rock Mechanics and Mining Sciences, vol. 100, pp. 28–39, 2017.
View at: Publisher Site | Google Scholar
Z.-L. Wang, Y.-C. Li, and R. F. Shen, “Numerical simulation of tensile damage and blast crater in brittle rock due to underground explosion,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 5, pp. 730–738, 2007.
View at: Publisher Site | Google Scholar
G. W. Ma and X. M. An, “Numerical simulation of blasting-induced rock fractures,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 6, pp. 966–975, 2008.
View at: Publisher Site | Google Scholar
Y. Gou, X. Shi, X. Qiu, J. Zhou, H. Chen, and X. Huo, “Propagation characteristics of blast-induced vibration in parallel jointed rock mass,” International Journal of Geomechanics, vol. 19, no. 5, Article ID 04019025, 2019.
View at: Publisher Site | Google Scholar
Z. Zhu, B. Mohanty, and H. Xie, “Numerical investigation of blasting-induced crack initiation and propagation in rocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 3, pp. 412–424, 2007.
View at: Publisher Site | Google Scholar
M. M. Dehghan Banadaki and B. Mohanty, “Numerical simulation of stress wave induced fractures in rock,” International Journal of Impact Engineering, vol. 40-41, pp. 16–25, 2012.
View at: Publisher Site | Google Scholar
J. D. Chen, W. G. Yu, and S. Lu, “Performance analysis about high-power sparker source,” Journal of Geophysical Equipment, vol. 23, pp. 165–168, 2013, (in Chinese).
View at: Google Scholar
H. R. Nicholls and W. I. Duvall, “Pre-splitting rock in the presence of a static stress field,” Tech. Rep., U.S. Dept. of the Interior, Bureau of Mines, Washington, DC, USA, 1966, Technical Report.
View at: Google Scholar
K. R. Y. Simha and W. L. Fourney, “Studies on explosively driven cracks under confining in-situ stresses: in: rock mechanics in productivity and protection,” in Proceedings of the 25th Symposium on Rock Mechanics, Evanston, IL, USA, June 1984.
View at: Google Scholar
W. L. Fourney, R. D. Dick, X. J. Wang, and Y. Wei, “Fragmentation mechanism in crater blasting,” International Journal of Rock Mechanics and Mining Science & Geomechanics Abstracts, vol. 30, no. 4, pp. 413–429, 1993.
View at: Publisher Site | Google Scholar
S. Aziznejad, K. Esmaieli, J. Hadjigeorgiou, and D. Labrie, “Responses of jointed rock masses subjected to impact loading,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 10, no. 4, pp. 624–634, 2018.
View at: Publisher Site | Google Scholar
H. B. Li, X. Xiang, J. C. Li, J. Zhao, B. Liu, and Y. Q. Liu, “Rock damage control in bedrock blasting excavation for a nuclear power plant,” International Journal of Rock Mechanics and Mining Sciences, vol. 48, pp. 210–218, 2011.
View at: Google Scholar
P. Yan, W.-b. Lu, M. Chen, Y.-g. Hu, C.-b. Zhou, and X.-x. Wu, “Contributions of in-situ stress transient redistribution to blasting excavation damage zone of deep tunnels,” Rock Mechanics and Rock Engineering, vol. 48, no. 2, pp. 715–726, 2015.
View at: Publisher Site | Google Scholar
J. Yang, W. Lu, Y. Hu, M. Chen, and P. Yan, “Numerical simulation of rock mass damage evolution during deep-buried tunnel excavation by drill and blast,” Rock Mechanics and Rock Engineering, vol. 48, no. 5, pp. 2045–2059, 2015.
View at: Publisher Site | Google Scholar
LSTC, LS-DYNA Theoretical Manual, Livermore Software Technology Corporation, Livermore, CA, USA, 2010.
M. Tao, X. Li, and C. Wu, “Characteristics of the unloading process of rocks under high initial stress,” Computers and Geotechnics, vol. 45, pp. 83–92, 2012.
View at: Publisher Site | Google Scholar
S.-Y. Xiao, L.-J. Su, Y.-J. Jiang, and Z.-X. Liu, “Numerical analysis of hard rock blasting unloading effects in high in situ stress fields,” Bulletin of Engineering Geology and the Environment, vol. 78, no. 2, pp. 867–875, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Zhanfeng Fan and Jianhua Cai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies