Abstract

Deformation and fracture of brittle materials, especially crack propagation, have drawn wide attention in recent years. But dynamic crack propagation under impact loading was not well understood. In this paper, we experimentally tested Brazilian disk (BD) fine sandstone specimens containing pre-cracks under cyclic impact loading by the Φ 74 mm diameter split Hopkinson pressure bar (SHPB) test device. The pre-cracked specimens were named central straight through crack flattened Brazilian disk (CSCFBD). By using the low air-pressure loading conditions (0.1 MPa, equal to the impact velocity of 3.76 m/s), a series of dynamic impact tests were detected successfully, and the effects of pre-cracks on dynamic properties were analyzed. Experimental results show that the multiple cracks mostly initiate at/or near the pre-crack tips and then propagate in different paths and directions varying by inclination angles, leading to the ultimate failure. Compared to static or quasi-static loading, dynamic crack propagation and fracture behavior are obviously different. Furthermore, we characterized the crack propagation paths, directions, and fracture patterns and discussed the influences of the pre-cracks during the breakage process. We concluded that the results obtained are significant in investigating the failure mechanism and mechanical properties of brittle materials under impact loading.

1. Introduction

Dynamic deformation and fracture of brittle materials are complex processes. Mining, tunnel excavation, and natural disasters such as landslides and earthquakes are all involved in the problems of dynamic damage. The failure of brittle materials is usually associated with crack propagation initiated from natural or artificial pre-existing defects. Moreover, the mechanical properties of brittle materials are closely related to external loading conditions, such as loading rate and load magnitude [13]. Rock is one of the most complex brittle materials containing different scale voids, cracks, and other defects, as shown in Figure 1, which are the main mechanical factors that affect the rock deformation and failure [4]. Crack growth and catastrophic failures initiated from pre-existing defects subjected to multiaxial loads are the main concerns for geotechnical engineers and designers of underground structures. The defects in rock can promote the initiation of new defects, which in turn may propagate and coalesce with other defects, and then can further decrease the strength of the rock mass [58]. The presence of pre-cracks may obviously reduce the fracture toughness, dynamic uniaxial compressive strength, and dynamic tensile strength and lead to fragmentation and multiple crack interactions, branching, and coalescence [9].

Due to the difficulties of in-situ tests, the laboratory experiment is an important and effective research method to investigate rock failure modes and fracture mechanisms. Over the past few decades, many experiments have been devoted to use semi-circular core in the three-point bending (SCB) specimen [10], Brazilian disk (BD) specimen with chevron flaws or other pre-existing flaws [11], radial cracked ring (RCR) specimen [12], and modified ring (MR) specimen [13] to investigate fracture toughness and crack propagation. Note that, in previous studies, the fracture behavior and crack propagation of brittle materials were mainly investigated under static/quasi-static loading or used in intact specimens. Irwin et al. [1416] divided the simple cracks into three types in the basic failure process: Mode I (the tension/opening mode) crack, Mode II (the sliding mode), and Mode III (the tearing mode) crack. In engineering, Mode I cracks are of prime importance. More complicated cracks can form from these simple cracks, called mixed-mode cracks. Many experimental studies have been conducted to explain the crack initiation, propagation, and coalescence in pre-cracked brittle materials under static or quasi-static loading [1721]. The BD test is an effective way to study crack propagation and fracture behavior of brittle materials [22, 23]. Al-Shayea [24] studied crack propagation paths in pre-cracked limestone central straight though crack Brazilian disk (CSCBD) specimens loaded with diametrical compression, and they investigated the possibility of using outcrop specimens to estimate the fracture toughness behavior of the reservoir rock at in-situ conditions of temperature and confining pressure. Haeri et al. [2527] used Portland Pozzolana Cement (PPC) experimentally and studied crack propagation and crack coalescence of BD pre-cracked and pre-holed specimens under quasi-static compressive and tensile loading. Aliha and Bahmani [28] investigated fracture toughness under mixed-mode loading using different cylindrical and disc shapes for the brittle material. However, in deep strata, the damage and destruction of rock mass are often caused by dynamic even cyclic dynamic loading. None of these investigations listed above has ever captured the dynamic multiple crack propagation. In such cases, it is of great significance and essential to investigate the dynamic crack propagation and fracture patterns of rock in deep strata under impact/cyclic impact loading.

