# Crack Initiation Criteria and Fracture Simulation for Precracked Sandstones

Authors: R. Q. Huang; L. Z. Wu; B. Li
Journal: Advances in Materials Science and Engineering
 (2019)

---

## Abstract

The friction coefficient, tip curvature, and different-width crack state influence the stress intensity factor (SIF). The maximum circumferential tensile stress (MTS) and minimum strain energy density criterion (S) face challenges in explaining the mode-II fracture propagation of cracks. The maximum radial shear stress (MSS) and modified twin shear stress factor (ITS) criteria are proposed as the brittle mode-II fracture criteria. The experiments and numerical analysis are also performed. The results indicate that the fracture angles of the MSS and ITS were similar and different from the results of MTS and S. The equivalent stress intensity factors (ESIFs) from the mixed mode I-II are proposed to determine the fracture mode. There are different fracture models for different cracks under tensile and compressive stresses. The ratio of the tensile strength to uniaxial compressive strength influenced the fracture angle of ITS. The lateral pressure coefficient (k) had a significant effect on the mode-II fracture angle when the angle between the crack and the vertical direction is less than 40° and the lateral pressure coefficient is more than 0. Because the same fracture mode k (k > 0) can inhibit mode-I fracturing, conversely, it can also promote mode-I fracturing. Experimental results and numerical simulations of fracture propagation under uniaxial compression confirmed that the theoretical results were correct.


---

## Body

## 1. Introduction

Rock develops in complex geological environments and includes various defects or flaws. These flaws weaken the mechanical properties of the rock mass and modify the stress distributions. Internal cracks in rock mass exert an important effect on determining the fracture mode, initiation, propagation, and rock strength [1–9].The cracks in a rock mass can be divided into three basic modes: mode-I crack (normal load), mode-II crack (shear load), and mode-III crack (antiplane shear load). A tensile, shear, and tear crack fracture is defined as a mode I, II, and III fracture, respectively [10]. The mixed mode I-II is the most common type of a mixed mode fracture. After fracture mechanics was introduced into rock mechanics, the SIF was employed to describe the stress state at the crack tip [11]. When the shear stress acting on the main crack exceeds the friction stress between the cracks, the stress will concentrate at the crack tip. The cracks continue to grow and curve toward the direction of the maximum principal stress, when the stress strength factor meets or exceeds KIC. Based on the maximum shear stress theory, a method has been proposed to compute the crack initiation angle under a mixed mode I-II fracture [12]. The open-type crack geometry influences crack propagation [12, 13]. The crack initiation and growth on a rock specimen subjected to compressive stress has been investigated experimentally [1, 3, 7, 14–26]. Many numerical methods have been used to analyze the fracture and crack propagation in rocks, such as the finite element method (FEM), boundary element method (BEM), and discrete element method (DEM). By using the finite element fracture software called Franc2D, the energy release rate (G), crack propagation, fracturing time, and static tensile and normal-distributed stresses were calculated to represent the crack initiation and growth in a rock specimen [11, 27]. A realistic failure process analysis has been developed to simulate the cracked rock failure [28]. The PFC2D software has been employed to discuss the effect of the initial flaw orientation in the specimen’s failure mode under compressive loading [29]. A fracture and crack propagation analysis system has been employed in the investigation of crack growth [30]. The maximum tangential stress [31], maximum energy release rate [32], and minimum energy density criterion [33, 34] have typically been considered as the fracture initiation criteria to identify the crack growth mechanism of brittle rocks. The F-criterion and modified energy release rate criterion have also been used to investigate quasi-brittle fracture characteristics [35–37].These criteria are based on the assumption of a mode-I fracture. However, the mode-II fracture extension of compression-shear cracks has rarely been investigated. Approaches toward identifying a fracture mode that can determine suitable fracture criteria are lacking. Therefore, it is difficult to determine the fracture mode of a crack and the modes I, II, or III fracture toughness according to experimental results.The main objective of this paper was to propose a mode-II fracture criterion and conduct a numerical analysis of an open and closed brittle rock cracks. Mixed mode I-II ESIFs are proposed to determine the fracture mode. The relationships between the fracture angle and model of crack propagation with a crack angle and thickness, lateral pressure coefficient, and ratio of tensile strength to compressive strength are discussed. The consistency between the theoretical and numerical results was verified.


## 2. SIF of Mixed Mode I-II Crack

