The Scientific World Journal

Volume 2014 (2014), Article ID 735497, 7 pages

http://dx.doi.org/10.1155/2014/735497

## Numerical Built-In Method for the Nonlinear JRC/JCS Model in Rock Joint

Institute of Mineral Resources, Chinese Academy of Geological Sciences, Beijing 100037, China

Received 15 November 2013; Accepted 1 January 2014; Published 11 February 2014

Academic Editors: A. Esmaeily and A. Rodríguez-Castellanos

Copyright © 2014 Qunyi 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

The joint surface is widely distributed in the rock, thus leading to the nonlinear characteristics of rock mass strength and limiting the effectiveness of the linear model in reflecting characteristics. The JRC/JCS model is the nonlinear failure criterion and generally believed to describe the characteristics of joints better than other models. In order to develop the numerical program for JRC/JCS model, this paper established the relationship between the parameters of the JRC/JCS and Mohr-Coulomb models. Thereafter, the numerical implement method and implementation process of the JRC/JCS model were discussed and the reliability of the numerical method was verified by the shear tests of jointed rock mass. Finally, the effect of the JRC/JCS model parameters on the shear strength of the joint was analyzed.

#### 1. Introduction

The Mohr-Coulomb linear model is generally used to study the characteristics of rock mass strength [1–3]. Lin et al. [1] used numerical algorithm to simulate the mechanical behavior of a layered rock mass under true triaxial compression by Mohr-Coulomb linear model. Most calculation software are based on the Mohr-Coulomb model, which represents rock strength with cohesion and internal friction angle [4–7]. However, the joint surface is widely distributed in the rock, thus leading to the nonlinear characteristics of rock mass strength and limiting the effectiveness of the linear model in reflecting characteristics [8]. Therefore, scholars introduced nonlinear models, such as the Hoek-Brown model [9, 10] and JRC/JCS model [11–13], which can describe rock mass strength. Halakatevakis and Sofianos [14] investigated the Hoek-Brown criterion analytically through an extended plane of weakness theory. Babanouri et al. [13] developed extended Barton’s shear failure criterion for rock joints to consider the effect of various paths of normal loading/unloading before shearing and overconsolidation ratio in a fracture. Park and Song [15] produced a rough joint with a joint roughness coefficient (JRC) value ranging from 10 to 20 in an intact sample by defining the jointcontacts along a predefined joint surface. Studies have contributed significantly to the nonlinear description of rock mass strength. However, the JRC/JCS model, which was proposed by Barton et al. on the basis of a large number of shear tests on jointed rock mass, is generally believed to describe the characteristics of joints better than other models [11, 16]. The JRC/JCS model is mainly studied by theoretical and experimental approaches and not by numerical approaches. On the basis of the above considerations, this paper established the relationship between the parameters of the JRC/JCS and Mohr-Coulomb models. Thereafter, the numerical implement method and implementation process of the JRC/JCS model were discussed and the reliability of the numerical method was verified by the shear tests of jointed rock mass. Finally, the effect of the JRC/JCS model parameters on the shear strength of the joint was analyzed.

#### 2. Theoretical Derivation

##### 2.1. Relationship between the Parameters of the JRC/JCS Model and Mohr-Coulomb Model

The JRC/JCS model was formulated as follows [11]: where is the rock shear strength; is the joint normal stress; is the basic friction angle of the rock, which can be fixed at 30°; JRC is the joint roughness coefficient, which has a value that is related to the joint rough shape; JCS is the rock compressive strength and has a slight effect on shear strength under low stress conditions. The influence of JCS increases with increasing normal stress.

The Mohr-Coulomb model was formulated as follows: Equation (2) shows that the strength of the joint surface depends on the following strength parameters: cohesion and internal friction angle in the Mohr-Coulomb model [17]. The derivation of (2) is as follows: By integrating (1) and (3), we obtained the following: where , .

The triangle transformation of (4) is as follows: The formula of became unacceptable when , . Therefore, Barton suggested that should be smaller than 70° in the actual project. The minimum of the normal stress could be obtained by :

##### 2.2. Relationship Analysis of the Model Parameters

Assuming that the rock basic friction angle was 30° and that MPa to 105 MPa and to 18 were changed, Figures 1 and 2 could be determined. The cohesion and internal friction angle increase gradually with increasing JRC (Figure 1), and the slope of the relation curve of and JRC increases gradually with increasing JRC. By contrast, the slope of the relation curve of and JRC decreases gradually with increasing JRC. The curve of the effect of JRC on the internal friction angle (Figure 1(b)) demonstrates that the internal friction angle that corresponds to the same JRC increases with increasing rock compressive strength. Furthermore, for the same increment of JCS, the increment of the internal friction angle is greater with increasing JRC. However, for the same increment of JRC, the range ability of the internal friction angle decreases with increasing JCS. This observation highlights that the effect of JRC on the internal friction angle is greater than the effect of JCS. Cohesion increases linearly with increasing JCS, whereas the internal friction angle increases nonlinearly with increasing JCS (Figure 2). The slope of the curve of JCS and cohesion increases gradually with increasing JRC; thus, a rougher joint surface corresponds to the greater influence of JCS on cohesion (Figure 2(a)). When , that is, the roughness of the joint surface is high, the internal friction angle of joint surface is the basic friction angle with a value of 30° (Figure 2(b)).

