Compression and Shear Fracture Analysis of Boundary Cracks Containing Water in Rock

In order to conserve the water resource during undergroundmining, the fracture andmechanical properties of rock are important for the stability of water-resisting layers, especially for the fracture behavior of boundary cracks containing water in rock. Considering the swelling of rock under water environment and the influence of water on rock, the stress intensity factors of modes I and II are derived for boundary cracks in rock under compressive and shear stresses. +e cracks are divided into the closed and open states. +e effects of the crack inclination angle, friction coefficient between crack surfaces, and initial crack length on stress intensity factors are also taken into account. +e stress intensity factors for closed and open boundary cracks are verified by numerical and physical experiments, respectively, and the deviation of the results is within 5%. It is shown that pore pressure has different effects on the relationship between stress intensity factor and friction coefficient under different lateral pressures. +e effect of water on crack propagation is mainly due to the deterioration of the fracture toughness of the rock. It is found that the critical coefficient λc is a key parameter to determine whether the boundary crack propagates in rock under compression-shear stress. Further studies should be performed to apply the present fracture theory to rock mass or water-resisting layers.


Introduction
Water has a significant effect on the fracture and mechanical properties of rock materials [1][2][3][4]. And a host of geological hazards and engineering disasters are closely related to water, such as landslides, karst collapse, fault activation, water inrush in mine, instability of water-resisting layer, and so on [5,6]. As coal mining moves westward in China, it has significant importance to conserve the water resource in the arid and semiarid areas of Northwest China. Due to the disturbance of underground mining or construction, the initial cracks are produced in the water-resisting strata, which are mainly composed of sandstone, shale, and clay. e water intrusion becomes an important factor that leads to the crack propagation until ultimately penetrating the whole rock stratum. In the condition of water environment, the geometric parameters of cracks and the water content will have a significant impact on the fracture and mechanical properties of rock. Yang et al. [7] performed a set of conventional triaxial compression tests on dry and brine-saturated sandstone specimens containing two preexisting three-dimensional flaws. It was shown that the peak strength and crack damage stress of the sandstone specimens increased as the flaw angle increased. Zuo et al. [8] conducted three-point bending tests on semicircular bend (SCB) sandstone specimens; the results showed that the peak load of the specimen exponentially increases with the increase of crack inclination angle. Guha Roy et al. [9] investigated the mechanical and fracture properties of sedimentary rocks, which were three types of sandstones and one type of shale saturated in water for different periods of time. e degree of saturation has significant effect on both the strength and fracture properties of sedimentary rocks. Yao [10] developed a 3D pore pressure cohesive zone model to predict hydraulic fracturing with different rock properties under fluid injection. Huang and Liu [11] studied the influence of shale bedding on the propagation of cracks containing water and summarized the propagation law of cracks for hydraulic fracture under true triaxial stresses when they meet the bedding plane.
Fracture toughness is a rock's property in the presence of a preexisting flaw, which governs the critical stress at which a crack in rock becomes unstable and propagates. e semicircular bend (SCB) specimens have been widely used for fracture toughness determination of rocks [8,[12][13][14][15][16][17][18]. Lim et al. [13,14] proved that the water saturation is a dominant factor affecting the fracture toughness of soft rocks and the fracture toughness increases with strain rate. Furthermore, Mahanta et al. [15] indicated that the fracture toughness and the energy release rate for mode I, mode II, and mixed mode are a function of strain rate. Kataoka et al. [16] investigated the influence of water vapor pressure in the surrounding environment on mode I fracture toughness of rocks. Combining semicircular bend specimens with acoustic emission monitoring, Zhou et al. [17] found that the fracture toughness and crack growth rate of sandstone would decrease with the increase of water content. For a clay-bearing sandstone, the fracture toughness decreases as relative humidity increases depending on the orientation of the microcracking and the toughness anisotropy appears during humidification [18]. Moreover, Hua et al. [19] used centrally cracked Brazilian disk specimens to study the effect of cyclic wetting and drying on the fracture toughness and crack propagation of rock. Using nanoindentation, Gupta et al. [20] measured both tensile and shear dominated fracture toughness in shales, where crack length method and energy method are suitable to mode I (tensile) and mode II (shear) fracture toughness, respectively.
In recent years, several criteria were presented for the fracture of rock, which determines whether a crack propagates or nucleates. Savitski and Detournay [21] obtained the asymptotic solution of the propagation criterion of the elliptical center crack containing fluid in the low permeable rock by analytical method. It is considered that three crack propagation regimes are identified: viscosity and fracture toughness dominated, as well as mixed regimes. e relationship between initial fracture geometry, rock compressive strength, and fracture toughness is explored from experiment and theory, which provides theoretical guidance for engineering practice [22,23]. Dong et al. [24] developed a unified interfacial crack initiation criterion expressed by the stress intensity factors, by performing three-point bending and four-point shear tests on the composite rock-concrete specimens. Sih [25] presented the K I -K II curve that governs the mixed mode fracture of cracks under remote compression. en, according to the special stress condition of underground rock mass, Zhou [26] proposed the fracture criterion of compression-shear crack of I-II mixed mode by the unilateral crack test. Based on the basic energy release rate criterion, Shen and Stephansson [27] put forward an F criterion which accords with the rock compression-shear fracture model. In order to improve the prediction accuracy of mode II stress intensity factor, Ayatollahi and Stananinia [28] proposed the maximum tangential stress criterion, which takes into account the first three terms of William's series expansion for elastic stresses at crack tip. Also, Zheng and Luo [29] deduced the criterion expression of compression-shear crack theoretically in combination with different failure criteria.
However, the current research mainly focused on the mechanical and fracture properties of central crack in rock under compression and shear state during underground mining or construction. e fracture behavior of boundary cracks in the compression-shear state has been rarely investigated in the literature, especially under the water action after the initial cracks in rock are produced due to the influence of mining or construction. And the crack opening state in rock is usually not considered strictly. erefore, the mechanical and fracture analysis of boundary cracks containing water in rock is carried out in the present study. e expressions of modes I and II stress intensity factors of the boundary cracks in the compression-shear stress state are obtained for both open and closed cracks. Subsequently, the corresponding stress intensity factors are verified by physical and numerical experiments, respectively. It provides a theoretical basis for estimating the instability law of boundary cracks in rock by compression-shear coefficient, which is a key parameter of rock.

