MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2016/4517360 4517360 Research Article Numerical Analysis of Flattened Brazilian Disc Test Based on the Cusp Catastrophe Theory http://orcid.org/0000-0002-4282-9963 Wang Min 1,2 2 Cao Ping 1 He Xiao-Qiao 1 School of Resources and Safety Engineering Central South University Changsha Hunan 410083 China csu.edu.cn 2 Hunan Provincial Key Laboratory of Shale Gas Resource Utilization Hunan University of Science and Technology Xiangtan Hunan 411201 China hnust.edu.cn 2016 3172016 2016 20 04 2016 11 06 2016 26 06 2016 2016 Copyright © 2016 Min Wang and Ping Cao. 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.

The Brazilian disc test is a simple and useful technique to determine the tensile strength of rock materials. By using FLAC3D, 63 numerical simulations in total were performed when flattened Brazilian disc coefficient and Poisson’s ratio were different. Based on Griffith theory, the corresponding FISH language was compiled to record the Griffith equivalent stress. Through analysis of numerical simulation results, it is indicated that fracture plane was not the plane going through center of the Brazilian disc, which was in good agreement with the references. In addition, the flattened Brazilian disc coefficients had greater influence on tensile strength than Poisson’s ratio. Based on cusp catastrophe theory, the flattened Brazilian disc coefficient should not exceed 0.035 for the flattened Brazilian disc tests. Consequently, a tensile strength empirical formula considering flattened Brazilian disc coefficient by utilizing the flattened Brazilian disc test was established, which was σ t = 0.9993 e x p ( - 11.65 ε ) 2 p / π D t , ε 0.035 .

Hunan University of Science and Technology E21527 National Natural Science Foundation of China 51174088 51174228 National Basic Research Program of China 013CB035401 Fundamental Research Funds for the Central Universities of Central South University 2015zzts077 2014zzts055
1. Introduction

The Brazilian disc test is a useful technique to determine the tensile strength of rock materials ; the Brazilian disc test has attracted the attention of a large number of scholars due to its extensive practical importance in engineering application. There are many studies on the Brazilian disc test in theory analysis, laboratory test, and numerical simulations. Jianhong et al.  analyzed the plane stress distribution of the Brazilian disc test; the tensile elastic modulus of a rock can be obtained by the results of theory analysis. Serati et al.  investigated the double Fourier series technique to construct an elastic stress filed in a cylindrical subject to lateral boundary loads, and this method was verified against available elastic solutions for axisymmetric and nonaxisymmetric engineering problems such as Brazilian tensile strength and point load strength. With regard to the experiments for the Brazilian test, Tan et al.  conducted a series of Brazilian tests on Monsel slate considering different foliation-loading angles, fracture patterns and strength of samples were analyzed, and how the microparameters influence the bearing capacity and failure modes of Brazilian disc tests for anisotropic rocks was revealed. Liu et al.  investigated the mechanical tensile behavior of Beisha granite with different test methods, especially the damage evolution process. It was revealed that, in direct tensile testing, the AE events accumulated mainly along the failure surface, whereas in indirect tests the location and variation of AE events were more complex. In the field of numerical simulations, Riera et al.  utilize the discrete element method (DEM) to simulate the Brazilian test, and it was concluded that specifications for the splitting test of quasi-fragile materials should be examined under the light of recent advances in fracture crack propagation and dynamic rupture, properly considering the influence of the not-yet-clearly-defined material properties and geometry of the strip. Yang and Huang  constructed a model for Brazilian disc splitting test using PFC2D, and then it was used to simulate Brazilian splitting test for jointed rock mass specimens and specimen containing a central straight notch, and the influences of joint angles on tensile strength and failure modes were analyzed.

However, three-dimension effects of Brazilian test can influence the results , and hence the thickness ( t )/diameter ( D ) ratio should be considered in the Brazilian disc. ASTM suggested that the thickness ( t )/diameter ( D ) ratio ranged between 0.2 and 0.75 , while ISRM recommended that t / D ratio should be 0.5 . In addition, the deformation on loading location would occur, and the loading mode and stress distribution also changed, which had given rise to the large error by utilizing Brazilian tests to obtain the tensile strength of rock materials. Obviously, the deformation on loading location is inevitable. Hence, reducing deformation on loading location would be the key to improving the results accuracy utilizing the Brazilian test, and the flattened Brazilian test would be a good alternative.

