Parameter Calibration and Experimental Study of a Discrete Element Simulation Model for Yellow Cinnamon Soil in Henan, China

Abstract

To investigate the interaction mechanism between agricultural tillage machinery and soil, this study established a precise simulation model by integrating physical and numerical experiments using typical yellow cinnamon soil collected from western Henan Province, China. The discrete element parameters for soils with varying moisture contents were calibrated based on the Hertz–Mindlin (no slip) contact model. Through Plackett–Burman screening, steepest ascent optimization, and Box–Behnken response surface methodology, a predictive model correlating moisture content, parameters, and repose angle was developed, yielding the optimal contact parameter combination: interparticle static friction coefficient (0.6), soil–65Mn static friction coefficient (0.69), and interparticle rolling friction coefficient (0.358). For the Bonding model, orthogonal experiments coupled with NSGA-II multi-objective optimization determined the optimal cohesive parameters targeting maximum load (673.845 N) and displacement (9.765 mm): normal stiffness per unit area (8.8 × 107 N/m3), tangential stiffness per unit area (6.85 × 107 N/m3), critical normal stress (6 × 104 Pa), critical tangential stress (3.15 × 104 Pa), and bonding radius (5.2 mm). Field validation using rotary tillers and power harrows demonstrated less than 6% deviation in soil fragmentation rates between simulations and actual operations, confirming parameter reliability and providing theoretical foundations for constructing soil-tillage machinery interaction models.