Stress Intensity Factors for Cracks Containing Water
Under the disturbance of underground engineering such as mining, the occurrence of cracks in underground aquifer strata can be divided into the upward and downward cracks. e downward cracks are generally in the state of compression-shear stress due to the existence of in situ stress and the upper water pressure. e crack propagation can be described as the I-II compression-shear combined type. And the fracture criterion for the compression-shear rock crack can be formulated as follows [25,26,29]: where K I and K ΙΙ are the stress intensity factors of modes I and II, respectively. K ΙΙc is the mode II fracture toughness of rock and λ is the compression-shear coefficient of rock. Considering the crack opening state, it can be divided into the closed and open cracks. e criterion for a crack to be open or closed was derived by David and Zimmerman [30]. e two-dimensional model is shown in Figure 1.

Stress Intensity Factors for Closed Cracks.
Due to the influence of groundwater in the upper part of the crack, the mechanical action and the swelling effect of rock within water mainly affect the stress intensity factors at the crack tip. e mechanical action mainly includes the far-field stress, the pore pressure, and the friction between crack surfaces. erefore, the stress intensity factor can be obtained by superimposing those produced by these two factors, respectively. For the closed cracks, the stress intensity factor is expressed as 2 Advances in Materials Science and Engineering where K b I and K b ΙΙ are the stress intensity factors of modes I and II closed cracks, respectively. K p I and K p ΙΙ are the modes I and II stress intensity factors caused by the mechanical action, and K h I and K h ΙΙ are the modes I and II stress intensity factors caused by the swelling effect of rock within water.
For the stress intensity factors caused by the mechanical action, the basic form of the modes I and II closed cracks is described as follows, which is according to the stress intensity factor of the boundary crack in a semi-infinite plate obtained by using the theories of complex variable function and rational mapping function [31]: where σ n and τ p are the stresses of crack tip in the normal and tangential directions, respectively. 2a is the length of crack. e negative sign indicates that the mode I stress intensity factor hinders the tensile fracture of the crack under compression-shear condition. e analysis of crack growth depends mainly on the stress field at the crack tip. e magnitude of normal stress σ n and tangential stress τ at the boundary crack tip can be determined as Taking into account that the pore pressure caused by the existence of groundwater above the rock layer and the friction effect caused by rough surface when the closed crack surface is dislocated relatively, equations (4) and (5) change to where f is the friction coefficient of crack surface. It is found that f ∈ [0.18, 0.8], which was measured by using the electronic universal testing machine CSS-55000 and the selfmade test box. P w is the pore pressure of rocks. us, equations (6) and (7) can be substituted for equation (3). e expression of stress intensity factor for closed crack without considering the swelling is as follows: When considering the swelling effect of water on rock immersion, Liang et al. [32] performed experiments and theoretical analyses to indicate that the hydration swelling stress mainly affects the stress intensity factor of mode I crack. At this moment, K h II is equal to 0. en, the stress intensity factor of rock produced by the swelling is obtained as where P h is the hydration stress produced by the swelling, which should be measured by experiments on rocks. b is the length of the hydration stress applied to the crack surface. Finally, the mode I stress intensity factor of closed crack containing water can be expressed as At the same time, the compression-shear fracture criterion for closed cracks containing water is also obtained:

