Research Article  Open Access
Diyuan Li, Feihong Gao, Zhenyu Han, Quanqi Zhu, "Full and LocalField Strain Evolution and Fracture Behavior of Precracked Granite under Coupled Static and Dynamic Loads", Shock and Vibration, vol. 2020, Article ID 8866673, 15 pages, 2020. https://doi.org/10.1155/2020/8866673
Full and LocalField Strain Evolution and Fracture Behavior of Precracked Granite under Coupled Static and Dynamic Loads
Abstract
Flaws and discontinuities play a crucial role in the failure process of rocks. To investigate the fracturing mechanism of rock with combined flaws composed of crack and hole, the digital image correlation (DIC) method is used to record and analyze the rock failure behavior. Coupled static and dynamic loads are applied on granite specimens with prefabricated flaws by a modified split Hopkinson pressure bar (SHPB) device. The dynamic mechanical properties of the granite specimens are affected by the flaw inclinations with the loading directions. With the inclination angle increasing, the combined strength and peak strain both decrease first and then increase. Full and localfield strain evolution of the granite specimens is analyzed in a quantitative way by using DIC technique. The specimens with a flaw angle of 45° are broken relatively evenly with homogenous small particle sizes. The variation trend of fragment sizes is consistent with that of combined strength and absorption energy of the specimens.
1. Introduction
The geometry, size, and spatial distribution of flaws incubating within the rock are of particular importance to the mechanical properties of rock [1]. Under the longterm effect of geological action and external force, numerous micro and macroflaws, such as faults, discontinuities, holes, and microfissures, have grown inside the rock. The existence of these natural weak surfaces and flaws is the most fundamental cause of rock failure. In general, exploring the fracturing mechanism of rock with flaws has been one of the focuses of rock mechanics in recent decades [2, 3].
In engineering practice, rock media exhibit very complex characteristics due to the heterogeneity, nonlinearity, and uncertainty properties [4], and usually contain both cracklike and holelike flaws. The mechanical properties and fracturing behavior of specimens with one of them have been systematically investigated by using different experimental and numerical methods. It has been concluded that the fracture initiation, propagation, and coalescence behavior of specimens with holelike flaws are influenced by a multitude of factors, such as hole size [5], hole shape [6, 7], number of holes [8, 9], and stress conditions [10, 11]. Meanwhile, those of specimens with cracklike flaws are likewise closely related to the geometry of cracks [12], crack numbers [13, 14], and crack size [15].
Generally, rock mass contains massive flaws in irregular arrangements and of different shapes in rock engineering. The crackhole combined with flaws is also a common type of flaws in rock mass. Yang et al. [16] conducted experiments to investigate the effect of fissure angle on the strength and deformation behavior of sandstone specimens with combined flaws. By using photographic monitoring observation, crack coalescence mechanism (tensile crack, hole collapse, farfield crack, and surface spalling) and nine different crack types for the combined flaws are identified. Liu et al. [17] carried out a number of laboratory tests and numerical simulations based on the particle flow code method to analyze the cracking process and stress field evolution in specimens containing combined flaws. The arrangement of the flaws in their studies is interconnected cracklike and holelike flaws. Yin et al. [18] conducted experiments and numerical simulations to investigate the cracking characteristics of specimens with the flaws in a noncontact geometry of combination of single hole and an inclined fissure underneath. However, most of these tests were usually loaded by static uniaxial compression for those studies aforementioned.
In recent decades, scholars have carried out numerous tests based on the SHPB system, to investigate the mechanism of crack initiation and propagation processes on brittle rock. Under individual dynamic loading, the mechanical properties and crack propagation of rock specimens containing a single flaw with different inclination angles are investigated [19–21]. Under coupled staticdynamic loading, the influence of dynamic strain rate and the axial static stress on dynamic behaviors of rock specimens with a hole or fissures are considered by some researches [22–24], while the study focusing on holecrack flaws is rather limited. In view of this, further study of the fracturing properties of brittle rock materials containing prefabricated holecrack flaws under coupled staticdynamic loading is essential.
Conventional sensors such as strain gauges, extensometers, and accelerometers are widely used to monitor the deformation of rock specimens. However, these sensors can only measure the strain at a few discrete locations [25–27], and it is obviously not applicable to obtain the local strain surrounding the flaws and cracks in the specimen. To overcome the limitations of traditional contact measurement technology, the digital image correlation (DIC) method, a typical representative of optical technology, has performed its great advantages [28–30]. Fracture is of great importance for the behavior of rock engineering structures [31–33]. Hence, the crack initiation and fracturing process of specimens with prefabricated flaws have been investigated with different auxiliary experimental means, such as acoustic emission (AE) technology [34], computerized tomography (CT) scanning [35, 36], and 3D Xray computed tomography [37]. All those mentioned methods are beneficial to provide inspections for rock failure, especially under static compression tests. However, these methods have limitations in recording the strainfield variation in the dynamic experiment.
Besides, there are also many numerical simulation methods. For example, Li et al. [38] adopt the damage smear method (DSM) to investigate transscale progressive rock failure process. Damage evolution and accumulation of representative volume elements (RVEs) are used to reveal the characteristics of macroscopic rock failure. In addition, the finite element method (FEM) [39], discrete element method (DEM) [9, 40], rock failure process analysis code (RFPA) [41], and particle flow code (PFC) [8, 42] are also widely used to reveal the failure mechanism of rock. However, to conduct numerical modeling, reliable and reasonable constitutive relationship must be built as the basic information. Besides, the parameter selection and validation of these numerical models also require experimental results. It is a challenge in directly obtaining the full and localfield response, especially measuring true dynamic response of the specimen or structure. Based on the calculation of the gray matrix of the digital speckle images, the surface strain evolution of a specimen can be extracted by using the DIC technique without the need to develop a numerical model.
In the current work, a modified split Hopkinson pressure bar (SHPB) equipment is used to apply the coupled static and dynamic loads to a rock specimen. Each rock specimen is precracked with crackhole flaws in a noncontact geometry. Digital image correlation (DIC) apparatus is adopted to record the realtime fracturing process of specimens. The fullfield cracking process and strain localization evolution are analyzed. The energy dissipation and rock fragmentation size characteristics are summarized. The main purpose of this paper is to obtain the mechanical properties and the strain evolution of brittle rock containing combined flaws under dynamic loading by using the DIC technique. Evidently, the dynamic rock failure analysis by the DIC method will provide some visual observations and proofs for the studies in rock mechanics, particularly for the rocks with natural or artificial flaws.
2. Methodology
2.1. Specimen Preparation
The specimens adopted in the experiment were made of granite with high strength and brittleness, which are extracted from Miluo county of Hunan province in China. To determine the basic mechanical properties of the test granite, cylindrical specimens with dimensions of Φ50 mm × 100 mm (diameter × height) were prepared and tested under uniaxial compression. The basic mechanical parameters are listed in Table 1.