#### 3. Numerical Implementation

##### 3.1. Numerical Computation Method

To establish the numerical computation method in the JRC/JCS model, the incremental formula of stress and deformation was obtained by elastic incremental theory: where , , is the initial stress, and is the stress increment. is the elastic incremental strain in the direction.

The stress components could be calculated by elastic theory: where and are the normal stress and shear strength updated by calculation, respectively.

Thereafter, new stress components , , and could be deduced by the theory that incremental strain is the sum of the elastic incremental strain and plastic incremental strain under the plastic state: By integrating (7) into (11), we obtained the following: because of where and are the normal stress and shear strength updated by calculation, respectively. is the angle of and normal direction of the joint surface.

##### 3.2. Verification of the Numerical Method

The VC++ language was used for the secondary development of the numerical calculation module of FLAC3D, and the corresponding program was developed [18–20]. During the calculation, the strain and the stress of every unit were calculated by elastic theory. Thereafter, to judge if the yield condition was achieved, the corresponding stress should be adjusted to meet the yield function by adding the JRC/JCS model. Test data that comprised the shear strength and normal stress of the joint obtained by joint shear tests in the laboratory were used to verify the correctness of the program, as shown in Table 1. The is taken for calculation, and the numerical model is shown in Figure 3. As described by Lin et al. [21], FLAC3D is difficult to use for building large, complex, and three-dimensional mining models. Lin et al. [21] combined the advantages of FLAC3D in numerical calculation and those of SURPAC in three-dimensional modeling and compiled the interface program. In this paper, we used ANSYS to build the model then transformed it to FLAC3D. The comparison between the calculated and test results is shown in Table 1 and Figure 4. The joint shear strength obtained from the test increases non-linearly with increasing normal stress. Moreover, the results from the Mohr-Coulomb model and JRC/JCS model are in good agreement with the test results. However, when the normal stress of the joint surface is small, the nonlinear characteristics of the shear strength are obvious and the Mohr-Coulomb model cannot describe the nonlinear characteristics. By contrast, the JRC/JCS model can describe the characteristics well. The correlation coefficient of the results from the JRC/JCS model and test model is 0.99315, which is higher than the correlation coefficient of the results of the Mohr-Coulomb model at 0.99072. Therefore, the reliability of the numerical calculation method of the JRC/JCS model is verified. The calculations show that the parameters of the JRC/JCS model are and MPa, whereas the parameters of the Mohr-Coulomb model are MPa and .

##### 3.3. Analysis of Parameter Effects

To study further how JRC and JCS affect the shear strength of the joint surface, MPa was fixed through calculation and the relationship of JRC and JCS with the shear strength of the joint surface was obtained (Figures 5 and 6). In Figures 5 and 6 the shear strength of the joint surface increases non-linearly with increasing JRC and JCS. The slope of the relation curve of the shear strength and JRC increases with increasing JRC, whereas the slope of the relation curve of shear strength and JCS decreases with increasing JCS. Thus, changes in JRC affect the shear strength more than changes in JCS.

#### 4. Conclusions

(1) The cohesion and internal friction angle increase gradually with increasing JRC. The internal friction angle that corresponds to the same JRC increases with increasing rock compressive strength. Furthermore, for the same increment of JCS, the increment of the internal friction angle is greater with increasing JRC. The effect of JRC on the internal friction angle is greater than the effect of JCS. Cohesion increases linearly with increasing JCS, whereas the internal friction angle increases non-linearly with increasing JCS.

(2) The results from the Mohr-Coulomb model and JRC/JCS model are in good agreement with the test results. However, when the normal stress of the joint surface is small, the nonlinear characteristics of the shear strength are obvious and the Mohr-Coulomb model cannot describe the nonlinear characteristics. By contrast, the JRC/JCS model can describe the characteristics well.

(3) The shear strength of the joint surface increases non-linearly with increasing JRC and JCS. Changes in JRC affect the shear strength more than changes in JCS.

#### Conflict of Interests

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

#### Acknowledgment

The research project was supported by the National Natural Science Foundation of China (41202057).

#### References

