An analytical method and numerical simulation were developed to investigate the shear performance of fiberglass rock bolts (20-tonne and 30-tonne) by conducting sixteen double-shearing tests with both clean and infilled shear interfaces. Following the preparation of the required samples, each test set-up was subjected to different ranges of pretension values. The infilled scenario involved 5 mm thick sandy clay infilled shear interfaces. The results of the double shearing tests unveiled that as pretension increased, so did the confining pressures at the shear interfaces for both clean and infilled joints. Also, an analytical model was developed utilising the Fourier transform, energy balance theory, and linear elastic theory. The result was an empirical relationship that could determine the double shear performance of fibreglass rock bolts in close agreement with the experimental data. Coefficients were incorporated to facilitate model calibration and tuning. Eventually, fast Lagrangian analysis of continua (FLAC) three-dimensional (3D) modelling was utilised to conduct numerical simulations of fibreglass rock bolts subjected to double shearing scenarios. The numerical model was calibrated against experimental data and then extended to conduct a sensitivity analysis on fibreglass rock bolts subjected to double shear test setup variations. Scenarios included rock bolt installation angles, shearing rates, and various host rock strengths. The results revealed that increasing the shear speed from the experimental test baseline yielded substantial displacement increases in the post-failure residual performance of the rock bolts. Changing the installation angle resulted in greater peak shear forces and extended residual zones. The least significant impacts were observed when changing the host rock UCS, suggesting neither rock bolt was drastically impacted by weak or strong host rocks.