According to the linear elastic method, the stress components of the mixed mode I-II at the crack tip in the polar coordinates can be expressed as follows [38]:(1a)σr=122πrKIcosθ23−cosθ+KIIsinθ23cosθ−1,(1b)σθ=122πrcosθ2KI1+cosθ−3KIIsinθ,(1c)τrθ=122πrcosθ2KIsinθ+KII3cosθ−1,where KI is the mode-I SIF, KII is the mode-II SIF, r is the distance from the crack tip, and θ is the angle by which the surface deviates from the original crack tip direction.An infinite plate with a central crack under biaxial loading (σy∞ and σx∞, σx∞=kσy∞) is shown in Figure 1, where the crack length is 2a. The stress state along the crack plane can be determined as follows:(2a)σT=σy∞cos2β+σx∞sin2β,(2b)σN=σy∞sin2β+σx∞cos2β,(2c)τ=σy∞−σx∞sinβcosβ,where σT, σN, and τ are the tangential stress, normal stress, and shear stress. The tensile stress is positive, while the compressive stress is negative.Figure 1
External and internal stresses for mixed mode I-II crack in infinite plate.The SIF for cracks with different widths, which causes friction and no friction along the crack plane, is different.KIT and KIN are SIF generated by the transverse compressive stress and normal stress, respectively. KIT is only considered when the transverse compressive stress σT<0, and KIT=−1/2σTρ/aπa when ρ/a⟶0. In addition, KIN=σNπa. Mode-I SIF (KI) is only considered when the crack is tensile (σN > 0) for a closed crack (causing friction), and KI = 0 when σN ≤ 0. Mode-I SIF is affected by the crack tip radius of curvature (ρ) and the transverse compressive stress (σT) for a nonclosed crack (no friction along the crack plane) [39] (Figure 1), but KI≤0 for σN ≤ 0. KI exerts an inhibitory effect on the circumferential stress of mode-II SIF (KII). Table 1 lists the results of calculating SIF for closed and nonclosed cracks.Table 1
SIFs for closed and nonclosed cracks.
Crack thickness (mm)
Loading conditions
K
IN
K
IT
K
I
K
II
0.5
σ
N
>
0
σ
N
π
a
0
K
IN
+
K
IT
τ
π
a
σ
N
<
0
,
τ
>
u
σ
N
0
0
K
IN
+
K
IT
τ
−
u
σ
N
π
a
σ
N
<
0
,
τ
≤
u
σ
N
0
0
K
IN
+
K
IT
0
2.0
σ
T
<
0
σ
N
π
a
−
1
/
2
σ
T
ρ
/
a
π
a
ρ
/
a
⟶
0,
K
IN
+
K
IT
τ
π
a
σ
T
>
0
σ
N
π
a
0
K
IN
+
K
IT
τ
π
a


## 3. Fracture Criterion of Brittle Material

### 3.1. Maximum Circumferential Tensile Stress Theory (MTS)

The maximum circumferential tensile stress criterion can effectively explain the tensile fracture of brittle rock. The crack propagation direction at the crack tip can be obtained by using the maximum circumferential stress [31]. The equivalent mode-I stress intensity factor (KIe), which is transformed from the mixed mode I-II fracture, is defined as follows:(3)KIe=12cosθ02KI1+cosθ0−3KIIsinθ0,where the initiation angle θ0=2arctan1−1+8KI/KII2/4KII/KI.WhenKIe equals the mode-I fracture toughness (KIC), the crack initiates, and KI is 0 for a closed crack. Thus, θ0 = 70.5° (when KII/KI⟶∞, limKII/KII⟶∞1−1+8KII/KI2/4KII/KI=1/2, thus θ0 = 70.5°), and KIe becomes maximum for closed cracks.


### 3.2. Minimum Strain Energy Density Criterion (S)

The fracture angle of the crack is determined by the direction of the minimum strain energy density, and the strain energy density of a near crack tip element is expressed as follows [40]:(4)S=a11KI2+2a12KIKII+a22KII2+a33KIII2,where KIII is the mode-III SIF (KIII = 0) and a11, a12, a22, and a33 can be obtained as follows:(5a)a11=116πμ3−4v−cosθ1+cosθ,(5b)a12=18πμcosθ−1+2vsinθ,(5c)a22=116πμ41−v1−cosθ+1+cosθ3cosθ−1,(5d)a33=14πμ,where μ is the shear modulus and ν is Poisson’s ratio. The initiation angle (θ0) can be obtained by ∂S/∂θ=0 and ∂2S/∂θ2>0.There exist more fracture criteria, including the maximum potential energy release rate criterion and energy-momentum tensor criterion. However, these fracture criteria require a mode-I fracture. Therefore, the radial shear stress criterion and modified twin shear stress factor criterion for a mode-II fracture are proposed.


### 3.3. Maximum Radial Shear Stress Criterion (MSS)

The maximum radial shear stress at the crack tip should satisfy the following [41]:(6)∂τrθ∂θ=0,∂2τrθ∂θ2<0or∂2τrθ∂θ2>0,τrθθ=θ0max.θ
0 is expressed as follows:(7)θ0=2arctan−2+Acosα0/3−3sinα0/33k0,where A=4+42KII/KI2, B=−4KII/KI, α0=arccosT, T=−4A−3k0B/2A3, k0=2KII/KI, and α0∈0,π,−1<T<1.


