Abstract: Abstract: In order to obtain the interaction parameters which can be used for discrete element simulation between the soil-engaging components and clayey black soil in northeast China, the Hertz-Mindlin with JKR Cohesion contact model in EDEM2.7 was used to simulate the clayey black soil. Aiming at the actual working environment with moisture content of 10%-20%, 2 kinds of clayey black soil with moisture content of 12.46%±1% and 17.15%±1% were prepared respectively. The rolling friction coefficient, recovery coefficient, JKR surface energy parameter and the static friction coefficient between the soil particles were used as the model calibrated parameters. The soil repose angle simulated under the simulation parameters was set as the response value. The regression model of the soil repose angle was established based on the Box-Behnken response surface method, and the regression model was optimized by using the soil repose angles under 2 moisture contents obtained by physical experiments. The optimal solution of the contact model parameters of the clayey black soil particles with 2 moisture contents was obtained. The parameter range of the discrete element contact model of clayey black soil with moisture content between 10%-20% were given. The soil repose angle and stacked shape obtained by the simulation experiment with the optimized solution had a high similarity with that of the physical test. At the same time, the static friction coefficients of 3 kinds of soil contact materials such as 65Mn plate, UHMW-PE plate and PTFE plate were tested under the conditions of 4 kinds of moisture content, which were 0, 12.46%, 17.15% and 23.5% respectively. Then for the 2 typical soil contact materials such as 65Mn and PTFE, the slope physical test bench were constructed for the clayey black soil with moisture content of 17.15%±1%. The slope simulation models of 65Mn and PTFE materials were established under the EDEM2.7 software platform. Based on the critical simulation parameters such as JKR surface energy parameter, static friction coefficient, rolling friction coefficient and recovery coefficient between soil particle and soil contact material, the simulated rolling distance of the soil ball on different soil contact materials was set as the response value. The excellent fit regression model of rolling distance for the 2 kinds of typical soil contact materials were established based on the response surface optimization method of Box-Behnkende. The 2 regression models were optimized by the sliding distance obtained by the physical experiment. Then the optimized solution of the contact model parameters of different materials was obtained. The simulation results showed that when the surface energy of JKR between soil and 65Mn was 5.5 J/m2, the recovery coefficient was 0.61, the static friction coefficient was 0.57 and the rolling friction coefficient was 0.056, the rolling distance simulation result was 153.56 mm, which was close to the average rolling distance obtained from physical experiments of 155.93 mm, and the relative error was 1.52%; When the surface energy of JKR between soil and PTFE was 4.08 J/m2, the recovery coefficient was 0.6, the static friction coefficient was 0.52 and the rolling friction coefficient was 0.045, the simulation result of rolling distance was 269.35 mm, which was close to the average rolling distance of 269.55 mm obtained from physical experiment, and the relative error was 0.07%. And the optimized parameters can be used to simulate the discrete parameters between the clayey black soil and the soil-engaging components. The study provides credible basic data for the design and simulation of agricultural machinery under clayey black soil conditions.