Abstract

The effect of water on fracture toughness in various modes of mudstones was investigated using semicircular bend (SCB) specimens exposed to three-point bendings. Natural mudstone specimens are obtained using a special coring method and are then classified into three types (A, B, and C) corresponding to three directions/configurations of bedding planes (divider, arrester, and transverse). The results show that the Type A (divider configuration) specimens possess the largest fracture toughness value for all tested modes and same soaking time, whereas the Type C (transverse configuration) specimens have the smallest one. By increasing soaking time, the fracture toughness in all three modes decreases and the fracture mechanism changes from brittle failure to ductile failure. Among them, the Type C specimens have the highest degree of degradation for each soaking time period. Regarding the fracture modes, the degradation degree of KIc is higher than that of KIIc for all three types of mudstone specimens. In addition, KIIc/KIc ratio increases when soaking time is extended. Furthermore, in the initial and short soaking time stages, the experimental KIIc/KIc results are consistent with theoretical findings from modified MTS criterion. However, after being soaked 300 minutes for three types of specimens, the test curves deviate from the theoretical curves. Analogously, the mixed-mode I/II ratio of Keff to KIc is consistent with the theoretical values in the initial stage when the degree of damage is low. With soaking time increasing, the experimental curve is gradually deviated from the theoretical curve. When soaking time reaches 300 minutes, the deviation is substantial. And the test data for the Type-A specimens are observed to provide better agreements with theoretical predictions by modified MTS criterion than those for the other two types of specimens.

1. Introduction

Mudstones, a sort of sedimentary rocks containing clay minerals, are found extensively in a variety of engineering geological environments, such as mining, hydraulic fracturing, tunneling, blasting, and geological CO2 storage. Under these natural conditions, water is a significant factor that alters the physical and mechanical properties of argillaceous swelling rocks [1, 2]. Some investigations of the effect of water on clay-bearing rocks have been conducted. For instance, Valès et al. [3] evaluated the influence of water content on anisotropic mechanical characteristics of clay-bearing shales, finding that the strength of these rocks and their cohesion and elastic modulus under uniaxial or triaxial stress states reduce with the increase of saturation degree. Pham et al. [4] and Erguler and Ulusay [5] drew out similar conclusions for mudstones. Gautam and Shakoor [6] observed that rocks containing higher amounts of clay slake more extensively and rapidly under natural climatic conditions and that the water content plays a key role in reducing the strength of clay-bearing rocks. Several experimental studies on wave velocities in clay-bearing rocks with rock-water interaction have been conducted [7].

Stress intensity factor (SIF) was defined by Irwin [8] based on linear elastic fracture mechanics (LEFM) to explain stress fields around crack tips. The critical value of SIF corresponding to crack propagation is referred to as the fracture toughness, which is a key mechanical parameter that characterizes rock resistance against crack propagation [912]. Experimental results have investigated the fracture toughness of synthetic rocks containing clay (which are regarded as isotropic materials) under the influence of water content [13, 14]. Chen et al. [15] investigated the fracture toughness of anisotropic Hualien marble; their study showed that the angle of bedding influences pure mode II fracture toughness of rock samples. Wu et al. [16] obtained different values of pure mode I fracture toughness of coals and sandy mudstones for three types of specimens defined by the bedding direction of the rocks. To date, however, little has been reported on water influence on mudstones fracture toughness with bedding. The aim of the present paper is to address this shortcoming.

Multiple experimental methods and specimens have been developed to obtain the mixed-mode I/II fracture toughness of brittle materials, including a four-point shear specimen [17, 18], semicircular specimens under three-point bending (SCB) [19], compact tension-shear specimens [20], and cracked specimens of chevron notched Brazilian disk [21, 22]. SCB is the best experimental configuration among these specimens and methods for evaluating modes I and II and mixed-mode I/II fracture toughness of rock materials owing to its simple machining process, ease of testing by simple application of three-point compressive loadings, and the minimal use of rock materials. Since Chong et al. [2325] advocated the use of SCB to evaluate different modes of fracture toughness in rock, Lim et al. [26] performed the simulation of the SIFs of SCB specimens in different geometries of specimens and showed that SCB technique could be applied to a range of mixed-mode I/II configurations. A test standard suggested by International Society for Rock Mechanics (ISRM) for mode I static fracture toughness of rock employing SCB specimen has been presented [19]. Ayatollahi et al. [27, 28] introduced the fracture parameter T-stress to calculate fracture toughness in SCB specimens using finite element method, and modified maximum tangential stress (MMTS) criterion has been used calculate mode II and I/II mixed-mode fracture toughness in terms of mode I fracture toughness. This test uses a specific coring method to acquire natural mudstone cores for machining SCB specimens. The mudstone SCB specimens are classified into three types based on orientation of the bedding plane and used to conduct experimental research on the influence of water on the pure modes I and II and mixed-mode I/II fracture toughness of mudstones. Furthermore, based on the MMTS criterion, a comparison is performed between test results and theoretical values.

