Failure characteristics and the law of the energy evolution of granite with different pre-crack inclination angles under uniaxial compression loading

Cracks affect the strength of rock masses and eventually threaten their stability in engineering. In order to study the fracture characteristics and mechanical properties of cracked rocks, uniaxial compression tests of pre-cracked granite samples with a central circular through hole were carried out by using MTS816 rock mechanics testing system. The inclination angles of different pre-cracks are 0°, 15°, 30°, 45°, 60°, 75°, and 90° respectively, and the influence of the crack stop hole near the crack tip on the failure behavior of pre-cracked samples is also considered. The results indicate that, compared with the intact sample, the peak strength of pre-cracked samples decreases significantly and is related to the pre-crack inclination angle. The failure mode of the sample varies with the pre-crack inclination angle, and the crack stop hole near the crack tip also has a certain influence on the crack growth to a certain extent. However, in terms of failure mode and its transformation law, the influence of central circular through hole and crack stop hole can be ignored. Generally speaking, the larger the inclination angle of the pre-crack, the more the total energy required for failure of the sample, and the more the stored elastic strain energy. Before the peak strength, the elastic strain energy of the sample is greater than the dissipated energy, after the peak strength, the dissipated energy gradually exceeds the elastic strain energy due to energy conversion. It is found that the pre-crack reduces the energy storage capacity of the sample, and the total energy is ultimately dominated by sample integrity. The dissipated energy rate increases first, then decreases, and finally increases again, the inflection points are the end of micro-crack closure and the peak strength, respectively. The crack stop hole changes the law of energy evolution to a certain extent, which can improve the ability of rocks to accumulate energy when designed at an appropriate position, so as to improve its load-bearing capacity in a certain range. The results display the mechanical properties of pre-cracked granite samples under uniaxial compression and are conducive to its application in engineering.


