Abstract

In this paper, the instantaneous rock burst test of Beishan granite is carried out by using a deep rock burst simulation test system and an acoustic emission monitoring system. The acoustic emission data were monitored in real time during the test. The variation of the number and energy of acoustic emission events was studied, and the distribution characteristics of rock burst debris were analyzed. Based on plate and shell mechanics, the failure process of surrounding rock is discussed from the perspective of structural stability. The results show that (1) when the vertical stress reaches 171.31 MPa, the specimen is destroyed and the number of acoustic emission events and cumulative absolute energy before the specimen is destroyed increase sharply. (2) The debris generated by rock burst is mainly composed of slab debris, flaky debris, and thin flaky debris, accounting for 93.53% of the total debris. (3) When the length or height of the rock slab is constant, the maximum tensile stress in the rock slab decreases nonlinearly with the increase of rock slab thickness. For the same size of the rock slab, the farther away from the roadway wall, the greater the maximum tensile stress in the rock slab. (4) When the thickness of the rock slab is constant, the maximum tensile stress in the rock slab increases nonlinearly with the increase of height to thickness ratio K. When the ratio of height to thickness K is constant, the maximum tensile stress in the rock slab increases with the increase of rock slab thickness h. (5) With the increase of covering depth, the critical failure thickness of the rock slab decreases nonlinearly and the surplus energy increases nonlinearly.

1. Introduction

At present, China produces about 2000 tons of high-level radioactive waste every year. The storage of nuclear waste is close to saturation, and the demand for high-level radioactive waste disposal libraries is very urgent. Geological treatment of rock layers in deep underground areas is an internationally recognized method for proper disposal of high-level radioactive waste [13]. Many countries, such as France, Finland, Italy, and South Korea, have started to build underground nuclear waste repositories [411]. Granite has the characteristics of high strength and low permeability, so it has been selected by many countries as a potential bearing rock for the construction of underground nuclear waste repository. In China, the Beishan area of Gansu Province is a large-scale granite underground, which has the advantages of good mechanical properties, stable geological structure, and less biological activity. After years of repeated assessment, the Beishan area has been identified as one of the best preselected sites for China’s nuclear waste repository [12, 13]. However, when drilling exploration in the Beishan area, it is found that there is high ground stress in the local area, such as core cake phenomenon. The nuclear waste repository requires tens of thousands of years of safety, and the local high ground stress area is likely to induce rock burst during excavation unloading, which will not only affect the stability of the underground repository but also may pose a major threat to the construction of the underground repository and the safety of the staff. Therefore, it is necessary to conduct in-depth research on the characteristics of rock burst and the initiation mechanism of Beishan granite.

Rock burst refers to a dynamic failure phenomenon caused by the instantaneous release of strain energy accumulated inside the rock due to excavation unloading under high ground stress conditions. At present, the definition of rock burst has not been unified. Cook [14] gave the definition of rock burst for the first time, that is, rock burst is the uncontrolled destruction of rock mass after the release of violent energy, while Russnes [15] believes that the new fracture surface of rock in the process of failure is accompanied by sound, spalling, spalling ejection, and other phenomena, which is rock burst. As the core of the research on rock burst, the mechanism of rock burst is related to the accuracy of rock burst prediction and the effectiveness of prevention and control measures. Therefore, it has become the focus of scholars at home and abroad. Chen et al. [16], based on uniaxial and triaxial rock tests under thermomechanical coupling, discussed the brittle failure mechanism of granite and evaluated the occurrence trend of rock burst at different temperatures. Based on the assumption of static and dynamic stress coupling, Cai et al. [17] studied the mechanism of rock burst and proposed four quantitative evaluation methods, which provided a method for understanding the mechanism of rock burst and monitoring rock burst. Based on a rock burst accident in a deep mine, Zheng et al. [18] analyzed the crack field, stress field, displacement field, and kinetic energy of roadway surrounding rock under cyclic dynamic load by using a particle discrete element method and discussed the instability mechanism of roadway rock burst roof. Si and Gong [19] conducted triaxial unloading compression tests on granite specimens to investigate the effects of the three-dimensional stress state and unloading rate on the rock burst mechanism. Tan et al. [20] established a mechanical model for energy release in surrounding rock of roadway, analyzed the law of energy release in surrounding rock of roadway, and revealed the mechanism of rock burst caused by kinetic energy release in surrounding rock. Gong et al. [21, 22] conducted true triaxial tests on granite and red sandstone cube samples containing prefabricated holes, reproduced the process of rock burst induced by spalling damage in deep-buried hard rock tunnels, and revealed the mechanism of rock burst induced by spalling damage. Based on the theory of slab deflection and energy principle, Guo et al. [23] established a mechanical model for roadway floor under high-level stress and studied the mechanism of rock burst in the floor of south track roadway in Xing’an Coal Mine.

