Research Article  Open Access
An Improved Statistical Damage Constitutive Model for Granite under Impact Loading
Abstract
To study the impact properties of granite, the parameters (including the stressstrain curve, elasticity modulus, peak strength, and peak strain) of the test pieces in each group were determined via standard splitHopkinson pressure bar tests. The results revealed that the prepeak stressstrain curves are approximately linear; the postpeak stressstrain curve declined sharply and exhibited the characteristics of brittle material failure after the stress exceeded the peak strength. In terms of the specimen form following failure, for increasing strain rate, the granite specimen became increasingly fragmented after failure. In addition, the singleparameter statistical damage constitutive model was improved, and a doubleparameter statistical damage constitutive model for describing the total stressstrain curve of granite under the action of impact loading was proposed. The parameters of the statistical damage model, m and a, were obtained via fitting. The results revealed that the parameter m decreases with increasing elasticity modulus, whereas the parameter a increases. Similarly, the peak strength and the peak strain increased (in general) with increasing strain rate.
1. Introduction
In many practical engineering problems, the applied load is dynamic. For example, rock blasting and excavation involve mechanical properties associated with the dynamic failure of rocks. The mechanical characteristics of impact have been extensively investigated. Vostretsov and Yakovitskaya [1] studied the mechanical properties of a preloaded rock specimen under the action of an impact load. Li et al. [2] proposed a novel technique for combined dynamic and static loading of a rock. The results revealed that, for this type of loading, the strength of the rock decreased significantly when the axial preload exceeds 70% of the static load strength of the rock. Wang et al. [3] studied the nonlinear damage of brittle rock under the action of an impact load. Furthermore, Millon et al. [4] investigated the mechanical characteristics of sedimentary rocks under the action of an impact load and the fragmentation of rock under different strain rates. Wu et al. [5] proposed a method for quantifying the fatigue damage of brittle rocks under an impact load. Peng et al. [6] assessed the failure process and the governing mechanism in sandstone under biaxial isostatic pressure and found that microcrack propagation is more likely under a combined load than under separate loads. Through numerical simulation, Bi et al. [7] evaluated the failure process of the rock material under the action of an impact load. Saksala [8] used a viscoplastic damage model to investigate the lowfrequency impact loading of rocks. Wang and Yang [9] studied the dynamic fracture toughness of coal rocks under the action of an impact load. Liu et al. [10] assessed the acoustic emission wavelet denoising threshold of rocks under a dynamic load. Scherbakov and Chmel [11] investigated the impact fracture characteristics of granite at different temperatures. In addition, Wang et al. [12–14] and Meng et al. [15] studied the mechanical behavior of clayey soil under different conditions and the weakening mechanisms of gypsum interlayers from Yunying salt, respectively. Tu et al. [16] proposed a new criterion for defining slope failure.
Previous studies have shown that the mechanical behavior of rocks under an impact load differs significantly from that under a static load. However, theoretical analyses and experimental studies focused on the dynamic failure of rock materials are rare, and thus, the dynamic characteristics of rocks under impact loading have received an increasing attention. Therefore, in the present work, granite specimens are subjected to a splitHopkinson pressure bar (SHPB) test. The model parameters of specimens in each group are obtained via fitting based on the improved statistical damage constitutive model. The determined impact properties of granite are important from both theoretical and practical points of view.
2. Tests
2.1. Testing Process
Cylindrical 63 mm × 31.5 mm black granite specimens with a diameter (D)toheight (H) ratio of D : H = 2 : 1 were tested. To guarantee the flatness of the upper and lower surfaces, these surfaces were symmetrically ground with a grinder. In the standard SHPB test, specimens with a diametertoheight ratio of 2 were investigated (Figure 1). The SHPB consisted of three parts (Figure 2 for the loading device): bullets, incident bar, and transmission bar. When a bullet strikes the incident bar, an incident wave will form, then propagate to the specimen along the bar, and interact with the specimen, resulting in a reflected wave and a transmission wave. The three types of wave can be measured by the strain gauge on the bar and acquired by an oscilloscope via an electrical bridge and ultrahigh dynamic strainometer.
2.2. Test Signal
Figure 3 presents the waveform signal of specimen #6 during testing. The black line and the red line denote the incident bar signal and the transmission bar signal, respectively. Incident waves, reflected waves, and transmission waves are typical triangular waves. The waveform suggests that the dynamic impact loading is a linear mode of loading and unloading that occurs immediately after the failure strength is reached, consistent with the dynamic failure characteristics of brittle materials.
2.3. Test Result
Figure 4 presents the dependence of the dynamic compression strength on the time. As shown in the figure, the dynamic failure process of the rock lasts for several dozens of microseconds. The minimum and maximum strengths are 189 MPa and 280 MPa, respectively. Most specimens have compression strengths of 260 MPa, except for specimens #2, #13, and #19, which have a slightly lower strength of 190 MPa–200 MPa.
Figure 5 presents the stressstrain curve of dynamic compression. The prepeak stressstrain curves are approximately linear. However, the postpeak stressstrain curve declines rapidly after the stress exceeds the peak strength. An average dynamic compression strength of 240 MPa (average error: 18.75%) is determined. Most other failures occur at strain ranging from 7% to 9%, except in the case of specimen #8 where failure occurs at a strain of >11%, consistent with the failure characteristics of brittle materials. The fitted stressstrain curve is described by a linear elastic relationship, and an average dynamic modulus E_{d} of 31.5 GPa is determined.
Figure 6 shows the dependence of the dynamic compression strength on the strain rate ranging from 81 to 210 s^{−1}. The strength increases significantly with increasing strain rate, indicating that the granite is hardened when withstanding the dynamic compression. Table 1 presents the test result of the granite under the action of the impact load.