### 3.4. Modified Twin Shear Stress Factor Criterion (ITS)

The twin shear stress factor criterion can be used to analyze the initiation angle of a pure mode-I fracture. However, there is significant deviation in the investigation of a mixed mode I-II fracture. Moreover, an improved twin shear stress factor criterion is proposed to predict the mode-II fracture angle.For 2D plane strain problems, the principal stressesσ1 and σ3 can be determined as follows:(8)σ1,3=12σr+σθ±12σr−σθ2+4τrθ2.Based on equation (8), σ2 can be obtained by assuming plane strain (ε2=0). This can satisfy σ1>σ2>σ3.The twin shear stressf can be determined by the shear stresses τ12, τ13, and τ23, as follows:(9a)f=τ13+ατ12=σ121+α−12ασ2+σ3,τ12>τ23,(9b)f=τ13+ατ12=12σ1+ασ2−σ32α+1,τ12<τ23,where the principle shear stresses can be expressed as τ12=0.5σ1−σ2, τ23=0.5σ2−σ3, and τ13=0.5σ1−σ3 and α is the ratio of σt to σc of the rocks such that α=σt/σc.By substituting equations (1a)–(1c) and (8) into equations (9a) and (9b), f can be expressed as follows:(10)f=122πrTsKI,KII,θ,where TsKI,KII,θ is the twin shear stress factor.The shear stress is constant (C=τrθ) on the radial shear stress line, and f on the equal radial shear stress line stress can be expressed as follows:(11)f=CFKI,KII,θ,where FKI,KII,θ is given by(12)FKI,KII,θ=TsKI,KII,θcosθ/2KIsinθ+KII3cosθ−1.The crack fractures along the direction of the twin shear stress minimum value, and the fracture angle are given by(13)∂F∂θ=0,∂2F∂θ2<0,F<0,C<0.The equivalent mode-II SIF that transforms from the mixed mode I-II fracture can be expressed as follows:(14)KIIe=12cosθ02KIsinθ0+KII3cosθ0−1.


## 3.1. Maximum Circumferential Tensile Stress Theory (MTS)

The maximum circumferential tensile stress criterion can effectively explain the tensile fracture of brittle rock. The crack propagation direction at the crack tip can be obtained by using the maximum circumferential stress [31]. The equivalent mode-I stress intensity factor (KIe), which is transformed from the mixed mode I-II fracture, is defined as follows:(3)KIe=12cosθ02KI1+cosθ0−3KIIsinθ0,where the initiation angle θ0=2arctan1−1+8KI/KII2/4KII/KI.WhenKIe equals the mode-I fracture toughness (KIC), the crack initiates, and KI is 0 for a closed crack. Thus, θ0 = 70.5° (when KII/KI⟶∞, limKII/KII⟶∞1−1+8KII/KI2/4KII/KI=1/2, thus θ0 = 70.5°), and KIe becomes maximum for closed cracks.


## 3.2. Minimum Strain Energy Density Criterion (S)

The fracture angle of the crack is determined by the direction of the minimum strain energy density, and the strain energy density of a near crack tip element is expressed as follows [40]:(4)S=a11KI2+2a12KIKII+a22KII2+a33KIII2,where KIII is the mode-III SIF (KIII = 0) and a11, a12, a22, and a33 can be obtained as follows:(5a)a11=116πμ3−4v−cosθ1+cosθ,(5b)a12=18πμcosθ−1+2vsinθ,(5c)a22=116πμ41−v1−cosθ+1+cosθ3cosθ−1,(5d)a33=14πμ,where μ is the shear modulus and ν is Poisson’s ratio. The initiation angle (θ0) can be obtained by ∂S/∂θ=0 and ∂2S/∂θ2>0.There exist more fracture criteria, including the maximum potential energy release rate criterion and energy-momentum tensor criterion. However, these fracture criteria require a mode-I fracture. Therefore, the radial shear stress criterion and modified twin shear stress factor criterion for a mode-II fracture are proposed.


## 3.3. Maximum Radial Shear Stress Criterion (MSS)

The maximum radial shear stress at the crack tip should satisfy the following [41]:(6)∂τrθ∂θ=0,∂2τrθ∂θ2<0or∂2τrθ∂θ2>0,τrθθ=θ0max.θ
0 is expressed as follows:(7)θ0=2arctan−2+Acosα0/3−3sinα0/33k0,where A=4+42KII/KI2, B=−4KII/KI, α0=arccosT, T=−4A−3k0B/2A3, k0=2KII/KI, and α0∈0,π,−1<T<1.


## 3.4. Modified Twin Shear Stress Factor Criterion (ITS)

