Crack Initiation Criteria and Fracture Simulation for Precracked Sandstones
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.
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 . 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 . 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 . 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 . The PFC2D software has been employed to discuss the effect of the initial flaw orientation in the specimen’s failure mode under compressive loading . A fracture and crack propagation analysis system has been employed in the investigation of crack growth . The maximum tangential stress , maximum energy release rate , 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 :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 ( and , ) is shown in Figure 1, where the crack length is 2a. The stress state along the crack plane can be determined as follows:where , , and are the tangential stress, normal stress, and shear stress. The tensile stress is positive, while the compressive stress is negative.
The SIF for cracks with different widths, which causes friction and no friction along the crack plane, is different. and are SIF generated by the transverse compressive stress and normal stress, respectively. is only considered when the transverse compressive stress , and when . In addition, . 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 () for a nonclosed crack (no friction along the crack plane)  (Figure 1), but 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.
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 . The equivalent mode-I stress intensity factor (KIe), which is transformed from the mixed mode I-II fracture, is defined as follows:where the initiation angle .
When KIe equals the mode-I fracture toughness (KIC), the crack initiates, and KI is 0 for a closed crack. Thus, θ0 = 70.5° (when , , 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 :where KIII is the mode-III SIF (KIII = 0) and , , , and can be obtained as follows:where μ is the shear modulus and ν is Poisson’s ratio. The initiation angle (θ0) can be obtained by and .
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 :
θ0 is expressed as follows:where , , , , , and .
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 and can be determined as follows:
Based on equation (8), can be obtained by assuming plane strain (). This can satisfy .
The twin shear stress f can be determined by the shear stresses , , and , as follows:where the principle shear stresses can be expressed as , , and and α is the ratio of to of the rocks such that .
The shear stress is constant () on the radial shear stress line, and f on the equal radial shear stress line stress can be expressed as follows:where is given by
The crack fractures along the direction of the twin shear stress minimum value, and the fracture angle are given by
The equivalent mode-II SIF that transforms from the mixed mode I-II fracture can be expressed as follows:
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 :
For most rocks, KIC is less than KIIC . 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 Figure 2. We can obtain for the modes I () and II crack shown in Figures 2(a) and 2(c), i.e., . Thus, a mode-I fracture occurs. However, there exists for the mode-I crack without friction along the crack plane (), that is, , which occurs with the mode-II fracture.
For a mixed mode I-II crack, and are superimposed. The circumferential stress field occurs according to Figures 2(a) and 2(c). If , , a mode-I fracture occurs. Otherwise, the circumferential stress decreases according to Figures 2(b) and 2(c), when . The circumferential stress field is the same as that shown in Figure 2(c) for a mixed mode I-II closed crack, when . Thus, and the fracture is mode-I. For a mixed mode I-II crack without friction, when , the fracture mode can be expressed as follows:
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 ; 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 . The lateral pressure coefficient (k > 0) exerted a significant effect on the mode-II fracture angle when β < 40° (Figures 3(a) and 3(b)).
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)).
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 and S, 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 :where is the material parameter and and are the energy release rates of the mode I-II fracture, , , and , is the stress threshold of the crack fracture, and is the displacement of the cohesive model when fracture occurs.
The computational models of the rock samples are described in Figure 5. 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.
Table 3 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.
5.2. Analysis of Closed Crack Fracture
When the crack angle was smaller, the crack extended with a small angle (1–5 in Table 3), 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 (). 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 and . 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 and . The testing results indicate that fracture criteria for predicting the crack propagation is closely related to precracked angle.
The mode-I fracture occurred in a sample with a nonclosed crack, when the crack angle was less than 45° (Figure 4(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 if and were satisfied.
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. For and , the fracture angle of the nonclosed crack was determined by MSS.
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
The authors thank the project supported by the Funds for Creative Research Groups of China (no. 41521002) and the National Natural Science Foundation of China (no. 41672282). Huang thanks the Innovative Team of the Chengdu University of Technology.
L. Z. Wu, H. Deng, R. Q. Huang, L. M. Zhang, X. G. Guo, and Y. Zhou, “Evolution of lakes created by landslide dams and the role of dam erosion: a case study of the Jiajun landslide on the Dadu River, China,” Quaternary International, vol. 503, pp. 41–50, 2019.View at: Publisher Site | Google Scholar
G. R. Irwin, “Analysis of stresses and strains near the end of a crack traversing a plate,” Journal of Applied Mechanics-Transactions, vol. 22, pp. 361–364, 1957.View at: Google Scholar
Y. P. Li and C. H. Yang, “Influence of geometric characteristics of pre-existing cracks on mixed mode fractures under compression-shear loading,” Chinese Journal of Rock Mechanics and Engineering, vol. 25, pp. 462–466, 2006.View at: Google Scholar
Y. Fujii and Y. Ishijima, “Consideration of fracture growth from an inclined slit and inclined initial fracture at the surface of rock and mortar in compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 6, pp. 1035–1041, 2004.View at: Publisher Site | Google Scholar
C. H. Park and A. Bobet, “The initiation of slip on frictional fractures,” in Proceedings of the 41st US Rock Mechanics Symposium, pp. 06–923, Golden, CO, USA, June 2006.View at: Google Scholar
C. H. Park and A. Bobet, “Crack initiation and propagation from frictional fractures,” in Proceedings of the 1st Canada—US Rock Mechanics Symposium, pp. 557–564, Vancouver, Canada, May 2007.View at: Google Scholar
A. Manouchehrian, M. Sharifzadeh, M. F. Marji, and J. Gholamnejad, “A bonded particle model for analysis of the flaw orientation effect on crack propagation mechanism in brittle materials under compression,” Archives of Civil and Mechanical Engineering, vol. 14, no. 1, pp. 40–52, 2014.View at: Publisher Site | Google Scholar
H. Dündar and A. O. Ayhan, “Multiple and non-planar crack propagation analyses in thin structures using FCPAS,” Frattura ed Integrità Strutturale, vol. 10, no. 35, pp. 360–367, 2016.View at: Google Scholar
M. A. Hussian, E. L. Pu, and J. H. Underwood, “Strain energy release rate for a crack under combined mode I and mode II,” in Proceedings of the Fracture Analysis, pp. 2–28, ASTM STP 560, West Conshohocken, PA, USA, 1974.View at: Google Scholar
B. R. Lawn and T. R. Wilshaw, Fracture of Brittle Solids, Cambridge University Press, London, UK, 1975.
N. I. Muskhelishvili, Some Basic Problems of the Mathematical Theory of Elasticity, Noordhoff International Publishing, Leyden, The Netherlands, 1953.
G. C. Sih, Mechanics of Fracture, Noordhoff International Publishing, Leyden, The Netherlands, 1973.
ABAQUS Inc, ABUQUS/Standard User’s Manual Version 6.9, ABAQUS Inc, Waltham, MA, USA, 2009.
R. Q. Huang, L. Z. Wu, Q. He, and J. H. Li, “Stress intensity factor analysis and the stability of overhanging rock,” Rock Mechanics and Rock Engineering, vol. 50, no. 8, pp. 2135–2142, 2017.View at: Google Scholar