Experimental Study on the Mechanical Properties and Damage Evolution of Hollow Cylindrical Granite Specimens Subjected to Cyclic Coupled Static-Dynamic Loads

To study the characteristics of roadway surrounding rock damage caused by frequent disturbances under different static stress conditions, cyclic impact tests on granite with vertical holes under different axial prestress conditions were performed by a modified split Hopkinson pressure bar test, and the damage of the specimens was recorded with a high-speed camera process. The test results show that under the same air pressure cyclic impact, the rock specimens mainly undergo the compactionfatigue-failure transition. As the axial prestress increases, the compaction-fatigue phase gradually weakens, and the dynamic compressive strength decreases. When the axial prestress is 42% of the UCS and 62% of the UCS, the rock specimen shows a certain “strengthening” effect during the initial cyclic impact stage. During the failure of the rock specimens, the axial prestressing effect limited the initiation of some transverse cracks, and a mixed tensile-shear failure mode appeared. The rock specimens with an axial prestress of 62% of the UCS showed energy release during cyclic impact. To some extent, the probability of “rock bursts” has been induced. Based on the one-dimensional stress wave theory, the damage variables of wave impedance during the cyclic impact loading of the rock with vertical holes are defined. It is found that when the rock specimen is in the stage of compaction and fatigue damage, the damage is small, and the damage is even reduced.


Introduction
As the main carrier of loads in underground engineering applications, the structure and structural stability of rock are closely related to engineering safety. In recent years, with the continuous lack of shallow resources and the need for economic development for various countries, an increasing number of underground spaces have been developed at greater depths, such as mining, highway tunnel construction, and roadway excavation engineering projects [1][2][3], which trigger the continuous change of the stress state of the deep surrounding rock. In the process of underground rock excavation, the surrounding rock structure will inevitably be subjected to ground stress and mechanical drilling and blasting. The common effect of frequent disturbance [4,5] and various factors lead to deformation, failure, and instability of the rock mass structure, which leads to a series of engineering disasters, such as intensified instability of the roadway [6][7][8] and rock burst [9,10]. Therefore, one-dimensional dynamic static coupling can better reflect the actual stress state of roadway excavation. The study of the dynamic characteristics of the surrounding rock under dynamic and static coupling has important reference significance for the safe construction [11,12], efficient operation [13], and long-term stability of underground rock engineering.
Therefore, domestic and foreign scholars have performed much research on the mechanical properties of rocks under coupled static-dynamic loads. Zhou et al. [14][15][16][17][18] used an improved Hopkinson pressure bar system to analyse the mechanical properties of rock under the combined action of one-dimensional dynamic and static combinations and tested dried sandstone and water-saturated sandstone with different axial prestresses to study the influence of water saturation on the dynamic characteristics of rock under static prestress. Wang et al. [19] carried out experimental studies on different axial pressures and different impact loads on red sandstone and analysed the effects of axial pressure and impact loads on their mechanical properties and energy loss. Gong et al. [20][21][22] studied the dynamic characteristics of sandstone under coupled static-dynamic loads. The study found that when sandstone is loaded and destroyed within a certain range of axial pressure ratios, as the incident energy increases, the sandstone successively undergoes the three stages of "absorbed energy-release energy-absorbed energy". The above studies are based on the experimental study of intact rock samples under impact loading and have obtained abundant research results. However, rock failure and instability usually begin with primitive defects such as rock cracks and holes [23,24]. At present, there are few studies on defective rocks under impact loading. Li et al. [25][26][27] performed a one-dimensional dynamic and static combination test on samples containing circular and elliptical holes and found that the size, shape, and location of the holes have a certain effect on the dynamic strength of the rock, and the existence of holes induces shear failure of rocks. Wang et al. [28] used Hopkinson rods and combined digital image correlation technology to carry out impact experiments on rock-like samples with an angle of 0°to 90°between the prefabrication and loading directions and analysed samples under conditions of multiple loading rates, single cracks, and multiple angles. Li et al. [29] tested granite with square and circular holes to test the instability failure characteristics and obtained the relaxation time T2 spectrum curve, nuclear magnetic porosity, and nuclear magnetic resonance images under nuclear magnetic resonance (NMR) tests. In summary, most of the previous studies utilized impact tests based on prefabricated lateral cavity rock samples perpendicular to the loading direction, and few studies have been conducted on vertical holes. In the actual underground tunnel excavation process, frequent dynamic disturbances generated by mechanical drilling and blasting in the form of stress waves act on the surrounding rock of the excavated tunnel, which has a greater impact on its stability, and the surrounding rock of the tunnel when subjected to dynamic load disturbances. The rock is subjected to static loads such as in situ stress and tectonic stress [30], as shown in Figure 1. Therefore, understanding the failure mechanism of roadway surrounding rock under cyclic impact loading and certain static stress conditions is the essential basic information and basis for scientifically evaluating the stability of surrounding rock during the excavation of underground roadways.
Thick wall cylinders are the preferred model for simulating underground tunnels, tunnels, mine shafts, and other structures [31]. To correctly understand the mechanical properties, energy, and damage evolution of surrounding rock under the coupling action of static stress and dynamic disturbance in an actual project, this paper selects a granite cylinder with inner and outer diameters of 5 mm and 25 mm, respectively, and uses the modified Hopkinson bar system of Central South University to carry out cyclic impact tests on the granite under the same air pressure and different axial prestresses and uses high-speed photography. The image machine records the failure process of the specimen and analyses the mechanical behaviour, failure mode, energy evolution, and damage characteristics of the specimen under dynamic static coupling action.

