Abstract

To study the tensile and fracture properties of the specimen under the quasistatic loading, the Brazilian disc splitting method and the notched semicircular bend (NSCB) method were used to test the tensile properties of coal specimens, and the fracture properties of NSCB specimens with different notch depths were tested and analyzed. The applicability of plane strain fracture toughness KIC and J-integral fracture toughness in evaluating the fracture properties of coal specimens was discussed. The influence of notch depth on the fracture toughness measurement of the NSCB specimen was studied. Combined with the surface strain monitoring of specimens during loading and the industrial CT scanning image of damaged specimens, the deformation characteristics of coal specimen under loads and the distribution law of crack after failure were analyzed. The results show that the NSCB test is suitable for measuring the tensile strength of a coal specimen; when the dimensionless notch depth is β = 0.28, the dispersion of plane strain fracture toughness KIC of the NSCB specimen is the smallest. Besides, the plane strain fracture toughness of coal is obviously affected by the notch depth and dimensionless stress intensity factor. The J-integral fracture toughness can be used to effectively evaluate the fracture performance of specimens.

1. Introduction

The tensile and fracture properties of rock are of great significance for rock mass stability control [1], rock hydraulic/gas fracturing [2], and roadway roof support design [3, 4]. Many scholars have carried out extensive research on the mechanism of tunnel deformation instability and rock burst. Zhang et al. [5] established an energy integral model for the rockburst-inducing area and a friction work calculation model for the plastic area. If the remaining energy after the coal seam is broken in the rockburst-inducing area is greater than the friction work required for the coal to burst out, then a rockburst accident will occur. Zhang et al. [6] also propose a robust classifier ensemble to predict squeezing conditions in rock tunnels. Seven individual machine-learning classifiers were aggregated using weighted voting methods to establish the classifier ensemble. In the research field of seepage and failure of rock mass, Ma et al. [7, 8] quantitatively analyzed the characteristics of water-sediment flow in rock fractures by computational fluid dynamics (CFD). Based on the two-phase flow theory, a resistance model of water-sediment flow in fractures was established and verified by a laboratory-scale test. Among these studies, the testing of relevant mechanical properties of coal rocks is critical. The accurate and effective parameter measurement of tensile and fracture properties of coal rock has always been a research hotspot in rock mechanics and mining engineering. At present, the determination methods for the tensile strength of rock materials in the laboratory mainly include the direct tensile test, the Brazilian disc splitting test, beam bending test, and semicircular bend test [9]. The direct tensile test is the most direct method to obtain the tensile strength of rock materials. However, the direct tensile test requires special fixtures, and the specimen processing is difficult [9]. In the Brazilian disc splitting test, the internal stress distribution of the specimen is complex, and the compressive stress on the specimen axis is 3 times the tensile stress [10], leading to the tendency of fracture failure for specimens. The rectangular beam three-point bending test [11] has good effect in the tensile strength test, but it is difficult to process raw coal into a rectangular beam specimen, and the specimen is easily damaged by the notch and grinding. In recent years, the semicircular bend test has been widely used to test the tensile properties of rock or asphalt due to its simple operation and specimen processing. However, there are few reports about the semicircular bend test of specimens.

To measure the fracture toughness of mode I rock materials, the International Society for Rock Mechanics (ISRM) recommends four types of specimens: short rod (SR), chevron bend (CB), cracked chevron notched Brazilian disc (CCNBD), and notched semicircular bend (NSCB) specimens [12]. Fowell et al. [1315] compared SR, CB, and CCNBD specimens. It was found that there were many unfavorable factors in SR and CB tests, such as the lower load required for the specimen failure, larger core, complex loading system, and complicated preparation of short rod specimens; the test device and fixture required for CCNBD specimens were relatively simple, and CCNBD specimens were easier to prepare and can be used for the toughness test of various fracture modes. Kuruppu et al. [12, 16] believed that the NSCB test has the advantages of fewer specimens and simple loading device in testing mode I fracture toughness of rock materials under quasistatic loading. However, the analysis of crack propagation characteristics and the influence of notch depth in the fracture toughness test of rock materials using the NSCB specimens is still limited.