The flattened Brazilian test has been a popular method in recent years to determine the tensile strength of rock materials . Wang et al.  conducted theory analysis, and the stress distribution for flattened Brazilian test was obtained; moreover, the numerical simulations of flattened Brazilian test were performed, and the formulas for elastic modulus, tensile strength, and fracture toughness of rocks were derived. Lin et al.  analyzed the three-dimensional effects in the flattened Brazilian test through numerical simulation, and the empirical formula of the flattened Brazilian test was proposed.

In this paper, the flattened Brazilian disc coefficient was used to measure the size of the flattened Brazilian disc plane instead of the loading angle. The flattened Brazilian disc with different flattened Brazilian disc coefficient was constructed in FLAC3D to simulate the flattened Brazilian tests. The FISH language was compiled to illustrate the Griffith equivalent stress of the flattened Brazilian discs, and the possible fracture planes of the flattened Brazilian discs were analyzed. Consequently, the empirical formula of tensile strength using flattened Brazilian tests was put forward.

2. Numerical Simulation of Flattened Brazilian Disc Test 2.1. Flattened Brazilian Disc Coefficient <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M8"> <mml:mrow> <mml:mi>ε</mml:mi></mml:mrow> </mml:math></inline-formula>

In most references, the loading angle was used to describe the size of loading plane; however, the angle is much complicated to measure in practice. Hence, in this paper, the flattened Brazilian disc coefficient ε was employed to depict the size of loading plane, displayed in Figure 1.

The flattened Brazilian disc coefficient.

The flattened Brazilian disc coefficient can be expressed as the following equation: (1) ε = B C O B = 1 - D C A B , where ε is the flattened Brazilian disc coefficient; D C is the height of the flattened Brazilian disc; and A B is the diameter of the disc. Based on (1), it is obvious that the flattened Brazilian disc coefficient can be measured more easily than the loading angle.

2.2. Numerical Simulations Construction and Mechanical Parameters

To study the influence of the flattened Brazilian disc coefficient on the flattened Brazilian disc test, the commercial software FLAC3D was utilized. To construct the numerical simulation model, its modeling and meshing were accomplished in ANSYS with Solid185 zone type; subsequently, the constructed numerical simulation models were introduced into FLAC3D, the numerical simulation models with different flattened Brazilian disc coefficient were constructed (Figure 2), the flattened Brazilian disc coefficient ranged from 0 to 0.04 with spacing of 0.005, diameter of Brazilian discs was 0.05 m, and their thickness was 0.05 m.

Some flattened Brazilian discs with different flattened Brazilian disc coefficient.

0.01

0.02

0.03

0.04

The uniaxial compression loading was applied on the flattened plane of the flattened Brazilian disc directly, and its loading rate was 5 e - 9  m/step. Simultaneously, the related FISH program  was compiled to record the applied load p . When the mechanical ratio fell below 1 e - 5 , then numerical simulations ended. The numerical simulation model was elastic model. Due to having no influence of Young’s modulus E on the stress distribution [23, 24], Young’s modulus for all numerical simulation models was 100 GPa. Poisson’s ratio included 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, and 0.4, and shear modulus G of numerical simulation models can be calculated based on the following equation: (2) G = E 2 1 + υ , where E is Young’s modulus and υ is Poisson’s ratio.

Through combination of different flattened Brazilian disc coefficient (9 types) and Poisson’s ratio (7 types), 63 numerical simulations in total were performed.

3. Numerical Simulation Results Analysis 3.1. Fracture Plane Analysis of the Brazilian Disc Test

For most rock materials, they belong to brittle materials, and the fracture mechanism can be explained by Griffith theory [25, 26]. Moreover, Yu and Xu  suggested that the location of the maximum equivalent stress σ G was the fracture plane location, and the key of finding the fracture plane was determining location of the maximum equivalent stress. Based on Griffith criterion, the Griffith equivalent stress σ G can be expressed as follows: (3) σ G = - σ 1 - σ 3 2 σ 1 + σ 3 , σ 1 + 3 σ 3 0 , σ 3 , σ 1 + 3 σ 3 < 0 , where σ 1 is the maximum principal stress, σ 3 is the minimum principal stress, and σ G is the Griffith equivalent stress.

In order to explore the location of fracture plane of the Brazilian disc test, the corresponding FISH language was compiled to calculate the average Griffith equivalent stress σ G P of the selected plane, which was parallel to y z -plane and with x -coordinate, showed in Figure 3.

A selected plane parallel to y z -plane and with x -coordinate.

