The cohesion and internal friction angle were characterized as quadratic functions of strain and were assumed to follow the Mohr-Coulomb criterion after the yield of peak strength. These mechanical parameters and their variations in post-peak softening stage can be exactly ascertained through the simultaneous solution based on the data points of stress-strain curves of triaxial compression tests. Taking the influence of the fault into account, the variation of strata pressure and roadway convergence with coal advancement, the temporal and spatial distribution of axial bolt load were numerically simulated by FLAC3D (Fast Lagrangian Analysis of Continua) using the ascertained post-peak mechanical parameters according to the cohesion weakening and friction strengthening model. The change mechanism of axial load of single rock bolt as abutment pressure changes was analyzed, through the comparison analysis with the results of axial bolt load by field measurements at a coal mine face. The research results show that the simulated results such as the period of main roof weighting, temporal and spatial distribution of axial bolt load are in accordance with field measurement results, so the validity of the numerical model is testified. In front of the working face, the front abutment pressure increases first and then decreases, finally tends to be stable. A corresponding correlation exists between the variation of axial bolt load and rock deformation along the bolt body. When encountered by a fault, the maximum abutment pressure, the influential range of mining disturbance and the roadway convergence between roof and floor before the working face are all increased. In the roadways along the gob, axial bolt loads on the side of the working face decrease, while the other side one increases after the collapse of the roof. As superficial surrounding rock mass is damaged, the anchoring force of rock bolts will transfer to inner rock mass for balancing the tensile load of the bolts.

Full-length bonded bolts are widely used in deep mining engineering and an in-depth understanding of their mechanical characteristics under complex and high ground stress conditions is of great significance for deep roadway support systems. Based on a quantitative GSI rating system of surrounding rocks and rock nonlinear dilatancy angle model, a nonlinear dilatancy angle model suitable for jointed rocks was developed. The Hoek–Brown strain-softening model parameters were transformed into equivalent Mohr–Coulomb strength parameters, and a numerical model of the deep roadway was constructed using FLAC ^3�� numerical simulation software as a tool. The force characteristics of fulllength bonded anchors under different constitutive model and dilatancy angle model conditions were analyzed, and the effects of different lengths of anchors on the stability of the surrounding rock were studied. The obtained results revealed a big difference between the axial forces of bolts calculated by strain-softening and ideal elastic-plastic models. It was also found that bolt shear force was less influenced by the strain-softening behaviors of surrounding rocks. Dilatancy angle greatly affected bolt axial force. Therefore, if the dilatancy angle was neglected, great errors would be created in the calculation results of supporting structure designs. The nonlinear dilatancy angle model of jointed rock masses more accurately captured the stress properties of bolts after field monitoring and analysis. The findings of the study can serve as a guide for calculating the stability of surrounding rocks in deep mining engineering.