Advances in Civil Engineering

Volume 2010 (2010), Article ID 149810, 13 pages

http://dx.doi.org/10.1155/2010/149810

## Bolt-Grout Interactions in Elastoplastic Rock Mass Using Coupled FEM-FDM Techniques

^{1}Department of Mining Engineering, Indian Institute of Technology, Kharagpur 721302, India^{2}Department of Mathematics, Indian Institute of Technology, Kharagpur 721302, India

Received 28 January 2010; Revised 21 April 2010; Accepted 30 June 2010

Academic Editor: Ser Quek

Copyright © 2010 Debasis Deb and Kamal C. Das. 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

Numerical procedure based on finite element method (FEM) and finite difference method (FDM) for the analysis of bolt-grout interactions are introduced in this paper. The finite element procedure incorporates elasto-plastic concepts with Hoek and Brown yield criterion and has been applied for rock mass. Bolt-grout interactions are evaluated based on finite difference method and are embedded in the elasto-plastic procedures of FEM. The experimental validation of the proposed FEM-FDM procedures and numerical examples of a bolted tunnel are provided to demonstrate the efficacy of the proposed method for practical applications.

#### 1. Introduction

Rock bolts have been widely used as a primary support system to stabilize the rock masses around tunnel, underground mine galleries, slopes, and others structure made in rock masses. In general, rock bolts reinforce rock masses through restraining the deformation within rock masses [1] and reduces the yield region around the excavation boundary. During the last four decades, different types of rock bolts have been in practice, out of which fully grouted active/passive bolts were the most common types. For a fully grouted passive rock bolt installed in deformable rock masses, a neutral point exists on the bolt rod, where shear stress at the interface between the bolt and grout material vanishes. The pickup length is then defined as the length between free end from the tunnel boundary to the neutral point, and the shear stress along this bolt length drags the bolt towards the tunnel [1]. The bolt length between the neutral point and the other free end of the bolt (inside the rock mass) is designated as anchor length and the developed shear stress drags the bolt towards the rock mass or, in other words, anchors the bolt into the rock mass. Based on these concepts, shear stresses and axial loads developed in a bolt rod are analytically formulated by many researchers who consider the effect of faceplate at the free end of the bolt [1, 2]. Bolt grout interactions around a circular tunnel in Hoek and Brown medium have been formulated analytically considering a bolt density factor [3]. In their study, Indraratna and Kaiser [3] demonstrated that due to the installation of bolts, the Hoek and Brown parameters , , and uniaxial compressive strength of intact rock, , would change to new values , , and based on the bolt density factor which consists of bolt spacing, shear stress on the bolt grout interactions, tunnel radius, and bolt diameter. Considering different approaches to bolt performance, Stille [4] presented a closed form elasto-plastic analytical solutions of grouted bolts. Based on shear lag model, Cai et al. [5] derived an analytical solution of rock bolts for describing the interaction behaviours of rock bolt, grout material and rock mass. In addition, recently Torres [6] presented closed-form solution of elastic rock bolt in elasto-plastic rock mass for grouted and anchored rock bolt. In his work, stiffness and strength of grout material are assumed to be infinitely high. Hence, there is virtually no shear displacement in rock bolt interface. Brady and Lorig [7] numerically analyzed the interactions of bolt grout in Mohr Columb medium using finite difference method (FDM) technique and showed that radial displacement, peak tangential stress, and yielded region would reduce due to the installation of grouted bolts around a circular tunnel. FDM technique has also been used for solving differential equation, which governs the interactions between rock bolt and grout material [2].

In addition, numerous studies have been published on the analytical solution of stresses and displacements around a circular tunnel considering elasto-plastic rock mass with Hoek and Brown yield criterion [8–11]. Finite Element Method (FEM) and/or FDM-based procedures are also developed for the analysis of the said problem and have been presented in many references for solving geotechnical problems [12–15].

