Vegetation roots enhance soil stability by influencing saturation and pore structure, playing a pivotal role in stabilizing slopes, reducing erosion, and enhancing soil structure. However, current research on the hydraulic effects of roots on soil remains relatively limited. The micro-mechanisms of vegetation's impact on soil and the macro-level connections are not yet fully understood, which poses a challenge to the modeling of root-soil system. This study develops a three-dimensional (3D) finite element model of root-soil composites based on root computed tomography (CT) images and experimental results. Four different groups are modeled, including the rootless group, and those with Festuca arundinacea (FA) roots at various growth stages. The simulation results show that the saturation in the shallow layers significantly decreases in root-soil composite groups, and the rhizosphere water content is lower than that away from the roots, resulting in a net water flux toward the roots. The influence range of roots on suction is gradually amplified with increasing root growth process and root water uptake time. Higher levels of root development result in a stronger overall water uptake effect, leading to a more pronounced decrease in saturation. Closer proximity to the surface roots results in a more rapid increase in soil suction. Compared with one-dimensional root water uptake models, this model considers the effects of spatial heterogeneity of root structures on soil, which provides a comprehensive modeling basis for studying the effect of root system on soil.