Shock and Vibration
Volume 2015 (2015), Article ID 167580, 9 pages
http://dx.doi.org/10.1155/2015/167580

Research Article

Acoustic Emission Behavior of Rock-Like Material Containing Two Flaws in the Process of Deformation Failure

Quan-Sheng Liu,1,2 Jie Xu,¹ Bin Liu,² and Jing-Dong Jiang¹

¹School of Civil Engineering, Wuhan University, Wuhan, Hubei 430072, China
²Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China

Received 17 October 2014; Revised 14 February 2015; Accepted 17 February 2015

Academic Editor: Shimin Liu

Copyright © 2015 Quan-Sheng Liu 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.

Abstract.

Many sudden disasters (such as rock burst) by mining extraction originate in crack initiation and propagation. Meanwhile a large number of shock waves are produced by rock deformation and failure. With the purpose of investigating crack coalescence and failure mechanism in rock, experimental research of rock-like materials with two preexisting flaws was performed. Moreover, the AE technique and photographic monitoring were adopted to clarify further the procedure of the crack coalescence and failure. It reveals that AE location technique can record the moments of crack occurrences and follow the crack growth until final failure. Finally, the influence of different flaw geometries on crack initiation strength is analyzed in detail. This research provides increased understanding of the fracture mechanism of mining-induced disasters.

1. Introduction

With the development of underground engineering excavating technology, more and more new equipment and new strategies are applied into mining industry. However, the capacity of mining disaster prevention and mitigation cannot meet the productivity demands. Many sudden disasters (such as rock burst) by mining extraction originate in the linkage of preexisting flaws due to the initiation, propagation, and interaction of new cracks and existing fractures. Meanwhile a large number of shock waves are produced by rock deformation and failure. Therefore, it is very urgent and important to understand the mining-induced shock waves during crack initiation, coalescence, and failure of rock material.

Many experimental investigations [1–4] have been carried out on enhancing the understanding of crack coalescence. The experimental results indicated that the mechanical behaviors of rock failure were closely bound up with the geometries of flaws. Meanwhile, the coalescence behavior of cracks is of great importance and difficulty. Following the development of test research, rock-like materials gain more attention and have been widely studied.

Vásárhelyi and Bobet [5], Sagong and Bobet [6], Zhou et al. [7], and Zhu et al. [8] have investigated the crack coalescence behavior by loading rock-like materials with different flaw geometries. Wong and Einstein [9] carried out an experimental study on cracking and coalescence behavior with molded gypsum and marble samples, which showed that the cracking processes and macroscopic deformation were bound up with the underlying microscopic change. Based on the combination of numerical simulations and experiments, Haeri et al. [10] have investigated the crack coalescence mechanism and propagation paths according to Mode I and Mode II stress intensity factors. In order to research the effect of flaw angle and distribution density on the failure characteristics, Pu and Cao [11] have performed experimental study on rock-like samples with prefabricated closed multiflaws. In accordance with cracking behaviors in modeled gypsum and Carrara marble samples, Wong and Einstein [12] have obtained seven different crack types. The reason for choosing rock-like materials is threefold: first, a large number of samples can be easily prepared; second, a wide range of brittle rocks can be represented; third, many of previous experimental results can be compared with the current investigation.

Acoustic emission technique (AE) has the capability to efficiently reflect the initiation and coalescence of cracks by tracing the mining-induced shock waves, which has been widely applied to many studies [13, 14] of crack behavior in the rock material. With the purpose of better understanding frost weathering mechanisms of rocks, Duca et al. [15] conducted an experimental study on the stored elastic energy by AE technique, which was released accompanied by the growth of microcracks during freezing tests. A series of AE location tests were carried out by Ting et al. [16] under unloading conditions, which indicated that energy increasing rate and AE space-time evolvement rules during strength failure process were dependent on the unloading rates. Jouniaux et al. [17] studied the fracturing and slip behaviors of the samples containing a preexisting macroscopic healed joint, which showed that precursory localization of microfractures in the final rupture plane had been observed in the early stage of deformation by AE technique. Based on the analysis of AE parameters in microcracking activity, Vidya Sagar et al. [18] found the AE behavior characteristics of each stage were linked closely with the number and distribution of microcracks. In order to identify the characteristics of crack coalescence behaviors in sandstone samples, Yang et al. [19–21] adopted the AE technique to record the real-time crack coalescence process of the samples in the whole loading course. All studies show that AE has the capacity to efficiently highlight the moments of crack occurrences and follow the crack growth until final failure.

