# The Analytical Approach to Evaluate the Load-Displacement Relationship of Rock Bolts

Authors: Jianhang Chen; Shengli Yang; Hongbao Zhao; Junwen Zhang; Fulian He; Shuaifeng Yin
Journal: Advances in Civil Engineering
 (2019)

---

## Abstract

Fully grouted rock bolts are widely used in civil engineering and mining engineering, playing a significant role in keeping the stability and safety of excavations. In this paper, the load transfer mechanism of fully grouted rock bolts was studied with an analytical model. A trilinear model was used to depict the bond-slip behaviour of the bolt/grout interface. The displacement of the confining medium was involved in this analytical model. Then, the shear stress propagation along the bolt/grout interface was analysed in the elastic, elastic-softening, elastic-softening-debonding, softening-debonding, and debonding stages. Experimental pull-out tests were used to validate this analytical model. There was a good correlation between experimental and analytical results. A parametric study was conducted to evaluate the influence of Young’s modulus of the confining medium, the shear strength of the bolt/grout interface, and the residual shear strength of the bolt/grout interface on the load transfer performance of rock bolts. The results show that increasing Young’s modulus of the confining medium was beneficial for improving the load transfer performance of rock bolts. However, once Young’s modulus of the confining medium was beyond a critical limit, it had marginal effect on the peak load of rock bolts. Furthermore, increasing the shear strength of the bolt/grout interface and the residual shear strength of the bolt/grout interface led to rising of the peak load of rock bolts. However, compared with the residual shear strength of the bolt/grout interface, increasing the shear strength of the bolt/grout interface had more apparent effect in improving the peak load of rock bolts.


---

## Body

## 1. Introduction

A rock bolt is a rock tendon which can be either a rod or a strand, installed in a borehole to reinforce the surrounding rock mass and soil [1]. Rock bolts have been used in civil and mining for a long time [2]. Numerous field applications and computer simulation have proved that rock bolts play a significant role in improving the internal strength of rock mass and enhancing the safety of underground excavations [3–6].The rock bolt reinforcement system is basically composed of four components: the surrounding confining medium, e.g., rocks or soils, the rock reinforcement element, the internal fixture, and the external fixture [7–9]. It is found that all these four components have an influence on the performance of a reinforcement system.To better understand the performance of the reinforcement system, Thompson et al. [10] classified the rock reinforcement system as three types: continuously mechanically coupled (CMC), continuously frictionally coupled (CFC), and discrete mechanically or frictionally coupled (DMFC). This study deals with fully grouted rock bolts, and based on this classification, fully grouted rock bolts belong to the CMC system.For the CMC system, cementitious grout or polyester resin can be used to bond rock bolts with the surrounding confining medium [11–13]. Via the bonding material, the load from exterior unstable rock mass can be transferred to interior stable rock mass. It should be mentioned that in the in situ applications, rock bolts can encounter axial, shear, and combined axial and shear loading conditions [14]. However, this paper focused on the axial performance of rock bolts.Failure of the rock bolting system has different modes. Jeremic and Delaire [15] summarised that failure of the rock bolting system can be classified as bond failure of the bolt/grout interface, bond failure of the grout/rock interface, failure of the grout annulus, and failure in the surrounding rock mass. Then, Hutchinson and Diederichs [16] added that rupture of the rock tendon also occurred in the in situ scenario. Nevertheless, among these failure modes, bond failure of the bolt/grout interface and the grout/rock interface are more common to encounter [17].Jalalifar [18] stated that when fully grouted rock bolts were loaded axially, the shear stress along the bolt/grout interface was larger than the shear stress along the grout/rock interface because the bolt/grout interface had a smaller contact area. When the grout and the rock had similar strength, failure was likely to occur along the bolt/grout interface. Previous experimental pull-out tests also showed that bond failure of the bolt/grout interface occurred in the tests [19, 20]. Therefore, this paper deals with the bond failure of the bolt/grout interface.Gambarova [21] stated that the bonding performance of the bolt/grout interface mainly relies on chemical adhesion, mechanical interlock, and friction. Chemical adhesion plays a negligible role on the bonding performance and can be easily removed after a small relative movement between the bolt and grout [22].To understand the load transfer performance of rock bolts, extensive research has been conducted via experimental pull-out tests [23–25]. Compared with that, much less work has been conducted with analytical modelling. Hawkes and Evans [26] conducted the pioneering work in analytically modelling the load transfer of rock bolts. Specifically, an analytical model was proposed to evaluate the load distribution along a fully grouted rock bolt. However, this assumption was only valid in the coupling stage. In the modelling process, a linear relationship between the axial stress in the bolt and the shear stress along the bolt/grout interface was assumed. However, this assumption was credible when the bolt was fully coupled with the surrounding medium. With the axial load in the bolt increasing, debonding occurred along the bolt surface, and this assumption may not be valid again. Farmer [27] deduced an analytical model to study the load distribution along a rock bolt. In this analytical model, it was assumed that the shear stress along the bolt/grout interface was linear with the grout shear modulus. Therefore, this model was only applicable when the bolt/grout interface underwent elastic behaviour. Once debonding occurred along the bolt/grout interface, the credibility of this analytical model decreased. Aydan et al. [28] deduced an analytical solution to study the load transfer of rock bolts. However, no experimental work was conducted to validate the credibility of the model. Rajaie [29] studied the load transfer of cable bolts, and an energy approach was adopted to analyse the loading process of cable bolts. There was a good correlation between laboratory results and analytical results. However, this model mainly focused on the peak load calculation, and no equation was deduced to calculate the load variation with displacement in the pull-out process. Therefore, this model cannot be used to study the load-displacement relationship of cable bolts.Yazici and Kaiser [30] proposed a bond strength model to study the load transfer of plain cable bolts. Later, this model was modified to study the influence of stress change on load transfer of cable bolts [31]. This model treated the bolt/grout interface as a zigzag geometry, similar to a rock joint. For the rock joint, strength of the rock is usually weak along two sides. However, for the cable bolting system, strength of the steel is much higher along the bolt side. Following that, Hyett et al. [32] developed a constitutive model to study the load-displacement performance of cable bolts. Two different boundary conditions were considered: constant normal stress and constant normal stiffness. This model assumed that the shear stress along the bolt/grout interface was uniform. However, experimental work proved that when the bolt had a long embedment length, the shear stress along the bolt/grout interface was not uniform [22]. Therefore, this model was applicable to calculate the load-displacement relationship for short embedded bolts. However, it may not be suitable for long embedded bolts. A load distribution model was proposed by Hyett et al. [33]. This model can be used to study the load distribution in the embedment section and the anchor section. Moosavi [34] developed an analytical model for modified cable bolts. However, when modified cable bolts were used to reinforce rock masses, the modified geometry can be installed at different positions in the embedment section. And, the influence of the modified geometry position on the load transfer performance of cable bolts was not considered. Li and Stillborg [35] developed an analytical model to study the load distribution of rock bolts. The bolt/grout interface was assumed to be composed of four sections: completely decoupled, partially decoupled, softening, and elastic.Cai et al. [36] used a modified shear-lag model to study the load transfer of rock bolts. However, the softening behaviour of the bolt/grout interface was neglected. Xiao and Chen [37] assumed that the bolt/grout interface underwent elastic, softening, and debonding stages and developed a model for rock bolts. Ren et al. [38] deduced the load transfer process of rock bolts with an analytical model. Experimental work was used to validate the credibility of this analytical model. However, when deducing the analytical model, it was assumed that the shear slippage of the bolt/grout interface was equal to the axial displacement of the rock bolt. Then, the displacement of the surrounding confining medium was neglected. Consequently, this model cannot be used to study the influence of the surrounding confining medium on the performance of rock bolts. Blanco Martín et al. [39] proposed an analytical solution to evaluate the load-displacement performance of rock bolts. Jalalifar [40] evaluated the performance of rock bolts installed in an elasto-plastic rock mass. In this model, it was assumed that the shear force in the grout annulus had a linear relationship with the shear slippage. Therefore, this model was only valid when the bolt/grout interface underwent elastic deformation. After the bolt/grout interface started plastic deformation, this model was not valid again. Cao et al. [41] studied the influence of rib geometry on the bonding capacity of rock bolts. It was predicted that rock bolts with small rib face angles were more effective in hard rock conditions.Ma et al. [42] used a nonlinear equation to depict the bond-slip behaviour of the bolt/grout interface and studied the load transfer of rock bolts. Later, the model was modified to study the influence of free end slip on the performance of rock bolts [43]. However, in this model, an empirical equation was used to depict the bond-slip relationship of the bolt/grout interface. Therefore, the properties of the rock bolting system may modify the coefficients in the bond-slip model. Consequently, this model cannot be used to predict the influence of rock mass confinement on the performance of rock bolts. Chen et al. [13] used an analytical model to study the axial performance of cable bolts. However, the model was only appropriate for laboratory pull-out scenario. Li et al. [44] developed an analytical model to study the load transfer of modified cable bolts. Cable bolts with short embedment length were used to validate the model. However, in situ tests proved that the installed cable bolts can be very long for the field scenario [45]. And, whether this model can be used to study the performance of cable bolts with long embedment length was not discussed.The literature review shows that for the rock bolting system, when the shear behaviour of the bolt/grout interface was depicted with a bond-slip equation, different analytical models can be used, such as the exponential model [42] or the trilinear model [46]. Among them, the trilinear model considers the elastic, softening, and debonding behaviour of the bolt/grout interface and thus is more commonly used. However, when the trilinear model was incorporated into the rock bolting system to deduce the load-displacement relationship of rock bolts, the involvement of the surrounding confining medium was usually not considered.It should be mentioned that Vaculik et al. [47] conducted an analytical study to evaluate the load-slip relationship of the FRP-to-substrate joints with the bond-slip model. A framework was proposed to depict the performance of FRP-to-substrate joints. However, for the rock bolting system, cylindrical rock bolts are always used to reinforce the excavated rock masses, such as tunnels, roadways, and stopes [48, 49]. Therefore, this paper applies the method proposed by Vaculik et al. [47] into the fully grouted rock bolting system and further deduced the specific analytical solution for fully grouted rock bolts. Additionally, they evaluated the performance of FRP-to-substrate joints with the embedment length ranging from short to long. However, for fully grouted rock bolts, the embedment length is usually extremely long [45, 50]. For example, Cao et al. [51] summarised that the embedment length for a rock bolt can be around 3 m. Furthermore, when the cable bolt is used, the embedment length may range from 3 m to 15 m. Considering this, this paper especially focused on the long embedment length scenario for rock bolts.First, the relationship between the shear stress of the bolt/grout interface and the shear slippage was simulated with a trilinear model. Then, the pull-out process of rock bolts was analysed with five stages: elastic, elastic-softening, elastic-softening-debonding, softening-debonding, and debonding stages. After that, the analytical model was validated with experimental pull-out tests. Finally, a parametric study was conducted to evaluate the influence of Young’s modulus of the confining medium, the shear strength of the bolt/grout interface, and the residual shear strength of the bolt/grout interface on the load transfer performance of rock bolts.


