Primary Investigation on Equivalent Anchoring Method of Jointed Rock Mass Tunnel and Its Engineering Application

Rock bolts, one of the main support structures of the tunnel, can improve the stress state and mechanical properties of the surrounding rocks. +e rock bolts are simulated by bar or beam elements in present numerical calculations for most 2D tunnel models. However, the methods of simulating rock bolt in three-dimensional models are rarely studied. Moreover, there are too many rock bolts in the long-span tunnel, which are hardly applied in the 3D numerical model.+erefore, an equivalent anchoring method for bolted rock masses needs to be further investigated. First, the jointed material model is modified to simulate the anisotropic properties of surrounding rock masses. +en, based on the theoretical analysis of rock bolts in reinforcing mechanical properties of the surrounding rock masses, the equivalent anchoring method of the jointed rock mass tunnel is numerically studied. +e equivalent anchoring method is applied to the stability analysis of a diversion tunnel in Western China. From the calculation results, it could be found that the reinforcement effect of rock bolts could be equivalently simulated by increasing the mechanical parameter value of surrounding rocks. For the jointed rock mass tunnel, the cohesion and internal friction angle of the surrounding rocks are improved as 1.7 times and 1.2 times of the initial value, which can simulate the reinforcement effect of rock bolts. Comparing with analytical results, the improved internal friction angle is nearly consistent with analytical result. +e reinforcement effect of rock bolts is simulated obviously when the mechanical parameters of surrounding rocks are increased simultaneously. +e engineering application shows that the equivalent anchoring method can reasonably simulate the effect of rock bolts, which can provide reference for stability analysis of three-dimensional tunnel simulations.


Introduction
Rock bolts are widely used in tunneling, mining, and other underground engineering projects for more than 100 years, which can improve the strength and stiffness of surrounding rock masses. e finite element method has been extensively employed to study the mechanism of bolting, and the rock bolts are mostly simulated by bar or beam elements for twodimensional models. However, the three-dimensional rock bolt element of FEM is rarely studied and used in 3D tunnel simulations. Moreover, for the long-span tunnel, the rock bolts are numerous and hardly installed in the 3D tunneling models.
erefore, an equivalent simulation method for bolted rock masses that will simplify the difficulty of the numerical calculation needs to be further investigated.
Rock bolt is a part of a series of support methods, denoted as "rock reinforcement system." e mechanism of the bolted rock masses is widely investigated by analytical analysis, laboratory and field tests, and numerical simulations, which includes the bonding capacity of rock bolts and the interfaces properties between the bolts and rocks [1][2][3][4][5][6][7]. Some investigations focus on the constitutive relationships of the bolted rock masses. Chen proposed an equivalent constitutive model for jointed rock masses reinforced by fully grouted bolts, and the model was implemented in a finite element program [8]. Densie regarded the bolted rock masses as a homogenized anisotropic medium and proposed an elastoplastic constitutive model for bolted rock masses [9]. Carranza presented a closed-form solution for stress and displacement distributions around a circular tunnel excavated in elastic material and reinforced by rock bolts [10]. Maghous proposed a three-dimensionally theoretical and numerical model for the behavior of tunnels reinforced by rock bolts [11]. Many similar studies have been published on the constitutive model of the bolted rock masses from different viewpoints [12][13][14][15][16]. ough the constitutive models can accurately reflect the mechanical properties of the bolted rock masses, the implementation and subsequent computational process are quite difficult and time-consuming.
To simplify the numerical calculations of the bolted rock masses, many investigations assume that the constitutive model of the rocks is not changed after the bolts installed, but the mechanical parameters of the bolted rock masses are strengthened comparing with the rock masses before rock bolts are applied [17][18][19]. Most published literature indicated that the rock reinforcement effect could be equivalently simulated by improving the parameters of surrounding rock masses [3,15,20,21] but the increased value of the parameters such as elastic modulus, cohesive, and internal friction angles of the surrounding rock masses needs to be further investigated.
e existing experimental results have indicated that the layered rocks display a strong anisotropy of strength and deformation, which is usually controlled by the orientation of bedding planes [22][23][24][25]. e anisotropic properties of the layered rocks were investigated by experiments, analytical solutions, and numerical simulations in the past [26][27][28][29][30][31][32][33][34]. However, the surrounding rock masses are assumed to be isotropic in most numerical modeling. erefore, a reasonable constitutive relationship of the layered rock masses should also be needed before investigating the equivalent anchoring method.
In our study, the jointed material model is modified by introducing the anisotropic elastic matrix into the model, and the modified model is numerically implemented and verified. Based on the modified jointed material model, the reinforcement effect of rock bolts installed during the excavation process of the jointed rock mass tunnel is investigated. e equivalent simulation method of improving the parameters of the surrounding rocks for the bolted rock masses is studied through a series of numerical simulations. en, the equivalent anchoring method is adopted for the stability analysis of a division tunnel in Western China. Figure 1 shows the five elastic constants, E 1 and ] 1 are the elastic modulus and Poisson's ratio parallel to transversely isotropic surfaces, and E 2 , ] 2 , and G 2 are the elastic modulus, Poisson's ratio, and shear modulus perpendicular to transversely isotropic surfaces, respectively. e five elastic constants can be determined by laboratory tests, such as uniaxial compression, triaxial compression, and Brazilian tests. Amadei and Pan [30] suggested that the stress matrix (σ x , σ y , σ z , τ yz , τ zx , τ xy ) and strain matrix (ε x , ε y , ε z , c yz , c zx , c xy ) of the layered rock mass were related not only to the five elastic parameters (E 1 , ] 1 , E 2 , ] 2 , G 2 ) but also to the dip angle of the joint set (the angle between the bedding planes and horizontal line: θ):

