Factors influencing the anti-impact performance of a “roadway rock support” system

A mechanical model of a circular section of a tunnel roadway considering damage is established to improve the impact protection performance of the “roadway rock-support” system and provide a theoretical basis for designing coal mine impact ground pressure roadway support. The formula of the critical rock burst load of a circular roadway is derived according to the instability theory of rock burst disturbance response. The influence of mechanical properties of surrounding rock and roadway support strength on the critical rock burst load of a roadway rock-support system is studied using the control variable method. The research shows that 1) under the support condition of a roadway, with the increase of uniaxial compressive strength, softening modulus, and internal friction angle of surrounding rock, the critical rock burst load of a roadway has an increasing trend; the critical rock burst load of roadway decreases with the increase of the elastic modulus of the surrounding rock. 2) Under the condition of no support, with the increase of uniaxial compressive strength, elastic modulus, and internal friction angle of surrounding rock, the critical rock burst load of roadway tends to increase. With the decrease of the surrounding rock’s softening modulus, the critical rock burst load of the roadway decreases. When the aforementioned four kinds of surrounding rock influence factors are the same, the critical load of rock burst under the supporting condition is much larger than that under the non-supporting condition. 3) A new impact tendency index K is defined as the ratio of the softening modulus and the elastic modulus. 4) The critical load of rock burst increases approximately linearly with the increase of support stress. The critical load of a rock burst is about 400 times as large as the supporting stress. Increasing roadway support strength can greatly improve the stability of the supporting and roadway surrounding rock system, and the stability of the supporting and roadway surrounding rock system can be improved by reasonably changing the mechanical properties of the surrounding rock.