Various experiments devices have been used to explore a wide range of strain rates. Split Hopkinson pressure bar (SHPB) technique, which decouples cleverly the inertia effect in structures and strain rate effect in materials, has been widely used to characterize the dynamic performance of various engineering materials at high strain rate, such as rock [2933], concrete [3438], and ceramics [39, 40] at high strain rates (102∼104 s−1). The strain rate sensitive behavior of brittle materials has been under investigation for several decades. The strain rate sensitivities are mainly measured by strength or the strains at the maximum stress [41, 42]. In recent years, a variety of researchers have investigated and demonstrated the dynamic properties of natural or artificial brittle materials [4347], and they found that the dynamic strength (including dynamic compressive strength and dynamic tensile strength) and impact toughness increases with strain rate and the strain rate sensitivity of brittle materials.

In this paper, we defined some central straight through crack flattened Brazilian disk (CSCFBD) specimens to investigate dynamic crack propagation and fracture patterns under impact loading by the SHPB device, and the term “pre-crack(s)” is used to describe the artificially created crack. The motivation for this work will focus on the following two points: (1) investigating the dynamic fracture patterns, multiple crack propagation paths, and directions in pre-cracked specimens; (2) characterizing and analyzing the crack types initiated from pre-cracks under cyclic impact loading. This paper is organized as follows: in Section 2, the preparation of tested CSCFBD specimens and experimental procedure are discussed. The SHPB experimental scheme is also introduced. In Section 3, we present the experimental results. We look at both the multiple crack propagation paths and directions. In Section 4, we discuss the experimental results and characterize the dynamic fracture patterns and failure modes. The conclusions concluded upon the foregoing results are given in Section 5.

2. Preparation of Disk Specimens and Experimental Procedure

Due to the crystalline and blocky structures, fine sandstone is widely used to investigate the fracture behavior and crack propagation of defected brittle materials. In this paper, we used fine sandstone to study the dynamic crack propagation under impact loading. The fine sandstone samples were excavated by geologic drilling from about 900-meter depth underground strata in Juye coalfield, whose Cenozoic formation is very thick. The average thickness of strata in the fourth system is 158.43 m, and the average thickness of the upper tertiary strata is 497.01 m. And the thickness of the new boundary layer is 530∼720 m, mainly composed of clay, sandy clay, sand, fine sand, and gravel. Main coal seam roof and floor sandstone thickness is 4.80∼75.65 m, mainly fine sandstone, local sandstone, and siltstone [4]. The pre-cracked specimens’ preparation and experimental procedures will be explained as follows.

2.1. Preparation of CSCFBD Specimens

For manufacturing CSCFBD specimens, the whole tests used six samples with each one cored into 30∼50 cm long, diameter D = 62 mm cylindrical columns in the construction site. In order to avoid external environmental influence, the surface of the samples was wrapped in the multilayer food preservation film after removing from the formation. The specimens were cut into BD shapes with the size of Φ 62 × 30 mm. Three specimens with same pre-crack geometry were prepared to guarantee the reproducibility of these experimental tests, and a total number of 21 CSCFBD specimens were manufactured in this work. The physical and mechanical properties of the tested fine sandstone are listed in Table 1.

Various Brazilian tests were conducted on CSCFBD specimens containing a single crack with different inclination angles. These pre-cracks were created by the high-speed water jet cutting machine [4], as shown in Figure 2. Figure 3 shows pre-cracked specimens with different inclination angles and the crack inclination angle β (the angle between the normal line of the pre-crack surface and the vertical direction of the impact load), β = 0°, 15°, 30°, 45°, 60°, 75°, and 90°. Figure 4 shows a schematic view of CSCFBD specimens. The pre-cracks length, 2b (b is the half of pre-crack length), is equal to 10 mm, and the crack width is 1 mm. The radius and thickness of the CSCFBD specimen are R = 31 mm and H = 30 mm. Crack length ratio is an important parameter for the pattern, trajectory, and the number of fractures; in this work, the ratio b/R is 0.16.

2.2. SHPB Experimental Procedure

In this work, the dynamic tests were taken in the structure laboratory of Hohai University, adopting Φ 74 mm diameter straight taper variable cross sections SHPB device. The Φ 74 mm diameter SHPB device mainly composes of the power system (which propelled by a gas gun), elastic bars (an incident bar, a transmitter bar, and an absorbing bar), a damper, energy absorbing setup, high dynamic strain indicator and data processing systems. The power system consists of an air compressor and pressure vessel. The impact velocity is measured by the light electric tachometer. Schematic of the SHPB test device is shown in Figure 5.

