A novel creep contact model for rock and its implement in discrete element simulation

The creep behavior of rocks significantly impacts projects' structural safety, including mining, tunneling, and hydraulic engineering. Numerous creep models used in numerical calculations incorporate basic elements, such as Burger's model. However, relying solely on series or parallel combinations of these basic components fails to account for the accelerated creep stage, prompting efforts toward developing enhanced creep models. Despite noted progress, a limited number of studies on the full-process simulation of graded creep loading and the consequent generation of cracks. To overcome this problem, an innovative creep model is proposed based on both the Parallel Bond and the Kelvin-Voigt contact models. This proposed creep model simulates the graded loading creep test of rocks under different load conditions and also conducts a parameter sensitivity analysis for the impact of the proportionate discrepancy between the two contact models. The results indicate that this new creep model could successfully simulate the comprehensive creep process, effectively capture crack progression, and represent ultimate failure modes during the creep process. Notably, to ensure the precision of simulations, it is recommended that the proportion of the Kelvin-Voigt contact model within the total contact is at most 25%.