A tensor-type capillary stress, instead of a scalar suction, has been proposed to serve as a stress-like state variable to capture the effects of capillarity in the mechanics of unsaturated granular soils. However, the influence of water content on the evolution of capillary stress in such soils remains insufficiently understood. This study performs numerical simulations of unsaturated granular soils in the pendular regime using the Discrete Element Method (DEM) involving a volume-controlled capillary bridge model. In these simulations, water content is maintained constant by redistributing the water from ruptured capillary bridges to adjacent ones. The evolution of capillary stress with varying water contents during triaxial and biaxial loading conditions is systematically examined. The DEM simulation results show that, under both loading conditions, the mean component of the capillary stress generally decreases, while its deviatoricity gradually develops. These changes are observed to become less significant as the initial degree of saturation increases. At low saturation levels, capillary bridges between non-contacting particle pairs rupture due to soil deformations, and the water from these ruptured bridges redistributes to existing contacts. This redistribution leads to an anisotropic distribution of pore water aligned with the contact network. At higher saturation levels, non-contacting capillary bridges persist due to their ability to sustain large relative displacements between particles, allowing the spatial distribution of pore fluids to remain less constrained by the solid contact network. Additionally, at higher water contents, relative sliding and particle rearrangement are the primary factors influencing the directional distribution of capillary bridges.

In this paper, a state-dependent, bounding surface plasticity model that simulates the behavior of unsaturated granular soils is presented. An unsaturated, soil mechanics-compatible elastoplastic response is adopted in which no part of the response occurs in a purely elastic fashion. To create an appropriate hydro-mechanical coupling, a newer generation stress framework, consisting of the Bishop-type effective stress and a second stress variable, is used in conjunction with a soil-water characteristic curve function. Details regarding the model development, parameter estimation, and assessment of the model's predictive capabilities are outlined. With a single set of parameter values, the model realistically simulates the main features that characterize the shear and volumetric behavior of unsaturated granular soils over a wide range of matric suction, density, and net confining pressure.