2. Experimental Methods and Progress

2.1. SCB Specimens under Three-Point Bendings

Chong and Kuruppu [23] suggested using an SCB specimen exposed to three-point bending to perform all modes of fracture toughness tests on rock samples, including pure modes I and II and mix mode I/II. Figure 1 shows a semicircular disc specimen with the radius of R and a radial edge crack with the length of α previously prepared from the center. Fracture toughness tests are conducted using a vertical compressive force () along specimen axis on a three-point-bending machine. KI and KII denoted SIFs of modes I and II for SCB specimen, respectively, and are written as follows [29]:where is compressive force. Crack propagation begins when compressive force reaches peak load, Pcr. Fracture toughness KIc and KIIc values are obtained by equations (1) and (2). Also, B is specimen thickness, β is the departure angle of preexisting crack from vertical direction, α/R is the ratio of crack length to semicircular radius, and S/R is the ratio of half distance between two supporting cylindrical rollers to radius. YI and YII are geometric factors of modes I and II, respectively, which are functions of S/R, α/R, and β geometric parameters. Test methods for various fracture toughness modes are achieved by the variation of the values of these geometric parameters. For instance, at β = 0, pure mode I loading is under SCB specimen. Ayatollahi and Aliha [29] used finite element method to determine YI and YII values for pure modes I and II loading under the conditions of different geometries and determined the value of β of pure mode II loadings for various S/R and α/R values.

2.2. Materials and Sample Preparation

Mudstone samples used in this work were obtained from 150–200 m depth of Xiegou Coal Mine in Lvliang City, Shanxi Province, China. The mineral component of the rock was determined semiquantitatively using X-ray diffraction (XRD), and Figure 2 shows a typical XRD pattern. Clay content in the specimens is 34%, including kaolinite and illite. The rest is composed predominantly of silicon dioxide and potash feldspar. It is almost impossible to obtain cylinder mudstone specimens by traditional drilling cores due to limitations of the brittle bedding architecture, the likelihood of mechanical damage, and the influence of the water. Several investigators have used synthetic mudstone instead of natural mudstone for test specimen preparation [30, 31]. However, in this test, a special coring method was used to acquire natural mudstone specimens in different directions of stratification. A 1.5 mm carborundum wire saw in a numerical control machine with lubricating oil coolant was used to collect core samples. Sedimentary rocks generally have layered structures and are considered to have transversely isotropic properties. For mudstone in particular, the influence of water on clay-bearing sedimentary rock specimens increases the degree of anisotropy. The ISRM has suggested that the sedimentary rock SCB specimens whose notches are oriented along three mutually perpendicular directions defined as arrester, short transverse, and divider directions should be tested to determine mode I fracture toughness [19].

This test defined three types of the SCB specimens designated as Types A (divider), B (arrester), and C (transverse) corresponding to three bedding directions (see Figure 3). Consequently, the semicircular specimens were cut into pieces with radius 25 mm (R = 25 mm) and thickness 20 mm (B = 20 mm). On this basis, a straight notch of each semicircle was previously prepared with a 0.5 mm diamond-impregnated fine wire saw. α/R ratio value was 0.35, and the values of β were 0°, 5°, 15°, 30°, 45°, and 54°. Figure 4 shows the geometrical characteristics of the SCB specimens (Types A, B, and C).

