Solution for strain-softening surrounding rock reinforced by grouted bolts

To evaluate the reinforcement effectiveness of the grouted rock bolts, a dimensionless parameter’s approach (i.e. bolt density β ) is adopted. Following the Mohr–Coulomb and generalized Hoek–Brown failure criteria, the concept of the neutral point for the grouted rock bolts is introduced, and the determination of the equivalent plastic zone is divided into three categories (i.e. minimal yielding, major yielding, and excessive yielding). The potential plastic zone is sufficiently subdivided into a large number of concentric annuli through the finite difference approximation of equilibrium and compatibility equations. A simple approach for circular tunnel in strain-softening surrounding rock incorporating the neutral point of the grout rock bolts is proposed to derivate the solutions of stress and displacement. The proposed approach is validated with the published results. With the increase in bolt density β , the plastic radius would reduce 10%–30%, and the radial displacements would decrease by 10%–40%.


Introduction
Grouted rock bolts have been widely used in all civil engineering fields because of its simple structure, convenient construction, remarkable reinforcement effectiveness, low cost, and engineering adaptability, especially in cracked and softening rock masses. The grouted rock bolts can effectively improve the strength and stabilization of the surrounding rock through convergence control. 1 Although most researchers have contributed to the stress and displacement of tunnels and have solved many technical problems in practice, the effects of strain-softening characteristic on the stress and displacement of a tunnel, reinforced with grouted rock bolts, have yet to be defined clearly. And the method of ground reaction curve of circular tunnel excavated in strain-softening rock mass was initially proposed by Brown et al. 2 based on finite difference method.
The effects of the grouted rock bolts on rock mass have been studied by many researchers. Stille 3 and Stille et al. 4 presented a closed-form elastic-plastic solution of the grouted rock bolts. Five different approaches of bolt performance are considered in this model, and the results agree well with the measured data. Indraratna and Kaiser 5 developed an analytical model in which the elastic-plastic concept was adopted and a proper interaction mechanism between the ground and the grouted bolts was considered. Peila and Oreste 6 presented a new convergence-confinement approach able to model a reinforced rock mass zone surrounding an opening zone with new improved properties. Oreste and Peila 7 presented a new procedure for the computation of the convergence-confinement curves of a bolted tunnel. Fahimifar and Soroush 8 developed solutions for the stress and displacement of the tunnels that were excavated in softening rock mass and reinforced by the grouted rock bolts. Ranjbarnia et al. 9 proposed an analytical solution to determine the stress and displacement of a tunnel that was excavated in a brittle and strainsoftening rock mass and reinforced with active grouted rock bolts. This solution is based on the nonlinear peak strength criterion. Osgoui and Oreste 10 developed an elastic-plastic analytical solution of an axisymmetrical problem for a circular tunnel reinforced by the grouted rock bolts. This improved model was based on the approach of Indraratna and Kaiser, 5

and the nonlinear
Hoek-Brown (H-B) failure criterion was adopted. B Liu et al. 11 presented an analytical design method for the truss-bolt system in reinforcing underground fractured rock roofs in coal mines. Carranza-Torres 12 proposed the method of quantifying the mechanical contribution of rock bolts installed systematically around tunnels. Reinforcement of the effects of using rock bolt on rock mass with cross-flaws was explored by B Zhang et al. 13 X Li et al. 14 proposed an analytical model of shear behavior of a fully grouted cable bolt subjected to shearing to effectively evaluate the cable shear behavior and the influence of fully grouted cable bolts on joint shear resistance. A comprehensive study including laboratory tests and numerical modeling was performed to investigate factors of rock bolt fracture under varying loading conditions by H Kang et al. 15 In the analyses of engineering problems, the Mohr-Coulomb (M-C) and H-B failure criteria have been used widely by most scholars such as Carranza-Torres, 16 Park and Kim, 17 and Yang and Yin. 18 Lee and Pietruszczak 19 developed a numerical procedure for calculating the stress and radial displacements of a circular tunnel excavated in strain-softening rock mass, which obeyed the M-C and generalized H-B failure criteria. Park et al. 20 proposed the pressure-induced inflatable pipe method. The concept of cavity expansion for tunnel reinforcement was used. The interaction between tunnel supports and ground convergence was also investigated analytically and numerically in this approach. Zou and Li 21 proposed an improved numerical approach to analyze the stability of the strainsoftening surrounding rock with the consideration of the hydraulic-mechanical coupling and the variation in elastic strain in the plastic region. In those literatures, the analytical and numerical solutions based on the plane strain assumption are discussed with different constitutive models (e.g. elastic-plastic model, elastic-brittle-plastic model, and strain-softening model) combined with the linear M-C or generalized H-B failure criteria. Although the solutions of the strain-softening surrounding rock based on M-C or H-B failure criterion are obtained by most researchers, the effect of the grouted rock bolts and the position of the neutral point was seldom considered.
The main objective of this article is to explore solutions for the stress and displacement of circular openings excavated in strain-softening surrounding rock. In the presented solutions, a numerical stepwise procedure that considers the reinforcement effectiveness and the neutral point of the grouted rock bolts is adopted and modified, where the deterioration of strength in the plastic region is considered. The improved model which considers the strain-softening behavior and is compatible with M-C and generalized H-B failure criteria was proposed based on model of Indraratna and Kaiser. 5 Theory and methodologies