Acoustic emission (AE) is a physical phenomenon in which part of strain energy is released as transient elastic wave under loading during crack propagation in rock. Because of its advantages of dynamic nondestructive monitoring and real-time analysis, it can record the whole process of rock failure. Therefore, many scholars have studied rock failure based on acoustic emission technology. Meng et al. [24] studied the energy accumulation, evolution, and dissipation characteristics of sandstone specimens during uniaxial cyclic loading and unloading compression at 6 different loading rates. The stress-strain relationship and acoustic emission characteristics of rock specimen deformation and failure are analyzed, and the law of energy evolution in the process of rock deformation and failure is revealed. Zhang et al. [25] conducted triaxial tests on granite specimens at a depth of 1150 m using an acoustic emission monitoring system to study the crack evolution and damage precursor characteristics during the deformation and damage process of deep granite under high peripheral pressure. Akdag et al. [26] studied the time-domain and frequency-domain response of strain rupture by using acoustic emission and kinetic energy analysis methods and quantified the thermal damage evolution process by the b value, cumulative acoustic emission energy, and event rate. Dong et al. [27] conducted tests on soft siltstone specimens under combined dynamic and static loads to quantitatively reveal the mechanical properties and acoustic emission characteristics of soft rock under different dynamic disturbances. Rasskazov et al. [28] proved a new method to solve the problem of local rock burst danger and rock burst prediction through acoustic measurement, spectrum correlation analysis based on acoustic emission signals, and the Griffith–Irwin [29, 30] fracture model. Moradian et al. [31] for a full understanding of the fracture process in brittle rocks tested prismatic specimens of granite containing pre-existing flaws in uniaxial compression tests, and their cracking process was monitored with AE. The characteristics of the AE parameters and the evolution of cracking sequences were analyzed for every cracking level. Based on microcrack and macrocrack damage, a classification of cracking levels is introduced. Ren et al. [32] evaluated the temporal and spatial evolution and damage of microcracks in schist during true triaxial compression and strain rock burst tests using acoustic emission localization techniques and moment tensor analysis. Qin et al. [33] carried out three-point bending tests of sandstone under different spans, analyzed the acoustic emission characteristic parameters of sandstone fracture process, and obtained the damage evolution law of sandstone based on acoustic emission characteristics. Based on triaxial compression and acoustic emission tests, Wang et al. [34] studied the damage characteristics and mechanical properties of Beishan deep granite under different confining pressure levels and discussed the evolution of acoustic emission parameters and strain energy of deep rock at different damage stages.

In the investigation and study of hard rock failure, a large number of scholars have found that a large number of slab cracking failures exist, ranging from shear slip failure of geological faults to small rock specimens in laboratory tests [35, 36]. Slabbing failure, as a general rule and phenomenon of hard and brittle surrounding rock in deep-buried high-stress roadway, is shown as the surrounding rock cut by cracks approximately parallel to the excavation face, forming layers of thin slab approximately parallel to the excavation face [37]. It is generally believed that these rock slabs are the key to rock burst initiation. When the rock slab remains stable, the surrounding rock remains stable, while when the surface rock slab is suddenly unstable, the surrounding rock will spontaneously fail from the surface to the depth. At present, the rock slab is often used to describe the failure phenomenon of rock burst or the initiation process of rock burst qualitatively. This paper intends to quantitatively describe the mechanical characteristics of rock burst initiation based on the failure of the rock slab.

In summary, scholars have conducted in-depth studies on the rock burst mechanism of different research objects, such as thermodynamic coupling, static and dynamic coupling, strength, and energy. However, there has been limited research on the initiation mechanisms of rock burst mechanics using structural stability theory. Extensive engineering practice has shown that slab structure failure is a common failure form of underground engineering surrounding rock. In view of this, this paper carried out an indoor instantaneous rock burst test on Beishan granite; during the test, acoustic emission data were monitored synchronously. The variation law of energy and acoustic emission events during the failure process of Beishan granite was studied, and the distribution characteristics of rock burst debris were analyzed. Based on plate shell mechanics, the initiation mechanisms of rock burst mechanics in Beishan granite were explored from the perspective of structural stability.

