Shock and Vibration
Volume 2018 (2018), Article ID 3105384, 13 pages
https://doi.org/10.1155/2018/3105384

Research Article

Rock Dynamic Fracture Characteristics Based on NSCB Impact Method

Yanbing Wang1,2,3

¹School of Mechanics and Architecture Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China
²State Key Laboratory for Geomechanics and Deep Underground Engineering, Beijing 100083, China
³Department of Mechanical Engineering, University of South Carolina, Columbia, SC 29208, USA

Correspondence should be addressed to Yanbing Wang; ceowyb818@163.com

Received 2 December 2017; Revised 12 March 2018; Accepted 20 March 2018; Published 21 May 2018

Academic Editor: Sergio De Rosa

Copyright © 2018 Yanbing Wang. 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.

Abstract.

In 2012, the International Society for Rock Mechanics (ISRM) recommended a new Notched Semicircular Bend (NSCB) method for the determination of dynamic fracture toughness of rock materials, but it did not consider the effect of some uncontrollable factors in the course of the experiment on the test result. This thesis firstly carried out dynamic fracture toughness experiments on several typical rock materials such as sandstone using the modified Split-Hopkinson Pressure Bar (SHPB) experimental system with high-speed photography, directly compared the dynamic fracture failure characteristics of several rock materials, and examined the loading rate dependency of the dynamic fracture toughness of rock materials. Based on the numerical analysis method of Discrete Lattice Spring Model (DLSM), it focused on the effect of bullet impact loading rate, loading area of incident bar, support restraints of clamping specimen, and other uncontrollable factors in the course of SHPB experiment on test results. The findings can be referenced for the improvement of NSCB method.

1. Introduction

In rock engineering practice, the action of dynamic loading is a commonplace, including engineering blasting, earthquake, or rock burst [1–4]. For the engineering safety and the maximization of economic benefits, the mechanical properties of rock such as dynamic fracture toughness have to be accurate in the design and construction of these projects [5, 6]. As the dynamic experiment is more complex than static experiment with respect to the rock materials, in 2012, the Committee of Rock Dynamics of International Society for Rock Mechanics (ISRM-CRD) first recommended some methods for rock dynamic experiment, including Notched Semicircular Bend (NSCB) for the testing on dynamic fracture toughness of rock materials [7]. Wei et al. [8] estimated theoretically the lengths of fracture process zones (FPZs) in the NSCB and CCNSCB tests with representative specimen geometries using the maximum normal stress criterion, and their effects on these fracture toughness tests were evaluated according to the effective crack model. Wang and Yang [9] tested the dynamic fracture toughness of coal samples based on the Notched Semicircular Bend (NSCB) method under an impact load and discussed the influence of the impact velocity and bedding angle on the dynamic fracture toughness, failure strain, and strain rate of coal samples. Yao et al. [10] studied the effect of thermally induced damage on the dynamic fracture toughness of two mortars. The Notched Semicircular Bend (NSCB) specimens were heat-treated at 150 degrees C, 250 degrees C, 350 degrees C, 450 degrees C, and 600 degrees C and then tested at room temperature of 25 degrees C with a Split-Hopkinson Pressure Bar (SHPB) apparatus. Chen et al. [11] selected the NSCB method to investigate the fracture initiation toughness of rocks subjected to preload. Wei et al. [12] determined the critical SIFs for a wide range of CCNSCB (cracked chevron notched semicircular bending method) specimen geometries via three-dimensional finite element analysis. Gao et al. [13] measured more fracture parameters of the NSCB specimen loaded by a SHPB system using the digital image correlation (DIC) method combined with ultra-high-speed photography. Dai et al. [14] studied the micromechanical model for the rate dependence of the fracture toughness anisotropy of Barre granite based on NSCB method. Yin et al. [15] studied the effect of thermal treatment on the dynamic fracture toughness of Laurentian granite using NSCB method. This NSCB method is a test unit based on the modified Split-Hopkinson Pressure Bar (SHPB), but it does not consider the effect of bullet impact loading rate, loading area of incident bar, support restraints of clamping specimen, and other uncontrollable factors in the course of SHPB experiment on test results.