The twin shear stress factor criterion can be used to analyze the initiation angle of a pure mode-I fracture. However, there is significant deviation in the investigation of a mixed mode I-II fracture. Moreover, an improved twin shear stress factor criterion is proposed to predict the mode-II fracture angle.For 2D plane strain problems, the principal stressesσ1 and σ3 can be determined as follows:(8)σ1,3=12σr+σθ±12σr−σθ2+4τrθ2.Based on equation (8), σ2 can be obtained by assuming plane strain (ε2=0). This can satisfy σ1>σ2>σ3.The twin shear stressf can be determined by the shear stresses τ12, τ13, and τ23, as follows:(9a)f=τ13+ατ12=σ121+α−12ασ2+σ3,τ12>τ23,(9b)f=τ13+ατ12=12σ1+ασ2−σ32α+1,τ12<τ23,where the principle shear stresses can be expressed as τ12=0.5σ1−σ2, τ23=0.5σ2−σ3, and τ13=0.5σ1−σ3 and α is the ratio of σt to σc of the rocks such that α=σt/σc.By substituting equations (1a)–(1c) and (8) into equations (9a) and (9b), f can be expressed as follows:(10)f=122πrTsKI,KII,θ,where TsKI,KII,θ is the twin shear stress factor.The shear stress is constant (C=τrθ) on the radial shear stress line, and f on the equal radial shear stress line stress can be expressed as follows:(11)f=CFKI,KII,θ,where FKI,KII,θ is given by(12)FKI,KII,θ=TsKI,KII,θcosθ/2KIsinθ+KII3cosθ−1.The crack fractures along the direction of the twin shear stress minimum value, and the fracture angle are given by(13)∂F∂θ=0,∂2F∂θ2<0,F<0,C<0.The equivalent mode-II SIF that transforms from the mixed mode I-II fracture can be expressed as follows:(14)KIIe=12cosθ02KIsinθ0+KII3cosθ0−1.


## 4. Fracture Mode and Fracture Criterion

### 4.1. Identification of Fracture Mode

The dimensionless stress field at the crack tip can be expressed as follows [42]:(15a)fθ=σθ2raσy∞+σx∞2,(15b)frθ=τrθ2raσy∞+σx∞2.For most rocks,KIC is less than KIIC [43]. The fracture mode can be determined according to the relationship between the stress field and fracture toughness.The dimensionless stress field of the mixed mode I-II crack is shown in Figure2. We can obtain fθmax>frθmax for the modes I (σN>0) and II crack shown in Figures 2(a) and 2(c), i.e., KIe/KIIe>KIC/KIIC. Thus, a mode-I fracture occurs. However, there exists fθmax<frθmax for the mode-I crack without friction along the crack plane (σN<0), that is, KIe/KIIe<0<KIC/KIIC, which occurs with the mode-II fracture.Figure 2
Dimensionless stress field of mode I-II crack: (a) mode-I crack with friction,σN>0; (b) mode-I crack without friction,σN<0; (c) mode-II crack.
(a)
(b)
(c)For a mixed mode I-II crack,KI and KII are superimposed. The circumferential stress field occurs according to Figures 2(a) and 2(c). If σN>0, KIe/KIIe>KIC/KIIC, a mode-I fracture occurs. Otherwise, the circumferential stress decreases according to Figures 2(b) and 2(c), when σN<0. The circumferential stress field is the same as that shown in Figure 2(c) for a mixed mode I-II closed crack, when σN<0. Thus, KIe/KIIe>KIC/KIIC and the fracture is mode-I. For a mixed mode I-II crack without friction, when σN<0, the fracture mode can be expressed as follows:(16a)KIe<0,for mode−II−fracture,(16b)KIeKIIe>1,for mode−I−fracture,(16c)0<KIeKIIe<1,KIeKIIe>KICKIIC,for mode−I fracture,(16d)0<KIeKIIe<1,KIeKIIe<KICKIIC,for mode−II fracture.


### 4.2. Example: Mixed Mode I-II Nonclosed Crack under Compression

By assuming that the crack aperture was 2 mm, the crack length was 2a = 10 mm, the curvature radius was ρ = 1 mm at the crack tip, and the crack was nonclosed during loading. The relationship between the crack angles and the fracture angles is shown in Figure 3, where the mixed mode I-II nonclosed crack was analyzed according to MTS, S, MSS, and ITS. The fracture angles based on MTS and S were similar. The fracture angles based on S were influenced by v; however, the results were quite different when β < 40° and k < 0. Otherwise, the fracture angles between MSS and ITS were similar and the fracture angles from ITS were influenced by α=σt/σc. The lateral pressure coefficient (k > 0) exerted a significant effect on the mode-II fracture angle when β < 40° (Figures 3(a) and 3(b)).Figure 3
Relationship between crack angle and fracture angle for mixed mode I-II nonclosed crack: (a)σy∞<0, k=0; (b) σy∞<0, k=0.4; (c) σy∞<0, k=−0.4.
(a)
(b)
(c)The mode I-II fracture regions for the mixed mode I-II nonclosed cracks (equations (16a)–(16d)) are shown in Figure 4. k exerted a significant effect on the fracture mode. The positive k (k > 0) value inhibited the mode-I fracture (Figures 4(a) and 4(b)), while the negative k value promoted the mode-I fracture (Figures 4(a) and 4(c)).Figure 4
Mode I-II fracture regions for mixed mode I-II nonclosed cracks: (a)σy∞<0, k=0; (b) σy∞<0, k=0.4; (c) σy∞<0, k=−0.4.
(a)
(b)
(c)The criteria selected to analyze the fracture angles of the crack were based on the fracture mode. In the mode-I fracture, the fracture criteria, such as MTS andS, could forecast the fracturing angle. Moreover, the MSS and ITS criteria are suggested for the mode-II fracture.