Specimen Preparation and Test
System. The test material is granite with relatively good integrity and uniformity, processed into a 50 mm × 100 mm cylindrical rock sample. Then, a vertical hole with a diameter of 10 mm is prefabricated in the centre of the circular cross-section, and the two cross sections at both ends of the sample, which exhibit nonparallelism and nonverticality, are less than 0.02 mm, as shown in Figure 2. To reduce the discreteness of the rock samples, a rock acoustic wave parameter test system (HS-YS4A type) is used to test the P-wave velocity of the rock samples. The average density of the rock samples is 2:79 × 10 3 kg/m 3 , the average P-wave velocity is 5345 m/s, and the average compressive strength is 140 MPa.

Testing Equipment.
The test equipment uses the dynamic and static combined loading test system based on the SHPB device modified by Central South University. It consists of a punch device, an incident rod, a transmission rod, a buffer rod, and an axial compression device, as shown in Figure 3. The rods of the test system are made of 40Cr alloy steel with a diameter of 50 mm. The longitudinal wave velocity is 5400 m/s, the density is 7810 kg/m 3 , and the elastic modulus is 240 GPa. The data acquisition instrument and display equipment are a CS-1D super dynamic strain gauge and DL-750 oscilloscope, respectively. A spindle-shaped bullet is used in the launch cavity to eliminate wave oscillations and achieve stable half-sine wave loading to achieve a constant strain rate loading effect. To facilitate the observation of the microdestruction process of the test piece, a FAS-TCAM SA1.1 high-speed camera is used to synchronously record the crack propagation process of the test piece. The frame number of the high-speed camera is set to 43,200 fps

Test Procedure.
To test the effectiveness, before the start of the test, the incident rod and the transmission rod are closely pressed for empty punching, and the trend of the incident wave and the transmitted wave on the oscilloscope is observed to ensure whether the two rods are normally secured.
To achieve the purpose of cyclic impact of the rock sample under different prestressing effects and to avoid the rock sample being damaged, it is necessary to test the rock sample and finally set the impact pressure to 0.8 MPa.
During the test, lubricating oil should be applied between the incident rod and the transmission rod to reduce the end effect of the sample, the rock should be sandwiched between the incident rod and the transmission rod, and then the manual oil pump should be slowly prestressed to the set value. The strain gauge, oscilloscope, and high-speed camera are initiated to record data.