2. Test Conditions and Test Scheme

2.1. Test Conditions

Test rock samples were taken from the preselection area of Beishan, Gansu Province. The rocks were tonalite, complete, and fresh, and the samples were processed into cuboid specimens of 150 mm × 60 mm × 30 mm. The six surfaces of the specimen were highly polished to a flatness of 0.0025 mm per surface and a surface roughness of Raw 1.6.

Test equipment adopts the deep rock burst process simulation test system (Figure 1) independently developed by the China University of Mining and Technology (Beijing). The system has a test host, loading system, and data acquisition system. The acquisition system consists of a stress acquisition system, an acoustic emission acquisition system, and a camera system (Figure 2). The system can not only realize conventional tests such as uniaxial tension and compression, triaxial compression, tension and compression, tension, and shear but also realize true triaxial loading single-sided independent sudden unloading tests.

Stress acquisition uses a DSG9803 strain amplifier and a USB8516 portable data acquisition instrument, which can realize 8-channel independent acquisition. Acoustic emission adopts a PCI-II system provided by American Physical Acoustics Company. The acquisition frequency of AE sensors was set to 200 kHz, and the threshold value for signal acquisition was set to 26 dB. The amplification factor for signal acquisition was set to 40. A high-speed photography system can take and store images at a resolution of 1024 × 1024 at 1000 frames/s.

2.2. Test Scheme

The instantaneous rock burst test of Beishan granite specimens was carried out through laboratory tests to simulate the possibility of instantaneous rock burst in surrounding rock under different excavation depths. Because the mechanical properties of Beishan granite are stable and the texture is uniform, we carried out a rock burst simulation test on one specimen. As shown in Figure 3, the stress loading process is divided into three stages: in the first stage, at the beginning of the test, the specimen was initially loaded to the original stress corresponding to a burial depth of 500 m, that is, σH = 15.14 MPa, σh = 10.57 MPa, and  = 13.40 MPa, and held for 5 minutes. According to the Kirsch equation, the excavation stress concentration is calculated; the vertical stress σV is 29.63 MPa after excavation along the minimum horizontal stress, and the whole process is recorded as level 1 loading. In the second stage, due to the small stress at this time, the possibility of instantaneous rock burst is less. To minimize the impact of continuous unloading on the rock properties, no unloading was performed under this stress state. After maintaining one stress state for 5 minutes, the next load is continued to be applied. Each level of loading is equivalent to an increase of 200 m in buried depth, and the whole process is 2 to 4 levels of loading in turn. In the third stage, when the burial depth increases to 1300 m, the original stresses are σH = 32.74 MPa, σh = 23.37 MPa, and  = 34.84 MPa and the vertical stress σV is 81.15 MPa. This stress state had the potential to induce instantaneous rock burst, so instantaneous unloading was performed along the σh direction while keeping σH and σV constant to observe whether rock burst occurred. If rock burst occurs, the test is stopped; otherwise, the next level of load is applied, that is, σH = 37.14 MPa, σh = 26.57 MPa, and σV = 94.03 MPa. Then, it is unloaded instantaneously along the direction of σh, keeping σH and σV unchanged and observing whether rock burst occurs. If rock burst occurs, the test is stopped. If there is no rock burst, the next level of load is applied. The whole process is 5 to 12 levels of loading in turn. When the load reaches level 12, rock burst occurs, and the test is terminated.

3. Test Results and Analysis

3.1. Rock Burst Stress and Acoustic Emission Event Number Analysis

After multiple loading and unloading, when the original stress reaches 72.36 MPa, σH reaches 63.54 MPa, σh reaches 45.77 MPa, and the vertical stress σV is 171.31 MPa, the specimen has rock burst. At the moment of rock burst, the vertical stress σV dropped to 109.6 MPa and the horizontal stress σH dropped to 43.32 MPa. As shown in Figure 4, from the beginning of loading to the occurrence of rock burst, the generation of acoustic emission events is mainly concentrated in the unloading process of stage. During level 11, the number of acoustic emission events increases greatly and 1244 acoustic emission events are generated and with a slight crack sound, indicating that internal cracks of the specimen increase and crack propagation intensifies. During level 12, the number of acoustic emission events increases sharply and 13,855 acoustic emission events are generated in a short time, indicating that a large number of cracks expand, aggregate, and penetrate at this stage, which eventually leads to rock burst of the specimen with a violent noise.