Based on the ISRM-recommended NSCB, the dynamic fracture toughness of several typical rock materials including sandstone was determined here, and the effect of uncontrollable factors in the course of the experiment on test results was analyzed using Discrete Lattice Spring Model (DLSM).

2. Test Unit and Principle

The Φ50 mm Split-Hopkinson Pressure Bar (SHPB) experimental system from State Key Laboratory for Geomechanics & Deep Underground Engineering was used. It comprises pressure control system, bullet, input bar, output bar, absorbing bar, damper, laser speed measurement system, measurement recording system (ultra-dynamic strain gauge, waveform memory, and data processing system), and high-speed camera, as shown in Figure 1. This experiment was based on the following assumptions: (1) one-dimensional stress wave and (2) uniformity.

Figure 1: Diagram of SHPB experimental system.

The specimen was clamped between the input and output bar, the incident wave was generated when the bullet impacted the input bar, and the action between the incident wave and the specimen generated the reflected wave and the transmitted wave , where the waveforms will be recorded by strain gauges attached on the input and output bar. In the course of the conventional SHPB experiment, the majority of incident waves were rectangular with high-frequency oscillations, which indicates that the dynamic force did not reach equilibrium in early period of the loading. The contact force between the input bar and the specimen and the contact force between the output bar and the specimen were, respectively, as follows:

where is the elastic modulus of bar material and is the cross-sectional area of the bar.

To make the dynamic loading on the ends of the specimen at equilibrium, 1 mm thick round rubber gasket was used as the pulse shaper to be attached to the end of the incident bar, resulting in a half-sinusoidal incident wave, as shown in Figure 2. The specimen was kept at the state of dynamic equilibrium before the failure, that is, ≈ (or ), which is the key to the success of dynamic loading experiment. The calibration of stress time-history curves at the ends of the specimen can be analyzed through comparison.

Figure 2: Output signals.

The experiment procedure was based on [16], and the dynamic fracture process of the rock specimen was captured with a high-speed camera. The test data was processed once the experiment was completed.

The dimension of NSCB specimen is indicated in Figure 3; according to the ISRM-recommended dimension [7], the diameter , the thickness , and the notch length for each specimen were 50, 25, and 4 mm, respectively, and the distance between two supports was 30 mm. Thus, the dimensionless quantity , , and .

Figure 3: Dimension of NSCB specimen.

The dynamic loading forces at the ends of specimen must be kept in balance. Figure 4 shows a typical dynamic balance check curve, and from (1), the dynamic loading force at the ends of the specimen was determined to be balanced.

Figure 4: Dynamic force balance check.

The type stress intensity factor (SIF), , can be expressed as [7]

where is the radius of a specimen and is a dimensionless quantity, and its value is related to the geometric dimension of the notch; for and in this experiment, can be expressed as

The peak point of curve was used as the dynamic fracture toughness, , which has a significant loading rate dependency [7], and the slope of the straight segment before the peak of SIF time-history curve was defined as the loading rate, and the loading rate as shown in Figure 5 was 90 GPa·m1/2/s.

Figure 5: SIF time-history curve.

3. Dynamic Fracture Toughness Test for Several Typical Rocks

3.1. SHPB-Based Dynamic Fracture Roughness Test

Five typical rock materials were selected and the basic physical and mechanical parameters were given in Table 1.

Table 1: Physical and mechanical parameters of typical rocks.

Density (kg/m³)
Elastic modulus (GPa)
Poisson’s ratio
Uniaxial compressive strength (MPa)
Place of Origin

Sandstone
2150
4.6
0.24
41
Immediate roof of concentrated track tunnel, the Second Mining Area, Shuguang Mine of Shanxi Coking Coal

Mudstone
2480
6.5
0.25
28
Immediate floor of concentrated track tunnel, the Second Mining Area, Shuguang Mine of Shanxi Coking Coal

