Microscopic contact model of lunar regolith for high efficiency discrete element analyses

The grains of lunar regolith are characterized with rough surfaces, angular shapes and mutual adhesions due to short-range interactions. These features control the macroscopic mechanical behavior of lunar regolith but have not been completely captured by contact models in previous Discrete Element Method (DEM) analyses. In this paper, a simplified two-dimensional microscopic contact model is proposed for high efficiency DEM analyses of lunar regolith. The model consists of three components in the normal, tangential and rolling directions respectively, plus two new parameters. A shape parameter is used to control the rolling resistance ability at the contact area between two particles to capture the features of grain shape and interlocking. The second parameter, micro-separation, which denotes the nominal minimum distance between the molecules of the two contacting particles, is introduced to account for van der Waals force as the major component of the short-range interactions that contribute to the adhesion of regolith grains in lunar environment conditions. The novel model has been implemented in a two-dimensional DEM code for numerical simulations of biaxial compression tests on lunar regolith. The effects of interparticle friction, grain shape, lunar environment conditions and void ratio on the strength of lunar regolith were numerically investigated. The results show that soils in the simulated lunar environment exhibit greater strength and more apparent strain-softening and shear dilatancy than on the Earth. The proposed model can capture the main features of the mechanical behavior of lunar regolith (apparent cohesion and high peak friction angle) and a wide range of strength indices can be obtained by the contact model.