The process of rock burst occurrence is shown in Figure 5. The failure process of the rock sample can be roughly divided into four stages: calm period, small particle ejection, flake stripping accompanied by particle mixing ejection, and recalm period. In the calm stage, as shown in Figure 5(a), the rock sample accumulates certain elastic strain energy through elastic deformation under loading. With the increase of σV, the initial microcracks in the rock sample develop and gradually converge to form larger cracks. Then, they enter the small particle ejection stage, and cracks transition from stable expansion to unstable expansion stage. Affected by material heterogeneity, several local fracture zones randomly form on the surface of the free surface, and the ejection of small particle cuttings occurs, as shown in Figure 5(b). In flake stripping accompanied by particle mixing ejection stage, the rock units in a certain depth range of the free surface are split under the action of compressive tensile cracking of σV, forming lamellar or plate-like rock slabs. With the continuous increase of σV, the rock slab on the surface of the free surface bulges and bends towards the free space, and stability gradually decreases. As shown in Figure 5(c), when σV increases to the peak strength of the rock sample, the rock slabs on the surface of the free surface will bend and break, and the residual elastic strain energy stored in the rock sample will be released sharply, forcing the broken thin rock slabs to eject quickly into free space, thus forming severe rock burst failure. Then, as shown in Figure 5(d), the rock sample returned to the calm period after experiencing severe rock burst.

3.2. Analysis of Acoustic Emission Energy Characteristics during Rock Burst

The change trend of the cumulative absolute energy of acoustic emission is basically the same as that of acoustic emission events; that is, there is an obvious increase when rock burst occurs, but the increase of absolute energy during destruction is more obvious than the number of acoustic emission events (Figure 6). Absolute energy is derived from the integral of the squared voltage signal divided by the reference resistance (10 k ohm) over the duration of the AE waveform packet. With the gradual loading of σV, the absolute acoustic emission energy increases slowly, the cumulative acoustic emission energy per level is small, and the average acoustic emission event energy is relatively low. When the load reaches level 12, the acoustic emission energy increases abruptly and the average acoustic emission energy increases sharply. At the moment of rock burst, multiple high-energy acoustic emission events are generated, indicating that large cracks expand rapidly at this stage, and rock burst occurs along with a violent sound.

3.3. Study on the Characteristics of Rock Burst Debris

The macroscopic failure of rock is the result of continuous initiation, development, expansion, aggregation, and penetration of internal cracks. The debris after rock failure is the direct embodiment of fracture expansion, which reflects the fracture form and splitting process inside the rock. Therefore, it is necessary to analyze the geometric characteristics of rock burst debris. The categorization of rock burst debris is shown in Figure 7. According to the debris size, rock burst debris can be divided into four types: coarse debris, medium debris, fine debris, and microdebris. The particle size of coarse debris is larger than 30 mm, the particle size of medium debris is located in the range of 5–30 mm, the particle size of fine debris is located in the range of 0.075–5.00 mm, and the particle size of microdebris is smaller than 0.075 mm. Because it is difficult to obtain the geometric characteristics of debris with a small particle size, the debris with particle size greater than 5 mm is mainly studied, and the height (H), length (L), and thickness (T) of debris are measured by using a vernier caliper. According to the height-to-thickness ratio of debris, the specimens were divided into four categories: block debris (H/T < 3), slab debris (3 < H/T < 6), flake debris (6 < H/T < 9), and thin flake debris (H/T > 9).

Figure 8 shows the distribution map of the debris size ratio of the specimen. The total amount of debris larger than 5 mm is 279, the maximum value of H/T is 16.92, the minimum value is 1.66, and the average value is 6.20. The maximum value of L/T is 10.16, the minimum value is 1.05, and the average value is 4.02. So, the characteristic size of H/L/T is 6.20 : 4.02 : 1. The number of massive debris, tabular debris, flake debris, and thin flake debris is divided into 18, 135, 91, and 34, accounting for 6.47%, 48.56%, 32.73%, and 12.23% of the total. The debris is mainly composed of slab debris, flake debris, and thin flake debris, accounting for 93.53% of the total debris. For slab debris, flake debris, and thin flake debris, the maximum value of H/L is 16.92, the minimum value is 3, and the average value is 6.46; the maximum value of L/T is 10.16, the minimum value is 3.03, and the average value is 5.01; and the characteristic size of H/L/T is 6.45 : 5.01 : 1. Considering the destruction of numerous sheet structures during the rock burst process, it is inferred that a substantial number of rock slabs had likely formed within the specimen’s interior before rock burst occurred. The disruption of these rock slabs is highly likely to have triggered rock burst occurrence.

