Experimental Investigation on Uniaxial Compression Mechanical Behavior and Damage Evolution of Pre-Damaged Granite after Cyclic Loading

: Failure behavior of pillars in deep mines is affected by various cyclic loads that cause initial pre-damage. Pillars will be further damaged and developed in the long-term compressive stress until they are destroyed. To reveal the strength characteristics and crack damage fracture laws after rock pre-damage, uniaxial compression tests were carried out on granite specimens damaged by cyclic loading using the digital speckle correlation method. The experimental results indicate that the mechanical properties of pre-damaged specimens show large damage differences for different cycles. The damage variable of the pre-damaged specimens increases with the increase of cycle number and conﬁning pressure. The damage of specimens is primarily due to the strength weakening effect caused by cycle numbers, and the conﬁning pressure restriction effect is not obvious. The evolution laws of uniaxial compression damage propagation in the pre-damaged specimens show differences and obvious localization phenomenon. Pre-damaged specimens experienced three failure modes in the uniaxial compression test, namely tensile shear failure (Mode I), quasi-coplanar shear failure (Mode II), and stepped path failure (Mode III), and under different pre-damage stress environments with high conﬁning pressures, the failure modes are dominated by Mode II and Mode III, respectively. stress-resistance tests conducted on


Introduction
Instability in rock engineering is primarily due to the further expansion of joints, faults, and cracks [1]. There are studies that show the effect of mechanical excavation as a cyclic load in deep mining, especially with roadheader excavation [2], which will lead to stress redistribution on pillars or surrounding rocks. At the same time, the rock masses are subjected to cyclic loading, causing meso-damage of the surrounding rock or pillars [3]. Mechanical properties for permanently supported pillars are further damaged and deteriorated under compressive stress after excavation, which introduces safety hazards to the engineering structure [4].
In recent years, rock mass damage evolution research has made significant achievements. Fan et al. [5] conducted uniaxial compression tests on specimens containing multinon-persistent joints with PFC 3D , and they found that the joint specimens exhibited four failure modes. Hu et al. [6] proposed an interval value for the damage initial stress range of intact rocks under uniaxial compression. Xiao et al. [7] presented the fatigue damage evolution laws of rocks under cyclic loading and believed that it primarily depended on the maximum stress, amplitude, and initial fatigue damage. Bagde et al. [8] studied the effects of amplitude and frequency on the mechanical properties of rock under uniaxial cyclic 1. Four groups of rock specimens were vacuumed and saturated with water using the vacuum water saturation equipment, first for 240 min dry pumping and then for 120 min wet pumping to completely saturate the specimens; 2. Uniaxial compressive strength of specimens in Group D were measured with the MTS815 testing machine loaded at a rate of 2 kN/s (Group D was the control group); 3. The remaining rock specimens of Group A, B, and C were cyclically loaded and unloaded for 10, 20, and 30 times under confining pressures of 50 MPa, 80 MPa, 100 MPa, and 120 MPa, respectively, without destroying the specimens, to realize the pre-damage treatment inside the granite specimens (see Table 2 for the specific scheme). The load cycle path used constant confining pressure and unloading axial pressure; both loading and unloading rates were 2 kN/s, loading stress point was 153 MPa (about 93% of the average uniaxial compressive strength), and unloading stress point was 5 MPa (about 3% of the average uniaxial compressive strength).

Pre-Damage Treatment of Specimens
The granite specimen was a standard cylinder of 50 mm × 100 mm, with uniform dense texture and no damage. To obtain the pre-damaged granite specimens, the intact specimens were first subjected to a triaxial compression test. The cyclic loading scheme is shown in Figure 2. The granites were classified into four groups (A, B, C, and D), and the basic physical parameters of group A3, B3, C3, and D are shown in Table 1. The procedures are mainly as follows: 1.
Four groups of rock specimens were vacuumed and saturated with water using the vacuum water saturation equipment, first for 240 min dry pumping and then for 120 min wet pumping to completely saturate the specimens; 2.
Uniaxial compressive strength of specimens in Group D were measured with the MTS815 testing machine loaded at a rate of 2 kN/s (Group D was the control group); 3.
The remaining rock specimens of Group A, B, and C were cyclically loaded and unloaded for 10, 20, and 30 times under confining pressures of 50 MPa, 80 MPa, 100 MPa, and 120 MPa, respectively, without destroying the specimens, to realize the pre-damage treatment inside the granite specimens (see Table 2 for the specific scheme). The load cycle path used constant confining pressure and unloading axial pressure; both loading and unloading rates were 2 kN/s, loading stress point was 153 MPa (about 93% of the average uniaxial compressive strength), and unloading stress point was 5 MPa (about 3% of the average uniaxial compressive strength).