Rectangular cuboid specimens (60 mm in length, 45 mm in width, and 15 mm in thickness) were tested under the coupled static and dynamic loads by a modified split Hopkinson pressure bar (SHPB) apparatus. A circular hole (10 mm in diameter) was cut by using a highpressure water jet apparatus at each specimen center. A contrast group with only circular hole (Group C) and an experimental group with crack unexposed to the hole (Group U) were prepared and tested afterwards. The specific geometrical layout about the test specimen can be seen in Figure 1. The prefabricated combined flaws contain a holelike flaw and a cracklike flaw 15 mm in length in different inclination angles β (the angle between the cracklike flaws and the loading direction, including 0°, 30°, 45°, 60°, and 90°). The artificial crack has a minimum distance d of 7 mm from the edge of the internal hole, and the layout is displayed in Figure 2. The prefabricated crack is always tangent to a circle with a radius of 12 mm (R = d + r). Tangent point, which is also the midpoint of the crack, rotates counterclockwise in the second quadrant. The angles aimed to investigate in this paper are limited in only one quadrant, and larger angles are not discussed due to symmetry.
(a)
(b)
(c)
Moreover, the dimensions and mechanical parameters of each type of specimens are listed in Table 2. It is noteworthy that the specimen no. in Table 2 describes the number of crack, crack inclination angle, and the specimen number (three specimens are prepared for each type). For example, U302 refers to the second specimen (no. 2) containing a single crack with β = 30°.
 
Note: is the axial static pressure, is the dynamic strength, and is the combined strength. The symbol “—” indicates that the data are eliminated for the specimen owing to large deviation. 
Before the tests, artificial speckle patterns should be prepared on the specimen surface. At first, a layer of white paint was sprayed thinly and evenly on the surface of specimens. If the paint is too thick, the crack initiation and propagation of the specimens are difficult to be reflected on the surface in time, which would result in the inaccurate acquisition of images captured by the highspeed digital camera (Phantom V711, produced by Vision Research Inc.). After the paint was dry, speckles were coated manually and randomly, and each speckle was set to cover approximately 4–6 pixels. It is noteworthy that a speckle particle size is the same as the distance between two adjacent speckles to ensure that it can be caught and optimized by computer simulation.
2.2. Testing Equipment
2.2.1. Loading System
Rectangular cuboid granite specimens were tested by using a modified split Hopkinson pressure bar (SHPB) system with coupled static and dynamic loads, as shown in Figure 1. The modified SHPB system consists of a gas gun, a coneshaped striker, an incident bar, a transmitted bar, a momentum bar, a damper, and a set of axial static pressure device [44]. The basic parameters of elastic bars in SHPB device are shown in Table 3.