Anisotropic Elastic Parameters.
e relationship between axial stress (σ y ) and strain (ε x , ε y , ε z ) can be expressed as where a 12 , a 22 , and a 23 are the functions of the jointed surface direction θ and the five elastic parameters (E 1 , ] 1 , E 2 , ] 2 , G 2 ), through the following relationships: en, the elastic parameters of each joint dip angle are given by Both the elastic modulus and Poisson's ratio vary as functions of the joint dip angles: e anisotropic elastic parameters of the layered sandstone were measured through the uniaxial tests by Zhang [35]; then, the elastic parameters of the layered sandstone with different inclination angles can be predicted by using equation (4) theoretically. Table 1 shows the comparison of the experimental and theoretical results of elastic parameters of the layered sandstone. It can be found that the difference of the elastic modulus of the layered sandstone with different inclination angles between theoretical and experimental value is negligible, and the biggest deviation is 1.68%. However, the difference of Poisson's ratio is a little bigger; the biggest deviation reaches 25.13%. As a whole, the elastic parameters of the layered rock mass with different inclination angles predicted by using equation (4) are close to the experimental results.

Modified Joint Material Model.
e jointed material model employed allows for the opening and frictional sliding of the joints. Additionally, the model includes a bulk material failure mechanism based on the Drucker-Prager failure criterion. e elastic matrix in jointed material model is isotropic, just having two independent elastic parameters (elastic modulus and Poisson's ratio). To consider the anisotropic nature of the jointed rock mass, the jointed material model has been modified by introducing the anisotropic parameters of the layered rock mass (equations (1)- (7)). e elastic matrix should be modified by equation (1), which can simulate the elastic properties of the layered rock mass with different inclination angles. e FORTRAN language is utilized for the model implementation. Figure 2 presents the geometry and basic condition of the strip foundation on layered rock mass discretized into a 2D model (10 m × 10 m). e width of the strip foundation is 2 m, and the pressure on the layered rock is 50 kPa. e basic properties of the rock matrix are elastic modulus of E 1 � 20 GPa and E 2 � 30 GPa, Poisson's ratio of υ 1 � 0.20 and υ 2 � 0.30, shear modulus of G 2 � 9.68 GPa, and density of 2.3 × 10 3 kg/m 3 . Soft interfaces have a normal stiffness of 7.5 GPa/m, and tangential stiffness of 5 GPa/m. e joint connectivity rate is 1, and joint spacing is 0.1 m. Figures 3 and 4 present the results of von Mises stress and displacement of the strip foundation which is related to the changes of dip angle. e distribution of stress and deformation is symmetric for horizontal layered rock mass and vertical layered rock mass and asymmetrical for other dip angles of soft interfaces. e curves of stress and displacement of the selected measuring point and their relation with the dip angle are shown in Figure 5. At the dividing dip angle of 90°, the curves are symmetric within the range [0°, 180°], and node displacement tends to be minimum in the vicinity of 90°. us, the modified model can effectively depict the anisotropic characteristics of the layered rock masses.

Primary Investigations on the Equivalent
Anchoring Method It assumes that the surrounding rocks are well integrated, incomplete development of bedding planes, isotropic, and axially symmetric. en, the surrounding rock masses of circular tunnel will be spiral failure after excavation. Two sets of slip lines (β � ±((π/4) + (φ/2))) will be generated in the plastic zone. e body of rock bolt will act as a "pin effect" to increase the shear strength of failure surfaces because of the generated slip lines. From the experimental results by Dulacska [36], it can be obtained that where F is the shear stress of rock bolt and F max is the maximum shear stress of rock bolt. According to the von Mises failure criterion, F max can be expressed as ) for the condition of pure failure. d is the diameter of rock bolt and σ s is the tensile strength of rock bolt. e shear strength of the failure surfaces is enhanced due to the effect of rock bolt; therefore, the equivalent cohesion of the bolted rock masses can be expressed as where c 0 and φ 0 are the initial cohesion and internal friction angle of surrounding rock masses and s l and s c are the axial and circumferential spacing of rock bolts along the tunnel, which can be seen in Figure 6.

