Abstract

In rock engineering of the cold region, there are a lot of rock joints. The shear characteristics of joints play a decisive role in the stability of rock engineering in the cold area. In this paper, based on the numerical simulation method of particle flow, reasonable microscopic parameters are selected for the numerical simulation of the direct shear test of bolted joints. The results show that the shear stiffness and contact modulus are linearly and positively correlated. The greater the contact modulus, the greater the residual stress, the better the synergetic effect between rock and bolt, and the more developed the microcrack. The smaller the contact stiffness ratio, the greater the residual stress. The shear stiffness decreases with the increase in the contact stiffness ratio, and the larger the contact stiffness ratio, the slower the shear stiffness decreases, while the shear strength does not change with the contact stiffness ratio. The contact stiffness ratio has a weak effect on the number of cracks in the model. The shear stiffness increases with the increase in the parallel bond modulus, and the shear strength decreases with the increase in the parallel bond modulus. The binding stiffness is independent of the shear stiffness, and the peak shear stress decreases with the increase in the binding stiffness ratio. The greater the bond stiffness ratio, the greater the number of cracks.

1. Introduction

There are various types of geological defects and discontinuities, such as faults, joints, slices, and unconformities. To a large extent, the movement and deformation of rock mass are mainly subject to these discontinuities, which have a controlling influence on the stability of rock mass engineering [15]. For a long time, in order to facilitate the research, rock mass is often idealized as a uniform continuous medium [611] in the research of exploring the mechanical properties of rock mass, and the theory describing the deformation and failure process of rock mass is the macroscopic elastic-plastic theory. Such a theoretical model has a certain significance for the research on the mechanical properties of rock mass [1217]. The laboratory test is the main method to study jointed rock mass [16, 1821]. Based on mass shear tests, Barton [22] proposed a JRC-JCS model for estimating the shear strength of rock joints. Lee et al. [23] studied and discussed in depth the characteristics of shear deformation, failure mechanism, and attenuation law of strength of specimens subjected to cyclic shearing. Jiang et al. [24] studied the influence of joint roughness on shear mechanical characteristics. Ghazvinian et al. [25] studied the variation rules of peak shear strength and shear strength of different sawtooth joints. Jahanian and Sadaghiani [26] studied joint shear curves and peak shear strength at different undulation angles. With the development of computer technology and the appearance of a large computer, the numerical simulation method is more and more abundantly used to study the mechanical properties of the rock [2731]. Numerical methods for studying the mechanical properties of the rock mass, such as particle flow code (PFC), have been widely used in the design and construction of rock mass engineering [3235]. Park and Song [36] carried out a series of direct shear tests of joints based on PFC3D and studied the effects of geometric characteristics and microscopic properties of joints on joint shear performance.

The relationship between shear mechanical characteristics of jointed rock mass and normal stress and joint surface topography has been studied from a variety of perspectives [3740]. However, it is not enough just to fit the mechanical properties of the rock mass from the macroscopic phenomenon without studying the intrinsic nature of the mechanical properties produced by rock mass. At present, there are few studies on the variation rules of shear mechanical characteristics of bolted jointed rock mass under different rock microscopic parameters [41, 42]. Due to the difficulties in obtaining and researching methods and means of bolted joint samples with different microscopic parameters, this paper adopts the PFC numerical calculation method to study in detail the effects of different microscopic parameters on the macroscopic shear mechanical behavior of jointed rock mass and the development and evolution characteristics of microcracks. The research conclusion has a certain reference value for the laboratory test.

2. Modelling and Parameter Calibration

2.1. Numerical Model

The numerical model is composed of homogeneous particles which are linked by parallel bonding [4346]. In the parallel bond contact model, the bond is equivalent to a number of springs acting side by side, which can transmit both force and torque and which is more in line with the mechanical properties of rock materials, and it is the most widely used model to simulate the mechanical properties of the rock. The final model is established in this paper. The direct shear boundary conditions and particle generation are shown in Figure 1.

Different from other methods of continuum analysis, the simulation of the macroscopic mechanical behavior of the particle discrete element model must be constructed by assigning values to a series of mesoscopic parameters related to material properties [4751]. The selection of model mesoscopic parameters is related to the accuracy of the final simulation results. For the microscopic parameters of the particle flow model of rock materials, it is usually necessary to select the macroscopic mechanical parameters (shear mechanical parameters and elastic modulus) and the stress-strain relationship in the laboratory physical tests [52]. In this paper, the macroscopic shear mechanical behavior of jointed rock mass under different joint microparameters is studied. Relevant physical and mechanical parameters obtained through laboratory tests are shown in Table 1. The microscopic parameters of rock were obtained through repeated debugging, as shown in Table 2. The peak shear strength of the rock and joints obtained by numerical calculation and the laboratory test was further compared by about 150 times, and the corresponding shear strength parameters, adhesive force and friction angle, were obtained through the Mohr–Coulomb criterion fitting [5355], as shown in Table 3, indicating that the values of the two are relatively consistent.