2.2.2. Digital Image Correlation (DIC) System
During the test, DIC apparatus is used to monitor the deformation evolution of the specimens and it mainly consists of the highspeed digital camera and software system. The installation of this apparatus is shown in Figure 3. A gooseneck lamp was installed perpendicularly to the specimen to provide illumination. Once triggered, all images were capture and simultaneously stored in a computer for the subsequent processing.
The lengthpixel ratio of the imaging system is 0.234 mm/pixel. The frame rate of the highspeed camera is set to 79161 fps (frames per second) with a resolution of 256 × 256 pixels so that it can capture an image each 12.63 μs. After impact tests, the captured speckle images are analyzed by twodimensional visual image correlation (VIC2D) software from Correlated Solutions, Inc. (CSI). (Please refer to reference [46] for more details about the software). Generally, VIC2D is advantageous for monitoring inplane deformation and strain measurements of any planar surface on account of its simple setup. However, there are some limitations for the DIC system: (a) measurement accuracy mostly depends on the quality of speckle particle. (b) Measurements are limited to inplane surface deformation and strain. (c) Once the strain is less than 100 µε, DIC is not recommended for any applications [47].
2.3. SHPB Test
The main loading stages of the test are as follows. First, static loads of 20 MPa were applied manually on the specimen in a slow and uniform speed by the axial static pressure loading device, and the static loading rate is 0.17 MPa/s. Second, dynamic impact loads were applied to the granite specimen by launching the coneshaped striker in the SHPB system. The dynamic strain rate is in the range of 50∼60 s^{−1}. In the test, once a transistorlogic (TTL) signal was generated by the strain gauge attached on the incident bar, the highspeed camera and the oscilloscope were both triggered simultaneously.
3. Testing Results and Analysis
3.1. Testing Validity
During the test, a slowly rising half sine wave is generated by launching the coneshaped striker to impact the incident bar. When the incident wave propagates to the barspecimen interface, a part of the wave is reflected (reflected wave), and the remainder passes through the specimen and is transmitted into the transmitted bar (the transmitted wave). Hence, stress equilibrium at the two ends of the specimen should be achieved in the test to ensure that results are convinced and reasonable. It means before the rock specimen reaches failure, the transmitted stress wave is almost equal to the superposition of the incident wave and reflected wave.
Taking specimen U03 as an example, the stress waves of the specimen are shown in Figure 4. The transmitted stress wave (Tra) and the superposition of incident wave and reflected wave (Inc + Re) in the middle part almost coincide with each other, indicating that the stress equilibrium of the specimen is satisfied during the dynamic test.
3.2. Dynamic StressStrain Curve
The time when impact stress wave propagates to the specimen end is set as the starting point in the test. Figure 5 shows the dynamic stressstrain curves of the granite specimens with prefabricated flaws. It can be approximately classified into three typical stages: (I) linear elastic deformation stage, (II) nonlinear deformation stage, and (III) macroscopic failure stage.
Evidently, the compaction stage which should be concave upward cannot be observed in this curve. Because the initial compaction of microfissures in the interior of rock specimens has achieved in prestatic loading, the curve goes directly to the stage of linear elasticity. The first stage has a quite short duration. With further loading, the stressstrain curve exhibits nonlinear property. New crack initiation and propagation occur in this stage.
When reaching the peak stress, the stressstrain curve shows a downward trend, which may result from the propagation of macrofractures. Hence, the postpeak stage is defined as macroscopic failure stage. In this stage, macroscopic fractures form and propagate to the specimen ends, finally causing the brittle failure.
3.3. Mechanical Properties
According to the onewave propagation theory in the SHPB test, the axial stress , strain , and strain rate can be derived by the following equations [48]:where are the crosssectional area, Young’s modulus, and Pwave velocity of the elastic bar, respectively. and are the crosssectional area and length of the specimen, respectively. , , and are the incident, reflected, and transmitted strain pulses, respectively.
Thus, the dynamic mechanical properties of the granite specimens can be obtained. The results are listed in Table 2. The strain rate of the specimens is maintained at about 50–60 S^{−1}, which can be regarded as a constant strain rate range. Dynamic strength reflects the impact resistance of the specimen, and the peak strain refers to the strain corresponding to the dynamic strength. Moreover, the combined strength is the sum of the prestatic pressure and the dynamic strength, which reflects the actual bearing capacity of the specimens [44]. The average value and the standard deviation of the testing results are provided in Figure 6.
The existence of artificial flaws are responsible for strength weakening [49]. Combined strength and peak strain of the specimens in Group U both decrease first and then increase with the inclination angle increasing gradually. The strength reduction of specimens with an artificial cracklike flaw of 0° is only 0.66% lower than that of specimen C0, while the strength reduction of specimens with a prefabricated flaw of 30°, 45°, 60°, and 90° is 24.40%, 30.79%, 25.17%, and 16.57%, respectively. The peak strain of specimen C0 is 0.32%, whereas it ranges between 0.177% and 0.204% for the specimens with combined flaws, which has a significant reduction.
3.4. Fracturing Process
During the past few decades, white patches are regarded as a precursor to observe new cracks in many studies [50]. Moreover, the initiation, extension, and coalescence of fractures are described in a qualitative way. In this study, based on DIC technique, the major principal strain contours of the specimen are obtained during the coupled staticdynamic loading, as shown in Table 4. It can reflect the evolution of the surface strain field in the process of loading intuitively and quantitatively.

