Related Rule Study of Subcritical Crack Growth and Threshold Values in Transversely Isotropic Slates

Ma, Gang; Li, Jiangteng; Wang, Huiwen

doi:https://doi.org/10.1155/2020/8843381

Geofluids

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Thermal-Hydraulic-Mechanical (THM) Coupling Behaviour of Fractured Rock Masses

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 8843381 | https://doi.org/10.1155/2020/8843381

Related Rule Study of Subcritical Crack Growth and Threshold Values in Transversely Isotropic Slates

Gang Ma,¹Jiangteng Li,¹and Huiwen Wang²

Academic Editor: Yanlin Zhao

Received09 Jun 2020

Revised31 Jul 2020

Accepted05 Oct 2020

Published30 Oct 2020

Abstract

Elastic parameters and the subcritical crack growth of different bedding angle slate specimens were studied using uniaxial compression testing and the double torsion constant displacement load relaxation method using SANS and MTS Insight machines. To study the relations of the mode-I stress intensity factor versus the subcritical crack growth velocity , the fracture toughness , the stagnation speed, and the threshold values, the double torsion constant displacement load relaxation method was carried out. The related rules between the bedding angles (β) and the uniaxial compressive strength, fracture toughness, and threshold values were investigated. Experimental results show that the uniaxial compression, the fracture toughness, and the threshold value curves move to the bottom then increase with the increase of the β angle. In addition, its fracture toughness is minimal when the β angle of the slates is 45°, and crack initiation and crack propagation are generated under load, which can lead to the failure of the slate. - relations of transversely isotropic slates measured by this method are in accordance with linear rules, which is in good agreement with the Charles theory. The range of for these different bedding angle slates is from 0.511 to 0.789. The test results would provide the basis for studying seepage and time-dependency of rock engineering stability.

1. Introduction

Rocks are mineral aggregates after a long geological process; rocks of layered structure can be simplified into a transversely isotropic body, whose physical and mechanical properties in parallel to the bedding plane are different from those in perpendicular to the bedding plane. The stability analysis of the layered rock is complicated by the isotropic characteristics of the bedding surface and the coupling effect of the seepage field and stress field [1]. Thus, to reveal the transversely isotropic permeability mechanism, it is necessary to investigate the transversely isotropic fracture characteristics.

At present, the researches on anisotropic rock properties are mainly concentrated on the following areas: (1) the measurement method study of elastic deformation parameters [2], (2) the anisotropy strength and yield criterion study [3, 4], (3) anisotropic rock constitutive model study [5, 6], and (4) anisotropic rock deformation characteristics and mechanical property study [7, 8]. Transverse isotropy is a special case of anisotropy; these studies have a great significance on helping us have a better understanding of the permeability and mechanical properties of transversely isotropic rocks. The instability instances in situ and rock laboratory examinations show that the rock failure is due to crack growth and propagation in rock mass [9], and the crack is the main seepage channel of the rock mass, which controls the seepage characteristics of the rock mass. According to the fracture mechanics and stress corrosion theory [10, 11], when a crack tip stress intensity factor exceeds its fracture toughness , the crack will quickly expand; when the stress intensity factor is less than its fracture toughness, the crack growth will process at the speed of stress corrosion subcritical crack growth. The study shows that the crack in certain circumstances propagates until the material fractures with a collective effect of tensile stress and corrosion medium, which is called stress corrosion [12]. Usually, the mode-I stress intensity factor at the crack tip is regarded as a parameter to control the crack growth. The stress corrosion lower limit is ; when the stress intensity factor is less than , the crack will not extension; when the stress intensity factor is between and , the crack grows at a certain velocity; the reasons can be listed as follows: (1) the stresses at the crack tip increases in the process of forming a microcrack and interspace and (2) the chemical bond fracture stresses decrease by the effects of the environmental medium, so the crack grows at a stable and static state in the process, which is called subcritical crack growth. The crack growth in the rock mass happens after the subcritical cracks propagate to a certain extent. Therefore, studies of the subcritical crack growth laws and fracture properties are important to rock engineering long-term stability. The elastic parameters and the principle of subcritical crack growth can be obtained in the paper. The related laws between the bedding angles (β) and compressive strength, fracture toughness, and threshold values had been analyzed. The relations of the mode-I stress intensity factor versus the subcritical crack growth velocity and the fracture toughness will be able to predict when the crack will be extended to the critical crack length under certain stress conditions, the result of which can cause structural damage, which gives the permeability and mechanics study a new approach.