Besides, the NSCB testing method has been successfully applied to the fracture performance test of asphalt mixture [17, 18], mineral mixture [19], marble [20, 21], concrete [22], limestone [23], and other materials. However, there are relatively few experimental studies on the fracture toughness of coal specimens under the quasistatic condition. Coal is a complex heterogeneous anisotropic material, and the complex spatial distribution of its components and bedding increases the dispersion of NSCB test results for coal. Therefore, it is necessary to perform semicircular bend tests and NSCB tests to study its tensile and fracture properties.

In this paper, the quasistatic BD tests and SCB tests of specimens were carried out to compare the measured tensile strength, and then, the applicability of the NSCB method in measuring the tensile strength of specimens was analyzed. The fracture toughness of NSCB specimens with different notch depths was measured, and the influence of notch depth on the fracture toughness test results was studied. Combined with the surface strain monitoring of specimens under loads and CT scanning images of damaged specimens, the deformation and crack distribution characteristics of specimens under loads were analyzed.

2. Specimen Preparation and Test Method

2.1. Specimen Collection and Processing

The specimens, belonging to bituminous coal, were taken from the coal seam 11 of Xinzhouyao mine to the northeast of Datong coalfield, Shanxi Province. The moisture content of the specimen was 4.13%, ash content was 2.04%, volatile content was 25.81%, fixed carbon content was 69.17%, density was 1.301 g/cm3, elastic modulus was 5.41 GPa, and Poisson’s ratio was 0.24. Table 1 shows the maceral composition of Datong coal. According to the Brazilian disc specimen preparation method recommended by ISRM [24], 15 disc specimens (4 of which were BD specimens for testing) of Φ 50 mm × 25 mm were drilled parallel to the coal and rock bedding plane and processed. Within the angle of 40–50° between the notch surface and the bedding surface of the disk specimen, 10 disc specimens were cut in half to form 20 semicircular SCB specimens. Then, 12 specimens were taken and slit perpendicular to the bottom surface at the center of the bottom; then, 3 groups of NSCB specimens with the notch width of 1 mm and the notch depth a of 4 mm, 7 mm, and 10 mm, respectively, were taken. There were 4 specimens in each group. Table 2 shows the parameters of BD, SCB, and NSCB specimens.

ISRM suggests that when the NSCB method is used to test the mode I quasistatic fracture toughness of rock materials, the support point spacing S and specimen diameter D should be 0.5 ≤ S/D ≤ 0.8 [12]. In this paper, S/D = 0.8 was selected for testing, and the geometric dimension of NSCB specimens is shown in Figure 1. In the SCB tests of specimen, the influence of cutting loss of the specimen on tensile strength was considered, and the spacing between support points and specimen diameter S/D = 0.8 was selected referring to previous studies [25].

2.2. Test Equipment and Process

The tensile and fracture tests of semicircular specimens were carried out by using the Shimazu AGS-H5KN precision electronic universal test loading device [11]. The loading was controlled by displacement, and the loading rate was 0.2 mm/min. During the loading, the load and displacement data system were automatically recorded, and the sampling interval was 0.1 s. The SCB and NSCB specimens were placed on the top of the two supports, respectively, and the support spacing was 4 cm. The damaged specimens were scanned by using the ACTIS5.0-225X-ray Industrial CT testing system to observe the internal crack distribution characteristics.

In the fracture toughness test of BD, SCB, and NSCB specimens, strain gauges were set at different positions on the surface of specimens for deformation monitoring. Figure 2 shows the specific position and number of strain gauges. In BD specimens, 1#–4# strain gauges were arranged at equal intervals (6 mm) from the center of the specimen to the lower boundary; in SCB specimens, only one group of strain gauges were set radially at the center of the bottom surface to record the tensile strain changes; in NSCB specimens, two groups of strain gauges were arranged on the line connecting the upper loading point and the notch, the distance between the center of 1# strain gauge and the loading point was 2 mm, and the distance between the center of 2# strain gauge and the end of the notch was 1.5 mm.