important indexes to characterize the problem of a rock burst when the roadway is disturbed in a specific stress environment. It is also the theoretical basis of stress control theory in rock burst prevention theory.
Scholars have studied and advanced the understanding of the critical load of rock bursts in roadways (Pan, 1999;Pan et al., 2014). Li et al. examined the rock burst instability theory and concluded that the modulus ratio was an important parameter affecting the critical load of a rock burst . Ma et al. used the roof shear beam model to deduce that the critical load decreased with an increase in the widthheight ratio of a rectangular section roadway but increased with the increase of roof thickness. The critical load decreased with an increase in the stiffness ratio of the coal seam to the roof but increased with the increase in the modulus ratio (Ma et al., 2011). Several studies found that the internal friction angle had a relatively strong influence on the damage radius, stress, and displacement of the surrounding rock, while the influence of cohesion is small (Pan and Wang, 2004a;Pan and Wang, 2004b;Pan et al., 2006;Pan et al., 2007). Guo et al. proposed that the rock burst occurrence conditions were related to the ratio of the ascending and descending modulus of rock mass and the development degree of rock fracture (Guo et al., 2011). Chen et al. conducted a theoretical analysis on the critical rock burst of a circular section of a roadway and concluded that the critical rock burst increases with the increase of unified strength theoretical parameters, supporting force, modulus ratio, and internal friction angle and decreases with the increase of dilatancy parameters (Chen et al., 2019;Chen et al., 2020). Meng et al. concluded that part of the input energy in the surrounding rock of a deep circular roadway is released in the form of work conducted by surrounding rock pressure and is used to resist support pressure. When there is no support, this part of the energy is converted into surrounding rock kinetic energy (Meng et al., 2014). Several studies proposed the shape method of the airfoil arch section roadway to improve the bearing capacity and optimize the support parameters (Xiao et al., 2016;Xiao et al., 2017;Ding et al., 2020). Wang et al. used the theory of elastic-plastic mechanics, defined the ratio of pre-peak plastic strain to elastic strain as the pre-peak damage coefficient (K), and established the constitutive relationship of coal and rock based on the Scott model. The results show that the smaller the K value is, the greater the impact tendency is . Yin et al. found that the lateral pressure coefficient, roadway radius, and mechanical parameters of surrounding rock (cohesion and internal friction angle) were the main factors affecting the shape and size of the plastic zone in a circular roadway (Yin et al., 2018a;Yin et al., 2018b). Huang and Gao proposed that the critical condition of a roadway rock burst was closely related to changes in the surrounding rock support mode, initial strength parameters, and strength parameters of the plastic zone. When the critical plastic zone depth and critical load are small, the rock burst occurs easily, the occurrence frequency is high, and the damage is small. When the critical plastic zone depth and critical load are large, a rock burst does not occur easily, and its occurrence frequency is low. Once it occurs, its strength is destructive (Huang and Gao, 2001). Liu et al. analyzed the airfoil tension crack at the preexisting crack tip in the roadway coal wall, determined the critical stress of rock burst, and analyzed its influencing factors (Liu et al., 2018;Liu et al., 2019;Liu et al., 2020;Liu et al., 2021). Several studies used acoustic emission and infrared monitoring to conduct experimental research on rock failures, which was of great significance in detecting rock mass failure on the engineering scale (Dong et al., 2019;Dong et al., 2020a;Dong et al., 2020b;Dong et al., 2021;Dong and Luo, 2022).
Due to the initial stress state of the surrounding rock before roadway excavation, the stress redistribution after excavation reaches the secondary stress state. Therefore, the surrounding rock is not only the main load source but also an important part of the bearing structure. An analytical solution is the most powerful means to analyze such roadway engineering problems. Li et al. believed that once the stress in a certain area of the surrounding rock exceeds the strength of the rock mass, the surrounding rock in this part will enter the plastic or failure state. The elastic solutions of a circular tunnel under a nonaxisymmetric external load and radial and shear internal load were derived. Using these solutions, the new elastic-plastic solution and plastic zone radius equation of a circular tunnel under non-axisymmetric external load and radial internal load were derived (Li et al., 2022). Guo et al. thought that the form of a plastic zone of roadway surrounding rock determines its failure form and degree. Based on the implicit equation of a plastic zone boundary, the variation law of the general form of the plastic zone of rock surrounding a circular roadway was comprehensively studied (Guo et al., 2016;Guo et al., 2018).
In conclusion, the critical load of roadway rock burst is one of the main research topics of tunnel anti-shock support design, but the influence trend of the properties of the surrounding rock on the critical load of a rock burst is not clear. This article considers the small disturbance and three-dimensional equivalent strain caused by the external load of the surrounding rock system in the equilibrium state and adopts the bilinear constitutive relationship. In other words, the overall stress is characterized by a linear climb before the peak and a rapid decline after the peak, which is mainly reflected in the elastic stage and the softening stage after the peak. Based on the Coulomb yield criterion and the rock burst initiation criterion, the critical load formulas with and without support are obtained. According to the analytical solution of the critical rock burst load of a circular roadway, the influence of surrounding rock properties such as uniaxial compressive strength, softening modulus, elastic modulus, and internal friction angle on the critical rock burst load of a roadway was studied by controlling single and double variables. FLAC3D software was used for numerical simulation, and the theoretical and numerical solutions were compared and analyzed to provide a theoretical basis for analyzing roadway anti-impact engineering problems.