## 4.1. Identification of Fracture Mode

The dimensionless stress field at the crack tip can be expressed as follows [42]:(15a)fθ=σθ2raσy∞+σx∞2,(15b)frθ=τrθ2raσy∞+σx∞2.For most rocks,KIC is less than KIIC [43]. The fracture mode can be determined according to the relationship between the stress field and fracture toughness.The dimensionless stress field of the mixed mode I-II crack is shown in Figure2. We can obtain fθmax>frθmax for the modes I (σN>0) and II crack shown in Figures 2(a) and 2(c), i.e., KIe/KIIe>KIC/KIIC. Thus, a mode-I fracture occurs. However, there exists fθmax<frθmax for the mode-I crack without friction along the crack plane (σN<0), that is, KIe/KIIe<0<KIC/KIIC, which occurs with the mode-II fracture.Figure 2
Dimensionless stress field of mode I-II crack: (a) mode-I crack with friction,σN>0; (b) mode-I crack without friction,σN<0; (c) mode-II crack.
(a)
(b)
(c)For a mixed mode I-II crack,KI and KII are superimposed. The circumferential stress field occurs according to Figures 2(a) and 2(c). If σN>0, KIe/KIIe>KIC/KIIC, a mode-I fracture occurs. Otherwise, the circumferential stress decreases according to Figures 2(b) and 2(c), when σN<0. The circumferential stress field is the same as that shown in Figure 2(c) for a mixed mode I-II closed crack, when σN<0. Thus, KIe/KIIe>KIC/KIIC and the fracture is mode-I. For a mixed mode I-II crack without friction, when σN<0, the fracture mode can be expressed as follows:(16a)KIe<0,for mode−II−fracture,(16b)KIeKIIe>1,for mode−I−fracture,(16c)0<KIeKIIe<1,KIeKIIe>KICKIIC,for mode−I fracture,(16d)0<KIeKIIe<1,KIeKIIe<KICKIIC,for mode−II fracture.


## 4.2. Example: Mixed Mode I-II Nonclosed Crack under Compression

By assuming that the crack aperture was 2 mm, the crack length was 2a = 10 mm, the curvature radius was ρ = 1 mm at the crack tip, and the crack was nonclosed during loading. The relationship between the crack angles and the fracture angles is shown in Figure 3, where the mixed mode I-II nonclosed crack was analyzed according to MTS, S, MSS, and ITS. The fracture angles based on MTS and S were similar. The fracture angles based on S were influenced by v; however, the results were quite different when β < 40° and k < 0. Otherwise, the fracture angles between MSS and ITS were similar and the fracture angles from ITS were influenced by α=σt/σc. The lateral pressure coefficient (k > 0) exerted a significant effect on the mode-II fracture angle when β < 40° (Figures 3(a) and 3(b)).Figure 3
Relationship between crack angle and fracture angle for mixed mode I-II nonclosed crack: (a)σy∞<0, k=0; (b) σy∞<0, k=0.4; (c) σy∞<0, k=−0.4.
(a)
(b)
(c)The mode I-II fracture regions for the mixed mode I-II nonclosed cracks (equations (16a)–(16d)) are shown in Figure 4. k exerted a significant effect on the fracture mode. The positive k (k > 0) value inhibited the mode-I fracture (Figures 4(a) and 4(b)), while the negative k value promoted the mode-I fracture (Figures 4(a) and 4(c)).Figure 4
Mode I-II fracture regions for mixed mode I-II nonclosed cracks: (a)σy∞<0, k=0; (b) σy∞<0, k=0.4; (c) σy∞<0, k=−0.4.
(a)
(b)
(c)The criteria selected to analyze the fracture angles of the crack were based on the fracture mode. In the mode-I fracture, the fracture criteria, such as MTS andS, could forecast the fracturing angle. Moreover, the MSS and ITS criteria are suggested for the mode-II fracture.


## 5. Numerical Analysis and Experimental Results of Single-Crack Sample Failure

### 5.1. Numerical Model