3. Tensile Properties of Specimens

To compare the differences of tensile strength obtained by the BD test and SCB test, 8 specimens in two groups were carried out for BD and SCB tests. The loading-displacement data were recorded during the test, and the strain-time curve was obtained by using strain gauges. After the test, specimens were scanned by CT to analyze the internal fracture characteristics of the specimen.

3.1. BD Test Results

Table 3 shows the peak load, peak displacement, and peak time of four BD specimens. The tensile strength of specimens can be obtained from the following equation:where σt is the tensile strength, MPa; Pmax is the ultimate load, N; and D and t are the diameter and thickness of specimens, mm. The tensile strength and mean value of specimens can be calculated by the abovementioned equation.

Figure 3 shows the loading-displacement curves of four BD specimens. It can be found that the peak displacement of four BD specimens is in the range of 0.5–0.8 mm, that is, the specimens fail rapidly after the radial strain reaches 1%–1.6%. The average tensile strength of the BD specimen is 1.749 MPa, and the dispersion coefficient is 0.18. Figure 4 shows the real image and CT scanning image of each specimen after failure. It is found that the main crack of BD specimens is along the loading direction, but the main crack is not straight due to the influence of specimen heterogeneity. The secondary cracks are propagating in the specimen, and most of them are distributed laterally near the middle region of the specimen.

Figure 5 shows the variation curve of transverse strain at four monitoring points on the surface of BD specimens with time. The strain near the center of the BD-1 specimen is firstly in tension and then suddenly changes to compression after about 30 s of loading. This is mainly due to the transverse fracture of the BD-1-1# strain gauge and the change of peripheral strain field. Both BD-2 and BD-3 specimens begin to be compressed in the corresponding 1# strain gauge region and then gradually turn to be in tension. In general, the closer to the bottom of the disc specimen, the greater the tensile strain and strain growth rate. The strain and strain increment of the BD-4 specimen at 4# strain gauge are smaller than those of other measuring points. This is mainly due to the abnormal evolution of strain caused by the fracture zone near the 4# strain gauge at the beginning of loading.

3.2. SCB Test Results

At present, there are many methods to calculate the tensile strength of specimens by the SCB test [2528]. In this paper, equation (2) is used to calculate the tensile strength [25]:where σ is the tensile strength of the specimen, MPa; P is the ultimate load, N; and D and t are the diameter and thickness of the specimen, mm.

Table 3 also shows the peak load, peak displacement, peak time, and tensile strength of four SCB specimens. The average tensile strength measured by the SCB test is 2.102 MPa, and the dispersion coefficient is 0.28. Through the comparison of data in Table 3, it is seen that the tensile strength of the specimen measured by BD tests is 0.83 times of that measured by the SCB method. Previous studies have found that the tensile strength of rock specimens measured by BD tests is about 0.63 times of that measured by the direct tensile method [29]. Therefore, the tensile strength of specimens measured by SCB tests in this paper is about 0.76 times of that measured by the direct tensile method. In other words, the tensile strength of specimens measured by SCB tests is closer to that measured by the direct tensile test.

Figure 6(a) and 6(b) show loading-displacement curves and bottom transverse strain-time curves of specimens in SCB tests, respectively. Except for the abnormal loading-displacement and bottom strain data of the SCB-0-3 specimen, the loading-displacement curves of other specimens are basically the same in terms of development trend. The deformation corresponding to the peak load of the specimen is about 0.2–0.3 mm, which is obviously smaller than the radial peak displacement of BD specimens. The loading-displacement curve and strain-time curve show two stages with different growth rates before the peak value. In other words, the initial stress and strain increment are small, but when the deformation is greater than 0.08 mm, the load and strain increase rapidly.

Figure 7 shows the crack distribution on the front of typical specimens after failure in the SCB test. The difference between the loading-displacement curve and bottom strain-time curve of SCB-0-3 and the other three specimens is mainly caused by the crack initiation position, which is not near the center of the bottom of the specimen. Therefore, the heterogeneity of coal has a certain influence on the SCB test results.

