The discrete element computer simulation method is an effective tool that enables the study of the interaction mechanism between the pulping components and the paddy soil during the paddy field pulping process. The findings are valuable in optimizing the parameters of the paddy beating device to improve its working quality and efficiency. However, the lack of accurate soil models for paddy soil has limited the application and development of the discrete element method in paddy pulping research. This study selected the Hertz-Mindlin with Johnson-Kendall-Roberts discrete element model for the pre-pulping paddy loam soil and used the slump error as the test index to select nine parameters, including soil Poisson’s ratio and surface energy, as test factors to calibrate the model parameters. The Plackett-Burman test identified soil shear modulus, surface energy, and soil-iron plate static friction coefficient as significant factors affecting the test index. The steepest ascent test results determined the test range of the above parameters. The Box-Behnken test obtained the regression model between the significant factors and the test index, and the regression model was optimized using the slump error as the target. The optimal combination of parameters was surface energy of 3.257 J/m2, soil shear modulus of 0.709 MPa, and static friction coefficient between soil and iron plate of 0.701. The slump simulation test using this combination of parameters yielded an average slump error of 2.04%. The collective results indicate the accuracy of the calibrated discrete element simulation parameters for paddy loam soil. These parameters can be used for discrete element simulation analysis of the paddy pulping process after paddy field soaking.