To determine the possible fracture plane during loading process, four numerical simulation results were taken as examples to illustrate the possible fracture plane. For convenient comparison of the numerical simulation results, contour of the dimensionless ratio σ G / ( 2 p / π D t ) was displayed in FLAC3D by compiling the FISH language, and the dimensionless ratio σ G P / ( 2 p / π D t ) of the selected plane with different x -coordinate is displayed in Figure 4.

Contour of the dimensionless ratio σ G / ( 2 p / π D t ) and ratio σ G P / ( 2 p / π D t ) with different x -coordinate when Poisson’s ratio υ is 0.1.

Flattened Brazilian disc coefficient is 0.01

Flattened Brazilian disc coefficient is 0.02

Flattened Brazilian disc coefficient is 0.03

Flattened Brazilian disc coefficient is 0.04

As shown in Figure 4, the maximum dimensionless ratio σ G P / ( 2 p / π D t ) of the selected plane was not in the y z -plane, which meant that the possible fracture plane was not the plane going through the center of Brazilian discs. For instance, when the flattened Brazilian disc coefficient was 0.04 (Figure 4(d)), the maximum ratio σ G P / ( 2 p / π D t ) of the selected plane was not the plane going through the center of Brazilian disc ( x -coordinate of the selected plane was 0 mm, and σ G P / ( 2 p / π D t ) was 0.4545), but the plane going through the compression plane ends ( x -coordinate of the selected plane was −4 mm, and ratio σ G P / ( 2 p / π D t ) was 0.4721). The larger the ratio σ G P / ( 2 p / π D t ) was, the larger the average Griffith equivalent stress σ G P was; consequently, the possible fracture plane would be the plane with −4 mm x -coordinate when the flattened Brazilian disc coefficient was 0.04 and Poisson’s ratio was 0.1. Based on numerical simulation results, it was implied that the possible fracture plane possibly occurred on the plane going through the compression plane ends ( A , B and A , B ), displayed in Figure 5.

Possible fracture plane.

The numerical simulation results indicated that the location of the maximum Griffith equivalent stress is not the center of the Brazilian disc; that is, the cracks would not initiate from the center of the flattened Brazilian discs, which cannot guarantee the accuracy of the tensile strength by utilizing the Brazilian discs test. Through analysis of the numerical simulation results, it is suggested that the accuracy of the tensile strength difference decreases with increasing of flattened Brazilian disc coefficient.

Moreover, these results were also verified by the study by You and Su . It is observed that the fracture plane is not going through the center of Brazilian discs, and these results were in good agreement with the observation of [29, 30]. Changing of the flattened Brazilian disc coefficient influenced stress distribution of the flattened Brazilian disc, which had given rise to the fracture plane away from the center of Brazilian discs.

Both numerical simulation results and laboratory tests  had shown that fracture plane location would not occur along the plane going through the center of Brazilian disc, but it had more possibly gone through the compression plane ends. Thus, when adopting the Brazilian disc test to determine the tensile strength, taking the Griffith equivalent stress of Brazilian disc center as the tensile strength of rock materials would not be so accurate in the flattened Brazilian disc tests. Hence, it was necessary to calibrate the empirical formula of the tensile strength by adopting the flattened Brazilian disc test.

3.2. Empirical Formula of Flattened Brazilian Disc Tests

In accordance with the aforementioned numerical simulations results, the fracture plane deviated from the center of Brazilian discs; hence, taking the Griffith equivalent stress of Brazilian discs’ center as the tensile strength of rock materials was not reasonable. However, taking the maximum Griffith equivalent stress σ G - m a x of the flattened Brazilian disc would be more accurate.

As discussed previously, the stress distribution changed with increasing the flattened Brazilian disc coefficient, and the possible fracture plane was not the plane going through the Brazilian discs center; in other words, the maximum Griffith equivalent stress location changed as well, and the correlation of the tensile strength and the applied load p should be calibrated. The correlation between the maximum Griffith equivalent stress σ G - m a x of the center plane and the applied load p can be expressed as follows: (4) σ t = σ G - m a x = k 2 p π D t , where σ t is the tensile strength of rock materials by using the flattened Brazilian disc test, σ G - m a x is the maximum Griffith equivalent stress of the Brazilian disc, p is the applied load, D is the diameter of Brazilian discs, t is the thickness of Brazilian discs, and k is the correction coefficient.

In order to explore influence of Poisson’s ratio and the flattened Brazilian disc coefficient on the correction coefficient k , its changing trend with increasing Poisson’s ratio and the Brazilian disc coefficient are displayed in Figure 6, respectively.

Influence of Poisson’s ratio and flattened Brazilian disc coefficient on the correction coefficient k .