Formula derivation of the critical rock burst load of the circular roadway
Assuming that the coal-rock mass system is subjected to external force P (Pa), the characteristic depth of the plastic softening deformation zone is ρ (m). If the coal-rock mass system in the equilibrium state produces a small disturbance ΔP due to the external load, the corresponding plastic softening zone of the coal-rock mass is ρ + Δρ. If the system has a stable equilibrium state, the disturbance ΔP must be bounded. If the equilibrium state of the coal-rock mass system is unstable, even in the case of a small disturbance, it will lead to the infinite expansion of the plastic softening zone. The discriminant relationship of the extreme point of the disturbance response when the rock burst occurs is derived: Δρ ΔP dρ dP ∞ .
( 1 ) Assuming that the roadway is a circular cavity with radius a (m), under the action of hydrostatic pressure P (Pa), the softening zone is generated; that is, the depth of the damage zone is ρ (m), and the inner surface of the roadway is subjected to a constant support stress P 0 (Pa). The roadway is simplified as an axisymmetric plane strain problem, as shown in Figure 1. We define a polar coordinate system r, θ through the roadway center. The Z-axis perpendicular to r, θ plane, takes unit thickness in this study. The radial stress and circumferential stress at the interface between the elastic zone and softening zone are σ p r (Pa) and σ p θ (Pa), and the stress at the interface can be obtained according to the stress solution in the elastic zone.
The coal rock at the junction of softening zone and the elastic zone has just entered the yield state, so the damage variable D 0 is obtained. According to the Coulomb yield criterion, In Eq. 4, σ c (Pa) is the uniaxial compressive strength of coal and rock; in the limit equilibrium of rock failure, the angle between the resultant force formed by the normal stress and the internal friction on the shear plane and the normal stress is called the internal friction angle (Pan and Dai, 2021). It is expressed by φ (°), which reflects the size of the internal friction of the rock, and the internal friction angle variable q 1+sin φ 1−sin φ is defined. According to the deformation geometry equation, In Eq. 5, B is the integral constant, u(m) is radial displacement, ε r is the radial strain component, and ε θ is the circumferential strain component. According to the geometric relationship and the incompressible volume property, the equivalent effect in the softening region becomes The von Mises yield criterion indicates that under certain deformation conditions when the equivalent effect force at a point within a stressed object reaches a certain value, that point begins to enter a plastic state. The equivalent effect variation is a physical quantity used to determine the location of the yield surface of a material after strengthening and is calculated using the same formula as the fourth strength theory formula for calculating the equivalent effect force, except that the stress is changed to a strain. For coal and rock materials, it is considered that the strain ε is equal to that under the uniaxial condition, and the equivalent strain ε under the three-dimensional condition is In Eq. 7, ε r , ε θ , and ε Z are the positive strain component and γ rθ , γ θZ , and γ Zr are the shear strain component.
In the area where the elastic zone and the failure damage zone intersect, the coal rock reaches peak strength and peak strain. Based on the plastic mechanics theory, the uniaxial strain is transformed into the triaxial strain, and the equivalent strain ε ε c is obtained. Let r ρ, and In Eq. 8, ε c is the peak strain. Without considering the volume force, the equilibrium equation in the softening zone is Based on the equilibrium equation and the yield condition, the stress partial differential equation of the damage zone is In Eq. 10, λ (Pa) is the softening modulus. By solving this equation, the stress is In Eq. 11, Q is the integral constant. Boundary conditions on the roadway surface are as follows: In the area where the elastic zone and the failure damage zone intersect, the continuous condition of radial stress is According to Eq. 11, combined with boundary conditions (Chen et al., 2019;Chen et al., 2020), Frontiers in Earth Science frontiersin.org By solving the aforementioned formula dρ dp , it is concluded that the initiation criterion of a rock burst is In Eq. 15, ρ* (m) is the critical softening zone depth of the roadway when a rock burst occurs; E (Pa) is the elastic modulus of coal and rock materials. Before the triaxial compression stress-strain curve of coal and rock reaches the peak, the coal and rock are in the elastic deformation stage, and the elastic modulus is E. When the curve reaches the peak, it is simplified to a linear deformation. The softening modulus is defined as λ (Pa), and its expression is λ ?σ f − σ ′ r ). ε f is the strain corresponding to the peak strength σ f (Pa) of the stress-strain curve under triaxial compression; σ ′ r (Pa) and ε ′ r are the termination strength and termination strain at the end of the triaxial compression test, respectively.
When support stress is not considered: Corresponding to Eq. (15), it is deduced that the critical load P* (Pa) of a rock burst when the section shape of the roadway is circular is equal to When support stress is not considered: 3 Analysis of factors that influence the critical rock burst load in a circular roadway According to Eq. (18), the critical rock burst load of a circular roadway is related to the mechanical properties of the surrounding rock, including uniaxial compressive strength σ c , softening modulus λ, elastic modulus E, internal friction angle φ, and support stress P 0 . According to the Mining Engineering Design Manual (Zhang, 2003), the rock grades of the surrounding rock are soft, medium-hard, and hard. The rock mass test parameters are as follows: the uniaxial compressive strength is 95.6 MPa, the softening modulus is 25 GPa, the elastic modulus is 24 GPa, and the internal friction angle is 30°. The mechanical properties of surrounding rock, such as uniaxial compressive strength, softening modulus, elastic modulus, internal friction angle, and the influence of support stress on the critical load of roadway rock burst, were studied by the control variable method.
The control variable method means that when a physical quantity is affected and restricted by multiple factors at the same time, the multi-factor problem is transformed into multiple single-factor problems. When studying the relationship between this physical quantity and one of the factors, only one factor is allowed to change, and the other factors are controlled to remain unchanged to determine the relationship between the relevant physical quantities.

