Numerical Investigation of Micromechanical Failure Evolution in Rocky High Slopes Under Multistage Excavation

High rock slopes are extensively distributed in areas of major engineering constructions, such as transportation infrastructure, hydraulic projects, and mining operations. The stability and failure evolution mechanism during their multi-stage excavation process have consistently been a crucial research topic in geotechnical engineering. In this paper, a series of two-dimensional rock slope models, incorporating various combinations of slope height and slope angle, were established utilizing the Discrete Element Method (DEM) software PFC2D. This systematic investigation delves into the meso-mechanical response of the slopes during multi-stage excavation. The Parallel Bond Model (PBM) was employed to simulate the contact and fracture behavior between particles. Parameter calibration was performed to ensure that the simulation results align with the actual mechanical properties of the rock mass. The research primarily focuses on analyzing the evolution of displacement, the failure modes, and the changing characteristics of the force chain structure under different geometric conditions. The results indicate that as both the slope height and slope angle increase, the inter-particle deformation of the slope intensifies significantly, and the shear band progressively extends deeper into the slope mass. The failure mode transitions from shallow localized sliding to deep-seated overall failure. Prior to instability, the force chain system exhibits an evolutionary pattern characterized by “bundling–reconfiguration–fracturing,” serving as a critical indicator for characterizing the micro-scale failure mechanism of the slope body.