Background: Rainfall-induced landslides are a widespread and destructive geological hazard that resist precise prediction. They pose serious threats to human lives and property, ecological stability, and socioeconomic development. Methods: To address the challenges in mitigating rainfall-induced landslides in high-altitude mountainous regions, this study proposes a digital twin framework that couples multiple physical fields and is based on the spherical discrete element method. Results: Two-dimensional simulations identify a trapezoidal stress distribution with inward-increasing stress. The stress increases uniformly from 0 kPa at the surface to 210 kPa in the interior. The crest stress remains constant at 1.8 kPa under gravity, whereas the toe stress rises from 6.5 to 14.8 kPa with the slope gradient. While the stress pattern persists post-failure, specific magnitudes alter significantly. This study pioneers a three-dimensional close-packed spherical discrete element method, achieving enhanced computational efficiency and stability through streamlined contact mechanics. Conclusions: The proposed framework utilizes point-contact mechanics to simplify friction modeling, enhancing computational efficiency and numerical stability. By integrating stress, rainfall, and seepage fields, we establish a coupled hydro-mechanical model that enables real-time digital twin mapping of landslide evolution through dynamic parameter adjustments.