Granite
2730
18.4
0.2
103
Tongtang Town, Tong’an District, Xiamen

Coal
1475
7.9
0.23
18
The Second Mining Area, Shuguang Mine of Shanxi Coking Coal

PMMA
1127
6.1
0.28

Beijing Xuanguang Acrylic Factory

3.2. Test Results and Analysis

Fastcam-SA5(16G) color high-speed digital camera produced by Japan Photron was used in the experiment. Figure 6 shows the process of rock specimen failure recorded by the high-speed camera. In view of the low gloss of the rock specimen, we have increased the brightness and contrast of some images by 40% and 40%, respectively, for example, coal and PMMA. The cracks initiated from the notch and ran through the entire specimen within 10 μs, and the semicircular rock specimen was symmetrically split into two equal parts. The mineral particles in the sandstone were relatively smaller, and the trajectory of its crack propagation was relatively straight to almost split the specimen into two equal parts; the mudstone had good integrity and the trajectory of its crack propagation was bent a little; the mineral composition in granite had a relatively large volume, and the cracks initiated from the notch and then propagated to the loading point along the boundary of minerals; its trajectory was distorted, with an obvious feature of intergranular fracture; the looseness of coal was relatively high, the internal joint fissure was fully developed, and the trajectory of its crack propagation was bent; PMMA is one of the transparent media; the trajectory of its crack propagation was straight but showed features of shear failure along the thickness direction.

(a) Sandstone

(b) Mudstone

(c) Granite

(d) Coal

(e) PMMA

Figure 6: Process of rock specimen failure.

Figure 7 shows the calculated results for dynamic fracture toughness. By comparing the five rock specimens in the experiment, given the same loading rate, the fracture toughness of granite was the highest, while that of coal was the lowest; sandstone and mudstone had similar results, and that of PMMA was lower than any of these samples, but higher than coal. The dynamic fracture toughness for five rock materials shows a significant loading rate dependency; as the loading rate increases, the dynamic fracture toughness of rock rises gradually. The fracture toughness and loading rate were fit into the linear equation

where is the type of rock specimen and and are the linear correlation coefficients. Table 2 shows the fitting results.

Table 2: Linear correlation coefficients for the fitting curve of fracture toughness-loading rate.

Rock Type
(s)
Standard Deviation
(Mpa⋅m1/2)
Standard Deviation

Sandstone
0.04049
0.00319
0.37337
0.29452

Mudstone
0.03785
0.00156
0.51833
0.13249

Granite
0.02928
0.00506
2.19798
0.45141

Coal
0.01765
0.0026
0.28858
0.22856

PMMA
0.01566
0.00103
1.71827
0.10356

Figure 7: Test results of the dynamic fracture toughness for several typical rocks.

From Table 2, the slope of the fitting curve for fracture toughness of sandstone, , was the biggest, but the difference was relatively smaller, which suggests that the change of the fracture toughness of sandstone was the quickest with the increase of loading rate; mudstone and sandstone showed similar characteristics; the difference of fracture toughness fitting curve for granite was the greatest, and the standard deviation for the slope and the difference of its fitting curve were the maximum among all rock specimens. And from Figure 7, when the loading rate was low, the increase of fracture toughness was relatively fast; when the loading rate was high, the curve was relatively stable and the increase of fracture toughness was slow, but its linear characteristics were not obvious, which is closely related to the mineral composition of granite; the difference of the fitting curve for fracture toughness of coal, , was the smallest, which is related to the characteristics of coal itself; the slope of the fitting curve for fracture toughness of PMMA, , was the smallest, but the difference was relatively higher, which suggests that the fracture toughness of such material was less sensitive to the loading rate, and its fracture toughness would increase linearly and gradually with the loading rate.

4. DLSM Numeral Simulation-Based Analysis of the Dynamic Fracture Characteristics for Typical Rocks