## 2. Analytical Modelling

### 2.1. Bond-Slip Model

A trilinear model was used to depict the bond-slip behaviour of the bolt/grout interface. The trilinear model is composed of three parts: ascending section, descending section, and horizontal section [13, 38, 47], as shown in Figure 1. Specifically, with the shear slippage of the bolt/grout interface increasing from 0 to δp, the shear stress increases linearly from zero to the shear strength of the bolt/grout interface. Then, with δ increasing to δr, the shear stress decreases linearly to the residual shear strength of the bolt/grout interface. After δr, the shear stress remains constant and equals τr.Figure 1
A trilinear model [13, 38, 47].Therefore, the trilinear model can be expressed as [13, 38](1a)τ=τpδpδ,0≤δ≤δp,(1b)τ=τpδr−τrδpδr−δp−τp−τrδr−δpδ,δp<δ≤δr,(1c)τ=τr,δ>δr,where τ is the shear stress of the bolt/grout interface, Pa; τp is the shear strength of the bolt/grout interface, Pa; δp is the shear slippage when the shear strength of the bolt/grout interface reaches, m; τr is the residual shear strength of the bolt/grout interface, Pa; and δr is the shear slippage when the residual shear strength of the bolt/grout interface reaches, m.


### 2.2. Governing Equation

In previous research, the shear slippage of the bolt/grout interface was usually simplified as the axial displacement of the rock bolt [38, 42]. However, previous experimental work proved that the shear slippage of the bolt/grout interface was different from the axial displacement of the rock bolt because of the displacement of the confining medium [52]. Furthermore, Hyett et al. [33] found that the displacement of the confining medium had an apparent influence on the load transfer performance of rock bolts. Therefore, neglecting the displacement of the confining medium in the analytical model could not properly reflect the load transfer performance of rock bolts. Considering this issue, in this paper, the displacement of the confining medium was included.Since this study aims at studying the bond failure of the bolt/grout interface, the grout/medium interface was assumed to be intact in the pull-out process of rock bolts. Then, the shear slippage of the bolt/grout interface was defined as the relative movement between the rock bolt and the confining medium [53, 54]:(2)δx=ubx−umx,where δ(x) is the shear slippage of the bolt/grout interface at the position of x, m; ub(x) is the rock bolt displacement at the position of x, m; and um(x) is the confining medium displacement at the position of x, m.Differentiating equation (2) twice with respect to x leads to(3)d2δxdx2=d2ubxdx2−d2umxdx2.The axial stress in the rock bolt can be expressed with the displacement of the rock bolt. As for the confining medium, this is also applicable. Then, the following equations can be acquired [40]:(4)σbx=Ebdubxdx,(5)σmx=Emdumxdx,where Eb is Young’s modulus of the rock bolt, Pa; σb(x) is the axial stress in the rock bolt at the position of x, Pa; Em is Young’s modulus of the confining medium, Pa; and σm(x) is the axial stress applied on the confining medium, Pa.Substituting equations (4) and (5) into (3) leads to(6)d2δxdx2=1Ebdσbxdx−1Emdσmxdx.According to the mechanical equilibrium, the relationship between the axial stress in the bolt and the axial stress applied on the confining medium can be expressed as [55](7)σm=−σbπrb2Am,where Am is the cross-sectional area of the confining medium, m2 and rb is the rock bolt radius, m.Substituting equation (7) into equation (6) leads to(8)d2δxdx2=1Eb+πrb2EmAmdσbxdx.The relationship between the axial stress in the rock bolt and the shear stress along the bolt/grout interface can be seen in Figure2.Figure 2
Cross section view of a rock bolting system.Then, the following equation can be acquired [42]:(9)τx=rb2dσbxdx,where τ(x) is the shear stress of the bolt/grout interface at the position of x, Pa.Substituting equation (9) into equation (8) leads to(10)d2δxdx2−2rb1Eb+πrb2EmAmτx=0.It is defined that(11)λ2=2rb1Eb+πrb2EmAm.Then, equation (10) can be written as(12)d2δxdx2−λ2τx=0.Equation (12) is the governing equation of this analytical model, and it can be solved with the predefined shear stress equations.


### 2.3. Pull-Out Stages of Rock Bolts

According to previous research, the pull-out stages of the bolt/grout interface can be classified as the elastic stage, elastic-softening stage, elastic-softening-debonding stage, softening-debonding stage, and debonding stage [38, 53, 56]. Therefore, in the following contents, the pull-out process of a rock bolt was analysed with those five stages.


#### 2.3.1. Elastic Stage

When the pull-out load is small, the bolt/grout interface undergoes elastic deformation. And, the full bolt/grout interface is elastic, as shown in Figure3.Figure 3
The full bolt/grout interface is elastic in the elastic stage whereDb is the diameter of the rock bolt, m.Then, equation (1a) can be substituted into equation (12), leading to(13)d2δxdx2−τpδpλ2δx=0.It is defined that(14)λ12=τpδpλ2.Then, equation (13) can be written as(15)d2δxdx2−λ12δx=0.Equation (15) can be solved with boundary conditions. First, the axial stress in the bolt should be deduced. According to equation (2), the following equation can be acquired:(16)dδxdx=dubxdx−dumxdx.It can be further written as(17)dδxdx=σbxEb−σmxEm.Substituting equation (7) into equation (17) leads to(18)dδxdx=1Eb+πrb2EmAmσbx.Substituting equation (11) into equation (18), the axial stress in the bolt can be written as(19)σbx=2rbλ2dδxdx.The boundary conditions are(20)σbx=0=0,(21)σbx=L=Fπrb2,where F is the pull-out load of the rock bolt, N.Combining equations (19)–(21) with equation (15) leads to(22)δx=Fλ1δp2πrbτpsinhλ1Lcoshλ1x,(23)σbx=Fπrb2sinhλ1Lsinhλ1x,(24)τx=Fλ12πrbsinhλ1Lcoshλ1x.It is defined that at the loaded end, the shear slippage of the bolt/grout interface is equal to the pull-out displacement of the rock bolt [57]. Substituting x = L into equation (22), the pull-out load of the rock bolt can be expressed as(25)F=2πrbτptanhλ1Lλ1δpub,where ub is the pull-out displacement of the rock bolt, m.At the end of this stage, the shear stress at the position ofx = L increases to the shear strength of the bolt/grout interface. Therefore, substituting τ(x = L) = τp into equation (24) leads to(26)Fsof=2πrbτptanhλ1Lλ1,where Fsof is the pull-out load when the bolt/grout interface starts softening, N.


#### 2.3.2. Elastic-Softening Stage