4. Study on Rock Burst Initiation Mechanism of Beishan Granite

4.1. Establishment and Solution of the Mechanical Calculation Model of Slabbing Failure

In some typical deep hard rock projects at home and abroad, many scholars have observed tensile slabbing failure parallel to the excavation face in the engineering site (Figure 9). In addition, from the distribution of rock burst debris, it can be seen that the debris generated by the instantaneous rock burst test is also dominated by sheet structure. Considering that many slab structures are broken during the rock burst process, it is speculated that a large number of slab structures are formed inside the specimen before rock burst occurs. There is a strong correlation and essential connection between this slabbing phenomenon and rock burst. It can usually be considered as a precursory feature of strain-type rock burst in hard rock. Therefore, rock burst is likely to be caused by the destruction of these slab structures.

Assuming that a series of rock slab structures are formed inside the surrounding rock of the Beishan underground nuclear waste disposal repository, the stress distribution of the surrounding rock is shown in Figure 10. The surrounding rock forms a plastic bearing zone within a certain radius, and in the plastic bearing zone, both radial and tangential compressive stresses increase with the increase of radius and the tangential compressive stress reaches the maximum value at the junction of the plastic bearing zone and the elastic bearing zone and then decreases with the increase of radius in the elastic bearing region, while the radial compressive stress keeps increasing with the increase of radius. Based on the stress distribution in the surrounding rock, it can be observed that both the tangential compressive stress and radial compressive stress encourage the rock slab structures to bend inward towards the interior of the roadway, resulting in unstable failure with ejection into the roadway’s interior.

For ease of analysis and research, considering the effects of blasting and stress concentration, it is assumed that the plastic bearing zone radius of the roadway’s surrounding rock is 1 m. Both the tangential compressive stress and radial compressive stress increase linearly. At the junction of the elastic bearing zone and the plastic bearing zone, the concentration factor for tangential stress is 2, while the concentration factor for radial stress is 0.5. After roadway excavation, it can be assumed that the surrounding rock remains unchanged in strain along the direction of roadway advancement. Under the stress condition of a 500 m underground laboratory,  = 13.4 MPa, σH = 15.14 MPa, and σh = 10.57 MPa, the analysis is carried out according to the plane stress-strain problem. The stress change equation in the plastic bearing zone is shown as follows:

In the equation, σ1 is the tangential normal stress, σ2 is the normal stress along the roadway, σ3 is the radial stress, and d is the distance from the side of the roadway.

The stress condition and the calculation model of the rock slab formed by the surrounding rock are shown in Figure 11. If the length of the rock slab in the x direction is a, the length in the y direction is b, and the thickness is h, then

In the equation, Mx and Fx are the bending moment and pressure equivalent to the eccentric distribution stress on the upper and lower boundaries of the rock slab, My and Fy are the bending moment and pressure equivalent to the eccentric distribution stress on the left and right boundaries of the rock slab, and q is the lateral load difference on both sides of the rock slab.

For a four-sided simply supported thin slab, the boundary conditions are

The displacement satisfying the above boundary conditions can be taken as

The Rayleigh–Ritz method is used to calculate the deformation and stress of the rock slab. The increase of the deformation potential energy Vε is

In the equation, D is the bending stiffness of the rock slab, E is the elastic modulus of the rock slab, and μ is Poisson’s ratio.

The work done by the external load is

Substituting FTx = -Fx, FTy = -Fy, and q(x, y) =q into (7), we obtain

On the upper and lower boundaries of the rock slab,

According to (5)–(8), the total energy of the system is

According to equations (5)–(8), the total energy of the system is

For the equilibrium system, the partial derivative of the total potential energy to Amn is 0:

From (11) and (12),

Substituting (13) into (4) yields

The stress of each point in the rock slab is