In this study, FDM technique is applied to evaluate the interactions between bolt and grout considering neutral point, pickup and anchor length concepts. Displacements and stress distributions in the rock mass are computed using FEM technique, considering material nonlinearity based on Hoek and Brown yield criterion. The analysis procedure of bolt grout interactions are evaluated based on FDM and embedded into the elasto-plastic procedure of FEM. The proposed FEM-FDM procedure is validated with the experimental results of pull-out tests conducted by Farmer [16]. It is found that the observed shear strains along the bolt rod are close to those obtained from numerical models. The coupled FEM-FDM technique is then applied to evaluate rock mass behaviour and bolt performance for a circular bolted tunnel. Results are obtained by varying bolt densities and shear stiffness at the bolt grout interface considering the faceplate option.

#### 2. Implementation of Hoek-Brown Yield Criterion in FEM

The generalized Hoek and Brown rock mass failure criterion for drained condition is expressed as follows [17], where are principal stresses and should be arranged as (considering compressive stress as negative), = uniaxial compressive strength of intact rock and negative value is considered, is reduced material constant based on intact property of rock, and Geological Strength Index (GSI), and and are rock mass properties based on GSI [18].

For numerical computations, it is convenient to express yield surface using mean hydrostatic stress, * σ_{m}*, and deviatoric stress invariants, namely, second invariant of deviatoric stress, and Lode angle [13, 18]. This is particularly helpful since computations of principal stresses and their directions in three dimensions can be avoided. In this study, (1) has been reformulated using invariant of stresses as and is given as
where, the expression for parameters , and = = = = = .

The angle is called Lode angle ranging between ±30° and is the Kronecker delta, whose value is unity when and zero when From (2), it is clear that the generalized Hoek and Brown yield function is dependent on all three invariants of stresses. This yield function is plotted in three mutually perpendicular principal stress axes () as shown in Figure 1. In principal stress space, the yield locus is conical whose normal section at any point is an irregular hexagon, just as Mohr-Coulomb yield surface [11]. The conical shape of this yield surface is the evidence that hydrostatic stress does influence the yielding. In addition, the apex of the hexagonal pyramid is zero (0) when mean hydrostatic stress is .

##### 2.1. Properties of the Yield Surface

Figure 2 shows the two dimensional deviatoric or plane () represent of the generalized Hoek and Brown yield surface. The yield surface is convex and has six sharp corners or points of singularity where tangent of yield locus is indeterminate. In all other points, the yield surface is regular. The singular points occur when approaches to ±30 degrees. As a result, the direction of plastic straining is also indeterminate and thus elasto-plastic constitutive matrix cannot be determined at those points on the yield surface.

##### 2.2. Flow Rule

A nonassociative flow rule assumes that the normality condition, the direction of plastic strain increment is normal to the yield surface, is not satisfied. Hence an alternative plastic potential function, must be assumed, such that . On the contrary, associative flow rule assumes that The linear plastic potential function in plane as given in (3a) has been considered in conjunction with Hoek and Brown yield criterion for the analysis of ground reaction curve and elasto-plastic analysis of circular tunnel [8, 9]. Equation (3a) is modified to obtain the expression for plastic potential function using mean stress and deviatoric stress invariants, as given in (3b). In deviatoric plane, this plastic potential surface is convex and resembles an irregular hexagon, except when zero dilation angle, it is regular, having six points of singularity when approaches ±30 degrees,where and = dilation angle.

##### 2.3. Elastic-Plastic Analysis in Plane Strain Condition

The elasto-plastic constitutive matrix for elastic-perfectly plastic (strain hardening parameter taken to be zero) is represented by the following equation [12]:
where for plane strain condition, stresses and strains are expressed in Cartesian coordinate system as , . The matrix **D** is the elastic constitutive matrix estimated based on elastic modulus, and Poisson’s ratio, which gives
For elastic stress increment, first three rows and columns are considered. Representing and , for plane strain condition, **a _{i}** (

*i*= 1,, 3) can be obtained in Cartesian coordinate system as where and .

The constants and are the derivatives of yield and plastic potential functions, respectively, with respect to , and and given in following equations: Finite element implementation of the above yield functions is performed, based on the procedure outlined in various literatures [12, 14, 19]. In this case, Newton-Raphson solution method is applied with initial stiffness conditions. Hence, large numbers of load steps are taken for satisfying convergence criterion. The increments of plastic strain and stress are estimated using the following equations:where plastic multiplier