In the elastic-softening stage, the full bolt/grout interface is composed of the elastic section and softening section, as shown in Figure4.Figure 4
The full bolt/grout interface is composed of the elastic section and softening section in the elastic-softening stage.It is defined thatas is the softening length. In the elastic-softening stage, the interface in the domain [0, L − as] is still elastic. One boundary condition is equation (20) and the other boundary condition is(27)τx=L−as=τp.Combining equations (9), (19), (20), and (27) with equation (15) leads to(28)δx=δpcoshλ1L−ascoshλ1x,(29)σbx=2τprbλ1coshλ1L−assinhλ1x,(30)τx=τpcoshλ1L−ascoshλ1x.Substitutingx = L − as into equation (29), the axial stress can be acquired:(31)σbx=L−as=2τprbλ1tanhλ1L−as.In the domain [L − as, L], the bolt/grout interface softens. Substituting equation (1b) into equation (12) leads to(32)d2δxdx2+λ2τp−τrδr−δpδx−λ2τpδr−τrδpδr−δp=0.It is defined that(33)λ22=τp−τrδr−δpλ2.Therefore, equation (32) can be further written as(34)d2δxdx2+λ22δx−λ22τpδr−τrδpτp−τr=0.The boundary conditions are equations (27) and (31). Combining equations (9), (19), (27), and (31) with equation (34) leads to(35)δx=λ1δpλ2tanhλ1L−assinλ2x+as−L−λ12δpλ22cosλ2L−as−x+τpδr−τrδpτp−τr,(36)σbx=2τprb1λ1tanhλ1L−ascosλ2L−as−x+1λ2sinλ2x+as−L,(37)τx=τpλ2λ1tanhλ1L−assinλ2L−as−x+cosλ2L−as−x.Substitutingx = L and σbx=L=F/πrb2 into equations (35) and (36), respectively, the pull-out displacement and pull-out load of the rock bolt can be acquired:(38)F=2πrbτp1λ1tanhλ1L−ascosλ2as+1λ2sinλ2as,(39)ub=λ1δpλ2tanhλ1L−assinλ2as−λ12δpλ22cosλ2as+τpδr−τrδpτp−τr.At the end of this stage, the shear stress at the position ofx = L reduces to τr. Therefore, substituting τ(x = L) = τr into equation (37) leads to(40)cosλ2as−λ2λ1tanhλ1L−assinλ2as=τrτp.For an infinite embedment length, tanh(λ1(L − as)) is equal to 1. Then, the softening length when the elastic-softening stage ends can be acquired:(41)as1=1λ2arcsinλ1λ12+λ22−arcsinτrλ1τpλ12+λ22,where as1 is the softening length when the elastic-softening stage ends, m.


#### 2.3.3. Elastic-Softening-Debonding Stage

Afterτ(x = L) reduces to τr, the bolt/grout interface enters the elastic-softening-debonding stage. In this stage, the full bolt/grout interface is composed of elastic section, softening section, and debonding section, as shown in Figure 5.Figure 5
The full bolt/grout interface is composed of the elastic section, softening section, and debonding section in the elastic-softening-debonding stage.Assuming that the debonding length isad, the bolt/grout interface in the domain [0, L − ad − as] is still elastic. Therefore, equations (28)–(30) are still valid if L is replaced with L − ad:(42)δx=δ1coshλ1L−ad−ascoshλ1x,(43)σbx=2τprbλ1coshλ1L−ad−assinhλ1x,(44)τx=τpcoshλ1L−ad−ascoshλ1x.As for the domain [L − ad − as, L − ad], the bolt/grout interface softens. Therefore, equations (35)–(37) are still valid if L is replaced with (L − ad):(45)δx=λ1δpλ2tanhλ1L−ad−assinλ2x+ad+as−L−λ12δpλ22cosλ2L−ad−as−x+τpδr−τrδpτp−τr,(46)σbx=2τprb1λ1tanhλ1L−ad−ascosλ2L−ad−as−x+1λ2sinλ2x+ad+as−L,(47)τx=τpλ2λ1tanhλ1L−ad−assinλ2L−ad−as−x+cosλ2L−ad−as−x.Substitutingτ(L − ad) = τr into equation (47) leads to(48)cosλ2as−λ2λ1tanhλ1L−ad−assinλ2as−τrτp=0.Equation (48) depicts the relationship between ad and as. It shows that in this stage, the softening length keeps varying with the debonding length increasing. Combining equations (45) and (48) together, the shear slippage of the bolt/grout interface at the position of x = L − ad can be acquired:(49)δx=L−ad=δr.Substitutingx = L − ad into equation (46), the axial stress can be acquired:(50)σbx=L−ad=2τprb1λ1tanhλ1L−ad−ascosλ2as+1λ2sinλ2as.For the domain [L − ad, L], the bolt/grout interface debonds and the shear stress equals τr:(51)τx=τr.Also, equation (1c) can be substituted to equation (12), leading to(52)d2δxdx2=λ2τr.Substituting equations (19), (49), and (50) into equation (52) leads to(53)δx=δr+12λ2τrx2−L−ad2+τpλ2λ1tanhλ1L−ad−ascosλ2as+τpλ2λ2sinλ2as−λ2τrL−adx+ad−L,(54)σbx=2τprb1λ1tanhλ1L−ad−ascosλ2as+1λ2sinλ2as+2τrrbx+ad−L.Substitutingx = L and σbx=L=F/πrb2 into equations (53) and (54), the pull-out displacement and load of the rock bolt can be acquired:(55)ub=δr+λ2ad2τr2+τpλ2λ1tanhλ1L−ad−ascosλ2as+τpλ2λ2sinλ2asad,(56)F=2πrbτpλ1tanhλ1L−ad−ascosλ2as+τpλ2sinλ2as+τrad.The maximum pull-out load occurs in this stage and can be acquired by solving equations (48) and (56) together.At the end of this stage, the elastic length disappears. Therefore,L − ad = as. Substituting it into equation (48) leads to(57)as2=1λ2arccosτrτp,where as2 is the softening length when the elastic-softening-debonding stage ends, m.


#### 2.3.4. Softening-Debonding Stage

After the elastic length disappears, the bolt/grout interface enters the softening-debonding stage. In the domain [0,L − ad], the bolt/grout interface softens. Substituting equations (9), (19), and (20) and τ(x = L − ad) = τr into equation (34) leads to(58)δx=τrδp−δrτp−τrcosλ2L−adcosλ2x+τpδr−τrδpτp−τr,(59)σbx=2τrrbλ2cosλ2L−adsinλ2x,(60)τx=τrcosλ2L−adcosλ2x.Substitutingx = L − ad into equations (58) and (59) leads to(61)δx=L−ad=δr,(62)σbx=L−ad=2τrrbλ2tanλ2L−ad.For the domain [L − ad, L], the bolt/grout interface debonds. Substituting equations (19), (61), and (62) into equation (52) leads to(63)δx=δr+λ2τr2x2−L−ad2+λ2τrλ2tanλ2L−ad−λ2τrL−adx+ad−L,(64)σbx=2τrrbtanλ2L−adλ2+x+ad−L.Substitutingx = L and σbx=L=F/πrb2 into equations (63) and (64), the pull-out displacement and load of the rock bolt can be acquired:(65)ub=δr+λ2τradad2+tanλ2L−adλ2,(66)F=2πrbτrtanλ2L−adλ2+ad.At the end of this stage, the debonding length equals the full embedment length. Therefore, substitutingad = L into equations (65) and (66), the pull-out displacement and load of the rock bolt at the end of the softening-debonding stage can be acquired:(67)ub0=δr+λ2τrL22,(68)F0=2πrbτrL,where ub0 is the pull-out displacement of the rock bolt at the end of the softening-debonding stage, m and F0 is the pull-out load of the rock bolt at the end of the softening-debonding stage, N.


#### 2.3.5. Debonding Stage

After the softening length disappears, the bolt/grout interface enters the debonding stage. Since the full bolt/grout interface undergoes debonding, the relationship between the pull-out load and displacement can be expressed as(69)F0−F=2πrbτrub−ub0.Substituting equations (67) and (68) into equation (69) leads to(70)F=2πrbτrL+δr+λ2τrL22−ub.


## 2.1. Bond-Slip Model

A trilinear model was used to depict the bond-slip behaviour of the bolt/grout interface. The trilinear model is composed of three parts: ascending section, descending section, and horizontal section [13, 38, 47], as shown in Figure 1. Specifically, with the shear slippage of the bolt/grout interface increasing from 0 to δp, the shear stress increases linearly from zero to the shear strength of the bolt/grout interface. Then, with δ increasing to δr, the shear stress decreases linearly to the residual shear strength of the bolt/grout interface. After δr, the shear stress remains constant and equals τr.Figure 1
A trilinear model [13, 38, 47].Therefore, the trilinear model can be expressed as [13, 38](1a)τ=τpδpδ,0≤δ≤δp,(1b)τ=τpδr−τrδpδr−δp−τp−τrδr−δpδ,δp<δ≤δr,(1c)τ=τr,δ>δr,where τ is the shear stress of the bolt/grout interface, Pa; τp is the shear strength of the bolt/grout interface, Pa; δp is the shear slippage when the shear strength of the bolt/grout interface reaches, m; τr is the residual shear strength of the bolt/grout interface, Pa; and δr is the shear slippage when the residual shear strength of the bolt/grout interface reaches, m.


## 2.2. Governing Equation