Effect of the uniaxial compressive strength on the critical load
The roadway surrounding rock softening modulus is 25 GPa, the elastic modulus is 24.9 GPa, the internal friction angle is 30°, the uniaxial compressive strength is 36.7-231.7 MPa, and roadway support stress is 5 MPa. These parameters are used in Eqs. 17 and 18. Figure 2A shows the relationship curve between uniaxial compressive strength and critical load under supporting conditions; Figure 2B represents the relationship curve between uniaxial compressive strength and critical load under supporting conditions. It can be seen from Figure 2 that when the support stress is 5 MPa and the uniaxial compressive strength is 36.7 and 231.7 MPa, the critical loads of rock burst are 782.1 and 4884 Mpa, respectively (Kang et al., 2015). The critical load of rock burst increases with increasing uniaxial compressive strength. The fitting curve equation is y 21x + 12 (Li and Wu, 2019). Without support, when the uniaxial compressive strength is 36.7 and 231.7 MPa, the critical loads of rock burst are 36.63 and 231.2 MPa, indicating that the critical load of rock burst increases with the increase of uniaxial compressive strength, and the critical load of rock burst is approximately equal to the uniaxial compressive strength. The fitting curve equation is y x (Pan, 1999). Under the condition of support and no support, the critical load of rock burst increases with the increase of uniaxial compressive strength.

Effect of the softening modulus on the critical load
The uniaxial compressive strength of the surrounding rock is 95.6 MPa, the elastic modulus is 24.9 GPa, the internal friction angle is 30°, the softening modulus is 10-40 GPa, and the roadway support stress is 5 MPa. When these parameters are used in Eqs. 17 and 18, Figure 3A represents the relationship curve between the softening modulus and the critical load under supporting conditions, and Figure 3B represents the relationship curve between the softening modulus and the critical load without supporting conditions. It can be seen from Figure 3 that when the support stress is 5 MPa and the softening modulus is 10 and 40 GPa, the critical loads of rock burst are 1524 and 2577 MPa, respectively (Kang et al., 2015). The critical load of rock burst increases with the increase of softening modulus. The fitting curve equation is y 1.2 × 10 −7 x 2 + 0.029x + 1.2 × 10 3 (Li and Wu, 2019). Without support, when the softening modulus is 10 and 40 GPa, the critical loads of rock bursts are 166.8 and 77.56 Mpa, respectively, which shows that the critical load of the rock burst decreases with the increase of softening modulus. The fitting curve equation is y −5.1 × 10 −12 x 3 + 4.9 × 10 −7 x 2 − 0.017x + 2.9 × 10 2 (Pan, 1999). Under the conditions of support and no support, the Frontiers in Earth Science frontiersin.org softening modulus of the roadway has the opposite influence of law on the critical load of the rock burst.

