To investigate the impact of rainfall on the stability of granite residual soil slopes, indoor model box tests were conducted at three rainfall intensities (30, 60, 90 mm/h) and two rainfall durations (3. 12 h). The variations in wetting front and vertical displacement were monitored. PFC discrete element software was used to simulate direct shear tests of granite residual soil, calibrate the mesoscopic parameters of granite residual soil for varying moisture contents, and develop a discrete element slope model. The analysis concentrated on the displacement and rotation fields, instability indicators, force chains, and fabric anisotropy to reveal the mesoscopic deformation and mechanical mechanisms underlying slope instability in the model box tests. The results show that when the rainfall intensity reaches 60 mm/h or above, the slip and disturbance range of the slope expand significantly, and the slip body exhibits a circular are shape along the slope face. The slip loss rate of the slope initially decreases and then increases with prolonged rainfall; short-term low-intensity rainfall can stabilize the slope, but continuous rainfall significantly increases the slip loss rate. After 9 hours of rainfall, the displacement and rotation angle of soil particles in the slope increase markedly, forming a distinct circular are slip failure surface. Furthermore, after 9 hours of rainfall, the distributions of force chains and contact force anisotropy within the slope change significantly, with force chains on the slip surface breaking and densely concentrating in stable regions.

For a binary structure slope with a soil layer on the top and a rock layer on the bottom, during the rainfall process, surface runoff will cause soil and water loss on the slope surface and damage to the slope environment. When rainwater infiltrates into the slope, the pore water pressure in the soil gradually increases, the shear strength of the soil decreases, and a weak zone is formed at the soil-rock interface, which has a significant impact on the stability of the slope. Therefore, to study the soil and water loss on the slope surface and the stability of the slope under rainfall conditions, we used theoretical analysis, indoor model tests, and numerical simulations to conduct a comprehensive exploration of this issue, and the following conclusions were formed: the pore water pressure in the shallow layer is greater than that in the deep layer, and the pore water pressure at the toe of the slope is greater than that at the top of the slope; as the slope gradient increases, the time when the pore water pressure at the toe of the slope begins to respond gradually speeds up; the slope displacement first occurs at the lower part of the slope, then in the middle, and finally at the upper part; the time when the displacement at each point on the slope surface begins to respond gradually speeds up with the increase in the slope; the damage form at a small slope gradient is mainly flow sliding, and the damage process is continuous; the damage form at a large slope gradient is mainly flow sliding and overall sliding, and the damage process is continuous and sudden; when the binary structure slope fails, the sliding surface includes the internal sliding surface of the soil and the sliding surface at the soil-rock interface, but when the slope gradient is small, the relative sliding at the soil-rock interface is small, and a continuous sliding surface cannot be formed; and when the slope gradients are small (30 degrees and 40 degrees), the displacement decreases continuously from top to bottom, and no overall sliding surface is formed. The larger values of plastic strain mainly occur in the upper and middle parts of the slope, there is no formation of a continuous plastic strain zone, and the damage mode is flow sliding; when the slope gradients are large (50 degrees and 60 degrees), the displacement is the largest in the upper part, and a large displacement also occurs in the lower part, forming a sliding surface that penetrates through the soil-soil and rock-soil layers. The larger values of plastic strain occur in the upper, middle, and lower parts of the slope, a continuous plastic strain zone is formed, and the damage modes are flow sliding and overall sliding; numerical simulations were carried out on a typical actual slope, and consistent results were obtained.