From the failure load of specimens in the two tests, it is found that, under the same loading rate and specimen thickness, the peak load of specimens in SCB tests is less than that of specimens in BD tests (about 15.14%). Therefore, the SCB test can effectively simulate the tensile state of the specimen, and the crack propagation mode in the specimen is simpler; the SCB test is more suitable for measuring the tensile strength of coal.

4. Fracture Toughness Test of NSCB Specimens

There were 12 NSCB specimens in 3 groups for the fracture toughness test, namely, 4 specimens with notch depths of 4 mm, 7 mm, and 10 mm. Table 2 shows the specific parameters. The test results of 11 specimens were obtained (the results of the NSCB-10-4 specimen were not obtained). Figure 8 shows the loading-displacement curve of NSCB specimens with different notch depths in three groups during the fracture toughness test. With the increase in the notch depth, the average peak load of the specimen decreases. The average peak load of the specimen with the notch depth of 4 mm, 7 mm, and 10 mm are 554 N, 363 N, and 255 N, and the corresponding dispersion coefficients are 0.32, 0.26, and 0.35, respectively, as shown in Table 4. The variation trend of these curves is similar, that is, they all increase rapidly before the peak load, then slow down, and then, increase rapidly when the displacement of the loading point reaches about 0.1 mm. Affected by the heterogeneity of coal and the difference of notch depth, the loading-displacement curve and surface strain change of specimens with different notch depths are also different. Figure 9 shows the change of strain at two monitoring points on the surface of three groups of NSCB specimens with different notch depths. It can be found that, with the increase in the notch depth, the strain near the loading point changes from compressive strain to tensile strain, while most of strain is still in the tensile state above the crack end.

At present, the commonly used equation for calculating plane strain fracture toughness KΙC of NSCB specimens is as follows [16]:where the dimensionless notch depth is and R is the specimen radius, mm. Pmax is the peak load of specimen failure, N; Y′ is the dimensionless stress intensity factor. In general, the plane strain fracture toughness of coal measured by the NSCB test decreases with the increase of notch depth. If β = 0.16, the standard deviation of KIC is 2.11, which is significantly higher than that of 1.45 and 1.71 if β = 0.28 and β = 0.40. This also confirms the conclusion of previous studies that the fracture toughness test of NSCB specimens will produce relatively large error if β < 0.2 [12].

The crack growth is affected by the notch depth and the bedding angle. When the notch depth is 4 mm, the crack propagation is easily affected by the bedding; when the notch depth is 7 mm and 10 mm, the crack propagation is approximately linear and is less affected by the bedding. Figure 10 shows the relationship between KIC measured by the fracture toughness test of NSCB specimens and notch depth. It is seen that the depth of notch depth also has a certain influence on the dispersion of plane strain fracture toughness of specimens. By fitting the upper boundary, lower boundary, and median value of KIC data with the change of notch depth a, the following equation can be obtained.

The upper boundary fitting curve is

The lower boundary fitting curve is

The fitting curve of average value is

The results show that the KIC variation of NSCB specimens with the notch depth shows a V-shaped feature, that is, when the notch depth is 7 mm, the dispersion of the plane strain fracture toughness test value is small.

In the 1960s, Rice put forward J-integral theory as the basic index to describe the elastic-plastic fracture parameters of materials. Subsequently, Liu et al. established relevant criteria to replace the stress intensity factor. The main advantages of J-integral theory are as follows: (1) due to the conservation property of J-integral [30], the specific calculation of stress and displacement near the crack tip is not involved, thus avoiding the tedious process of solving the stress field; (2) the influence of the high-stress gradient region near the crack tip is much smaller than that of the stress intensity factor method on the calculation accuracy of J-integral [16, 31]. Therefore, it is more reasonable to use the theory of elastic-plastic fracture mechanics to evaluate the fracture properties of specimens with cracks.

For mode I fracture of brittle materials, the J-integral fracture toughness is [30]where U is the fracture energy of the specimen, 10−3 J.