The SHPB device is based on the one-dimensional theory of elastic wave and uses Lagrangian coordinates to describe all the physical parameters. The control equations of waves in the bars are on the following assumptions. One-dimensional assumption: speed and the strain ε are only the function of point X and time t and the assumption that the strain rate is independent. Assuming axial wave propagation and homogeneous stress distribution in the specimen, the resulting stress , strain , and strain rate of the specimen are obtained by the following equation [48]:where , , and are the cross-sectional area (mm2), Young’s modulus (GPa), and the wave velocity (km/s) of the bar material and and are the length (mm) and cross-sectional area (mm2) of the specimen. , , and are the strain singles in specimen. The diameter, Young’s modulus, and density of the elastic bars are 74 mm, 210 GPa, and 7850 kg/m3, respectively.

In fact, during the experimental tests by BD specimens, the cracks of samples are first produced by the center and then expanded along the radial direction. Wang et al. [49] presented that if the samples between two planes parallel to the plane of the degree of smoothness and not less than 0.05 mm, then they can ensure that the fracture initiated from the specimen center. The loading areas corresponding to the center angle 2α to meet 20° ≤ 2α ≤ 30°. Figure 6 shows the flattened BD specimen under radial loading. Before the installation of the specimens, a layer of Vaseline is evenly applied to the contact between the end of the specimen and the end of the compression bar, and the specimens need to be tightened between the incident bar and the transmitter bar. During SHPB tests, the brittle materials may fail before stress uniformly is achieved within the specimens. Modification of the incident pulse to closely match the elastic response is required. The pulse shaping technique has been widely applied in SHPB testing of engineering materials, and it is especially used for investigating the dynamic response of brittle materials. In this paper, we used a pulse shaping technique, a thin copper disk (12 mm-diameter and 1 mm-thickness) which was placed on the impact side of the incident bar. The pulse shape can attenuate high-frequency oscillations of the incident stress wave to improve the stress wave shape. Pulse shaping technique allows for controlling the damage of brittle materials [29, 30].

3. Experimental Results

3.1. Impact Time Analysis of CSCFBD Specimens

During testing, the air pressure in the gas gun chamber was kept constant 0.1 MPa, equal to the impact velocity 3.76 m/s, and the bullet must be brought back to its original position before the next loading. The times of impacts were recorded for each specimen until final failure. As shown in Figure 7, to reach the final crack propagation paths forms, 6 times impact loading are needed when β = 0° and the impact times are 3, 4, 6, 4, 2, and 2 when β = 15°, 30°, 45°, 60°, 75°, and 90° respectively. It is obvious that the impact times in final failure are related to the pre-existing inclination angles. Because pre-cracks can make the specimen strength to reduce, so when β = 15°, 75°, and 90°, the impact times are usually less than the cases of β = 0°, 30°, 45°, and 90°.

3.2. Dynamic Crack Propagation Paths and Directions of CSCFBD Specimens

In this paper, we investigated the dynamic crack propagation paths and directions in CSCFBD specimens. In Figure 8, the crack propagation paths in CSCFBD specimens with different inclination angles are indicated, and for all figures, the upper plane edges of the specimens are contacted with the incident bar. As shown in Figures 8(b)8(g), cracks initiated from the tips of pre-cracks and approximately propagated towards the direction of the maximum stress. It should be noted that, in Figure 8(a), there is a crack initiated from the middle portion of the pre-crack when β = 0°.

Figure 8(a) shows that five crack propagation paths appeared: two cracks started from the pre-existing crack’s left-end tip and propagated to the specimen’s lower plane edge; two cracks started from the pre-existing crack’s right-end tip; the rest one crack propagated to the upper and lower plane edges. However, because of the high values of crack orientation with respect to the loading direction, there is one crack which did not propagate from the tip of the pre-crack. In Figures 8(c) and 8(d), β = 30°, 45° respectively, five crack propagation paths appeared. Figure 8(c) shows that there are two cracks that started from the pre-crack upper tip and propagated to the upper plane edge, two cracks started from the pre-crack lower tip and propagated to the lower plane edge, and the rest one crack started from the pre-crack upper tip and propagated to the lower plane edge. While in Figure 8(d), there is one crack that started from the pre-crack lower tip and propagated to the upper plane edge. In Figures 8(b), 8(e), and 8(f), β = 15°, 60°, and 75°, respectively, four main crack propagation paths appeared, and two of them started from the pre-crack upper tip and propagated to the upper plane edge. The other two started from the pre-crack lower tip and propagated to the lower plane edge. In Figure 8(b), an intermittent crack appeared, starting from the pre-crack lower tip, and has a downward trend to the lower plane edge. Figure 8(e) shows that the scatter of the four cracks presents regular symmetrical characteristic. As shown in Figure 8(g), β = 90°, and the upper pre-crack tip propagated two crack propagation paths to the upper plane edge and one to the lower plane edge.