2.4. Analysis of Failure Modes
The failure modes of some specimens are shown in Figure 7. Overall, the specimens are broken into pieces of different sizes under dynamic loading at different strain rates. The strain rates of specimens #2, #4, #6, #8, #13, and #19 are 81.90 s^{−1}, 81.26 s^{−1}, 107.02 s^{−1}, 185.98 s^{−1}, 97.34 s^{−1}, and 99.29 s^{−1}, respectively. With increasing strain rate, the fragment size and impact compressive strength of the granite specimens decrease and increase, respectively, after failure.
(a)
(b)
(c)
3. Improved Statistical Damage Constitutive Model of Granite under Impact Load
3.1. Improved Statistical Damage Constitutive Model
According to Krajcinovic and Silva [17], the Weibull distribution can be used to describe the internal damage of materials. Therefore, in this work, the microunit strength of granite is described by the Weibull distribution, with a probability density distribution function, which is given aswhere and are the probability density distribution function and strength distribution variable of the microunit, respectively, while and are the Weibull distribution parameters.
The damage degree can then be expressed as
The microunit strength based on the Drucker–Prager failure criteria iswhere is the internal frictional angle of the material, is the first stress invariant of the stress tensor, and is the second stress invariant of the stress tensor.
Li et al. [18] have reported that the microunit strength under a onedimensional stress state is given aswhere is the axial strain.
The improved timedependent damage models based on the Kelvin model includewhere , , and are the axial stress, elasticity modulus, and viscosity coefficient, respectively.
Substituting formula (2) into formula (5) gives
At the peak of the stressstrain curve , and . Taking the derivative of formula (6) yields
By substituting formula (7) into formula (6), we obtain
Substituting formulas (4) and (7) into formula (2) yields
The statistical damage computational formula, expressed as equation (9), consists of only one parameter, i.e., m. The prepeak damage and postpeak stressstrain curve slope are inadequately reflected by a singleparameter statistical damage constitutive model. In other words, the specimens with similar prepeak damage may have different postpeak stressstrain curve slopes, suggesting that the singleparameter statistical damage constitutive model has limitations. For this reason, the present work proposes an improved computational formula for statistical damage:where and are statistical damage parameters. Parameter a reflects the peak strength of the specimen; i.e., the peak strength of the specimen increases with decreasing . The parameter reflects the decline slope of the postpeak stressstrain curve; i.e., the decline slope of the curve increases with increasing .
By substituting formulas (8) and (10) into formula (5), we obtain the expression describing the improved statistical damage constitutive model:
This formula describes the improved statistical damage constitutive model based on the Weibull distribution.
3.2. Fitting Result of the Improved Statistical Damage Constitutive Model
The fitting result of the improved statistical damage constitutive model is shown in Figure 8 and Table 2. As the table shows, good fitting results are obtained. A correlation coefficient of R^{2} > 0.9390 is obtained for all specimens, except for specimen #13. The poor fitting result of this specimen results from the occurrence of significant plastic deformation and downward bending of the prepeak stressstrain curve, as shown in Figure 8(e). This trend differs significantly from the characteristics of other specimens, which are described by approximately linear prepeak curves. Based on Tables 1 and 2, parameter m decreases with increasing elasticity modulus, whereas parameter a increases. Similarly, the peak strength and the peak strain increase (in general) with increasing strain rate.
(a)
(b)
(c)
(d)
(e)
(f)