Test
Principle. The SHPB system can simultaneously realize dynamic and static combined tests with an axial prestress of 0-200 MPa and an impact dynamic load of 0-500 MPa. The punch impacts the incident rod at a certain rate under the action of air pressure, and an incident stress wave is formed in the incident rod. When the incident stress wave is transmitted to the end face of the rock sample, transmission stress is formed on the end face of the rock sample due to the different wave impedances of the rock sample and the rod. The incident and transmitted waves can be measured from incident rod strain gauge 1 and transmitted rod strain gauge 2, respectively, as shown in Figure 2. Using the "three wave method" [32] to process the collected wave   3 Geofluids forms, we can obtain the changes in the strain, stress, and strain rate with time as follows: In the SHPB test, it is assumed that the energy loss generated by the end face of the specimen is negligible. The energy before the test is equal to the energy after the test. According to the law of the conservation of energy [32], the incident energy, reflection energy, transmitted energy, and formula for absorbed energy are as follows: where A e and A s denote the cross-sectional area of the elastic rod and the sample, respectively; E e represents the elastic modulus of the elastic rod; C e and L s represent the longitudinal wave velocity of the elastic rod and the length of the sample; ε I ðtÞ, ε R ðtÞ , and ε T ðtÞ represent incident, reflected, and transmitted strain wave signals, respectively; εðtÞ and σðtÞ are the strain and stress, respectively; W I , W R , W T , and W S are the incident energy, reflected energy, transmitted energy, and absorbed energy, respectively.

Test Results and Discussion
The posttest data are processed based on the three-wave method. Due to space limitations, some test specimens are selected for analysis in this article. The test results are shown in Table 2.

Mechanical Behaviour of Granite under Dynamic and Static Coupling
3.1.1. Dynamic Stress-Strain Curve. Figure 5 shows the stressstrain curve of a rock sample during cyclic impact under different axial prestressing forces. Figure 5 shows that the stressstrain curve can be roughly divided into a compaction stage, a crack propagation stage, and a postpeak unloading stage. In the crack growth stage, when a relatively smooth curve appears before the peak stress is approached, the stress growth is slow, and the strain growth rate is faster as the yield stage continues; when the strain increases to a certain level, the peak stress rapidly increases and is defined as the strengthening stage. As the axial prestress increases, the yield stage becomes less and less obvious. Before reaching the peak stress, the change trends of the stress-strain curves are similar, but there are obvious differences in the unloading phases after the peak. Before the failure of the rock sample, the maximum strain at the unloading stage after the peak increases first and then decreases with the decrease in the stress, showing a rebound phenomenon; when the rock sample is impacted and destroyed, the rebound phenomenon gradually disappears.

Effect of the Axial Prestress on the Peak Stress.
In the process of cyclic impact loading under different axial prestresses, the peak stress and the number of cyclic impact changes are shown in Figure 6. As the number of cycles increases, the peak stresses of Z-0-1 and Z-C-1 continue to decrease, while the peak stresses of rock samples Z-A-3 and Z-B-2 have an increasing trend during the initial impact stage, and their peak stress performance thereafter gradually decreases. According to the analysis, in the initial impact stage of samples Z-A-3 and Z-B-2, internal microcracks are closed under the impact load, which means that the wave  impedance of the rock sample increases and the peak stress increases [33]. In the subsequent impact process, the internal cracks of the rock samples gradually increase, the wave impedance of the samples gradually decreases, and the ability to resist external loads continues to weaken, which reflects the strength degradation characteristics of the rock sample during cyclic impact. However, due to the large axial prestress, Z-C-1 causes many internal microcracks, which greatly weakens the resistance of the rock sample to external loads. Under dynamic loading, the peak stress shows a decrease with an increase in the number of cycles. Figure 6 can further illustrate the effect of the axial prestress on the peak stress of the rock samples under the same air pressure. It can be seen from Figure 6 that the dynamic peak strength of unstressed rock samples is higher than that of prestressed axial forces, and the dynamic peak strength decreases as the axial prestress increases. When the axial prestress is 42% of the UCS and 62% of the UCS, the number of impacts increases compared to the case without axial prestress, and when the axial prestress is 83% of the UCS, the number of impacts is reduced. From the above phenomenon, it can be determined that the axial prestress has a certain effect on the dynamic strength and stiffness of the material. These phenomena can be explained by different stress stages under uniaxial compression [14]. Generally, the behaviour of brittle rock can be divided into three stages during uniaxial compression: (I) elastic deformation stage: approximately 0-30% of the UCS; (II) crack steady propagation stage: approximately 30-70% of the UCS; and (III) unstable crack growth stage: approximately 70-100% of the UCS, as shown in Figure 7.
The axial prestress is in the stage of stable crack propagation, and activated cracks regenerated inside the rock sample.
When the peak strength of the dynamic load is lower than a certain value, during the initial cyclic impact of the rock sample, the damage caused by the dynamic load is generally lower than the crack closure amount, and the rock sample shows a compaction process [34]. When the axial prestress is 42% of the UCS and 62% of the UCS, the dynamic strength of the rock sample is lower than that without prestressing, indicating that the prestress has caused some damage to the rock sample, but due to the axial prestressing effect, the local damage has an impact on the specimen failure and has a strengthening effect on the rock material. Therefore, the stiffness is enhanced compared with the prestress. At an axial prestress of 83% of the UCS, the rock sample is in the stage of unstable crack propagation. The internal cracks greatly deteriorate the mechanical properties of the rock, showing its minimum dynamic strength and minimum number of impacts.