2. Transverse Isotropy of Different Bedding Angle Slates

2.1. Constitutive Relations of Transversely Isotropic Elastic Body

Transversely isotropic elastic body is made of an elastic symmetry plane, whose physical and mechanical properties in parallel to the elastic symmetry plane are different from those in perpendicular to the elastic symmetry plane. For the layered rock, engineering generally regards the bedding plane as an elastic symmetry plane; it is believed that it has the same physical and mechanical properties in parallel to the bedding plane which are different from those that are vertical to the elastic symmetry plane. As shown in Figure 1, the XOY plane is assumed as the elastic symmetry plane; the transversely isotropic elastic parameters should satisfy the following conditions [13]:

The transversely isotropic elastic body constitutive relation can be simplified as follows [14]: where is the compliance matrix, which can be can be calculated by where and are Young’s modulus of rock and Poisson’s ratio in the direction parallel to the transversely isotropic plane, respectively; and is Young’s modulus of rock and Poisson’s ratio in the direction perpendicular to the transversely isotropic plane, respectively; and is the shear modulus for the plane normal to the transversely isotropic plane.

Because the stiffness matrix and the compliance matrix are inverse matrix, so the stiffness matrix can be calculated as follows:

where .

When an orthotropic body degenerates into a transversely isotropic elastic body, the nine independent elastic parameters are reduced to the five independent elastic parameters, respectively; , , , , , and are not independent elastic parameters, because of .

2.2. Coordinate Conversion of the Anisotropic Elastic Constants

Due to the different mechanical properties of the rock in all directions in anisotropic media, the coordinate system often needs to be converted to solve the elastic parameters, and coordinate transformation of anisotropic elastic parameters generally follows the tensor rule.

The global coordinate system () and the local coordinate system () of a transversely isotropic body are shown in Figure 2; the local coordinate system around the -axis rotates the β angle to get the global coordinate system. Stiffness matrix in the local coordinate system can be converted to the stiffness matrix in the global coordinate system. It can be expressed as [13] where is the transformation matrix, which can be expressed as follows:

Transversely isotropic elastic parameters in the local coordinate system can be transferred into the elastic parameters in the global coordinate system, which are shown as follows: where can be calculated by the following equations:

2.3. Thermodynamic Constraint of Transversely Isotropic Constant

To make the rock strain energy positive in this experiment, the 5 elastic parameters of the transversely isotropic rock should meet the related inequality of hot mechanical constraints, which are expressed as follows [15]:

3. Elastic Parameters of Transversely Isotropic Slate Specimen

Elastic modulus, Poisson’s ratio, and uniaxial compressive strength values of the cylinder slates whose diameter is 50 mm and height is 100 mm with different bedding angles (0°, 30°, 45°, 70°, 80°, and 90°) can be obtained by uniaxial compression test using the SANS machine, which are listed in Table 1.

According to the translation of rock mechanics by Milou and the relation between the constant in compliance matrix and elastic constants of rock, the following relationship can be obtained [16]: where and are the elastic modulus and Poisson’s ratio of the slate with a bedding angle of 90°, respectively; and are the elastic modulus and Poisson’s ratio of the slate with a bedding angle of 0°, respectively; and is the elastic modulus of the slate with β bedding angle.

For the slates with a certain angle between the transversely isotropic plane and horizontal plane, , , and of the slate specimen in three orthogonal directions can be obtained by Equation (10), which are listed in Table 2. Based on the transversely isotropic and loading direction of the slate, we select as the shear modulus of the slate specimen in the calculation of its fracture toughness.

4. The Relations between Bedding Angles of Slate Specimen and Uniaxial Compressive Strength and Elastic Modulus

For the layered rock with a single layer, its bedding angles have a great influence on its destructive strength. For the slates in this experiment, it can be considered as the transversely isotropic constitutive model which has only one set of parallel beddings with different angles loaded. We can obtain the related law between the bedding angles of the slate specimen and the uniaxial compressive strength and elastic modulus from Table 1, which can be expressed as Figure 3, from which we can know that the bedding angles have a great effect on the uniaxial compressive strength and elastic modulus of the slates. What is more, the uniaxial compression and elastic modulus decreases then increases with the increase of β angle and the relation curves show a “U” type. The results show that when the β angle of the slate is 45°, its uniaxial compression and elastic modulus are minimal. Therefore, it is most likely to be damaged when the bedding angle of slates is 45°.