Equivalent Internal Friction Angle of Surrounding Rock
Masses. e mechanical parameters of rock masses and rock bolts both have an impact on the equivalent strength of the bolt-rock composite structure. According to the bolt density factor proposed by Indraratna [17], there is where c is the circumferential angle between the rock bolts, D is the diameter of tunnel, and η is the friction coefficient between rock bolt and rock mass, which relates to the surface roughness of rock bolt. When the threaded rock bolts are adopted, η � tan φ 0 ; otherwise, η � tan(φ 0 /2).

Advances in Civil Engineering
Based on the theory of fully grouted anchor, the same deformation will be generated by rock bolts and surrounding rocks after rock bolts installed. At this time, the shear stress generated at the surface of rock bolts prevents the surface deformation of surrounding rocks. is kind of reinforcement performs as the compressive strength (σ c * ) of surrounding rocks increased. erefore, the equivalent compressive strength of the bolt-rock composite structure can be expressed as Assume that the bolt-rock composite structure satisfies the Mohr-Coulomb yield criterion. e gradient of the yield locus in principle stress coordinate system before and after the bolts supporting can be defined as f 0 and f 0 ′ , expressed as e intercept of gradient f 0 ′ increases after rock bolts are applied; combining with density factor of rock bolts, this variation relationship can be expressed as Equation (11) has been verified by Intraratna's experiments [17], the equivalent internal friction angle of the boltrock composite structure can be obtained through solving equation (11):

Equivalent Elastic Modulus of Surrounding Rock
Masses. Since the elastic modulus (E b ) of the rock bolt is much larger than the elastic modulus (E 0 ) of surrounding rocks, therefore, the codeformation between rock bolts and surrounding rocks will make a difference. For hard rock condition, the rock bolts have little impact on deformation parameters of surrounding rocks, and the deformation parameters have little impact on the evaluation of surrounding rocks stability. It can be assumed that Poisson's ratio is unchanged after rock bolts are applied, and the bolted rock masses are isotropic. According to the experimental results of Carranza [10], the equivalent elastic modulus of surrounding rock masses can be expressed as