#### 3. Analysis of Bolt-Grout Interactions Using FDM

A bolt installed in deformable rock mass is subjected to axial loading and it provides resistance to the movement of rock mass through shear stresses developed axially in the bolt-grout interfaces. The bolt is represented by a one-dimensional bar and hence the relationship between the shear force and the axial load can be established considering a small (differential) section of the bolt as shown in Figure 3. The force equilibrium in the axial direction leads to the following expression: or where shear force due to bond per unit length, is the area of the cross-section of the bolt , = length of the bolt, axial stress in the bolt, and radius of the bolt rod.

Now, for an elastic bolt with elastic modulus of , Hence, substituting (11) into (10) gives Farmer [16] proposed and justified that shear force developed due to bond between bolt and grout for a unit length of bolt. This relationship is a linear function (presumed) of the relative movement between rock and bolt [2, 20]. The relationship is expressed as where is the shear stiffness of bolt grout interface and can be determined from stress-displacement response during laboratory pullout tests. In general, depends on the bolt and rock displacements due to debonding between bolt rod and grout. Experimental studies of pull-out tests conducted by Farmer [16] showed that load-displacement curve is non-linear, which is evident from de-bonding between bolt rod and grout. However, in this present study, de-bonding or elasto-plastic behaviour of grout material has not been considered and hence shear stiffness is kept constant in (13). For yielding of bolt steel material, Tresca criterion has been implemented. In the example problem, given in the later part of the paper, the axial stresses are found to be far less than the yield strength of bolt steel and hence bolts do not yield for any condition. As a result, theoretical background of plasticity of bolt steel is omitted.

Displacement function describes the movement of rock mass over bolt length In general, will be a monotonically decreasing function with radial distance *r, *measured from a circular tunnel boundary. The form and rate of decrease of will depend on, namely. the size and shape of the opening, the strength and structure of the rock mass, stress redistribution occurring around excavation due to the yield of rock mass, and number of bolts installed at the excavation boundary.

Now, combining (12) and (13), the distribution of displacement along bolt length is described by the second order inhomogeneous linear ordinary differential equation as Equation (14) suggests that the solution of can be obtained if is known beforehand. In this study, the solution is obtained from the nodal displacements of a FEM model with bolts installed at the excavation boundary.

For this purpose, bolt is discretized along its length for implementation of finite difference scheme as shown in Figure 4. Hyett et al. [2] provided the basic formulation of (14) using the finite difference method. In the following, the significant part of that concept is mentioned for implementation of (14) in FEM procedure. In discretized form, (14) can be written with central difference form as
where the superscript *i* refers to the th discretized point on the bolt and denotes step length as . Equation (15) can be rewritten as
Equation (16) can be rearranged in more compact form as
where .

Now, using force-free boundary condition for the free end (without faceplate bolt), the first equation transforms: For a faceplated bolt, the displacements of rock mass and bolt are same at . Thus, reaction force is developed at that end. Then the unknown reaction force is determined as It may be noted that installation of faceplate ensures that the entire length of bolt acts as anchor length and the maximum tensile stress occurs at the tunnel boundary end.

For a passive grouted bolt, the other end of bolt, for both with or without faceplated, th nodal point is free. Hence, the *n*th equation becomes
Equations (16), (19), and (20) corporate linear simultaneous equations, which become a tri-diagonal matrix with *n* unknown values of . The matrix can be solved directly by Thomas algorithm or Gauss elimination method. In addition, axial tensile load on bolts is calculated from the equation
or in discretized form of above equation becomes The shear force at bolt grout interface is obtained from (13).

#### 4. Implementation of Bolt-Grout Interactions into Elasto-Plastic FEM Procedures