As shown in Figure 6, the influence of Poisson’s ratio on the correction coefficient k was not obvious (in Figure 6(a)). When the flattened Brazilian disc coefficient was constant, the correction coefficient k remained at a stable level with variation of Poisson’s ratio. For example, the flattened Brazilian disc coefficient was 0.02, and the correction coefficients k were 0.7467, 0.7499, 0.7731, 0.7833, 0.7903, 0.7638, and 0.7765 when Poisson’s ratios were 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, and 0.4, separately. However, the flattened Brazilian disc coefficient influenced the correction coefficient k a lot (shown in Figure 6(b)), and k noticeably changed with increasing the flattened Brazilian disc coefficient.

Through analysis of Figure 6(b), the tensile strength was influenced strongly by the flattened Brazilian disc coefficient; therefore, it was necessary to determine the maximum flattened Brazilian disc coefficient for avoiding big error when using (4). In order to determine the maximum flattened Brazilian disc coefficient, cusp catastrophe theory  was adopted. Because Poisson’s ratio had little influence on k , Poisson’s ratio was neglected in analysis of determining the maximum flattened Brazilian disc coefficient. Subsequently, the average error e t a of the tensile strength with different flattened Brazilian disc coefficient was introduced to express the error of adopting (4) in the Brazilian disc test, and the average error e t a can be defined as follows: (5) e t a = i = 1 n σ G - m a x i - 2 p i / π D t / 2 p i / π D t n = i = 1 n k i - 1 n , where e t a is the average error of tensile strength with different flattened Brazilian disc coefficient, p i is the maximum applied load, D is the diameter of Brazilian discs, t is the thickness of Brazilian discs, k i is the correction coefficient, and n is the number of Poisson’s ratio types ( n is 7 in this paper).

Through analysis of Figure 6 and (5), it is found that the calculation error would increase with increasing of the flattened Brazilian disc coefficient. The large flattened Brazilian disc coefficient would result in the crack initiation location far away from the flattened Brazilian disc center; thereafter, a large error would occur when using the flattened Brazilian disc test to determine the tensile strength of rock materials. In this paper, the cusp catastrophe theory was utilized to determine the maximum flattened Brazilian disc coefficient. The correlation of the error e t a and the flattened Brazilian disc coefficient ε can be expressed as the cusp catastrophe function e t a ( ε ) . According to cusp catastrophe theory, the function e t a ( ε ) can be written as follows: (6) e t a ε = ε 4 + u ε 2 + v ε .

The potential function e t a ( ε ) can be plotted as shown in Figure 7.

Cusp catastrophe function e t a ( ε ) .

As shown in Figure 7, cusp catastrophe function e t a ( ε ) consists of three parts: up lobe, middle lobe, and down lobe. Based on cusp catastrophe theory, the characteristics of cusp catastrophe function e t a ( ε ) were controlled by u and v , which formed the control variable plane. The stability of cusp catastrophe function e t a ( ε ) was controlled by the bifurcation set Δ = 8 u 3 + 27 v 2 . If 8 u 3 + 27 v 2 0 , the function e t a ( ε ) was stable. While 8 u 3 + 27 v 2 < 0 , the function e t a ( ε ) was unstable. In other words, if 8 u 3 + 27 v 2 < 0 , (6) describing the correlation of the tensile strength σ t and the correction coefficient k was not stable, and catastrophe change would occur. Thus, a big error would occur when using (6) to determine the tensile strength of rock materials in the flattened Brazilian disc test.

To determine the maximum flattened Brazilian disc coefficient, a series of data sets { ( e t a ( ε ) 1 , ε 1 ) , ( e t a ( ε ) 2 , ε 2 ) , , ( e t a ( ε ) i , ε i ) } ( i = 4,5 , , 9 ) were formed firstly. For example, when the maximum flattened Brazilian disc coefficient was 0.015, the corresponding data set would be { ( e t a ( ε ) 1 , ε 1 ) , ( e t a ( ε ) 2 , ε 2 ) , ( e t a ( ε ) 3 , ε 3 ) , ( e t a ( ε ) 4 , ε 4 ) } , which meant the flattened Brazilian disc coefficients ε 1 , ε 2 , ε 3 , and ε 4 were 0, 0.005, 0.01, and 0.015, respectively. Then, the least square fittings of the data sets { ( e t a ( ε ) 1 , ε 1 ) , ( e t a ( ε ) 2 , ε 2 ) , , ( e t a ( ε ) i , ε i ) } ( i = 4,5 , , 9 ) were conducted adopting (4). Subsequently, different groups of u and v were obtained when the maximum flattened Brazilian disc coefficient was different, and the corresponding bifurcation sets Δ = 8 u 3 + 27 v 2 were calculated. When Δ = 8 u 3 + 27 v 2 < 0 , then the catastrophe change happened, and then the maximum flattened Brazilian disc coefficient was determined.