4. Conclusion
A statistical damage constitutive model is used to investigate the behavior of granite under impact loading. The major conclusions are summarized as follows:(1)Under the impact load, the average dynamic modulus and average peak strength are 31.5 GPa and 240 MPa, respectively; the average peak strain ranges from 7% to 9%, consistent with the failure characteristics of brittle materials.(2)The failure morphology of the granite samples indicates that, for increasing strain rate, the samples become increasingly fragmented after failure, and the impact compressive strength increases.(3)The present work proposes an improved statistical damage constitutive model for granite subjected to an impact load and reveals the physical significance of the model parameters, m and a. The improved statistical damage constitutive model is used to fit the test data. A good fitting result is obtained, and the association coefficients are largely larger than 0.9390. Investigation of the parameters obtained via fitting revealed that parameter m decreases with increasing elasticity modulus, whereas parameter a increases. Similarly, the peak strength and the peak strain increase (in general) with increasing strain rate.
The impact behavior of rocks (including granite) is a complex scientific problem. In this work, the statistical damage theory is used to study this behavior from a macroscopic point of view and, hence, the present study has certain limitations. Therefore, future work will also consider the behavior from a microscopic point of view.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the Fujian Provincial Natural Science Foundation Projects (grant no. 2018J01624), National Natural Science Foundation of China (grant nos. 51478118 and 51678164), Key Projects of Guangxi Natural Science Foundation (grant no. 2017JJD150035), Guangxi Science and Technology Base and Special Foundation For Talents (grant no. 2017AD23050), Open Fund of Key Laboratory of Underground Engineering (grant no. KFT18015), and Zhejiang Provincial Natural Science Foundation of China (grant no. LY16E040002). The authors gratefully acknowledge these supports.
References
 A. G. Vostretsov and G. E. Yakovitskaya, “Effect of external pulsed lowenergy impact on destruction of preloaded rock specimens,” Journal of Mining Science, vol. 50, no. 6, pp. 1033–1039, 2014. View at: Publisher Site  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
 Y. X. Wang, P. Cao, Y. H. Huang, R. Chen, and J. T. Li, “Nonlinear damage and failure behavior of brittle rock subjected to impact loading,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 13, no. 1, pp. 61–68, 2012. View at: Publisher Site  Google Scholar
 O. Millon, M. L. RuizRipoll, and T. Hoerth, “Analysis of the behavior of sedimentary rocks under impact loading,” Rock Mechanics and Rock Engineering, vol. 49, no. 11, pp. 4257–4272, 2016. View at: Publisher Site  Google Scholar
 B. Wu, P. Kanopoulos, X. Luo, and K. Xia, “An experimental method to quantify the impact fatigue behavior of rocks,” Measurement Science and Technology, vol. 25, no. 7, Article ID 075002, 2014. View at: Publisher Site  Google Scholar
 J. Y. Peng, Y. H. Li, F. P. Zhang, and Z. G. Qiu, “Failure process and mechanism of sandstone under combined equal biaxial static compression and impact loading,” Strain, vol. 54, no. 2, Article ID e12267, 2018. View at: Publisher Site  Google Scholar
 J. Bi, X. P. Zhou, and X. M. Xu, “Numerical simulation of failure process of rocklike materials subjected to impact loads,” International Journal of Geomechanics, vol. 17, no. 3, Article ID 04016073, 2017. View at: Publisher Site  Google Scholar
 T. Saksala, “Damageviscoplastic consistency model with a parabolic cap for rocks with brittle and ductile behavior under lowvelocity impact loading,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 34, no. 13, pp. 1362–1386, 2010. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 X.L. Liu, Z. Liu, X.B. Li, M. Rao, and L.J. Dong, “Wavelet threshold denoising of rock acoustic emission signals subjected to dynamic loads,” Journal of Geophysics and Engineering, vol. 15, no. 4, pp. 1160–1170, 2018. View at: Publisher Site  Google Scholar
 I. P. Scherbakov and A. E. Chmel, “Impact fracture of granite at temperatures from 20 to 500°C,” Russian Geology and Geophysics, vol. 55, no. 10, pp. 1223–1228, 2014. View at: Publisher Site  Google Scholar
 Y. X. Wang, S. B. Shan, C. Zhang, and P. P. Guo, “Seismic response of tunnel lining structure in a thick expansive soil stratum,” Tunnelling and Underground Space Technology, vol. 88, pp. 250–259, 2019. View at: Publisher Site  Google Scholar
 Y. Wang, P. Guo, X. Li, H. Lin, Y. Liu, and H. Yuan, “Behavior of fiberreinforced and limestabilized clayey soil in triaxial tests,” Applied Sciences, vol. 9, no. 5, p. 900, 2019. View at: Publisher Site  Google Scholar
 Y. X. Wang, P. P. Guo, H. Lin et al., “Numerical analysis of fiberreinforced soils based on the equivalent additional stress concept,” International Journal of Geomechanics, 2019. View at: Google Scholar
 T. Meng, Y. Hu, R. Fang, Q. Fu, and W. Yu, “Weakening mechanisms of gypsum interlayers from Yunying salt cavern subjected to a coupled thermohydrochemical environment,” Journal of Natural Gas Science and Engineering, vol. 30, pp. 77–89, 2016. View at: Publisher Site  Google Scholar
 Y. Tu, X. Liu, Z. Zhong, and Y. Li, “New criteria for defining slope failure using the strength reduction method,” Engineering Geology, vol. 212, pp. 63–71, 2016. View at: Publisher Site  Google Scholar
 D. Krajcinovic and M. A. G. Silva, “Statistical aspects of the continuous damage theory,” International Journal of Solids and Structures, vol. 18, no. 7, pp. 551–562, 1982. View at: Publisher Site  Google Scholar
 X. B. Li, S. M. Wang, F. Q. Gong, H. P. Ma, and F. P. Zhong, “Experimental study of damage properties of different ages concrete under multiple impact loads,” Chinese Journal of Rock Mechanics and Engineering, vol. 31, no. 12, pp. 2464–2472, 2012, in Chinese. View at: Google Scholar
Copyright
Copyright © 2019 Zhenwei Zhao 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.