This study examines the stability of the Huangyukou landslide in Yanqing District, Beijing, under varying rainfall conditions, focusing on the effects of rainfall infiltration and surface runoff on slope stability. Using a combination of field surveys, geophysical methods, drone photogrammetry, and laboratory testing, a high-precision 2D and 3D numerical model was developed. A hydrological-soil-structure coupling model was employed to simulate rainfall-induced infiltration and runoff processes, revealing that increased saturation and pore water pressure significantly reduce shear strength, enhancing the risk of slope failure. Stability analysis, using a reduction factor method, yielded stability coefficients of 1.06 and 1.04 for 20-year and 100-year return period rainfall scenarios, respectively. The results highlight the critical role of rainfall in destabilizing the upper layers of dolomite and shale, with significant deformation observed in the middle and rear slope sections. This research provides a comprehensive framework for assessing landslide risk under extreme rainfall events, offering practical implications for risk mitigation in similar geological contexts.