Evolution of anisotropic capillarity in unsaturated granular media within the pendular regime

While the shear behavior of granular soils is directly related to the microstructure of contacts which often leads to the coaxiality between Cauchy stress and Satake fabric tensors, it is generally accepted by the geomechanics and geotechnical engineering community that the capillary effects are isotropic. At low saturation levels, however, the pore fluid tends to form interparticle menisci that can also manifest an anisotropic structure, which may result in the development of anisotropic capillarity in wetted granular media. To study the interplay between the solid grain contacts and the liquid bridges at the micro-scales, this study adopts a coupled discrete element method that utilizes a linear contact model combined with a capillary model, and explores their effects by conducting a series of numerical experiments. The distributions of contact and capillary force orientations during the experiment are further investigated to better understand how their alignments affect the global response of the granular assembly subjected to a deviatoric loading. The results indicate that the global shear stress response is not only affected by the contact fabric but also by the network of liquid bridges, and we also observe that the particles may lose contact while the pendular menisci may not be destroyed during the elastic unloading.