In the evolution of strain field, different types of specimen behave quite differently. It demonstrates the geometries of combined flaws significantly affect the deformation characteristics [18], which are summarized as follows:(1)In static loading pattern, the whole specimen is in the state of compression. However, some local high strains can be observed prior to observable cracks, which usually occur around the combined flaws. These local high strains are similar to the white patches mentioned earlier [51], which could establish foundation to reveal the new fracture trajectories and initiation mechanisms in the subsequent dynamic loading.(2)In dynamic loading patterns, at early stage of dynamic loading, the high strain shows obvious migration characteristics with the propagation direction of the stress wave (specimen C0 is not involved). The incident wave enters from the barspecimen interface, and soon afterwards high strain concentrations are formed in the right half of specimens. Next, they move forward along the propagation of the stress wave and then gradually distribute near the combined flaws of the specimens. With further loading, the high strain zones are widened and propagated, accompanied with some secondary cracks.(3)Two types of fractures are commonly observed: tensile fractures and shear fractures. Tensile fracture is nearly parallel to the loading direction, while shear fracture is oblique with respect to the loading direction [52–54]. The tensile fractures mostly occur at the flaw tip, while a few of them occur at the side wall of the circular hole. Shear fractures which formed at the top and bottom of the hole extend along the diagonal directions at an angle with the loading direction, generating an X or a half Xshaped fracture, which dominate the rock failure process during the dynamic loading.(4)Fractures emanating from the top or bottom of the hole coalesce with those initiating from crack tips. The crack tip which is near with the internal hole has the precedence. The coalescence modes of five types of flaw geometry will be discussed in Section 4.1 in detail.
3.5. Strain Evolution Analysis of the Measurement Points
In this paper, in order to obtain the strain localization evolution processes, some key measurement points are selected at discrete locations on the specimen surface, as shown in Figure 7. DIC technology is used to extract the strain data during the failure process, and then strain analysis and comparison about these points are conducted. In specimen C0, A, B, C, and D are, respectively, located at the top, left side, bottom, and right side of the hole. E, F, G, and H are, respectively, located at the midpoints of the four lines from measurement points A, B, C, and D to the edges of the specimen. In the specimens of U0, U30, U45, U60, and U90, A and B are, respectively, the left and right tips of the crack; C is the triangular vertex with line AB as the base; D and E are the top and bottom of the hole; F and G are the left and right wall of the hole; H, I, and J are equidistant on the median of the triangle. Taking typical specimens as examples from six groups, the straintime curves of measurement points are plotted in Table 5, including the maximum principal strain e1 and the minimum principal strain e2.

(a)
(b)
(c)
(d)
(e)
(f)
For specimen C0, it can be seen from Table 5(a) that when t = 88 μs, the peak strength is 170.3 MPa, and the local strain of all measurement points is in a low fluctuation. After that, the maximum principal strain at the top and bottom of the hole begins to increase significantly, and tensile cracks appear along the loading direction (Table 4). The straintime (e2) curve in Table 5(b) shows that there are shear fractures perpendicular to the loading direction appearing firstly at the top and bottom of the hole (A, C), and then they are restrained gradually. At t = 200 μs, the value of e2 drops sharply to a small value but later on increases again. At this moment, the specimen is completely destroyed, and shear fractures dominate the failure process.
For the specimens with prefabricated cracks, the statistic results are summarized in Table 6, and the variation rules are summarized as follows:(1)The initiation of new fractures mainly occurs at crack tips (A, B), and the top and bottom of holes (D, E). However, the strain of C has an obvious rise only at the inclination angle of 0°, 30°, and 45°. That is mainly because point C is in a quite close position to point D (hole top) at those angles. The strain of F, G, H, I, and J is in low fluctuation.(2)The e2 value at point A increases significantly, while the e1 value increases in a relatively low level. The e1 value of point B increases rapidly, while the fluctuation of e2 value is weak. It illustrates that shear stress is mainly concentrated at the left crack tip (A), while tensile stress is mainly concentrated at the right crack tip (B). At the angle of 90°, the strain at the top and bottom of the hole is significantly higher than that at the two tips of the crack. For this condition, new cracks are more easily occurred around the hole rather than at the crack tips. As the circular hole is completely located in front of crack and exposed to the stress wave earlier and larger than the preset crack. Therefore, the roof and floor of the internal hole could become the preferable location for microcracks to gather and develop, rather than at the tip of the crack.(3)The e1 and e2 values of point D (hole top) increase significantly at the angles of 0° and 90°, while those two values are relatively low at the angles of 30°, 45°, and 60°. The e1 and e2 values of E (hole bottom) are both large in the angle range of 0°∼90°. When these specimens are impacted to fail, the inverted Vshaped fractures appear at the bottom of the hole.
 