In previous research, the shear slippage of the bolt/grout interface was usually simplified as the axial displacement of the rock bolt [38, 42]. However, previous experimental work proved that the shear slippage of the bolt/grout interface was different from the axial displacement of the rock bolt because of the displacement of the confining medium [52]. Furthermore, Hyett et al. [33] found that the displacement of the confining medium had an apparent influence on the load transfer performance of rock bolts. Therefore, neglecting the displacement of the confining medium in the analytical model could not properly reflect the load transfer performance of rock bolts. Considering this issue, in this paper, the displacement of the confining medium was included.Since this study aims at studying the bond failure of the bolt/grout interface, the grout/medium interface was assumed to be intact in the pull-out process of rock bolts. Then, the shear slippage of the bolt/grout interface was defined as the relative movement between the rock bolt and the confining medium [53, 54]:(2)δx=ubx−umx,where δ(x) is the shear slippage of the bolt/grout interface at the position of x, m; ub(x) is the rock bolt displacement at the position of x, m; and um(x) is the confining medium displacement at the position of x, m.Differentiating equation (2) twice with respect to x leads to(3)d2δxdx2=d2ubxdx2−d2umxdx2.The axial stress in the rock bolt can be expressed with the displacement of the rock bolt. As for the confining medium, this is also applicable. Then, the following equations can be acquired [40]:(4)σbx=Ebdubxdx,(5)σmx=Emdumxdx,where Eb is Young’s modulus of the rock bolt, Pa; σb(x) is the axial stress in the rock bolt at the position of x, Pa; Em is Young’s modulus of the confining medium, Pa; and σm(x) is the axial stress applied on the confining medium, Pa.Substituting equations (4) and (5) into (3) leads to(6)d2δxdx2=1Ebdσbxdx−1Emdσmxdx.According to the mechanical equilibrium, the relationship between the axial stress in the bolt and the axial stress applied on the confining medium can be expressed as [55](7)σm=−σbπrb2Am,where Am is the cross-sectional area of the confining medium, m2 and rb is the rock bolt radius, m.Substituting equation (7) into equation (6) leads to(8)d2δxdx2=1Eb+πrb2EmAmdσbxdx.The relationship between the axial stress in the rock bolt and the shear stress along the bolt/grout interface can be seen in Figure2.Figure 2
Cross section view of a rock bolting system.Then, the following equation can be acquired [42]:(9)τx=rb2dσbxdx,where τ(x) is the shear stress of the bolt/grout interface at the position of x, Pa.Substituting equation (9) into equation (8) leads to(10)d2δxdx2−2rb1Eb+πrb2EmAmτx=0.It is defined that(11)λ2=2rb1Eb+πrb2EmAm.Then, equation (10) can be written as(12)d2δxdx2−λ2τx=0.Equation (12) is the governing equation of this analytical model, and it can be solved with the predefined shear stress equations.


## 2.3. Pull-Out Stages of Rock Bolts

According to previous research, the pull-out stages of the bolt/grout interface can be classified as the elastic stage, elastic-softening stage, elastic-softening-debonding stage, softening-debonding stage, and debonding stage [38, 53, 56]. Therefore, in the following contents, the pull-out process of a rock bolt was analysed with those five stages.


### 2.3.1. Elastic Stage

When the pull-out load is small, the bolt/grout interface undergoes elastic deformation. And, the full bolt/grout interface is elastic, as shown in Figure3.Figure 3
The full bolt/grout interface is elastic in the elastic stage whereDb is the diameter of the rock bolt, m.Then, equation (1a) can be substituted into equation (12), leading to(13)d2δxdx2−τpδpλ2δx=0.It is defined that(14)λ12=τpδpλ2.Then, equation (13) can be written as(15)d2δxdx2−λ12δx=0.Equation (15) can be solved with boundary conditions. First, the axial stress in the bolt should be deduced. According to equation (2), the following equation can be acquired:(16)dδxdx=dubxdx−dumxdx.It can be further written as(17)dδxdx=σbxEb−σmxEm.Substituting equation (7) into equation (17) leads to(18)dδxdx=1Eb+πrb2EmAmσbx.Substituting equation (11) into equation (18), the axial stress in the bolt can be written as(19)σbx=2rbλ2dδxdx.The boundary conditions are(20)σbx=0=0,(21)σbx=L=Fπrb2,where F is the pull-out load of the rock bolt, N.Combining equations (19)–(21) with equation (15) leads to(22)δx=Fλ1δp2πrbτpsinhλ1Lcoshλ1x,(23)σbx=Fπrb2sinhλ1Lsinhλ1x,(24)τx=Fλ12πrbsinhλ1Lcoshλ1x.It is defined that at the loaded end, the shear slippage of the bolt/grout interface is equal to the pull-out displacement of the rock bolt [57]. Substituting x = L into equation (22), the pull-out load of the rock bolt can be expressed as(25)F=2πrbτptanhλ1Lλ1δpub,where ub is the pull-out displacement of the rock bolt, m.At the end of this stage, the shear stress at the position ofx = L increases to the shear strength of the bolt/grout interface. Therefore, substituting τ(x = L) = τp into equation (24) leads to(26)Fsof=2πrbτptanhλ1Lλ1,where Fsof is the pull-out load when the bolt/grout interface starts softening, N.


### 2.3.2. Elastic-Softening Stage

In the elastic-softening stage, the full bolt/grout interface is composed of the elastic section and softening section, as shown in Figure4.Figure 4
The full bolt/grout interface is composed of the elastic section and softening section in the elastic-softening stage.It is defined thatas is the softening length. In the elastic-softening stage, the interface in the domain [0, L − as] is still elastic. One boundary condition is equation (20) and the other boundary condition is(27)τx=L−as=τp.Combining equations (9), (19), (20), and (27) with equation (15) leads to(28)δx=δpcoshλ1L−ascoshλ1x,(29)σbx=2τprbλ1coshλ1L−assinhλ1x,(30)τx=τpcoshλ1L−ascoshλ1x.Substitutingx = L − as into equation (29), the axial stress can be acquired:(31)σbx=L−as=2τprbλ1tanhλ1L−as.In the domain [L − as, L], the bolt/grout interface softens. Substituting equation (1b) into equation (12) leads to(32)d2δxdx2+λ2τp−τrδr−δpδx−λ2τpδr−τrδpδr−δp=0.It is defined that(33)λ22=τp−τrδr−δpλ2.Therefore, equation (32) can be further written as(34)d2δxdx2+λ22δx−λ22τpδr−τrδpτp−τr=0.The boundary conditions are equations (27) and (31). Combining equations (9), (19), (27), and (31) with equation (34) leads to(35)δx=λ1δpλ2tanhλ1L−assinλ2x+as−L−λ12δpλ22cosλ2L−as−x+τpδr−τrδpτp−τr,(36)σbx=2τprb1λ1tanhλ1L−ascosλ2L−as−x+1λ2sinλ2x+as−L,(37)τx=τpλ2λ1tanhλ1L−assinλ2L−as−x+cosλ2L−as−x.Substitutingx = L and σbx=L=F/πrb2 into equations (35) and (36), respectively, the pull-out displacement and pull-out load of the rock bolt can be acquired:(38)F=2πrbτp1λ1tanhλ1L−ascosλ2as+1λ2sinλ2as,(39)ub=λ1δpλ2tanhλ1L−assinλ2as−λ12δpλ22cosλ2as+τpδr−τrδpτp−τr.At the end of this stage, the shear stress at the position ofx = L reduces to τr. Therefore, substituting τ(x = L) = τr into equation (37) leads to(40)cosλ2as−λ2λ1tanhλ1L−assinλ2as=τrτp.For an infinite embedment length, tanh(λ1(L − as)) is equal to 1. Then, the softening length when the elastic-softening stage ends can be acquired:(41)as1=1λ2arcsinλ1λ12+λ22−arcsinτrλ1τpλ12+λ22,where as1 is the softening length when the elastic-softening stage ends, m.


### 2.3.3. Elastic-Softening-Debonding Stage