In order to examine the dynamic fracture toughness for different rock types, it is necessary to conduct a three-point bending experiment given the same impact rate. However, the SHPB used here was in pressure control; it would be too difficult to control the bullet impact rate at a given value. Furthermore, the total factor test would add to the cost of trials. To address this, DLSM numerical simulation with physical experiment was an alternative approach.

4.1. Discrete Lattice Spring Model (DLSM)

In the Discrete Lattice Spring Model (DLSM), the solid medium was represented by a set of particles and connected with a spring [17], and its mechanical responses were demonstrated by the appropriate lattice structure. When the gap among particles was less than the given threshold, the adjacent particles would be connected by the dummy “spring ligament”, as shown in Figure 8. Different assembly forms of particles and threshold values would result in different structures of lattice spring. A “ligament” comprises a normal spring and a set of multibody tangential springs; DLSM can simulate the short- and long-distance interaction among the particles through a ligament, as shown in Figure 8(a). The specific details about DLSM and the associated validation and application are stated in [17–19].

(a) Normal, tangential spring

(b) Forces applied on the particle

Figure 8: Particles and spring ligament.

4.2. Establishment of the Test Model for Rock Fracture Toughness

Figure 9 shows the DLSM and its size is the same as the actual specimen. The diameter of particle is 0.25 mm and the whole model is composed of 300,000 particles. The physical and mechanical parameters selected here are the actual parameters of rock, as listed in Table 1, with a spring damage parameter of 0.005. The dynamic loads in the course of the experiment were characterized by the rate loading applied above the middle part of the specimen, and its loading parameters are

where is the applied dynamic rate loading, in mm/s; is the required time for such rate changing from 0 to the given value, which was selected as 20 μs for all models. At the initial stage after the loads were applied, the specimen would easily reach a dynamic equilibrium by increasing the amplitude of loads. Such time of increase was usually the time taken for 5 round travels of the stress wave within the specimen, and the loading rate induced by different loading rates may be varied.

Figure 9: Model of NSCB specimen.

4.3. Numerical Calculation Results

Take sandstone as an example; Figure 10 shows the propagation of stress wave from the impact load in the sandstone specimen and the dynamic fracture process of the semicircular three-point bending sample when  m/s; the stress wave propagated from the impact loading point to the lower boundary and reached the notch at 36.25 μs, when the stress concentration at the notch showed obvious features; at 73.75 μs, the crack initiated from the tip of the notch. In the DLSM simulation, the propagation morphology of cracks and their displacement characteristics are similar to the actual crack characteristics as shown in the high-speed captures of Figure 6, which also makes up for the deficiency that the process of crack initiation and propagation is incomplete due to low configuration of the high-speed camera. Figure 11 gives the numerical results for mudstone, granite, coal, and PMMA, respectively.

Figure 10: Simulation of the rock specimen fracture (sandstone).

(a) Mudstone

(b) Granite

(c) Coal

(d) PMMA

Figure 11: DLSM calculated results for the other four rock specimens.

Figure 12 shows the curve of the position of crack tip with time, and the crack propagation can be roughly divided into two stages, namely, the stable propagation at the early stage of crack propagation and the unstable propagation at its late stage. The slope of curve at the stable propagation stage was used as the average velocity for the crack propagation; the majority of average velocities for stable crack propagation in numerical calculation were 350 m/s–400 m/s, while that of the cracks in coal was about 250 m/s. The occurrence of unstable propagation stage is attributed to the fact that a large deformation appeared on the specimen and the damage of particles was not regular anymore at the late stage of numerical calculation.

Figure 12: Displacement-time curve for the crack propagation,  m/s.