Introduction
It is well known that the stability of rock masses is often determined by different factors , in nature, the internal defects in rock masses, such as micro-cracks, joints, intergranular pores, are significant factors affecting their mechanical properties and stability (Zhu et al., 2018;Ban et al., 2021;Chen et al., 2022;Sonmez et al., 2022;Zhao et al., 2022;Ban et al., 2023). For example, in underground roadways, tunnels and stopes, unfavorable joint surfaces may lead to rock bursts, and even threaten the safety of underground mining operations. Therefore, the study of the influence of cracks on the mechanical properties of rocks is of great significance in mining and engineering construction.
The mechanical properties of cracked rocks have been a hot topic in the field of rock mechanics for at least half a century. Advances in the mechanics of fracture development have been paralleled by related yet largely independent work on the statistical physics of fracture systems (Griffith, 1924;Schultz, 2000). In order to better study the fracture behavior of rocks, the International Society of Rock Mechanics (ISRM) has recommended some standard test methods (Xie et al., 2023), such as the cracked chevron notched Brazilian disc (CCNBD) test (Fowell, 1995), the short rod (SR) and chevron bend (CB) test (Ouchterlony, 1988), the notch semi-circular bend (NSCB) test (Kuruppu et al., 2014), etc. In recent years, many achievements have also been continuously reported. For example, Aliha and Ayatollahi (2014) observed that the average mode II fracture toughness was considerably larger than that of mode I fracture toughness. Mahdavi et al. (2020) provided comprehensive data for determining the normalized critical crack length of semicircular bend specimen with chevron notch and under pure mode I loading. In order to characterize the stress field near the crack tip under external stress, many studies on stress intensity factors have been carried out. Ayatollahi et al. (2016) calculated the minimum dimensionless stress intensity factors of semi-circular bend specimen (SCB) with straight-through and chevron notches. Combined with the finite element analysis of three kinds of Brazilian disc samples with prefabricated widths less than 2 mm, it is found that the dimensionless stress intensity factor increases with crack width increasing (Li C. et al., 2019).
Currently, some researchers have focused on the failure characteristics of different fractured rocks. Zhong et al. (2020) observed that fractures were initiated at the tips of pre-existing cracks and rapidly propagated to upper mortargranite bonded surfaces. Xiong et al. (2022) found two main fracture evolution processes of sandstone containing intersecting fissures under uniaxial compression. Based on the experimental results, Liu and Lv (2019) proposed a new model to calculate the growth path of wing cracks under compressive shear stress. Li S. et al. (2019) stated that the flaw inclination angle had little effect on the optimum index, but had a certain effect on sample fracture modes. Nara et al. (2017) stated that the crack velocity increased with increasing temperature and/or relative humidity. The crack propagation in rock sample with inclined crack under uniaxial compression had been investigated by Eftekhari et al. (2016). The progressive failure and mechanical properties of pre-cracked marble with planar or non-planar cracks were studied by using the latticespring-based synthetic rock mass (LS-SRM) modeling method (Bastola and Cai, 2019).
Moreover, digital image correlation (DIC) and acoustic emission (AE) technology were used to study the influence of bolt on fracture behavior of infilled jointed rock under uniaxial loading (Li et al., 2023). The freeze-thaw failure characteristics and strength loss of non-penetrating fractured rock mass with different fracture densities were studied by means of theoretical derivation and laboratory tests (Chen et al., 2023a). Wang et al. (2021) declared that the peak strength and elastic modulus first decrease and then increase as the crack inclination angle increases. Zhao et al. (2017) investigated the mechanical and permeability properties of fractured limestone during full-range stress-strain processes. With the development of computer technology, advanced numerical simulation technology has been widely applied in geotechnical engineering (Lei et al., 2022), taking the study of fractured rocks as an example, a grain-based distinct element model (GBM) was used to investigate the influence of grain size on the fracturing response of pre-cracked granite (Saadat and Taheri, 2019). Chen et al. (2023b) discussed the internal relations among compressionshear angle, fracture distribution and crack propagation, and the changes in mechanical properties of cracked rock samples were explored combining with the discrete element method. Haeri et al. (2015) carried out a coupled numerical-experimental analysis on crack propagation, crack coalescence and failure process of jointed rock slopes.
The deformation and failure of rocks are caused by the continuous input of external energy (Tang et al., 2021), therefore, the study of the failure process of rocks from the perspective of energy is conducive to exploring the mechanical mechanism of its failure. Yu et al. (2022) discussed the relationship among the strain energy, volume fraction, confining pressure and failure mode of the backfill-encased-rock through triaxial compression test, acoustic emission, and microscopic tomography. Lin H. et al. (2022) concluded that the dissipated energy difference increases with the increase of confining pressure when the strain increases to a certain value. Gong et al. (2022a) and Gong et al. (2022b) put forward and verified the damage constitutive model of brittle rock based on linear energy dissipation law. Xu et al. (2021) found that the total input energy density, elastic strain energy density, and dissipated energy density of the red sandstone with five different fracture dip angles show obvious quadratic growth in response to the actual unloading stress. There is a good consistency between the experimentally obtained fracture loads and the theoretical predictions based on the constancy of the mean strain energy density over the material volume (Razavi et al., 2018). Bertram and Kalthoff (2003) found that the energy dissipated by a crack propagating at a high velocity was almost two orders of magnitude larger than the energy to initiate a crack. The relationship between energy dissipation and fracture mechanism in graphite fracture zone was evaluated by a tri-linear softening model (Tang et al., 2021).
The above discussion indicates that the mechanical properties of cracked rocks have been one of the hot topics concerned by scholars all over the world, and their achievements also provide an important reference for this study. In this experiment, the pre-cracked granite samples with seven different inclination angles were prepared, different Frontiers in Earth Science frontiersin.org from the traditional oblique straight pre-crack, a central circular through hole with a diameter of 5 mm was drilled at the center of the sample. In addition, for samples with pre-crack inclination angle of 45°, the influence of the crack stop hole near the crack tip is considered. In brief, we discussed the failure behavior of precracked samples under uniaxial compression, and analyzed the laws of their energy evolution. The results verify the mechanical properties of this type of granite, which is conducive to its application in engineering.