The SIF of a single-crack sandstone was investigated by using a finite element software ABAQUS. Using the commercial finite element software, the extended FEM and a cohesive model were employed to simulate the crack propagation for a single-crack rock without considering the progressive process [44, 45]. The Benzeggagh–Kenane (B-K) model can be expressed as follows [46]:(17)GTC=GIC+GIIC−GICGsGTη,where η is the material parameter and GIC and GIIC are the energy release rates of the mode I-II fracture, Gs=GII+GIII, GT=GI+Gs, and GTC=0.5σ0uf, σ0 is the stress threshold of the crack fracture, and uf is the displacement of the cohesive model when fracture occurs.The computational models of the rock samples are described in Figure5. The mechanical simulation analysis parameters are listed in Table 2. The SIF of a single crack was obtained numerically based on the elasticity theory, by applying a compressive stress of 10 MPa and k = 0. A collapse element was employed to simulate the singularity of the crack tip.Figure 5
Calculation model diagram of single-crack sample (r = 0): (a) closed crack: crack angle β = 15°, crack thickness d = 0.5 mm; (b) nonclosed crack: crack angle β = 45°, crack thickness d = 2 mm.
(a)
(b)Table 2
Numerical parameters for analysis.
Material
Elastic modulusE (GPa)
Friction coefficientu
Poisson’s ratiov
G
IC (N/mm)
G
IIC (N/mm)
η
Maximum nominal stress of type I failureNmax (MPa)
Maximum nominal stress of type II failureSmax (MPa)
Maximum nominal stress of type III failureTmax (MPa)
Brittle material
1.354
0.12
0.21
0.62
0.76
1.0
6.63
8.09
8.09Table3 depicts the changes in the maximum circumferential stress and the radial shear stress at the crack tip. Figures 6–8 show the computational results of a single-crack brittle rock fracture with different precracked angles.Table 3
Numerical results of stress near the crack tip.
Condition
Crack thicknessd (mm)
Crack angleβ (°)
Maximum circumferential stress at crack tip (MPa)
Direction of the maximum circumferential stress (°)
Maximum radial shear stress at crack tip (MPa)
Direction of the maximum radial shear stress (°)
1
0.5
15
10.258
13
8.206
0
2
30
13.057
47
10.435
0
3
45
14.984
63
11.981
0
4
60
12.731
65
10.193
0
5
75
11.726
73
9.437
0
6
2.0
15
13.127
34
8.533
−10
7
30
15.424
45
14.677
2
8
45
12.387
81
17.749
7
9
60
7.742
101
17.066
12
10
75
1.548
125
14.678
41
The tip radius of curvature is 0.25 mm for conditions 1–5 and 1.0 mm for conditions 6–10.Figure 6
Numerical simulation of fracture propagation withβ = 15°: (a) 0.5 mm thick crack (ρ = 0.25 mm); (b) 2.0 mm thick crack (ρ = 1.0 mm).
(a)
(b)Figure 7
Numerical simulation results of single-crack fracture propagation withβ = 45°: (a) 0.5 mm thick crack (ρ = 0.25 mm); (b) 2.0 mm thick crack (ρ = 1.0 mm).
(a)
(b)Figure 8
Numerical simulation results of single-crack fracture propagation withβ = 75°: (a) 0.5 mm thick crack (ρ = 0.25 mm); (b) 2.0 mm thick crack (ρ = 1.0 mm).
(a)
(b)


### 5.2. Analysis of Closed Crack Fracture

When the crack angle was smaller, the crack extended with a small angle (1–5 in Table3), rather than with an angle of 70.5°, from the maximum circumferential stress theory. This was caused by the strong influence of material expansion, which was in turn caused by compression. Therefore, equation (1b) cannot be used to calculate the circumferential stress with a small crack angle. However, the material expansion had a smaller effect on the stress field at the crack tip, when the crack angle was greater than 15°. By using the maximum circumferential stress criterion, the fracture initiation angles in the experiments were approximately 70.5°. Therefore, equation (1b) could be used to calculate the circumferential stress in a rock with crack angles greater than 15°.The results revealed that the material expansion exerted a considerable effect on the stress at the crack tips, when the crack angle was small and the initiation angles gradually approached 70.5°, according to the maximum circumferential stress theory. Therefore, the maximum circumferential stress theory can be used to estimate the fracture criterion and calculate the compression-induced initiation angles for the closed cracks, if the crack angle is greater than 15°.


### 5.3. Fracture Analysis of Nonclosed Crack