Stress Intensity Factors for Open Cracks.
Unlike closed cracks, there is no friction between crack surfaces because the cracks are not closed, and there should be a curvature ρ at the crack tip. If the curvature ρ is zero, the stress singularity at the crack tip is not affected by the shear stress at the same region, which is caused by the far-field compressive stress and the flow in the crack. If the curvature ρ is not zero, the far-field compressive stress and the flow will produce a local additional stress at the crack tip. As the swelling effect, the mode I stress intensity factor induced by the additional stress acts with K p Ι on the contrary. e open crack is regarded as a blunt crack, and the maximum additional stress at the crack tip is equal to the compressive stress produced by the far-field compressive stress, rock swelling, and water flow: e mode I stress intensity factor K T Ι caused by additional stress can be calculated by the following formula [33]: en, Finally, the expressions of modes I and II stress intensity factors for open crack are as follows: en, the compression-shear fracture criterion for open cracks containing water is obtained:

Verification of Physical Experiments for Open Cracks.
For open cracks, a set of physical experiments were performed on rock specimens with boundary cracks in different water-bearing states to verify the accuracy of the stress intensity factors. In the experiments, red sandstone was selected as the test material, and the rock specimens were processed into the prismatic specimens with the size of 60 mm × 60 mm × 120 mm. An inclined crack was prefabricated at the top center of the specimen, with a length of 30 mm and an angle of 80°with the horizontal direction. e boundary cracks were cut with a wire saw, and the crack tip was polished with steel wire. e size of the prismatic specimen is shown in Figure 2(a).
By soaking the rock, rock specimens of four water-bearing states were prepared, which are natural, water content of 2.3%, water content of 2.5%, and natural saturation states. Compared with the weight of the dried sample, the water content of the rock is about 1.8% in the natural state and 3.3% in the natural saturated state. Four specimens were made for each water-bearing state. e uniaxial compression experiments were carried out on the MTS-816 test machine. In the tests, the compression rate was 0.002 mm/s. e specimens during the tests are shown in Figures 2(b) and 2(c). Figure 3 shows the stress-strain relations of rocks with boundary open crack under four water-bearing conditions, where the stress and strain for each water-bearing state were averaged the experimental results in the same water-bearing states. It can be found that water intrusion greatly weakens the ultimate stress of rock. From the natural state to the saturated state, the ultimate stresses of rock for crack initiation decreased by 47.7%. e elastic modulus of rock slowly decreases with the increase of water content, indicating that the property of rock exhibits the trend from brittleness to plasticity under the action of water. Table 1 compares the ultimate stresses obtained from the physical experiments and the theoretical values calculated by using equation (18). e ultimate stress from test for each water-bearing state was averaged the experimental results in the same water-bearing states. e relative errors are also shown in Table 1. By comparing the experimental and theoretical results, it can be seen that the relative errors are less than 5%, indicating that the calculated ultimate stresses are in good agreement with the measured ones. Furthermore, the theoretical results are usually smaller than the experimental ones of the rock. It means that the ultimate stress of rock with open crack calculated by equation (18) is conservative and has a certain safety factor.