It is assumed that the shear force developed at the bolt-grout interface, resists the movement of rock mass along bolt axis direction. In order to implement this concept numerically, shear force at discretized points (FDM) is transformed to nodal reaction forces (FEM) of a rock element and is added to the reaction forces generated by that element. These reaction forces are applied into the next iteration in Newton-Raphson solution scheme (see (22)). In the following section, this concept is discussed in detail. Given
where global stiffness matrix at load step and iteration **i** = , if yielding of rock does not occur at the Gauss point. The above integration is performed based on Gauss-Legendre quadrature rule.:volume of the element *e*;:number of elements;:incremental displacement vector;:applied load at load step ;:reaction force vector (rock element) at load step at iteration ;:reaction force vector (bolt element) at load step at iteration .

Let us consider a 4-noded rock element is bisected by a bolt element as shown in Figure 5. The discretized points of the bolt element are *i* and *j*, respectively, and numerals 1 to 4 denote the node numbers of the rock element. It may be noted that points *i* and *j* are automatically calculated as the intersection points between bolt and edges of solid elements. These points are taken as nodal points for FDM procedure. From (13), the shear force per unit length at points *i* and *j* can be estimated after the bolt displacements are computed using FDM procedures. Let and denote the shear forces per unit length at points *i* and *j* respectively. If the values of shape functions at the midpoint (denoted by *m*) between the points *i* and *j* are then shear force acting on the bolt segment can be transferred to nodal reaction forces as

Nodal reaction force is then transformed into global coordinate system as and . In (23), the negative sign denotes shear force acting on rock mass. The Newton-Raphson method is applied to obtain solution for displacements. In this procedure, for each load step and for every iteration within a load step, Hoek-Brown yield criterion (see (2)) is evaluated at each Gauss point of a rock element with the updated trial stress, as follows. If the criterion is satisfied, plastic correction is applied to estimate the increment of plastic strain and corrected stress level (8b). If the previous state of stress was in elastic domain then transition from elastic to plastic regime is estimated as elasto-plastic ratio (a value between 0 and 1.0) for correction of total strain increment in plastic regime. With the updated stress (24), nodal reaction force for each element is computed using For each bolt element, bolt displacements are computed at each discretized points using (17). Then, shear force per unit length is obtained at each point using (13). Shear force is transformed into nodal reaction forces using (23) and added to the above nodal reaction forces. The global internal out-of-balance force vector, as given in right-hand side of (22) is applied as input parameter to N-R procedure for the next iteration. In this way, iterations are continued within a load step until convergence criterion is satisfied. Then load step is incremented and the same process continues until the total load is applied. The algorithm for coupled FEM-FDM procedure for bolt-grout interactions in elasto-plastic rock mass is given in Algorithm 1.

#### 5. Experimental Validations

The proposed FEM-FDM procedure has been validated using the results obtained from “pull-out” tests of rock bolts conducted by Farmer [16]. In general, the “pull-out” test is intended to measure the short-term strength of rock bolt anchor [21]. The strength is measured by a pull test in which bolt-head displacement is measured as a function of the applied load. In his pioneering work, Farmer [16] conducted pull-out tests of resin grouted rock bolts anchored in elastic concrete, limestone, and chalk. In his study, grouted lengths of bolts were 350 and 500 mm in concrete. The tests were conducted after 24 hours from the installation of the bolts. Five strain gauges were installed on the rod at regular interval along the length of the bolt. During testing, the anchors were loaded against surface bearing plate of 300 , at a force rate of 5 kN/min considering conditions of rigid boundary. Farmer [16] also derived an ordinary differential equation of bolt displacement, similar to (14), ignoring the displacement of rock or concrete.

Here, a concrete block of 1.0 m 1.5 m has been discretized with 4 noded quadrilateral elements of 100 mm size and grouted rock bolt of length 500 mm is loaded axially, as shown in Figure 6. In pull-out tests, since axial load is applied only on the bolt rod, axial force of the first discretized point (*i *= 1) on bolt can be rewritten using (21a) and (21b) as
or