3. Influence of Rock Mesoscopic Parameters on Macroscopic Shear Mechanical Behavior

In order to study the influence of rock mesoscopic parameters on the macroscopic shear mechanical behavior of bolted jointed rock mass, the control variable method was adopted to conduct the direct shear test numerical simulation and result analysis for the bolted jointed rock mass specimens with different rock mesoscopic parameters. Without changing other parameters, only the microscopic parameters of the research object were changed to observe the change of the shear mechanical response.

3.1. Contact Modulus

The contact modulus is the parameter that controls the elastic modulus of the particle contact. The larger the is, the greater the elastic modulus of the particle contact is. The contact modulus can be combined with the contact stiffness ratio . By obtaining the contact area conversion between the two particles, the normal stiffness kn and tangential stiffness ks of the two particles can be obtained, respectively. Only the contact modulus was changed, which was set as 0.1, 0.3, 1, and 3 (GPa), respectively. The shear deformation curve of the bolted jointed rock mass was obtained, as shown in Figure 2. By analyzing the elastic modulus corresponding to different contact modulus, the shear stiffness and shear strength corresponding to different contact modulus are obtained, as shown in Table 4 and Figure 3. Although in general, the contact modulus and the shear strength increase, the contact modulus has a major impact on the stiffness of the model, so only the relationship between the contact modulus and the shear stiffness is emphatically analyzed.

Since the contact modulus mainly controls the elastic modulus between particle contacts, the larger the is, the harder the lithology is. It can be seen from Figure 2 that the slope of the linear segment of the shear stress-shear displacement curve increases with the increase in the contact modulus. From Figure 3, it can be seen that the relationship between shear stiffness and contact modulus is also basically linear. The larger the contact modulus, the greater the rock stiffness and the smaller the shear displacement corresponding to the shear stress peak without changing the other parameters. In addition, the larger the contact modulus , the greater the postpeak curve fluctuation, the greater the postpeak residual stress, and the better the synergetic effect between rock bolt and rock. The smaller the contact modulus, the smaller the residual stress. This is because the bolt has a regulating effect on the comprehensive shear stiffness of the structural plane. When the shear modulus of the rock is small, the bolt can increase the comprehensive shear stiffness of the structural plane. When the shear modulus of the rock is large, the bolt can reduce the comprehensive shear stiffness of the structural plane.

3.2. Contact Stiffness Ratio

The contact stiffness ratio is the ratio of the normal contact stiffness kn and the tangential contact stiffness ks in the parallel bonding linear group, namely, kn/ks, which is a dimensionless quantity. The contact stiffness ratio needs to be assigned only when the contact method command is used, and kn and ks can be directly assigned when the contact property command is directly used, so it is not necessary to specify the contact stiffness ratio . The PFC program uses the contact stiffness ratio to automatically calculate the normal stiffness kn and the tangential stiffness ks of two contact particles in the software through formula (1) and formula (2):where A is the contact surface area and L is the distance between the centers of the two particles.

The contact stiffness ratio means that the contact modulus is used to distribute the normal stiffness kn and the tangential stiffness ks. The larger the contact modulus , the larger the normal stiffness and the smaller the ks. The smaller the is, the smaller the normal stiffness kn is and the larger the tangential stiffness ks is. Through the control variable method, only was changed, and the contact stiffness ratio was set as 0.5, 1.0, 1.5, and 2.2. The direct shear test was conducted on the structural plane model with bolt and without bolt, and the influence of the contact stiffness ratio on the mechanical behavior of the structural plane was observed. The straight shear-simulated shear deformation curve of rock mass with bolted joints is shown in Figure 4. The corresponding shear stiffness and shear strength of each contact stiffness ratio were obtained as shown in Table 5. The relationship between contact stiffness ratio and shear stiffness and shear strength was plotted as shown in Figure 5.

According to the shear deformation curve of bolted joints in Figure 4, the larger the contact stiffness ratio, the smaller the slope of the linear segment before the peak of the curve, that is, the smaller the shear stiffness. The smaller the contact stiffness ratio, the higher the slope of the linear section of the stress deformation curve. This is because the larger the contact stiffness ratio is, the smaller the tangential stiffness ks of the particles will be, and the decreasing tangential stiffness of the particles will reflect the decreasing shear stiffness of the structural plane on the macrolevel. The smaller the contact stiffness ratio, the greater the residual stress and the more obvious the strain hardening characteristics.

