不同含水率黏重黑土与触土部件互作的离散元仿真参数标定

为了获得可用于东北地区黏重黑土与触土部件相互作用的离散元仿真模拟参数,该文利用EDEM中Hertz-Mindlin with JKR Cohesion接触模型对不同含水率的东北地区黏重黑土进行相关参数标定,针对含水率在10%～20%的实际作业环境,分别配置含水率为12.46%±1%和17.15%±1%的2种黏重黑土,以土壤颗粒间的滚动摩擦系数、恢复系数、表面能参数及静摩擦系数为标定对象,并以土壤颗粒的仿真堆积角为响应值,基于Box-Behnken的响应面优化方法得到堆积角回归模型,并对回归模型进行寻优,得到2种含水率的模型参数优化解,并给出了模型参数范围。测定了4种含水率下黏重黑土对3种触土部件材料(65Mn、UHMW-PE和PTFE)的静摩擦系数,并以此为基础分别对65Mn(典型铁基材料)和PTFE(典型低表面能材料)板进行斜面试验,以含水率为17.15%±1%的黏重黑土为试验对象,分别搭建斜面物理试验平台和仿真模型,以土壤颗粒与触土部件材料之间的滚动摩擦系数、恢复系数、表面能参数及静摩擦系数为标定对象,以仿真得到的土球在65Mn和PTFE板上的滚动距离为响应值,基于响应面优化法得到滚动距离回归模型,以实测的滚动距离为目标对回归模型进行寻优,得到黏重黑土对2类典型触土部件材料接触模型参数的优化解。研究结果表明,标定优化后的土壤模型能够近似代替真实的东北地区黏重黑土进行仿真,可利用标定后的参数进行黏重黑土与触土部件间的离散元仿真,可为东北黏重黑土作业条件下的农业机械触土部件仿生减阻设计与优化提供基础数据。

Calibration of parameters of interaction between clayey black soil with different moisture content and soil-engaging component in northeast China

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.