In (26b), the axial force (pull-out load) is a known parameter. It is incorporated into (17) for *i* =1 keeping all other conditions same. shear strain along the bolt is determined using (27) and compared with those of experimental observations recorded at gauge stations,
The initial shear stiffness of bolt grout interface is calculated from load-displacement curve recorded by Farmer [16] and it is found to be 6 GN/m per unit length. In numerical models, the elastic moduli of concrete, and bolt rod are considered to be 20 GPa and 180 GPa, respectively, as given by Farmer [16]. The bolt and hole diameters are 20 mm and 28 mm, respectively.

Farmer [16] had recorded the gauge strains for different pull-out loads ranging between 20 kN and the breaking load. Shear strains along bolt length of the pull-out test are shown in Figure 7. Farmer [16] mentioned that experimental shear strain/stress distribution for pull force of 20 kN is close to the theoretical distribution which he had proposed. However, he observed substantial differences in shear stress at 40 kN pull force and this difference increased with increased force. He concluded that de-bonding between bolt rod and grout at higher level of pull force might be the reason for the differences and it was also evident from the non-linear load-displacement curve.

Numerical study conducted on the same elastic rock medium suggests that during pull-out load, the displacement of rock medium is negligible as compared to bolt displacement Hence, ignoring (27) provides the axial strains which are similar to the measured values. In this study, shear strains obtained from numerical models are compared with those of experimental results of different pull-out loads. Shear strains are estimated at discretized points along bolt length. It is found that up to 80 kN pull-out force, shear strain distribution along the bolt rod is very close to the values obtained from the experimental results (Figure 7). However, numerical results differ from those of experimental observations if the pull-out force exceeds 80 kN. The difference in results are contributed due to de-bonding between bolt rod and grout which is not considered in FEM-FDM procedure ( being constant). From this study, it is clear that FEM-FDM procedure proposed in this paper simulates pull-out test results and hence can be applied to other applications, one of which is elaborated in the next section.

#### 6. Numerical Examples of a Circular Tunnel Having Rock Bolts

A circular tunnel made in rock mass, having 4 m in diameter, is considered for numerical studies using the above FEM-FDM procedures. The finite element and finite difference codes are developed by the authors and are used for the numerical computations of the above problem. The rock mass is subjected to yield, based on Hoek and Brown criterion, as given in (1). Bolts are installed around the tunnel boundary (tunnel and bolt-grout parameters are given in Table 2). The geometry, finite element mesh and boundary conditions are shown in Figure 8(a). Finer mesh is developed near the tunnel boundary for improving the accuracy of displacement and stress distributions. Since the problem is axisymmetric, it is sufficient to analyze only one quarter of the tunnel. The models are made up with 352 quadrilateral elements having 405 nodes and are analyzed in plane strain condition. It is also considered that faceplate may be installed at the free end of the bolt (tunnel boundary side). Bolt density and shear stiffness at bolt grout interface are also varied to evaluate the performance of bolts and response of the rock mass for a uniform far field stress of 8.0 MPa.

Figure 8(b) shows the enlarged view near the tunnel boundary depicting the 9 installed bolts. It may be noted that, bolts can be installed at any part of rock element. However, in this study, they are installed along the boundary of the rock elements for symmetry. Material properties of intact rock and Hoek-Brown parameters are listed in Table 1. It may be seen that residual properties of Hoek and Brown parameters are kept same as the peak properties. Results of bolt models are compared with those of no-support model. Again, results of no-support model are compared with analytical solutions having elastic perfectly plastic conditions in Hoek-Brown media, Hence, in this study, all numerical models are analyzed considering elastic perfectly plastic conditions. The initial stiffness scheme is adopted for solving displacements using Newton-Raphson method.

##### 6.1. Results and Discussions

In the following subsections, stress and displacement distributions around the tunnel boundary, and load and shear response of bolts are discussed in details. Influence of faceplate on the bolt performance, significance of different shear stiffness of bolt grout interface, and effect of bolt density on rock mass are mainly focused in the discussion. For comparison purpose, a model without bolts is also analyzed and the corresponding displacements and stresses in the rock mass are plotted with those bolt models. From here on, the model without support will be termed as “no-support” model. The models having and not having faceplated bolts will be called as “with faceplate” and “without faceplate” models, respectively.