Uniaxial Compression Test with DSCM
The fundamentals of DSCM is to match the same pixel areas (sub-images) betw two digital speckle images of the specimen's surface before deformation (reference im

Uniaxial Compression Test with DSCM
The fundamentals of DSCM is to match the same pixel areas (sub-images) between two digital speckle images of the specimen's surface before deformation (reference image) and after deformation (deformed image) [22]. Its main purpose is to search the displacement component of the pixels in the speckle sub-image ( Figure 3). To identify the sub-images in the two images, the properties of the sub-images are represented by the gray function. Thus, by tracking the grayscale value of each pixel in the reference image and the deformation image and evaluating their image correlation, the displacement field on the specimen's surface can be calculated, and then the strain field can be obtained by further calculation. The correlation coefficient is defined as follows [23]: where f (x i , y i ) is the grayscale value of the reference image at the point (x i , y i ) and g x i , y i is the grayscale value of the deformed image at the point x i , y i ; f , g are the average grayscale values of the sub-images in the reference image and the deformed image, respectively. and after deformation (deformed image) [22]. Its main purpose is to search the displacement component of the pixels in the speckle sub-image ( Figure 3). To identify the subimages in the two images, the properties of the sub-images are represented by the gray function. Thus, by tracking the grayscale value of each pixel in the reference image and the deformation image and evaluating their image correlation, the displacement field on the specimen's surface can be calculated, and then the strain field can be obtained by further calculation. The correlation coefficient is defined as follows [23]: where , is the grayscale value of the reference image at the point , and , is the grayscale value of the deformed image at the point , ; ̅ , ̅ are the average grayscale values of the sub-images in the reference image and the deformed image, respectively. Combined with the DSCM method, secondary uniaxial compression tests were performed on the pre-damaged granite to explore the damage evolution and failure laws. The observation surface of each specimen was uniformly sprayed with a layer of matte white paint, and they were then randomly sprayed with matte black paint. After air drying, a uniaxial compression rigidity testing machine was used for displacement loading at 0.02 mm/min, and a CCD camera was used for shooting until the specimen was destroyed. The experimental scheme is shown in Figure 4. Combined with the DSCM method, secondary uniaxial compression tests were performed on the pre-damaged granite to explore the damage evolution and failure laws. The observation surface of each specimen was uniformly sprayed with a layer of matte white paint, and they were then randomly sprayed with matte black paint. After air drying, a uniaxial compression rigidity testing machine was used for displacement loading at 0.02 mm/min, and a CCD camera was used for shooting until the specimen was destroyed. The experimental scheme is shown in Figure 4.

Uniaxial Compressive Strength Results for Intact Specimens
The stress-strain curves of specimens (control group) in group D are illustrated in

Uniaxial Compressive Strength Results for Intact Specimens
The stress-strain curves of specimens (control group) in group D are illustrated in Figure 5, and the uniaxial compressive strength is 164.20, 162.57, and 167.87 MPa, respectively. To study the damage laws of the other groups of pre-damaged granite specimens during uniaxial compression, the uniaxial compressive strength of intact specimens is taken as the average value of 164.88 MPa. In addition, it can be seen from the stress-strain curve of uniaxial compression that intact granite belongs to brittle rock, and there is almost no plastic yield stage.

Uniaxial Compressive Strength Results for Intact Specimens
The stress-strain curves of specimens (control group) in group D are illustrated in Figure 5, and the uniaxial compressive strength is 164.20, 162.57, and 167.87 MPa, respectively. To study the damage laws of the other groups of pre-damaged granite specimens during uniaxial compression, the uniaxial compressive strength of intact specimens is taken as the average value of 164.88 MPa. In addition, it can be seen from the stress-strain curve of uniaxial compression that intact granite belongs to brittle rock, and there is almost no plastic yield stage.

Uniaxial Compressive Strength Results for Pre-Damaged Specimens
The stress-strain curves for pre-damaged specimens under 50 MPa confining pressure are illustrated in Figure 6. There are significant differences in the axial stress-strain curves after different cycles. The uniaxial compressive strength of pre-damaged specimens with 10 cycles is the largest, and the stress-strain curves showed obvious brittleness,