Verification of Numerical Experiments for Closed Cracks.
For closed cracks, taking account of the difficulty in processing closed cracks in physical experiments, numerical models are established by ABAQUS finite element analysis software. Figure 4 illustrates the numerical model of 120 mm height, 60 mm width, and 20 mm thickness. e crack length is 30 mm, with a horizontal angle of 80°. Eight-node hexahedral elements with reduced integration were adopted in the FEM. e size of element was 2 mm × 2 mm × 2 mm and the total number of the elements was 18000. e fixed displacement constraints were set at the bottom of the FEM model. e uniform pressure of 30 MPa was loaded on the top of the model. e models with boundary cracks are numerically analyzed by XFEM module. e crack tip element is selected as the object in the FEM analysis. Due to the existence of friction between crack surfaces during the compression-shear fracture of closed cracks, the influence of friction coefficient between crack surfaces on stress intensity factor is dominant. Consequently, the pore pressure and swelling characteristics of rocks are not considered in the present simulation for simplicity. e elastic modulus E of 800 MPa and Poisson's ratio μ of 0.25 were used in the FEM. e maximum principal stress was used as the initial damage criterion of the material, and the value was set as 25 MPa. e parameter of damage evolution was 0.0005. e viscosity coefficient of 0.1 for damage stabilization was to facilitate convergence of the results. Table 2 compares the theoretical and numerical results, and the relative errors are also shown in the table. It is found that the relative errors of modes I and II stress intensity factors are less than 5% when the friction between closed cracks is taken into account. e theoretical calculation results of mode II stress intensity factor are in good agreement with the numerical results.

Analysis of Influencing Factors on Stress Intensity Factors
Whether closed cracks (equations (9) and (11)) or open cracks (equations (16) and (17)), in addition to stress conditions, the factors affecting stress intensity factors mainly  include crack inclination, friction coefficient of crack surface, crack length, and pore pressure in water-bearing rocks. In the following section, the sensitivity of parameters of compression-shear boundary cracks containing water under different pore pressures and lateral pressures is analyzed.