Figure 13 shows the loading force-displacement curves which are typical unimodal smooth. The peak point in each curve was used as the value for , which was substituted into (2) for the calculation of dynamic fracture toughness. Figure 14 shows when  m/s, 6 m/s, 8 m/s, 12 m/s, 14 m/s, 16 m/s, 18 m/s, or 20 m/s, the dynamic fracture toughness corresponding to each type of rock. When the impact velocity exceeded a certain value, the damage on rock specimen firstly occurred at the part where the load was applied, rather than at the tip of the notch. For example, Figure 15 shows that when  m/s, the loading point of coal specimen was firstly damaged; when  m/s, the sandstone and mudstone specimens also showed a similar case. This suggests that, with such impact velocity, the dynamic semicircular three-point bending experiment will be unsuitable to obtain the dynamic fracture toughness of rock. Numerical results show that the dynamic fracture toughness of rock specimen increased linearly with the loading rate; given the same impact rate, the dynamic fracture toughness of granite was the greatest, while it was the smallest for coal; that of the sandstone and mudstone was between the values and their fracture toughness was similar; the fracture toughness of PMMA was greater than coal but less than rock.

Figure 13: Loading force-displacement curve,  m/s.

Figure 14: DLSM calculated results of the fracture toughness.

(a) Coal:  m/s

(b) Granite:  m/s

Figure 15: Damage morphology of partial rock specimens.

4.3.1. Effect of the Loading Area

In the course of SHPB experiment, the part where the input bar impacted the specimen can be considered as the linear loading, namely, the top of the upper boundary of the arc contacted with the input bar. However, in the course of the actual experiment, when the input bar was applied to the specimen, the upper boundary of the arc was firstly under pressure, which means that the dynamic loading was actually applied to a very small plane, but the exact area of such plane remains to be identified.

With respect to the problems stated above, take the sandstone as an example; the DLSM for each of the four different loading areas was established as shown in Figure 16, where (a) is the free face loading, namely, the above simulations used this method; (b) is the approximately linear loading, in which the loading face is only 1 mm long; (c) surface loading, the loading surface was 2 mm long; (d) surface loading, the loading surface was 5 mm long;  m/s.

(a)

(b)

(c)

(d)

Figure 16: Four DLSMs for each of the four different loading areas.

Figure 17 shows the DLSM calculated results. The calculated results in case of the free face loading and the 5 mm surface loading are in good consistency with the experiment; the specimen was symmetrically split into two parts; only the particles near the left support were damaged when the splitting process was about to end, which is consistent with the reality. At the late stage of SHPB experiment, the input bar was continuously pressed to the broken specimen at a high speed, and the support became the point of stress concentration, which will be absolutely damaged. In case of the 2 mm surface loading, the fracture began from the notch, but as the crack propagated, the particles at the loading point were also damaged. Because the basic specimen had been damaged, due to the relatively small area of the loading end, the stress concentration was still very obvious. In case of the linear loading, although the fracture began from the notch, as soon as the crack initiated, the particles at the loading point were also damaged, and as the crack propagated, the crack tip generated two branches; the crack finally became distorted, until it ended at the right end of the loading point.

(a)

(b)

(c)

(d)

Figure 17: DLSM calculated results.

From the displacement-loading curve in Figure 18, as the displacement increased, the load at the time of linear loading increased to a certain value and then entered the stage of “plastic platform”, followed by the rapid reduction to zero; when the displacement continued to increase, the load increased again, which indicates that the semicircular three-point bending dynamic fracture experiment with the linear loading will be unsuitable to obtain the dynamic fracture toughness of rock. From the curve, it also shows that, with the decrease of loading area, the time when the loading peak occurs gradually delayed, along with the gradual reduction of peak value. From the displacement-time curve in Figure 19, once the loading point at the upper boundary was damaged, the displacement curve will “jump”; the software failed to induce the curve of crack tip position automatically; for instance, in case of the 2 mm loading face, at 116 μs, the loading point had a fracture and the displacement curve would exhibit an error, but it does not affect the calculation of fracture toughness, because the stress value as soon as the crack initiated at the notch was regarded as the peak value of loading force, and then the appropriate fracture toughness can be calculated.

Figure 18: Displacement-loading curve.

Figure 19: Displacement-time curve.

4.3.2. Effect of the Support Restraints