Uniaxial Compressive Strength Results for Pre-Damaged Specimens
The stress-strain curves for pre-damaged specimens under 50 MPa confining pressure are illustrated in Figure 6. There are significant differences in the axial stress-strain curves after different cycles. The uniaxial compressive strength of pre-damaged specimens with 10 cycles is the largest, and the stress-strain curves showed obvious brittleness, which is very similar to that of the intact granite specimens. With the increase of cycle numbers, the uniaxial compressive strength tended to decrease, and there was strain softening behavior after the peak, which exhibited significant ductility. This indicates that the mechanical properties of pre-damaged specimens in high cycle numbers (20 cycles and 30 cycles) were significantly different from those in low cycle numbers (10 cycles); especially under high cycle numbers, the mechanical properties of pre-damaged specimens decreased significantly. In addition, the stress-strain curve of the pre-damaged specimens has many micro-stress drops, which show pre-peak and post-peak fluctuations, indicating that the pre-damaged specimens were not destroyed once, but destroyed step by step. Because the pre-damage specimen was cyclically loaded and unloaded during triaxial compression pre-damage, it resulted in the repeated initiation and closure of cracks or micro-holes inside the specimen. When cycle numbers increased, the interaction of defects inside the specimen and its damage accumulation effect lead to the expansion of the local deformation area of the specimen under uniaxial compression, greatly reducing the mechanical properties of the specimen. Moreover, when the pre-damaged specimen had good bearing capacity, it is possible that yield failure could have first occurred locally in the specimen. Therefore, a slight stress drop could occur before the specimen fails completely. This result was observed in many studies [24][25][26].
cracks or micro-holes inside the specimen. When cycle numbers increased, the interaction of defects inside the specimen and its damage accumulation effect lead to the expansion of the local deformation area of the specimen under uniaxial compression, greatly reducing the mechanical properties of the specimen. Moreover, when the pre-damaged specimen had good bearing capacity, it is possible that yield failure could have first occurred locally in the specimen. Therefore, a slight stress drop could occur before the specimen fails completely. This result was observed in many studies [24][25][26].  Table 3 lists the average uniaxial compressive strength of pre-damaged specimens under different damage conditions, and its relationship with the cycle numbers and confining pressures is shown in Figure 7.

Strength Characteristics of Pre-Damaged Specimens
Axial strain ε (%) 10 cycles Figure 6. Axial stress-strain curves of pre-damaged specimens with σ 3 = 50 MPa. Table 3 lists the average uniaxial compressive strength of pre-damaged specimens under different damage conditions, and its relationship with the cycle numbers and confining pressures is shown in Figure 7.  The uniaxial compressive strength of pre-damaged specimens with less than 20 cycles decreased with the increase of cycle numbers under different confining pressures ( Figure 7a). Due to cyclic loading and unloading, cracks were generated inside the specimen, which reduced the bearing capacity of the specimen. Therefore, the uniaxial com-  The uniaxial compressive strength of pre-damaged specimens with less than 20 cycles decreased with the increase of cycle numbers under different confining pressures ( Figure 7a). Due to cyclic loading and unloading, cracks were generated inside the specimen, which reduced the bearing capacity of the specimen. Therefore, the uniaxial compressive strength of the specimen was weakened. However, the pre-damaged specimens with more than 20 cycles showed complex strength laws. Under the confining pressures of 50 MPa and 80 MPa, the uniaxial compressive strength of the pre-damaged specimens increased with the increase of cycle numbers and decreased with the increase of cycle numbers at 100 MPa and 120 MPa. In the high cyclic loading and unloading numbers (more than 20 cycles) of pre-damaged specimen, a large number of cracks inside the specimen was initiated and accumulated. Under the confining pressures of 50 MPa and 80 MPa, it is shown that the uniaxial compressive strength of the pre-damaged granite specimen did not always decrease with the increase of the number of cycles. This result differs from that revealed in He et al. [27]. This difference may be due to Griffith cracks, which affect the strength of the specimen according to its size. It has also been shown that if the deviatoric stress of cyclic loading is lower than the fatigue threshold (generally 94% of rock peak strength), the peak strength obtained by post-monotonic loading may increase [28]. In this test, the upper limit of deviatoric stress cycle was about 93% of the peak strength of granite, so a slight increase in strength was also possible. Under the confining pressure of 100 MPa and 120 MPa, due to the excessive confining pressure, the closed damage cracks inside the specimen started and propagated again, so the uniaxial compressive strength of pre-damaged specimen would continue to weaken. In general, the uniaxial compressive strength of pre-damaged granite specimens is much lower than that of intact granite specimens, which further proves that the bearing capacity of rock decreases after fatigue damage.

Strength Characteristics of Pre-Damaged Specimens
The relationship between uniaxial compressive strength and confining pressure is nonlinear (Figure 7b). Under different cycle numbers (10 cycles, 20 cycles, and 30 cycles), the uniaxial compressive strength of pre-damaged specimens decreased with the increase of confining pressure, which decreased by 28.9%, 14.8%, and 52.9%, respectively. This shows that the restriction effect of confining pressure on the internal microcrack evolution of the specimen was weakened with the increase of cycle numbers, and the damage of the specimen was mainly dominated by cycle numbers.