Effect of Crack Inclination on Modes I and II Stress Intensity Factors.
When rock saturates naturally, the friction coefficient between crack surfaces is assumed to be 0.8. e relationship between stress intensity factor of closed crack and crack inclination is drawn as follows. Figures 5 and 6 show the influence of crack inclination on stress intensity factors of compression-shear boundary cracks under different pore pressures. Compared with Figures 5 and 6, it can be found that the variation trend of stress intensity factor with crack inclination is basically the same under the same pore pressure, whether the crack is open or closed. When the crack inclination is constant, the stress intensity factor increases with the increase of pore pressure. According to equation (17) and Figures 5 and 6, it can be seen that the tendency of the mode I fracture of the crack becomes more obvious with the increase of the inclination angle of the crack. At the same time, the more the increase of pore pressure, the more the increase of the mode I stress intensity factor. e existence of pore pressure is also conducive to crack propagation in mode I. It shows that the contribution of water to crack propagation is mainly due to the weakening effect on the physical and mechanical properties of the rock, thus affecting the fracture toughness of the rock. Moreover, as the crack inclination increases, the mode I stress intensity factor increases gradually in a shape of "S." And the rate of change is faster at the beginning, then tends to be stable, and finally slows down. At the stage where the rate of change tends to be stable, the inclination varies between 25°and 65°, and the increase of stress intensity factor is more than 75%. Figures 5(b) and 6(b) show that the mode II stress intensity factor decreases first and then increases as the crack inclination increases. When the closed crack inclination is 65 ∘ , K b II reaches the minimum value; when the open crack inclination is 45 ∘ , K b II reaches the minimum value. Due to the friction between the crack surfaces, the closed crack has a larger absolute value of the mode II stress intensity factor when the crack inclination angle is close to horizontal. For the open crack, the absolute value of the mode II stress intensity factor first increases and then decreases as the crack inclination angle increases, and it reaches the maximum at θ � 45 ∘ , showing a parabolic form. is indicates that if the crack inclination angle is too large or too small, it is not conducive to the crack propagation for the open cracks. Namely, when the boundary crack is close to the boundary or perpendicular to the boundary, it is not easy to expand.  Figures 7(a) and 8(a), the influence of crack inclination on mode I stress intensity factor is suppressed as the lateral pressure increases. e variation range of mode I stress intensity factor of crack under the relatively higher lateral pressure is less than that under the relatively lower lateral pressure. When the inclination angle is less than 10°, it can be considered that the lateral pressure has no effect on K b Ι . Later, with the increase of angle, the inhibition of lateral pressure on K b Ι becomes more and more obvious. Since the inclination of the downward crack is generally above 60°, the compression-shear crack is less prone to mode I crack propagation when the lateral pressure increases. It can be found from Figures 7(b) and 8(b) that the effect of lateral pressure on K b II is consistent with that of pore pressure when the inclination of crack is constant for mode II stress intensity factor. at is, the higher the lateral pressure is, the   6 Advances in Materials Science and Engineering greater the stress intensity factor of mode II is. In addition, as the crack inclination angle increases, the variation range of K b II at a relatively higher lateral pressure is smaller than that at a relatively lower lateral pressure. Combined with Figures 5 and 6, it is shown that the variation range of modes I and II stress intensity factor decreases significantly as the lateral pressure increases gradually, and the change of crack inclination is no longer the main factor affecting the crack growth.

Influence of Friction Coefficient of Crack Surface on Stress
Intensity Factor. Figure 9 shows the relationship between mode II stress intensity factor and friction coefficient under the high (σ x � 25 MPa) and low (σ x � 0 MPa) lateral pressures.
e friction coefficient f between the crack surfaces mainly affects the mode II stress intensity factor K b II of the closed crack. As shown in Figure 9, K b II is proportional to the friction coefficient f. When the lateral pressure is very small, the mode II stress intensity factor is negative. With the increase of friction coefficient between crack surfaces, the absolute value of mode II stress intensity factor decreases, which indicates that the increase of friction coefficient is not conducive to crack growth. At this time, the larger the pore pressure is, the smaller the decline rate of |K b ΙΙ | is. It can be concluded that the pore pressure greatly weakens the determination of friction coefficient on mode II stress intensity factor. If the pore pressure reaches a larger value, the effect of friction coefficient on K b II can be neglected, indicating that the crack is in the open state at this time, which is consistent with the actual situation. When the lateral pressure is high,  the mode II stress intensity factor is positive. With the increase of friction coefficient, the mode II stress intensity factor and its absolute value increase. is shows that the increase of friction coefficient is beneficial to crack propagation under the high lateral pressure. e smaller the pore pressure is, the larger the growth rate of |K b ΙΙ | is and the smaller variation of curve slope is. is indicates that the small variation of pore pressure has little effect on the mode II stress intensity factor with friction coefficient at high lateral pressure. In addition, when f � 0.2, pore pressure P w changes from 0 MPa to 3 MPa, and mode II stress intensity factor increases by 72.6% and 30.9% at the high and low lateral pressures, respectively. When f � 0.8, the mode II stress intensity factor increases by 3.4% and 9.4%, respectively, under the above two lateral pressures. It is shown that the greater the friction coefficient is, the smaller the effect of the pore pressure on the mode II stress intensity factor is, whether at the high lateral pressure or the low lateral pressure. On the contrary, Figure 10 shows that, considering pore pressure is constant, the greater the friction coefficient is, the greater the influence of lateral pressure on mode II stress intensity factor is. Figures 11 and 12 show the relationships between crack length and stress intensity factor for the closed and open cracks, respectively. It can be found that the stress

