Abstract
Under natural conditions, joint failure usually occurs after long-term loading rather than short-term loading. In order to prevent rock mass failure caused by creep behavior, bolt is widely used as a mature and effective reinforcement method. Therefore, it is necessary to consider the rheological effect of bolted joint. In this paper, the particle flow code (PFC) was used to study the effect of mesoparameters on the rheological shear mechanical behavior of joints without bolts and with bolts. The effects of mesoparameters such as Maxwell elastic coefficient , Maxwell viscosity coefficient , Kelvin elastic coefficient , Kelvin viscosity coefficient , and friction coefficient were analyzed. The results show that and mainly affect the instantaneous shear displacement, but have little effect on the rheological shear displacement. However, , , and have little effect on the instantaneous shear displacement, and the sensitivity is mainly reflected in the rheological shear displacement. and mainly affect the instantaneous shear displacement, but have little effect on the rheological shear displacement. However, , , and have little effect on the instantaneous shear displacement, and the sensitivity is mainly reflected in the rheological shear displacement. Compared with the parameter sensitivity characteristics without bolt, the instantaneous shear displacement, shear stiffness, and steady rheological displacement under anchor condition are less affected by parameters, indicating that anchor can reduce the influence of joint parameters.
1. Introduction
The mechanical properties of joints are important factors affecting the stability of rock mass engineering [1–3]. Under natural conditions, joint failure usually occurs after long-term loading rather than short-term loading. Under the condition of constant external stress, the displacement and stress field of joint are constantly adjusted and reorganized with time, resulting in the creep phenomenon that the strain increases continuously with time, which leads to the failure of jointed rock mass [4–6]. In order to prevent rock mass failure caused by creep behavior, bolt is widely used as a mature and effective reinforcement method [7, 8]. The mechanical properties of bolted joints can be studied by experiment, theoretical analysis, and numerical simulation [9–13]. At present, numerical calculation methods have been widely used in the field of rock mechanics, forming a variety of numerical calculation methods, including finite element method, discrete element method, motion element method, and meshless method [14–16]. Each method has its own advantages and disadvantages, which can solve geotechnical engineering-related problems. In the process of joint shear, due to the fluctuation of rock rupture, the cracks will occur in the rock mass and lead to stress and strain change [17–19]. However, the numerical calculation method based on continuum cannot easily simulate crack development, joint shear slip of continuous deformation of the line. Therefore, it is necessary to consider the material discontinuity numerical calculation method, such as particle flow code (PFC) [20, 21]. In order to study the mesomechanical characteristics of bolted joints, some scholars adopted PFC to simulate the macroscopic mechanical response of bolted joints through the interaction between mesoparticles. Shang et al. [22] established a numerical calculation model of discontinuous joints, studied the characteristics of joint mesoscopic parameters, and simulated the failure mode and strength characteristics of joints under equal normal stress and equal normal stiffness. The previous work mainly studied the shear mechanical properties of joints under conventional conditions, and the rheological effects of joint shear were seldom considered. Therefore, this paper adopts discrete element numerical calculation method to conduct rheological shear simulation tests on both anchored and unbolted joints and adopts control variable method in parameter sensitivity analysis. The effects of rheological mesoscopic parameters such as Maxwell body elastic coefficient , Maxwell body viscosity coefficient , Kelvin body elastic coefficient , Kelvin body viscosity coefficient , and friction coefficient on shear mechanical behavior of joints were investigated.
2. Numerical Model
The numerical calculation model of joint was established using PFC, and the size of the model was . Based on the established intact rock mass geometry model, the corresponding joint geometry is drawn in CAD. The geometry is imported into PFC by the geometry import command, and it is set to DFN fracture to establish the rock mass geometry model with joints. In the model, different particle parameters are set to simulate the bolt, as shown in Figure 1. In this paper, a flat-joint model is used to simulate the upper and lower rock blocks. Due to shear box in the rheological shear test parts in addition to a constant load, loading wall outside wall is consistent with the only wall loading rate of unconstrained free wall.
In the study of shear mechanical behavior, the experimental phenomenon at the moment of peak shear displacement is an important research object, and it is necessary to keep the simulated peak shear displacement consistent with the real test as much as possible [23–25]. Therefore, in this paper, the results of numerical simulation when normal stress is 1 MPa and laboratory test are the same as the standard, and the trial-and-error method is adopted to calibrate the mesoscopic parameters of rock. Firstly, a complete sample numerical model of flat joints is established to simulate the shear stress-shear displacement curve of flat joints. Finally, the mesoscopic parameters matching the actual situation are calibrated. The microscopic parameters of the numerical model are shown in Table 1. According to the bolt parameters in reference [26], a parallel bond model is used to describe the bolt. The friction coefficient is 0.5, which is smaller than that of rock. The tensile strength of bolt particles is 4.0 GPa, the bond force is 0.9 GPa, and the internal friction angle is 20°. Set a 5 cm long bolt with a diameter of 2 mm, as shown in Table 2 for specific parameters.
3. Calculation Results and Analysis
3.1. Relationship between Maxwell Elastic Coefficient and Macroscopic Shear Stiffness
In the rheological contact model, Maxwell volume elastic coefficient controls the instantaneous displacement of intergranular contact. It is of important reference value to analyze the influence of , the microscopic parameter of intergranular contact, on the macroscopic shear mechanical properties of joints. The control variable method was adopted, set to , , , and (unit: Pa), respectively, to analyze the sensitivity of elastic coefficient of Maxwell body. Shear displacement-time curves of bolt-free joints under different conditions are shown in Figure 2. As can be seen from the figure, the instantaneous shear displacement decreases with the increase of , and has a significant effect on the instantaneous shear displacement. Different corresponds to the same rheological shear displacement, and has no effect on the rheological shear displacement. By analyzing the elastic moduli corresponding to different , the macroscopic shear stiffness corresponding to different is obtained, and the relationship between and shear stiffness is shown in Figure 3. The shear stiffness is obtained from the curve of relation between shear stress and shear displacement; the curve slope represents the shear stiffness. It can be seen that there is a nonlinear relationship between and macroscopic shear stiffness, and the shear stiffness changes greatly and converges to a value, indicating that has a significant influence on macroscopic shear stiffness. The nonlinear expression is used to fit it, and the results are shown in Equation (1); is 0.98, indicating that the formula can well predict the macroscopic shear stiffness under different conditions.
In the case of bolted joint, shear displacement-time curves under different conditions were recorded, as shown in Figure 4. It can be seen that, similar to the case without bolt, the instantaneous shear displacement decreases with the increase of , and has a significant effect on the instantaneous shear displacement. Different corresponds to the same rheological shear displacement, and has no effect on the rheological shear displacement.
By analyzing the elastic moduli corresponding to different , the macroscopic shear stiffness corresponding to different is obtained, and the relationship between and shear stiffness is shown in Figure 5. It can be seen that there is a nonlinear relationship between and macroscopic shear stiffness, and the shear stiffness changes greatly and converges to a value, indicating that has a significant influence on macroscopic shear stiffness. Nonlinear expression was used to fit them, and the results are shown in Equation (2). is 0.96, indicating that this formula can well predict the macroscopic shear stiffness under different . Compared with Formula (1) without bolt, the shear stiffness with bolt is less affected by , and the convergence value is also smaller, indicating that the bolt can reduce the influence brought by the change of joint parameters.
3.2. Influence of Maxwell Volume Viscosity Coefficient on the Relationship between Joint Shear Displacement and Time
In Burger’s contact model, Maxwell viscosity coefficient controls the creep displacement of interparticle contact, and the displacement is closely related to time. Analyzing the microscopic parameter of interparticle contact can effectively reflect the macroscopic mechanical characteristics of Burger’s model parameters in the process of joint shear. The control variable method was adopted, set to , , , and , respectively, to analyze the sensitivity of Maxwell body viscosity coefficient . The shear displacement-time curves of bolt-free joints under different are shown in Figure 6. It can be seen that has little influence on instantaneous shear displacement, but significant influence on rheological shear displacement. The larger is, the smaller the rheological shear displacement is, and there is a negative correlation between the two. The rheological shear simulation results of unbolted joints obtained by changing show that the adjustment of has no significant effect on the joint numerical model displacement, crack. The results show that affects the shear mechanical properties under rheological conditions, but the effect is small when no damage occurs. is a factor affecting particle contact, but different from the displacement affected by . The displacement affected by is controlled by time. Therefore, when the time scale is not very large and is not very small, the displacement affected by is not particularly obvious.
In the case of bolted joint, shear displacement-time curves under different conditions were recorded, as shown in Figure 7. It can be seen that has little influence on instantaneous shear displacement, but significant influence on rheological shear displacement. The larger is, the smaller the rheological shear displacement is, and there is a negative correlation between the two. Compared with Figure 6 without bolt, the discretization of shear displacements and time curves of different in Figure 7 with bolt is smaller, indicating that anchor can reduce the influence of joint parameter on joint aging deformation. The rheological shear simulation results of bolted joints obtained by changing show that adjusting has no significant effect on the joint numerical model displacement and crack. The results show that affects the shear mechanical properties under rheological conditions, but the effect is small when no damage occurs. Compared with the condition without bolt, the influence of on displacement is less obvious in the condition with bolt, and the influence of on displacement can be weakened by bolt.
3.3. Influence of Nonlinear Kelvin Elastic Coefficient on Shear Mechanical Behavior of Joints
In the rheological contact model, the Kelvin elastic coefficient () affects the aging displacement of interparticle contact. The analysis of the micro parameter () of interparticle contact can effectively study the reflection of Burger’s model parameters in macroscopic shear mechanical properties of joints. The control variable method was adopted, set to , , , and (unit: Pa), respectively, to analyze the sensitivity characteristics of the elastic coefficient of Kelvin body . The shear displacement-time curves of boltless joints under different conditions are shown in Figure 8. It can be seen that has little effect on the instantaneous shear displacement, but has an obvious effect on the rheological shear displacement. The larger is, the smaller the rheological shear displacement is, and there is a negative correlation between them.
The rheological shear simulation results of unbolted joints obtained by changing show that adjusting has no obvious effect on the numerical model displacement, crack of the joints. The results show that affects the shear mechanical properties under rheological conditions, but the effect is small when no damage occurs. affects particle displacement, but the displacement is controlled by time. Therefore, when the time range is not very large and is not very small, the displacement affected by is not particularly obvious. The relationship between shear displacement and time at different was recorded, as shown in Figure 9. It can be seen from Figure 9 that has little influence on instantaneous shear displacement, but has obvious influence on rheological shear displacement. The larger is, the smaller the rheological shear displacement is, and there is a negative correlation between them. The rheological shear simulation results of bolted joints obtained by changing show that adjusting has no obvious effect on the numerical model displacement and crack of the joints. The results show that affects the shear mechanical properties under rheological conditions, but the effect is small when no damage occurs. Compared with the condition without bolt, the rock mass displacement with bolt is less affected by .
3.4. Relationship between Joint Shear Displacement and Time under Different Kelvin Body Viscosity Coefficients
The Kelvin viscosity coefficient affects the aging displacement of the contact between particles. The analysis of the microscopic parameters of the contact between particles can effectively study the reflection of the parameters of Burger’s model on the macroscopic shear mechanical properties of joints. Control variable method was adopted, set as , , , and , respectively, to analyze the sensitivity characteristics of Kelvin body viscosity coefficient . Shear displacement-time curves under different conditions are shown in Figure 10. It can be seen that has little influence on instantaneous shear displacement, but significant influence on rheological shear displacement. The larger is, the smaller the rheological shear displacement is, and there is a negative correlation between the two.
The rheological shear simulation results of unbolted joints obtained by changing show that adjusting has no significant effect on the joint numerical model displacement, crack. The results show that affects the shear mechanical properties under rheological conditions, but the effect is small when no damage occurs. is a factor affecting particle contact, but different from the displacement affected by , the displacement affected by is controlled by time. Therefore, when the time scale is not very large and is not very small, the displacement affected by is not particularly obvious. Shear displacement-time curves of bolted joints under different conditions are shown in Figure 11. It can be seen that has little influence on instantaneous shear displacement, but significant influence on rheological shear displacement. The larger is, the smaller the rheological shear displacement is, and there is a negative correlation. The rheological shear simulation results of bolted joints obtained by changing show that adjusting has no significant effect on the joint numerical model displacement and crack. The results show that affects the shear mechanical properties under rheological conditions, but the effect is small when no damage occurs. Compared with the condition without bolt, the rock mass displacement with bolt is less affected by .
3.5. Influence of Joint Friction Coefficient on Joint Mechanical Behavior
In the rheological contact model, the friction coefficient affects the contact force between particles, and the analysis of the friction coefficient can effectively study the reflection of the rheological model parameters on the macroscopic shear mechanical properties of joints. The control variable method was adopted, set to 0.1, 0.3, 0.5, and 0.7, respectively, to analyze the sensitivity characteristics of friction coefficient . Shear displacement-time curves under different conditions are shown in Figure 12. As can be seen from the figure, the instantaneous shear displacement decreases with the increase of , and has a significant impact on the instantaneous shear displacement. The rheological shear displacement corresponding to different is roughly the same, and has no effect on the rheological shear displacement. By analyzing the elastic moduli corresponding to different , the macroscopic shear stiffness corresponding to different is obtained, and the relationship between and shear stiffness is shown in Figure 13. It can be seen that there is a nonlinear relationship between and the macroscopic shear stiffness, and the shear stiffness changes greatly, indicating that has a significant impact on the macroscopic shear stiffness. The nonlinear expression is used to fit it, and the results are shown in Equation (3); is 0.99, indicating that the formula can well predict the macroscopic shear stiffness under different conditions.
Respectively, change the interparticle friction coefficient of , by calculation, the joint can be obtained corresponding crack number, displacement, as shown in Figure 14. It shows that, in the process of shearing, joint displacement and number of crack are significantly affected by . The greater the shear displacement is smaller, the less crack. has a negative relationship with shear displacement and crack number. In the rheological model, the shear strength is controlled by the normal stress and friction coefficient . Therefore, the larger is, the greater the static friction force to overcome under the condition of the same displacement caused by particle contact is, and it is less likely to produce sliding and cracks.
(a)
(b)
In the case of bolted joint, shear displacement-time curves under different conditions are shown in Figure 15. As can be seen from the figure, the instantaneous shear displacement decreases with the increase of , and has a significant impact on the instantaneous shear displacement. The rheological shear displacement corresponding to different is roughly the same, and has no effect on the rheological shear displacement. By analyzing the elastic moduli corresponding to different , the macroscopic shear stiffness corresponding to different is obtained, and the relationship between and shear stiffness is shown in Figure 16. It can be seen that there is a nonlinear relationship between and the macroscopic shear stiffness, and the shear stiffness changes greatly, indicating that has a significant impact on the macroscopic shear stiffness. The nonlinear expression is used to fit it, and the results are shown in Equation (4); is 0.99, indicating that the formula can well predict the macroscopic shear stiffness under different conditions. Compared with Formula (3) without bolt, the shear stiffness with bolt is less affected by , and the anchor can reduce the influence of joint parameter on mechanical properties.
For bolted joints, is changed separately, and it is found that the larger is, the smaller the shear displacement is, the fewer the cracks are, and the less obvious the tension chain concentration on the joints is. has a negative relationship with shear displacement and crack number. Shear strength is controlled by normal stress and friction coefficient and is proportional to . Therefore, the larger is, the more static friction force to be overcome under the condition of the same displacement caused by particle contact and the less prone to sliding and cracking of rock mass. However, the rock mass under the condition of anchor is less affected by than that without bolt, and the anchor can reduce the influence brought by the change of rock mass mesoscopic parameters, which is conducive to enhancing the stability of rock mass.
4. Conclusions
(1)For unbolted joint, and mainly affect the instantaneous shear displacement, but have little effect on the rheological shear displacement. However, , , and have little effect on the instantaneous shear displacement, and the sensitivity is mainly reflected in the rheological shear displacement(2)For bolted joint, and mainly affect the instantaneous shear displacement, but have little effect on the rheological shear displacement. However, , , and have little effect on the instantaneous shear displacement, and the sensitivity is mainly reflected in the rheological shear displacement. It should be noted that the instantaneous shear displacement, shear stiffness, and steady rheological displacement with bolt are less affected by parameters than the parameter sensitivity without bolt, indicating that the anchor can reduce the influence of joint parameters
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.