Relationship of Damage Variables to Cycle Numbers and Confining Pressures
Cyclic damage of rocks refers to the existence of internal cracks, micro-holes, and other defects in rocks in a certain pressure field when subjected to a certain cyclic load, resulting in micro-stress concentrations at their ends, and the interaction between microdefects results in the activation of primary cracks or holes and dislocation between rock particles [29]. The meso-damage evolution within the rock is related to cycle numbers and confining pressure, among which the existence of initial confining pressure restricts the generation of micro-defects and particle dislocation. The damage variable D is the description of the meso-damage evolution within the rock [30]. Therefore, the damage variable D for pre-damaged specimens in this paper is defined as: where σ c is the uniaxial compressive strength of the pre-damaged specimen, σ c0 is the uniaxial compressive strength of the intact specimen, and D is the damage variable of the pre-damaged specimen. It is generally believed that when D = 0, the rock is in a nondestructive state, at which time σ c0 = σ c ; when 0 < D < 100, the rock is in the process of deformation failure and is damaged to varying degrees. When D = 100, it reaches rock material failure, corresponding to a total loss condition. According to Equation (2), the damage law of pre-damaged specimens under different cycle numbers and confining pressures is obtained (Figure 8).
pre-damaged specimens increased, and the magnitude of the damage changes varied with different confining pressure. At 50 MPa and 80 MPa, the increasing rate of the specimen's damage gradually slows down, but there is a certain degree of damage in each cycle. At 100 MPa and 120 MPa, the increasing rate of the specimen's damage accelerated with the increase of cycle numbers, with an increase of 167% and 87%, respectively. The damage characteristics of pre-damaged specimens were not only related to cycle numbers, but also related to the confining pressure. With the increase of confining pressure, the damage variable of pre-damaged specimens after the same number of cycles became higher. When the specimen was subjected to a high number of cyclic loads, the confining pressure effect was not obvious.  According to Section 3.1, the uniaxial compressive strength of the intact specimen σ c0 was 164.88 MPa. In Figure 8, as cycle numbers increased, the damage variables for the pre-damaged specimens increased, and the magnitude of the damage changes varied with different confining pressure. At 50 MPa and 80 MPa, the increasing rate of the specimen's damage gradually slows down, but there is a certain degree of damage in each cycle. At 100 MPa and 120 MPa, the increasing rate of the specimen's damage accelerated with the increase of cycle numbers, with an increase of 167% and 87%, respectively. The damage characteristics of pre-damaged specimens were not only related to cycle numbers, but also related to the confining pressure. With the increase of confining pressure, the damage variable of pre-damaged specimens after the same number of cycles became higher. When the specimen was subjected to a high number of cyclic loads, the confining pressure effect was not obvious.