Statement of the problem
As shown in Figure 1, in a continuous, homogeneous, isotropic, initially elastic rock mass, a circular opening with a radius (a) is excavated. The cross section of the tunnel is subjected to uniform hydrostatic pressure (s 0 ). After excavating, an internal support pressure (p in ) acts on the tunnel wall surface in the radial direction. Grouted rock bolts are installed systematically over the periphery of the opening. In the system of the cylindrical polar coordinates, the stress and displacement of the surrounding rock are only related to radius, so the gravity field can be ignored. During excavation of tunnel, p in is gradually reduced. When the internal support pressure (p in ) is less than a critical support pressure, plastic region will appear around the surrounding rock. After this initial yielding, the rock strength drops gradually following the post-yield softening behavior.
The strain-softening model is used, and the numerical stepwise procedure is improved for the solution for stress and strain in the plastic region.
The stress-strain relationship of the strain-softening rock mass is shown in Figure 2. Figure 2 displays the idealized relationships of s 1 À s 3 versus e 1 . Both the decrease in strength from the peak and the continued deformation at the residual region are accompanied by plastic dilation.
Several assumptions below have been made to determine the solution of stress and displacement of the strain-softening surrounding rock incorporating the reinforcement effectiveness and the neutral point of the grouted rock bolts. The rock mass around the circular opening is regarded as an isotropic, continuous, and permeable medium which obeys the M-C or generalized H-B failure criteria under the plane strain condition. The strain-softening constitutive model that follows a non-associated flow rule is employed for formula derivation. The elastic strain in the plastic and softening regions of the surrounding rock obeys Hooke's law.
Considering the reinforcement effectiveness and the neutral point of the grouted rock bolts, some assumptions are also needed. Bolt density parameter and neutral point As shown in Figure 3, a circular tunnel is reinforced by grouted rock bolts. The bolt density parameter proposed by Indraratna and Kaiser 5 is adopted to describe the global intensity parameter of the reinforced surrounding rock as follows where b is dimensionless. It reflects the relative density of the rock bolts and takes into consideration the shear stress on the bolt surface, which opposes the rock mass displacements near the tunnel wall. The diameter of grouted rock bolts is d and the tunnel radius is a. l 0 is the friction factor which is analogous to the coefficient of friction. The number of the grouted rock bolts is presumably distributed uniformly on the tunnel surface in both the circumferential and longitudinal directions to maintain model symmetry. The horizontal distance between the two adjacent bolts, in the inner surface of the tunnel, is denoted by S T . The longitudinal space along the axis of the tunnel is represented by S L .
The shear stress is related to the first derivative of the axial stress. Hence, one point where the value of t Z is 0 is defined as the neutral point that exists where the axial stress attains the maximum, as shown in Figure 3. The relationship between the shear stress and the axial stress is presented by The shear stress is characterized by the division of the rock bolts, a pick-up length and an anchor length, on either side of the neutral point. By the division of the rock bolts, pick-up side and anchor side, the shear stress is characterized as presented in Figure 4. The location of the neutral point is given by Indraratna and Kaiser 5 r is also the anchor length.

Three categories of the equivalent plastic zone
The plastic zone around the surrounding rock reinforced by the grouted rock bolts is redefined as the equivalent plastic zone. Because the location of the plastic zone is related to the neutral point and the bolt length, three cases (as shown in Figure 5) are analyzed and calculated for the determination of the equivalent plastic zone radius (R) and corresponding classifications are shown below.

