Abstract Based on the stress redistribution analysis of rock mass during the deep underground excavation, the unloading process of pre-flawed rock material was simulated by distinct element method (DEM). The effects of unloading rate and flaw inclination angle on unloading strengths and cracking properties of pre-flawed rock specimens are numerically revealed. The results indicate that the unloading failure strength of pre-flawed specimen exhibits a power-function increase trend with the increase of unloading period. Moreover, combined with the stress state analysis on the flaws, it is found that the unloading failure strength increases with the increase of flaw inclination angle. The cracking distribution of pre-flawed specimens under the unloading condition closely depends on the flaw inclination angle, and three typical types of flaw coalescence are observed. Furthermore, at a faster unloading rate, the pre-flawed specimen experiences a sharper and quicker unloading failure process, resulting in more splitting cracks in the specimens.