Based on the method above, a series of u , v , and 8 u 3 + 27 v 2 were calculated when the maximum flattened Brazilian disc coefficient was different (Table 1), and Figure 8 displays the variation of Δ = 8 u 3 + 27 v 2 when i is 8 and 9.

Different u , v , and 8 u 3 + 27 v 2 when maximum flattened Brazilian disc coefficients are different.

i 4 5 6 7 8 9
Maximum flattened Brazilian disc coefficient 0.015 0.020 0.025 0.030 0.035 0.040
u 0.548 141.100 23.990 125.600 48.030 −16.580
v 8.834 7.2100 8.922 7.091 8.697 10.230
Δ = 8 u 3 + 27 v 2 >0 >0 >0 >0 >0 <0

8 u 3 + 27 v 2 changing trend when maximum flattened Brazilian disc coefficient is 0.035 and 0.04.

Through analysis of Figure 8, Δ = 8 u 3 + 27 v 2 was positive when i was 8. While i was 9, Δ = 8 u 3 + 27 v 2 was negative; that is to say, cusp catastrophe function e t a ( ε ) for data set { ( e t a ( ε ) 1 , ε 1 ) , ( e t a ( ε ) 2 , ε 2 ) , , ( e t a ( ε ) 8 , ε 8 ) } was stable, which suggested that catastrophe change of error will not occur when using (4) to determine tensile strength. While data set { ( e t a ( ε ) 1 , ε 1 ) , ( e t a ( ε ) 2 , ε 2 ) , , ( e t a ( ε ) 9 , ε 9 ) } was not stable, the catastrophe change of error occurred; hence, the data ( e t a ( 0.04 ) , 0.04 ) caused the catastrophe change, and the flattened Brazilian disc coefficient should be less than 0.04. However, for accuracy of the tensile strength when using (4), flattened Brazilian disc coefficient should not exceed 0.035, which was in good agreement with the study by Lin et al. .

Because of little influence of Poisson’s ratio on the correction coefficient k , the tensile strength determined by (4) without consideration of Poisson’s ratio was acceptable. To establish correlation of flattened Brazilian disc coefficient ε and the correction coefficient k , the data was used when flattened Brazilian disc coefficient ranged from 0 to 0.035, and the average correction coefficient k a was calculated when flattened Brazilian disc coefficients were different, listed in Table 2.

Average correction coefficient k a when flattened Brazilian disc coefficients are different.

Correction coefficient (k ) Flattened Brazilian disc coefficient ( ε )
0 0.005 0.010 0.015 0.020 0.025 0.030 0.035
Poisson’s ratio ( υ )
0.10 0.9398 0.9188 0.9057 0.8288 0.7772 0.7467 0.6393 0.6393
0.15 0.9404 0.9473 0.9080 0.8398 0.7939 0.7499 0.6537 0.6537
0.20 0.9402 0.9350 0.9231 0.8391 0.7879 0.7731 0.6461 0.6461
0.25 0.9428 0.9650 0.9316 0.8607 0.8086 0.7833 0.6723 0.6723
0.30 0.9586 0.9558 0.9231 0.8750 0.8126 0.7903 0.6596 0.6596
0.35 0.9712 0.9685 0.9256 0.8977 0.7934 0.7638 0.6823 0.6823
0.40 1.0002 0.9911 0.9440 0.8813 0.8014 0.7765 0.6874 0.6874

Average correction coefficient ( k a ) 0.9562 0.9542 0.9230 0.8603 0.7962 0.7691 0.6630 0.6545

Based on the data in Table 2, correlation of k a and flattened Brazilian disc coefficient ε was constructed via the least square fitting method, shown in Figure 9.

Correlation of the average correction coefficient and flattened Brazilian disc coefficient.

Consequently, correlation of the average correction coefficient and flattened Brazilian correction coefficient can be expressed as follows: (7) k a = 0.9993 e x p - 11.65 ε .

Due to little influence of Poisson’s ratio on the correction coefficient k , k in (4) can be replaced by k a ; combining (4) and (7), the tensile strength formulation of the Brazilian disc test can be denoted as follows: (8) σ t = 0.9993 e x p - 11.65 ε · 2 p π D t , where ε is the flattened Brazilian disc coefficient ( ε 0.035 ), σ t is the tensile strength, p is the applied load, D is the diameter of Brazilian disc, and t is the thickness of Brazilian disc.