Test sample and equipment
Experiment is an important research method in the field of geotechnical engineering (Kang et al., 2022;Lin S. Q. et al., 2022;Wang et al., 2022;Zhao et al., 2023). In this study, the common granite was selected as the research object. The granite samples obtained after processing have no obvious joints and cracks, and have a very good compositional homogeneity. Its shape is a thin cuboid with length, width and thickness equal to 160, 80, and 10 mm respectively. In addition, there is a circular through hole with a diameter of 5 mm in the center of the sample, an oblique straight pre-crack passes through the hole, and their centers coincide with each other. The length and width of the pre-crack are 10 and 1 mm respectively. The pre-crack inclination angles α are 0°, 15°, 30°, 45°, 60°, 75°, and 90°respectively. The left side of Figure 1 shows a typical sample, its pre-crack inclination angle α is equal to 45°.
In order to study the influence of crack stop holes on the failure behavior of granite samples, for some samples with pre-crack inclination angle equal to 45°, a crack stop hole was drilled near each end of the pre-crack for comparative analysis. The function of the crack stop hole is to prevent further propagation of cracks in areas where stress is easily concentrated . Nevertheless, the role of crack stop hole in this study is to change the crack growth path, therefore, it is arranged near the route of wing crack propagation, as shown in Figure 1. The vertical distance d 1 from the center of the crack stop hole to the tip of the pre-crack is 10 mm, while the corresponding horizontal distance d 2 from their center to the corresponding pre-crack tip is 0, 5, and 10 mm respectively, and the radius of the crack stop hole is 2 mm. Figure 1 shows two typical crack stop holes with a horizontal distance of 5 mm and a vertical distance of 10 mm.
Uniaxial compression test was carried out by MTS816 rock mechanics testing system, as shown in Figure 2. During the test, the axial stress was applied to samples by displacement control, and the loading rate was 0.002 mm/s. The system can automatically acquire parameters such as time, axial load and displacement, etc., and the collection interval was set to 0.1 s.

Principle of energy calculation
Energy evolution plays an important role in the analysis of mechanical characteristics, there are mainly four types of energy changes in the process of rock failure under external load, namely, energy input, energy accumulation, energy dissipation and energy release. According to the laws of thermodynamics, based on the assumption that there is no heat exchange with the external environment during loading, and ignoring the energy released through thermal radiation, heat exchange and other forms, the total energy input from the outside is completely converted into elastic strain energy and dissipated energy inside the sample (Peng et al., 2020). Therefore, the energy of the pre-cracked granite sample under uniaxial compression is calculated as follows Qiao et al., 2022). Frontiers in Earth Science frontiersin.org where U is the total energy input from the outside, MJ·mm −3 ; U e is the elastic strain energy, MJ·mm −3 ; U d is the dissipated energy, MJ·mm −3 . In the principal stress space, the total energy and elastic strain energy absorbed by sample element can be calculated by the following formula : where σ i , σ j , σ k (i, j, k =1, 2, 3) are the principal stresses in each direction; ε i , ε e i are the total strain and elastic strain in each principal stress direction; μ i , E i are Poisson's ratio and elastic modulus during the unloading stage, which can be replaced by μ 0 , E 0 at elastic stage (Xie et al., 2008).
In this study, it should be noted here that the elastic modulus of pre-cracked granite samples should be called apparent stiffness. Under uniaxial compression, the circumferential stress is 0 MPa, and only the axial stress does work. Therefore, according to Eqs 2-4, the energy can be calculated as follows Qiao et al., 2022): where σ 1(i) , ε 1(i) are the stress and strain values corresponding to each point on the stress-strain curve. According to Eqs 1, 5, 6, the dissipated energy of the sample is the difference between the total energy and the elastic strain energy (Peng et al., 2020).

Variation of mechanical parameters
It can be seen from Figure 3A that the trend of stress-strain curves of samples with different pre-crack inclination angles is basically consistent with that of samples without cracks, and they all go through four stages, namely, micro-crack compaction stage, elastic deformation stage, unstable fracture development stage and post failure stage. In addition, the peak strength of the precracked samples is always lower than that of the sample without pre-cracks, indicating that pre-crack reduces the peak strength of the granite sample, and make it more vulnerable to damage.
As shown in Figure 3B, for granite samples with pre-crack inclination angle of 45°, crack stop hole has a certain influence on the stress-strain curve. When the horizontal distance d 2 between the center of the crack stop hole and the tip of the pre-crack is equal to 0 mm, the peak strength of the sample is slightly lower than that of FIGURE 2 MTS816 rock mechanics testing system and its main components.