The mode-I SIF of the nonclosed crack was negative when the crack was under compression (σN<0). According to equation (1), the mode-I circumferential compressive stress with KI < 0 at the crack tip will restrain the circumferential tensile stress caused by the mode-II SIF. The maximum circumferential tensile stress at the crack tip was less than the radial shearing stress under certain conditions. Thus, the mode-II fracture in a nonclosed crack occurred only if 0<KIe/KIIe<1 and KIe/KIIe<KIC/KIIC. Therefore, the mode domain could be divided into two regions, namely, the mode-I and mode-II fracture regions (Figure 4). Figure 9 represents the crack propagation between numerical and experimental results where crack thickness is 2 mm. The tests were carried out on an MTS815 test system (MTS Systems Corporation, Eden Prairie, MN, USA). The sample was from rock-like material with 2 mm thick crack with a width of 50 mm. The testing and numerical results show that the initiation angle of the crack and fracture propagation are similar. Figure 10 compares the initiation angles between the theoretical analysis and numerical results when σy∞<0 and k=0. The testing results indicate that fracture criteria for predicting the crack propagation is closely related to precracked angle.Figure 9
Numerical and experimental results of fracture propagation with 2 mm thick crack andβ = 75°: (a) numerical analysis; (b) experimental results.
(a)
(b)Figure 10
Comparison of initiation angles between theoretical and experimental results.The mode-I fracture occurred in a sample with a nonclosed crack, when the crack angle was less than 45° (Figure4(a)). The crack rupture angles were small, which does not agree with the initiation angles (70.5°) obtained by the MTS and listed in Table 3. For crack angles greater than 30°, there existed considerable differences between the fracture angles of the closed and nonclosed cracks (Figure 8), and the rupture angles of the open cracks were similar to the results obtained by the MSS and ITS (Figure 10). Therefore, the identification of the fracture mode and proposed mode-II fracture criteria is correct.The analysis revealed that the maximum circumferential stress criterion faces challenges in describing the fracture growth of open cracks subject to compression. The initiation angles could be determined by the radial shear stress criterion if0<KIe/KIIe<1 and KIe/KIIe<KIC/KIIC were satisfied.


## 5.1. Numerical Model

The SIF of a single-crack sandstone was investigated by using a finite element software ABAQUS. Using the commercial finite element software, the extended FEM and a cohesive model were employed to simulate the crack propagation for a single-crack rock without considering the progressive process [44, 45]. The Benzeggagh–Kenane (B-K) model can be expressed as follows [46]:(17)GTC=GIC+GIIC−GICGsGTη,where η is the material parameter and GIC and GIIC are the energy release rates of the mode I-II fracture, Gs=GII+GIII, GT=GI+Gs, and GTC=0.5σ0uf, σ0 is the stress threshold of the crack fracture, and uf is the displacement of the cohesive model when fracture occurs.The computational models of the rock samples are described in Figure5. The mechanical simulation analysis parameters are listed in Table 2. The SIF of a single crack was obtained numerically based on the elasticity theory, by applying a compressive stress of 10 MPa and k = 0. A collapse element was employed to simulate the singularity of the crack tip.Figure 5
Calculation model diagram of single-crack sample (r = 0): (a) closed crack: crack angle β = 15°, crack thickness d = 0.5 mm; (b) nonclosed crack: crack angle β = 45°, crack thickness d = 2 mm.
(a)
(b)Table 2
Numerical parameters for analysis.
Material
Elastic modulusE (GPa)
Friction coefficientu
Poisson’s ratiov
G
IC (N/mm)
G
IIC (N/mm)
η
Maximum nominal stress of type I failureNmax (MPa)
Maximum nominal stress of type II failureSmax (MPa)
Maximum nominal stress of type III failureTmax (MPa)
Brittle material
1.354
0.12
0.21
0.62
0.76
1.0
6.63
8.09
8.09Table3 depicts the changes in the maximum circumferential stress and the radial shear stress at the crack tip. Figures 6–8 show the computational results of a single-crack brittle rock fracture with different precracked angles.Table 3
Numerical results of stress near the crack tip.
Condition
Crack thicknessd (mm)
Crack angleβ (°)
Maximum circumferential stress at crack tip (MPa)
Direction of the maximum circumferential stress (°)
Maximum radial shear stress at crack tip (MPa)
Direction of the maximum radial shear stress (°)
1
0.5
15
10.258
13
8.206
0
2
30
13.057
47
10.435
0
3
45
14.984
63
11.981
0
4
60
12.731
65
10.193
0
5
75
11.726
73
9.437
0
6
2.0
15
13.127
34
8.533
−10
7
30
15.424
45
14.677
2
8
45
12.387
81
17.749
7
9
60
7.742
101
17.066
12
10
75
1.548
125
14.678
41
The tip radius of curvature is 0.25 mm for conditions 1–5 and 1.0 mm for conditions 6–10.Figure 6
Numerical simulation of fracture propagation withβ = 15°: (a) 0.5 mm thick crack (ρ = 0.25 mm); (b) 2.0 mm thick crack (ρ = 1.0 mm).
(a)
(b)Figure 7
Numerical simulation results of single-crack fracture propagation withβ = 45°: (a) 0.5 mm thick crack (ρ = 0.25 mm); (b) 2.0 mm thick crack (ρ = 1.0 mm).
(a)
(b)Figure 8
Numerical simulation results of single-crack fracture propagation withβ = 75°: (a) 0.5 mm thick crack (ρ = 0.25 mm); (b) 2.0 mm thick crack (ρ = 1.0 mm).
(a)
(b)