Based on the numerical simulation results, the tensile strength adopting the flattened Brazilian disc test with variation flattened Brazilian disc coefficient is proposed. Moreover, flattened Brazilian disc coefficient is easily obtained during the Brazilian disc test compared to loading angle; hence, the empirical tensile strength equation is workable, and it could be adopted to estimate the tensile strength in the flattened Brazilian disc test.

4. Discussion

The tensile strength is a key mechanical parameter for rock materials, which plays a dominant role in rock engineering. Both direct method and indirect method were used to determine the tensile strength of rock materials; however, the direct method is time-consuming and hard-conducting. Thereafter, the indirect method was more widely used to obtain the tensile strength of rock materials compared to the direct method.

The Brazilian disc test is one of the most widely used indirect methods, whereas some scholars found that the loading location was smashed and the crack did not initiate from the center of the Brazilian disc due to the stress concentration in the loading location, which decrease the accuracy of the tensile strength by using the Brazilian disc test. It is because the most important point of the Brazilian disc test is the crack initiation from the center of the Brazilian disc center. To increase the accuracy of the tensile strength of the Brazilian disc test, some new Brazilian disc tests were proposed, and the flattened Brazilian disc test was one of them.

In order to obtain accurate tensile strength of rock materials by utilizing the flattened Brazilian disc test, the shape parameter of the flattened Brazilian disc was defined as the flattened Brazilian disc coefficient. Then, 63 numerical simulations in total were performed, and the related FISH language was applied to obtain the Griffith equivalent stress based on the Griffith theory. Through analysis of the numerical simulation results, it is found that the crack initiation location (the maximum Griffith equivalent stress location) deviates from the center of the flattened Brazilian disc center. Moreover, the corresponding tensile strength difference becomes larger with increasing of the flattened Brazilian disc coefficient. Thus, it is necessary to determine the appropriate value of the flattened Brazilian disc coefficient to avoid large error between the experimental results and the true value. In this paper, the cusp catastrophe theory was adopted to determine the maximum flattened Brazilian disc coefficient. Based on the cusp catastrophe theory, the flattened Brazilian disc coefficient should be less than 0.035. Additionally, an empirical formula was proposed evaluating the tensile strength of rock materials by using the flattened Brazilian disc test.

The numerical simulation is mainly used in this paper, the numerical simulation result was not verified by the experimental results, and it lacks further experimental study. To validate the numerical simulation results, the experiments study on the flattened Brazilian disc test would be our next tasks.

5. Conclusions

In this paper, the flattened Brazilian disc tests were conducted, through numerical simulations results analysis, the fracture plane locations of flattened Brazilian disc were determined, and, finally, the tensile strength empirical equation considering flattened Brazilian disc coefficient was obtained; the main conclusions of this paper were as follows:

Through numerical simulations analysis, it was suggested that the fracture plane location was not the plane going through the center of Brazilian discs with increasing flattened Brazilian disc coefficient; it was the plane going through the compression plane ends.

Combining numerical simulation results, an equation for the tensile strength adopting flattened Brazilian disc test was established. Because the maximum Griffith equivalent stress location was not the Brazilian disc center any longer, which was verified by both laboratory tests and numerical simulations, it is necessary to calibrate the empirical formula. The empirical formula for the tensile strength determined by the Brazilian disc test was established, and the empirical equation was σ t = 0.9993 e x p ( - 11.65 ε ) 2 p / π D t , ε 0.035 .

Competing Interests

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

Acknowledgments

This study was funded by the Open Research Fund Program of Hunan Provincial Key Laboratory of Shale Gas Resource Utilization; Hunan University of Science and Technology (Grant no. E21527); the National Natural Science Foundation of China (Grants nos. 51174088, 51174228); the National Basic Research Program of China (Grant no. 013CB035401); and the Fundamental Research Funds for the Central Universities of Central South University (Grants nos. 2015zzts077, 2014zzts055).