Frontiers in Earth Science
frontiersin.org the one without crack stop hole. However, when the horizontal distance d 2 is equal to 5 and 10 mm, the peak strength is increased slightly. Moreover, the overall strain decreases under the influence of crack stop hole, indicating that the crack stop hole reduces the ductility of the sample before failure. Although there are differences in the position of crack stop holes, the distance between curves is closer than that of the sample without crack stop holes, and the slope of the stress-strain curve is larger than that of the sample without crack stop holes, indicating that crack stop holes increase the apparent stiffness of granite samples. Finally, in Figures 3A, B,  Frontiers in Earth Science frontiersin.org the stress-strain curve drops rapidly after the peak strength, reflecting the brittle fracture characteristics of granite. Figure 4 shows the peak strength of different pre-cracked granite samples, it can be seen that the peak strength increases with the increase of the inclination angle, as shown in Figure 4A. When the inclination angle of the pre-crack is 0°, the peak strength of the sample reaches the lowest value equal to 49.22 MPa, and when the inclination angle rises to 90°, the peak strength increases to 151.14 MPa. According to fracture mechanics theory and Figure 5, the stress field near the crack tip in polar coordinate system can be obtained , as shown in Eq. 8.
where σ r is the radial stress; σ θ is the circumferential stress; τ rθ is the shear stress; θ is the angle in polar coordinate system; r is the distance from the crack tip; K I and K II are the stress intensity factors of type I and II cracks, respectively. Under uniaxial compression, the pre-cracked granite sample is only subjected to axial stress. Therefore, for an open crack, K I and K II are determined by Eq. 9.
where a is half the length of the long axis of the crack and equal to the distance from the crack tip to the crack center; β is the angle between crack strike and the direction of axial stress, as shown in Figure 1B. According to Eq. 9, when β is equal to 90°, that is, the precrack inclination angle α (see Figure 1) is 0°, in this case, K I obtains the maximum value and K II = 0, therefore, it only shows the characteristics of type I crack. Considering that the tensile strength of rocks is far lower than the shear strength and compressive strength, the granite sample with pre-crack inclination angle of 0°does have the lowest peak strength. As pre-crack inclination angle increases, K I gradually decreases, K II is no longer zero, and the sample shows the characteristics of both type I and II cracks. Moreover, the width of the pre-crack in this study is 1 mm, according to the results obtained by Zhang et al. (2020), when the crack width is greater than 0.8 mm, the peak strength decreases with the increase of β, and the minimum value is obtained when β = 90°. When β = 0°, the pre-crack inclination angle α is 90°, both K I and K II are equal to zero, indicating that there is no stress concentration at the crack tip, in this case, the pre-cracked granite sample obtains the highest value of peak strength. Combined with the experimental results in this paper, it can be seen that the theoretical results mentioned above have been verified again.
For pre-cracked samples with crack stop holes, the peak strength changes with the horizontal distance d 2 between the crack stop hole and the crack tip from 0 to 10 mm, as shown in Figures 3B, 4B. In detail, the peak strength is 80.68 MPa when the horizontal distance d 2 is 0 mm, which is smaller than the peak strength of the sample without crack stop holes (87.45 MPa). The peak strength increases to 105.33 MPa when the horizontal distance d 2 is equal to 5 mm, and decreases to 101.02 MPa when d 2 is equal to 10 mm. It is estimated that the crack stop hole changes the law of crack propagation, resulting in changes in stress distribution. Therefore, there is a reasonable location of the crack stop hole to obtain the extreme value of the peak strength of the pre-cracked granite sample. In summary, the arrangement of the crack stop hole plays an important role in improving the stability of the structure of cracked rocks.
The peak strength of pre-cracked granite samples under uniaxial compression is positively and linearly correlated with the inclination angle of the pre-crack, and the fitting equation is as follows: σ p 1.04α + 45.54 R 2 0.92 where σ p represents the peak strength of the pre-cracked granite samples, MPa; α is the pre-crack inclination angle,°, as shown in Figure 1B. The ratio of the peak strength of the pre-cracked sample to the peak strength of the intact sample is defined as strength weakening coefficient to characterize the weakening degree of the peak strength of granite caused by pre-crack inclination angle, as shown in Figure 4A. When the pre-crack inclination angle is 0°, the peak strength of the sample is only 26% of that of the intact sample. When the pre-crack inclination angle increases from 0°to 60°, the strength weakening coefficient increases from 0.26 to 0.46, indicating that the peak strength of the pre-cracked granite sample increases relatively slowly. However, when the pre-crack inclination angle increases from 60°to 90°, the peak strength increases relatively quickly, for example, the peak strength is 78.73% of that of intact sample when the inclination angle is 90°.
The slope of the straight line segment of the stress-strain curve is defined as the elastic modulus of the rock. However, considering that the elastic modulus is only related to the performance of the rock, such as lithology, composition, particle size, water content and pore size, therefore, the so-called elastic modulus here should actually be regarded as apparent stiffness. According to the