Note: “e1, e2” denotes an obvious rise process as shown in its straintime curves. The symbol “—” denotes a relatively low fluctuation of the straintime curves, generally less than 0.2. 
4. Discussion
4.1. Crack Coalescence Modes
Figure 8 shows the crack coalescence modes of precracked specimens with different flaw inclination angles. There are three types of crack coalescence modes in the present study: regional failure (Figure 8(a)), monotype fracture (Figures 8(c)–8(d)), and mixed mode fractures (Figures 8(b)–8(e)). The detailed coalescence categories are described as follows:(i)Regional failure: from Figure 8(a), fractures derived from the bottom of horizontal crack are all connected with the fissures initiating from crown area of the circular hole. Then, the damage zone is formed as an inverted triangle, and the damage intensity in this area is relatively high.(ii)Monotype fracture: in Figures 8(c) and 8(d), a shear fracture of long and narrow shape is formed and it connects the lower crack tip and the crown area of hole.(iii)Mixed mode fractures: in Figure 8(b), the tensile fracture derived from the lower tip of crack is connected with the shear fracture above the hole. Furthermore, in Figure 8(e), tensile fractures generated at both tips of crack coalesce with the shear fractures emerged at the top and bottom of the hole, respectively, forming a symmetrical failure.
(a)
(b)
(c)
(d)
(e)
It can be seen that tensile fractures are usually occurred from the crack tip, while shear fractures are often formed at the boundary of circular hole. However, this phenomenon is not applicable in the situation that the crack tip is relatively close to the hole. In this condition, the microfissures which initiate from the crack tip and crown area of the hole propagate toward each other rapidly, and then a visible fracture is formed afterwards (Figures 8(c) and 8(d)).
4.2. Average Strain Evolution of Failure Zone
Based on the above analysis about the crack coalescence modes between the combined flaws, the main failure zones which involve combined flaws are selected to analyze the average strain in specimens. As shown in Figure 9, taking the specimen with a flaw angle of 45° as an example, a circle with the center of point C and a radius of 12 mm is delineated. The average strain (only e1 is analyzed) of the area is calculated by DIC technique, and the straintime curves of five typical specimens are plotted in Figure 10.
It can be seen from Figure 10 that the average strain of the failure zone in specimen U903 (symmetrical mixed fractures along the loading direction) is most prominent. Then, specimen U03 (regional failure) follows. In specimens U301, U453, and U601, there is a dominant fracture connecting the prefabricated crack with the hole. The order of the average strain within the three failure zones is 45° > 30° > 60°. The characteristics of this trend are consistent with the crack coalescence modes in Figure 8.
4.3. Energy Characteristics
Rock failure is actually a process of energy dissipation. The specimens applied prestatic loading are in the compressive state with considerable elastic energy accumulated. Once it is subjected to sufficiently large dynamic disturbances, the stored elastic energy is rapidly released and may result in severe damage. Therefore, it is of great importance to analyze the law of energy consumption. Based on the onedimensional stress wave theory and data recorded in the SHPB tests, the average value of incident energy E_{I}, the reflected energy E_{R}, the transmitted energy E_{T}, and the absorption energy E_{A} can be calculated [48] by the following equations. The results are listed in Table 7:where , and are the incident, reflected, and transmitted stress at time t, respectively, and is the wave impedance of the bars.

At different flaw inclination angles, the incident energy E_{I}, the reflected energy E_{R}, the transmitted energy E_{T}, and the absorption energy E_{A} of the granite specimens have different values. Particularly, E_{A} is the energy absorbed by the rock specimens, and this part of energy has the main contribution for the failure of rock and also the new fractures in the specimen. Figure 11 plots the relationship between the absorption energy E_{A} and the flaw inclination angle. It is obvious that E_{A} first decreases then rises gradually with the inclination angle increasing, and it reaches the lowest value at 45°. This is because the precracked specimens with a flaw of 45° absorb only a relatively small amount of energy, fractures generate inside the specimens, and failure occurs, while other types of specimens require more energy to break. The variation trend of absorption energy is similar to that of combined strength of the specimens, as shown in Figure 6.
4.4. Rock Fragmentation Size Analysis
Macrofractures are formed in the rock specimens under coupled staticdynamic loading. This behavior inevitably leads to the fragmentation and energy dissipation in the rock failure process. The distribution characteristics of rock fragmentation size can reflect the extent of deformation and failure behavior of rocks. The greater the damage degree of rock specimen will result in better dispersion and more homogenous diameter of rock fragment. Figure 12(a) shows the fragment of the granite specimen after impact load. It is worth noting that in order to make the results of size analysis more reasonable, all fragments of the specimens in the same inclination angle are mixed together. Then, the split desktop, a rock fragmentation size measurement software, is used to analyze the size data in each group. The processed picture of Group U30 is shown in Figure 12(b), and all rock fragmentation size data are summarized in Table 8. It lists the percentage of rock fragment in each interval, as well as the average size and top size in each group. Almost all the fragment size is less than 5.08 cm, except for some fragment in Group U90, which accounts for 5.87%. Furthermore, it is clear that with the increase of inclination angle, the average size and top size of the specimens show a trend of decreasing firstly and then increasing, and they both reach the lowest at the angle of 45°, which is consistent with the variation trend of combined strength and absorption energy.
(a)
(b)