Variation in the Average Strain Rate and Maximum
Strain. The average strain rate represents the amount of strain per unit time of the rock during a single impact [34]. Figures 8 and 9 show the relationship between the average strain rate and the maximum strain with the number of cycles when the rock samples are subjected to different axial prestresses. It can be seen that the average strain rate and the maximum strain increase as a whole. Z-0-1 shows a relatively large average strain rate and maximum strain during the first impact. According to the analysis, some initial cracks cannot be avoided inside the rock sample, and more stress may be generated around the hole wall during the impact due to the stress concentration. The crack shows a larger strain rate and maximum strain during the initial impact; the average strain rate and maximum strain of Z-A-3 during 5 Geofluids the first impact are lower than those of Z-0-1, which is due to the axial prestress. Effected by the initial crack closure, the first time the impact is initiated with the same air pressure, the crack closure amount is relatively small. In the same way, with the increase in the axial prestress, the internal crack closure of the rock samples increases, and the average strain rate and maximum strain of Z-B-2 and Z-C-1 are substantially reduced compared with Z-A-3 at the first impact. As the number of impacts accumulates, Z-A-3 and Z-B-2 show an average strain rate and a maximum strain that first decrease and then increase, and the strain of Z-C-1 increases as a whole. The analysis shows that during the initial impact of Z-A-3, the crack closure amount of the rock sample is larger than the crack initiation amount, showing a decreasing trend of the average strain rate and that the maximum strain decreases; with the increase in the number of cycles, the rock sample is in the third stage. During the fifth impact process, the average strain rate and the maximum strain change range decrease, indicating that the microcracks at this stage of the rock specimens develop slowly. The analysis indicates that  6 Geofluids the rock specimens transition from the initial compaction stage to the fatigue damage stage and finally collapse. At this time, cracks in the rock samples rapidly expand, and the average strain rate and maximum strain increase substantially. During the initial impact of Z-B-2, the average strain rate and the maximum strain change are small, and the crack growth rate of the rock sample is very slow or does not expand. The average strain rate and the maximum value during the third impact should reach a minimum, indicating that the rock sample is in the first three stages. The crack closure stage is at the initial stage of fatigue damage; when the rock sample is impacted for the fifth time, the average strain rate and maximum strain of the rock sample increase substantially, indicating that the crack growth speed of the rock sample accelerates, which represents the transition from the fatigue damage stage to the failure stage. The Z-C-1 strain rate and average strain rate increase with increasing number of cycles. According to the analysis, a large number of microcracks appear in the sample due to the prestressing effect. When the impact load is applied, the initiation cracks are promoted to develop into the macrocrack direction. No compaction stage or fatigue damage stage is shown. In summary, when the axial prestress is low, the rock sample undergoes three stages during the cyclic impact process: the compaction phase, fatigue damage phase, and failure phase; as the prestress increases, the rock sample compaction phase weakens. There are two phases: the fatigue damage phase and failure phase; when the prestress is greater than a certain value, the rock sample only shows the failure phase.

