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)).
| (a) Effect of JRC on |
| (b) Effect of JRC on |
| (a) Effect of JCS on |
| (b) Effect of JCS on |
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).