With the purpose of figuring out the strength failure mechanism and crack coalescence behavior in rock materials, AE technique and photographic monitoring were used in the uniaxial compression tests for samples containing two parallel flaws. Moreover, the whole process of the real-time crack coalescence was also recorded by tracing mining-induced shock waves in the whole loading course, which was not carried out in rock-like materials previously. In this paper, the researching emphasis centralizes in investigating the effect of the geometry (flaw angle and ligament angle) on the strength failure mechanism of rock-like materials and in studying real-time crack coalescence process under uniaxial compression on the foundation of photographic monitoring and AE technique combination.

2. Sample Materials and Experimental Techniques

2.1. Sample Preparation

The samples were made of high-strength gypsum and water-cement ratio is 1. The size and flaw geometry of the samples are shown in Figure 1. Each sample was cast in a mold with internal dimensions of 60 mm × 40 mm × 120 mm. Meanwhile, the open flaws were produced with different-sized metallic shims, which are 0.2 mm in thickness, 60 mm in length, and 10–30 mm in width. In order to further investigate the influence of flaw geometry on the crack coalescence behavior, different geometries of two parallel flaws were chosen, which were set by varying inclination angle ( and ) while keeping two constants ( mm and  mm). Two flaw inclination angles were used: 45° and 75°; six ligament angles were used: 45°, 60°, 75°, 90°, 105°, and 120°.

(a)

(b)

Figure 1: Sample with two flaws. (a) Overall view and (b) geometry of the flaws.

To obtain the physicomechanical parameters of the rock-like materials, uniaxial and biaxial compression tests were done with the standard cylindrical samples of 50 mm in diameter and 100 mm in length and Brazilian splitting tests were done with the samples of 50 mm in diameter and 50 mm in length. The physicomechanical parameters of intact samples are shown in Table 1.

Table 1: The physicomechanical parameters of intact samples.

(MPa)
(MPa)
(GPa)

(°)
(g/cm³)
P-wave velocity (m/s)
S-wave velocity (m/s)

5.45
0.90
1.30
2.44
9.03
0.84
3123
2657

Note: = uniaxial compressive strength, = tensile strength, = internal friction angle, = cohesive force, and = density.

2.2. Test Equipment

A servo-controlled hydraulic testing machine, namely, RMT-150C (Figure 2(a)), was used in the experiment, which can perform well in load or displacement control. The maximum capacity of RMT-150C was 1 MN and all the real-time data were recorded and analyzed in time. The samples were loaded in the displacement-controlled condition and the stain rate was set as 1 × 10⁻³ mm/s. In addition, two rigid steel blocks were placed between the testing samples and loading frame. With the purpose of decreasing the influence of the friction effects and obtaining more accurate testing results, plastic cushions were placed between the end face of samples and rigid steel blocks, as shown in Figure 2(b).

(a) RMT-150C

(b) Experimental setup for the sample

Figure 2: RMT-150C testing system and experimental setup of samples.

A multichannel, high-speed AE testing system, namely, PCI-2, was used to acquire the real-time AE information, which was manufactured by American Physical Acoustic Corporation (Figure 3(a)). The AE system has become mature as a technique for monitoring the crack growth in rock materials, which is composed of AE transducers, preamplifier, signal acquisition, and so on. In this study, AE signals were amplified by 40 dB, and the amplitude threshold, the sampling frequency, and the sampling length were set at 18 mV, 2.5 MHz, and 8192, respectively. Eight Nano-30 type sensors are used to acquire AE signals, of which the operating frequency is 100–400 kHz. Each sensor was equipped with a 1220A-AST type preamplifier. The sensors were fixed on rock faces by plastic tapes and Vaseline was used for coupling, as shown in Figure 3(b). Before the test, a pencil lead break (PLB) was conducted to check out coupling quality of the acoustic emission sensor with the specimen and correct the AE location system, shown in Figure 3(b). If the amplitude value tested by one AE sensor was above 95 dB, the sensor could be considered to have a good coupling quality with the specimen. If the location of the AE event tested by AE setup was close to the actual location, this system could achieve the demand of position control. Otherwise the PLB tests should be repeated and setup parameters should be modified until the location has met the requirements. Besides, the velocity of longitudinal wave and transverse wave is 3123 m/s and 2657 m/s, respectively.

