Simulation Analysis and Experiments for Blade-Soil-Straw Interaction under Deep Ploughing Based on the Discrete Element Method

: The desirable sowing period for winter wheat is very short in the rice-wheat rotation areas. There are also lots of straw left in harvested land. Deep rotary tillage can cover rice straw under the surface to increase soil organic matter. Clarifying the effect of the rotary tillage blade on the soil and straw, as well as analyzing the movement patterns and forces on the straw and soil, are essential to investigate the deep rotary tillage process in order to solve the problems of energy consumption and poor straw burial effect of deep tillage and deep burial machinery. In this study, we built the interaction model of rotary blade-soil-straw through the discrete element method to conduct simulation and identiﬁed the factors that affect the power consumption and operation quality of the rotary blade. The simulation process reﬂects the law of rotary blade-soil-straw interaction, and the accuracy of the simulation model has been veriﬁed by ﬁeld trials. The simulation test results show that the optimized structural parameters of the rotary tillage blade were 210 mm, 45 mm, 37 ◦ and 115 ◦ (R, H, α and β ) designed based on this theoretical model can cultivate to a depth of 200 mm. The operating parameters were 8 π rad/s for rotational speed and 0.56 m/s for forward speed, respectively; the simulated and ﬁeld comparison tests were conducted under the optimal combination of parameters, and the power, soil breaking rate, and straw burial rate were 1.73 kW, 71.34%, and 18.89%, respectively; the numerical error rates of simulated and ﬁeld test values were 6.36%, 5.42%, and 8.89%, respectively. The accuracy of the secondary model was veriﬁed. The simulation model had good accuracy at all factor levels. The model constructed in this study can provide a theoretical basis and technical reference for the interaction mechanism between rotary tillage and soil straw, the optimization of machine geometry, and the selection of operating parameters.


Introduction
Rotary tillers are widely used in the seed bed preparation segment of agricultural production [1]. The study demonstrated that a rotation tillage depth of 200 mm could fully mix the straw into the soil and significantly improve the sowing quality and wheat yield [2]. However, there are a significant number of challenges, such as insufficient tillage depth (less than 150 mm), unsatisfactory straw burying effect (low straw burying rate, insufficient straw burying depth), high power consumption, etc., in rice-wheat rotation areas. Moreover, long-term shallow rotational tillage results in a thinning and low-fertility soil tillage layer, which in turn leads to lower crop yields [3]. Simultaneously, inadequate straw mulching resulted in the distribution of straw in the upper part of the tillage layer, which allowed gaps between the top and bottom of the soil tillage layer, thus leading to segregation of the soil. It is easy to cause the next crop to fall dry and die, reducing the seeding rate in the cultivation area of two crops a year [4]. The concentration of a huge number of agricultural machines in the busy season, coupled with the increase in power consumption, is prone to air contamination [5,6], which is not in accordance with the concept of green development of agricultural machinery and decelerates the achievement of the carbon peaking and carbon neutral development goals.
To meet the above challenges and realize deep plowing, four solutions are applied in production practice as follows. (1) Increase the radius of the rotary blade. The test proved that the machine could achieve a tillage depth of 300 mm, but there are problems with the huge mechanism, heavy tillage load, and low reliability of the machine [7]. (2) Double-axis front shallow and behind deep rotary tillage. The machine sets two rows of knife rollers at the front and behind, the former for the first operation, the latter again, and finally, two operations to achieve deep tillage. In this solution, the soil accumulates in the gap between the front and rear blade axes of the machine, which interferes with the soil and prevents the rear blade axes from sinking into the ground. Moreover, the structure is complicated and heavy and consumes much energy [8,9]. (3) Submerged soil reversal rotary tillage. This solution is a tillage method in which the cutter shaft is reversed and submerged below the ground surface. This form of mechanism resulted in the transmission box and frame being difficult to enter the soil and congestion problems [10][11][12][13][14][15][16][17]. (4) Diagonal submerged soil reversal rotary tillage. This solution effectively solves the problem that the transmission box cannot enter the soil in the above two solutions. The blade roller is angled in the horizontal plane during the operation of the machine, and tillage is carried out obliquely, but there is a significant lateral force, which is not favorable for the operator to control the direction [18,19].
The structure of key components is unreasonable, and the equipment design lacks academic guidance. The above methods restrict the research progress of deep rotary tillage. The systematic study of rotary blade-soil-straw interaction in deep rotary tillage operation is insufficient, and the design of rotary tiller parameters and whole machine design lack theoretical support. This study synthesizes the research results in related fields of agricultural engineering and conducts research and discussion on the mechanism of rotary blade-soil-straw interaction in deep rotary tillage operation based on the discrete element method. The Discrete Element Method (DEM) is a numerical simulation method for basic research on solving discontinuous media problems [20]. It is different from the Finite Element Method (FEM) for solving continuous media problems. The basic principle of this method is to separate the research object into a set of rigid units. In order to make each unit conform to Newton's second law, the motion equation of each unit is calculated by means of central difference, and the whole motion state and parameters of the research object are obtained. This method is favored by most researchers in the field of agricultural engineering. In this paper, the discrete element method is used to deal with the fact that the soil model is a discontinuous medium [21,22].
We have developed a discrete element model of rotary blade-soil-straw for issues such as soil breakage rate, straw displacement, and power consumption during observation of rotary tillage operation. After that, the prototype was tested according to the 3D model of the machine designed on the basis of this study, and field verification tests were done. Finally, the validity of the model was determined by comparing the results of field tests and simulation analysis. The structural arrangement of the remainder of this article is as follows: Part 2 describes the analysis of the relevant phenomena during the simulation test. Part 3 shows the relevant results. Part 4 presents the discussion of the test. Part 5 concludes and summarizes the whole article. In this study, a blade-soil-straw interaction model was established for rice-wheat rotation. The model can provide a theoretical basis and technical reference for the interaction mechanism between rotary tillage and soil straw, the optimization of machine geometry, and the selection of operating parameters. The tractor is connected to the rotary tiller by a towing device and transmits power through a universal coupling. The schematic diagram of the rotary tiller is shown in Figure 1a. The rotary tiller moves forward with a tractor. The rotary tiller shaft rotates through the gearbox, following the tractor power output shaft to drive the rotary tiller. The scraper impacts and levels the scattered soil to achieve a flat cultivated land. In order to study the interaction rule of rotary tillage blade-soil-straw more intuitively, this paper first studies the operation mechanism of rotary tillage single blade-soil-straw and then optimizes the form and structure parameters of rotary tillage single blade.