In the course of DLSM simulation three-point bending experiment, the necessary boundary restraints should be added to the two supports, as shown in Figure 16(a) (L Support, R Support); different restraints may incur an unexpected effect on the calculated results. To examine the effect of support restraints on dynamic fracture toughness test, three different DLSMs were established for the sandstone experiment and Table 3 shows its boundary restraints, in which Scheme I is the support restraints as mentioned earlier.

Table 3: Support restraints.

No.
L Support
R Support
Remark

-direction
-direction
-direction
-direction

I
0
0
—
0
L Support was fixed, while R Support was smooth

II
0
0
0
0
Two supports were fully fixed

III
—
0
—
0
Two supports were fully smooth in -direction, while being restrained in -direction

“0” means that the velocity was zero; “—” means that there was no restraint.

Figure 20 shows the numerical results of the different restraints, in which the results of Scheme I are as shown in Figure 17(a) and they are not mentioned again. In Scheme II, although the crack initiated at the notch but did not propagate, there was a fracture at the loading point and the range of such fracture expanded continuously; the crack began to propagate from the notch in a “herringbone” shape. In Scheme III, the crack was gradually branched into two cracks to run through the specimen, after it propagated with certain distance. Since these fracture morphologies are inconsistent with the experiment, the fracture toughness cannot be calculated accordingly, which also indicates that the support restraints will have a significant effect on the calculated results.

(a) II

(b) III

Figure 20: Numerical results for different restraints.

5. Conclusion

The dynamic fracture toughness of several typical rock materials including sandstone in normal state was studied using NSCB-recommended test method and SHPB experiment system, with high-speed photography and DLSM numerical analysis. The following conclusions can be drawn.

The dynamic fracture toughness of rock materials showed an obvious loading rate dependency and the loading rate tended to be positively correlated with the dynamic fracture toughness of rock. Through this experiment, given the same loading rate, the fracture toughness of granite was the highest, while it was the lowest for coal; sandstone and mudstone had similar results, and that of PMMA was lower than any of these samples, but higher than coal.

DLSM numerical results are also similar to the test results. With the increase of the bullet impact loading rate, the dynamic fracture toughness increased linearly; when the impact velocity exceeded a certain value, the damage on rock specimen firstly occurred at the part where the load was applied, rather than at the tip of notch. In this case, the dynamic fracture toughness of rock cannot be calculated correctly. The loading face had a significant effect on the tested results; an ideal linear loading would be unsuitable to the dynamic fracture toughness test of rock materials, but the numerical results in case of the free face loading and the 5 mm surface loading were in good consistency with the experiment. The restraints at the bottom support will also affect the test. Since the numerical simulation results in case of “two supports are fully fixed” and “two supports are fully smooth in -direction, while restrained in -direction” were inconsistent with the experiment, the fracture toughness cannot be calculated accordingly.

As for the dynamic fracture toughness test of rock, the tips of prefabricated small notch should be as sharp as possible. We can try to change the shape of the end of the incident bar. Perhaps we can design an arc section.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author would like to thank Dr. Gao-feng Zhao for the help in performing the DLSM numerical simulation and the National Natural Science Foundation of China (no. 51704295) and Project funded by China Postdoctoral Science Foundation (2016M600147) for their support.

References

R.-S. Yang, Y.-B. Wang, Z.-W. Yue, and G.-Q. Liu, “Dynamic behaviors of crack propagation in directional fracture blasting with two holes,” Explosion and Shock Waves, vol. 33, no. 6, pp. 631–637, 2013.

R. Yang, Y. Wang, and C. Ding, “Laboratory study of wave propagation due to explosion in a jointed medium,” International Journal of Rock Mechanics and Mining Sciences, vol. 81, pp. 70–78, 2016.

N. Li, L.-J. Chen, and P. Zhang, “Dynamic analysis for fracturing progress by detonation gas,” Chinese Journal of Geotechnical Engineering, vol. 28, no. 4, pp. 460–463, 2006.