Effect of Crack Inclination and Friction Coefficient on Ratio of I and II Stress Intensity
Factor. e stress intensity factor of crack is the key index to influence whether crack propagates or not, but the propagating of compression-shear crack also depends on the compression-shear coefficient λ. In this section, the relationships between the critical compression-shear coefficient λ c and the crack inclination or the friction coefficient of the crack surface are discussed.
Equations (12) and (18) show that when the critical compression-shear coefficient is used to judge whether the rock is fractured or not, the criteria for judging the fracture are symbolically related to the mode I stress intensity factor. When K I < 0, if the crack propagates, λ < λ c ; when K I > 0, if the crack propagates, λ > λ c . Figure 13 shows the relationship between the critical compression-shear coefficient and the crack inclination of closed crack. It can be seen from Figure 13 that the critical compression-shear coefficient always changes near 0 when the crack inclination is less than 40°. At the moment, the mode I stress intensity factor is negative when the crack inclination is close to the horizontal direction. According to the aforementioned criteria, if the crack propagates in case of K I < 0, then λ < λ c . In fact, the compression-shear coefficient of rock is generally larger than the calculated critical value at this time, so the crack is not easy to propagate in this state. Subsequently, the growth rate of the compression-shear coefficient is accelerated until the crack inclination is 85°, where λ c is 38. When the crack inclination is greater than 85°, the critical compression-shear coefficient λ c decreases rapidly. When the crack inclination is 90°, the critical compression-shear coefficient decreases again to nearly 0. At this stage, mode I stress intensity factor continues to increase with the increase of crack inclination, and mode I fracture gradually becomes the main fracture mode. Figure 14 shows the relationship between critical compression-shear coefficient and friction coefficient when the crack inclination is 70°. At this time, pore pressure and lateral pressure are both 0 MPa. It is shown that the critical compression-shear coefficient decreases with the increase of friction coefficient between crack surfaces. e value of mode I stress intensity factor is negative at this time. When the crack propagates, the compression-shear coefficient λ should be less than the critical compression-shear coefficient λ c , that is, λ < λ c . is shows that the increase of friction coefficient is not conducive to crack propagation at this moment.

Conclusions
e fracture behavior of boundary cracks containing water in rock is discussed in this study. e detailed calculation methods of stress intensity factors for closed and open cracks are proposed by considering the influences of pore pressure, lateral pressure, friction coefficient between closed cracks, and the swelling effect.
Combining the stress intensity factors with the fracture criterion, the role of compression-shear coefficient in judging the crack initiation was highlighted. e compression-shear coefficient is related to the properties of rock materials. It can reveal the contribution of each aforementioned factor to crack initiation. e effect of water on crack propagation is mainly due to the deterioration of the physical and mechanical properties of the rock; that is, the fracture toughness is degraded. e small range change of pore pressure has little effect on the stress intensity factor. For the closed cracks, the effect of friction coefficient on mode II stress intensity factor is inhibited by the increase of pore pressure at the relatively lower lateral pressure.
is study provides a theoretical basis for estimating the instability law of boundary cracks in rock materials by compression-shear coefficient. Further studies should be performed to apply the proposed fracture theory to rock mass or water-resisting layers.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.