(a) AE PCI-2

(b) Setup for sensors

(c) PLB

Figure 3: AE PCI-2 acoustic emission system.

3. Strength and Deformation Behavior

Take flaw angles = 45° and = 75°, for example; the uniaxial compression curves for the samples containing two flaws are shown in Figure 4. It can be found that flaw angle and ligament angle make a great difference in the strength and deformation properties of rock-like materials, and the detailed influence will be analyzed in the following.

(a) = 45°

(b) = 75°

Figure 4: Axial stress-strain curves of the samples, in which (a) and (b), respectively, show the flaw angles = 45° and = 75° ( = axial stress and = axial strain) (black, red, blue, dark cyan, magenta, and dark yellow correspond to 45°, 60°, 75°, 90°, 105°, and 120° of ligament angle , resp.).

3.1. Relationship between Stress and Strain for the Samples with Different Flaw Geometries

The relationship between stress and strain for the samples containing two parallel flaws can go through approximately four classic stages, that is, flaw closure, elastic deformation, crack initiation and propagation, and strain-softening.

At the initial stage of the stress-strain curves, the downward concave and nonlinear characteristics are clearly observed, which result from the closure of some primary pores and voids at low stress levels and are defined as flaw closure. Moreover the better consistency of the rigidity for all samples at this stage lies in the better homogeneity of tested materials. With the increase of stress, the relationship of stress and strain is at the stage where the linear stress-strain behavior of flawed samples keeps unchanged although there are some irrecoverable processes, which is defined as elastic deformation. Moreover, the slope of the line at elastic deformation stage is dependent distinctly on ligament angle . In accordance with Figure 4, it is clear that the slope of the line increases with the increase of ligament angle . Nevertheless, at the stage of crack initiation and propagation, new cracks will initiate and propagate from the tips of the preexisting flaws, which result from stress concentration. Therefore the stress-strain curve shows distinctly nonlinear deformation before peak strength point. From Figure 4, it is clear that the peak strength of samples with ligament angle = 105° or 120° is obviously greater than that of others, so one can conclude that it is more difficult to fail for overlapping flaw geometries () than nonoverlapping ones (). When the relationship of stress and strain develops into the strain-softening stage, the macroscopic crack occurs and then the postpeak behavior shows a slower fall.

3.2. Effect of Different Flaw Geometries on Mechanical Parameters of Samples

In order to investigate the influence of different flaw geometries on the strength and deformation behaviors, the mechanical parameters of the samples under uniaxial compression are discussed in detail, which are shown in Figure 5.

(a) Uniaxial compressive strength

(b) Elastic modulus

Figure 5: Effect of flaw angle and ligament angle on mechanical parameters of samples.

The uniaxial compressive strength of intact sample is 5.450 MPa. Corresponding to the uniaxial strength of the integrity sample, for the samples with = 45° or = 75°, the average reduction extent of uniaxial compressive strength is 32.02% or 12.76%, respectively. According to Figure 5(a), the uniaxial compressive strength of the samples containing flaw angle = 75° is greater than that of samples containing flaw angle = 45°; therefore, one can conclude that the samples with flaw angle = 75° have the ability to bear greater pressure. Moreover, when the angle is increased from 45 to 75 degrees and from 45 to 90 degrees for angles of 75 and 45 degrees, respectively, the uniaxial compressive strength keeps almost unchanged. Then it begins to increase with increasing of ligament angle and the demarcation point is = 90° for = 45° and = 75° for = 75°.

In accordance with Figure 5(b), the change rule of the elastic modulus is different from that of uniaxial compressive strength shown in Figure 5(a). For the samples with = 45°, there is a distinct valley when = 75°; then the elastic modulus increases with the increasing ligament angle. However, the elastic modulus of the samples with = 75° keeps an increasing trend from beginning to end, which has a demarcation point at = 60° because of the different slope of increase.

From the above analysis, one can conclude that the mechanical parameters of the samples with = 75° are greater than those with = 45°, which shows that the flaw angle makes a significant difference on the strength and deformation behaviors. In addition, when , all mechanical parameters increase with increasing ligament angle, which shows that the overlapping flaw geometries have contributed more to the stability of the samples.