Effect of Different Prestressing Effects on the Dynamic
Deformation Modulus. Figure 10 shows the relationship between the dynamic deformation modulus of the rock sample and the number of cycles during different axial prestressed cyclic impacts. The dynamic deformation modulus reflects the resistance of the rock to deformation. Because the dynamic stress-strain curve has no obvious straight line segment, to better reflect the rock's resistance to deformation during impact, reduce errors, and reduce test dispersion, this   7 Geofluids paper adopts the following calculation method proposed by Tang et al. [35]: where E 1 is the first secant modulus, E 2 is the second type of secant modulus, E 3 is the deformation modulus of the loading section, E d is the dynamic deformation modulus, σ d50 is 50% of the peak stress, andε d50 is 50% of the peak stress. The corresponding strain σ d is the peak stress, ε d is the peak strain, α is the angle between the tangent and the axis at 50% of the peak stress, and E 0 is the elastic modulus of the material without damage.
In combination with the stress-strain curve in Figure 3, it can be seen that when the axial prestress is low and when the dynamic load initially acts on the rock sample, the internal cracks are further closed, showing that the dynamic deformation modulus increases first and then decreases as the number of cycles increases. When the axial prestress is greater than a certain degree, there are more microcracks in the rock sample, but the overall structure is in a state of stress equilibrium due to axial compression. When the dynamic load is applied, the internal cracks expand rapidly, showing that the dynamic deformation modulus changes. The number of cycles increases at a smaller level. The overall trend of the Z-0-1, Z-A-3, and Z-B-2 dynamic deformation moduli first increases and then decreases, while Z-C-1 shows a decreasing trend. The increase in the dynamic deformation modulus means that its resistance to deformation is enhanced. The dynamic deformation modulus of the Z-A-3 and Z-B-2 rock samples increases during the first three impacts, indicating that the rock samples are in the crack-closing stage under the impact load, showing a certain hardening effect. The dynamic deformation modulus followed by the cycle number decreases, and the resistance of the rock samples to deformation weakens.
It is further seen from Figure 8 that the maximum deformation modulus of the axial prestress of 42% of the UCS and 62% of the UCS is substantially higher than when the axial prestress is not applied. This is because the rock is enhanced when the axial prestress is low. The ability of the sample to resist deformation and the dynamic deformation modulus of the rock sample under impact failure are lower than when axial prestress is not applied. According to the analysis, the rock sample has sufficient internal crack expansion when subjected to multiple impact loads. The deformation modulus is low. When the axial prestress is 83% of the UCS, the maximum dynamic deformation modulus of the rock sample and the dynamic deformation modulus at failure are lower than those without axial prestress. This is because when the axial prestress is large, the prestress promotes the formation of microcracks in the rock sample, reducing its resistance to deformation.

Energy Evolution Characteristics of Granite under
Dynamic and Static Coupling. Figure 11 shows the relationship between the absorbed energy of the rock sample and the number of cycles under different axial prestressing forces. It can be seen from Figure 11 that the variation law of the absorbed energy varies with different axial prestresses. When axial prestress is not applied, the absorption energy increases with the number of cycles, and the absorption energy is positive, which means that during the cyclic impact process, the rock sample continuously absorbs energy for the initiation and expansion of its internal cracks. In addition, when the axial prestress is 42% of the UCS, the rock sample increases with the number of cycles, and the energy is converted from the released energy to the absorbed energy. This is because the rock sample has a certain amount of energy stored during the axial prestress, and the first two impact loads, and during the action, the rock samples exhibit energy release, and the released energy decreases as the number of impacts increases. The analysis shows that during the initial impact process, the internal structure of the rock sample is fine-tuned, and the microcracks are continuously closed, resulting in a decrease in the released energy and energy as the number of cycles increases. During the third impact, the rock sample begins to absorb energy, and the absorbed energy is 1.14 J, indicating that the rock sample has a relatively stable internal structure during the third impact, and there is basically no energy dissipation; the rock sample is subjected to the fifth impact. When the load is applied, the energy absorption value is 7.88 J, which is substantially larger than that during the third and fourth impacts, indicating that the crack growth rate and degree of the rock sample increase. When the rock sample is impacted and damaged, the energy absorption value is 28.4 J, which increases sharply, indicating that the specimen has more absorbed energy when it is destroyed, and the degree of damage increases. The axial prestress is 62% of the UCS, the rock sample releases energy during cyclic impact, and  8 Geofluids the overall trend is decreasing. The analysis shows that the rock sample stores substantial strain energy. As the number of cycles accumulates, the internally activated cracks will increase. Due to the release of energy, the internal cracks will continue to expand. During impact failure, the residual elastic energy inside the rock sample will be greater than the fracture. The surface requires energy, and under the action of dynamic loading, it triggers the release of its energy, which eventually leads to a "rock burst" in the rock sample. The axial prestress is 83% of the UCS, and more cracks have been induced inside the rock sample, resulting in less stored energy. During the first impact, the energy is released internally, causing more cracks to activate and propagate inside, and a certain macroscopic damage occurs on the surface. In this state, because the internal energy storage is low, the rock sample is under the impact of the impact load. "Rock bursts" are less likely to occur in this stage, so the specimens show absorbed energy when they finally fail.