Komurlu E. Kesimal A. Evaluation of indirect tensile strength of rocks using different types of jaws Rock Mechanics and Rock Engineering 2015 48 4 1723 1730 10.1007/s00603-014-0644-3 2-s2.0-84932642951 Razavi S. M. Gonzalez M. Cuitiño A. M. General and mechanistic optimal relationships for tensile strength of doubly convex tablets under diametrical compression International Journal of Pharmaceutics 2015 484 1-2 29 37 10.1016/j.ijpharm.2015.02.030 2-s2.0-84923586798 Karaman K. Cihangir F. Ercikdi B. Kesimal A. Demirel S. Utilization of the Brazilian test for estimating the Uniaxial compressive strength and shear strength parameters Journal of the Southern African Institute of Mining and Metallurgy 2015 115 3 185 192 10.17159/2411-9717/2015/v115n3a3 2-s2.0-84937153994 Khanlari G. Rafiei B. Abdilor Y. An experimental investigation of the Brazilian tensile strength and failure patterns of laminated sandstones Rock Mechanics and Rock Engineering 2014 48 2 843 852 10.1007/s00603-014-0576-y 2-s2.0-84939886419 Mikl-Resch M. J. Antretter T. Gimpel M. Kargl H. Pittino G. Tichy R. Ecker W. Galler R. Numerical calibration of a yield limit function for rock materials by means of the Brazilian test and the uniaxial compression test International Journal of Rock Mechanics and Mining Sciences 2015 74 24 29 10.1016/j.ijrmms.2014.12.001 Riera J. D. Miguel L. F. F. Iturrioz I. Assessment of Brazilian tensile test by means of the truss-like Discrete Element Method (DEM) with imperfect mesh Engineering Structures 2014 81 10 21 10.1016/j.engstruct.2014.09.036 2-s2.0-84908010471 Jianhong Y. Wu F. Q. Sun J. Z. Estimation of the tensile elastic modulus using Brazilian disc by applying diametrically opposed concentrated loads International Journal of Rock Mechanics and Mining Sciences 2009 46 3 568 576 10.1016/j.ijrmms.2008.08.004 2-s2.0-60749133488 Serati M. Alehossein H. Williams D. J. 3D elastic solutions for laterally loaded discs: generalised Brazilian and point load tests Rock Mechanics and Rock Engineering 2014 47 4 1087 1101 10.1007/s00603-013-0449-9 2-s2.0-84903533027 Tan X. Konietzky H. Frühwirt T. Dan D. Q. Brazilian tests on transversely isotropic rocks: laboratory testing and numerical simulations Rock Mechanics and Rock Engineering 2014 48 4 1341 1351 10.1007/s00603-014-0629-2 2-s2.0-84928813753 Liu J. F. Chen L. Wang C. P. Man K. Wang L. Wang J. Su R. Characterizing the mechanical tensile behavior of Beishan granite with different experimental methods International Journal of Rock Mechanics and Mining Sciences 2014 69 50 58 10.1016/j.ijrmms.2014.03.007 2-s2.0-84898876499 Yang S.-Q. Huang Y.-H. Particle flow study on strength and meso-mechanism of Brazilian splitting test for jointed rock mass Acta Mechanica Sinica 2014 30 4 547 558 10.1007/s10409-014-0076-z 2-s2.0-84906269370 Yu Y. Yin J. Zhong Z. Shape effects in the Brazilian tensile strength test and a 3D FEM correction International Journal of Rock Mechanics and Mining Sciences 2006 43 4 623 627 10.1016/j.ijrmms.2005.09.005 2-s2.0-33644627567 Standard test method for splitting tensile strength of cylindrical concrete specimens 1984 ASTM C496 Philadelphia, Pa, USA ASTM ISRM Suggested methods for determining tensile strength of rock materials International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts 1978 15 3 99 103 10.1016/0148-9062(78)90003-7 Zhang Q. B. Zhao J. A review of dynamic experimental techniques and mechanical behaviour of rock materials Rock Mechanics and Rock Engineering 2014 47 4 1411 1478 10.1007/s00603-013-0463-y 2-s2.0-84903515348 Chen R. Dai F. Qin J. Lu F. Flattened Brazilian disc method for determining the dynamic tensile stress-strain curve of low strength brittle solids Experimental Mechanics 2013 53 7 1153 1159 10.1007/s11340-013-9733-6 2-s2.0-84881542344 Wang Q. Z. Feng F. Ni M. Gou X. P. Measurement of mode I and mode II rock dynamic fracture toughness with cracked straight through flattened Brazilian disc impacted by split Hopkinson pressure bar Engineering Fracture Mechanics 2011 78 12 2455 2469 10.1016/j.engfracmech.2011.06.004 2-s2.0-79960311323 Keles C. Tutluoglu L. Investigation of proper specimen geometry for mode I fracture toughness testing with flattened Brazilian disc method International Journal of Fracture 2011 169 1 61 75 10.1007/s10704-011-9584-z 2-s2.0-79953187898 Wang Q.-Z. Xing L. Determination of fracture toughness KIC by using the flattened Brazilian disk specimen for rocks Engineering Fracture Mechanics 1999 64 2 193 201 10.1016/S0013-7944(99)00065-X 2-s2.0-0033333865 Wang Q. Z. Jia X. M. Kou S. Q. Zhang Z. X. Lindqvist P.-A. The flattened Brazilian disc specimen used for testing elastic modulus, tensile strength and fracture toughness of brittle rocks: analytical and numerical results International Journal of Rock Mechanics and Mining Sciences 2004 41 2 245 253 10.1016/s1365-1609(03)00093-5 2-s2.0-0742287770 Lin H. Xiong W. Xiong Z. Y. Gong F. Three-dimensional effects in a flattened Brazilian disk test International Journal of Rock Mechanics and Mining Sciences 2015 74 10 14 10.1016/j.ijrmms.2014.11.006 2-s2.0-84921038386 Itasca Consulting Group Fast Lagrangian Analysis of Continua in 3 Dimensions, User Manual, Version 3.1 2004 Minneapolis, Minn, USA Itasca Consulting Group Yu Y. Questioning the validity of the Brazilian test for determining tensile strength of rocks Chinese Journal of Rock Mechanics and Engineering 2005 24 7 1150 1157 2-s2.0-17944370287 Xu Z. L. Elasticity Beijing, China Higher Education Press 2006 (Chinese) Jaeger J. C. Cook N. G. W. Zimmerman R. W. Fundamentals of Rock Mechanics 2007 4th Oxford, UK Blackwell Griffith A. A. The phenomena of flow and rupture in solids Philosophical Transactions of the Royal Society of London A 1920 221 163 198 Yu Y. Xu Y. L. Method to determine tensile strength of rock using flattened Brazilian disk Chinese Journal of Rock Mechanics and Engineering 2006 25 7 1457 1462 2-s2.0-33748445326 You M. Q. Su C. D. Experimental study on split test with flatten disk and tensile strength of rock Chinese Journal of Rock Mechanics and Engineering 2004 23 18 3106 3112 Fairhurst C. On the validity of the ‘Brzilian’ test for brittle materials International Journal of Rock Mechanics and Mining Sciences 2009 46 421 425 Swab J. J. Yu J. Gamble R. Kilczewski S. Analysis of the diametral compression method for determining the tensile strength of transparent magnesium aluminate spinel International Journal of Fracture 2011 172 2 187 192 10.1007/s10704-011-9655-1 2-s2.0-84905505586 Liu Z.-X. Lan M. Xiao S.-Y. Guo H.-Q. Damage failure of cemented backfill and its reasonable match with rock mass Transactions of Nonferrous Metals Society of China 2015 25 3 954 959 63684 10.1016/s1003-6326(15)63684-6 2-s2.0-84927510554 Zhang C.-P. Han K.-H. Fang Q. Zhang D.-L. Functional catastrophe analysis of collapse mechanisms for deep tunnels based on the Hoek-Brown failure criterion Journal of Zhejiang University: Science A 2014 15 9 723 731 10.1631/jzus.A1400014 2-s2.0-84919910643 Pei-Jian J. En-Yuan W. Ning H. Si-Heng W. Catastrophe progression method on forecast of rock burst Disaster Advances 2013 6 7 34 38 2-s2.0-84883207448 Pan Y. Zhao Z.-G. Analysis of main shock of thrust fault earthquake by catastrophe theory Applied Mathematics and Mechanics-English Edition 2012 33 7 845 864 10.1007/s10483-012-1590-8 MR2984802 2-s2.0-84864553804 Yao W.-J. Qiu Y.-Z. Chen Z.-K. Ultimate bearing capacity analysis of a super-long rock-socketed filling pile based on catastrophe theory Advances in Structural Engineering 2010 13 2 331 338 10.1260/1369-4332.13.2.331 2-s2.0-77953353057 Pan Y. Li A.-W. Fold catastrophe model of strike-slip fault earthquake Applied Mathematics and Mechanics-English Edition 2010 31 3 349 362 10.1007/s10483-010-0308-9 MR2655470 2-s2.0-77951672970 Wang S. Y. Lam K. C. Au S. K. Tang C. A. Zhu W. C. Yang T. H. Analytical and numerical study on the pillar rockbursts mechanism Rock Mechanics and Rock Engineering 2006 39 5 445 467 10.1007/s00603-005-0075-2 2-s2.0-33750520952