- H. Lin, P. Cao, and Y. Wang, “Numerical simulation of a layered rock under triaxial compression,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 60, pp. 12–18, 2013. View at Google Scholar - L. C. Li, C. A. Tang, W. C. Zhu, and Z. Z. Liang, “Numerical analysis of slope stability based on the gravity increase method,”
*Computers and Geotechnics*, vol. 36, no. 7, pp. 1246–1258, 2009. View at Publisher · View at Google Scholar · View at Scopus - T. Unlu and H. Gercek, “Effect of Poisson's ratio on the normalized radial displacements occurring around the face of a circular tunnel,”
*Tunnelling and Underground Space Technology*, vol. 18, no. 5, pp. 547–553, 2003. View at Publisher · View at Google Scholar · View at Scopus - H. Lin, W. Xiong, and P. Cao, “Stability of soil nailed slope using strength reduction method,”
*European Journal of Environmental and Civil Engineering*, vol. 17, pp. 872–885, 2013. View at Google Scholar - T. Liu, P. Cao, and H. Lin, “Evolution procedure of multiple rock cracks under seepage pressure,”
*Mathematical Problems in Engineering*, vol. 19, pp. 1–12, 2013. View at Google Scholar - S. Y. Wang, S. W. Sloan, C. A. Tang, and W. C. Zhu, “Numerical simulation of the failure mechanism of circular tunnels in transversely isotropic rock masses,”
*Tunnelling and Underground Space Technology*, vol. 32, pp. 231–244, 2012. View at Google Scholar - D. Sun, Y.-P. Yao, and H. Matsuoka, “Modification of critical state models by Mohr-Coulomb criterion,”
*Mechanics Research Communications*, vol. 33, no. 2, pp. 217–232, 2006. View at Publisher · View at Google Scholar · View at Scopus - A. Khani, A. Baghbanan, S. Norouzi, and H. Hashemolhosseini, “Effects of fracture geometry and stress on the strength of a fractured rock mass,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 60, pp. 345–352, 2013. View at Google Scholar - E. Hoek and C. Carranza-Torres, “Corkum B. Hoek-Brown failure criterion-2002 edition,” in
*Proceedings of NARMS-TAC*, pp. 267–273, 2002. - E. Hoek and E. T. Brown, “Practical estimates of rock mass strength,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 34, no. 8, pp. 1165–1186, 1997. View at Google Scholar · View at Scopus - N. Barton and S. Bandis, “Review of predictive capabilities of JRC-JCS model in engineering practice,”
*Publikasjon*, vol. 182, pp. 1–8, 1991. View at Google Scholar · View at Scopus - S. Du, Y. Hu, X. Hu, and X. Guo, “Comparison between empirical estimation by JRC-JCS model and direct shear test for joint shear strength,”
*Journal of Earth Science*, vol. 22, no. 3, pp. 411–420, 2011. View at Publisher · View at Google Scholar · View at Scopus - N. Babanouri, S. Karimi Nasab, A. Baghbanan, and H. R. Mohamadi, “Over-consolidation effect on shear behavior of rock joints,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 48, no. 8, pp. 1283–1291, 2011. View at Publisher · View at Google Scholar · View at Scopus - N. Halakatevakis and A. I. Sofianos, “Correlation of the Hoek-Brown failure criterion for a sparsely jointed rock mass with an extended plane of weakness theory,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 47, no. 7, pp. 1166–1179, 2010. View at Publisher · View at Google Scholar · View at Scopus - J.-W. Park and J.-J. Song, “Numerical simulation of a direct shear test on a rock joint using a bonded-particle model,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 46, no. 8, pp. 1315–1328, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Zhao, “Joint surface matching and shear strength. Part B: JRC-JMC shear strength criterion,”
*International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts*, vol. 34, no. 2, pp. 179–185, 1997. View at Publisher · View at Google Scholar · View at Scopus - H. Lin, W. Zhong, P. Cao, and T. Liu, “Variational safety factors and slip surfaces of slope using three-dimensional strength reduction analysis,”
*Journal of the Geological Society of India*, vol. 82, pp. 545–552, 2013. View at Google Scholar - M. Cai, “Influence of stress path on tunnel excavation response—numerical tool selection and modeling strategy,”
*Tunnelling and Underground Space Technology*, vol. 23, no. 6, pp. 618–628, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. Rutqvist, “Status of the TOUGH-FLAC simulator and recent applications related to coupled fluid flow and crustal deformations,”
*Computers and Geosciences*, vol. 37, no. 6, pp. 739–750, 2011. View at Publisher · View at Google Scholar · View at Scopus - Q. Zhang, B. S. Jiang, X. S. Wu, H. Q. Zhang, and L. J. Han, “Elasto-plastic coupling analysis of circular openings in elasto-brittle-plastic rock mass,”
*Theoretical and Applied Fracture Mechanics*, vol. 60, pp. 60–67, 2012. View at Google Scholar - H. Lin, T. Liu, J. Li, and P. Cao, “A simple generation technique of complex geotechnical computational model,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 863104, 8 pages, 2013. View at Publisher · View at Google Scholar