Dynamic and Static Coupling Failure Modes of Granite.
Understanding the rock failure mode is of practical engineering significance for underground mining processes and the prevention of hazards such as rock bursts. Figure 12 shows the fracture process of specimens under different axial prestressing forces under the same air pressure impact, and Figure 11 shows the corresponding failure mode. When the axial prestress is not applied, under the impact of the dynamic load cycle of the sample, the reflected wave and the transmitted wave in the rock sample cannot be transmitted to the two ends of the sample in time, resulting in a superposition at a certain position inside the rock sample and a greater tensile wave. The existence of holes further weakens the tensile strength of the rock sample so that microcracks in other positions of the rock sample will not be able to propagate and penetrate in the future, and a through tensile crack perpendicular to the loading direction appears at a certain position, triggering the damage of the rock sample. The spalling damage leads to a larger size of the final damage and less damage. Figure 12(a) shows a diagram of the failure process of the rock sample without axial prestress failure. At 1104 μs, cracks perpendicular to the loading direction begin to appear in the centre of the rock sample, and debris appears at the incident end of the rock sample, which eventually leads to spalling failure of the rock sample at 3335 μs.
When the rock sample is under axial prestress, a potential compression-shear surface is formed inside. During the initial cyclic impact process, because the impact strength is lower than the maximum compressive strength of the prestress, the rock sample does not undergo macroscopic damage. As the cyclic impact accumulates to a certain degree, the maximum strength of the prestress is continuously weakened. When the impact strength of the prestress is higher than the minimum strength to which the prestress is subjected, the rock sample will fail.
When the axial prestressing force is 42% of the UCS and 62% of the UCS, the axial prestress and dynamic load act simultaneously, resulting in mixed tensile-shear failure. When the axial prestress is 42% of the UCS, Figure 12(b) shows that when the transmitting end of the rock sample is at 276 μs, compression cracks begin to appear, and at 345 μs, tensile cracks appear at the ends while the end face of the rock sample is sheared and fragments began to peel off. At 598 μs, a tensile crack appears on the surface of the rock sample. When the axial prestress is 62% of the UCS, Figure 12(c) shows that the rock sample begins to undergo compression-shear failure when the rock sample is broken. At 161 μs, compression shear cracks appear from the entrance side of the specimen. At 368 μs, the fragments on the compression shear surface begin to eject from the side and make a loud noise. This phenomenon is similar to a "rock burst" in rock engineering.
When the axial prestressing force is 83% of the UCS, the axial prestress plays a major role, and shear tensile mixed failure occurs. It can be observed from Figure 12(d) that the rock sample cracks at 253 μs, and at 414 μs, chipping occurs on the side of the rock. The analysis shows that the prestress is 83% of the UCS, a more obvious compression zone is formed at the end of the sample, and the internal cracks parallel to the loading direction expand to the inside [21], forming a potential "shear failure surface", resulting in more obvious compression on both sides. The shear rupture surface, when impact load is applied, further promotes this formation trend and triggers rock chipping on the side of the rock.
Sliding loads act on different axial prestress failure modes. Figure 13 shows that as the axial prestress increases, the failure block size of the rock samples does not decrease substantially. The main reason for this behaviour is that in this study, granite was not severely damaged. When the prestress is in the stage of stable crack growth, which is lower than the crack damage stress (70% of the UCS) [14], under cyclic impact loading, the rock sample will undergo certain fatigue damage, and its internal microcracks gradually change from the initial compaction stage. The crack propagation stage is transitional, and the crack propagation speed is slow. When the final impact is broken, the crack is full, and the fragmentation is small. When the prestress is higher than the crack damage stress, the number of cyclic impacts of the rock sample is lower, which indicates that the main crack is propagating faster, and the internal microcracks have time 9 Geofluids to grow. The main crack expands and penetrates, resulting in larger fragments during failure.