Damage Evolution Mechanism of Pre-Damaged Specimens
To further analyze the evolution mechanism (initiation, propagation, and coalescence) of surface cracks in pre-damaged granites, representative specimens were selected and analyzed in detail with the digital speckle correlation method. The surface principal strain cloud maps of the pre-damaged specimen under uniaxial compression are illustrated in Figures 9-11, in which different colors represent different principal strain regions; blue represents low principal strain regions and red represents high principal strain regions. From the principal strain cloud map corresponding to the marks on the stress-strain curve, it can be intuitively seen that there is a one-to-one correspondence between the macroscopic mechanical behavior and the principal strain field during the loading process of pre-damaged specimens.
curve is concave, obviously showing the phenomenon of strain increase. Entering the plastic phase means that the rock will have irreversible plastic deformation. Finally, there is the failure stage of the rock. It should be noted that the plastic stage of the granite specimens after pre-damage in this test was very short, and in some cases almost non-existent, as shown in Figures 10b and 11a, which exhibit the characteristics of plastic-brittle failure [31]. In order to combine the processed speckle images, we took the typical representative images of each stage of the granite specimens for analysis.  was not obvious. Because the surface expanded by the tensile crack was peeled off, the speckle surface could not be captured by the system during image processing, which lead to the failure of recognition and calculation. In the whole process of crack evolution, tensile crack was generated first, then shear crack was generated. In the final failure of the specimen, the tensile crack and shear crack were independent of each other, and there was no coalescence between the cracks. This damage evolution was concentrated in the prepeak stage, which started from the initial elastic stage to generate cracks, and then propagated crack in the elastic stage and plastic stage, until the specimen was destroyed at the peak, suggesting that the overall evolution was relatively slow.
At 80 MPa and 10 cycles, the stress at points a-h were 11.8%σc, 52.8%σc, 70.5%σc, 93.8σc, 95.6%σc, 100%σc, 89.1%σc, and 86.7%σc, respectively (Figure 9b). Point a corresponds to the compaction stage of the curve, at which time there was a small strain concentration in the upper left corner of the specimen. Points b and c correspond to the elastic stage of the curve, and when loaded to point c, a tensile band began to form at the upper right side of the specimen. Points d and e correspond to the plastic stage of the curve; when loaded to point d, the stretch band in the upper right corner deepened further, while the shear band in the upper left corner did not change significantly. When loaded to point e, a wide range of shear bands were formed in the upper left corner of the specimen. Points f-h correspond to the failure stage of the curve. When loaded to point f, cracks on the upper left corner of the specimen propagated and penetrate, leading to fracture. After point f, the stress drops rapidly. When the stress drops to points g and h, the specimen cracks propagated rapidly and coalesced into and penetrated each other, eventually destroying a large area in the upper left corner of the specimen. During the whole evolution, numerous tensile cracks and partial shear cracks were initiated, all of which propagated along the direction of compressive stress, and the various cracks interacted with each other, eventually penetrating to form stepped shear cracks, on one side of which the stepped shear cracks generally rupture. This is similar at 100 MPa and 120 MPa confining pressures. At 50 MPa and 20 cycles, the stress at points a-h were 2.8%σc, 12.2%σc, 70.8%σc, 81.2%σc, 95.8%σc, 100%σc, 97.5%σc, and 98.9%σc, respectively (Figure 10a). Point a and The deformation characteristics of the granite specimen after pre-damage show four stages. First, the compaction stage, which has an upward-concave stress-strain curve, is more significant compared to the undamaged granite specimens of group D ( Figure 5). This is followed by the elastic stage; the stress-strain curve in this stage is basically straight, as shown in the red line part in Figure 9a. Next is the plastic stage, when the curve is concave, obviously showing the phenomenon of strain increase. Entering the plastic phase means that the rock will have irreversible plastic deformation. Finally, there is the failure stage of the rock. It should be noted that the plastic stage of the granite specimens after pre-damage in this test was very short, and in some cases almost non-existent, as shown in Figures 10b and 11a, which exhibit the characteristics of plastic-brittle failure [31]. In order to combine the processed speckle images, we took the typical representative images of each stage of the granite specimens for analysis. connected with each other to form two crossed macroscopic shear cracks, namely quasicoplanar shear cracks. This damage evolution was concentrated near the peak, and it started to initiate cracks in the middle elastic stage, then crack propagation and penetration occurred in the plastic failure stage, and finally destroyed the specimen. The evolution of cracks under other confining pressures (Figure 10b) is similar to Figure 9b in that tensile cracks were initiated first, and then the interactions between the tensile cracks merged to form stepped shear cracks. At 50 MPa and 30 cycles, the stress at points a-h were 20.2%σc, 67.6%σc, 86.9%σc, 100%σc, 90.4%σc, 90.2%σc, 86.9%σc, and 72.6%σc, respectively (Figure 11a). Point a corresponds to the compaction stage of the curve, points b and c correspond to the elastic stage of the curve, and the deformation at points a and b is uniform. When loaded to point c, a strain concentration began to form in the lower left corner of the specimen. Points d-h correspond to the failure stage of the curve. When loaded to points d-f, an obvious tensile band formed on the left side of the specimen and began to deform unevenly. When loaded to points g-h, a shear crack through the specimen formed, which interacted with the tensile crack, forming a "Y" type quasi-coplanar shear crack. This damage evolution was concentrated in the post-peak stage, and the evolution speed was rapid.
At 80 MPa and 30 cycles, the stress at points a-h were 3.8%σc, 64.2%σc, 96.0%σc, 98.1%σc, 99.1%σc, 100%σc, 89.3%σc, and 73.6%σc, respectively (Figure 11b). Point a corresponds to the compaction stage of the curve, and the deformation was uniform. Point b corresponds to the initial elastic phase of the curve. When loaded to point b, the deformation was still uniform. Points c-e correspond to the plastic stage of the curve. When loaded to point c, a stress concentration began to appear on the left side of the specimen with a slight rupture. When loaded to points d-e, the strain area on the left side of the At 50 MPa and 10 cycles, the stress at points a-h were 0.6%σ c , 1.6%σ c , 23.0%σ c , 37.8% σ c , 72.6%σ c , 85.7%σ c , 88.6%σ c , and 100%σ c , respectively (Figure 9a). Points a, b, c, and d correspond to the compaction stage of the curve; the specimen deformation was small at this time. When loaded to point c, an obvious tensile zone was formed in the upper left corner of the specimen, and crack initiation began. Points e, f, and g correspond to the elastic stage of the curve. When loaded to point d-g, the tensile band in the upper left corner extended further along the direction of compressive stress until it reached both ends of the specimen. Point h corresponds to the failure stage of the curve, and it is the peak stress point. When loaded to point h, a penetrating shear band rapidly formed in the upper right corner of the specimen. At this time, the tensile crack zone on the left side of the specimen was not obvious. Because the surface expanded by the tensile crack was peeled off, the speckle surface could not be captured by the system during image processing, which lead to the failure of recognition and calculation. In the whole process of crack evolution, tensile crack was generated first, then shear crack was generated. In the final failure of the specimen, the tensile crack and shear crack were independent of each other, and there was no coalescence between the cracks. This damage evolution was concentrated in the pre-peak stage, which started from the initial elastic stage to generate cracks, and then propagated crack in the elastic stage and plastic stage, until the specimen was destroyed at the peak, suggesting that the overall evolution was relatively slow.
At 80 MPa and 10 cycles, the stress at points a-h were 11.8%σ c , 52.8%σ c , 70.5%σ c , 93.8% σ c , 95.6%σ c , 100%σ c , 89.1%σ c , and 86.7%σ c , respectively (Figure 9b). Point a corresponds to the compaction stage of the curve, at which time there was a small strain concentration in the upper left corner of the specimen. Points b and c correspond to the elastic stage of the curve, and when loaded to point c, a tensile band began to form at the upper right side of the specimen. Points d and e correspond to the plastic stage of the curve; when loaded to point d, the stretch band in the upper right corner deepened further, while the shear band in the upper left corner did not change significantly. When loaded to point e, a wide range of shear bands were formed in the upper left corner of the specimen. Points f-h correspond to the failure stage of the curve. When loaded to point f, cracks on the upper left corner of the specimen propagated and penetrate, leading to fracture. After point f, the stress drops rapidly. When the stress drops to points g and h, the specimen cracks propagated rapidly and coalesced into and penetrated each other, eventually destroying a large area in the upper left corner of the specimen. During the whole evolution, numerous tensile cracks and partial shear cracks were initiated, all of which propagated along the direction of compressive stress, and the various cracks interacted with each other, eventually penetrating to form stepped shear cracks, on one side of which the stepped shear cracks generally rupture. This is similar at 100 MPa and 120 MPa confining pressures.
At 50 MPa and 20 cycles, the stress at points a-h were 2.8%σ c , 12.2%σ c , 70.8%σ c , 81.2%σ c , 95.8%σ c , 100%σ c , 97.5%σ c , and 98.9%σ c , respectively (Figure 10a). Point a and point b correspond to the compaction stage of the curve, at which time the deformation was relatively uniform. Points c-d correspond to the elastic stage of the curve, and when loaded to point c, a strain concentration occurred in the upper left corner of the specimen and the specimen began to deform unevenly. Point e corresponds to the plastic stage of the curve. When loaded to point e, multiple tensile bands initiated above the specimen, all of which propagated along the direction of compressive stress. Points f-h correspond to the failure stage of the curve, and when loaded to points f-h, the tensile crack terminals connected with each other to form two crossed macroscopic shear cracks, namely quasi-coplanar shear cracks. This damage evolution was concentrated near the peak, and it started to initiate cracks in the middle elastic stage, then crack propagation and penetration occurred in the plastic failure stage, and finally destroyed the specimen. The evolution of cracks under other confining pressures (Figure 10b) is similar to Figure 9b in that tensile cracks were initiated first, and then the interactions between the tensile cracks merged to form stepped shear cracks.
At 50 MPa and 30 cycles, the stress at points a-h were 20.2%σ c , 67.6%σ c , 86.9%σ c , 100%σ c , 90.4%σ c , 90.2%σ c , 86.9%σ c , and 72.6%σ c , respectively (Figure 11a). Point a corresponds to the compaction stage of the curve, points b and c correspond to the elastic stage of the curve, and the deformation at points a and b is uniform. When loaded to point c, a strain concentration began to form in the lower left corner of the specimen. Points d-h correspond to the failure stage of the curve. When loaded to points d-f, an obvious tensile band formed on the left side of the specimen and began to deform unevenly. When loaded to points g-h, a shear crack through the specimen formed, which interacted with the tensile crack, forming a "Y" type quasi-coplanar shear crack. This damage evolution was concentrated in the post-peak stage, and the evolution speed was rapid.
At 80 MPa and 30 cycles, the stress at points a-h were 3.8%σ c , 64.2%σ c , 96.0%σ c , 98.1%σ c , 99.1%σ c , 100%σ c , 89.3%σ c , and 73.6%σ c , respectively (Figure 11b). Point a corresponds to the compaction stage of the curve, and the deformation was uniform. Point b corresponds to the initial elastic phase of the curve. When loaded to point b, the deformation was still uniform. Points c-e correspond to the plastic stage of the curve. When loaded to point c, a stress concentration began to appear on the left side of the specimen with a slight rupture. When loaded to points d-e, the strain area on the left side of the specimen gradually expanded to form multiple shear cracks, and the strain area was concentrated toward the center of the specimen. Points f-h correspond to the failure stage of the curve. Point f was the peak stress point, and when the stress dropped to point h, the specimen strain zone exhibited an "X" type quasi-coplanar shear crack. Under the confining pressure of 100 MPa and 30 cycles, the crack evolution process was similar to that of the stepped crack formed above. Under the confining pressure of 120 MPa and 30 cycles, the crack evolution process was a typical quasi-coplanar shear crack evolution process, but its evolution began very early.
DSCM results showed that the surface crack initiation period of the specimen was related to the cycle numbers. The larger the cycle numbers, the surface crack initiation period changed from the elastic stage to the plastic stage, respectively. During the uniaxial compression failure process of pre-damaged specimens, the final macroscopic failure of all specimens was the result of the continuous accumulation of early damage [32]. Due to the influence of the boundary, stress concentration phenomenon usually occurred at the end of the specimen, so the crack initiation was mostly concentrated at the end and always expanded along the direction of the maximum principal strain, leading to the occurrence of early damage. In addition, tensile crack was always the earliest crack in the pre-damaged granite specimen, and the initiation of secondary cracks around the tensile crack was mostly the result of the joint action of tension and shear. The shear gradually played a leading role in the crack propagation, and the crack evolution presented a localized phenomenon. These results are consistent with those of crack evolution in [16].