###### 6.1.1. Influence of Faceplated Bolts in Rock Mass Behaviour

Figure 9 shows the radial displacement of rock mass with radial distance from tunnel boundary. It can be seen that the peak radial displacement is reduced by about 15% and 25% for without and with faceplate models, respectively, as compared to no-support model. In other words, faceplated bolts reduce the peak radial displacements by about 11% as compared to bolts without faceplate. The influence of faceplate can also be observed in the stress distribution plots shown in Figure 10. It is clear that the extent of yield region around the tunnel has decreased by at least 0.25 m from no-support model due to the installation of bolts. The peak tangential stress is lower and the overall radial (confining) stress is higher in the rock mass, if faceplate is installed in the bolts. These results suggest that bolts especially with faceplates can provide effective support to the yielding rock mass around a tunnel boundary, provided sufficient shear stiffness exists between bolt and grout interface.

###### 6.1.2. Effect of Faceplate on Bolt Performance

Figure 11 depicts the bolt load and shear load per unit length obtained from with and without faceplate models. The maximum tensile load carried by a faceplated bolt is almost twice than that of without faceplated bolt. The maximum shear load in the anchored portion of the bolt length is found to be 3.5 times in faceplated bolt as compared to nonfaceplated bolts. Moreover, for non-faceplated bolts, the anchor length is about 1.9 m as compared to 2.0 m for faceplated bolts. These results clearly demonstrate the need for faceplated bolts for enhancing the performance of a 2 m rock bolt.

###### 6.1.3. Significance of Shear Stiffness of Bolt-Grout Interface on Rock Mass Behaviour

As the shear stiffness at the bolt grout interface increases, overall stresses and displacements around the tunnel boundary reduce, due to increase in shear resistance along the bolt length. In this study, three shear stiffness values are considered, namely. 10, and 100 GPa/m. Shear stiffness of 1 GPa/m may represent a weak grout mix whereas 100 GPa/m may represent very strong grout (probably resin mix). If shear stiffness is increased from 1 GPa/m to 100 GPa/m, the radial displacement at the boundary of the tunnel is reduced by 22% for models with faceplate bolted tunnel (Figure 12). Similar trend is observed in models without faceplate. It is also found that the faceplated bolts have significantly reduced the rock movement all along the bolt length and thus the stability of the rock mass have increased around the tunnel boundary.

Stress distributions around the tunnel boundary are shown in Figure 13 for with faceplate models. As the shear stiffness increases, the radius of elasto-plastic zone decreases. For these models having 9 bolts installed equidistant from each other, it is found that the radius of elasto-plastic zone (measured from the centre of the tunnel) is reduced from 2.438 m to 2.313 m once is increased from 1 GPa/m to 100 GPa/m. These results emphasize upon bolts as being primary support system would be effective if the bond (shear stiffness) between the bolt and grout is sufficiently stronger.

###### 6.1.4. Effect of Shear Stiffness on Bolt Performance

As before, bolt performance is analyzed in terms of bolt load and shear load developed in a rock bolt. As expected, grout having stronger bond strength or higher shear stiffness produces higher axial load and shear load in the bolt grout interface providing better support against the rock movement. The peak shear load is 0.042 MN/m for GPa/m and 0.9217 MN/m for GPa/m. Hence, the peak anchoring shear load is about 22 times higher for GPa/m as compared to GPa/m for faceplated bolts (Figure 14). In the same token, the peak axial tensile load increases by 290 kN for the same conditions. These results are the manifestation for the need of higher bond strength (shear stiffness) grout material, for better performance of bolts.

###### 6.1.5. Influence of Bolt Density on Rock Mass Behaviour

As mentioned earlier, bolt density is defined as the number of bolts installed per unit area of the tunnel boundary. Mathematically, for a circular tunnel of radius , bolt density, , is expressed as or where is number of bolts, = Transverse distance (m), = Longitudinal distance between two bolts (m).

