Research Article  Open Access
Debasis Deb, Kamal C. Das, "BoltGrout Interactions in Elastoplastic Rock Mass Using Coupled FEMFDM Techniques", Advances in Civil Engineering, vol. 2010, Article ID 149810, 13 pages, 2010. https://doi.org/10.1155/2010/149810
BoltGrout Interactions in Elastoplastic Rock Mass Using Coupled FEMFDM Techniques
Abstract
Numerical procedure based on finite element method (FEM) and finite difference method (FDM) for the analysis of boltgrout interactions are introduced in this paper. The finite element procedure incorporates elastoplastic concepts with Hoek and Brown yield criterion and has been applied for rock mass. Boltgrout interactions are evaluated based on finite difference method and are embedded in the elastoplastic procedures of FEM. The experimental validation of the proposed FEMFDM 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 elastoplastic 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 closedform solution of elastic rock bolt in elastoplastic 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 elastoplastic rock mass with Hoek and Brown yield criterion [8–11]. Finite Element Method (FEM) and/or FDMbased 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 elastoplastic procedure of FEM. The proposed FEMFDM procedure is validated with the experimental results of pullout 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 FEMFDM 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 HoekBrown 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 MohrCoulomb 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 elastoplastic 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 elastoplastic 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. ElasticPlastic Analysis in Plane Strain Condition
The elastoplastic constitutive matrix for elasticperfectly 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, NewtonRaphson 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 BoltGrout 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 boltgrout interfaces. The bolt is represented by a onedimensional 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 crosssection 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 stressdisplacement 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 pullout tests conducted by Farmer [16] showed that loaddisplacement curve is nonlinear, which is evident from debonding between bolt rod and grout. However, in this present study, debonding or elastoplastic 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 forcefree 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 nth equation becomes Equations (16), (19), and (20) corporate linear simultaneous equations, which become a tridiagonal 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 BoltGrout Interactions into ElastoPlastic FEM Procedures
It is assumed that the shear force developed at the boltgrout 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 NewtonRaphson 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 GaussLegendre 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 4noded 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 NewtonRaphson method is applied to obtain solution for displacements. In this procedure, for each load step and for every iteration within a load step, HoekBrown 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 elastoplastic 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 outofbalance force vector, as given in righthand side of (22) is applied as input parameter to NR 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 FEMFDM procedure for boltgrout interactions in elastoplastic rock mass is given in Algorithm 1.

5. Experimental Validations
The proposed FEMFDM procedure has been validated using the results obtained from “pullout” tests of rock bolts conducted by Farmer [16]. In general, the “pullout” test is intended to measure the shortterm strength of rock bolt anchor [21]. The strength is measured by a pull test in which bolthead displacement is measured as a function of the applied load. In his pioneering work, Farmer [16] conducted pullout 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 pullout 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 (pullout 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 loaddisplacement 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 pullout loads ranging between 20 kN and the breaking load. Shear strains along bolt length of the pullout 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 debonding 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 nonlinear loaddisplacement curve.
Numerical study conducted on the same elastic rock medium suggests that during pullout 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 pullout loads. Shear strains are estimated at discretized points along bolt length. It is found that up to 80 kN pullout 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 pullout force exceeds 80 kN. The difference in results are contributed due to debonding between bolt rod and grout which is not considered in FEMFDM procedure ( being constant). From this study, it is clear that FEMFDM procedure proposed in this paper simulates pullout 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 FEMFDM 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 boltgrout 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 HoekBrown 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 nosupport model. Again, results of nosupport model are compared with analytical solutions having elastic perfectly plastic conditions in HoekBrown 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 NewtonRaphson method.

(a) Tunnel parameters.  
 
(b) Bolt and grout parameters.  

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 “nosupport” 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 nosupport 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 nosupport 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 nonfaceplated 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 BoltGrout 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 elastoplastic zone decreases. For these models having 9 bolts installed equidistant from each other, it is found that the radius of elastoplastic 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 nosupport 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 elastoplastic 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 FEMFDM procedure has been validated with results obtained from pullout tests. These procedures have been implemented in elastoplastic finite element analysis using HoekBrown yield criterion. The developed technique is then applied for analyzing a bolted circular tunnel with variations in shear stiffness in the boltgrout 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 nonfaceplated 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 nonfaceplated bolts and nosupport 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 elastoplastic zone reduces by 0.25 m from nosupport 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 nosupport 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
 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
 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
 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
 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 Site  Google Scholar
 C. CarranzaTorres, “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 Site  Google Scholar
 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
 C. CarranzaTorres and C. Fairhurst, “The elastoplastic response of underground excavations in rock masses that satisfy the HoekBrown failure criterion,” International Journal of Rock Mechanics and Mining Sciences, vol. 36, no. 6, pp. 777–809, 1999. View at: Publisher Site  Google Scholar
 C. CarranzaTorres, “Elastoplastic solution of tunnel problems using the generalized form of the HoekBrown failure criterion,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 1, pp. 629–639, 2004. View at: Publisher Site  Google Scholar
 S. K. Sharan, “Elasticbrittleplastic analysis of circular openings in HoekBrown media,” International Journal of Rock Mechanics and Mining Sciences, vol. 40, no. 6, pp. 817–824, 2003. View at: Publisher Site  Google Scholar
 D. Deb and S. O. Choi, “Singularity adjustment on generalized Hoek and Brown yield surface for elasticplastic analysis,” in Proceedings of the 40th US Rock Mechanics Symposium, Anchorage, Alaska, June 2005. View at: Google Scholar
 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., Minneapolis2005. FLAC (Fast Lagrangian Analysis of Continua), Version5.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
 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 (NARMSTAC '02), pp. 267–271, 2002. View at: Google Scholar
 J. Chakrabarty, Theory of Plasticity, vol. 1, McGrawHill, 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. View at: Google Scholar
 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.
Copyright
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.