Selection and Establishment of the Dynamic Rock
Damage Model. At present, the commonly used damage variables are defined as the elastic modulus, maximum strain, and ultrasonic wave speed. The ultrasonic wave speed is widely used because it is easy to operate and can effectively reflect the crack propagation trend in rock, but it has a cycle under a certain axial prestress. Impact tests limit the application of this method. For this reason, Jin Jiefang established a method based on the one-dimensional stress wave theory. Under the same cross-sectional area of the test specimen and the elastic rod, the wave impedance is used to define the damage variable for the first time, as shown in Equation (4). The elastic modulus, energy dissipation, and maximum strain are not suitable for defining the damage variables of rock during cyclic impact [36]. In this paper, the vertical centre hole is prefabricated in the centre of the sample. Therefore, Equation (4) cannot be directly applied to this test. Assuming that the discontinuous contact surface between the rock sample and the elastic rod is regarded as the contact surface between two elastic halfspaces with displacement discontinuity, the propagation of the stress wave on the displacement discontinuity surface 10 Geofluids can be transformed into a boundary value problem for solving the wave equation [32], as shown in Figure 14.
where D is the damage variable, (ρC) ̅ is the wave impedance at a certain impact, and ρC is the initial wave impedance of the rock. As shown in the figure, based on the one-dimensional wave theory, when the stress wave propagates from incident rod a to test piece b, its wave impedance changes from m a to m b . According to the one-dimensional stress wave theory, the interface has a continuous force and velocity on O 1 O 2 with the following available conditions: In the same way, σ T ′ ðt i Þ inside the rock enters the transmission rod at interface O 3 O 4 to generate a transmitted wave: The relationship between the transmitted wave and the incident wave at any time can be written as follows: If m a and m b are known, Equation (11) can be rearranged, and Equations (7), (8), and (9) can be simultaneously used to obtain the actual wave impedance of the rock specimen at a certain moment during the i time impact: By taking Equation (12) into Equation (4), the expression of the damage degree of the rock is as follows: where σ I , σ R , σ T , and σ ′ T are the incident wave, reflected wave, transmitted wave, and transmitted wave in the specimen, respectively; A a and A b are the cross-sectional areas of the elastic rod and the rock sample, respectively; ρ a and C a , and ρ b and C b are the density and P-wave velocity of the elas-tic rod and rock sample, respectively; λ is the reflection coefficient of the wave entering the specimen from the elastic rod; and i is the ith time impact order. The above inference formula assumes that the stress wave in the rock specimen generates a transreflective reflection. After the incident wave propagates into the specimen for a period of time, multiple transmissive reflections occur in the specimen at each instant. From the research in [36], it is evident that when t ≤ ð2L b /C b Þ, the reflected and transmitted waves obtained in the test are not affected by multiple transmissions and reflections. Figure 15 shows a trend chart of the change in the wave impedance when a certain specimen is impacted. The impedance of section AB basically remains unchanged. To produce comparable results, a certain time in section AB is selected as the reference point for each cycle of the test specimen. Figure 16 shows the change in the damage degree and impact number of the rock samples under different axial prestressing forces. It can be seen from the figure that the rock samples show different damage evolution trends with the increase in the number of cycles under different prestresses. The initial damage of Z-0-1 is lower than that of the other rock samples, and the initial damage of Z-A-3, Z-B-2, and Z-C-1 increases with the increase in the axial prestress. The rock samples have different degrees of damage when the prestress is 42% of the UCS, 62% of the UCS, and 83% of the UCS. Z-C-1 has a large number of microcracks because the axial prestress is higher than the crack damage stress, resulting in a substantial reduction in the number of cycles of rock samples subjected to cyclic impact. The cumulative growth rate of the damage of Z-0-1 continues to increase, with growth rates of 60% and 150%, respectively. During the initial cyclic impact process, Z-A-3 is dominated by compaction, resulting in the reduction in the first three damage levels, and then with the increase in the number of cycles, the cumulative damage growth rate of the rock sample increases with the acceleration of the internal crack growth. Z-B-2 increases the cumulative damage rate during the first four cycles of slow impact, with a range of 0.24-0.25. From the establishment of the dynamic damage model, it can be seen that the damage D is a function of the strain [37], which is consistent with the maximum strain in the previous analysis. Since the rock sample belongs to the fatigue damage stage at this stage, the amplitude of the strain change is small, which means that during this impact process, the damage to the rock sample is small and can almost be ignored. This shows that the cumulative damage growth of the rock sample during the first four impacts is slow, and the damage accumulation during the subsequent cyclic impact process increases with the increase in the maximum strain.