5. Double Torsion Method

5.1. Principle of Double Torsion Testing

The double torsion specimen was first adopted by Williams and Evans [17] to decide the fracture properties of glass and ceramics, which is shown in Figure 4. Since 1977, Ciccotti et al. [18], Saadaoui et al. [19], and Nara and Kaneko [20] applied it to rock research. The double torsion specimen can be assumed as two elastic torsion bars, of which each cross-section is a rectangle.

(a)

(b)

(c)

For the small deformation, when the width of the torsion bar is much larger than the thickness of the torsion bar, the torsion strain can be calculated by the following equation [21]: where is the torque , is the applied load on the torsion bar, is the shear modulus of the rock, is the crack length, is the specimen thickness, of the torsion bar’s width, and is the length of twisted arm, and Equation (11) can be rewritten as where is the elastic compliance of the double torsion specimen. The relation among crack strain energy release rate and the flexibility of the specimens satisfies the following relationship: where is the area of the crack and is the thickness of plane where crack exists. For the double torsion specimen, strain energy release rate can be defined as

Based on the relations between stress intensity factor and the strain energy release rate and the relationship between elastic modulus and shear modulus , a lot of research has been done, such as by Cao et al. [22], Andreas et al. [23], Ciccotti et al. [18, 24], and Zhao et al. [25]. The stress intensity factor is introduced in this study.

When the load reaches a critical value of the on the specimen, the crack quickly grows, So reaches a critical value , which is shown as follows:

Tests showed that the compliance change rate can be obtained by the displacement change rate with the given loads or by the load change rate with the given displacement. Under the constant displacement conditions, through the derivative of Equation (12), it can be obtained that

Equation (17) shows that if the specimen’s size and displacement are known, the subcritical crack growth velocity is related to the load relaxation rate with the given displacement state.

5.2. Stress Intensity Factor and Crack Velocity of the Double Torsion Method

The stress intensity factor was independent of the crack length in this study, and the double torsion test is convenient when it uses opaque materials, such as rocks for which crack length measurements are inconvenient. There are three methods to obtain the crack velocity in the double torsion tests, each of which has a different loading process. The constant load method was used by Li et al. for the original double torsion test [26]. Two other methods are the constant displacement rate and the constant displacement load relaxation method which were introduced by Evans [27]. In the constant displacement load relaxation method, the displacement of the loading point is kept constant during the experiment and the load relaxation is measured with the crack growth. Because the stress intensity factor is the function of the loading and the crack growth rate is the function of the temporal load and decreasing rate of the loading, the relation between the stress intensity factor and the crack growth rate can be obtained over a wide range of loads using a single experimental run. For this reason, the constant displacement load relaxation method was adopted in this paper.

5.3. Experimental

5.3.1. Specimen Preparation

Slates of different bedding angles (0°, 30°, 45°, 70°, 80°, and 90°) were machined as a rectangle plate double torsion specimen whose standard size was and width-to-thickness ratio was 15 : 1 ~ 12 : 1. There was a longitudinal groove on the bottom surface, which was 1 mm in width and in depth, and a notch was cut at the end. The crack front proceeded with a straight line starting from a machined initial notch and guided by a side groove thus producing a relaxation of the load. The parallel error between the top surface and the bottom surface of specimens was less than 0.025 mm. The mechanical parameters are listed in Table 3.

5.3.2. Experimental Procedure

The constant displacement load relaxation method only needs to monitor the decrease of the load with time to measure both the stress intensity factor and the subcritical crack growth velocity , regardless of the crack length measurement. The double torsion test was processed using an MTS Insight machine in the testing center of the Central South University. Firstly, the slate specimen should be precracked to obtain reasonable and - data [28, 29], and precrack velocity was 0.05 mm/min. The load-time curve of a typical specimen in precrack is shown in Figure 5. Secondly, the precracked specimen was loaded in 5 mm/min displacement velocity during constant displacement load relaxation testing [19]. Then, the displacement of the slate specimen was kept constant when the load comes up to 95% of the precrack load and the relaxation testing was done [30, 31]. The load-time data was recorded until the load relaxation was stabilized. The load-time curve and the displacement-time curve of a typical specimen are shown in Figures 6 and 7, respectively.

The precracked specimen was loaded in 20 mm/min displacement velocity during the fracture toughness testing, and the critical load was recorded. The fracture toughness can be obtained by Equation (l7). The load-time curve of a typical specimen in testing is shown in Figure 5.