Figure 13 shows the fragment size distribution of granite specimens with different angles after crush. It reflects the fragmentation degree of the specimens. It can be seen that when the crack inclination angle is 45°, the specimens are broken into smaller particle sizes with homogenous size distribution. The specimens with an angle of 30° rank second; nevertheless, the amount of rock fragment with the size less than 0.5 cm is the maximum at this angle. The specimens with flaw angles of 0° and 60° have less small size fragments but more large size fragments. The uniformity of fragment size distribution for the specimen with 90° flaw angle is the worst.
5. Conclusions
A series of impact tests are conducted on granite specimens with combined hole and crack flaws. The DIC technique is used to monitor the surface strain evolution and failure behavior of the specimens. The conclusions are drawn as follows:(1)Mechanical properties of granite with crackhole flaws are susceptible to the crack inclination angles. The combined strength and the peak strain of the specimens both decrease first and then increase with the flaw inclination angle increasing, and both of them reach the lowest value at the angle of 45°.(2)The influence of flaw inclination angle on the deformation of the specimen reflects the fracturing behavior. Two crack tips (A, B), and the top and bottom of the internal hole (D, E) are the main stress concentration zones, and their local maximum strain is basically more than 0.5 at postpeak stage.(3)The crack coalescence modes between combined flaws are summarized, i.e., regional failure, monotype fracture, and mixed mode fractures. Tensile fractures tend to occur at the crack tip, while shear fractures usually emerge at the top and bottom of the hole and dominate the ultimate rock failure process under dynamic impact loading.(4)The absorption energy E_{A} first increases and then decreases with the flaw inclination angle. Furthermore, the size distribution of rock fragments is analyzed. The specimens with crack of 45° are broken relatively evenly with homogenous particle size distribution. The characteristics of absorption energy and size distribution are consistent with the variation trend of combined strength.
Data Availability
The relevant data could be found in the manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to acknowledge the financial support from the Outstanding Youth Science Foundations of Hunan Province of China (no. 2019JJ20028), the Key Projects of the National Natural Science Foundation of China (no. 41630642), and the Graduate Student Research and Innovation Fund of Central South University (no. 1053320183421).
References
 E. Hoek and C. D. Martin, “Fracture initiation and propagation in intact rocka review,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 6, no. 4, pp. 287–300, 2014. View at: Publisher Site  Google Scholar
 F. Feng, X. B. Li, J. Rostami, D. X. Peng, D. Y. Li, and K. Du, “Numerical investigation of hard rock strength and fracturing under polyaxial compression based on MogiCoulomb failure criterion,” International Journal of Geomechanics, vol. 19, no. 4, Article ID 04019005, 2019. View at: Publisher Site  Google Scholar
 C. P. Zhang, P. G. Ranjith, and Z. Y. Zhang, “Experimental study of fracturing behaviour in ultralow permeability formations: a comparison between CO_{2} and water fracturing,” Engineering Fracture Mechanics, vol. 217, Article ID 106541, 2019. View at: Publisher Site  Google Scholar
 F. Dai, S. Huang, K. Xia, and Z. Tan, “Some fundamental issues in dynamic compression and tension tests of rocks using split Hopkinson pressure bar,” Rock Mechanics and Rock Engineering, vol. 43, no. 6, pp. 657–666, 2010. View at: Publisher Site  Google Scholar
 Y. L. Gui, J. L. Shang, J. J. Ma, and Z. Y. Zhao, “Numerical study of the circular opening effect on mechanical behaviour of rock under confinement,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 11, no. 6, pp. 1201–1210, 2019. View at: Publisher Site  Google Scholar
 X.R. Liu, S.Q. Yang, Y.H. Huang, and J.L. Cheng, “Experimental study on the strength and fracture mechanism of sandstone containing elliptical holes and fissures under uniaxial compression,” Engineering Fracture Mechanics, vol. 205, pp. 205–217, 2019. View at: Publisher Site  Google Scholar
 W. Zeng, S.Q. Yang, and W.L. Tian, “Experimental and numerical investigation of brittle sandstone specimens containing different shapes of holes under uniaxial compression,” Engineering Fracture Mechanics, vol. 200, pp. 430–450, 2018. View at: Publisher Site  Google Scholar
 Q. Zhang, X. Wang, L. G. Tian, and D. M. Huang, “Analysis of mechanical and acoustic emission characteristics of rock materials with doublehole defects based on particle flow code,” Shock and Vibration, Article ID 7065029, 11 pages, 2018. View at: Publisher Site  Google Scholar
 Q.S. Bai, S.H. Tu, and C. Zhang, “DEM investigation of the fracture mechanism of rock disc containing hole(s) and its influence on tensile strength,” Theoretical and Applied Fracture Mechanics, vol. 86, pp. 197–216, 2016. View at: Publisher Site  Google Scholar
 Y. Liu, F. Dai, L. Dong, N. Xu, and P. Feng, “Experimental investigation on the fatigue mechanical properties of intermittently jointed rock models under cyclic uniaxial compression with different loading parameters,” Rock Mechanics and Rock Engineering, vol. 51, no. 1, pp. 47–68, 2018. View at: Publisher Site  Google Scholar
 Y. Gao, X.T. Feng, Z. Wang, and X. Zhang, “Strength and failure characteristics of jointed marble under true triaxial compression,” Bulletin of Engineering Geology and the Environment, vol. 79, no. 2, pp. 891–905, 2020. View at: Publisher Site  Google Scholar
 S.Q. Yang, X.R. Liu, and H.W. Jing, “Experimental investigation on fracture coalescence behavior of red sandstone containing two unparallel fissures under uniaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 63, pp. 82–92, 2013. View at: Publisher Site  Google Scholar
 L. N. Y. Wong and H. H. Einstein, “Systematic evaluation of cracking behavior in specimens containing single flaws under uniaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 46, no. 2, pp. 239–249, 2009. View at: Publisher Site  Google Scholar
 D. Huang, D. Gu, C. Yang, R. Huang, and G. Fu, “Investigation on mechanical behaviors of sandstone with two preexisting flaws under triaxial compression,” Rock Mechanics and Rock Engineering, vol. 49, no. 2, pp. 375–399, 2016. View at: Publisher Site  Google Scholar
 J. Wu, M. Feng, B. Yu, and G. Han, “The length of preexisting fissures effects on the mechanical properties of cracked red sandstone and strength design in engineering,” Ultrasonics, vol. 82, pp. 188–199, 2018. View at: Publisher Site  Google Scholar
 S.q. Yang, H.w. Jing, and T. Xu, “Mechanical behavior and failure analysis of brittle sandstone specimens containing combined flaws under uniaxial compression,” Journal of Central South University, vol. 21, no. 5, pp. 2059–2073, 2014. View at: Publisher Site  Google Scholar
 T. Liu, B. Lin, W. Yang, Q. Zou, J. Kong, and F. Yan, “Cracking process and stress field evolution in specimen containing combined flaw under uniaxial compression,” Rock Mechanics and Rock Engineering, vol. 49, no. 8, pp. 3095–3113, 2016. View at: Publisher Site  Google Scholar
 Q. Yin, H. Jing, and H. Su, “Investigation on mechanical behavior and crack coalescence of sandstone specimens containing fissurehole combined flaws under uniaxial compression,” Geosciences Journal, vol. 22, no. 5, pp. 825–842, 2018. View at: Publisher Site  Google Scholar
 C. Zou and L. N. Y. Wong, “Experimental studies on cracking processes and failure in marble under dynamic loading,” Engineering Geology, vol. 173, pp. 19–23, 2014. View at: Google Scholar
 C. Zou, L. N. Y. Wong, J. J. Loo, and B. S. Gan, “Different mechanical and cracking behaviors of singleflawed brittle gypsum specimens under dynamic and quasistatic loadings,” Engineering Geology, vol. 201, pp. 71–84, 2016. View at: Publisher Site  Google Scholar
 X. B. Li, T. Zhou, and D. Y. Li, “Dynamic strength and fracturing behavior of singleflawed prismatic marble specimens under impact loading with a splithopkinson pressure bar,” Rock Mechanics and Rock Engineering, vol. 50, no. 1, pp. 29–44, 2017. View at: Publisher Site  Google Scholar
 Z. Jing, W. Dong, and W. L. Chang, Experimental Study on Postpeak Mechanical Properties of Hard Rocks under Combined Static and Dynamic Loading, Sichuan University, Chengdu, China, 2015.
 P. Feng, F. Dai, Y. Liu, and H.b. Du, “Mechanical behaviors of rocklike specimens with two noncoplanar fissures subjected to coupled staticdynamic loads,” Engineering Fracture Mechanics, vol. 199, pp. 692–704, 2018. View at: Publisher Site  Google Scholar
 W. C. Zhu, Y. Bai, X. B. Li, and L. L. Niu, “Numerical simulation on rock failure under combined static and dynamic loading during SHPB tests,” International Journal of Impact Engineering, vol. 49, pp. 142–157, 2012. View at: Publisher Site  Google Scholar
 K. Bharadwaj, A. Sheidaei, A. Afshar, and J. Baqersad, “Fullfield strain prediction using mode shapes measured with digital image correlation,” Measurement, vol. 139, pp. 326–333, 2019. View at: Publisher Site  Google Scholar
 J. Baqersad, J. Carr, T. Lundstrom, C. Niezrecki, P. Avitabile, and M. Slattery, “Dynamic characteristics of a wind turbine blade using 3D digital image correlation,” International Society Optics and Photonics, vol. 8348, pp. 83482I–83489I, 2012. View at: Google Scholar
 D. Li, L. N. Y. Wong, G. Liu, and X. Zhang, “Influence of water content and anisotropy on the strength and deformability of low porosity metasedimentary rocks under triaxial compression,” Engineering Geology, vol. 126, pp. 46–66, 2012. View at: Publisher Site  Google Scholar
 K. Muniandy, Z. Mohamad Ariff, and A. Abu Bakar, “Digital image correlation utilization in measuring displacement and strain during plastic film blowing process: a feasibility study,” Measurement, vol. 136, pp. 487–500, 2019. View at: Publisher Site  Google Scholar
 G. L. Golewski, “Measurement of fracture mechanics parameters of concrete containing fly ash thanks to use of Digital Image Correlation (DIC) method,” Measurement, vol. 135, pp. 96–105, 2019. View at: Publisher Site  Google Scholar
 W. W. Ji, P. Z. Pan, S. T. Miao, F. S. Su, and M. P. Du, “Fracture characteristics of two types of rocks based on digital image correlation,” Rock and Soil Mechanics, vol. 37, no. 8, pp. 2299–2305, 2016. View at: Google Scholar
 X. Li, Z. Zhou, T.S. Lok, L. Hong, and T. Yin, “Innovative testing technique of rock subjected to coupled static and dynamic loads,” International Journal of Rock Mechanics and Mining Sciences, vol. 45, no. 5, pp. 739–748, 2008. View at: Publisher Site  Google Scholar
 W.W. Ji, P.Z. Pan, Q. Lin, X.T. Feng, and M.P. Du, “Do disktype specimens generate a mode II fracture without confinement?” International Journal of Rock Mechanics and Mining Sciences, vol. 87, pp. 48–54, 2016. View at: Publisher Site  Google Scholar
 D. Y. Li, Z. Y. Han, X. L. Sun, T. Zhou, and X. B. Li, “Dynamic mechanical properties and fracturing behavior of marble specimens containing single and double flaws in SHPB tests,” Rock Mechanics Rock Engineering, vol. 52, no. 6, pp. 1623–2164, 2019. View at: Google Scholar
 H. W. Zhou, Z. H. Wang, C. S. Wang, and J. F. Liu, “On acoustic emission and postpeak energy evolution in beishan granite under cyclic loading,” Rock Mechanics and Rock Engineering, vol. 52, no. 1, pp. 283–288, 2019. View at: Publisher Site  Google Scholar
 W. G. P. Kumari, P. G. Ranjith, M. S. A. Perera et al., “Hydraulic fracturing under high temperature and pressure conditions with micro CT applications: geothermal energy from hot dry rocks,” Fuel, vol. 230, pp. 138–154, 2018. View at: Publisher Site  Google Scholar
 X. P. Zhou, Y. X. Zhang, and Q. L. Ha, “Realtime computerized tomography (CT) experiments on limestone damage evolution during unloading,” Theoretical and Applied Fracture Mechanics, vol. 50, no. 1, pp. 49–56, 2008. View at: Publisher Site  Google Scholar
 M. M. Kou, X. R. Liu, S. D. Tang, and Y. T. Wang, “3D Xray computed tomography on failure characteristics of rocklike materials under coupled hydromechanical loading,” Theoretical and Applied Fracture Mechanics, vol. 104, Article ID 102396, 2019. View at: Publisher Site  Google Scholar
 G. Li, X. F. Cheng, H. Pu, and C. A. Tang, “Damage smear method for rock failure process analysis,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 11, no. 6, pp. 1151–1165, 2019. View at: Publisher Site  Google Scholar
 Y. Q. Zhou, Q. Sheng, N. N. Li, and X. D. Fu, “Numerical investigation of the deformation properties of rock materials subjected to cyclic compression by the finite element method,” Soil Dynamics and Earthquake Engineering, vol. 126, Article ID 105795, 2019. View at: Publisher Site  Google Scholar
 Y. Ju, Y. L. Wang, C. S. Su, D. S. Zhang, and Z. Y. Ren, “Numerical analysis of the dynamic evolution of mininginduced stresses and fractures in multilayered rock strata using continuumbased discrete element methods,” International Journal of Rock Mechanics and Mining Sciences, vol. 113, no. 1, pp. 91–210, 2019. View at: Publisher Site  Google Scholar
 C. A. Tang and S. Tang, “Applications of rock failure process analysis (RFPA) method,” Journal of Rock Mechanics and Geotechnical Engineering, vol. 3, no. 4, pp. 352–372, 2011. View at: Publisher Site  Google Scholar
 W. Zhao, R. Huang, and M. Yan, “Study on the deformation and failure modes of rock mass containing concentrated parallel joints with different spacing and number based on smooth joint model in PFC,” Arabian Journal of Geosciences, vol. 8, no. 10, pp. 7887–7897, 2015. View at: Publisher Site  Google Scholar
 M. Tao, A. Ma, W. Cao, X. Li, and F. Gong, “Dynamic response of prestressed rock with a circular cavity subject to transient loading,” International Journal of Rock Mechanics and Mining Sciences, vol. 99, pp. 1–8, 2017. View at: Publisher Site  Google Scholar
 D. Y. Li, P. Xiao, Z. Y. Han, and Q. Q. Zhu, “Mechanical and failure properties of rocks with a cavity under coupled static and dynamic loads,” Engineering Fracture Mechanics, vol. 225, Article ID 106195, 2020. View at: Publisher Site  Google Scholar
 Q. Zhu, D. Li, Z. Han, X. Li, and Z. Zhou, “Mechanical properties and fracture evolution of sandstone specimens containing different inclusions under uniaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 115, pp. 33–47, 2019. View at: Publisher Site  Google Scholar
 Correlated Solutions Vic2D User’s Manual, 4th edition, SC, West Columbia, 2006.
 B. Pan, A. Asundi, H. M. Xie, and J. X. Gao, “Digital image correlation using iterative least squares and pointwise least squares for displacement field and strain field measurements,” Optics and Lasers in Engineering, vol. 47, no. 78, pp. 865–874, 2009. View at: Publisher Site  Google Scholar
 Z. Y. Han, D. Y. Li, T. Zhou, Q. Q. Zhu, and P. G. Ranjith, “Experimental study of stress wave propagation and energy characteristics across rock specimens containing cemented mortar joint with various thicknesses,” International Journal of Rock Mechanics and Mining Sciences, vol. 131, Article ID 104352, 2020. View at: Publisher Site  Google Scholar
 C. Zou, L. N. Y. Wong, J. J. Loo, and B. S. Gan, “Different mechanical and cracking behaviors of singleflawed brittle gypsum specimens under dynamic and quasistatic loadings,” Engineering Geology, vol. 201, pp. 71–84, 2016. View at: Publisher Site  Google Scholar
 Q. B. Zhang and J. Zhao, “Determination of mechanical properties and fullfield strain measurements of rock material under dynamic loads,” International Journal of Rock Mechanics and Mining Sciences, vol. 60, no. 8, pp. 423–439, 2013. View at: Publisher Site  Google Scholar
 S. P. Morgan, C. A. Johnson, and H. H. Einstein, “Cracking processes in Barre granite: fracture process zones and crack coalescence,” International Journal of Fracture, vol. 180, no. 2, pp. 177–204, 2013. View at: Publisher Site  Google Scholar
 J. J. F. Bonnen and T. H. Topper, “Threshold crack growth behavior of shear and tensile cracks,” International Journal of Fatigue, vol. 42, pp. 122–130, 2012. View at: Publisher Site  Google Scholar
 M. G. Tarantino, S. Beretta, S. Foletti, and I. Papadopoulos, “Experiments under pure shear and rolling contact fatigue conditions: competition between tensile and shear mode crack growth,” International Journal of Fatigue, vol. 46, pp. 67–80, 2013. View at: Publisher Site  Google Scholar
 S. Liu, Y. J. Chao, and X. Zhu, “Tensileshear transition in mixed mode I/III fracture,” International Journal of Solids and Structures, vol. 41, no. 2223, pp. 6147–6172, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Diyuan Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.