R. Yang, Y. Wang, L. Hou, and L. Yang, “DLDC experiment on crack propagation in defective medium under impact loading,” Chinese Journal of Rock Mechanics and Engineering, vol. 33, no. 10, pp. 1971–1976, 2014.

R.-S. Yang, Y.-B. Wang, H.-J. Xue, and K. Liu, “SEM experiment of rock crack cross section morphology after explosion fracturing with slotted cartridge,” Journal of China University of Mining and Technology, vol. 42, no. 3, pp. 337–341, 2013.

L. Yang, R. Yang, X. Zhao et al., “Experimental investigation on dynamic fractures of some metamorphic rock using reflected caustics,” Chinese Journal of Rock Mechanics and Engineering, vol. 34, no. 4, pp. 779–786, 2015.

Y. X. Zhou, K. Xia, X. B. Li et al., “Suggested methods for determining the dynamic strength parameters and mode-I fracture toughness of rock materials,” International Journal of Rock Mechanics and Mining Sciences, vol. 49, pp. 105–112, 2012.

M.-D. Wei, F. Dai, N.-W. Xu, T. Zhao, and Y. Liu, “An experimental and theoretical assessment of semi-circular bend specimens with chevron and straight-through notches for mode I fracture toughness testing of rocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 99, pp. 28–38, 2017.

Y. Wang and R. Yang, “Study of the dynamic fracture characteristics of coal with a bedding structure based on the NSCB impact test,” Engineering Fracture Mechanics, vol. 184, pp. 319–338, 2017.

W. Yao, Y. Xu, H.-W. Liu, and K. Xia, “Quantification of thermally induced damage and its effect on dynamic fracture toughness of two mortars,” Engineering Fracture Mechanics, vol. 169, pp. 74–88, 2017.

R. Chen, K. Li, K. Xia, Y. Lin, W. Yao, and F. Lu, “Dynamic Fracture Properties of Rocks Subjected to Static Pre-load Using Notched Semi-circular Bend Method,” Rock Mechanics and Rock Engineering, vol. 49, no. 10, pp. 3865–3872, 2016.

M.-D. Wei, F. Dai, N.-W. Xu, J.-F. Liu, and Y. Xu, “Experimental and numerical study on the cracked chevron notched semi-circular bend method for characterizing the mode I fracture toughness of rocks,” Rock Mechanics and Rock Engineering, vol. 49, no. 5, pp. 1595–1609, 2016.

G. Gao, S. Huang, K. Xia, and Z. Li, “Application of Digital Image Correlation (DIC) in Dynamic Notched Semi-Circular Bend (NSCB) tests,” Experimental Mechanics, vol. 55, no. 1, pp. 95–104, 2015.

F. Dai, K. Xia, and M. H. B. Nasseri, “Micromechanical model for the rate dependence of the fracture toughness anisotropy of Barre granite,” International Journal of Rock Mechanics and Mining Sciences, vol. 63, pp. 113–121, 2013.

T. Yin, X. Li, K. Xia, and S. Huang, “Effect of thermal treatment on the dynamic fracture toughness of laurentian granite,” Rock Mechanics and Rock Engineering, vol. 45, no. 6, pp. 1087–1094, 2012.

Y. Zhao, G.-F. Zhao, Y. Jiang, D. Elsworth, and Y. Huang, “Effects of bedding on the dynamic indirect tensile strength of coal: Laboratory experiments and numerical simulation,” International Journal of Coal Geology, vol. 132, pp. 81–93, 2014.

G. F. Zhao, Development of micro-macro continuum-discontinuum coupled numerical method [Ph.D. thesis], EPFL, Switzerland, 2010.

G.-F. Zhao, J. Fang, and J. Zhao, “A 3D distinct lattice spring model for elasticity and dynamic failure,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 35, no. 8, pp. 859–885, 2011.

Y. Zhao, G.-F. Zhao, and Y. Jiang, “Experimental and numerical modelling investigation on fracturing in coal under impact loads,” International Journal of Fracture, vol. 183, no. 1, pp. 63–80, 2013.