Failure criterion
The M-C failure criterion is given by where s 1 is the major principal stress; s 3 is the minor principal stress; Y and N are the strength parameters of the rock mass (Y = 2c cos f=(1 À sin f), where c and f are the cohesion and internal friction angle, respectively.
The following generalized H-B failure criterion is adopted by where s 1 and s 3 are the major and minor principal stresses, respectively. s c is the uniaxial compressive strength of the rock mass. a, m, and s are the strength parameters of the generalized H-B failure criterion. The strength and deformation parameters of the strain-softening rock mass are evaluated based on plastic deformation and are controlled by the plastic deviatoric strain 19,22 where e p 1 and e p 3 are the major and minor plastic strains, respectively.
The bilinear function that adopts the plastic shear strain to describe the physical parameters of the surrounding rock mass 19,22 where v represents a strength parameter, such as c, u, m, s, and a; g p r is the critical deviatoric plastic strain from which the residual behavior is observed first, and which should be identified through experimentation. The subscripts p and r represent the peak and residual values, respectively.

Stress and strain solutions in elastic region
The solutions of stress, strain, and displacement in the elastic zone considering the reinforced rock are given by Indraratna and Kaiser 5 as follows  The relationship of the radial stress and the tangential stress at the interface between the elastic rock and the equivalent plastic zone is presented by Combining equations (4) and (11) leads to where s rR is the radial stress at the interface between the elastic zone and the equivalent plastic zone.
Stress and strain solutions of excessive yielding (R Ã .(a + L)) Zone I: a + L\r\R Ã In zone I, the plastic zone is beyond the reinforced zone. The stress equilibrium equation of an element near an unsupported opening can be represented by In a bolted element, as shown in Figure 6, the stress equilibrium equation is given by Equations (13) and (14) describe the stress equilibrium condition of the unreinforced and reinforced segments, respectively. If the terms N (1 + b) and Y (1 + b) are replaced by the equivalent strength parameters N Ã and Y Ã , respectively, equation (14) for the bolted parts can be simplified to where If the bolt density parameter is introduced into the stress equilibrium equation, the radial stress at the interface between the elastic rock and the equivalent plastic zone can be modified by Considering the reinforcement effectiveness of the grouted rock bolts, the residual strength parameters of rock mass can be represented by The total plastic region is presumably divided into n connect annuli which are bounded by two annuli of the radius r (iÀ1) and r (i) . The radius of the first ring is r (0) = R which is at the interface between the elastic region and the equivalent plastic zone. To satisfy the equilibrium condition, the thickness of the annuli is determined automatically by the numerical process. In Figure 7, the ith annulus which is determined by the outer radius (r (iÀ1) ) and the inner radius r (i) adjacent to the elastic region is given. The slight increment of the radial stress is regarded as the known increment. So the (n 2 1) iterations of brittle-plastic analysis are conducted on the inner plastic annulus until the residual strength is reached. s r on both the inner and outer boundaries of the plastic zone are assigned, and s r decreases monotonically from s rR at r = R to p in at r = r 0 . The increment of radial stress is as follows So the radial stress components for ith radius can be For a large magnitude n, the corresponding tangential stress is given by where N Ã (i) and Y Ã (i) are the equivalent strength parameters.
So the tangential stress components can be represented by Under the plane strain condition, the elastic strain increments can be obtained using Hooke's law as follows The stress equilibrium differential equation that considers the reinforcement effectiveness of the grout rock bolts for the ith annulus is derived using equation (15) and expressed as equation (22) s where The normalized inner radius r (i) = r (i) =R can be expressed as where X = (((1 À N (i) ) s (i) + Y (i) )=( s (i) )). The compatibility equation can be written in the general form as follows where e r = du=dr and e u = u=r. The total strain in the plastic zone that is the sum of the elastic strain and the plastic strain is given by Combining equations (22) and (23) leads to where A 1 = (1=Dr (i) ) + (1 + k (iÀ1) (1= r (i) )), A 2 = (De e u(i) = Dr (i) ) À((De p u(iÀ1) À De p r(iÀ1) )=( r (i) )), A 3 = ((e e u(i) À e e r(i) )= ( r (i) )), and k (iÀ1) = ((1 + sin u (iÀ1) )=(1 À sin u (iÀ1) )).
Using equation (6), the plastic deviatoric plastic strain can be expressed by Through the relationship between strain and displacement, the displacement at ith annulus can be obtained as Thus, the stress and displacement solutions of the reinforced rock in this plastic region are obtained.
Zone II: r\r\a + L This is the middle plastic zone which is confined between the neutral point and the bolt end. Consequently, the rock mass is subjected to the negative shear stress of the rock bolts. So the rock mass has weak values of the strength parameters as follows If N 0 is greater than the residual strength of rock mass N Ã r , the rock mass strength first drops to N 0 and then continues to decrease until N Ã (i) = N Ã r . If N 0 is less than the residual strength of rock mass, as this plastic zone becomes a residual plastic zone, the rock mass strength no longer changes. So the rock mass strength parameters maintain the residual strength N r .
The first annulus in zone II is the interface between the plastic zones I and II, and the boundary conditions of stress and strain in this region can be obtained. Then as with zone I the stress and strain can be calculated by stepwise method.

Zone III: a\r\r
This yielded zone is confined by a range between the tunnel boundary and the neutral point. It falls within the pick-up length of the rock bolts and is stabilized by the positive shear stress. The residual plastic zone is reinforced and its strength parameters are effectively increased, and the strength parameters maintain the residual strength as follows The boundary condition of the stress and strain in region II is used as the initial condition of region III. Then the calculation steps are the same as region II.
Since e u = u=r, the displacement at each ring can be obtained by The procedure explained above can be repeated for n times until s r(i) reaches the internal support pressure. Therefore, we can obtain the equivalent plastic radius of zone III as follows Stress and strain solutions of major yielding (r\R Ã \(a + L)) Zone I: r.a + L Zone I is in elastic state and beyond the rock bolts. The stresses in elastic zone are given by where s L = s 0 ½1 À ((R)=(a + L)) 2 + s rR ((R)=(a + L)) 2 .
Zone II: R Ã \r\a + L This is the inner elastic zone within the reinforced zone. The stress in this zone is given by Zone III: r\r\R Ã This is the outer plastic zone beyond the neutral point. The strength parameters of rock mass begin to decrease from maximum values N Ã and Y Ã . If N 0 is greater than the residual strength of rock mass N Ã r , the rock mass strength first drops to N 0 and then continues to decrease until N Ã (i) = N Ã r . If N 0 is less than the residual strength of rock mass, as this plastic zone becomes a residual plastic zone, the rock mass strength no longer changes. Then the rock mass strength parameters maintain the residual value N Ã r . The stress and strain of rock mass are solved by the same method explained in section ''Stress and strain solutions of excessive yielding (R Ã .(a + L)).''

Zone IV: a\r\r
This yielded zone is confined between the neutral point and tunnel boundary, and the same approach in section ''Stress and strain solutions of excessive yielding (R Ã .(a + L))'' can be used.
Stress and strain solutions of minimal yielding (R Ã \r\(a + L)) Zone I: r.(a + L) This is the outermost elastic zone, and the stress is given by where s L = s 0 ½1 À ((r)=(a + L)) 2 + s r ((r)=(a + L)) 2 .
Zone II: r\r\(a + L) The anchor length of the rock bolts is contained in this elastic zone. The stress can be represented by Zone III: R Ã \r\r This elastic zone is confined to the end of the pickup length of the rock bolts. The stress is given by Zone IV: a\r\R Ã This region is the plastic zone. The stress and strain of rock mass are solved by the same method explained in section ''Stress and strain solutions of excessive yielding (R Ã .(a + L)).''

Stress and strain solutions in elastic region
Taking excessive yielding (R Ã .a + L) as an example, its equivalent plastic radius exceeds the length of the rock bolts. In the plastic zone, the strength parameters of the reinforced rock mass are increased to equivalent parameters and are as follows The relationships of the radial and tangential stresses at the interface between the elastic rock and the equivalent plastic zone are presented by Considering the reinforcement effectiveness of the grouted rock bolts, the radial stress at the interface between the elastic rock and the equivalent plastic zone can be modified as follows Equation (41) cannot be solved linearly, and the value of s rR can be solved by Newton iteration method.
The stress, strain, and displacement in the elastic zone are represented as follows, respectively

Stress and strain solutions in plastic region
Similar to the solutions based on the M-C failure criterion mentioned in section ''Stress and strain solutions of excessive yielding (R Ã .(a + L)),''s r on both the inner and outer boundaries of the plastic zone are known a priori and s r decreases monotonically from s rR at r = R Ã to p in at r = a.
The stress equilibrium differential equation that considers rock bolts' effectiveness for the ith annulus is expressed as equation (45) s r(i) À s r(iÀ1) where r (i) = r (i) =R Ã , r (i) = (1=2)(r (i) + r (iÀ1) ), and . The normalized inner radius (r (i) = r (i) =R Ã ) can be given by The increment of the tangential strain is represented by where . The displacement at each ring can be obtained using the relationship between the strain and displacement.
We can use the same method in section ''Stress and strain solutions in plastic region'' to calculate the stress and displacements of another two categories: the major yielding (r\R Ã \a + L) and the minimal yielding (R Ã \r\a + L).

Validations
To validate the accuracy of the proposed approach based on the M-C and the generalized H-B failure criteria, the results of the proposed approach for g p = 0 are compared with those by Indraratna and Kaiser 5 and Osgoui and Oreste, 10   As is demonstrated in Figure 8, for the cases of g p = 0, the results of the proposed approach are consistent with those presented by Indraratna and Kaiser 5 and Osgoui and Oreste, 10 because the elastic-brittleplastic model (g p = 0) is just a special case of strainsoftening model. Hence, the solutions of Indraratna and Kaiser 5 and Osgoui and Oreste 10 are the special cases of the proposed approach based on the M-C and generalized H-B failure criteria, respectively.

Numerical analysis and discussion
Comparison with elastic-brittle-plastic model To study the effect of strain-softening in combination with the grouted rock bolts on the stress and displacement, the stress and strain curves are plotted to compare the strain-softening model with the elastic-brittle-plastic model proposed by Indraratna and Kaiser 5 and Osgoui and Oreste, 10 respectively. Figures 9-11 depict the differences between the predicted results of the strain-softening model and elasticbrittle-plastic model with different bolt densities. It is shown from Figures 9 and 11 that stresses presented by the proposed model are larger than those based on the elastic-brittle-plastic model at the same point in the tunnel. Moreover, as is shown in Figures 14 and 15, the plastic radius of the proposed model would be reduced by 10% and 12%, respectively, for b = 0.145 and 0.291. The probable reason is that the strength parameters of strain-softening model are larger than those of the elastic-brittle-plastic model in some annuli of the plastic zone. As can be seen from Figure 9(c) that with the increase in the bolt densities, the stress curves of the proposed model are more agreeable with the curves based on the ideal elastic-plastic model in comparison with the elastic-brittle-plastic model.
It can be found from Figure 10(a) and (b) that the tangential strains presented by the proposed model are  It can be noted that the displacement in Figure 12 presented by the proposed model is smaller than those presented by Osgoui and Oreste. 10 For example, the radial displacement decreases from 18.2 to 15.3 mm if b = 0. The reduction is due to the fact that the strength parameters of strain-softening model are larger than those of the elastic-brittle-plastic model in some annuli of the plastic zone. Therefore, this model is more accurate than the elastic-brittle-plastic model. The tendency that the displacement of the reinforced surrounding rock decreases with the b increasing is provided in Figure 14. For example, the displacement reduction is 25% as b increases from 0 to 0.145 for s 0 = 12 MPa. The displacement of the reinforced surrounding rock corresponding to a larger b is smaller. One change that the displacement decreases about 55% as b increases from 0 to 0.291 for p in = 0.4 MPa can been seen. As shown in Figure 14, the radial displacement decreases with the increasing internal support pressure. So the bolt density b has a significant effect on the surrounding rock reinforcement.

Effects of the bolt density parameters
Based on the H-B failure criterion, two cases are carried out corresponding to (1) b = 0 and g p = 0.0047   The results are displayed in Figure 15(a)-(c). It can be seen from Figure 15(a) that stress increases with the increasing bolt densities at the same point. The plastic radius would be reduced correspondingly for a larger b. For example, the plastic radius of the proposed model would be reduced from 4.35 to 3.75 m if b increases from 0 to 0.132. Therefore, the effects of the bolt density parameters are significant in the surrounding rock reinforcement. With different hydrostatic pressures, displacements in the tunnel wall and the equivalent plastic radius vary with different bolt density parameters as shown in Figures 15(b) and 15(c). It can be seen that the greater the bolt density parameters, the smaller the radial displacements in the tunnel wall and equivalent plastic radius. For example, the displacement decreases from 25 to 15.5 mm as b increases from 0 to 0.325 for s 0 = 14 MPa.
As shown in Figure 15(d), the radial displacement in the tunnel wall decreases when the increasing internal support pressure increases. The displacement of the reinforced surrounding rock corresponding to a larger b is smaller. For instance, the displacement decreases about 30% as b increases from 0 to 0.25 for p in = 0 MPa. Therefore, the effect of the bolt density on the surrounding rock reinforcement is great.

Effects of the strain-softening parameters
To identify the effect of g p on the stress, strain, and displacement field under M-C and H-B failure criteria, three different g p are analyzed, respectively. Meanwhile, several examples are taken to explain the effect.
It can be seen from Figures 16(a) and 17(a) that the stress increases with the increase in g p at the same point. The plastic radius would be reduced correspondingly for a larger g p . Moreover, the plastic radius would be reduced by 10% for b = 0.291. Figures 16(a) and 17(a) and (b) illustrate the reduction in tangential strain and radial displacement for different g p , respectively. Figure 17(b) reflects that the radial displacement decreases from 14 to 11.5 mm if g p increased from 0 to 0.008.
Hence, it would be easy to note that g p has a significant effect on the surrounding rock reinforcement.
The effects of g p on the stress are significant, and the equivalent plastic radius decreases with the increase in g p . For example, the equivalent plastic radius for g p = 0 decreases from 4.05 to 3.5 m if g p is equal to 0.008.

Comparison with in situ monitoring data
General situation of engineering. Sanchahe Tunnel is located in the slope transition zone between Yunnan Guizhou Plateau and Guangxi hilly. The pile number of the left tunnel of Sanchahe Tunnel is ZK74 + 277-ZK75 + 245, the surrounding rock of the tunnel is mainly the plastic gravel containing silty clay, the length is 968 m, and the maximum buried depth is about 91.3. The zone ZK74 + 277-ZK75 + 615 is the key monitoring object; because the surrounding rock is soft and the stability is pretty poor, the surrounding rock properties are poor, such as joint and fissure, rock mass fragmentation, and the shape of the loose soil. The vertical section profile of ZK74 + 277-ZK75 + 615 is shown in Figure 18.
Measurement of radial displacement. The ground deformation of this tunnel is monitored and measured by multipoint borehole extensometers. Multi-point borehole extensometers are installed within the tunnel ground Based on the M-C failure criterion, the results of the improved model that considering the bolt density b and the neutral point of the grout rock bolts and the strainsoftening model are compared with the monitoring data of Sanchahe Tunnel. As shown in the figure, for the relationships between the radial displacements and the radius for the case of b = 0:18, when r=a.2, the calculation result is slightly larger than the monitoring result, which represents that the model is reliable. However, when r=a\2, the monitoring result is much smaller than the calculation result, most probably because of the following: s 0 = 15 MPa and P in = 1:5 MPa are given by substantial engineering experience; not the real value, but there are some temporary supports for tunnel wall to reduce the radial displacement; the deformation of multi-point borehole extensometers is not coordinated with the surrounding rock deformation.
Based on the generalized H-B failure criterion, for the relationships between the radial displacements and the radius for the case of g p = 0:08, when r=a\2, the tendency of the calculation curve is nearly the same as the monitoring curve, which represents that the model is also reliable. Comparing Figure 20 with Figure 21, the curve of the calculation result based on the M-C failure criterion is more consistent with the monitoring data curve. This is because the strength parameters of the generalized (H-B) failure criteria are estimated approximately, which can cause slight differences with the actual value.

Conclusion
Incorporating the strain-softening characteristic of the surrounding rock and the neutral point of the grouted rock bolts, the equivalent plastic zone is divided into three categories (i.e. minimal yielding, major yielding, and excessive yielding). The potential plastic zone is sufficiently subdivided into a large number of  concentric annuli; the bolt density parameters are introduced, the solutions of stress and displacement of circular tunnel are presented, and the corresponding calculation approach is proposed as well. Based on the elastic-brittle-plastic model, the proposed model is validated by the special case as g p = 0. Moreover, the stress presented by the proposed model is larger than those presented by elastic-brittle-plastic model at the same point. The comparison results show that the radial displacements in the tunnel wall and the plastic radius will be effectively reduced with the increase in the bolt density parameter b and strain-softening parameter g p . In practical tunnel engineering, the proposed approach can be widely used to evaluate the reinforcement effectiveness due to its simple practicality. In our   future work, we will compare our results to in situ monitoring data. Furthermore, the model can be extended to a similar problem in other projects, such as slope engineering, foundation pit engineering, and mining engineering.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The authors are grateful to the 973 Program (2013CB036004) and National Natural Science Foundation of China (no. 51208523).