In this case, since quarter tunnel is considered, the above equation can be expressed as for m and m. For maintaining symmetry, 2, 3, 5, and 9 bolts are included in the numerical models corresponding to 4, 8, 16, and 32 bolts for the entire tunnel boundary, respectively. It is obvious that, as the number of bolt increases, value of also increases. For no-support condition, is considered to be 0.

In general, if the frequency of bolting (bolt density) is increased then the displacements, the peak tangential stress, and the radius of elasto-plastic boundary may reduce. On the contrary, the overall radial stress will increase with increasing values of . Results obtained from numerical models support the above concepts as mentioned below. Figure 15 depicts the variation in the peak radial displacement at tunnel boundary with bolt density. From this figure, it is clear that (i) the peak radial displacement decreases with bolt density or number of bolt installed, and (ii) bolts with faceplate performs better as compared to bolts without faceplate as far as reduction in radial displacement is concerned. Figures 16 and 17 show the distribution of radial displacements for different bolt densities for faceplated and without faceplated model. It is found that in faceplated bolt model, the peak normalized radial displacement in percentage is reduced from 0.225 to 0.175 for GPa/m, if is increased from 0.0 to 2.55. Similar reduction in radial displacement is observed for without faceplated bolt model with increased bolt densities (Figure 16).

###### 6.1.6. Effect of Bolt Density on Bolt Performance

Axial tensile loads on the bolt are plotted with different bolt densities for without and with faceplate models (Figures 18 and 19, resp.). Load on the bolts has marginally decreased with increasing . This suggests that as the number of bolts increases, loads are distributed equally among all installed bolts. It may also be noticed that with faceplate, bolts carry almost twice the load as compared to without faceplate bolts, and hence the performance of the former is superior than the later.

#### 7. Conclusions

This paper proposes a methodology for the analysis of interaction between rock bolt and grout material and behaviour of rock mass using coupled FEM and FDM procedures. The FDM procedures are applied for evaluation of displacements, axial load of a bolt, and shear load per unit length at the bolt grout interface. The developed shear load at the interface is transformed as nodal reaction forces of rock element in FEM analysis. The FEM-FDM procedure has been validated with results obtained from pull-out tests. These procedures have been implemented in elasto-plastic finite element analysis using Hoek-Brown yield criterion. The developed technique is then applied for analyzing a bolted circular tunnel with variations in shear stiffness in the bolt-grout interface, bolt density, and faceplate options. This study shows the efficacy of the proposed method for practical applications. The specific conclusions of this study are summarized as follows.(i)The performance of faceplated bolts in terms of axial tensile load and shear load in the interface is superior to non-faceplated bolts. The study finds that for faceplate bolts the peak radial displacements in rock mass have reduced by 11% and 25% from in place of non-faceplated bolts and no-support conditions, respectively.(ii)The peak tangential stress is lower and the overall radial stress (confining stress) is higher if faceplates are attached in bolts.(iii)The radius of elasto-plastic zone reduces by 0.25 m from no-support conditions if bolting is considered along the 2 m radius tunnel boundary.(iv)The shear stiffness (bond strength grout material) plays an important role in bolt performance. The study finds that the peak axial tensile load increases from 110 to 400 kN if shear stiffness increases from 1 to 100 GPa/m. Similarly, shear load per unit length of the bolt also increases 22 times if the same condition exists.(v)Bolt density (frequency of bolting) around the tunnel also has an effect on the rock mass movements. It is found that the peak radial displacement in rock mass reduces by about 1.08 mm (with faceplated) and by 1.0 mm (without faceplated) from no-support condition if the number of bolts in a quarter tunnel of 2 m radius is increased from 0 to 9.

#### References