4. Real-Time Recording and Analysis of Crack Coalescence Process for Flawed Samples

Besides AE energy technique, the AE location-based technology and photographic monitoring were also adopted to illuminate further the characteristics of crack coalescence behaviors in flawed samples under uniaxial compression. It is generally known that the located AE events directly reflect the time-space distribution of the catastrophic fracture within the samples. Based on the combined results, the real-time recording and analysis of crack coalescence process for flawed samples are shown in detail as in Figure 6.

Figure 6: Real-time recording of crack coalescence process for flawed sample ( = 45° and = 75°).

In accordance with Figures 6, 7, and 8, the crack coalescence process of the sample containing two preexisting flaws ( = 45° and = 75°) in real time was analyzed in detail. Before point A () as shown in Figures 7 and 8, the stress concentration is too small to reach the material strength, so that there is no crack initiation at the tips of flaws. For the moment, the average modulus is 1.09 GPa. But a wing crack begins to initiate from the outer tips of the preexisting flaws as soon as the sample is loaded to point A, which is identified by AE energy shown in Figure 7. With the axial stress increasing, the wing crack further lengthens and continues to expand along the direction of loading. Due to some minor damage of structural supports inside the material, the average modulus (1.0 GPa) becomes a little lower than the elastic modulus 1.09 GPa. When the sample is loaded to the first yielding point B (), a vertical tensile crack occurs rapidly in the bridge in one second, which can be observed from the AE energy shown in Figure 7. At this time, the axial stress undergoes a minor drop to 3.074 MPa from 3.152 MPa. However it is more and more difficult for the wing crack to continue propagating because of the limitation of sample boundary. As time grows, the axial supporting capacity of the sample becomes increasingly great. At the moment of axial stress loaded to point C (), the second wing crack initiates from the outer tip of the upper flaw. What is more, the average modulus (~0.262 GPa) becomes significantly lower than previous elastic modulus, which results from the serious damage of structural supports inside the materials. When the sample was loaded to the second yielding point D (), two shear cracks emanate concurrently from the internal tips of two flaws, which indicates that the coalescence occurs. This moment is captured by the AE technology, as shown in Figure 7. During this loading process, it is much in evidence that the axial supporting capacity turns into growing slowly compared to previous period. In addition, the critical decline of average modulus (~0.096 GPa) is largely due to serious damage of structural supports inside the materials. Afterwards, with the sample being loaded to 3.465 MPa (98.5% ), surface spalling starts to appear on the surface of the sample, which can be obviously observed from Figure 6. When the sample is loaded to 3.518 MPa (100% ), at this moment, the first wing crack develops rapidly to the edge of the sample in one second, which leads to the strength failure. At the same time, the axial supporting capacity reduces to 2.035 MPa.

Figure 7: Relation between AE energy rate and time of flawed sample ( = 45° and = 75°).

Figure 8: Several important points in stress-strain curve for flawed sample ( = 45° and = 75°).

The crack coalescence process of flawed sample in real time can be confirmed and characterized by virtue of AE location technique as well, as shown in Figure 9. At the beginning of load, the AE events are so little that no rules can be obtained. When the sample is loaded to 2.25 MPa (point A), more AE events are observed at the outer tips of preexisting flaws, indicating that wing crack occurs. At the moment, the number of AE events only accounts for 6.9% of the total. Afterwards, with the increasing axial stress, a red zone can be distinctly observed in the bridge, which indicates the vertical tensile crack occurs. However the continuous increasing stress leads the AE events to increase rapidly. At this moment of about 406 s (corresponding axial stress is 3.357 MPa), the coalescence occurs, which is shown in Figure 9(e). Meanwhile, the proportion of AE events has risen to 79.3%. When the sample is loaded to 3.518 MPa (corresponding time is 425 s), a vertical red zone can be observed through the whole sample, which results from the failure of the sample (Figure 9(g)).

(a) = 1.69 MPa ( = 207 s)

(b) = 2.25 MPa ( = 274 s)

(c) = 3.15 MPa ( = 381 s)

(d) = 3.178 MPa ( = 385 s)

(e) = 3.357 MPa ( = 406 s)

(f) = 3.465 MPa ( = 419 s)