Structure of Rotary Blades
The tractor is connected to the rotary tiller by a towing device and transmits power through a universal coupling. The schematic diagram of the rotary tiller is shown in Figure 1a. The rotary tiller moves forward with a tractor. The rotary tiller shaft rotates through the gearbox, following the tractor power output shaft to drive the rotary tiller. The scraper impacts and levels the scattered soil to achieve a flat cultivated land. In order to study the interaction rule of rotary tillage blade-soil-straw more intuitively, this paper first studies the operation mechanism of rotary tillage single blade-soil-straw and then optimizes the form and structure parameters of rotary tillage single blade.
The blade is the key working part of the rotary tiller. The shape and structural parameters of the blade badly affected the quality and power consumption of the rotary tiller. The blade in this paper is designed and researched based on a standard rotary tillage cutter according to the research objectives, as shown in Figure 1b. The remaining are implemented in accordance with GBT 5669-2017 Rotary Tillage Machine Blades and Holders. The main structural parameters are shown in Table 1.  The values of forward velocity and linear velocity of the blade endpoint determine the normal progress of soil cutting operation. The blade endpoint presents a composite motion; therefore, arrangements of velocity parameters will have diverse trajectories.
Let Vd-the tangential velocity of the blade endpoint, m/s;  Table 1. The values of forward velocity and linear velocity of the blade endpoint determine the normal progress of soil cutting operation. The blade endpoint presents a composite motion; therefore, arrangements of velocity parameters will have diverse trajectories.
Let V d -the tangential velocity of the blade endpoint, m/s; The blade velocity ratio λ and the motion trajectory are shown in Figure 2. (1) When λ < 1, that is, the forward velocity V m is greater than the blade linear velocity V d , its motion trajectory is a short cycloid, the direction of the horizontal component of the linear velocity at any point on the cycloid is consistent with the forward direction, and it is hardly difficult to realize the operation of post throwing soil blocks.
(2) When λ = 1, that is, the two values are equal, its motion trajectory is a cycloid curve, The horizontal velocity at any point is zero, and the same as above cannot realize soil throwing.
(3) When λ > 1, that is, the forward velocity V m is less than the blade's linear velocity V d , its motion trajectory is a complementary cycloid. The horizontal partial velocity direction of any point below the maximum cross chord is opposite to the forward direction, exact converse above, the soil block can be reprojected effectively at present. The blade velocity ratio λ and the motion trajectory are shown in Figure 2.
(1) When λ < 1, that is, the forward velocity Vm is greater than the blade linear velocity Vd, its motion trajectory is a short cycloid, the direction of the horizontal component of the linear velocity at any point on the cycloid is consistent with the forward direction, and it is hardly difficult to realize the operation of post throwing soil blocks.
(2) When λ = 1, that is, the two values are equal, its motion trajectory is a cycloid curve, The horizontal velocity at any point is zero, and the same as above cannot realize soil throwing.
(3) When λ > 1, that is, the forward velocity Vm is less than the blade's linear velocity Vd, its motion trajectory is a complementary cycloid. The horizontal partial velocity direction of any point below the maximum cross chord is opposite to the forward direction, exact converse above, the soil block can be reprojected effectively at present.
where R is the radius of gyration of a blade, with the unit of m, Vm is forward velocity, with the unit of m/s, ω is angular velocity, with the unit of rad/s, t is time, with the unit of s, h is tillage depth, with the unit of m.
To meet the condition of throwing soil backward, the horizontal partial velocity Vx of the absolute velocity at any point on the absolute motion track of the blade is less than 0, according to the above equation, let: The movement locus of the blade endpoint is shown in Figure 2b. As can be seen from the above equation: R increases, H increases, but the increased torque leads to increased power consumption, which leads to a decrease in speed and, therefore, productivity. Therefore, generally, R = 215, 245, 260 mm, and in this study, a range of 150 to 240 Set point M is set as the blade cutting point; it will meet the normal working conditions of the rotary cultivator from the beginning of soil penetration to the end of soil throwing and leaving the ground; there are: where R is the radius of gyration of a blade, with the unit of m, V m is forward velocity, with the unit of m/s, ω is angular velocity, with the unit of rad/s, t is time, with the unit of s, h is tillage depth, with the unit of m.
To meet the condition of throwing soil backward, the horizontal partial velocity V x of the absolute velocity at any point on the absolute motion track of the blade is less than 0, according to the above equation, let: The movement locus of the blade endpoint is shown in Figure 2b. As can be seen from the above equation: R increases, H increases, but the increased torque leads to increased power consumption, which leads to a decrease in speed and, therefore, productivity. Therefore, generally, R = 215, 245, 260 mm, and in this study, a range of 150 to 240 mm was chosen in order to explore whether a small radius could be used for rototilling under ground level. The discrete element parameters consist of the intrinsic parameters (material density, Poisson's ratio, shear modulus), material contact parameters (coefficient of restitution, coefficient of static friction, dynamic friction coefficient), and contact model parameters. The parameters employed in this article are taken from the reference [21]. The relevant parameters are shown in Table 2.
The increase in the number of particles in the software will lead to a decline in the geometric level of the simulation calculation speed [22]. In addition, too small particles will lead to an uncertain impact on the model under the micro force. Therefore, the size of the soil model is generally larger than the real soil particles in the simulation calculation. Without affecting the simulation test, the soil particle radius selected in this paper is 5 mm, 10 mm, 15 mm, and 30 mm, as shown in Figure 3. mm was chosen in order to explore whether a small radius could be used for rototilling under ground level.

Soil and Straw Model Parameters
The discrete element parameters consist of the intrinsic parameters (material density, Poisson's ratio, shear modulus), material contact parameters (coefficient of restitution, coefficient of static friction, dynamic friction coefficient), and contact model parameters. The parameters employed in this article are taken from the reference [21]. The relevant parameters are shown in Table 2.
The increase in the number of particles in the software will lead to a decline in the geometric level of the simulation calculation speed [22]. In addition, too small particles will lead to an uncertain impact on the model under the micro force. Therefore, the size of the soil model is generally larger than the real soil particles in the simulation calculation. Without affecting the simulation test, the soil particle radius selected in this paper is 5 mm, 10 mm, 15 mm, and 30 mm, as shown in Figure 3.  In this paper, the rotary tillage blade mainly covers the straws scattered on the ground during operation and will not crush and cut the straws. For the purpose of improving the simulation efficiency, the straws established in this paper are also simplified to a certain extent. It is mainly composed of 10 particles with a radius of 5 mm. The center distance between adjacent particles is 5 mm. The total length of the straw model is 55 mm. The simulation model of straw is shown in Figure 4.

Contact Model Settings
Edem software 2018 integrates many options to apply to different material types of contact. We selected Hertz-Mindlin with Bonding along with the characteristics that the soil particles are difficult to reconnect after fracture in this paper. At the same time, the In this paper, the rotary tillage blade mainly covers the straws scattered on the ground during operation and will not crush and cut the straws. For the purpose of improving the simulation efficiency, the straws established in this paper are also simplified to a certain extent. It is mainly composed of 10 particles with a radius of 5 mm. The center distance between adjacent particles is 5 mm. The total length of the straw model is 55 mm. The simulation model of straw is shown in Figure 4.

Soil and Straw Model Parameters
The discrete element parameters consist of the intrinsic parameters (mater Poisson's ratio, shear modulus), material contact parameters (coefficient of rest efficient of static friction, dynamic friction coefficient), and contact model param parameters employed in this article are taken from the reference [21]. The releva eters are shown in Table 2.
The increase in the number of particles in the software will lead to a dec geometric level of the simulation calculation speed [22]. In addition, too sma will lead to an uncertain impact on the model under the micro force. Therefor of the soil model is generally larger than the real soil particles in the simulatio tion. Without affecting the simulation test, the soil particle radius selected in th 5 mm, 10 mm, 15 mm, and 30 mm, as shown in Figure 3.  In this paper, the rotary tillage blade mainly covers the straws scatter ground during operation and will not crush and cut the straws. For the purp proving the simulation efficiency, the straws established in this paper are also to a certain extent. It is mainly composed of 10 particles with a radius of 5 mm. distance between adjacent particles is 5 mm. The total length of the straw mode The simulation model of straw is shown in Figure 4.  Edem software 2018 integrates many options to apply to different materi contact. We selected Hertz-Mindlin with Bonding along with the characteristi soil particles are difficult to reconnect after fracture in this paper. At the sam

Contact Model Settings
Edem software 2018 integrates many options to apply to different material types of contact. We selected Hertz-Mindlin with Bonding along with the characteristics that the soil particles are difficult to reconnect after fracture in this paper. At the same time, the most commonly used contact theoretical model Hertz-Mindlin (no slip), is selected for ordinary contact calculation. Both normal and tangential forces have damping components, where the damping coefficient is related to the coefficient of restitution. The tangential friction force follows the Coulomb friction law, and the rolling friction adopts the research results of Sakaguchi. The relevant calculation principles are described in the literature [23][24][25][26][27][28].
The normal force F n is expressed as follows: where F n is the normal force, with the unit of N, E* is equivalent to Young's modulus, with the unit of Pa, R* is the equivalent radius, with the unit of m, δ n is the amount of normal overlap, with the unit of m, v i,j is Poisson's ratio, E i,j is Particle Young's modulus, with the unit of Pa, R i,j is Particle radius, with the unit of m. The normal damping force F n d is expressed as follows: ln 2 e + π n 2 (12) where F n d is Normal damping force, with the unit of N, β is viscous damping coefficient, e is coefficient of restitution, Sn _ is normal stiffness, with the unit of N/m, m* is equivalent mass, with the unit of g, m i,j is Particle mass, with the unit of g, v → rel n is Normal relative velocity, with the unit of m/s. The expression for the tangential force F t is as follows: where F t is the normal force, N,S t is tangential stiffness, with the unit of N/m, δ t is tangential overlap, with the unit of m, G* is equivalent shear modulus, with the unit of Pa. The tangential damping force F t d is expressed as follows: where F t d is angential damping force, with the unit of N, S t is tangential stiffness, with the is Tangential relative velocity, with the unit of m/s. The total tangential force is limited by the Coulomb friction force. When the tangential force reaches the Coulomb friction force, that is, F t > µ s F n , relative sliding occurs between the particles, where s is the static friction coefficient. To simplify the contact model, The model does not take into account the sliding friction of the contact model. The rolling friction is expressed in terms of the contact torque of the particle model, and the expression is as follows: where τ i is torque, with the unit of N·m, µ r is coefficient of rolling friction, R i is the distance from the contact point to the centroid, with the unit of m, ω i -angular velocity of particles at the contact point, with the unit of r/min.

(b) Hertz-Mindlin with Bonding
Hertz-Mindlin with Bonding contact model bonds particles together. Since this model includes Hertz-Mindlin (no slip), the default Hertz-Mindlin (no slip) contact model needs to be deleted from the software Creator Tree list when generating the bond.
Hertz-Mindlin with Bonding contact model uses finite-sized "bonding" bonds to bond particles. The bond can withstand the normal and tangential motion between the contacting particles to generate normal and tangential shear stress until the bond breaks at the critical maximum stress. The particles interact as rigid spheres, and no new bonds are generated, which is in line with the soil stress mechanical behavior characteristics of remaining loose after crushing [29][30][31].
In this model, the first generated particles interact with each other according to the Hertz-Mindlin contact model, and the connection starts after the set bond generation time point. After the particles are connected, the force and torque between the particles are both set to 0, according to the following Stepwise adjustment at each time step: where F n,t is normal and tangential force, with the unit of N, v n,t is normal and tangential velocity, with the unit of m/s, S n,t is normal and tangential stiffness, with the unit of N/m, A is the cross-sectional area of bonding bond, with the unit of m 2 , δ t is time step, with the unit of s, M n,t is normal and tangential moments, with the unit of N·m, ω n,t is normal and tangential angular velocities, with the unit of r/min, J is the moment of inertia, with the unit of m 4 , R B is the bond-bond radius. When the normal and tangential stresses exceed preset critical values, the bonded bonds will break, the particles will interact as rigid spheres, and no new bonds will be formed. Spawn bonds in this model can create cohesive forces when particles are generated, even if the particles are not in contact, so the contact radius should be set larger than the radius of the particle body. Table 3 lists the parameters adopted in this study. In this article, we use Edem's dynamic generation method to generate soil particle beds layer by layer, and the model is shown in Figure 5.
The particle contact model can quantify the interaction between discrete elements within a given contact range. Adding a contact model to the particle model is a key step in the test. Appropriate contact model to accelerate simulation test.   Figure 5. The discrete element model.

Initial Conditions for Simulation
In this simulation test, in order to ensure that the tillage depth was equal t the rotation center was set to the corresponding position into the soil accord radius of rotation of each blade, as shown in Figure 6. According to the research objectives of this article, combined with the e search experience and actual working conditions, the main indicators were dete the power (P), the rate of crushed soil (RC), and the rate of straw burial (RB).

Initial Conditions for Simulation
In this simulation test, in order to ensure that the tillage depth was equal to 200 mm, the rotation center was set to the corresponding position into the soil according to the radius of rotation of each blade, as shown in Figure 6. within a given contact range. Adding a contact model to the particle model is a ke in the test. Appropriate contact model to accelerate simulation test.

Initial Conditions for Simulation
In this simulation test, in order to ensure that the tillage depth was equal to 20 the rotation center was set to the corresponding position into the soil according radius of rotation of each blade, as shown in Figure 6.  According to the research objectives of this article, combined with the existing research experience and actual working conditions, the main indicators were determined as the power (P), the rate of crushed soil (RC), and the rate of straw burial (RB). Figure 1b shows that the value of B is fixed when R and β are certain, so the factor B is removed. The initial conditions of the experiment are shown in Table 4, and the factor level design is shown in Table 5. In this study, for convenience of calculation, the ratio of the number of bond-bond breaks per unit volume to the number of bond bonds at the beginning is used as the index of soil breakage. The sampling position is set 500 mm from the center of rotation after rotary tiller operation, and 1000 mm is set by grid-grid function after EDEM treatment. 1000 mm × 500 mm × 250 mm grid solver calculates the bonding bond between particles as a problem to be solved.
The value of power (P) is selected from the data within 2 s of the operation stabilization stage. Each blade is tested three times, and the average value is obtained, which is directly exported after EDEM software post-processing.
The crushed soil rate (RC) is measured by the ratio of bond breakage between particles in the tillage area. In this study, for the convenience of calculation, the ratio of the number of soil bond breakage per unit volume to the number of bonded bonds at the initial time was used as the indicator of crushed soil rate. The sampling position was set at 500 mm from the center of rotation after the rotary tillage knife operation, and the grid solver of 1000 mm × 500 mm × 250 mm was set using the grid function of EDEM post-processing, and the particle-to-particle bond was calculated as the problem to be solved.
Straw burial rate (RB) is the percentage of straw residual mass on the ground posttillage to the pre-tillage. In the post-processing of the EDEM software, L0 is the straw layer. The grid solver was divided into four layers (L1, L2, L3, L4) in the vertical direction, and the percentage of the total straw mass from L1 to L4 to the initial L0 mass was calculated as the straw burial rate.

Simulation
To better observe the movement of soil and straw in the simulation, we selected the velocity clouds from the simulation with different angles of blade rotation and placed them in Figure 7.
The velocity of the model is represented by the color-changing legend. The velocity ranges from green to red is 0-1 m/s, which means the velocity is the same as the forward direction. The velocity ranges from green to blue is 0-−1 m/s, which means the velocity is opposite to the forward direction. L refers to the whole tillage layer, L0 represents the straw, and L1-L4 represents the soil from 0-50 m, 50-100 mm, 100-150 mm, and 150-200 mm below the horizontal surface, respectively.
It can be seen from L: At 0 • the rotary blade touches the soil, the soil and straw on the side edge are first disturbed by the force of the rotary blade, L1 soil, and L0 straw splash. The disturbed area has a definite initial velocity and appears yellow, representing a low value; At 90 • , the rotary blade has fully entered the soil. The color of the soil and straw starts to change, which indicates that the number of particles gaining velocity starts to increase. At L, L0, and L1, it is clearly observed that the soil and straw near the rototiller knife are red and dark blue in color, with the majority of blue particles, which indicates that the velocity starts to increase and most of the soil acquires the initial velocity backwards.
Concurrently, the color of the soil and straw on the other side was light blue, indicating that the latter was significantly more disturbed than the former, and the velocity was negative. At 180 • , the rotary blade was about to break out of the soil. The soil near the front section of the rotary tiller is dark blue and blue, which means that the rotary tiller throws the soil in the lower part of the plow layer backwards and upwards. The color of the soil near the side edge of the rototiller is red, illustrating that the soil in this position moves forward with the rototiller. At this time, the color of the soil with greater disturbance on the left side changes to yellow, indicating that this part of the soil starts to move forward and gradually returns to its original state. The color of the thrown soil and straw particles changed after the rotary tillage knife came out of the soil, indicating that the soil and straw continued to move in the direction of the speed obtained when thrown and finally fell into the ground. The color of the disturbed area is dark green, demonstrating that the velocity returns to 0.
L0 and L1 received a significant disturbance. The soil close to the center is taken to L2 and L3 when the blade starts to cut the soil, and the soil far from the center is thrown upward and backward and finally falls into L2 and L4. L2 and L3 are disturbed similarly to L1 and L0. The soil far from the center is thrown out and finally falls into L1 and L3; For L4, the soil disturbance is minimal, but a small part of the soil far from the center is still affected by the force and acquires a certain initial velocity and flow. starts to change, which indicates that the number of particles gaining velocity starts to increase. At L, L0, and L1, it is clearly observed that the soil and straw near the rototiller knife are red and dark blue in color, with the majority of blue particles, which indicates that the velocity starts to increase and most of the soil acquires the initial velocity backwards. Concurrently, the color of the soil and straw on the other side was light blue, indicating that the latter was significantly more disturbed than the former, and the velocity was negative. At 180°, the rotary blade was about to break out of the soil. The soil near the front section of the rotary tiller is dark blue and blue, which means that the rotary tiller throws the soil in the lower part of the plow layer backwards and upwards. The color of the soil near the side edge of the rototiller is red, illustrating that the soil in this position moves forward with the rototiller. At this time, the color of the soil with greater disturbance on the left side changes to yellow, indicating that this part of the soil starts to move forward and gradually returns to its original state. The color of the thrown soil and straw particles changed after the rotary tillage knife came out of the soil, indicating that the soil and straw continued to move in the direction of the speed obtained when thrown and finally fell into the ground. The color of the disturbed area is dark green, demonstrating that the velocity returns to 0.
L0 and L1 received a significant disturbance. The soil close to the center is taken to L2 and L3 when the blade starts to cut the soil, and the soil far from the center is thrown upward and backward and finally falls into L2 and L4. L2 and L3 are disturbed similarly to L1 and L0. The soil far from the center is thrown out and finally falls into L1 and L3; For L4, the soil disturbance is minimal, but a small part of the soil far from the center is still affected by the force and acquires a certain initial velocity and flow. The migration of the soil between L1 and L4 gives a side view of the mixing state of the tillage layer after rotary tillage. The more the layers migrate, the further they are from The migration of the soil between L1 and L4 gives a side view of the mixing state of the tillage layer after rotary tillage. The more the layers migrate, the further they are from their original position and the better the mixing of the tillage layer. When the rototiller knife enters the soil again, it continues to disturb the soil after tillage, causing the straw to rebound and thus causing it to continue to migrate in the direction of the initial velocity. The soil is significantly displaced by the action of the rotary cutter, and the movement of the soil particles causes more significant displacement of the straw as the mechanical properties of the straw material are lower than those of the soil particles. The movement of the soil particles causes the straw to move more significantly due to the lower mechanical properties of the straw material than the soil particles. When the force is greater than the ultimate strength of the bond, the soil particles are separated and fractured, which means that the soil particles that were originally connected by the bond are separated and move to a different layer; as the tillage continues to advance, some of the particles will be separated by secondary cutting. During the tillage process, the soil is disturbed in different areas of the rototiller, with a clear range of disturbance away from the center. The surface of the tilled soil shows some undulations, mainly due to the fact that the tillage process has changed the surface form of the original soil, and there is a certain friction effect between the soil particles, which makes it difficult to recover the initial form.

Field Trials
In a comparison of the power consumption and operational effects of EDEM simulation and actual soil entry conditions, in field trials, the length of the experiment area is 200 m, and the width is 15 m. The experiment area belongs to the plain area with two crops a year, with an average annual precipitation of about 600 mm, and the soil texture is loam. The first crop of the experiment field is wheat, the vegetation coverage before tillage is 302.4 g/m 2 , the soil firmness of 0-250 mm thick is 2098.0 kPa, the absolute soil moisture content is 21.63%, the average volume mass is 1.67 g/cm 3 , and the temperature during the experiment is 28 • C. We connected the torque sensors and the farm machinery test box to the rotary tiller and the tractor and debugged them to ensure that the test system could complete the work normally, as shown in Figure 8.
the ultimate strength of the bond, the soil particles are separated an means that the soil particles that were originally connected by the bond move to a different layer; as the tillage continues to advance, some of t separated by secondary cutting. During the tillage process, the soil is ent areas of the rototiller, with a clear range of disturbance away fr surface of the tilled soil shows some undulations, mainly due to the process has changed the surface form of the original soil, and there effect between the soil particles, which makes it difficult to recover the

Field Trials
In a comparison of the power consumption and operational effect tion and actual soil entry conditions, in field trials, the length of the exp m, and the width is 15 m. The experiment area belongs to the plain are year, with an average annual precipitation of about 600 mm, and the s The first crop of the experiment field is wheat, the vegetation covera 302.4 g/m 2 , the soil firmness of 0-250 mm thick is 2098.0 kPa, the abs content is 21.63%, the average volume mass is 1.67 g/cm 3 , and the temp experiment is 28 °C. We connected the torque sensors and the farm ma the rotary tiller and the tractor and debugged them to ensure that th complete the work normally, as shown in Figure 8. The power (P) is measured by a torque sensor mounted between coupling. The power (P) is measured by a torque sensor mounted between the tractor and the coupling.
The rate of soil crushed (RC) is determined by measuring a 200 mm × 200 mm × 200 mm full tilled layer of soil on ploughed land. The size of the sampled soil is divided into three classes according to their side lengths: 0-4 cm, 4-8 cm, and 8 cm or more, and the percentage of soil mass from 0-4 cm to the total mass is the rate of soil fragmentation, measured one point per stroke.
The straw burial rate (RB)was determined by separating and weighing the straw from the soil at 0-5 cm, 5-10 cm, 10-15 cm, 15-20 cm, and 20-25 cm. During the sampling process, samples were taken strictly in accordance with the sampling frame and the boundaries of the tillage layer. Straw beyond the boundaries was cut off with scissors, keeping only the straw within the boundaries of the sampling area, and finally, the weight of each layer was used for statistical analysis.
The rotary tiller was connected to the tractor through a three-point suspension, and the torque sensor was installed at the end of the power take-off shaft. The rotary tiller was connected to the tractor through a three-point suspension, and the torque sensor was installed at the end of the power take-off shaft. The rototiller is mounted in the middle of the blade shaft, and to facilitate the same experimental conditions as the simulation, the soil is treated as a monopoly with a height of 250 mm and a width of 200 mm to facilitate the entry of a 150 mm radius rototiller stick into the soil. The structure diagram is shown in Figure 9.
The rate of soil crushed (RC) is determined by measuring a mm full tilled layer of soil on ploughed land. The size of the sam three classes according to their side lengths: 0-4 cm, 4-8 cm, an percentage of soil mass from 0-4 cm to the total mass is the ra measured one point per stroke.
The straw burial rate (RB)was determined by separating a from the soil at 0-5 cm, 5-10 cm, 10-15 cm, 15-20 cm, and 20-25 c process, samples were taken strictly in accordance with the s boundaries of the tillage layer. Straw beyond the boundaries w keeping only the straw within the boundaries of the sampling are of each layer was used for statistical analysis.
The rotary tiller was connected to the tractor through a thre the torque sensor was installed at the end of the power take-off sh connected to the tractor through a three-point suspension, and t stalled at the end of the power take-off shaft. The rototiller is m the blade shaft, and to facilitate the same experimental condition soil is treated as a monopoly with a height of 250 mm and a widt the entry of a 150 mm radius rototiller stick into the soil. The stru in Figure 9.

Results of Box-Behnken Design
Simulation tests were conducted according to the Box-Be power (P) and the rate of crushed soil (RC), and the rate of straw ation indexes, and the test results are shown in Table 6.

Results of Box-Behnken Design
Simulation tests were conducted according to the Box-Behnken design, with the power (P) and the rate of crushed soil (RC), and the rate of straw burial (RB) as the evaluation indexes, and the test results are shown in Table 6. According to the data samples in Table 6, the quadratic model of power (P), rate of crushed soil (RC), and straw burial rate (RB) was calculated by Design-Expert 12.0 software: According to the ANOVA results in Table 7, the Model F-value of 135.01 implies the model is significant. There is only a 0.01% chance that an F-value this large could occur due to noise. p-values less than 0.0500 indicate model terms are significant. In this case, A, B, D, A 2 , B 2 are significant model terms. Values greater than 0.1000 indicate the model terms are not significant. Removing the insignificant term, the regression equation of the power consumption model in this study.
The ANOVA results showed that the interaction of the factors in this model was not significant, and only the radius (R), the end cutter height (H), and the bending angle (β) were significant among the factors. It can be found from Figure 10a that when H = 45, α = 45, β = 120, the power consumption value is the lowest when R is 180, which is 1.12 kW, and as the radius R gradually increases, the power consumption increases significantly, and the analysis of the reason is due to the increase of torque, which leads to the increase of power consumption; when the radius R is lower than 180 and starts to decrease, the reason for the increase of power consumption is that the center of the blade sinks too deep into the tillage layer, which causes the blade to repeatedly Cutting soil, unable to throw out the soil fully, causing an increase in power consumption. Figure 8b,c show that the power consumption decreases with the increase of H and increases with the decrease of β, respectively. The ANOVA results showed that the interaction of the factors in this model was not significant, and only the radius (R), the end cutter height (H), and the bending angle (β) were significant among the factors. It can be found from Figure 10a that when H = 45, α = 45, β = 120, the power consumption value is the lowest when R is 180, which is 1.12 kW, and as the radius R gradually increases, the power consumption increases significantly, and the analysis of the reason is due to the increase of torque, which leads to the increase of power consumption; when the radius R is lower than 180 and starts to decrease, the reason for the increase of power consumption is that the center of the blade sinks too deep into the tillage layer, which causes the blade to repeatedly Cutting soil, unable to throw out the soil fully, causing an increase in power consumption. Figure 8b,c show that the power consumption decreases with the increase of H and increases with the decrease of β, respectively. •

Rate of crushed soil (RC)
According to the ANOVA results in Table 8, the Model F-value of 4.00 implies the model is significant. There is only a 0.70% chance that an F-value this large could occur due to noise. p-values less than 0.0500 indicate model terms are significant. In this case, A, B, A 2 , B 2 are significant model terms. Values greater than 0.1000 indicate the model terms are not significant. Removing the insignificant term, the regression equation of the rate of crushed soil model in this study. •

Rate of crushed soil (RC)
According to the ANOVA results in Table 8, the Model F-value of 4.00 implies the model is significant. There is only a 0.70% chance that an F-value this large could occur due to noise. p-values less than 0.0500 indicate model terms are significant. In this case, A, B, A 2 , B 2 are significant model terms. Values greater than 0.1000 indicate the model terms are not significant. Removing the insignificant term, the regression equation of the rate of crushed soil model in this study. RC = 67.31 + 4.82A + 5.97B − 5.52A 2 − 6.22B 2 (25)  The ANOVA results showed that the interaction of the factors in this model was not significant. The radius (R) and the end cutter height (H) were significant among the factors. It can be found from Figure 11 that when H = 45, α = 45, β = 120, RC gradually increases with the increase of R and H. The curve stops growing at R = 195 and H = 45, respectively. The value of RC hovers around 68%. It indicates that the space for disturbing the soil at the same tillage depth under the action of a single blade is limited, and the scanning trajectory formed by its rotation is approximated. The ANOVA results showed that the interaction of the factors in this model was not significant. The radius (R) and the end cutter height (H) were significant among the factors. It can be found from Figure 11 that when H = 45, α = 45, β = 120, RC gradually increases with the increase of R and H. The curve stops growing at R = 195 and H = 45, respectively. The value of RC hovers around 68%. It indicates that the space for disturbing the soil at the same tillage depth under the action of a single blade is limited, and the scanning trajectory formed by its rotation is approximated. •

Rate of straw burial (RB)
According to the ANOVA results in Table 9, the Model F-value of 11.95 implies the model is significant. There is only a 0.01% chance that an F-value this large could occur •

Rate of straw burial (RB)
According to the ANOVA results in Table 9, the Model F-value of 11.95 implies the model is significant. There is only a 0.01% chance that an F-value this large could occur due to noise. p-values less than 0.0500 indicate model terms are significant. In this case, A, AD, BC, A 2 , D 2 are significant model terms. Values greater than 0.1000 indicate the model terms are not significant. Removing the insignificant term, the regression equation of the rate of crushed soil model in this study is: The ANOVA results showed that the interaction of the factors in this model was significant. The R, Rβ, and Hα were significant among the factors. The trend of RB with R is shown in Figure 12a, which shows that the RB value increases and then decreases, with a maximum at R = 200 and a model prediction of 17.78%. This is due to the upper limit of the maximum straw burial capacity of a single blade role at the same tillage depth. Too small a blade radius causes the straw not to be fully tilled into the ground, and conversely, too large a radius causes the straw to be thrown farther by the blade, also preventing effective mulching. Figure 12b presents the interaction between R and β. It is clear that the effect of changing R on RB is similar to that shown in Figure 12a. Changing the value of β, when β is 110, the maximum value of RB is slightly less than that at 120. When β increases to 130, the extreme value of the variation curve decreases further, and the value of R is taken to 217.5. Figure 12c presents the interaction curve of H with α. It can be seen that the magnitude of α affects the trend of RB with H. This is because α affects the position of the bending line and, thus, the direction of the end line. α is 35, and RB decreases with increasing H. It can be seen that the larger H is at this angle, the more straw can be buried by a single blade. On the contrary, for α of 55, RB behaves in the opposite way to the above. The intersection points of the two are near H = 45.

Parameter Optimization and Comparison Test
The optimization function in Design-Expert 12.0 software was applied for the quadratic model in this study. The optimal parameters of the blade were obtained by solving the regression model with the conditions that the broken soil rate was taken as the minimum value of 65%, the minimum straw burial rate was 15%, and the power consumption was the lowest, as shown in Figure 13

Parameter Optimization and Comparison Test
The optimization function in Design-Expert 12.0 software was applied for the quadratic model in this study. The optimal parameters of the blade were obtained by solving the regression model with the conditions that the broken soil rate was taken as the minimum value of 65%, the minimum straw burial rate was 15%, and the power consumption was the lowest, as shown in Figure 13. The theoretical optimal parameters were obtained as the combination of radius (R) 209.761 mm, end face height (H) 45.8743 mm, pinch angle (α) 37.607°, and bending angle (β) 113.209°. In accordance with the actual fabrication process and practice, the values were taken as R = 210 mm, H = 45 mm, α = 37°, and β = 115°, etc. A comparison test between the EDEM simulation and field was conducted under the optimal parameters to verify the reliability of the optimized parameters. Table 10 shows the results of this comparative test. The simulation test results showed that the power, soil

Parameter Optimization and Comparison Test
The optimization function in Design-Expert 12.0 software was applied for the quadratic model in this study. The optimal parameters of the blade were obtained by solving the regression model with the conditions that the broken soil rate was taken as the minimum value of 65%, the minimum straw burial rate was 15%, and the power consumption was the lowest, as shown in Figure 13. The theoretical optimal parameters were obtained as the combination of radius (R) 209.761 mm, end face height (H) 45.8743 mm, pinch angle (α) 37.607°, and bending angle (β) 113.209°. In accordance with the actual fabrication process and practice, the values were taken as R = 210 mm, H = 45 mm, α = 37°, and β = 115°, etc. A comparison test between the EDEM simulation and field was conducted under the optimal parameters to verify the reliability of the optimized parameters. Table 10 shows the results of this comparative test. The simulation test results showed that the power, soil A comparison test between the EDEM simulation and field was conducted under the optimal parameters to verify the reliability of the optimized parameters. Table 10 shows the results of this comparative test. The simulation test results showed that the power, soil breakage rate, and straw burial rate were 1.62 kW, 67.47%, and 17.21%, respectively; the field results showed that the power, soil breakage rate, and straw burial rate were 1.73 kW, 71.34% and 18.89%, respectively. The average error rates of the simulated and field test values were 6.36%, 5.42%, and 8.89%, respectively, which indicated that the model was a decent fit.

Discussion
(a) Rotary blade radius R has the most obvious influence on the power consumption of single blade operation, while other indexes have no significant influence. If deep tillage is to be realized with small radius blades, the whole blade roll must be trapped in the inner part of the tillage layer in structural design, causing congestion and increasing power consumption; Conversely, using a blade with a large radius results in an increase in torque and a sharp increase in power consumption due to blade wear. The rotary tillage radius used in this paper is 210 mm, which can meet the purpose of 200 mm tillage depth and reduce power consumption at the same time. The reason for this is that a good mechanical structure and excellent operating performance make it stand out. Half of the rotary tiller rollers of this length fall into the tillage layer, which not only causes a sudden increase in power consumption caused by congestion but also achieves a good soil-cutting effect and reduces blade life loss; (b) The H had the most obvious effect on the straw cover rate. The reason is that they are close to the end of the blade and come into contact with the straw earlier, which has the most obvious impact on spraying and burying of straw. One of the objectives of this paper is to increase the straw cover rate and depth in order to better meet the agronomic requirements, reduce the straw residue on the surface of cultivated land, and bury it as far as possible to a deeper position in the tillage layer, so as to prevent the phenomenon of shelving and airing, and provide a good seeding bed for seeding. Considering comprehensively, 45 mm can meet the agronomic demand: (c) The α-influenced trend of each operation is different, so we can know by analyzing it.
In the range of 35-55 with the increase of value, power consumption P first decreases and then increases and reaches the minimum value at 36. The reason for these changes is α. The size has different effects on the shape of the end of the blade. Too small a value will cause the forward cutting edge to turn outwards and the angle of penetration to become larger, resulting in flat cutting of the blade, which cannot effectively exert the soil cutting effect of the positive cutting edge of the blade. On the contrary, excessive value results in the form of a hook down on the positive cutting edge, which is easy to cause excavation of the positive cutting edge when cutting soil, and the soil carrying effect is improved, which is not conducive to the crushing of soil blocks; (d) Bending angle β affected trend of each index is basically the same; only the RB is greatly affected. Based on the research results, 115 degrees is selected as the optimal solution. The change of its value results in the change of bending degree of the forward cutting edge. The larger the value, the closer the blade approaches the straight blade, the higher the cutting performance, and the lower the throwing effect. The smaller the value, the opposite.
It is the first time, to our knowledge, that the effect of rototill blades on soil and straw after deep rototilling has been analyzed using tillage soil stratification to describe the effect of rototill blades on soil and straw in a rice-wheat rotation. This study reveals the interaction between the rototill blade and the soil and straw, a finding that will assist in the development of technologies for deep rototill straw burial operations. In addition, the design of deep rototill blades with improved blade-breaking energy and increased straw mulching rates could lead to the provision of good seedbeds.

Conclusions
In this study, a discrete meta-simulation model of rotary tillage knife-soil-straw interaction was developed. The interaction law of soil-touching components in deep tillage operation was explored, and the field experiment results also verified its usability. The optimized structural parameters of the rotary tillage blade were 210 mm-45 mm-37 • -115 • (R-H-α-β). The simulated and field comparison tests were conducted under the optimal combination of parameters, and the power, soil breaking rate, and straw burial rate were 1.73 kW, 71.34%, and 18.89%, respectively; the numerical error rates of simulated and field test values were 6.36%, 5.42%, and 8.89%, respectively. The accuracy of the secondary model was verified.
Note that the aim of the study was to investigate the interaction of the rotary single blade with the soil straw in deep rototilling operations, applicable to optimizing the structural parameters of the rotary blade, and our experiments found that the soil movement was more consistent with the simulation results, which confirmed the validity of the simulation model. Papers exploring the arrangement of the rototill knives and the operating parameters have been published in the literature and will be the focus of our future work. Deep rototill mulching of straw will also be a hot research issue in the tillage chain. Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available within the article.