FIGURE 5
Stress components at the crack tip in polar coordinate system.

Frontiers in Earth Science
frontiersin.org stress-strain curve in Figure 3, the apparent stiffness of different precracked samples can be obtained, as shown in Figure 6. It can be seen that the apparent stiffness increases linearly with the increase of the pre-crack inclination angle. The linear fitting equation is as follows: where E t is the apparent stiffness of different pre-cracked granite samples, MPa.

Failure mode
The failure mode of the sample is a key feature reflecting its failure mechanism . Generally speaking, the failure of the sample starts from the initial cracking, and then the crack propagates until it penetrates the sample. In this study, as the pre-crack inclination angle increases, the characteristics of tensile failure gradually weaken, while the characteristics of shear failure strengthen, and shear failure plays a dominant role during this process, as shown in Figure 7A, this characteristic can be clearly seen as the inclination angle increases from 0°to 90°. For example, when α = 0°, a tensile crack initially forms and propagates on the upper edges of the central circular through hole, however, as the precrack gradually closed during the loading process, the wing crack appears at the pre-crack tip, whose growth rate is gradually faster than that of the initial tensile crack, and finally the tensile failure mode is dominant, and the tensile crack still penetrates the upper part of the sample. Moreover, the number of cracks caused by shear increases with the increase of the pre-crack inclination angle, including wing cracks and secondary cracks, especially when α = 90°, the shear failure characteristics are more significant.
It can be seen from Figure 7B that, the crack propagation of samples is significantly affected by crack stop holes. The number of cracks after sample failure is relatively minimum when the horizontal distance d 2 is equal to 0 mm, is maximum when d 2 = 5 mm, and ranks second when d 2 = 10 mm. This is closely related to the variation of peak strength of granite samples, in general, the more crack propagation, the higher the peak strength.

Energy evolution
According to Eqs 1-3, the total energy, elastic strain energy and dissipation energy of different granite samples during uniaxial compression can be calculated respectively, the results are depicted in Figures 8, 9. It can be seen from Figure 8A that the change of total energy of different samples is roughly similar, that is, the total energy absorbed by each sample increases non-linearly. Except when α = 0°, the total energy increases with the increase of the pre-crack inclination angle, indicating that the larger the pre-crack inclination angle, the more energy the sample consumes during the failure process, which reflects the compression resistance of the sample macroscopically. In Figure 8B, the change of elastic strain energy is similar to the change of stress-strain curve, showing a non-linear increase at first and then a rapid decline, that is, the elastic strain energy is released rapidly after the peak strength, resulting in rapid failure of rock samples, especially when the pre-crack inclination angle is relatively large. In addition, similar to the total energy, the elastic strain energy increases with the increase of the precrack inclination angle, except when α = 0°, indicating that the energy storage capacity of the sample increases.
In Figure 8C, dissipated energy increases slowly and linearly at first, and then rapidly and non-linearly after reaching the peak

Frontiers in Earth Science
frontiersin.org strength. The reason why the dissipated energy increases slowly before the peak strength is that most of the energy is stored as elastic strain energy. After the peak strength, the elastic strain energy is released rapidly, leading to a rapid increase in dissipated energy. It should be noted that the data of elastic strain energy and dissipated energy in the residual stage cannot be collected due to the rapid failure of the sample after the peak strength. It can be seen from Figures 9A, B that, under the influence of the crack stop hole, the total energy and elastic strain energy stored in the sample increase under the same strain. In addition, compared with the sample without crack stop hole, the total energy-strain curve and elastic strain energy-strain curve of the sample are closer to each other respectively, and the crack stop hole reduces the total strain. When d 2 = 0 mm, the total energy and elastic strain energy are smaller than those without crack stop holes. When d 2 = 5 mm and 10 mm, both types of energy are greater than those without crack stop holes. The results show that if the horizontal distance d 2 is relatively small, the elastic strain energy of the sample will be reduced, therefore, the reasonable arrangement of crack stop holes can improve the energy storage capacity of the sample.
When the peak strength is reached, the elastic strain energy reaches the maximum value. If the loading continues, the elastic strain energy is released in the form of dissipated energy, resulting in a rapid growth of dissipated energy after the peak strength, as shown in Figure 9C.