Afterτ(x = L) reduces to τr, the bolt/grout interface enters the elastic-softening-debonding stage. In this stage, the full bolt/grout interface is composed of elastic section, softening section, and debonding section, as shown in Figure 5.Figure 5
The full bolt/grout interface is composed of the elastic section, softening section, and debonding section in the elastic-softening-debonding stage.Assuming that the debonding length isad, the bolt/grout interface in the domain [0, L − ad − as] is still elastic. Therefore, equations (28)–(30) are still valid if L is replaced with L − ad:(42)δx=δ1coshλ1L−ad−ascoshλ1x,(43)σbx=2τprbλ1coshλ1L−ad−assinhλ1x,(44)τx=τpcoshλ1L−ad−ascoshλ1x.As for the domain [L − ad − as, L − ad], the bolt/grout interface softens. Therefore, equations (35)–(37) are still valid if L is replaced with (L − ad):(45)δx=λ1δpλ2tanhλ1L−ad−assinλ2x+ad+as−L−λ12δpλ22cosλ2L−ad−as−x+τpδr−τrδpτp−τr,(46)σbx=2τprb1λ1tanhλ1L−ad−ascosλ2L−ad−as−x+1λ2sinλ2x+ad+as−L,(47)τx=τpλ2λ1tanhλ1L−ad−assinλ2L−ad−as−x+cosλ2L−ad−as−x.Substitutingτ(L − ad) = τr into equation (47) leads to(48)cosλ2as−λ2λ1tanhλ1L−ad−assinλ2as−τrτp=0.Equation (48) depicts the relationship between ad and as. It shows that in this stage, the softening length keeps varying with the debonding length increasing. Combining equations (45) and (48) together, the shear slippage of the bolt/grout interface at the position of x = L − ad can be acquired:(49)δx=L−ad=δr.Substitutingx = L − ad into equation (46), the axial stress can be acquired:(50)σbx=L−ad=2τprb1λ1tanhλ1L−ad−ascosλ2as+1λ2sinλ2as.For the domain [L − ad, L], the bolt/grout interface debonds and the shear stress equals τr:(51)τx=τr.Also, equation (1c) can be substituted to equation (12), leading to(52)d2δxdx2=λ2τr.Substituting equations (19), (49), and (50) into equation (52) leads to(53)δx=δr+12λ2τrx2−L−ad2+τpλ2λ1tanhλ1L−ad−ascosλ2as+τpλ2λ2sinλ2as−λ2τrL−adx+ad−L,(54)σbx=2τprb1λ1tanhλ1L−ad−ascosλ2as+1λ2sinλ2as+2τrrbx+ad−L.Substitutingx = L and σbx=L=F/πrb2 into equations (53) and (54), the pull-out displacement and load of the rock bolt can be acquired:(55)ub=δr+λ2ad2τr2+τpλ2λ1tanhλ1L−ad−ascosλ2as+τpλ2λ2sinλ2asad,(56)F=2πrbτpλ1tanhλ1L−ad−ascosλ2as+τpλ2sinλ2as+τrad.The maximum pull-out load occurs in this stage and can be acquired by solving equations (48) and (56) together.At the end of this stage, the elastic length disappears. Therefore,L − ad = as. Substituting it into equation (48) leads to(57)as2=1λ2arccosτrτp,where as2 is the softening length when the elastic-softening-debonding stage ends, m.


### 2.3.4. Softening-Debonding Stage

After the elastic length disappears, the bolt/grout interface enters the softening-debonding stage. In the domain [0,L − ad], the bolt/grout interface softens. Substituting equations (9), (19), and (20) and τ(x = L − ad) = τr into equation (34) leads to(58)δx=τrδp−δrτp−τrcosλ2L−adcosλ2x+τpδr−τrδpτp−τr,(59)σbx=2τrrbλ2cosλ2L−adsinλ2x,(60)τx=τrcosλ2L−adcosλ2x.Substitutingx = L − ad into equations (58) and (59) leads to(61)δx=L−ad=δr,(62)σbx=L−ad=2τrrbλ2tanλ2L−ad.For the domain [L − ad, L], the bolt/grout interface debonds. Substituting equations (19), (61), and (62) into equation (52) leads to(63)δx=δr+λ2τr2x2−L−ad2+λ2τrλ2tanλ2L−ad−λ2τrL−adx+ad−L,(64)σbx=2τrrbtanλ2L−adλ2+x+ad−L.Substitutingx = L and σbx=L=F/πrb2 into equations (63) and (64), the pull-out displacement and load of the rock bolt can be acquired:(65)ub=δr+λ2τradad2+tanλ2L−adλ2,(66)F=2πrbτrtanλ2L−adλ2+ad.At the end of this stage, the debonding length equals the full embedment length. Therefore, substitutingad = L into equations (65) and (66), the pull-out displacement and load of the rock bolt at the end of the softening-debonding stage can be acquired:(67)ub0=δr+λ2τrL22,(68)F0=2πrbτrL,where ub0 is the pull-out displacement of the rock bolt at the end of the softening-debonding stage, m and F0 is the pull-out load of the rock bolt at the end of the softening-debonding stage, N.


### 2.3.5. Debonding Stage

After the softening length disappears, the bolt/grout interface enters the debonding stage. Since the full bolt/grout interface undergoes debonding, the relationship between the pull-out load and displacement can be expressed as(69)F0−F=2πrbτrub−ub0.Substituting equations (67) and (68) into equation (69) leads to(70)F=2πrbτrL+δr+λ2τrL22−ub.


## 2.3.1. Elastic Stage

When the pull-out load is small, the bolt/grout interface undergoes elastic deformation. And, the full bolt/grout interface is elastic, as shown in Figure3.Figure 3
The full bolt/grout interface is elastic in the elastic stage whereDb is the diameter of the rock bolt, m.Then, equation (1a) can be substituted into equation (12), leading to(13)d2δxdx2−τpδpλ2δx=0.It is defined that(14)λ12=τpδpλ2.Then, equation (13) can be written as(15)d2δxdx2−λ12δx=0.Equation (15) can be solved with boundary conditions. First, the axial stress in the bolt should be deduced. According to equation (2), the following equation can be acquired:(16)dδxdx=dubxdx−dumxdx.It can be further written as(17)dδxdx=σbxEb−σmxEm.Substituting equation (7) into equation (17) leads to(18)dδxdx=1Eb+πrb2EmAmσbx.Substituting equation (11) into equation (18), the axial stress in the bolt can be written as(19)σbx=2rbλ2dδxdx.The boundary conditions are(20)σbx=0=0,(21)σbx=L=Fπrb2,where F is the pull-out load of the rock bolt, N.Combining equations (19)–(21) with equation (15) leads to(22)δx=Fλ1δp2πrbτpsinhλ1Lcoshλ1x,(23)σbx=Fπrb2sinhλ1Lsinhλ1x,(24)τx=Fλ12πrbsinhλ1Lcoshλ1x.It is defined that at the loaded end, the shear slippage of the bolt/grout interface is equal to the pull-out displacement of the rock bolt [57]. Substituting x = L into equation (22), the pull-out load of the rock bolt can be expressed as(25)F=2πrbτptanhλ1Lλ1δpub,where ub is the pull-out displacement of the rock bolt, m.At the end of this stage, the shear stress at the position ofx = L increases to the shear strength of the bolt/grout interface. Therefore, substituting τ(x = L) = τp into equation (24) leads to(26)Fsof=2πrbτptanhλ1Lλ1,where Fsof is the pull-out load when the bolt/grout interface starts softening, N.


## 2.3.2. Elastic-Softening Stage

In the elastic-softening stage, the full bolt/grout interface is composed of the elastic section and softening section, as shown in Figure4.Figure 4
The full bolt/grout interface is composed of the elastic section and softening section in the elastic-softening stage.It is defined thatas is the softening length. In the elastic-softening stage, the interface in the domain [0, L − as] is still elastic. One boundary condition is equation (20) and the other boundary condition is(27)τx=L−as=τp.Combining equations (9), (19), (20), and (27) with equation (15) leads to(28)δx=δpcoshλ1L−ascoshλ1x,(29)σbx=2τprbλ1coshλ1L−assinhλ1x,(30)τx=τpcoshλ1L−ascoshλ1x.Substitutingx = L − as into equation (29), the axial stress can be acquired:(31)σbx=L−as=2τprbλ1tanhλ1L−as.In the domain [L − as, L], the bolt/grout interface softens. Substituting equation (1b) into equation (12) leads to(32)d2δxdx2+λ2τp−τrδr−δpδx−λ2τpδr−τrδpδr−δp=0.It is defined that(33)λ22=τp−τrδr−δpλ2.Therefore, equation (32) can be further written as(34)d2δxdx2+λ22δx−λ22τpδr−τrδpτp−τr=0.The boundary conditions are equations (27) and (31). Combining equations (9), (19), (27), and (31) with equation (34) leads to(35)δx=λ1δpλ2tanhλ1L−assinλ2x+as−L−λ12δpλ22cosλ2L−as−x+τpδr−τrδpτp−τr,(36)σbx=2τprb1λ1tanhλ1L−ascosλ2L−as−x+1λ2sinλ2x+as−L,(37)τx=τpλ2λ1tanhλ1L−assinλ2L−as−x+cosλ2L−as−x.Substitutingx = L and σbx=L=F/πrb2 into equations (35) and (36), respectively, the pull-out displacement and pull-out load of the rock bolt can be acquired:(38)F=2πrbτp1λ1tanhλ1L−ascosλ2as+1λ2sinλ2as,(39)ub=λ1δpλ2tanhλ1L−assinλ2as−λ12δpλ22cosλ2as+τpδr−τrδpτp−τr.At the end of this stage, the shear stress at the position ofx = L reduces to τr. Therefore, substituting τ(x = L) = τr into equation (37) leads to(40)cosλ2as−λ2λ1tanhλ1L−assinλ2as=τrτp.For an infinite embedment length, tanh(λ1(L − as)) is equal to 1. Then, the softening length when the elastic-softening stage ends can be acquired:(41)as1=1λ2arcsinλ1λ12+λ22−arcsinτrλ1τpλ12+λ22,where as1 is the softening length when the elastic-softening stage ends, m.