Under the combined action of vertical and horizontal loads, along with lateral loads, rock slabs may experience two possible scenarios: tensile failure and compressive failure. According to (14) and (15), it can be seen that σx and σy take extreme values at (a/2, b/2, h/2) and (a/2, b/2, −h/2); that is, at the midpoint of the slab surface, the two sides of the slab will produce the maximum compressive stress and the maximum tensile stress. Since the compressive strength of granite is much greater than its tensile strength, the tensile stress is taken as the judgment basis in this paper, combined with the strength theory. If the following conditions are met, the rock slab will be destabilized and fail:

In the equation, [σt] is the stress corresponding to when the rock slab reaches tensile strength and σt,max is the maximum tensile stress in the rock slab.

It is very important to consider the size effect when simulating instantaneous rock burst in the laboratory. The size of the test sample will affect the rock burst effect, dynamic response, stress and energy release, and other parameters. Selecting the appropriate sample size can simulate and analyze the characteristics and process of instantaneous rock burst more accurately and provide reliable test data and basis for the research and prevention of rock burst phenomenon. Since the assumption of a “thin plate” in elastic-plastic mechanics is followed when solving the force model of the rock slab, that is, the thickness (h) of the slab is less than 1/8∼1/5 of the minimum size (a or b) in the middle plane. Therefore, in the following analysis of the influence of the rock slab shape such as the height, length, and height-to-thickness ratio on the internal stress of the rock slab, as well as the variation of the critical failure thickness and surplus energy with the buried depth of the rock slab, we also adopted the above assumption.

4.2. Analysis of the Influence of Rock Slab Height on the Internal Stress of Rock Slab

Assuming that the maximum thickness of the rock slab is 0.2 m and the length of the rock slab is 1 m, according to (15), the variation trend of the maximum tensile stress within the rock slab with respect to the thickness of the rock slab can be calculated for rock slab heights of 0.6 m, 0.8 m, 1.0 m, 1.2 m, and 1.4 m. The rock slabs with d = 0, d = 0.5 m, and d = 1.0 m are selected to analyze the influence of the distance from the roadway wall on the stability of the rock slab (Figure 12). When the length of the rock slab is the same, the maximum tensile stress in the rock slab decreases significantly nonlinearly with the increase of the thickness of the rock slab, and the larger the height of the rock slab, the greater the maximum tensile stress in the rock slab. In addition, the maximum tensile stress in the rock slab increases with the increase of the distance from the roadway wall.

4.3. Analysis of the Influence of Rock Slab Length on the Internal Stress of Rock Slab

Assuming that the maximum thickness of the rock slab is 0.2 m and the height of the rock slab is 1 m, according to (15), we can calculate the trend of the maximum tensile stress within the rock slab with respect to the thickness of the rock slab for rock slab lengths of 0.6 m, 0.8 m, 1.0 m, 1.2 m, and 1.4 m. The rock slabs with d = 0, d = 0.5 m, and d = 1.0 m are selected to analyze the influence of the distance from the roadway wall on the stability of the rock slab (Figure 13). It can be seen from Figure 13 that the variation trend of tensile stress inside the rock slab is the same as that in Figure 12. When the height of the rock slab is the same, it can be observed that the maximum tensile stress within the rock slab exhibits a significant nonlinear reduction as the thickness of the rock slab increases. Moreover, as the length of the rock slab increases, the maximum tensile stress within the rock slab also increases. In addition, for rock slabs of the same size, the maximum tensile stress within the rock slab increases with the distance from the roadway wall.

4.4. Analysis of the Influence of the Height-Thickness Ratio on the Internal Stress of Rock Slab

The stress in the rock slab under different height-to-thickness ratios is calculated according to (15). The rock slabs with a length of 1.0 m and a thickness of 0.05 m, 0.1 m, and 0.15 m were selected to compare and analyze the change trend of stress. In addition, according to the assumption of elastic mechanics on thin slabs, rock slabs with height-to-thickness ratios K equal to 6, 7, 8, 9, and 10 are selected for stress analysis. The results are shown in Figure 14. When the thickness of the rock slab remains constant, it is evident that the maximum tensile stress within the rock slab nonlinearly increases with an increase in the height-to-thickness ratio K. Furthermore, for rock slabs of the same size, as the distance from the roadway wall decreases, the maximum tensile stress within the rock slab decreases, while it increases as the distance from the roadway wall increases. When the height-to-thickness ratio K is constant, as the thickness h of the rock slab increases, the maximum tensile stress within the rock slab also increases. Moreover, with increasing rock slab thickness h, the rate of stress increase becomes less steep. This occurs because thicker rock slabs have greater load-bearing capacity, experience less deformation, and have reduced bending stress’s influence on overall stress.