Effect of the elastic modulus on the critical load
The uniaxial compressive strength of the roadway surrounding rock is 95.6 MPa, the softening modulus is 25 GPa, the internal friction angle is 30°, the elastic modulus is 5.2-79.8 GPa, and the roadway support stress is 5 MPa. When these parameters are used in Eqs. 17 and 18, Figure 4A represents the relationship curve between the elastic modulus and the critical load under supporting conditions, and Figure 4B represents the relationship curve between the elastic modulus and the critical load without supporting conditions. It can be seen from Figure 4 that under the condition of 5 MPa support stress, with the increase of elastic modulus, the critical rock burst load has experienced three stages of rapid decrease, slow decrease, and basic stability (Kang et al., 2015). The elastic modulus is in the stage of rapid decrease of the critical load from 5 to 10 GPa, and the slope of the curve is −0.5. The elastic modulus is in the stage of a slow decrease of the critical load from 10 to 50 GPa, and the slope of the curve is −0.0375. The elastic modulus of 50-79.8 GPa is the stable stage of the critical load. The fitting curve equation is y −3.4 × 10 −20 x 5 + 8.3 × 10 −15 x 4 − 7.7 × 10 −10 x 3 + 3.4 × 10 −5 x 2 − 0.74x + 8.3 × 10 3 (Li and Wu, 2019). Without support, when the elastic modulus is 20 and 79.8 GPa, the critical loads of a rock burst are 86.1 and 199.5 MPa, which shows that the critical load of a rock burst increases with the increase of elastic modulus. The fitting curve equation is y −1.6 × 10 −23 x 2 − 0.0019x + 48.

Effect of the internal friction angle on the critical load
The uniaxial compressive strength of the roadway surrounding rock is 95.6 MPa, the softening modulus is 25 GPa, the elastic modulus is 24.9 GPa, the internal friction angle is 20°-40°, and the roadway support stress is 5 MPa. When these parameters are used in Eqs. 17 and 18, Figure 5A represents the relationship curve between the internal friction angle and the critical load under supporting conditions, and Figure 5B represents the relationship curve between the internal friction angle and the critical load without supporting conditions. It can be seen from Figure 5 that under the condition of support stress of 5 MPa, the critical load of rock burst increases with the internal friction angle, and it has experienced three different change trend stages of 20°-30°, 30°-35°, and 35°-40° (Kang et al., 2015). In particular, when the internal friction angle is greater than 35°, the critical load of rock burst increases rapidly with the increase of the internal friction angle (Li and Wu, 2019). Without support, when

Effect of the support stress on the critical load
The uniaxial compressive strength of the roadway surrounding rock is 95.6 MPa, the softening modulus is 25 GPa, the elastic modulus is 24.9 GPa, the internal friction angle is 30°, and the roadway support stress is 0-10 MPa. When these parameters are used in Eq. 17, the relationship curve between the roadway support stress and the critical rock burst load of a circular roadway is obtained, as shown in Figure 6. According to Figure 6, when the support stress is 0 and 10 MPa, the critical loads of the rock burst are 95.41 and 3947 Mpa, respectively. The critical load of a rock burst increases with the increase of support stress, the slope of the curve is 3.9 × 10 2 , approximately, and the critical load of a rock burst is about 400 times higher than the support stress. The fitting curve equation is y 3.9 × 10 2 x + 95. It can be seen that the stability of the support and the roadway surrounding the rock system can be greatly improved by increasing the support strength.

Effect of K on the critical load
According to the analysis provided in Sections 3.1-3.5, under the condition of support, the critical load of a roadway rock burst increases with the increase of the uniaxial compressive strength, softening modulus, and internal friction angle of the surrounding rock and decreases with the increase of elastic modulus of the surrounding rock. Under the condition