Surrounding Rock Engineering Applications of Roadways under Different Prestresses
The in situ stress and dynamic load affect the safety and stability of the roadway. Among them, the magnitude of the ground stress (prestress) determines the strength and stiffness of the surrounding rock of the roadway. Figure 17 shows the maximum peak stress and dynamic deformation modulus, which indirectly reflect the rock mass stiffness. The test results show that when the prestress is higher than the crack initiation stress (30% of the UCS), the dynamic strength and stiffness are reduced and increased, respectively, compared with the absence of axial prestress, which indicates that the in situ stress improves the "toughness" of the surrounding rock. Under the same dynamic load disturbance action, the ability to resist the instability and damage of the surrounding rock of the roadway is enhanced. When the axial prestress is higher than the crack damage stress (70% of the UCS), the strength and stiffness are reduced compared with the axial prestress. It can be concluded that under high prestress conditions, the dynamic load disturbance is extremely detrimental to the stability of the surrounding rock of the roadway, and it is more likely to cause instability. The failure type and energy evolution law of the surrounding rock of the roadway are also important features that reveal the failure mechanism of the surrounding rock. For tunnels with shallow depths, it can be considered that there is no static stress. If the surrounding rock of the roadway is to be destabilized and destroyed, it is necessary to absorb the external energy, and the failure of the surrounding rock of the roadway is usually caused by the tensile spalling damage caused by the superposition of the stress waves due to the dynamic load at a certain position of the roadway surrounding rock. For the surrounding rock of the roadway where the prestress is in the stage of stable expansion of the crack, as the depth increases, a certain degree of shear surface is formed inside, and the stored energy increases; when the external load is disturbed, the surrounding rock of the roadway is transformed from absorbed energy to released energy when the surrounding rock is destroyed, which results in a greater probability of triggering a rock burst. For the prestressed roadway surrounding rock in the unstable stage of crack expansion, which is higher than the crack damage stress, the structure of the roadway surrounding rock is almost in an unstable state, there are many microcracks inside, the rock density is very low, and the stored energy is relatively small. When the load is disturbed, the internal absorbed energy is greater than the released energy. In this case, rock burst does not easily occur, and the roadway is easily induced to collapse and fall off.

Conclusion
(1) The axial prestress has a certain effect on the dynamic strength and stiffness of granite: the axial prestress is in the stage of stable crack propagation, and its dynamic strength is lower than that of the pure dynamic compressive strength, but its stiffness increases. The axial prestress is in the unstable stage of crack propagation, showing its weakest dynamic strength and stiffness (2) During the cyclic impact failure process of granite, the granite mainly undergoes a transition from the compaction stage to the fatigue damage stage to the failure stage. With the increase in the prestressed axial pressure, the compaction stage and the fatigue   damage stage continue to weaken. When the axial prestress is 42% of the UCS, the rock sample shows the compaction stage-fatigue damage stage-failure stage trend; when the axial prestress is 62% of the UCS, the rock sample compaction stage weakens, showing the fatigue damage stage-failure stage trend; and when the axial prestress is 83% of the UCS, the rock sample only shows the failure stage (3) The existence of static stress is a necessary condition for "rock burst" to occur. Granite under different axial prestresses, the energy dissipation of rocks varies with the number of impacts. In general, when the axial prestress and impact load are at a certain value, the probability of "rock burst" is large (4) The failure mode of granite under different axial prestresses is different. When the axial prestress is 0, spalling tensile failure occurs. As the axial prestress increases, mixed tensile-shear failure occurs. This is due to the effect of axial prestressing, and different levels of shear planes are formed inside, which limits the initiation of transverse cracks (5) Based on the one-dimensional stress wave theory, the method of defining the rock damage degree by wave impedance is further studied. When the rock sample is in the compaction stage and the fatigue damage stage, the damage of the specimen is small, and the damage reduction phenomenon occurs

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request. Dynamic deformation modulus (GPa) Figure 17: Relationship between the dynamic deformation modulus, peak stress, and axial prestress.