4.5. The Change in Critical Failure Thickness and Surplus Energy with the Depth of Rock Slab

The average characteristic size (a/b/h = 6.45 : 5.01 : 1) of the rock slab is selected to analyze the change in tensile stress with the thickness of the rock slab. As shown in Figure 15, with the increase of the thickness, the maximum tensile stress increases. When the thickness of the slab reaches a certain value, the maximum tensile stress may reach the tensile strength of the rock, leading to the failure of the rock slab. This thickness is defined as the critical failure thickness hcr. As the thickness of the rock slab is greater than hcr, the rock slab is destroyed; otherwise, the rock slab remains stable. According to (14) and (15), hcr is calculated with different depths (Figure 16). Further analysis shows that with the increase of depth, the critical failure thickness hcr of the rock slab decreases from 0.1237 m to 0.0547 m, showing a nonlinear decrease, which indicates that in the case of a certain rock slab size ratio, the deeper the rock slab is, the more likely it is to be damaged.

Energy dissipation in the process of rock burst can reflect the intensity of rock burst. It is assumed that a series of rock slabs with a thickness of 0.15 m are formed in the surrounding rock (a/b/h = 6.45 : 5.01 : 1), and its energy is E0. The energy corresponding to the critical failure thickness of the rock slab is defined as the critical failure energy Ecr. When E0 is greater than Ecr, the rock slab will be destroyed, and furthermore, the difference between input energy and critical energy is defined as surplus energy. Surplus energy will be released in other forms (such as kinetic energy). The higher the surplus energy, the more severe the damage. Taking the energy E0 as the reference object, the variation trend of the surplus energy of the rock slab at different depths is analyzed. As shown in Figure 16, with the increase of depth, the surplus energy increases from 8.61 J to 78.52 J, showing a nonlinear increase, which indicates that the same rock slab structure is more and more likely to be destroyed with the increase of depth, and the degree of damage is more and more serious.

5. Conclusions

In this paper, we carried out an instantaneous rock burst test on the Beishan granite using the deep rock burst simulation test system. The variation law of the number of acoustic emission events and acoustic emission energy during the test was studied, the distribution characteristics of rock burst debris were analyzed, and the damage process of the surrounding rock from the perspective of structural stability was further explored preliminarily based on plate and shell mechanics. The research results can provide theoretical references for the construction and structural optimization of the Beishan underground nuclear waste repository and similar deep underground projects around the world. In addition, because the underground nuclear waste repository is still under construction, the stress distribution of the surrounding rock cannot be measured; therefore, when analyzing the damage process of the surrounding rock through the proposed model, the selection of parameters such as the radius of the plastic bearing zone and the stress concentration coefficient is different from the real stress of the surrounding rock of the underground experiments. Specifically, the following conclusions are obtained:(1)When the vertical stress reaches 171.31 MPa, the specimen is destroyed, and the number of acoustic emission events and cumulative absolute energy before the specimen is destroyed increase sharply.(2)The debris generated by rock burst is mainly composed of slab debris, flaky debris, and thin flaky debris, accounting for 93.53% of the total debris.(3)When the length or height of the rock slab is constant, the maximum tensile stress in the rock slab decreases nonlinearly with the increase of rock slab thickness. For the same size of the rock slab, the farther away from the roadway wall, the greater the maximum tensile stress in the rock slab.(4)When the thickness of the rock slab is constant, the maximum tensile stress in the rock slab increases nonlinearly with the increase of the height-to-thickness ratio K. When the ratio of height to thickness K is constant, the maximum tensile stress in the rock slab increases with the increase of rock slab thickness h.(5)With the increase of covering depth, the critical failure thickness of the rock slab decreases nonlinearly, and the surplus energy increases nonlinearly.

Data Availability

The data used to support the findings of the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This article was funded by the National Natural Science Foundation of China Youth Foundation Project (grant no. 52104082) and the National Natural Science Foundation of China Key Foundation Project (Regional Joint Fund) (grant no. U20A20266).