Frontiers in Earth Science
frontiersin.org of no support, the critical rock burst load of the roadway increases with the increase of uniaxial compressive strength, elastic modulus, and internal friction angle of the surrounding rock, and decreases with the decrease of the surrounding rock softening modulus. Under the conditions of support and non-support, the softening modulus and elastic modulus have opposite effects on the critical rock burst load of the roadway. Therefore, to better obtain the influence law of the softening modulus and elastic modulus on the critical load of a roadway rock burst, we define a new index K of rock burst tendency as the ratio of the softening modulus and the elastic modulus: K λ E . The uniaxial compressive strength of the roadway surrounding rock is 95.6 MPa, the internal friction angle is 30°, the roadway support stress is 5 MPa, and K is 0.1-10. When these parameters are used in Eqs. 17 and 18, Figure 7A represents the relationship curve between the new index K of impact tendency and the critical load under supporting conditions, and Figure 7B represents the relationship curve between the new index K of impact tendency and the critical load without supporting conditions. It can be seen from Figure 7 that 1) when the support stress of the roadway is 5 MPa, and the new indexes of impact tendency K are 0.1 and 10, the critical loads of a rock burst are 1627 and 87000 MPa, indicating that the critical load of a rock burst increases with the increase of the new index of impact tendency K. The fitting curve equation is y −2.9x 3 + 53x 2 + 6.6 × 10 2 x + 1.5 × 10 3 . 2) Without support, with the increase of K, the critical load of a rock burst experiences two stages of rapid decrease and slow decrease. When K is between 0.1 and 1, the critical load decreases rapidly from 525.8 to 91.25 MPa; the critical load decreases slowly from 91.25 to 52.53 MPa when K is between 1 and 10.

Effect of double variables on the critical load of the rock burst
To study the influence of double variables on the critical load P*, consider the conditions when the uniaxial compressive strength σ c is 95.6 MPa, the softening modulus λ is 25 GPa, the elastic modulus E is 24.9 GPa, and the internal friction angle φ is 30°. The support stress P 0 is 5 MPa. According to Figure 8A, when the softening modulus and uniaxial compressive strength change, and other parameters remain unchanged, the critical load increases with the increase of softening modulus and uniaxial compressive strength. According to Figure 8B, when the elastic modulus and uniaxial compressive strength change, and other parameters remain unchanged, the critical load increases with the decrease of elastic modulus and the increase of uniaxial compressive strength. According to Figure 8C, when the support stress and uniaxial compressive strength change, and other parameters remain unchanged, the critical load increases with the increase of support stress and uniaxial compressive strength. Figure 8D shows that when the elastic modulus and softening modulus change, and other parameters remain unchanged, the critical load increases with the decrease of elastic modulus and the increase of softening modulus. It can be seen from Figure 8E that when the support stress and softening modulus change, and other parameters remain unchanged, the critical load increases with the increase of support stress and softening modulus. According to Figure 8F, when the elastic modulus and internal friction angle change, and other parameters remain unchanged, the critical load increases with the decrease of elastic modulus and the increase of internal friction angle.
4 Numerical solution of the uniaxial compressive strength, elastic modulus, and internal friction angle According to the results of the analysis of the factors affecting the critical load of the rock burst in Section 3, the uniaxial compressive strength, elastic modulus, and internal friction angle are selected for numerical simulation. Using Flac3D software, a circular roadway model is established that is 400 m underground, 12 m long, and has a radius of 3 m, as shown in Figure 9.
The constitutive model adopts the Mohr-Coulomb model, which requires six parameters: bulk modulus, shear modulus, internal friction angle, dilatancy angle, cohesion, and tensile strength. The bulk modulus and shear modulus are obtained from the elastic modulus and Poisson's ratio, as shown in Eq. 19. Frontiers in Earth Science frontiersin.org In the formula, K* is the bulk modulus, G is the shear modulus, and v is Poisson's ratio.
According to the Griffith strength criterion, the tensile strength of rock is generally one-tenth of the compressive strength. For soil or rock, the dilatancy angle ranges from 0°to 20°. In the Flac3D model, the default value is 0°. The outside of the circular roadway is homogeneous surrounding rock. The rock mass test parameters are as follows: volume weight = 25 KN/M 3 , cohesion = 0.8 Mpa, internal friction angle = 30°, Poisson's ratio = 0.3, and elastic modulus = 24.9 GPa.