Failure Mode of Pre-Damaged Specimens
Tensile failure is caused by the transverse tensile stress of the failure surface exceeding the maximum tensile stress of the rock, and the crack is called tensile crack. The path of propagation of the tensile crack is smooth. Shear failure is caused by the shear stress of the failure surface exceeding the maximum shear stress of the rock, and the crack is called shear crack. The surface of a shear crack is considerably rougher than that of a tensile crack [33][34][35].
The pre-damaged specimens showed different failure modes under uniaxial compression. Based on the two basic failure modes of tensile crack and shear crack, the failure modes of the pre-damaged specimen can be divided into three categories: (1) tensile shear failure (mode I), (2) quasi-coplanar shear failure (mode II), and (3) stepped path failure (mode III) ( Figure 12). (1) Mode I: The main cracks in this failure mode are tensile cracks and shear cracks, and these two main cracks are independent of each other at the time of final failure. No crack coalescence occurs, so this failure mode is called tensile shear failure ( Figure  12a). The red line is a tensile crack, the blue line is a shear crack, and the gray line is secondary shear cracks in Figure 12a. The final failure in this mode will cause spalling on the tensile crack or shear crack surface, and the overall degree of damage is small. This failure mode is also identified in [33,36], which is called mixed coalescence mode. (2) Mode II: The main cracks in this failure mode are two intersecting shear cracks, so this mode is called quasi-coplanar shear failure (Figure 12b). The blue line is a quasicoplanar shear crack, and the gray line is a shear secondary crack in Figure 12b. The damage in this mode is generally marked by peeling from its conjugate surface or collapse, and the degree of damage is large. Such quasi-conjugate or single shear plane failure results have been reported in [37]. (3) Mode III: The main cracks in this failure mode are tensile cracks. Because multiple tensile cracks interact, propagate and coalescence to form a shear crack in a stepped path, this mode is called a stepped path failure (Figure 12c). The purple line is a step path crack, and the blue line is a tensile crack in Figure 12c. The main feature of this (1) Mode I: The main cracks in this failure mode are tensile cracks and shear cracks, and these two main cracks are independent of each other at the time of final failure. No crack coalescence occurs, so this failure mode is called tensile shear failure (Figure 12a). The red line is a tensile crack, the blue line is a shear crack, and the gray line is secondary shear cracks in Figure 12a. The final failure in this mode will cause spalling on the tensile crack or shear crack surface, and the overall degree of damage is small. This failure mode is also identified in [33,36], which is called mixed coalescence mode. (2) Mode II: The main cracks in this failure mode are two intersecting shear cracks, so this mode is called quasi-coplanar shear failure (Figure 12b). The blue line is a quasicoplanar shear crack, and the gray line is a shear secondary crack in Figure 12b. The damage in this mode is generally marked by peeling from its conjugate surface or collapse, and the degree of damage is large. Such quasi-conjugate or single shear plane failure results have been reported in [37]. (3) Mode III: The main cracks in this failure mode are tensile cracks. Because multiple tensile cracks interact, propagate and coalescence to form a shear crack in a stepped path, this mode is called a stepped path failure (Figure 12c). The purple line is a step path crack, and the blue line is a tensile crack in Figure 12c. The main feature of this failure mode is that there are generally multiple groups of parallel tensile cracks above or below the step. Failure mode III in this mode generally shows layers of peeling in the tensile crack area, and the degree of damage is relatively large. This result is consistent with those of a granite specimen that contains multiple pre-existing holes in [33].
Both cycle numbers and confining pressure affect the failure mode of pre-damaged specimens to varying degrees. Table 4 summarizes failure modes of pre-damaged specimens under uniaxial compression, with red lines representing tensile cracks, blue lines representing shear cracks, and purple lines representing stepped cracks. When confining pressure is 50 MPa, pre-damaged specimens with 10 cycles belong to mode I, and specimens with 20 and 30 cycles belong to mode II. When confining pressure is 80 MPa, predamaged specimens with 10 cycles belong to the transition failure mode I-III, specimens with 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. When confining pressure is 100 MPa, all pre-damaged specimens belong to mode III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II. III. When confining pressure is 120 MPa, pre-damaged specimens with 10 and 20 cycles belong to failure mode III, and specimens with 30 cycles belong to failure mode II.