## 2.3.3. Elastic-Softening-Debonding Stage

Afterτ(x = L) reduces to τr, the bolt/grout interface enters the elastic-softening-debonding stage. In this stage, the full bolt/grout interface is composed of elastic section, softening section, and debonding section, as shown in Figure 5.Figure 5
The full bolt/grout interface is composed of the elastic section, softening section, and debonding section in the elastic-softening-debonding stage.Assuming that the debonding length isad, the bolt/grout interface in the domain [0, L − ad − as] is still elastic. Therefore, equations (28)–(30) are still valid if L is replaced with L − ad:(42)δx=δ1coshλ1L−ad−ascoshλ1x,(43)σbx=2τprbλ1coshλ1L−ad−assinhλ1x,(44)τx=τpcoshλ1L−ad−ascoshλ1x.As for the domain [L − ad − as, L − ad], the bolt/grout interface softens. Therefore, equations (35)–(37) are still valid if L is replaced with (L − ad):(45)δx=λ1δpλ2tanhλ1L−ad−assinλ2x+ad+as−L−λ12δpλ22cosλ2L−ad−as−x+τpδr−τrδpτp−τr,(46)σbx=2τprb1λ1tanhλ1L−ad−ascosλ2L−ad−as−x+1λ2sinλ2x+ad+as−L,(47)τx=τpλ2λ1tanhλ1L−ad−assinλ2L−ad−as−x+cosλ2L−ad−as−x.Substitutingτ(L − ad) = τr into equation (47) leads to(48)cosλ2as−λ2λ1tanhλ1L−ad−assinλ2as−τrτp=0.Equation (48) depicts the relationship between ad and as. It shows that in this stage, the softening length keeps varying with the debonding length increasing. Combining equations (45) and (48) together, the shear slippage of the bolt/grout interface at the position of x = L − ad can be acquired:(49)δx=L−ad=δr.Substitutingx = L − ad into equation (46), the axial stress can be acquired:(50)σbx=L−ad=2τprb1λ1tanhλ1L−ad−ascosλ2as+1λ2sinλ2as.For the domain [L − ad, L], the bolt/grout interface debonds and the shear stress equals τr:(51)τx=τr.Also, equation (1c) can be substituted to equation (12), leading to(52)d2δxdx2=λ2τr.Substituting equations (19), (49), and (50) into equation (52) leads to(53)δx=δr+12λ2τrx2−L−ad2+τpλ2λ1tanhλ1L−ad−ascosλ2as+τpλ2λ2sinλ2as−λ2τrL−adx+ad−L,(54)σbx=2τprb1λ1tanhλ1L−ad−ascosλ2as+1λ2sinλ2as+2τrrbx+ad−L.Substitutingx = L and σbx=L=F/πrb2 into equations (53) and (54), the pull-out displacement and load of the rock bolt can be acquired:(55)ub=δr+λ2ad2τr2+τpλ2λ1tanhλ1L−ad−ascosλ2as+τpλ2λ2sinλ2asad,(56)F=2πrbτpλ1tanhλ1L−ad−ascosλ2as+τpλ2sinλ2as+τrad.The maximum pull-out load occurs in this stage and can be acquired by solving equations (48) and (56) together.At the end of this stage, the elastic length disappears. Therefore,L − ad = as. Substituting it into equation (48) leads to(57)as2=1λ2arccosτrτp,where as2 is the softening length when the elastic-softening-debonding stage ends, m.


## 2.3.4. Softening-Debonding Stage

After the elastic length disappears, the bolt/grout interface enters the softening-debonding stage. In the domain [0,L − ad], the bolt/grout interface softens. Substituting equations (9), (19), and (20) and τ(x = L − ad) = τr into equation (34) leads to(58)δx=τrδp−δrτp−τrcosλ2L−adcosλ2x+τpδr−τrδpτp−τr,(59)σbx=2τrrbλ2cosλ2L−adsinλ2x,(60)τx=τrcosλ2L−adcosλ2x.Substitutingx = L − ad into equations (58) and (59) leads to(61)δx=L−ad=δr,(62)σbx=L−ad=2τrrbλ2tanλ2L−ad.For the domain [L − ad, L], the bolt/grout interface debonds. Substituting equations (19), (61), and (62) into equation (52) leads to(63)δx=δr+λ2τr2x2−L−ad2+λ2τrλ2tanλ2L−ad−λ2τrL−adx+ad−L,(64)σbx=2τrrbtanλ2L−adλ2+x+ad−L.Substitutingx = L and σbx=L=F/πrb2 into equations (63) and (64), the pull-out displacement and load of the rock bolt can be acquired:(65)ub=δr+λ2τradad2+tanλ2L−adλ2,(66)F=2πrbτrtanλ2L−adλ2+ad.At the end of this stage, the debonding length equals the full embedment length. Therefore, substitutingad = L into equations (65) and (66), the pull-out displacement and load of the rock bolt at the end of the softening-debonding stage can be acquired:(67)ub0=δr+λ2τrL22,(68)F0=2πrbτrL,where ub0 is the pull-out displacement of the rock bolt at the end of the softening-debonding stage, m and F0 is the pull-out load of the rock bolt at the end of the softening-debonding stage, N.


## 2.3.5. Debonding Stage

After the softening length disappears, the bolt/grout interface enters the debonding stage. Since the full bolt/grout interface undergoes debonding, the relationship between the pull-out load and displacement can be expressed as(69)F0−F=2πrbτrub−ub0.Substituting equations (67) and (68) into equation (69) leads to(70)F=2πrbτrL+δr+λ2τrL22−ub.


## 3. Calibration of the Input Parameters

In this analytical model, the mechanical properties of the rock bolt and the confining medium, such asrb, Eb, and Em, can be acquired based on the rock bolt pull-out scenario. Specifically, rb and Eb can be acquired from the rock bolt specification provided by the bolt manufacturer. Em can be acquired by conducting uniaxial compressive strength tests on samples cored from the confining medium.As for the mechanical properties of the bond-slip model, they can be calibrated from the experimental load-displacement results. Figure6 shows a typical load-displacement curve of rock bolts.Figure 6
Selecting critical points to calibrate the input bond-slip parameters.In the initial stage of the load-displacement curve, the pull-out load of the rock bolt generally had a linear relationship with the pull-out displacement. Therefore, Point 1 represents the point where the linear relationship of the pull-out load and displacement ends in the load-displacement curve. Then,δp equals the pull-out displacement at Point 1. Also, Fsof equals the pull-out load at Point 1. Substituting Fsof and δp into equation (26), τp can be acquired.When the rock bolt reaches the maximum pull-out load, that point can be selected as Point 2, as shown in Figure6. Then, δr can be selected in the range between the pull-out displacement at Point 1 and the pull-out displacement at Point 2. Then, τr can be calibrated until the maximum analytical pull-out load fits well with the maximum experimental pull-out load.


## 4. Validation of the Analytical Model

In this section, two experimental pull-out tests were used to validate the credibility of this analytical model.


### 4.1. Validation with a Pull-Out Test

Bai et al. [58] conducted experimental in situ pull-out tests on rock bolts. Specifically, rock bolts with a radius of 14 mm were tested, and the embedment length was 3 m. The surrounding rock mass had Young’s modulus of 6 GPa.The mechanical properties of the rock bolt and the confining medium used in this analytical model are tabulated in Table1.Table 1
Mechanical properties of the bolt and the confining medium used in this model.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
14
51
3
6
1As for the input parameters for the bond-slip model, they can be calibrated from experimental pull-out results, following the methods indicated in Section3. Then, the parameters tabulated in Table 2 are used in this paper to analytically model the load transfer performance of rock bolts.Table 2
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
2.2
3.57
1
8.91A comparison between the experimental result and the analytical result is shown in Figure7. There is a close match between the experimental result and the analytical result, confirming the credibility of this analytical model.Figure 7
Comparison between the experimental result acquired from Bai et al. [58] and the analytical result.


### 4.2. Validation with the Other Pull-Out Test

Stillborg [59] conducted a number of pull-out tests in the mine site to study the load transfer performance of cable bolts. Cable bolts with a long embedment length of 10 m were tested. The cable bolts had a radius of 19 mm and Young’s modulus of 83 GPa. The surrounding rock mass had Young’s modulus of 57 GPa. The mechanical properties of the cable bolt and the confining medium used in this analytical model are tabulated in Table 3.Table 3
Mechanical properties of the cable bolt and the confining medium used in this model.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
19
83
10
57
1As for the input parameters for the bond-slip model, they can be calibrated from the experimental pull-out results. In this paper, the parameters tabulated in Table4 are used to calculate the load-displacement relationship of the cable bolt.Table 4
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
1.34
10.37
0.47
35.02The comparison between the experimental result and the analytical result is shown in Figure8. It shows that the analytical result agrees well with the experimental result, further confirming that this analytical model can be used to simulate the load transfer performance of rock bolts.Figure 8
Comparison between the experimental result acquired from Stillborg [59] and the analytical result.


## 4.1. Validation with a Pull-Out Test