All the samples were placed in a drying oven for 24 h at 105°C and were then cooled to room temperature. To investigate early effect on mudstone, the SCB specimens were submerged in water for different duration. For characteristics and the size of this mudstone specimen, the specimens were placed in water for 0, 50, 100, 200, and 300 minutes to impose different degrees of damage. Here, S/R ratio was 0.5 and displacement loading rate was set at 0.02 mm/min to determine static fracture toughness [32]. Geometric factors YI and YII of SCB specimens at a/R and S/R values of 0.35 and 0.5, respectively, were calculated by finite element method under pure mode II and mixed-mode I/II loadings [33]. Pure mode I fracture toughness was calculated under YI ≠ 0 and YII = 0 conditions for β = 0°. Pure mode II fracture toughness was obtained under YI = 0 and YII ≠ 0 conditions for β = 54° (Ayatollahi and Aliha [33]). Finally, mixed-mode I/II fracture toughness was obtained under  ≠ 0 and  ≠ 0 conditions for β ≠ 0° and β ≠ 54° (β = 5°, 15°, 30°, and 45° were tested here). The parameter values of YI and YII are listed in Table 1. KI and KII of the mudstone are calculated by equations (1) and (2). In this work, at least three specimens were tested for each condition. In total, approximately 300 mudstone SCB specimens were tested.

3. Results and Discussion

3.1. Pure Mode I Fracture Toughness of Mudstones

Tables 24 summarize mode I fracture toughness of three mudstone specimens with various soaking durations. For the initial specimens, the average mode I fracture toughness was 0.85 MPa·m0.5, 0.82 MPa·m0.5, and 0.73 MPa·m0.5. For each soaking time, average mode I fracture toughness of different mudstone specimens also showed differences. Load-displacement curves of three types of mudstone samples after they were soaked in water for the pure mode I loading condition are shown in Figure 5. Mudstone peak load after being submerged in water is smaller than the peak load of intact mudstone, and peak load is decreased with increasing soaking time.

The load-displacement curve is concave in initial compaction loading stage because of microcracks inside rock closing exposed to external forces. In elastic compaction stage, the load shows an approximately linear relationship with the displacement up to the peak load. Increasing submersion time widens compaction stage and load-displacement curve slope in elastic stage is sharply decreased. This change is more apparent in Type B specimens, where the direction of loading is perpendicular to bedding plane, than in the other two types. Furthermore, brittle fracture occurs in the undamaged and low-damage mudstone specimens. However, at a soaking time of 300 minutes, 200 minutes, and 200 minutes for Type A, B, and C specimens, respectively, mudstone specimens become more softened and ductile. When the soaking time is far more than 300 minutes, the three types of these mudstone samples will disintegrate. Water has a prominent effect on mudstone failure characteristics, which transforms from being brittle to being ductile. In addition, the load-displacement curves of three different specimens in mode II or mixed-mode I/II loadings variation tendencies are approximately consistent with that of the mode I loading.

3.2. The Effect of Water on the Fracture Toughness of the Mudstones

Several pure modes I and II and mixed-mode I/II fracture toughness tests were performed on mudstone for various soaking times with the three types of SCB specimens. After soaking for 50, 100, 200, and 300 minutes in water, the average gravimetric water content of the specimens was 1.8%, 3.3%, 5.9%, and 7.8%, respectively. After a short soaking time, the interaction of rock and water is not strong, and water content is lower. With increasing soaking time, the water-rock interaction becomes stronger, and fracture toughness of the mudstone degrades rapidly. Fracture toughness degradation degree is obtained bywhere Kit = different mode fracture toughness of mudstone after t minutes of soaking in the water. Deterioration trends of three specimens are equal. Nevertheless, the degree of degradation of the Type C specimens is higher than that of the other types for each time point. During contact with water, because of the presence of clay minerals such as kaolinite and illite particles that exhibit strong hydrophilicity, a water molecular layer is formed and continuously extended as water penetrates between these particles. In addition, the water molecules enter illite; consequently, the clay mineral crystal cell layers cause the spacing to increase, and the binding force is reduced between the crystals. The chemical reaction process is as follows [35]:

The first mechanism causes external expansion of the clay minerals, and the second causes internal expansion. Nonuniform stress is produced due to uneven expansion within the rock. Furthermore, the coarse grains and their intermediate cementation become softened and dissolve, increasing internal microcracks and expansion. With increasing soaking time, the fracture toughness of mudstone gradually decreases.

3.3. The Pure Mode II Fracture Toughness of Mudstones