Numerical Model and Calculation
Parameters. e excavation process of a jointed rock masses tunnel is numerically simulated by different configurations including supporting with rock bolts and without rock bolts and improving the equivalent mechanical parameters of surrounding rocks (the zone of three times to five times of the tunnel diameter). e plastic zone of surrounding rocks and the displacement of key points (tunnel vault and tunnel floor) are compared with different configurations. e cohesion and internal friction angle of surrounding rocks are selected as the parameters for equivalent supporting by improving its initial value.
e size of numerical model is 60 m × 60 m, and the diameter of tunnel is 10 m. e inclination angle of the jointed rocks is θ � 30°. It assumes that the numerical model

Advances in Civil Engineering
is under isotropic stress field of 4 MPa along the X and Y directions of the model's boundary, and the bottom constraints of the model are fixed. Figure 7 is the diagram of the numerical model. In order to simplify the numerical model, the bar element and beam element are adopted to simulate the rock bolt and lining, respectively. e thickness of lining is 0.4 m. e diameter of rock bolts is 25 mm. e range interval distance of rock bolts is 1.3 m, and the rock bolts are applied perpendicular to the tunnel cycle.
e modified jointed material model is adopted for the surrounding rocks. Table 2 lists the mechanical parameters of the surrounding rocks, rock bolts, and lining.

Numerical Results and Analysis.
e plastic zone can reflect stability and failure of the surrounding rock masses. Figure 8 shows the plastic zone of surrounding rocks after excavation of different configurations. Comparing Figures 8(a) and 8(b), it can be found that the size and range of plastic zone supporting with rock bolts are obviously smaller than those without rock bolts installed, and the range of the plastic zone is directional which related to the bedding plane angle. It demonstrates that the rock bolts play an important role in controlling deformation and failure area of the surrounding rock masses. Figures 8(c)∼8(e) are the plastic zone of surrounding rock mass under equivalent supporting by taking mechanical parameters of surrounding rocks as 1.7c, 1.3φ, and 1.3(c, φ), respectively. Comparing with Figure 8(b), it can be found that the plastic zone of surrounding rocks by equivalent supporting is basically analogous with the plastic zone of Figure 8(b). erefore, the equivalent supporting method can simulate the reinforcement effect of rock bolts for jointed rock mass tunnel. e equivalent mechanical parameter of cohesion is taken as 1.7c to accurately simulate the reinforcement effect of rock bolts for jointed rock mass tunnel. However, the internal friction angle of surrounding rocks is only taken as 1.3φ. Obviously, it can be found that degree of sensitivity to simulate the reinforcement effect of rock bolts through improving the equivalent mechanical parameters of surrounding rocks is different, and the most sensitivity parameter is the internal friction angle and then cohesion.
e plastic zone of surrounding rocks cannot roundly reflect the reinforcement effect of rock bolts simulated by equivalent supporting. e vertical displacement of dome and horizontal displacement of arch foot are selected as the comparison objects. e numerical results can be seen in Tables 3-5. e numerical results of supporting with and without rock bolts installed are also given in all the tables.
In Tables 3-5, it can be seen that improving the parameters of surrounding rocks can simulate the reinforcement effect of rock bolts for jointed rock mass tunnel, but the improved value of each parameter is different. e value of cohesion and internal friction angle improved independently are 1.7 times and 1.2 times of the initial value, which can equivalently simulate the reinforcement effect of rock bolts. It also can be found that the degree of sensitivity to simulate the reinforcement of rock bolts by improving the mechanical parameters of surrounding rocks is different. e most sensitivity parameter is internal friction angle and then cohesion. is finding is consistent with the analysis of plastic zone distribution of surrounding rocks. When the equivalent parameter of cohesion is taken by 1.7c, the vertical displacement of dome is close to the result supporting by rock bolts, and the difference of the horizontal displacement of arch foot between equivalent supporting and supporting with rock bolts is 10.8%. With the increasing of cohesion, the displacement of lining is decreased. When the equivalent parameter of internal friction angle is taken by 1.2φ, both the vertical displacement of dome and the horizontal displacement of arch foot are close to the result calculated by supporting with rock bolts. erefore, it can simulate the reinforcement effect of rock bolts vividly. Table 5 shows that when the shear strength parameters of surrounding rocks are increased by 1.1 times the initial value, the result of vertical displacement of lining vault is very closed to the result calculated by supporting with rock bolts. Similarly, with the increasing shear strength parameters, the displacement of lining is decreased firstly and then increased because the internal friction angle is the most sensitivity parameter of surrounding rocks to simulate the reinforcement effect of rock bolts equivalently.