Generally, the J-integral fracture toughness can be measured by only two specimens with different notch depths. In this paper, the J-integral fracture toughness is determined by the mean value of the test results of three specimens with different notch depths. Firstly, the fracture energy can be obtained by integrating the loading-displacement curves, and then, the J-integral fracture toughness of the specimen can be calculated according to equation (7).

Figure 11 shows the energy absorbed by NSCB specimens with different notch depths during failure. The energy absorbed by the specimen with the small notch depth is greater than that of the specimen with small notch depth. The fracture energy per unit thickness decreases by 0.65 × 10−3 J/mm with the increase of 3 mm of notch depth. The dispersion coefficients of fracture energy per unit thickness of NSCB specimens with notch depths of 4, 7, and 10 mm are 0.36, 0.28, and 0.41, respectively. When the notch depth is 7 mm, the dispersion of plane strain fracture toughness and fracture energy per unit thickness is the smallest. Due to the limited number of specimens, further research is needed to determine the relationship between β and the discreteness of plane strain fracture toughness.

By fitting the upper limit, lower limit and average value of fracture energy per unit thickness in Figure 11, the relationship between fracture energy per unit thickness and notch depth of specimen can be obtained.

The upper boundary fitting curve is

The lower boundary fitting curve is

The fitting curve of average value is

The absolute value of the slope of the average fitting curve is the J-integral fracture toughness of the specimen in the test. The J-integral fracture toughness value of the specimen is 0.215 KJ/m2, while that of the asphalt mixture material is 1.5–3.0 KJ/m [31, 32]. The large difference of the two values indicates that the specimen is brittle and the cohesive force between particles is small [33]. At the same time, the two J-integral fracture toughness values obtained from specimens with three kinds of notch depth have little difference. Therefore, J-integral fracture toughness can be used to evaluate the fracture properties of coal more effectively.

5. Conclusions

(1)Compared with BD tests, the tensile strength measured by SCB tests is about 0.76 times of that measured by the direct tensile test, and it is closer to the result of the direct tensile test. Besides, the corresponding deformation of SCB specimens under peak load is about 0.2–0.3 mm, which is significantly less than the corresponding radial displacement (0.5–0.8 mm) of BD specimens under peak load. Before the peak load, the stress and strain of SCB specimens show two stages of different growth rates. In other words, the initial stress and strain increment are small, but when the deformation is greater than 0.08 mm, the load and strain increase rapidly.(2)Under the same loading rate and specimen thickness, the peak load of the SCB test is less than that of the BD test by 15.14%. Therefore, the SCB test can effectively simulate the tensile state of specimen, and the crack propagation mode in the specimen is simpler; the SCB test is more suitable for measuring the tensile strength of coal.(3)With the increase in notch depth, the average peak load and plane strain fracture toughness of the NSCB specimen decrease continuously, and the strain near the loading point of specimen gradually changes from compressive strain to tension strain, while most of them are in tension above the crack end. If β < 0.2, the dispersion of KIC value of plane strain fracture toughness is larger; while if β = 0.28, the dispersion of KIC value of specimens is the smallest.(4)The plane strain fracture toughness KIC values of the tested specimens are 6.30–7.72 , and the average J-integral fracture toughness is 0.215 KJ/m2. The plane strain fracture toughness KIC of coal is obviously affected by the notch depth and the dimensionless stress intensity factor Y′, while the J-integral fracture toughness calculation method is simple. The J-integral fracture toughness calculated by different notch depths has little change. Therefore, the J-integral fracture toughness can be used to effectively evaluate the fracture properties of coal rock and other brittle materials.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by the research fund of the National Natural Science Foundation of China (51604093), the Fundamental Research Funds for the Universities of Henan Province (Grant no. NSFRF200332), the Henan Key Laboratory for Green and Efficient Mining and Comprehensive Utilization of Mineral Resources (Henan Polytechnic University) (KCF201804), the Key Scientific Research Project Fund of Colleges and Universities of Henan Province (21A610005 and 20B440001), and the Doctoral Foundation of Henan Polytechnic University (B2019-22). All support is greatly appreciated.