(g) = 3.518 MPa ( = 425 s)

Figure 9: The AE location for cracks initiating and propagating in flawed sample ( = 45° and = 75°).

By adopting the combination of photographic monitoring and AE technique, the crack coalescence behavior and failure process of all other samples with different flaw geometries were also recorded in real time. One can conclude that each bigger elastic modulus decline in all the compression curves signifies a bigger crack initiation, which is also accompanied by greater amplifying AE energy and more AE location events. Figure 10 summarizes the stress of first crack initiation, which also is called crack initiation strength.

Figure 10: Crack initiation strength of the samples with different geometries.

According to Figure 10, we can draw the conclusion that the crack initiation strength is linked closely with flaw angle and ligament angle. To explain in concrete terms, the crack initiation strength in the sample with = 75° is higher than that with = 45° for each ligament angle except = 45°; that is, the initiation of first crack in the sample with = 45° is easier than that with = 75°. In addition, for the sample with = 45°, the crack initiation strength showed a trend of fluctuations when the ligament angle was lower than 90°. When the ligament angle increased from 45° to 60°, the crack initiation strength increased from 3.1 MPa to 3.25 MPa. Then it decreased to 3.17 MPa at = 75° and 3.12 MPa at = 90°. After that the crack initiation strength increased with the increase of ligament angle. Similarly, for the sample with = 75°, the crack initiation strength was 4.52 MPa at = 60°. Then it decreased to 3.8 MPa at = 75°. After that it increased to the maximum value of 5.15 MPa at = 120°. What is more, we can come to a conclusion that the first crack in the sample with smaller ligament angle is easier to initiate than that with larger ligament angle on the whole. Particularly for the sample with overlapping flaw geometries (), the crack initiation strength is distinctly higher than that of nonoverlapping flaw geometries ().

5. Conclusions

This study focused on the influence of the geometry (flaw angle and ligament angle) on the strength failure mechanism of rock-like materials and analyzed real-time crack coalescence process by tracing mining-induced shock waves on basis of AE technique. It provides better understanding of the fracture mechanism of mining-induced disasters. The main conclusions of this study are the following.

(1)Compared with intact sample, the flawed samples fail with lower strengths, which is closely related to the flaw angle and ligament . The mechanical parameters of the samples with = 75° are greater than those with = 45°. When , all mechanical parameters increase with increasing ligament angle, which shows that the overlapping flaw geometries have contributed more to the stability of the samples.(2)The crack coalescence behavior and strength failure process of the samples with different flaw geometries are recorded in real time by tracing induced shock waves. Each bigger elastic modulus decline in all the compression curves signifies a bigger crack initiation, which is also accompanied by greater amplifying AE energy and more AE location events.(3)The crack initiation strength is linked closely with flaw angle and ligament angle. The crack initiation strength in the sample with = 75° is higher than that with = 45° for each ligament angle except = 45°; that is, the initiation of first crack in the sample with = 45° is easier than that with = 75°. Moreover, the initiation of first crack in the sample with smaller ligament angle is easier than that with larger ligament angle on the whole. Particularly for the sample with overlapping flaw geometries (), the crack initiation strength is distinctly higher than that of nonoverlapping flaw geometries ().

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Financial supports from the National Basic Research Program of China (973 Program) Grant no. 2014CB046904 and the National Natural Science Foundation of China Grant nos. 41130742 and 51474205 are gratefully acknowledged.

References

G. Chen, J. Kemeny, and S. S. Harpalani, “Fracture propagation and coalescence in marble plates with pre-cut notches under compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 30, p. 279, 1993.

H. Jiefan, C. Ganglin, Z. Yonghong, and W. Ren, “An experimental study of the strain field development prior to failure of a marble plate under compression,” Tectonophysics, vol. 175, no. 1–3, pp. 269–284, 1990.

