Abstract
Existing analytical models usually fail to match with the actual conditions due to ignoring the nonlinear behavior of the surrounding material reaction force, which changes progressively with the joint shear displacement from elastic stage to yield stage. To tackle this problem, this study proposes a new analytical model to describe the bolt deformation and bolt contribution from elastic stage to plastic stage. The developed model is verified by available experimental direct shear tests of bolted joints and compared with existing models. Then, based on this model, the effects of the joint dilation angle, the bolt installation angle, the friction angle, and the surrounding material strength on bolt contribution are also analyzed and its implication is further discussed. Our results show that the proposed model can precisely describe the evolution of bolt contribution from elastic stage to plastic stage. Compared with surrounding material strength, the augmentation of the joint dilation angle and friction angle is more beneficial to increase the bolt contribution and the optimal installation angle. The work presented is to attempt to provide a reference for the understanding of bolting mechanism of jointed rock mass, the development of bolting theories and the practice of bolting engineering.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Fully grouted bolt is an important technique for fractured rock mass reinforcement, widely used in the slope, tunnel, and dam, and so on. Bolt anchoring is an effective method to improve the mechanical properties of jointed rock mass. When the unstable rock mass slides along the structural plane or has a sliding trend, the rock mass has a strong transverse restraint effect on the bolt body, and the fully grouted bolt has a large shear stiffness, which restricts the sliding deformation or trend of the unstable rock mass along the structural plane. In this process, the compression stress of surrounding grout material, which generates the bolt transverse shear effect (bolt dowel effect), changes gradually and eventually reaches the yield state [1]. Therefore, the transverse shear effect of bolt has a significant contribution to the anchorage strength of jointed rock mass.
Since Bjurstrom conducted experimental studies on the shear resistance of the rock bolting [2], the reinforcement mechanism of bolts has been widely concerned, and many achievements have been achieved, such as the mechanism of the bolt shear resistance [3,4,5], the interaction mechanism between rock and bolt [6, 7], the failure mechanism of bolts [8, 9], and the dowel shear effect of bolts [10, 11]. In addition to the laboratory test, numerical method and theoretical model have also been used to study the mechanical behavior of the bolt in joint. Haile conducted numerical simulations to analyze the effect of the contribution factors (i.e., bolt diameter, bolt inclination angle, and bolt types) on bolt contribution to joint shear strength [12]. Li et al. conducted numerical experiments on the influence of contribution factors (i.e., concrete strength, bolt installation angle and bolt diameter) on the bolt contribution [13]. Li et al. carried out the numerical experiments of the shear behaviour of cable bolt in FLAC3D [14]. Nie et al. simulated the shear behaviour of rock bolts based on beam theory and DDA code [15].
Meanwhile, the theory of beams was applied to the mechanical behavior of the bolted joint [16,17,18]. The interaction between bolt and surrounding material was simulated by a series of Winkler springs. It is supposed that the dowel effect of bolt is activated by a small displacement, which cannot absolutely reflect the large deformation behavior of the bolt and the surrounding material. Moreover, other researchers have proposed the structure mechanical model to simulate the bolt dowel effect and estimate the bolt contribution to the joint strength [14, 19]. The bolt is considered as a statically determinate beam with both ends fixed, and the bolt axial and transverse shear forces are calculated through the structural mechanics theory in this model. However, this model still cannot reflect the progressive deformation mechanism of bolt transverse shear deformation section.
In summary, the previous research has mainly focused on the two aspects: the influence of the factors like the bolt installation angle, the bolt diameter, the bolt pretension, the grout strength, the joint roughness, the joint opening stress etc., on the bolted joint shear resistance behavior; the calculation of joint shear strength contributed by bolt axial force and bolt transverse shear force. Although the existing work is helpful to explore the mechanical action of the bolted joints, the existing research does not provide a wide acceptance constitutive model. Most researchers pay close attention to the calculation of bolt axial and bolt transverse shear forces subjected to joint shear resistance, but they do not pay more attention to the gradual changes of nonlinear interactions between the bolt and the grout material.
As depicted in Fig. 1, a lot of field observations, experimental direct tests and numerical calculations show that the rock bolt will be "s" deformation shape after the bolt reaches its yield strength [20], which is mainly due to the bolt deformation caused by the transverse shear forces and bending moments, causing the reaction of grout material on the compression side of the bolt. However, in the existing analytical models, the reaction force is simply calculated under the condition of the elastic state or yield state of surrounding grout material, while the transition stage from elastic state to yield state is not considered [21]. That is, the reaction force is a sudden change from elastic state to yield state, as shown in Fig. 2. In fact, in the elastic deformation process, the reaction force of grout on bolt progressively augments with the transverse shear displacement, and is non-uniform distribution along the bolt [1]. As the bolt stress reaches its yield limit value, the reaction force density is evenly distributed along the compression side of the bar. In this process, the length of the surrounding grout reaction force applied to the bolt increases continuously. Therefore, the theoretical model is deduced on the basis of the elastic state or limit state of the reaction force of the surrounding grout on the compression side of the bolt, which will reduce the accuracy of the calculation results.
In addition, during the shearing process, the joint will have the dilatancy effect. Thus, it is must be considered that the perpendicular deformation component to the joint depend on the joint dilation angle. Unfortunately, the influence of this shear dilatancy effect on the axial and transverse deformation of bolt has not been reflected in the existing analytical models.
In this study, we have four main objectives: (1) The theoretical relationships between bolt axial force and bolt axial deformation, and bolt transverse shear force and bolt lateral deformation are established from elastic stage to plastic stage considering the progressive yield of surrounding material based on the structure mechanical model; (2) The theoretical model of the bolt contribution to the bolt axial force and the bolt transverse shear force with the shear displacement are proposed in both the elastic and the plastic stages; (3) The developed model is validated by available experimental direct shear tests of bolted joints; (4) According to the proposed model, the effect of the contribution factors (i.e., bolt inclination angle, joint dilation angle, friction angle, and grout strength) on contribution of bolt are systematically analyzed using this proposed analytical model.
Section 2 of this paper introduces the analytical model. Section 3 verifies the proposed analytical model. Section 4 presents the influence analysis of the contribution factors. Section 5 gives the discussion. Section 6 describes the conclusion.
2 The model
2.1 Mechanical behavior of the bolt in elastic domain
The shear displacement of the joint Uo is caused by the shear loads. Meanwhile, the joint dilatancy during shearing will cause the opening displacement of the joint Ud. Thus, the axial force, No and transverse shear force, Qo are developed in the bar. At the same time, the axial deformation u0 and the transverse deformation vo of the bolt are induced by axial force, No and transverse shear force, Qo respectively, which makes the bar similar to "s" deformation shape (see Fig. 3). Based on the geometrical relationship between the displacement of the joint and the deformation of the bolt in Fig. 3, the axial deformation u0 and transverse deformation vo on one half of the bolt are deduced.
The joint dilation angle is denoted as β, and then the joint opening displacement Ud can be expressed as:
The reaction force of surrounding grout on the compression side of the bolt is caused by transverse deformation of the bar. It is assumed that the reaction force of grout on the bolt is non-uniform distribution along the direction of the bolt axis, and the reaction force of grout on bolt increases gradually with the transverse deformation, and that the reaction force density is evenly distributed along the axis of the bolt until the rock bolt reaches the yield limit. Point O is the intersection point of the rock bolt and joint. The maximum bending moment of the bolt is at Point A, and the rotation and deflection of the bolt are both zero at this location. It also supposed that the length from the joint (point O) and the potential plastic hinge of the bolt (point A) progressively increases with the joint displacement. The simplified progressive deformation model of bolt pressure side and the bolt mechanical model are shown in Figs. 4 and 5, respectively. qx is the grout reaction force density on the pressure side of the rock bolt, and qo is the reaction force density at the point O. The reaction force density q(x) can be assumed to a parabolic distribution as follows:
where m is a coefficient varying with the joint shear displacement, m = (Uoe − Uo)/Uoe. When Uo = 0, m is equal to 1. When the bolt reaches the yield limit (Uo = Uoe), m finally approaches zero.
The reaction force Q(x) and the bending moment M(x) can be obtained by integrating Eq. (4).
where L is the length from the joint (point O) (i.e. x = L) and the potential plastic hinge of the bolt (point A) (i.e. x = 0). Thus, the transverse shear force Qo = Q(L) = qoL − 2mL3/3, the maximum bending moment MA = M(0) = QoL − (2qoL2 − mL4)/4.
To consider the influence of the localized curvature in bolt, the term of the curvature influencing zone, Lc, is introduced by Maekawa and Qureshi [23], as shown in Fig. 6.
where σc is the uniaxial compressive strength of grout. E is the bolt elastic modulus. I is the inertia moment of the bolt cross section, equal to πD4/64. D is the bolt diameter.
Thus, the location of the maximum bending moment L, from the joint, can be calculated by Eq. (8) [23]. The length of L is changed with the transverse deformation vo of the bolt in elastic conditions.
Based on the approximate differential equation of beam rotation and deflection, the distribution function of rotation and deflection on one side of bolt can be derived from Eq. (6).
where φ(x) and v(x) are, respectively, the rotation and deflection equation of the bolt on one side of the joint; when x = 0, φ(x) and v (x) are both 0; when x = L, φ(L) and v(L) are, respectively, the bolt rotation angle and the transverse shear displacement at the joint.
The axial displacement of the bolt at the joint u0 is caused by the axial force No. Therefore, the axial displacement of the bolt u0 can be expressed as:
where Ab is the area of the bolt cross section, which is equal to πD2/4.
Based on Eq. (9), the rotation angle of the bolt φo at the joint can be obtained:
According to Eq. (10), and considering the influence of shear deformation on v0, the transverse displacement v0 can be expressed by.
where k is equal to 4/3 [24].
According to Eq. (11) and combining Eq. (12), Eq. (13) and Qo = Q(L) = qoL-2mL3/3, the axial force, No, transverse shear force, Qo and reaction force density qo at the joint in the elastic deformation stage have the following forms:
There are two yielding criteria of the bolt [6]: one is the tensile-shear yield of the rock bolt at joint, which is caused by axial force, No and transverse shear force, Qo; the other is the bending yield of the rock bolt at point A, which is caused by axial force No and bending moment MA.
where Noe is the yield axial force at joint, Qoe is the yield transverse shear force at the point O; σe is the yield stress of the bar.
2.2 Mechanical behavior of the bolt after yielding
After the bolt reaches the yield, the reaction pressure density of the surrounding grout material on the pressure side of the bar quickly reaches the limit value Pu (the maximum pressure of the grout reaction per unit length), as shown in Fig. 7. Beyond the yield of the bolt, the section OA of the bolt has the mechanical properties of the truss, so only the bolt axial force, No, increases gradually, and transverse shear force, Qo, remains unchanged [25]. When the ultimate strain of the bolt εf is reaching, the bolt failure at the joint will be occur. Thus, the bolt failure criteria (point O) at the joint caused by the combination of axial force Nof and transverse shear force Qof can be expressed as:
where Qof is equal to Qoe, Nop is equal to Abσp, Qop is equal to \(A_{b} \sigma_{p} /\sqrt 3\), σp is the failure stress of the bolt.
According to the Eq. (19), Nof is expressed as,
In addition, based on the geometrical consideration in Fig. 7, we can obtain the following equation to correlate the bolt installation angle α, the plastic length of the bar Lpj and the shear displacement of joint Uopj at step j:
where Uopj and ΔUopj are the globe shear displacement and the shear displacement increment of joint in the plastic stage at step j. Leu and Δαeu are the distance between the joint and the plastic hinge (point A) and the increment rotation of the bolt at the yield, respectively. Lpj and Δαpj are, respectively, the distance from the joint to the plastic hinge (point A) and the increment rotation angle of the bolt in plastic stage at step j. εpj is the strain of the bolt at step j.
From Eq. (21), the increment rotation inclination angle of the bolt in plastic stage Δαpj at step j can be expressed as:
Due to the complex of the axial deformation characteristics of the bolt between yield and failure, the description of the development of the bolt axial force is rather difficulty. Based on Neal [26], the diminution of the bolt axial rigidity can be described by a power function. The bolt axial rigidity is related to the plastic strain of the bolt, thus
where Kj is the axial rigidity of the bolt at step j, A and B are the coefficient, and are discussed in the part 5.
Therefore, the axial force of the bolt after yielding can be expressed as,
where Nopj is the axial force at step j.
According to the Eq. (26), when the strain εpj = εf, the axial force of the bolt Nopj is equal to Nof, thus Kf can be expressed as:
2.3 Bolt contribution to joint shear resistance
According to the mechanical behavior of the bolted joint, the bolt contribution to the joint shear resistance T should be divided into two components: the bolt contribution to the joint shear resistance provided by the axial force, TN and the bolt contribution to the joint shear resistance provided by the transverse shear force, TQ. In this paper, the bolt contribution caused by the transverse shear force is regarded as the dowel effect. The values of T, TN and TQ can be expressed as:
where ϕb is the joint friction angle.
The superiority degree of the dowel effect, f can be express as the ratio of TQ to TN.
The f can quantify the advantage of TQ relative to TN. f > 1, the joint shear resistance is almost contributed by the transverse shear force of the bar; f < 1, the joint shear resistance is almost offered by the axial force of the bar; f = 1, the shear-resist effect of the two components is the same.
When calculating the shear strength of the bolt in plastic stage, the bolt inclination angle α in Eqs. (28), (29) and (30) should be modified as follow:
Excel Worksheet program is used to calculate the whole process curve of No, Qo and T against the joint shear displacement Uo respectively. The specific calculation procedure is shown in Fig. 8.
3 Model validation
A series of laboratory direct shear tests of the bolted joints were carried out by Spang [5, 27]. In order to demonstrate the proposed analytical model, some available direct shear tests of the reinforced joint are used. In his study, two blocks with the size of 150 × 150 × 130 mm are stacked up to form a specimen with a steel bolt installed across the joint center. The diameter of the bolt is 8 mm. The bolt inclination angles are varied from 60° to 90° for shear testing. Granite (σc = 105 MPa), concrete (σc = 40 MPa), sandstone (σc = 10 MPa) are used as block material to simulate the rock. The shear load, the bolt deformation, and the joint shear displacement are recorded during the loading. Based on his experiments, we can obtain the material properties of the bolt and joint directly, as shown in the Table 1. Due to the small thickness of mortar material in the laboratory shear tests, only the rock material is considered as the surrounding material of the bolt in the analytical model.
Figures 9, 10, 11 and 12 show the comparison of the complete curves of shear displacement versus shear load of the bolted joint from the experiments and those from other analytical model without considering progressive yield of surrounding material [27]. The results show that the proposed analytical model is able to precisely predict the bolt contribution evolution for the direct shear tests. Because Pellet did not consider progressive yield of surrounding material and the joint dilatancy effect, when the strength of the surrounding material is large or small, the theoretical calculation by Pellet has a large deviation from the experimental results. Therefore, comparison of the analytical model, the proposed analytical results are more consistent with experimental data. In addition, it can also be concluded in Figs. 9, 10, 11 and 12, that the increasing surrounding material strength could increase the bolt contribution, but reduce the displacement at failure.
The calculation results also show that the bolt yielding meets the tensile-shear yield criteria due to the small diameter of bolt and the hard surrounding material, which is consistent with the conclusion obtained by Ferrero [6]. Ferrero believes that the yielding form of bolt is related to the bolt diameter and the surrounding material strength. The small diameter bolt generally yields by the combination of tensile and shear, while a large diameter bolt may bend yield. The tensile-shear yield occurs in hard material, while the bending yield occurs easily in soft material.
4 Analysis of the bolt contribution
In order to analyze the influence of the different factors (i.e., joint dilation angle β, surrounding material strength σc, friction angle ϕb, and bolt installation angle α) affecting the bolt contribution at the yield limit and failure respectively, we carried out a parametric study. The bolt total contribution (T) at the yield limit, Ty, the bolt total contribution (T) at the failure, Tf, and the ratio of TQ to TN, f, at the yield limit were calculated. Thus, we obtained the complete curves of Ty, Tf, TQ/TN versus bolt installation angle α with different β, σc, and ϕb. The parameters required in the analytical model are shown in Table 2.
4.1 Influence of bolt installation angle on the ratio f
The effect of the bolt installation angle variation on the ratio f is shown for different values of joint dilation angle in Fig. 13. As can be seen in Fig. 13, for different joint dilation angle, the installation angle has the same width of fluctuation range, in which the dowel effect is more obvious and the dominant of the dowel effect is greatly affected by installation angle, but there is a delay in the position of fluctuation range. For example, when β = 0°, the fluctuation range of installation angle is 80° < α ≤ 130°. When the joint dilation angle increase to 10°, the fluctuation range of installation angle is delayed for 10°, that is, 90° < α ≤ 140°. The fluctuation range delay value is consistent with the increment of joint dilation angle. When α ≤ 80°, the ratio f is not influenced by the installation angle and joint dilation angle, and f is less than 1. In other words, the bolt contribution is almost generated by the axial force.
Figure 14 shows the influence of the bolt installation angle α on the ratio f for different values of surrounding material strength. As the surrounding material strength increases, the width and location of fluctuation range remains the same, that is, 80° < α ≤ 130°. When 110° < α ≤ 130°, the ratio f increases with the surrounding material strength rising, but under others bolt installation angle, the ratio f is not affected by the surrounding material strength. When α < 80°, and 130° < α ≤ 180°, the ratio f is less than 1, and not influenced by the installation angle and surrounding material strength.
The influence of the bolt installation angle α on the ratio f for different values of friction angle is shown in Fig. 15. The width of fluctuation range becomes lager when the friction angle increases. For example, when ϕb = 10°, the fluctuation range of installation angle is 80° < α ≤ 110°. When ϕb = 32°, the fluctuation range of installation angle is 80° < α ≤ 130°. When ϕb = 50°, the fluctuation range of installation angle is increased to 80° < α ≤ 150°. When α < 80°, and 150° < α ≤ 180°, the ratio f is less than 1, and not influenced by the installation angle and friction angle. In fluctuation range of installation angle, the dowel effect can be improved by increasing the friction angle.
From Figs. 13, 14 and 15, the evolution curve has two regions: weak impact range and fluctuation range. In the weak impact range, the axial force mobilized in the bolt is more obvious, and the ratio is almost not influenced by the joint dilation angle, friction angle, and surrounding material strength. However, in the fluctuation range, the dowel effect becomes obvious, but the width, location, and peak value of the fluctuation range influenced by the joint dilation angle, friction angle, and surrounding material strength. The joint dilation angle will affect the position of the fluctuation range, which leads to the delay of the fluctuation range. The surrounding material strength will affect the peak value of the fluctuation range. The friction angle will affect the peak value and width of the fluctuation range.
4.2 Influence of bolt installation angle on the bolt contribution
The influence of the bolt installation angle α on the bolt contribution at yield limit Ty and failure Tf for different values of joint dilation angle is shown in Fig. 16. At the yield limit, the bolt contribution evolution with the installation angle can be divided into two regions: peak region (0° < α ≤ 120°) and rise region (120° < α ≤ 180°). In peak region, as the joint dilation angle increases, the bolt contribution gets lager. However, in rise region, the joint dilation angle bolt does not influence the bolt contribution. In peak region, the turning point (a maximum bolt contribution point) is not a fixed value but increases with the increase of joint dilation angle. The increment of the optimal angle is almost equal to that of the dilation angle. In rise region, the bolt contribution increases linearly with the augmentation of the bolt installation angle.
At failure, the bolt contribution evolution with the installation angle can be divided into three regions: rise region, parabolic region and descending region. At the same installation angle, with the augmentation of the joint dilation angle, the bolt contribution is larger. In the rise region, the bolt contribution increases with the increase of the installation angle. In the parabolic region, when α = 90°, the bolt contribution reaches a minimum value, and the installation angle corresponding to the lowest point does not change with the joint dilation angle. In the descending region, the bolt contribution decreases with the augmentation of the installation angle. With the increase of joint dilation angle, the slope of the curve has a large change rate in the rise and descending regions, while the curve becomes smoother in the parabolic region.
The influence of the bolt installation angle on the bolt contribution at the yield limit Ty and the failure Tf for different values of surrounding material strength is shown in Fig. 17. At the yield limit, the bolt contribution evolution with the installation angle can also be divided into peak region and rise region. When 90° < α ≤ 140°, the bolt contribution is greatly influenced by the surrounding material strength, and increases with the increase of surrounding material strength. On the other hand, the bolt contribution is not obviously affected by surrounding material strength for other installation angles. The turning point (a maximum bolt contribution point) is a fixed value, not influenced by the surrounding material strength.
At the failure, the surrounding material strength σc does not affect the finial contribution of the bolt. It can be found from Fig. 17 that with the decrease of the surrounding material strength, the ultimate shear resistances of the bolt have the same value. The reason for this is that the compressive strength of surrounding material will affect the reaction force on the bolt per unit length Pu (Eq. 33). Pu will improve the bolt contribution provided by the transverse shear force. Based on Eq. 33, we can conclude that the greater the compressive strength of material is, the greater the reaction force on the bolt per unit length Pu is, that is, the greater the bolt contribution at the yield limit. However, in the plastic stage, when the reaction force reaches the maximum value, the transverse shear force in bolt will not be improved, so the surrounding material strength will not affect the finial bolt contribution.
The influence of the bolt installation angle α on the bolt contribution at the yield limit Ty and the failure Tf for different values of friction angle ϕb is shown in Fig. 18. At the yield limit, the curve of bolt contribution to the installation angle can also be divided into two regions. In the peak region, when α = 90°, the friction angle ϕb has little influence on the bolt contribution. When 0° < α ≤ 90°, the bolt contribution increases with the augmentation of the friction angle ϕb, but decreases with increasing the friction angle ϕb when the installation angle ranges from 90° to 180°. The turning point is also not a fixed value but increases with the augmentation of friction angle ϕb.
At failure, the bolt contribution evolution with the installation angle can also be divided into three regions: rise region, parabolic region and descending region. With the augmentation of friction angle ϕb, the installation angle range of rise region and descending region increases, while the installation angle range of region parabola region decreases. And the bolt contribution increases with the augmentation of the friction angle, whatever of the bolt orientation is.
According to Figs. 16, 17 and 18, the curves of the resistant contribution of the bolt over the change of joint dilation angle, surrounding material strength and friction angle have nearly the same shape characteristics in the cases of yield and failure limit of bar. Compared with the joint dilation angle and friction angle, the surrounding material strength has little influence on the bolt contribution at yield limit Ty and the failure Tf. But with the augmentation of joint dilation angle and friction angle, the contribution of bolt can be easily mobilized.
The change of surrounding material strength does not affect the optimal installation angle, but with the augmentation of friction angle and joint dilation angle, the optimal installation angle increases. The optimum installation angle between the joint and the bolt is approximately equal to the joint friction angle. The curves of failure shear resistance of the bolted joint over the change of joint dilation angle, surrounding material strength and friction angle have an almost symmetric distribution. The symmetry point is located near 90°.
5 Discussion
The proposed analytical model offers an effective methodology for predicting the bolt contribution to joint strength-shear displacement curves. The limitation of the model is that the joint shear strength is not considered. Due to the interaction between bolt and joint in the process of shear loading, the influence of joint shear strength on axial force and transverse shear force should be considered in the model.
Although there are many researchers have analyzed the hinge point location L, but they judge the value L (i.e., 2–3 times bolt diameter) based on experiments. Jalalifar and Aziz derived L according to the bolt yield criterion (Eq. 34) [20]. However, the value of L gradually decreases with the shear displacement, which is inconsistent with the actual situation. In this study, we introduce and modify the term of the curvature influencing zone, Lc and the hinge point location, L. Before yielding, the hinge point location, L, gradually increases with the shear displacement (Eqs. 7 and 8). However, there is an uncertain parameter a in the formula. According to the existing literature, a is taken as 0.8 [23]. After calculation, it is found that the value of a is related to the bolt diameter, installation angle, surrounding material strength and other factors, so it is necessary to further investigate.
In the previous models, the bolt axial deformation has not been considered properly. For instance, the influence of the joint dilation angle and the increment rotation on the bolt axial deformation is not properly considered. In addition, the bolt axial deformation is not considered after yielding in most analytical models. Therefore, we calculate the bolt axial deformation by taking into consideration the joint dilation angle and the increment rotation from elastic stage to plastic stage. However, the bolt–grout interface bond stress is ignored in our model, thus the bolt axial force in the bolt axis does not change.
In the plastic hardening stage, the diminution of the bolt axial rigidity is described by a power function, which is verified by the experimental tests results. However, the bolt axial rigidity (Eq. 25) is related to two uncertain parameter A and B. Since we can calculate Kf, if we know K1 (the axial rigidity of the bolt at step 1), it is easy to obtain A and B. However, there is very little discussion about K1. In our study, we use the programming solver to calculate A and B by setting constraint and objective function:
The following constraints should be satisfied:
where n is the number of intervals between elastic limit and failure, Tbeu and Tbp1 are, respectively, the bolt contribution at elastic limit and the bolt contribution in plastic stage at step 1.
6 Conclusion
In this paper, an analytical model of bolted rock joint is proposed to predict the bolt contribution evolution, the bolt axial force evolution and the bolt transverse shear force evolution. In our model, we take into account the progressive yield characteristic of surrounding material. In addition, the joint dilation angle and the increment rotation are considered from the elastic stage to the plastic stage. The proposed model agrees well with the available experimental results. The results show that the proposed model can predict accurately the evolution of the bolt contribution, the bolt axial force, and the transverse shear force from elastic stage to plastic stage. Moreover, parametric studies were conducted about four main contribution factors including joint dilation angle, bolt installation angle, friction angle and surrounding material strength, based on the proposed analytical model. Our results show that the above improvements of the model are more consistent with the actual conditions of the bolt deformation.
References
Holmberg M (1991) The mechanical behaviour of untensioned grouted rock bolts. Ph.D. thesis, Royal Institute of Technology, Stockholm
Bjurstrom S (1974) Shear strength of hard rock joints reinforced by grouted untensioned bolts. In: Proceedings of the 3rd international conference ISRM congress, Denver, USA
Haas CJ (1981) Analysis of rock bolting to prevent shear movement in fractured ground. Min Eng 33(6):698–704
Ludvig B (1983) Shear tests on rock bolts. In: rock bolting: theory and application in mining and underground construction, proceedings of the international symposium on rock bolting. Abisko, pp 193–203
Spang K, Egger P (1990) Action of fully-grouted bolts in jointed rock and factors of influence. Rock Mech Rock Eng 23(3):201–229. https://doi.org/10.1007/BF01022954
Jalalifar H, Aziz N (2010) Analytical behaviour of bolt-joint intersection under lateral loading conditions. Rock Mech Rock Eng 43:89–94. https://doi.org/10.1007/s00603-009-0032-6
Grasselli G (2005) 3D behaviour of bolted rock joints: experimental and numerical study. Int J Rock Mech Min Sci 42(1):13–24. https://doi.org/10.1016/j.ijrmms.2004.06.003
Li YZ, Liu CH (2019) Experimental study on the shear behavior of fully grouted bolts. Constr Build Mater 223:1123–1134. https://doi.org/10.1016/j.conbuildmat.2019.06.207
Jalalifar H, Aziz N, Hadi M (2006) The effect of surface profile, rock strength and pretension load on bending behaviour of fully grouted bolts. Geotech Geol Eng 24:1203–1227. https://doi.org/10.1007/s10706-005-1340-6
Ferrero AM (1995) The shear strength of reinforced rock joints. Int J Rock Mech Min Sci Geomech Abstr 32(6):595–605. https://doi.org/10.1016/0148-9062(95)00002-X
Ferrero AM, Giani GP, Migliazza M, Egger P (1997) Physical and numerical analysis of reinforced rock joint behaviour. Rivista Italiana di Geotecnica 3:63–89
Haile A (1999) A mechanistic evaluation and design of tunnel support systems for deep level South African mines. PhD thesis, University of Natal, Durban, South Africa
Li L, Hagan PC, Saydam S, Hebblewhite B (2016) Li Y (2016) Parametric study of rockbolt shear behaviour by double shear test rock. Rock Mech Rock Eng 49(12):4787–4797. https://doi.org/10.1007/s00603-016-1063-4
Li XW, Nemcik J, Mirzaghorbanali A, Aziz N, Rasekh H (2015) Analytical model of shear behaviour of a fully grouted cable bolt subjected to shearing. Int J Rock Mech Min Sci 80:31–39. https://doi.org/10.1016/j.ijrmms.2015.09.005
Nie W, Guo W, Ma SQ, Zhao ZY (2020) Numerical modelling of fully grouted rockbolts subjected to shear load. Rock Mech Rock Eng 53(5):2493–2503. https://doi.org/10.1007/s00603-020-02049-8
Stimpson B (1987) An analytical method for determining shear stiffness of an inclined grouted bolt installed across on open discontinuity. Int J Min Geo-Eng 5(3):299–305
Sq MA, Zhao ZY, Peng J, Gui Y (2018) Analyticalmodeling of shear behaviors of rockbolts perpendicular to joints. Constr Build Mater 175:286–295
Ma SQ, Zhao ZY, Shang JL (2019) An analytical model for shear behaviour of bolted rock joints. Int J Rock Mech Min Sci 121:104019. https://doi.org/10.1016/j.ijrmms.2019.04.005
Liu CH, Li YZ (2017) Analytical study of the mechanical behavior of fully grouted bolts in bedding rock slopes. Rock Mech Rock Eng 50(2):413–2423. https://doi.org/10.1007/s00603-017-1244-9
Jalalifar H, Aziz N (2010) Experimental and 3D numerical simulation of reinforced shear joints. Rock Mech Rock Eng 43:95–103. https://doi.org/10.1007/s00603-009-0031-7
Zhao ZH, Sun W, Chen SJ, Wang WM, Wang QB (2020) Coupling model of jointed rock mass and rock bolt in offshore LPG underground storage. Energy Sci Eng 8:1468–1483
Li CC (2010) Field observations of rock bolts in high stress rock masses. Rock Mech Rock Eng 43(4):491–496
Maekawa K, Qureshi J (1996) Computational model for reinforcing bar embedded in concrete under combined axial pullout and transverse displacement. Proc JSCE 31(538):227–239. https://doi.org/10.2208/jscej.1996.538_227
Gere JM, Timoshenko SP (1990) Mechanics of materials. Kent Publishing Company, Amsterdam
Pellet F, Egger P (1996) Analytical model for the mechanical behaviour of bolted rock joints subjected to shearing. Rock Mech Rock Eng 29:73–97. https://doi.org/10.1007/BF01079755
Neal BG (1977) The plastic methods of structural analysis. Wiley, New York
Pellet F (1994) Strength and deformability of jointed rock masses reinforced by rock bolts. Ph.D. thesis, Swiss Federal Institute of Technology, Lausanne, Switzerland
Acknowledgements
This study was financially supported by the National Natural Science Foundation of China through Grants 51809194 and Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin (China Institute of Water Resources and Hydropower Research), Grant No: IWHR-SKL-201817.
Author information
Authors and Affiliations
Contributions
W.C. designed this study, proposed the method, analyzed the results, and wrote this manuscript; Y.L. revised this manuscript. All authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wen-qiang, C., Yi-jia, L. Analytical model of bolt shear resistance considering progressive yield of surrounding material. SN Appl. Sci. 4, 42 (2022). https://doi.org/10.1007/s42452-021-04923-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s42452-021-04923-8