5.3.3. Results and Analysis

(1)KI-V Relations in Transversely Isotropic Slates. Based on the double torsion testing data, the relation between subcritical crack growth velocity and the stress intensity factor can be obtained by Equations (15) and (17). According to the Charles theory [32], to analyze the relation between the subcritical crack growth velocity and the stress intensity factor , a linear function of double logarithmic space can be used, and their - coordinate points in logarithm coordinates are shown in Figures 8–13. - data are regressed as follows: where and are the fitting coefficients.

The power function relation between and is defined as where and are the fitting coefficients.

Table 4 shows each specimen’s subcritical crack growth parameters. Researching rock seepage and stable time effect must be based on the and the relation -, from the perspective of fracture mechanics, and these obtained parameters provide basic information for further study of geotechnical engineering seepage and time-dependent stability.

(2)Fracture Toughness (KIC) and the Threshold Values of Different Bedding Angles of Slates. Based on the testing data, the fracture toughness () can be obtained by Equation (16), which is shown in Table 5.

Typically, during the testing process of the rock subcritical crack growth rate, with the expansion of the crack, the crack tip’s stress intensity factor decreases and the crack extended speed is reduced. When approaching subcritical crack extension threshold value, the crack growth rate is difficult to measure [21]. For rock-like material, the previous researchers did not detect the threshold value, at present. When Wilkins did the research of the subcritical crack expansion of granite, he still did not find the threshold value () with the expansion speed of 10^-12 m/s [17]. In this experiment, the lowest subcritical crack extension speed is in the order of 10^-5 to 10^-6. From the point of view of engineering significance, it can be considered that crack growth has been stopped at this time. In this article, the lowest subcritical crack extension velocity is regarded as the stagnation speed of the crack propagation which is shown in Table 6. According to Equation (19), the stress intensity factor related to the crack propagation can be obtained as the threshold values , which are shown in Table 7.

(3)The Relations between the KIC and the Bedding Angles of Slates. The relation between the bedding angles (β) and the fracture toughness () can be obtained, as shown in Figure 14. It is clear that the basic relation can be expressed as , , from which we can see that fracture toughness decreases then increases with the increase of the β angle and the relation curves show a “U” type; the result can well reflect the change rule of the with the increase of β angles, which is consistent with the theoretical result. Fracture toughness is an important parameter to characterize the rock to resist crack initiation and propagation. This experiment shows that when the bedding angle of the slates is 45°, its fracture toughness is minimal; it means that its capacity to resist crack propagation is minimal, and it is most likely to be damaged during loading.

(4)The Relation between the Threshold Values and the Bedding Angles of Slates. Figure 15 shows the relation between and β angles for different bedding angles’ slates. The basic relationship is , . According to Figure 15, it is found that threshold values () decrease then increase with the increase of the β angle. Because of the different bedding angles corresponding to different thresholds, we can know that the threshold values of the slate specimen have obvious anisotropy. When the bedding angle of the slates is 45°, its threshold value is minimal; it means that when the slate specimen is under load, crack initiation and expansion are most likely to happen, which can cause slate specimen failure.

6. Conclusion

Transversely isotropic elastic parameters and the subcritical crack growth of different bedding angle slate specimens were studied, and a series of uniaxial compression tests and the double torsion constant displacement load relaxation methods were carried out. The following conclusions can be drawn based on this study: (1)Based on the slate specimen elastic parameters obtained from the uniaxial compression test and the double torsion test, it is clear that the slate is significantly transversely isotropic(2)The relations of the mode-I stress intensity factor versus the subcritical crack growth velocity and the fracture toughness were obtained by the double torsion constant displacement load relaxation method; the threshold values were calculated. The calculation results indicate that - relations of different bedding angle slate specimens measured are in accordance with the linear rule, and they are in great agreement with the Charles theory(3)The threshold value of subcritical crack propagation was obtained from the perspective of engineering applications. In this experiment, the range of for these different bedding angle slates is from 0.511 to 0.789, and the lowest subcritical crack extension speed of the different bedding angle slates is in the order of 10^-5 to 10^-6(4)The slate’s bedding angle has an inevitable impact on its mechanical properties. The uniaxial compression, fracture toughness, and threshold values decrease then increases with the increase of the β angle and the relation curve is U-shaped. It is also observed that when the β angle of the slate is 45°, its fracture toughness is minimal. Therefore, crack initiation and crack propagation are generated under loading, resulting in the failure of the slate

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This paper receives funding from a project (No. 51979293) supported by the National Natural Science Foundation of China. The authors wish to acknowledge this support.

