Calibration of Discrete Element Method Parameters for a High-Fidelity Lunar Regolith Simulant Considering the Effects of Realistic Particle Shape

The Discrete Element Method (DEM) is an important tool for investigating the geotechnical properties of lunar regolith. The accuracy of DEM simulations largely depends on precise particle modeling and the appropriate selection of mesoscopic parameters. To enhance the reliability and accuracy of the DEM in lunar regolith studies, this paper utilized the high-fidelity IRSM-1 lunar regolith simulant to construct a DEM model with realistic particle shapes and conducted an angle of repose (AoR) simulation test. The optimal DEM parameters were calibrated using a combination of the Plackett-Burman test, steepest ascent test, and Box-Behnken design. The results indicate that the sliding friction coefficient, rolling friction coefficient, and surface energy significantly influence the simulation AoR. By optimizing against the measured AoR using a second-order regression model, the optimal parameter values were determined to be 0.633, 0.401, and 0.2, respectively. Under these optimal parameters, the error between the simulation and experimental AoR was 2.1%. Finally, the calibrated mesoscopic parameters were validated through a lifting cylinder test, showing an error of 6.3% between the simulation and experimental results. The high similarity in the shape of the AoR further confirms the accuracy and reliability of the parameter calibration method. This study provides a valuable reference for future DEM-based research on the mechanical and engineering properties of lunar regolith.