Conclusions
(1) The mechanical properties of granite specimens that undergo high cycle pre-damage are significantly different from those of their low cycle pre-damage specimens. The mechanical properties of the high cycle pre-damaged specimens decreased significantly. Under uniaxial compression, the pre-damaged specimens showed complex strength characteristics, and uniaxial compressive strength presents a non-

Conclusions
(1) The mechanical properties of granite specimens that undergo high cycle pre-damage are significantly different from those of their low cycle pre-damage specimens. The mechanical properties of the high cycle pre-damaged specimens decreased significantly. Under uniaxial compression, the pre-damaged specimens showed complex strength characteristics, and uniaxial compressive strength presents a non-

Conclusions
(1) The mechanical properties of granite specimens that undergo high cycle pre-damage are significantly different from those of their low cycle pre-damage specimens. The mechanical properties of the high cycle pre-damaged specimens decreased significantly. Under uniaxial compression, the pre-damaged specimens showed complex strength characteristics, and uniaxial compressive strength presents a non-

Conclusions
(1) The mechanical properties of granite specimens that undergo high cycle pre-damage are significantly different from those of their low cycle pre-damage specimens. The mechanical properties of the high cycle pre-damaged specimens decreased significantly. Under uniaxial compression, the pre-damaged specimens showed complex strength characteristics, and uniaxial compressive strength presents a non-Failure mode II

