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

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.


Introduction
The Mohr-Coulomb linear model is generally used to study the characteristics of rock mass strength [1][2][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][5][6][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][12][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.

Relationship between the Parameters of the JRC/JCS
Model and Mohr-Coulomb Model. The JRC/JCS model was formulated as follows [11]: The Scientific World Journal 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: The triangle transformation of (4) is as follows: The formula of became unacceptable when → 0, + JRC ⋅ log 10 (JCS/ ) → ∞. Therefore, Barton suggested that +JRC⋅log 10 JCS/ should be smaller than 70 ∘ in the actual project. The minimum of the normal stress could be obtained by + JRC ⋅ log (JCS/ ) = 70 ∘ :

Relationship Analysis of the Model Parameters.
Assuming that the rock basic friction angle was 30 ∘ and that JCS = 5 MPa to 105 MPa and JRC = 0 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 JRC = 0, 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)).

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 1 = + 4 /3, 2 = − 2 /3, 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 0 and Δ are the normal stress and shear strength updated by calculation, respectively. Thereafter, new stress components 1 , 2 , and 3 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 The Scientific World Journal where and are the normal stress and shear strength updated by calculation, respectively. is the angle of 1 and normal direction of the joint surface.

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][19][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 JRC = 16 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 4 The Scientific World Journal

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.
The Scientific World Journal 5  (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 nonlinearly with increasing JRC and JCS. Changes in JRC affect the shear strength more than changes in JCS.