Comparison between Numerical and Analytical Results.
To equivalently simulate the reinforcement effect of rock bolts, the cohesion and internal friction angle of the surrounding rocks need to improve as 1.7 times and 1.2 times of the initial value according to the numerical results. Based on the parameters of the numerical model, the equivalent cohesion and internal friction angle can be calculated by equations (7) and (12). e equivalent cohesion and internal friction angle of surround rocks are c 0 ′ � 0.699 MPa and φ 0 ′ � 44.231°, respectively. is means that the cohesion and internal friction angle of the surrounding rocks are improved by 7.5% and 26.3%. From the comparison between numerical and analytical results, it can be found that the deference of internal friction angle is small, and the deference of cohesion is a little bigger. It should point out that the analytical results calculated are based on the assumption that the surrounding rocks are isotropic. However, the surrounding rocks in the numerical simulation are assumed as anisotropic material. In the other words, the internal friction angle of surrounding rocks improved after rock bolts installed may be analogous for isotropic and anisotropic rock mass tunnel.

Geological Conditions.
A water diversion tunnel is in the west of China, in which terrain in a whole is slightly fluctuating, and the north of the tunnel area is lower than the south, and the west of the tunnel area is lower than the east. e terrain of tunnel is slowly inclined from the northwest to the southeast, and the bedrock is mostly bare. e section of the tunnel is 46.153 km long, which is constructed by TBM and drilling and blasting. According to the analysis of the borehole data along the tunnel, the bottom lithology revealed by the tunnel is mainly of Devonian tuffaceous sandstone, tuff, and carbonaceous siltstone; Carboniferous andesite, calcareous glutenite, calcareous siltstone; and Hualixi period Quartz diorite and granite. e stability of a cross caverns is analyzed during excavation process in the section of drilling and blasting in this paper. e rocks in this section are mainly of Devonian tuff; the inclination angle of the rocks is 30°. e rock masses are layered, fresh, and intact. e stability of the surrounding rocks is very well.