Conclusions
(1) The mechanical properties of granite specimens that undergo high cycle pre-damage are significantly different from those of their low cycle pre-damage specimens. The mechanical properties of the high cycle pre-damaged specimens decreased significantly. Under uniaxial compression, the pre-damaged specimens showed complex strength characteristics, and uniaxial compressive strength presents a non-linear relationship with cycle numbers and confining pressures. Damage to pre-damaged specimens was primarily due to the strength weakening effect caused by cycle numbers. (2) The damage variables of pre-damaged specimens positively correlate with the cycle numbers and confining pressure. Damage to the pre-damaged specimens continues to increase with increasing cycle numbers, and the magnitude of damage varied with various confining pressure. The greater the confining pressure, the greater the damage variable in pre-damaged specimens. When specimens were subjected to high cycle numbers, the confining pressure restriction effect was not obvious. (3) The evolution laws of uniaxial compression damage propagation in pre-damaged specimens showed differences and obvious localization phenomenon under different cycle numbers and confining pressures. The number of initial cracks was less at low cycle numbers. However, the number of initial cracks was large under high cycle numbers, and they coalesced with each other to form stepped shear cracks. Tensile cracks were always the earliest cracks to appear in pre-damaged granite specimens, and the greater the cycle numbers, the more concentrated the crack initiation was in the plastic deformation stage. (4) There are three uniaxial compression failure modes for pre-damaged specimens, namely tensile shear failure (mode I), quasi-coplanar shear failure (mode II), and stepped path failure (mode III). The failure modes of pre-damaged specimens are related to confining pressure. Under different pre-damage stress environments with high confining pressures, the failure modes are dominated by Mode II and Mode III, respectively.