L. N. Y. Wong and H. Einstein, “Fracturing behavior of prismatic specimens containing single flaws,” in Proceedings of the 41st US Rock Mechanics Symposium (Golden Rocks '06), Omni Press, Golden, Colo, USA, June 2006.

L. G. Zhao, J. Tong, and M. C. Hardy, “Prediction of crack growth in a nickel-based superalloy under fatigue-oxidation conditions,” Engineering Fracture Mechanics, vol. 77, no. 6, pp. 925–938, 2010.

B. Vásárhelyi and A. Bobet, “Modeling of crack initiation, propagation and coalescence in uniaxial compression,” Rock Mechanics and Rock Engineering, vol. 33, no. 2, pp. 119–139, 2000.

M. Sagong and A. Bobet, “Coalescence of multiple flaws in a rock-model material in uniaxial compression,” International Journal of Rock Mechanics and Mining Sciences, vol. 39, no. 2, pp. 229–241, 2002.

X. P. Zhou, H. Cheng, and Y. F. Feng, “An experimental study of crack coalescence behaviour in rock-like materials containing multiple flaws under uniaxial compression,” Rock Mechanics and Rock Engineering, vol. 47, no. 6, pp. 1961–1986, 2013.

W. S. Zhu, W. Z. Chen, and J. Shen, “Simulation experiment and fracture mechanism study on propagation of Echelon pattern cracks,” Acta Mechanica Solida Sinica, vol. 19, pp. 355–360, 1998.

L. N. Y. Wong and H. H. Einstein, “Crack coalescence in molded gypsum and Carrara marble. Part 1. Macroscopic observations and interpretation,” Rock Mechanics and Rock Engineering, vol. 42, no. 3, pp. 475–511, 2009.

H. Haeri, K. Shahriar, M. F. Marji, and P. Moarefvand, “Cracks coalescence mechanism and cracks propagation paths in rock-like specimens containing pre-existing random cracks under compression,” Journal of Central South University, vol. 21, no. 6, pp. 2404–2414, 2014.

C.-Z. Pu and P. Cao, “Failure characteristics and its influencing factors of rock-like material with multi-fissures under uniaxial compression,” Transactions of Nonferrous Metals Society of China, vol. 22, no. 1, pp. 185–191, 2012.

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.

A. Zang, F. C. Wagner, S. Stanchits, G. Dresen, R. Andresen, and M. A. Haidekker, “Source analysis of acoustic emissions in Aue granite cores under symmetric and asymmetric compressive loads,” Geophysical Journal International, vol. 135, no. 3, pp. 1113–1130, 1998.

X. Lei, K. Kusunose, M. V. M. S. Rao, O. Nishizawa, and T. Satoh, “Quasi-static fault growth and cracking in homogeneous brittle rock under triaxial compression using acoustic emission monitoring,” Journal of Geophysical Research B: Solid Earth, vol. 105, no. 3, pp. 6127–6139, 2000.

S. Duca, C. Occhiena, and L. Sambuelli, “A procedure to determine the optimal sensor positions for locating AE sources in rock samples,” Rock Mechanics and Rock Engineering, vol. 48, no. 2, pp. 481–493, 2014.

A. Ting, Z. Ru, L. Jianfeng, and R. Li, “Space-time evolution rules of acoustic emission location of unloaded coal sample at different loading rates,” International Journal of Mining Science and Technology, vol. 22, no. 6, pp. 847–854, 2012.

L. Jouniaux, K. Masuda, X. Lei et al., “Comparison of the microfracture localization in granite between fracturation and slip of a preexisting macroscopic healed joint by acoustic emission measurements,” Journal of Geophysical Research B: Solid Earth, vol. 106, no. 5, pp. 8687–8698, 2001.

R. Vidya Sagar, R. V. Prasad, B. K. Raghu Prasad, and M. V. M. S. Rao, “Microcracking and fracture process in cement mortar and concrete: a comparative study using acoustic emission technique,” Experimental Mechanics, vol. 53, no. 7, pp. 1161–1175, 2013.

S. Q. Yang and H. W. Jing, “Strength failure and crack coalescence behavior of brittle sandstone samples containing a single fissure under uniaxial compression,” International Journal of Fracture, vol. 168, no. 2, pp. 227–250, 2011.

S. Q. Yang, Y. H. Dai, L. J. Han, and Z. Q. Jin, “Experimental study on mechanical behavior of brittle marble samples containing different flaws under uniaxial compression,” Engineering Fracture Mechanics, vol. 76, no. 12, pp. 1833–1845, 2009.

S.-Q. Yang, “Crack coalescence behavior of brittle sandstone samples containing two coplanar fissures in the process of deformation failure,” Engineering Fracture Mechanics, vol. 78, no. 17, pp. 3059–3081, 2011.