Pure mode II fracture toughness values of three different types of initial mudstone specimens are 0.272 MPa·m0.5, 0.269 MPa·m0.5, and 0.255 MPa·m0.5, as shown in Tables 57. By increasing soaking time, mode I fracture toughness degradation degree is larger than the corresponding value for mode II. It demonstrates that the water sensitivity of pure mode II fracture toughness is less than that of mode I. When the mudstone interacts with water, the tensile (mode I) microcracks tend to appear around the crystal in the mudstone specimens [36]. It causes the mode I fracture toughness of mudstone specimens easier to reduce after soaking in water.

Along with experimental methods, many theoretical criteria such as maximum energy release rate [37], minimum strain energy density [38], and maximum tangential stress (MTS) [39] have been developed for fracture behavior prediction in mixed modes. Among the three criteria, the MTS criterion is the most prominent and widely used because it is convenient and quantifiable. The traditional MTS criterion states that fracture happens when tangential stress σθθ at a certain distance rc along the direction defined by θm is above the critical stress σθθc. The stresses near the crack tip depend on the KI and KII. The traditional MTS criterion is given by

According to conventional and modified MTS criteria developed by Smith et al. [40] considers the T-stress which is a nonsingular and constant stress term and depends on the sample shape and loading conditions. T-stress has been calculated by Ayatollahi et al. [41] through finite element method, which is expressed for SCB specimens using the dimensionless parameter [29]:

The dimensionless parameters hinges on the geometric parameters S/R, α/R, and β. Ayatollahi and Aliha [27] demonstrated in detail that the modified MTS criterion could determine mode II fracture toughness using SCB specimens. Fracture toughness ratio KIIc/KIc is obtained as follows:

According to equations (2) and (6), can be simplified as

Dimensionless parameters and YII are obtained by the finite element method [29]. In this test, is 4.678, corresponding to the geometrical parameter α/R being 0.35 and S/R being 0.5. Finally, equation (7) is rewritten aswhere rc is a constant material parameter independent from specimen geometry. The critical distance rc in rocks is defined as fracture process zone (FPZ) [42] in front of crack tip and is created in rocks because of microcrack initiation and coalescence by increasing the applied load. FPZ size in rock is estimated by maximum principal stress model suggested by Schmidt [43]:where KIc is mode I fracture toughness of mudstone and σt is rock tensile strength. In this test, three types of disc mudstone specimens, corresponding to the three types of the mudstone SCB specimens soaking for the same time, are used to calculate σt by the Brazilian test [44], as shown in Figure 6. The rock tensile strength and calculated rc values for the mudstone are summarized in Table 8.

In this test, the crack angle β of the SCB specimen is taken as 54° to calculate the pure mode II fracture toughness. Figure 7 depicts crack-growth trajectory of several typical SCB specimens with different degrees of damage under mode II loading. The Type A specimens AII0-1, AII50-1, AII100-1, AII200-1, and AII300-1 were separately immersed in water for 0, 50, 100, 200, and 300 minutes, respectively, and the other two types of specimens were exposed to the same soaking times as these Type A specimens. In the initial and low degree of damage stages, a distinct macrocrack starts at the crack tip up to axis loading location. With increasing soaking time, small cracks appear on the surface of the specimens, and the macrocrack initiation deviates gradually from the crack tip. The bedding structure of the mudstone influences the crack propagation, and the effect becomes stronger with the deepening degree of damage, particularly for the Type C specimens. With increasing soaking time, the crack propagation path of the Type C specimens overlaps with the direction of the weakened bedding. The crack spreading trends of the three types of specimens under mixed-mode I/II loadings are identical to mode II.

Based on the conventional MTS criterion, fracture toughness ratio KIIc/KIc constantly equals 0.87. Nevertheless, this test ratio is, in reality, smaller than 0.87. According to equation (9), the fracture toughness ratio KIIc/KIc is obtained in terms of using modified MTS criterion. Figure 8 shows the plot of relevant numerical results, and mode II fracture toughness is decreased as increases.

Figures 9(a)9(c), corresponding to the three types of the specimens, compares KIIc to KIc ratio achieved from test results with theoretical values obtained from modified MTS criterion. In the initial stage, when the degree of damage is low, the curve of the experimental results is in accordance with the theoretical curve. However, after soaking Types A, B, and C specimens for 300 minutes, test curves and theoretical curves clearly deviate. With increasing soaking time, the ratio of KIIc to KIc increases. In addition, the test results of the Type B specimens have less consistency with the theoretical values than the other two types of specimens. The bedding has an effect on mudstone mode II fracture toughness.