From Figures 8(a)8(g), it can be clearly seen that the final dynamic crack propagation paths and fracture patterns under impact loading are obviously different compared to those under static or quasi-static loading (e.g., Figure 9 [4]). In many studies, there is only one fracture that propagates from each tip of BD specimens under static or quasi-static loading. But under impact loading in the SHPB test, due to fractures with high strain rate, there are multiple crack propagation paths. The results also show the pre-existing inclination angles affect the multiple crack propagation paths and directions initiated from tips of the pre-cracks under cyclic dynamic loading.

4. Discussion

The crack propagation patterns were obtained in previous investigations of brittle materials with pre-cracks as shown in Figure 10 [1, 2]. From the experimental results above, two types of cracks were observed: wing cracks and secondary cracks. Usually, the wing cracks are tensile cracks and the secondary cracks are shear cracks (oblique shear cracks and coplanar or quasi-coplanar shear). Table 2 summarizes crack types initiated from the pre-cracks in CSCFBD specimens under cyclic impact loading. Most of the tensile cracks and shear cracks initiated from the tips of pre-cracks at an angle and then propagate to the parallel to the compressive direction. Shear cracks’ initiation patterns depend on the inclination angles of the pre-cracks. Under impact loading, shear cracks caused the failure of the tested specimens mostly.

Fracture and failure of pre-cracked CSCFBD specimens under cyclic impact loading involve tensile and shear crack types. Because of pre-cracks, the fracture patterns and failure modes are much more complex than those of intact brittle materials. Note that all the CSCFBD specimens present multiply crack propagation paths. But the geometries of pre-cracks appear to play limited effects on the crack type of CSCFBD specimens under cyclic impact loading.

In this paper, we studied the dynamic crack propagation and fracture patterns in pre-cracked CSCFBD specimens from deep underground strata under cyclic impact loading. The dynamic crack propagation paths and directions are obviously different from those under static or quasi-static loading in the previous studies, and the experimental results present some regular symmetrical characteristics. The dynamic crack propagation paths and propagation directions are not only determined by material related but also dependent on the geometries of the pre-cracks. However, the geometries of pre-cracks appear to play limited effects on the final dynamic fracture patterns and failure modes of CSCFBD specimens under cyclic impact loading. Further research may be a focus on the dynamic crack propagation mechanism of brittle materials studied by experimental and numerical simulation comprehensively. Various numerical simulation methods, e.g., finite element method (FEM), extended finite element method (XFEM), finite differential method (FDM), and different criteria, e.g., strain energy density (SED) criterion, cohesive zone model (CZM) [50] for theoretical analysis, have been developed to investigate crack propagation, fracture patterns, and failure modes in brittle materials but have not shown satisfactory effectiveness in modeling dynamic damage evolution and crack propagation [5153]. Using mesh-free particle methods to numerically study dynamic crack propagation may be an effective way. Especially in peridynamics, cracks are part of the solution, not part of the problem [54]. Therefore, using peridynamics to simulate dynamic multiple crack propagation and fracture patterns in brittle materials under impact or cyclic impact loading would be very meaningful which we plan for in the future.

5. Conclusions

In this paper, we provide experimental results of dynamic fracture tests carried out on fine sandstone CSCFBD specimens under cyclic impact loading by the Φ 74 mm-diameter SHPB device. Dynamic crack propagation and fracture mechanism are rather complicated processes. From the results presented and analyzed above, the following conclusions could be drawn:(1)Compared to static or quasi-static loading, the dynamic crack propagation and fracture behavior are much more complex, and it presents multiple crack propagation paths and directions in dynamic fracture.(2)Both tensile and shear cracks were observed, most of them mainly initiated from tips of the pre-cracks and propagated in a stable manner. According to the different geometries of pre-cracks, the tested CSCFBD specimens experienced tensile or shear crack propagation failure.(3)The natural or artificial pre-existing defects can change the crack propagation paths, especially directions under impact loading compared to static or quasi-static loading. However, the geometries of pre-cracks appear to play limited effects on cracks type of CSCFBD specimens under cyclic impact loading.(4)Dynamic crack propagation and fracture mechanism are rather complicated processes. This study revealed multiple dynamic crack propagation behavior that had not been observed previously. Numerical simulation is further needed to be summarized and explored which we plan for in the future.

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

The research described in this paper was financially supported by the National Key Technologies Research & Development Program (nos. 2018YFC0406703 and 2017YFC1502600), the National Natural Science Foundation of China (nos. 11672101 and 11372099), the Postgraduate Research & Practice Innovation Program of Jiangsu Province (no. KYLX16_0700), and the Fundamental Research Funds for the Central Universities (no. 2016B45214).