Bai et al. [58] conducted experimental in situ pull-out tests on rock bolts. Specifically, rock bolts with a radius of 14 mm were tested, and the embedment length was 3 m. The surrounding rock mass had Young’s modulus of 6 GPa.The mechanical properties of the rock bolt and the confining medium used in this analytical model are tabulated in Table1.Table 1
Mechanical properties of the bolt and the confining medium used in this model.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
14
51
3
6
1As for the input parameters for the bond-slip model, they can be calibrated from experimental pull-out results, following the methods indicated in Section3. Then, the parameters tabulated in Table 2 are used in this paper to analytically model the load transfer performance of rock bolts.Table 2
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
2.2
3.57
1
8.91A comparison between the experimental result and the analytical result is shown in Figure7. There is a close match between the experimental result and the analytical result, confirming the credibility of this analytical model.Figure 7
Comparison between the experimental result acquired from Bai et al. [58] and the analytical result.


## 4.2. Validation with the Other Pull-Out Test

Stillborg [59] conducted a number of pull-out tests in the mine site to study the load transfer performance of cable bolts. Cable bolts with a long embedment length of 10 m were tested. The cable bolts had a radius of 19 mm and Young’s modulus of 83 GPa. The surrounding rock mass had Young’s modulus of 57 GPa. The mechanical properties of the cable bolt and the confining medium used in this analytical model are tabulated in Table 3.Table 3
Mechanical properties of the cable bolt and the confining medium used in this model.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
19
83
10
57
1As for the input parameters for the bond-slip model, they can be calibrated from the experimental pull-out results. In this paper, the parameters tabulated in Table4 are used to calculate the load-displacement relationship of the cable bolt.Table 4
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
1.34
10.37
0.47
35.02The comparison between the experimental result and the analytical result is shown in Figure8. It shows that the analytical result agrees well with the experimental result, further confirming that this analytical model can be used to simulate the load transfer performance of rock bolts.Figure 8
Comparison between the experimental result acquired from Stillborg [59] and the analytical result.


## 5. Parametric Study

A parametric study was conducted to evaluate the input parameters on the load transfer performance of rock bolts. Specifically, the influence of Young’s modulus of the confining medium, the shear strength of the bolt/grout interface, and the shear slippage where shear strength of the bolt/grout interface reaches on the load transfer performance of rock bolts were studied.


### 5.1. Young’s Modulus of the Confining Medium

To study the influence of Young’s modulus of the confining medium on the load transfer performance of rock bolts, an analytical rock bolt pull-out scenario was created. A rock bolt with a radius of 10 mm was tested. Young’s modulus of the rock bolt and the embedment length were 196 GPa and 1.5 m. Then, the mechanical properties of the bolt and the confining medium are tabulated in Table5. As for the input parameters for the bond-slip model, they are tabulated in Table 6.Table 5
Mechanical properties of the bolt and the confining medium.
r
b (mm)
E
b (GPa)
L (m)
E
m (MPa)
A
m (m2)
10
196
1.5
10
1
50
90Table 6
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
2
1.5
0.8
3.5Young’s modulus of the confining medium was changed from 10 MPa to 90 MPa, and the load-displacement performance of rock bolts was calculated, as shown in Figure9. It shows that Young’s modulus of the confining medium had an apparent influence on the load transfer performance of rock bolts. The initial stiffness and the peak load of the rock bolting system increased with Young’s modulus of the confining medium.Figure 9
Influence of Young’s modulus of the confining medium on the performance of rock bolts.This finding agrees well with previous research. Hyett et al. [60] conducted laboratory and in situ pull-out tests on plain cable bolts. In the laboratory, PVC, aluminium, and steel pipes were used to confine cable bolts. The results showed that the peak load of cable bolts increased with Young’s modulus of the confining medium. As for the in situ case, cable bolts were pulled from shale, limestone, and granite. The results showed that when the confining medium changed from shale to granite, the peak load of cable bolts increased apparently.Figure9 also shows that with Young’s modulus of the confining medium increasing, the increasing rate of the peak load was not constant. To further study the influence of Young’s modulus of the confining medium on the peak load, a series of analytical pull-out tests was conducted. The relationship between the rock bolt peak load and Young’s modulus of the confining medium was acquired, as shown in Figure 10.Figure 10
Influence of Young’s modulus of the confining medium on the peak load of rock bolts.It shows that for this scenario, before Young’s modulus of the confining medium increased to 280 MPa, the peak load increases rapidly. After that, the peak load is still increasing but with a much lower slope. This indicates that when Young’s modulus of the confining medium was higher than a critical limit, it had marginal effect on the peak load of rock bolts.


### 5.2. Shear Strength of the Bolt/Grout Interface

To study the influence of the shear strength of the bolt/grout interface on the load transfer performance of rock bolts, three analytical pull-out tests were conducted. In this case, the shear strength of the bolt/grout interface was varied from 2 MPa to 6 MPa. As for the other parameters, they were kept constant.The mechanical properties for the rock bolt and the confining medium are tabulated in Table7. As for the input parameters for the bond-slip model, they are tabulated in Table 8.Table 7
Mechanical properties of the bolt and the confining medium.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
10
196
1.5
1
1Table 8
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
2
1.5
0.5
3.5
4
6A plot of the input bond-slip model is shown in Figure11. Based on the input bond-slip model, the calculated load transfer performance of rock bolts is shown in Figure 12.Figure 11
A plot of the input bond-slip model when the shear strength of the bolt/grout interface was different.Figure 12
Load transfer performance of rock bolts when the shear strength of the bolt/grout interface was different.The results show that the shear strength of the bolt/grout interface played a significant role in determining the load transfer performance of rock bolts. Specifically, when the shear strength of the bolt/grout interface was 2 MPa, the rock bolt had a maximum pull-out load of 146 kN when the pull-out displacement reached 3.1 mm.Compared with that, when the shear strength of the bolt/grout interface was 6 MPa, the rock bolt had a maximum pull-out load of 282 kN, increasing by 93.2%. And, the corresponding pull-out displacement reached 4.5 mm. Apparently, the larger the shear strength of the bolt/grout interface, the higher the maximum load transfer capacity of rock bolts.Furthermore, the results show that the shear strength of the bolt/grout interface also had an effect on the initial stiffness of the rock bolting system. With the shear strength of the bolt/grout interface increasing, the initial stiffness of the rock bolting system also increased, as shown in Figure12.


### 5.3. Residual Shear Strength of the Bolt/Grout Interface

To study the influence of the residual shear strength of the bolt/grout interface on the load transfer performance of rock bolts, another rock bolt pull-out scenario was created. In this case, the residual shear strength of the bolt/grout interface was varied from 0.5 MPa to 1.5 MPa. As for the other parameters, they were kept constant.The mechanical properties of the rock bolt and the confining medium are tabulated in Table9. As for the input parameters for the bond-slip model, they are tabulated in Table 10.Table 9
Mechanical properties of the bolt and the confining medium.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
10
196
1.5
1
1Table 10
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
4
1.5
0.5
3.5
1
1.5A plot of the input bond-slip model is shown in Figure13. As for the load transfer performance of rock bolts, it is shown in Figure 14.Figure 13
A plot of the input bond-slip model when the residual shear strength of the bolt/grout interface was different.Figure 14
Load transfer performance of rock bolts when the residual shear strength of the bolt/grout interface was different.The results show that the residual shear strength of the bolt/grout interface had no effect on the initial stiffness of the rock bolting system. For example, in this case, when the residual shear strength of the bolt/grout interface increased from 0.5 MPa to 1.5 MPa, the initial stiffness of the rock bolting system kept constant.Nevertheless, the residual shear strength of the bolt/grout interface had an effect on the maximum load transfer capacity of rock bolts. For example, in this scenario, when the residual shear strength of the bolt/grout interface was 0.5 MPa, the maximum pull-out load of rock bolts was 226 kN and the corresponding pull-out displacement was 3.6 mm. Then, when the residual shear strength of the bolt/grout interface increased to 1.5 MPa, the maximum pull-out load of rock bolts increased to 256 kN, increasing by 13.3%. And, the corresponding pull-out displacement was 4.9 mm. Apparently, increasing the residual shear strength of the bolt/grout interface was beneficial for improving the maximum load transfer capacity of rock bolts.Additionally, compared with the shear strength of the bolt/grout interface, increasing the residual shear strength of the bolt/grout interface had less apparent effect in improving the maximum load transfer capacity of rock bolts. For example, in the aforementioned analysis, when the residual shear strength of the bolt/grout interface increased by 200%, the maximum pull-out load of rock bolts only increased by 13.3%. However, when the shear strength of the bolt/grout interface increased by 200%, the maximum pull-out load of rock bolts largely increased by 93.2%. Therefore, the shear strength of the bolt/grout interface had more significant effect in deciding the maximum load transfer capacity of rock bolts.The contribution of this paper is that it concentrated on fully grouted rock bolts with long embedment length grouted in the surrounding confining medium. Specifically, based on the previous work [38, 47], this paper acquired the specific analytical solution for this scenario. Furthermore, based on the acquired specific analytical solution, it further studied the influence of Young’s modulus of the confining medium on the load transfer performance of rock bolts. It revealed that although the peak load of rock bolts increased with Young’s modulus of the confining medium, there was a critical limit. Once Young’s modulus of the confining medium was beyond that critical limit, it had marginal effect on the peak load of rock bolts.