3.4. Mixed-Mode I/II Fracture Toughness of Mudstones

We also measured mixed-mode I/II fracture toughness of three mudstone specimens with different degrees of water damage. The mixed-mode equivalent intensity factor Keff and the mixity parameter Me are [38]where the KIm (KIIm) is mode I (II) stress intensity factor. Me value is 1 and 0 for pure modes I and II, respectively. According to equation (5), considering T-stress and ignoring higher order terms of O (r1/2), mixed-mode fracture toughness ratio Keff/KIc is determined by Ayatollahi and Aliha [34] aswhere rc is the fracture progress zone of the three types of SCB specimens with different soaking times (refer to Table 8). By substituting equations (1), (2), and (6) into equations (12) and (13), the mixture parameter Me and mixed-mode fracture toughness ratio Keff/KIc are achieved for SCB specimens in terms of YI, YII, and as follows:

Through finite element analyses, Ayatollahi and Aliha [34] obtained the nondimensional parameters YI, YII, and of SCB specimens under various crack inclination angles (β). Based on equation (14), the mixity parameter Me can be obtained and is shown in Table 1.

Each type of mudstone SCB specimens includes four groups of mixed I/II fracture toughness experiments corresponding to 4 crack inclination angles (5°, 15°, 30°, and 45°) under varying degrees of water damage. There are at least three effective sets of experimental data for each experimental condition. Figure 10 describes the mixed I/II fracture toughness Keff/KIc of mudstone SCB specimens with different soaking duration by the modified MTS criterion. Figure 10 shows a curve fitted to the average of the test data employing least-squares technique. In the no- and low-damage stages, the experimental ratio of Keff to KIc is in good compliance with theoretical predictions through modified MTS criterion. With increasing soaking time, the degree of damage deepens, and microcracks arise inside the mudstone specimens. In addition, the bedding has a relatively large effect on mixed-mode I/II fracture toughness. Test data curve is gradually deviated from theoretical curve. When soaking time reaches 300 minutes for the three types of specimens, these deviations become easily detectable. By comparing the ratio Keff/KIc for the three types of mudstone SCB specimens, the test data for the Type A specimens is observed to provide better compliance with theoretical predictions by modified MTS criterion than those for the other two types of specimens because the SCB specimens are simplified as a plane stress model and because the Type A specimens are regarded as an isotropic material in which the influence of bedding is disregarded to a certain extent.

4. Conclusions

In this test, the effects of water on the fracture toughness of mudstones in pure modes I and II and mix mode I/II with bedding are investigated using three types of SCB specimens. Theoretical results and experimental findings according to the modified MTS criterion are compared here. The following conclusions are drawn in this work:(1)This test used a special coring method to obtain native mudstone SCB specimens, which are classified into three types corresponding to the three bedding directions. The fracture toughness of Type A (divider configuration) specimens has the largest value for all tested modes and same soaking time and that of the Type C (transverse configuration) specimens is considerably smaller than the other two types.(2)With increasing soaking time, fracture toughness for all tested modes and for all three types of mudstone specimens decreases, and failure characteristics are changed from being brittle to being ductile. In particular, the degree of degradation of the Type C specimens is higher than that of the other types for each soaking time period.(3)The degree of degradation of pure mode I fracture toughness is slightly higher than pure mode II, proving that the water sensitivity of pure mode I fracture toughness is higher. KIIc/KIc ratio increases when the soaking time is extended. In the initial and short soaking time stages, the experimental KIIc/KIc values are consistent with theoretical findings according to modified MTS criterion. However, after being soaked 300 minutes, the test curves deviate from the theoretical curves.(4)Based on the modified MTS criterion, mixed-mode I/II ratio of Keff to KIc is consistent with the theoretical values for a low degree of water damage. With soaking time increasing, test data curve is gradually deviated from the theoretical curve. When soaking time reaches 300 minutes for the three types of specimens, the deviation is substantial. The test data characterizing the Type A specimens are in better agreement with the theoretical prediction than the other two specimen types.

Data Availability

Data supporting this research article are available from the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the support by the National Natural Science Foundation for Young Scientists of China (grant no. 51904202) and China Scholarship Council (201703780108).