References

S. Wang, J. Li, and C. Wang, “Experimental study on subcritical crack propagation in transverse isotropic rocks,” Journal of Central South University(Science and Technology), vol. 48, no. 6, pp. 1607–1614, 2017.
View at: Google Scholar
M. Hakala, H. Kuula, and J. A. Hudson, “Estimating the transversely isotropic elastic intact rock properties for in situ stress measurement data reduction: a case study of the Olkiluoto mica gneiss, Finland,” International Journal of Rock Mechanics and Mining Sciences, vol. 44, no. 1, pp. 14–46, 2007.
View at: Publisher Site | Google Scholar
S. Xie, H. Lin, Y. Chen, R. Yong, W. Xiong, and S. Du, “A damage constitutive model for shear behavior of joints based on determination of the yield point,” International Journal of Rock Mechanics and Mining Sciences, vol. 128, p. 104269, 2020.
View at: Publisher Site | Google Scholar
S.-Q. Yang, P.-F. Yin, B. Li, and D.-S. Yang, “Behavior of transversely isotropic shale observed in triaxial tests and Brazilian disc tests,” International Journal of Rock Mechanics and Mining Sciences, vol. 133, p. 104435, 2020.
View at: Publisher Site | Google Scholar
P. Lade, “Modeling failure in cross-anisotropic frictional materials,” International Journal of Solids and Structures, vol. 44, no. 16, pp. 5146–5162, 2007.
View at: Publisher Site | Google Scholar
H. Lin, S. Xie, R. Yong, Y. Chen, and S. Du, “An empirical statistical constitutive relationship for rock joint shearing considering scale effect,” Comptes Rendus Mécanique, vol. 347, no. 8, pp. 561–575, 2019.
View at: Publisher Site | Google Scholar
X. Zhou, J. Li, and H. Lin, “Analysis of internal stress distribution and mechanics characteristics of pre-existing cavity in brittle rock under triaxial cyclic loading,” Frontiers in Earth Science, vol. 8, 2020.
View at: Publisher Site | Google Scholar
X.-P. Zhou, J.-Z. Zhang, and L. N. Y. Wong, “Experimental study on the growth, coalescence and wrapping behaviors of 3D cross-embedded flaws under uniaxial compression,” Rock Mechanics and Rock Engineering, vol. 51, no. 5, pp. 1379–1400, 2018.
View at: Publisher Site | Google Scholar
R. Cao, R. Yao, J. Meng, Q. Lin, H. Lin, and S. Li, “Failure mechanism of non-persistent jointed rock-like specimens under uniaxial loading: laboratory testing,” International Journal of Rock Mechanics and Mining Sciences, vol. 132, p. 104341, 2020.
View at: Publisher Site | Google Scholar
F. Zhang, X. Wang, and J. Lu, “translators,” in Fracture mechanics, China Architecture & Building Press, Beijing, 1981.
View at: Google Scholar
X. Zhou, H. Cheng, and Y. 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.
View at: Google Scholar
G. Liu and M. Cai, “Modeling time-dependent deformation behavior of brittle rock using grain-based stress corrosion method,” Computers and Geotechnics, vol. 118, p. 103323, 2020.
View at: Publisher Site | Google Scholar
R. Gholami and V. Rasouli, “Mechanical and elastic properties of transversely isotropic slate,” Rock Mechanics and Rock Engineering, vol. 47, no. 5, pp. 1763–1773, 2014.
View at: Publisher Site | Google Scholar
Z. Wang, Z. Zong, L. Qiao, and W. Li, “Elastoplastic model for transversely isotropic rocks,” International Journal of Geomechanics, vol. 18, no. 2, p. 04017149, 2018.
View at: Publisher Site | Google Scholar
B. Amadei, “Importance of anisotropy when estimating and measuring in situ stresses in rock,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 33, no. 3, pp. 293–325, 1996.
View at: Publisher Site | Google Scholar
L. Milou, Rock Mechanics, China Coal Industry Publishing House, Beijing, 1981.
D. Williams and A. Evans, “A simple method for studying slow crack growth,” Journal of Testing and Evaluation, vol. 1, no. 2, pp. 264–270, 1973.
View at: Google Scholar
M. Ciccotti, N. Negri, G. Gonzato, and F. Mulargia, “Practical application of an improved methodology for the double torsion load relaxation method,” International Journal of Rock Mechanics and Mining Sciences, vol. 38, no. 4, pp. 569–576, 2001.
View at: Publisher Site | Google Scholar
M. Saadaoui, P. Reynaud, G. Fantozzi, F. Peronnet, and J. P. Caspar, “Slow crack growth study of plaster using the double torsion method,” Ceramics International, vol. 26, no. 4, pp. 435–439, 2000.
View at: Publisher Site | Google Scholar
Y. Nara and K. Kaneko, “Study of subcritical crack growth in andesite using the double torsion test,” International Journal of Rock Mechanics and Mining Sciences, vol. 42, no. 4, pp. 521–530, 2005.
View at: Publisher Site | Google Scholar
H. Xiao, “Rheologic damage study of three gorges ship lock high slope and analysis on in-situ measured deformation,” Chinese Journal of Rock Mechanics and Engineering, vol. 18, no. 5, pp. 497–502, 1999.
View at: Google Scholar
P. CAO, L. I. Jiang-teng, and H.-p. YUAN, “Testing study of subcritical crack growth rate and fracture toughness in different rocks,” Transactions of Nonferrous Metals Society of China, vol. 16, no. 3, pp. 709–713, 2006.
View at: Publisher Site | Google Scholar
K. Andreas, P. Eckhard, W. Jorg, and B. Wolfgang, “Subcritical crack growth in AI₂O₃ with submicron grain size,” Journal of the European Ceramic Society, vol. 23, no. 1, pp. 81–89, 2003.
View at: Google Scholar
M. Ciccotti, G. Gonzato, and F. Mulargia, “The double torsion loading configuration for fracture propagation: an improved methodology for the load-relaxation at constant displacement,” International Journal of Rock Mechanics and Mining Sciences, vol. 37, no. 7, pp. 1103–1113, 2000.
View at: Publisher Site | Google Scholar
Y. Zhao, L. Zhang, J. Liao, W. Wang, Q. Liu, and L. Tang, “Experimental study of fracture toughness and subcritical crack growth of three rocks under different environments,” International Journal of Geomechanics, vol. 20, no. 8, p. 04020128, 2020.
View at: Publisher Site | Google Scholar
K. Li, Y. Cheng, Z. Y. Yin, D. Han, and J. Meng, “Size effects in a transversely isotropic rock under Brazilian tests: laboratory testing,” Rock Mechanics and Rock Engineering, vol. 53, no. 6, pp. 2623–2642, 2020.
View at: Publisher Site | Google Scholar
A. Evans, “A method for evaluating the time dependent failure characteristics of brittle materials and its applications to polycrystalline alumina,” Journal of Materials Science, vol. 7, no. 11, pp. 37–46, 1972.
View at: Google Scholar
B. Atiunson, “Fracture toughness of Tennessee sandstone and Carrara marble using the double torsion testing method,” International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, vol. 16, no. 1, pp. 49–53, 1979.
View at: Publisher Site | Google Scholar
X.-P. Zhou, Y.-T. Wang, J.-Z. Zhang, and F.-N. Liu, “Fracturing behavior study of three-flawed specimens by uniaxial compression and 3D digital image correlation: sensitivity to brittleness,” Rock Mechanics and Rock Engineering, vol. 52, no. 3, pp. 691–718, 2019.
View at: Publisher Site | Google Scholar
M. Ciccotti, N. Negri, L. Sassi, G. Gonzato, and F. Mulargia, “Elastic and fracture parameters of Etna, Stromboli, and Vulcano lava rocks,” Journal of Volcanology and Geothermal Research, vol. 98, no. 1-4, pp. 209–217, 2000.
View at: Publisher Site | Google Scholar
M. E. Ebrahimi, J. Chevalier, and G. Fantozzi, “Slow crack-growth behavior of alumina ceramics,” Journal of Materials Research, vol. 15, no. 1, pp. 142–147, 2000.
View at: Publisher Site | Google Scholar
E. Z. Lajtai and L. P. Bielus, “Stress corrosion cracking of Lac du Bonnet granite in tension and compression,” Rock Mechanics and Rock Engineering, vol. 19, no. 2, pp. 71–87, 1986.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Gang Ma 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.

PDF Download Citation

Download other formats

Order printed copies

Views

343

Downloads

837

Citations