Frontiers in Earth Science
frontiersin.org strength is far more than the dissipated energy, in addition, the storage capacity of elastic strain energy of the intact sample is stronger than that of the pre-cracked samples. The pre-crack significantly reduces the total energy absorbed by the sample, which may be due to the reduction of peak strength caused by the pre-crack, thus reducing its energy storage capacity. The total energy, elastic strain energy and their proportion can be increased within a limited range by crack stop holes, and both kinds of energy are ultimately dominated by the integrity of the sample. According to Table 1, the elastic strain energy and dissipated energy of each sample at peak strength are shown in Figure 10. Generally speaking, the elastic strain energy increases non-linearly with the increase of pre-crack inclination angle, and can be expressed by Eq. 12. The dissipated energy fluctuates between 0.032 and 0.073, and its distribution can be expressed mainly by Eq. 13. It should be noted that Eq. 13 does not consider the values of points where the pre-crack inclination angle is 30°and 60°in Figure 10, otherwise the correlation coefficient will fall to a relatively low value. The reason may be that the dissipation energy is very small at peak strength, and its slight fluctuation can easily lead to a larger change in the correlation coefficient. In addition, the dispersion of granite samples increases the difference of dissipated energy. U e −0.00164α + 4.26α 2 *10 −5 + 0.116 R 2 0.93 (12) U d −0.00134α + 1.501α 2 *10 −5 + 0.0547 R 2 0.79 (13)

Change in dissipated energy rate
The mechanical behavior of rock sample under uniaxial compression can be reflected by energy dissipation . The ratio of dissipated energy to total energy is defined as the Frontiers in Earth Science frontiersin.org dissipated energy rate η e , as shown in Eq. 14, which reflects the relationship between the dissipated energy and the total energy during loading.
The variation of dissipated energy rate of each sample with axial strain is shown in Figure 11. The dissipated energy rate increases first and then decreases during the microcracks compaction stage. The reason may be that the internal microcracks of the sample are compacted at the initial stage, and the energy is released as dissipated energy. After that, if the loading continues, strain hardening gradually takes a leading position, and the energy is mainly stored in the form of elastic strain energy. When approaching the peak strength, the decline of dissipated energy rate slows down, and then the dissipated energy rate gradually increases, which may be caused by the accelerated growth of microcracks. After the peak strength, the elastic strain energy is released rapidly due to the macroscopic failure of the sample, resulting in a sharp rise in the dissipated energy rate.
As shown in Figure 11B, all curves show rising characteristics at the initial stage, however, the dissipated energy rate of the sample without crack stop hole (α = 45°) keeps stable before ε = 0.002, and begins to decline after ε is greater than 0.002, indicating that the internal microcracks of the sample without crack stop hole are relatively more developed, so its dissipated energy rate remains at a higher value. On the other hand, the crack stop hole inhibits the crack growth from the initial stage of loading, thus reducing the dissipated energy rate of samples. In addition, the difference in the location of the crack stop hole leads to the difference of the dissipated energy rate. Regardless of the well-known achievements, the new findings will be further discussed in this section. First of all, different from the traditional oblique straight pre-crack, a central circular through hole with a diameter of 5 mm was drilled at the center of the sample. Although the central circular through hole may affect the strength of pre-cracked granite samples, the law of variation of mechanical parameters has not been changed correspondingly. For example, pre-crack inclination angle has a significant influence on the peak strength, which increases with the increase of the inclination angle Saadat and Taheri, 2019;Wu et al., 2021;Ma et al., 2022), furthermore, Ma et al. (2022) pointed out that the elastic modulus (apparent stiffness) increases linearly with the increase of the pre-crack inclination angle, their results are the same as those in Figures 4A, 6A. Despite the fact that the central circular through hole does not change the law of variation of mechanical parameters, its influence on the value of mechanical parameters of pre-cracked samples is still worthy of further study.
The failure modes of pre-cracked rock samples with different inclination angle are inevitably different from each other. Wu et al. (2021) stated that, with the increase of friction angle (inclination angle), the change of crack type occurred, and more transition from tensile cracks to shear cracks can be observed. Compared with Section 3, it can be seen that the addition of the central circular through hole does not change the law of transformation of the failure mode of pre-cracked granite samples. In addition, according to Figure 7, the crack stop hole also fails to change the failure mode of the granite sample with pre-crack inclination angle of 45°, namely, the failure of the sample is dominated by both tensile and shear modes during the loading process, and the crack stop hole mainly affects the crack growth path.
Crack stop holes are always used to prevent crack propagation in areas where stress is easily concentrated, especially in the field of metal materials. For example, the experiment carried out by Fu et al. (2017) indicated that the stress at the crack tip was reduced by a crack stop hole. Since stop-holes are regarded as an emergency measure to delay crack propagation, Okura and Ishikawa (2002) presented the conditions under which stop-holes prevent the re-initiation of fatigue cracks. A model was modified based on the short crack effect and the residual stress effect to estimate the fatigue crack growth with the stop-hole retardation (Lu et al., 2020). The excellent performance of the crack stop hole has also attracted the attention of rock mechanics researchers. Rao et al. (2009) and Rao et al. (2022) stated that the adoption of optimal double holes symmetrical to the original crack plane is an effective method to prevent shear crack in brittle rocks, and pointed out that the design of stop holes is of great significance to improve the stability of mining activities by calculating the stress intensity factor of anisotropic elastic solids with multiple oval-holes and cracks.