## 5.1. Young’s Modulus of the Confining Medium

To study the influence of Young’s modulus of the confining medium on the load transfer performance of rock bolts, an analytical rock bolt pull-out scenario was created. A rock bolt with a radius of 10 mm was tested. Young’s modulus of the rock bolt and the embedment length were 196 GPa and 1.5 m. Then, the mechanical properties of the bolt and the confining medium are tabulated in Table5. As for the input parameters for the bond-slip model, they are tabulated in Table 6.Table 5
Mechanical properties of the bolt and the confining medium.
r
b (mm)
E
b (GPa)
L (m)
E
m (MPa)
A
m (m2)
10
196
1.5
10
1
50
90Table 6
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
2
1.5
0.8
3.5Young’s modulus of the confining medium was changed from 10 MPa to 90 MPa, and the load-displacement performance of rock bolts was calculated, as shown in Figure9. It shows that Young’s modulus of the confining medium had an apparent influence on the load transfer performance of rock bolts. The initial stiffness and the peak load of the rock bolting system increased with Young’s modulus of the confining medium.Figure 9
Influence of Young’s modulus of the confining medium on the performance of rock bolts.This finding agrees well with previous research. Hyett et al. [60] conducted laboratory and in situ pull-out tests on plain cable bolts. In the laboratory, PVC, aluminium, and steel pipes were used to confine cable bolts. The results showed that the peak load of cable bolts increased with Young’s modulus of the confining medium. As for the in situ case, cable bolts were pulled from shale, limestone, and granite. The results showed that when the confining medium changed from shale to granite, the peak load of cable bolts increased apparently.Figure9 also shows that with Young’s modulus of the confining medium increasing, the increasing rate of the peak load was not constant. To further study the influence of Young’s modulus of the confining medium on the peak load, a series of analytical pull-out tests was conducted. The relationship between the rock bolt peak load and Young’s modulus of the confining medium was acquired, as shown in Figure 10.Figure 10
Influence of Young’s modulus of the confining medium on the peak load of rock bolts.It shows that for this scenario, before Young’s modulus of the confining medium increased to 280 MPa, the peak load increases rapidly. After that, the peak load is still increasing but with a much lower slope. This indicates that when Young’s modulus of the confining medium was higher than a critical limit, it had marginal effect on the peak load of rock bolts.


## 5.2. Shear Strength of the Bolt/Grout Interface

To study the influence of the shear strength of the bolt/grout interface on the load transfer performance of rock bolts, three analytical pull-out tests were conducted. In this case, the shear strength of the bolt/grout interface was varied from 2 MPa to 6 MPa. As for the other parameters, they were kept constant.The mechanical properties for the rock bolt and the confining medium are tabulated in Table7. As for the input parameters for the bond-slip model, they are tabulated in Table 8.Table 7
Mechanical properties of the bolt and the confining medium.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
10
196
1.5
1
1Table 8
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
2
1.5
0.5
3.5
4
6A plot of the input bond-slip model is shown in Figure11. Based on the input bond-slip model, the calculated load transfer performance of rock bolts is shown in Figure 12.Figure 11
A plot of the input bond-slip model when the shear strength of the bolt/grout interface was different.Figure 12
Load transfer performance of rock bolts when the shear strength of the bolt/grout interface was different.The results show that the shear strength of the bolt/grout interface played a significant role in determining the load transfer performance of rock bolts. Specifically, when the shear strength of the bolt/grout interface was 2 MPa, the rock bolt had a maximum pull-out load of 146 kN when the pull-out displacement reached 3.1 mm.Compared with that, when the shear strength of the bolt/grout interface was 6 MPa, the rock bolt had a maximum pull-out load of 282 kN, increasing by 93.2%. And, the corresponding pull-out displacement reached 4.5 mm. Apparently, the larger the shear strength of the bolt/grout interface, the higher the maximum load transfer capacity of rock bolts.Furthermore, the results show that the shear strength of the bolt/grout interface also had an effect on the initial stiffness of the rock bolting system. With the shear strength of the bolt/grout interface increasing, the initial stiffness of the rock bolting system also increased, as shown in Figure12.


## 5.3. Residual Shear Strength of the Bolt/Grout Interface

To study the influence of the residual shear strength of the bolt/grout interface on the load transfer performance of rock bolts, another rock bolt pull-out scenario was created. In this case, the residual shear strength of the bolt/grout interface was varied from 0.5 MPa to 1.5 MPa. As for the other parameters, they were kept constant.The mechanical properties of the rock bolt and the confining medium are tabulated in Table9. As for the input parameters for the bond-slip model, they are tabulated in Table 10.Table 9
Mechanical properties of the bolt and the confining medium.
r
b (mm)
E
b (GPa)
L (m)
E
m (GPa)
A
m (m2)
10
196
1.5
1
1Table 10
Input parameters for the bond-slip model.
τ
p (MPa)
δ
p (mm)
τ
r (MPa)
δ
r (mm)
4
1.5
0.5
3.5
1
1.5A plot of the input bond-slip model is shown in Figure13. As for the load transfer performance of rock bolts, it is shown in Figure 14.Figure 13
A plot of the input bond-slip model when the residual shear strength of the bolt/grout interface was different.Figure 14
Load transfer performance of rock bolts when the residual shear strength of the bolt/grout interface was different.The results show that the residual shear strength of the bolt/grout interface had no effect on the initial stiffness of the rock bolting system. For example, in this case, when the residual shear strength of the bolt/grout interface increased from 0.5 MPa to 1.5 MPa, the initial stiffness of the rock bolting system kept constant.Nevertheless, the residual shear strength of the bolt/grout interface had an effect on the maximum load transfer capacity of rock bolts. For example, in this scenario, when the residual shear strength of the bolt/grout interface was 0.5 MPa, the maximum pull-out load of rock bolts was 226 kN and the corresponding pull-out displacement was 3.6 mm. Then, when the residual shear strength of the bolt/grout interface increased to 1.5 MPa, the maximum pull-out load of rock bolts increased to 256 kN, increasing by 13.3%. And, the corresponding pull-out displacement was 4.9 mm. Apparently, increasing the residual shear strength of the bolt/grout interface was beneficial for improving the maximum load transfer capacity of rock bolts.Additionally, compared with the shear strength of the bolt/grout interface, increasing the residual shear strength of the bolt/grout interface had less apparent effect in improving the maximum load transfer capacity of rock bolts. For example, in the aforementioned analysis, when the residual shear strength of the bolt/grout interface increased by 200%, the maximum pull-out load of rock bolts only increased by 13.3%. However, when the shear strength of the bolt/grout interface increased by 200%, the maximum pull-out load of rock bolts largely increased by 93.2%. Therefore, the shear strength of the bolt/grout interface had more significant effect in deciding the maximum load transfer capacity of rock bolts.The contribution of this paper is that it concentrated on fully grouted rock bolts with long embedment length grouted in the surrounding confining medium. Specifically, based on the previous work [38, 47], this paper acquired the specific analytical solution for this scenario. Furthermore, based on the acquired specific analytical solution, it further studied the influence of Young’s modulus of the confining medium on the load transfer performance of rock bolts. It revealed that although the peak load of rock bolts increased with Young’s modulus of the confining medium, there was a critical limit. Once Young’s modulus of the confining medium was beyond that critical limit, it had marginal effect on the peak load of rock bolts.


## 6. Conclusions

A trilinear model was adopted to depict the bond-slip behaviour of the bolt/grout interface. It was assumed that the shear stress of the bolt/grout interface increased linearly to the shear strength of the bolt/grout interface with the shear slippage. Then, the shear stress decreased linearly to the residual shear strength of the bolt/grout interface with the shear slippage. After that, the shear stress kept constant.The pull-out stages of rock bolts were analysed. During the pull-out process, the bolt/grout interface underwent elastic, elastic-softening, elastic-softening-debonding, softening-debonding, and debonding stages.Experimental pull-out tests were used to validate the credibility of this analytical model. The results show that there was a good correlation between experimental and analytical results.A parametric study was conducted to evaluate Young’s modulus of the confining medium, the shear strength of the bolt/grout interface, and the residual shear strength of the bolt/grout interface on the load transfer performance of rock bolts. The results show that increasing Young’s modulus of the confining medium was beneficial for improving the load transfer performance of rock bolts. However, once Young’s modulus of the confining medium was beyond a critical limit, it had a marginal effect on the peak load of rock bolts. Additionally, increasing the shear strength of the bolt/grout interface and the residual shear strength of the bolt/grout interface led to rising of the peak load of rock bolts. However, compared with the residual shear strength of the bolt/grout interface, increasing the shear strength of the bolt/grout interface had more apparent effect in improving the load transfer capacity of rock bolts.



---

*Source: 2678905-2019-11-28.xml*