- C. Li and B. Stillborg, “Analytical models for rock bolts,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 36, no. 8, pp. 1013–1029, 1999. View at Google Scholar · View at Scopus - A. J. Hyett, M. Moosavi, and W. F. Bawden, “Load distribution along fully grouted bolts, with emphasis on cable bolt reinforcement,”
*International Journal for Numerical and Analytical Methods in Geomechanics*, vol. 20, no. 7, pp. 517–544, 1996. View at Google Scholar · View at Scopus - B. Indraratna and P. K. Kaiser, “Analytical model for the design of grouted rock bolts,”
*International Journal for Numerical & Analytical Methods in Geomechanics*, vol. 14, no. 4, pp. 227–251, 1990. View at Google Scholar · View at Scopus - H. Stille, M. Holmberg, and G. Nord, “Support of weak rock with grouted bolts and shotcrete,”
*International Journal of Rock Mechanics and Mining Sciences and Geomechanics*, vol. 26, no. 1, pp. 99–113, 1989. View at Google Scholar · View at Scopus - Y. Cai, T. Esaki, and Y. Jiang, “A rock bolt and rock mass interaction model,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 41, no. 7, pp. 1055–1067, 2004. View at Publisher · View at Google Scholar · View at Scopus - C. Carranza-Torres, “Analytical and numerical study of the mechanics of rockbolt reinforcement around tunnels in rock masses,”
*Rock Mechanics and Rock Engineering*, vol. 42, no. 2, pp. 175–228, 2009. View at Publisher · View at Google Scholar · View at Scopus - B. Brady and L. Lorig, “Analysis of rock reinforcement using finite difference methods,”
*Computers and Geotechnics*, vol. 5, no. 2, pp. 123–149, 1988. View at Google Scholar · View at Scopus - C. Carranza-Torres and C. Fairhurst, “The elasto-plastic response of underground excavations in rock masses that satisfy the Hoek-Brown failure criterion,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 36, no. 6, pp. 777–809, 1999. View at Publisher · View at Google Scholar · View at Scopus - C. Carranza-Torres, “Elasto-plastic solution of tunnel problems using the generalized form of the Hoek-Brown failure criterion,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 41, no. 1, pp. 629–639, 2004. View at Publisher · View at Google Scholar · View at Scopus - S. K. Sharan, “Elastic-brittle-plastic analysis of circular openings in Hoek-Brown media,”
*International Journal of Rock Mechanics and Mining Sciences*, vol. 40, no. 6, pp. 817–824, 2003. View at Publisher · View at Google Scholar · View at Scopus - D. Deb and S. O. Choi, “Singularity adjustment on generalized Hoek and Brown yield surface for elastic-plastic analysis,” in
*Proceedings of the 40th US Rock Mechanics Symposium*, Anchorage, Alaska, June 2005. - D. Deb,
*Finite Element Method Concepts and Applications in Geomechanics*, Prentice Hall, New Delhi, India, 2006. - G. N. Pande, G. Beer, and J. R. Williams,
*Numerical Methods in Rock Mechanics*, John Wiley & Sons, New York, NY, USA, 1990. - Rockscience Inc.,
*Phase2: Finite Element Analysis and Support Design*, Toronto, Canada, 1999. - Itasca Consulting Group, Inc., Minneapolis-2005. FLAC (Fast Lagrangian Analysis of Continua), Version-5.0.
- I. W. Farmer, “Stress distribution along a resin grouted rock anchor,”
*International Journal of Rock Mechanics and Mining Sciences and Geomechanics*, vol. 12, no. 11, pp. 347–351, 1975. View at Google Scholar · View at Scopus - E. Hoek, C. C. Torres, and B. Corkum, “Hoek and Brown failure criterion,” in
*Proceedings of the 5th North American Rock Mechanics Symposium and 17th Tunneling Association of Canada Conference (NARMS-TAC '02)*, pp. 267–271, 2002. - J. Chakrabarty,
*Theory of Plasticity*, vol. 1, McGraw-Hill, New York, NY, USA, 1987. - D. R. J. Owen and E. Hinton,
*Finite Elements in Plasticity: Theory and Practice*, Pineridge Press, 1980. - T. Zhen Yu and C. Jie Xian, “Behavior of bolting as tunnel support,” in
*Proceedings of the International Symposium on Rock Bolting*, O. Stephansson, Ed., pp. 87–98, Abisko, Sweden, 1983. - R. Ulusay and J. A. Hudson,
*The Complete ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 1974–2006*, ISRM Turkish National Group, Ankara, Turkey, 2007.