FIGURE 10
The variation of elastic strain energy and dissipated energy at peak strength.
Frontiers in Earth Science frontiersin.org At present, the application of crack stop holes in rock mechanics is relatively few compared with the field of metal materials, and the research on crack stop holes mainly focuses on stress concentration and crack initiation (Okura and Ishikawa, 2002;Rao et al., 2009;Fu et al., 2017;Liu et al., 2017;Saadat and Taheri, 2019;Lu et al., 2020;Wu et al., 2021;Ma et al., 2022;Rao et al., 2022). In this paper, we introduce the crack stop hole into the study of granite materials, nevertheless, the focus of attention is different from the existing researches, that is, we care about the effect of the crack stop hole on apparent stiffness, peak strength, energy evolution, stress-strain curve, failure mode and shape. Therefore, the layout of the crack stop hole in this paper is inevitably different from that of previous studies, we are concerned about the location of the crack stop hole where the rock sample has the optimal mechanical properties. However, in this study, we only considered the influence of crack stop hole on granite samples with pre-crack inclination angle of 45°, in the future, the research scope can be expanded, for example, the effect of the crack stop hole on rock samples under different factors, such as different crack inclination angle, sample shape and lithology, can also be considered. Finally, numerical simulation can also be carried out in the future to further reveal the mechanical properties of rocks under the influence of central circular through hole and crack stop hole.

Conclusion
The study in this paper is beneficial to verify the mechanical properties of fractured granite rocks with different inclination angles, in which there is a circular through hole in the center of the sample, and the effect of the crack stop hole on its mechanical behavior. The main conclusions are as follows: (1) The central circular through hole does not change the law of variation of mechanical parameters of pre-cracked granite samples. With the increase of the pre-crack inclination angle, the growth of tension cracks gradually becomes less obvious, indicating that the central circular through hole does not change the law of transformation of the failure mode of pre-cracked samples. The crack stop hole also failed to change the failure mode of the pre-cracked granite sample with inclination angle of 45°, but affected the crack propagation path and its peak strength to a certain extent. (2) In general, the larger the inclination angle of the pre-crack, the more total energy required for failure of the sample and the more stored elastic strain energy, and reflecting the compressive resistance of the sample. The crack stop hole can change the law of evolution of energy to a certain extent. The pre-crack reduces the energy storage capacity of the sample, and the total energy is ultimately dominated by its integrity. (3) The law of dissipated energy rate during loading is analyzed, which implies the energy evolution mechanism of deformation and failure of granite samples. In addition, it is also found that the ability of the sample to accumulate elastic strain energy can be improved within a certain range by using crack stop holes at appropriate locations. In a word, the mechanical properties of pre-cracked granite samples are further verified in this study, thus contributing to its application in engineering. The authors declare that this study received funding from Hefei Rail Transit Group Co., Ltd (2021BFFBZ02689). The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.

Conflict of interest
Authors EZ and XZ were employed by Lanzhou Engineering and Research Institute of Non-ferrous Metallurgy Co., Ltd. Author JL was employed by Cscec Aecom Consultants Co., Ltd. Author CZ was employed by Sinosteel Maanshan General Institute of Mining Research Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.