Understanding the mechanical behaviour of water ice-bearing lunar soil is essential for future lunar exploration and construction. This study employs discrete element method (DEM) simulations, incorporating realistic particle shapes and a flexible membrane, to investigate the effects of ice content, initial packing density, and gravitational conditions on lunar soil behaviour. Initially, we calibrated DEM model parameters by comparing triaxial tests on lunar soil without ice to physical experiments and the angle of repose simulations, validating the accuracy of our approach. Building on this, we conducted simulations on water ice-bearing lunar soil, examining stress-strain responses, shear strain, bond breakage, deviatoric fabric, and N-ring structures. DEM simulations demonstrate that increasing ice content from 0 % to 10 % elevates peak strength from 85 kPa to 240 kPa in loose samples and from 0.2 MPa to 1.62 MPa in dense samples. This strengthening aligns with microstructural stabilization evidenced by 5-ring configurations and narrowed branch vector distributions. Strain field analysis reveals greater deformation magnitudes in icy regolith, suggesting a trade-off between enhanced load-bearing capacity and reduced ductility. These quantified mechanical responses, including strength gain, structural stabilization, and strain localization, reveal the dual engineering implications of water ice in lunar soil.

来源平台：COMPUTERS AND GEOTECHNICS