## 5.2. Analysis of Closed Crack Fracture

When the crack angle was smaller, the crack extended with a small angle (1–5 in Table3), rather than with an angle of 70.5°, from the maximum circumferential stress theory. This was caused by the strong influence of material expansion, which was in turn caused by compression. Therefore, equation (1b) cannot be used to calculate the circumferential stress with a small crack angle. However, the material expansion had a smaller effect on the stress field at the crack tip, when the crack angle was greater than 15°. By using the maximum circumferential stress criterion, the fracture initiation angles in the experiments were approximately 70.5°. Therefore, equation (1b) could be used to calculate the circumferential stress in a rock with crack angles greater than 15°.The results revealed that the material expansion exerted a considerable effect on the stress at the crack tips, when the crack angle was small and the initiation angles gradually approached 70.5°, according to the maximum circumferential stress theory. Therefore, the maximum circumferential stress theory can be used to estimate the fracture criterion and calculate the compression-induced initiation angles for the closed cracks, if the crack angle is greater than 15°.


## 5.3. Fracture Analysis of Nonclosed Crack

The mode-I SIF of the nonclosed crack was negative when the crack was under compression (σN<0). According to equation (1), the mode-I circumferential compressive stress with KI < 0 at the crack tip will restrain the circumferential tensile stress caused by the mode-II SIF. The maximum circumferential tensile stress at the crack tip was less than the radial shearing stress under certain conditions. Thus, the mode-II fracture in a nonclosed crack occurred only if 0<KIe/KIIe<1 and KIe/KIIe<KIC/KIIC. Therefore, the mode domain could be divided into two regions, namely, the mode-I and mode-II fracture regions (Figure 4). Figure 9 represents the crack propagation between numerical and experimental results where crack thickness is 2 mm. The tests were carried out on an MTS815 test system (MTS Systems Corporation, Eden Prairie, MN, USA). The sample was from rock-like material with 2 mm thick crack with a width of 50 mm. The testing and numerical results show that the initiation angle of the crack and fracture propagation are similar. Figure 10 compares the initiation angles between the theoretical analysis and numerical results when σy∞<0 and k=0. The testing results indicate that fracture criteria for predicting the crack propagation is closely related to precracked angle.Figure 9
Numerical and experimental results of fracture propagation with 2 mm thick crack andβ = 75°: (a) numerical analysis; (b) experimental results.
(a)
(b)Figure 10
Comparison of initiation angles between theoretical and experimental results.The mode-I fracture occurred in a sample with a nonclosed crack, when the crack angle was less than 45° (Figure4(a)). The crack rupture angles were small, which does not agree with the initiation angles (70.5°) obtained by the MTS and listed in Table 3. For crack angles greater than 30°, there existed considerable differences between the fracture angles of the closed and nonclosed cracks (Figure 8), and the rupture angles of the open cracks were similar to the results obtained by the MSS and ITS (Figure 10). Therefore, the identification of the fracture mode and proposed mode-II fracture criteria is correct.The analysis revealed that the maximum circumferential stress criterion faces challenges in describing the fracture growth of open cracks subject to compression. The initiation angles could be determined by the radial shear stress criterion if0<KIe/KIIe<1 and KIe/KIIe<KIC/KIIC were satisfied.


## 6. Conclusions

The mode-II fracture criteria and a method of fracture mode identification were proposed. We conducted a uniaxial compression numerical analysis on brittle material containing closed and open cracks at various angles and thicknesses. The numerical and theoretical results were analyzed and the following conclusions were drawn:(1)
The fracture angles, based on the MSS and ITS, were close. The ITS results were influenced by the ratio of the tensile strength to the uniaxial compressive strength of the rocks.k exerted an important effect on the mode-II fracture angle. A positive k value inhibited the mode-I fracture, while a negative k value promoted the mode-I fracture.(2)
Friction caused by crack thickness has an important role in the fracture initiation angle. The initiation direction of the closed crack at the crack tip approximates the direction of the maximum circumferential stress. When the angle between the preloaded crack and the loading direction was large, the crack initiation behavior under loading conformed to the results of the maximum circumferential stress criterion.(3)
MTS faces challenges in interpreting the fracture propagation of open cracks subject to compressional stress. The circumferential compressive stress caused by the mode-I SIF withheld the circumferential tensile stress from the mode-II SIF. For0<KIe/KIIe<1 and KIe/KIIe<KIC/KIIC, the fracture angle of the nonclosed crack was determined by MSS.



---

*Source: 9359410-2019-08-04.xml*