Numerical Model and Calculation Parameters.
e drilling and blasting section of the tunnel has a main carven and a lateral carven, which cannot be simplified as a plane strain model. It needs to establish a three-dimensional model to analyze the stability during the excavation process.   1(c, φ)), and influence of the rock type is ignored. According to the geological data, mechanical parameters of the rocks and lining are listed in Table 6. e process of tunnel excavation can be simulated by two configurations: (1) only supported by lining and (2) Figures 10 and 11 show the maximum and minimum principle stress of the surrounding rocks after the main cavern is excavated. e maximum principle stress appears at the floor and the middle sidewall of the tunnel. For the equivalently bolted rock masses, the maximum principle stress is compressive strength; however, the tensile stress appears in the surrounding rocks without rock bolts are applied. Stress concentration may appear at the junction of the sidewall and bottom floor, and the maximum compressive stress is about 16 MPa. e compressive stress at the roof of the tunnel is a little bigger, and the stress is about 12 MPa. Figure 12 shows the vertical displacement of the surrounding rocks after the excavation of main cavern. For the equivalently bolted rock masses, the vertical displacements of the roof and floor are 8 mm and 3.45 mm, respectively. However, for the surrounding rocks without rock bolts being applied, the vertical displacement of the roof and floor is bigger. It is demonstrated that the equivalent anchoring method can control the deformation of the surrounding rock masses. Figure 13 presents the plastic zone of the surrounding rocks, for the equivalently bolted rock masses, the plastic zone only appears at the middle of the sidewall, and the range of the plastic zone is small. However, for the surrounding rocks without rock bolts being applied, the plastic zone appears at the middle of all sidewalls, and the floor of the tunnel. e range of the plastic zone is wide. Figures 14 and 15 show the maximum and minimum principle stress after the lateral cavern is excavated. For the zone of the lateral cavern, the maximum principle stress appears at the roof of the cross section between the main cavern and lateral cavern, which the tensile stress may emerge in this area. It should note that the tensile failure will be occurred. After the lateral cavern is excavated, the compressive stress at the sidewall of the lateral cavern is about 9 MPa, which the compressive failure may be occurred. Figure 16 shows the plastic zone of the surrounding rocks after the lateral cavern is excavated. For the equivalently bolted rock masses, the plastic zone only appeared at sidewall of the cross section between the main cavern and lateral cavern. It is closer to the actual construction of the tunnel. It demonstrates that the equivalent anchoring method can reflect the failure zone of the surrounding rocks during the tunnel excavation. Meanwhile, 17 shows that the vertical displacement of the lateral cavern roof is 1.39 cm for    the surrounding rocks without rock bolts being applied. However, for the equivalently bolted rock masses, the vertical displacement is 1.20 cm. Figures 18 and 19 are the principle stress of the lining. Figure 18 shows that the most area of the lining is under the condition of tension, and the maximum tensile stress is less than 1.78 MPa. Figure 19 shows that the compressive stress of the color area of the lining is bigger than 16.7 MPa, and the maximum compressive stress may be occurring at the cross section between the roof and sidewall of lining, and the cross section between the sidewall and floor of lining.        Advances in Civil Engineering 13

Conclusions
e reinforcement effect of surrounding rock masses installed with rock bolts is equivalently simulated based on the assumption that the constitutive model of the jointed rock masses is anisotropic. e equivalent anchoring method of the jointed rock mass tunnel is proposed, and the method is applied to the stability analysis of a division tunnel in Western China. Some conclusions can be obtained from the studies in this paper.
(1) Considering the anisotropic properties of the layered rocks, the elastic matrix of the jointed material model is modified. e model verification shows that the modified jointed material model can reflect the anisotropic properties of the layered rocks. (2) Based on the modified jointed material model, the equivalent anchoring method of the jointed rock mass tunnel is numerically studied. e results of a series of numerical simulations show that the reinforcement effect of rock bolts can be equivalently simulated by improving the shear strength parameters of surrounding rocks. In order to simulate the reinforcement effect of rock bolts, the cohesion and internal friction angle of the surrounding rocks need to improve as 1.7 times and 1.2 times of the initial value. e shear strength parameters of surrounding rocks, improved 1.1 times the initial value, can accurately simulate the reinforcement effect of rock bolts. (3) e excavation process of a division tunnel in Western China is numerically studied in this paper. e simulation results show that the equivalent anchoring method can reflect the reinforcement effect of the rock bolts. e numerical results are close to the actual conditions of the tunnel construction. It demonstrates that the equivalent anchoring method can be used to the simulation of 3D tunnel construction process. It solves the problem of the simulation method of rock bolts in three dimensions and many rock bolts applied in three-dimensional model.

Data Availability
e data used to support the findings of this study are included within the article.    14 Advances in Civil Engineering