Numerical solution of the uniaxial compressive strength
When the uniaxial compressive strength changes and other parameters remain unchanged, Figure 10 is obtained. Figure 10A represents the vertical displacement curve of a circular roadway when the uniaxial compressive strength is 36.7 MPa, and (B) represents the vertical displacement curve of a circular roadway when the uniaxial compressive strength is 231.7 MPa. It can be concluded that the vertical displacement increases with the increase of the load step. When the load step is the same, the vertical displacement value is small when the uniaxial compressive strength is large.

Numerical solution of the elastic modulus
When the elastic modulus changes and other parameters remain unchanged, Figure 11 is obtained. Figure 11A represents the vertical displacement curve of the circular roadway when the elastic modulus is 5.2 GPa, and (B) represents the vertical displacement curve of the circular roadway when the elastic modulus is 79.8 GPa. It can be concluded that the vertical displacement increases with the increase of the load step. When

FIGURE 11
Vertical displacement curve when the elastic modulus changes. (A) Vertical displacement curve when the elastic modulus is 5.2 GPa (with support). (B) Vertical displacement curve when the elastic modulus is 79.8 GPa (with support).

FIGURE 12
Vertical displacement curve when the internal friction angle changes. (A) Vertical displacement curve when the internal friction angle is 20°(with support). (B) Vertical displacement curve when the internal friction angle is 40°(with support).

Frontiers in Earth Science
frontiersin.org the load step is the same, the vertical displacement value is small when the elastic modulus is small.

Numerical solution of the internal friction angle
When the internal friction angle changes and other parameters remain unchanged, Figure 12 is obtained. Figure 12A represents the vertical displacement curve of the circular roadway when the internal friction angle is 20°, and (B) represents the vertical displacement curve of the circular roadway when the internal friction angle is 40°. It can be concluded that when the vertical displacement is the same, the load step is less when the internal friction angle is small, indicating that the displacement changes faster when the internal friction angle is small.

Conclusion
(1) A mechanical model of a circular "roadway rock-support" system considering damage is established. The critical load formula of the rock burst in the roadway with and without supporting conditions is derived according to the instability theory of rock burst disturbance response.
(2) Under the conditions of support and no support, the critical rock burst load of the roadway increases with the increase of uniaxial compressive strength. When the uniaxial compressive strength is the same, the critical load of the rock burst with support is 20 times larger than that without support. The critical load of rock burst increases with the increase of support stress, which is about 400 times larger than the support stress. The stability of the support and the roadway surrounding the rock system can be greatly improved by increasing the support strength.
(3) Under the condition of support, the critical load of rock burst increases with the internal friction angle. In particular, when the internal friction angle is greater than 35°, the critical load of rock burst increases rapidly with the increase of the internal friction angle. When the internal friction angle is the same, the critical load of the rock burst with support is much larger than that without support, and the difference is nearly three orders of magnitude. (4) Under the support condition of the roadway, the critical load of rock burst increases with the increase of the softening modulus. The critical load of rock burst experienced three stages of rapid decrease, slow decrease, and basic stability with the increased elastic modulus. When the values of the softening modulus and the elastic modulus are the same, the critical load with support is much larger than that without support.
(5) A new index of impact tendency K is defined as the ratio of the softening modulus to the elastic modulus. When the new index K of impact tendency is the same, the critical load of the rock burst with support is much larger than that without support. The conclusions have certain limitations due to the complex rock burst mechanism. The assumption of strain equivalence, isotropic, and homogeneity should be continuously reexamined in subsequent research into establishing a circular section roadway model.

Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors without undue reservation.