It can be seen from Table 5 and Figure 5 that the shear stiffness decreases with the increase in the contact stiffness ratio, and the model shear stiffness is negatively correlated with the contact stiffness ratio. The larger the contact stiffness ratio, the slower the shear stiffness decreases, which shows the nonlinear change. The shear strength does not change with the contact stiffness ratio.

3.3. Parallel Bonding Modulus

Parallel bond modulus is the cementation between two particles in each parallel bonding model; the parallel bond modulus is set to be 0.01 GPa, 0.3 GPa, 3 GPa, and 5 GPa while other mesoscopic parameters are kept the same, then the shear deformation of the bolted joint is shown in Figure 6. By analyzing the elastic modulus corresponding to different parallel bond modulus, the shear stiffness corresponding to different contact modulus is obtained, and the relationship between the parallel contact bond modulus, shear stiffness, and shear strength is obtained, as shown in Figure 7.

According to the shear deformation curve of bolted joints in Figure 6, the smaller the parallel bonding modulus is, the smaller the slope of the linear segment before the curve peak is, and the higher the peak strength is, the sharper the curve peak is. The higher the parallel bonding modulus, the higher the slope of the linear segment before the curve peak and the lower the peak strength. The strain hardening occurred after the peak, and the parallel bond modulus increased. When the shear displacement reaches to the vicinity of 4 mm, the four curves experience a sudden decrease in stress of different degrees, which is caused by the sudden embrittlement of the rock bite. The greater the parallel bonding modulus, the greater the stress reduction degree, and the shear displacement corresponding to the stress reduction decreases with the increase in the bonding modulus. It shows that the bond strength remains unchanged; the greater the parallel bond modulus, the smaller the displacement generated by the contact fracture, and the brittleness of the contact bond increases. According to Figure 7, the relationship between the parallel bonding modulus and shear stiffness and shear strength and the shear stiffness increases with the increase in the parallel bonding modulus, but the increase in the rate of the shear stiffness decreases with the increase in the parallel bonding modulus. The shear strength decreases with the increase in the parallel bond modulus, and the rate of shear strength decreases with the increase in the parallel bond modulus.

3.4. Parallel Bonding Stiffness Ratio

The meaning of parallel bond stiffness ratio and contact stiffness is similar, the parallel bond stiffness ratio is the ratio of the stiffness tothe caking tangential stiffness , namely ,. When the contact method is used, it is needed to assign value for the parallel bond stiffness ratio , PFC software uses formulas (3) and (4) to calculate the bonding stiffness and the tangential stiffness:where L is the distance between the centers of two particles.

By using the method of control variable to control other mesoscopic parameters as the same, only the parallel bond stiffness ratio is changed and is, respectively, set to 0.5, 1.0, 1.5, and 2.2. The peak shear stress corresponding to different bond stiffness ratios is obtained, and the relationship between the bond stiffness ratio and peak shear stress is shown in Table 6. Figure 8 shows that the bond stiffness ratio does not have much effect on the deformation and stress curve, namely, binder modulus has effect on bonding and tangential stiffness, but almost has no effect on the overall shear stiffness of the model, the stress deformation curves are similar at the elastic stage basic coincidence. As can be seen from the relationship between the bond stiffness ratio and peak shear stress in Figure 9, the peak shear stress decreases with the increase in the bond stiffness ratio, and the decreasing rate of peak shear stress decreases with the increase in the bond stiffness ratio, which tends to be flat. If the bond modulus remains the same, the bond stiffness ratio decreases, that is, the bond tangential stiffness increases. It shows that increasing the tangential stiffness of the bond can affect the shear strength and increase the shear strength.

4. Conclusions

(1)Shear stiffness and contact modulus are linearly and positively correlated. The smaller the contact stiffness ratio of bolted joints, the greater the residual stress and the more obvious the strain hardening characteristics. The shear stiffness decreases with the increase in the contact stiffness ratio. The influence of the contact stiffness ratio on the number of cracks in the model is positive, but the correlation is weak.(2)The shear stiffness increases with the increase in the parallel bond modulus, but the increased rate of the shear stiffness decreases with the increase in the parallel bond modulus. The strain hardening characteristic is enhanced with the increase in the parallel bond modulus.(3)The peak shear stress decreases with the increase in the bond stiffness ratio, and the reduction rate decreases with the increase in the bond stiffness ratio. The bond stiffness ratio also has an effect on the number of cracks, but the effect is far less than that of the parallel bond modulus.

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.

Acknowledgments

This paper was funded by the project supported by the Hunan Provincial Natural Science Foundation of China (2019JJ50031 and 2018JJ2500). The authors wish to acknowledge this support.