Design and Test of Soil–Fertilizer Collision Mixing and Mulching Device for Manure Deep Application Machine

: Aiming at the problems of low uniformity and utilization rate in the traditional deep application method of orchard manure, a soil–fertilizer collision mixing and mulching device was designed. The ditching mechanism, the structure of the soil divider, and the fertilizer delivery auger were analyzed and designed. Furthermore, with the rotational speed of the cutter and auger and the deﬂection angle of the soil divider as factors, as well as the uniformity of soil–fertilizer mixing and mulching as evaluation indexes, the discrete element simulation tests were conducted. The simulation results showed that when the turret speed, the stirrer speed and the soil separator deﬂection angle were 140 r/min, 146 r/min and 22 ◦ , the mixing uniformity and mulching uniformity were the highest, which was 88.35% and 96.86%, respectively. Based on the optimal parameters, the ﬁeld test was conducted, and the soil–fertilizer mixing uniformity was 87.02%, with a relative error of 1.33% compared with the simulation test results. The relative error of 94.37% of mulch uniformity is 2.49%, which indicates that the simulation optimization results are reliable and the mixing performance of the device is good and can meet the requirements of soil–fertilizer mixing operation. The results of this study can provide an important reference for the design of the soil and fertilizer mixing machine


Introduction
In the process of forest fruit cultivation, manure can improve the soil environment and soil fertility, and its reasonable application can increase the yield and improve the quality of forest fruits [1,2].In recent years, in order to promote sustainable land development and improve the quality of forest fruit production, relevant state departments have been issuing relevant documents and increasing the proportion of organic fertilizers in agricultural production, such as the "Agricultural Pollution Prevention Policy" formulated by the Ministry of Agriculture, which clearly requires reducing the use of chemical fertilizers and increasing the resource utilization of livestock and poultry manure [3][4][5].
At present, the manure applied by deep manure is larger, and the manure is not easily absorbed.As a result, most of the existing ditching fertilizer deep application machine mulching does not use soil-fertilizer mixing, making it difficult to play the effect of manure to regulate soil pH and improve the soil environment.With the increasing emphasis on green agriculture, the policy of "organic fertilizer instead of chemical fertilizer" and the agronomic requirements for the application of forest and fruit stall fertilizer are clear, and the application of manure in orchards should be fully mixed with soil before application, which can further improve the efficiency of fertilizer utilization, promote the effective absorption of fruit tree roots, and ensure the effect of fertilizer application and fertilizer utilization efficiency [6,7].
At present, scholars at home and abroad have conducted some related studies on soil and fertilizer mixing.Kwapinska et al. [8] used the discrete unit method to simulate the motion inside a horizontal rotating drum mixer, analyzed the mixing characteristics of the mixture at different rotational speeds, and explored the effects of rotational speed and filling rate on mixing time and mixing number.Ucgul et al. [9] studied a rotating shovel that was able to mix fertilizer from the surface into the soil and improve soil properties.Xiao Hongru et al. [10] designed a 1KS60-35X orchard-type double-spiral trencher and fertilizer spreader, which used double-spiral ditching and simultaneously mixed fertilizer thrown from the middle hollow shaft with soil mixing.The fertilizer application performance was stable in each soil layer and meets the agronomic requirements of orchard fertilization.Huang Yue [11] designed an orchard ditching and fertilizer mixing backfill device, which finished soil-fertilizer mixing in the trench after ditching and applicating and finally mulched the soil.The discrete element model of fertilizer mixing was established by EDEM (particle mechanics simulation software based on discrete element method), and the effect of blades with different curvature radii on the soil-fertilizer mixing uniformity was analyzed.The results showed that it achieved the technical requirement of efficient utilization of soilfertilizer mixing.Pingping Zhang et al. [12] developed a fruit tree soil-fertilizer mixing precision control fertilizer applicator, where the discharged fertilizer was pushed by a spiral auger into the soil collection box, then mixed with soil evenly, and finally sent by the auger into the ditch to be mulched by the mulching plate.It has high operational efficiency and reduces labor intensity and production cost.Yuan Quanchun et al. [13] designed and studied a soil-fertilizer mixing and layering backfill device.The soil was thrown to the soilfertilizer mixing and layered backfill device by the ditching device, after which the organic fertilizer was discharged in three ways by the fertilizer discharge mechanism, the backfill fertilizer was conveyed by the auger, and the paddles were spirally arranged between the auger to improve the soil-fertilizer mixing quality, and the soil-fertilizer mixture was backfilled to the fertilizer application trench in turn to realize layered backfill.At present, in the deep application of orchard fertilizer, the soil and fertilizer mixing methods are mainly mixed directly in the ditch and then backfilled to the ditch after being mixed, both of which are mixing soil and fertilizer by stirring, and the operation efficiency is not high.
In this paper, combining agronomy and farmers' needs, a soil-fertilizer collision mixing and mulching device, mixing soil and fertilizer while mulching in ditching operation, was designed.By designing the ditching mechanism, fertilizer discharge mechanism and soil dividing plate, establishing a discrete element simulation model, the best combination of device parameters was established, and field trials were conducted for verification.The hope is that we can realize soil-fertilizer mixing, improve the fertilizer utilization rate, and provide technical support for the fine and deep application of manure.

Overall Structure
The overall structure of the soil and fertilizer collision mixing and mulching device is shown in Figure 1, which is mainly composed of the ditching mechanism, flow-guided mechanism, soil dividing plate and fertilizer discharging mechanism.The ditching mechanism is mainly composed of a rotary tillage knife and knife plate, and the fertilizer discharging mechanism is mainly composed of auger, a fertilizer box and a scraper.

Working Process
When working, as shown in Figure 2, the ditching mechanism cut soil for ditching with device advancing, and the cut soil moves to the rear under the action of ditching knife and flow-guided mechanism.Then, the back-thrown soil was shunted and delaminated to form soil mixed with fertilizer and soil mulched with fertilizer.The soil mixed with fertilizer, in the back-throwing process, collides and mixes with the fertilizer conveyed by the auger to form soil-fertilizer mixture and then mulch to the trench.Finally, the soil mulched with fertilizer is back-thrown to cover the soil-fertilizer mixture, and the mixing and mulching work in the application of orchard manure is completed.

Main Technical Parameters
Combined with the agronomic requirements of autumn manure deep application in orchards, and to achieve soil and manure mixing and mulching, the main technical parameters of the manure deep application machine were determined, as shown in Table 1.

Working Process
When working, as shown in Figure 2, the ditching mechanism cut soil for ditching with device advancing, and the cut soil moves to the rear under the action of ditching knife and flow-guided mechanism.Then, the back-thrown soil was shunted and delaminated to form soil mixed with fertilizer and soil mulched with fertilizer.The soil mixed with fertilizer, in the back-throwing process, collides and mixes with the fertilizer conveyed by the auger to form soil-fertilizer mixture and then mulch to the trench.Finally, the soil mulched with fertilizer is back-thrown to cover the soil-fertilizer mixture, and the mixing and mulching work in the application of orchard manure is completed.

Working Process
When working, as shown in Figure 2, the ditching mechanism cut soil for ditching with device advancing, and the cut soil moves to the rear under the action of ditching knife and flow-guided mechanism.Then, the back-thrown soil was shunted and delaminated to form soil mixed with fertilizer and soil mulched with fertilizer.The soil mixed with fertilizer, in the back-throwing process, collides and mixes with the fertilizer conveyed by the auger to form soil-fertilizer mixture and then mulch to the trench.Finally, the soil mulched with fertilizer is back-thrown to cover the soil-fertilizer mixture, and the mixing and mulching work in the application of orchard manure is completed.

Main Technical Parameters
Combined with the agronomic requirements of autumn manure deep application in orchards, and to achieve soil and manure mixing and mulching, the main technical parameters of the manure deep application machine were determined, as shown in Table 1.

Main Technical Parameters
Combined with the agronomic requirements of autumn manure deep application in orchards, and to achieve soil and manure mixing and mulching, the main technical parameters of the manure deep application machine were determined, as shown in Table 1.The ditching mechanism cuts the soil and completes the furrowing operation while throwing the cut soil back to the inflow mulching component, including ditching the disc, furrowing the knife and driving the mechanism of the cutter disc.The disc-type trencher breaks soil evenly and works efficiently, and it is widely used in various orchards for ditching and fertilization and other agricultural production fields [13].During the ditching operation, the motion of the ditching blade is synthesized by the horizontal motion of the trencher and the rotational motion of the disc, and the motion trajectory is a coswing [14][15][16], as shown in Figure 3.The ditching mechanism cuts the soil and completes the furrowing operation while throwing the cut soil back to the inflow mulching component, including ditching the disc, furrowing the knife and driving the mechanism of the cutter disc.The disc-type trencher breaks soil evenly and works efficiently, and it is widely used in various orchards for ditching and fertilization and other agricultural production fields [13].During the ditching operation, the motion of the ditching blade is synthesized by the horizontal motion of the trencher and the rotational motion of the disc, and the motion trajectory is a coswing [14][15][16], as shown in Figure 3. Let the coordinates of any point on the trajectory of the end point of the ditching blade be F(x, y).Take the center of rotation of the ditching blade O as the coordinate origin, the forward direction of the machine is the positive direction of the x-axis, and the vertical downward direction of the y-axis, as shown in the figure.The equation of the motion trajectory of the point The tool endpoint at any moment x-axis and y-axis partial velocity is Let the coordinates of any point on the trajectory of the end point of the ditching blade be F(x, y).Take the center of rotation of the ditching blade O as the coordinate origin, the forward direction of the machine is the positive direction of the x-axis, and the vertical downward direction of the y-axis, as shown in the figure.The equation of the motion trajectory of the point F(x, y) is The tool endpoint at any moment x-axis and y-axis partial velocity is Then, the absolute velocity v at the endpoint of the tool is As shown in Equation ( 4), the higher the rotational speed of the cutter, the greater the linear speed of the cutter edge.However, the increase in cutting speed will also cause the rise in power consumption of soil cutting and throwing.When the cutter speed is low, the soil cannot be thrown out, but it can only be passively pushed upward by the soil below.It may damage the cutter.Combined with the pre-experimental result, the final setting of the cutter speed is 140-160 r/min.
Cutter diameter D is one of its important structural parameters, which has a certain influence on the soil spreading distance of the rotary ditching device, power consumption of the machine and the size of the model, etc.According to the formula where H-depth of ditching, m.The ditching depth designed in this paper is 30 cm, which is finally taken as D = 40 cm.

Fertilizer Delivery Auger
In this paper, a spiral auger was selected to convey manure, and the auger parameters were determined according to the operational requirements [17].According to the agronomic requirements, the auger fertilizer delivery rate was determined to be 10 t/h, and the auger outer diameter and delivery volume equations are as follows: where D-spiral diameter, m; Q-conveying volume per hour, t/h; K-material characteristic coefficient; ϕ-filling coefficient; ρ-material accumulation density, t/m 3 ; c-inclination coefficient; S-pitch; n-rotational speed, r/min.
The material characteristic coefficient K is 0.05, the filling coefficient is 1, the auger operation is horizontal, so the inclination angle coefficient C is 1, and the screw diameter D is 0.25 m.The pitch S is 0.20 m.According to the auger conveying capacity shown in Formula ( 7), the minimum winch speed n is 120 r/min to reach the theoretical conveying capacity.The spiral conveying component assembly is shown in Figure 4.During the operation of the spiral auger, with the propulsion of the spiral blades, the fertilizer moves along the circumferential direction of the spiral blades while making horizontal movements along the axial direction [18,19].The velocity analysis of a fertilizer particle on the spiral blade is shown in Figure 5, where the location of the fertilizer particle During the operation of the spiral auger, with the propulsion of the spiral blades, the fertilizer moves along the circumferential direction of the spiral blades while making horizontal movements along the axial direction [18,19].The velocity analysis of a fertilizer particle on the spiral blade is shown in Figure 5, where the location of the fertilizer particle is at a distance from the axis and the spiral lift angle is α.When the spiral blade rotates, the implication velocity of the fertilizer particle is equal to the linear velocity V N at that point, and the implication velocity direction is perpendicular to the wall spread line of the spiral blade.During the operation of the spiral auger, with the propulsion of the spiral blades, the fertilizer moves along the circumferential direction of the spiral blades while making horizontal movements along the axial direction [18,19].The velocity analysis of a fertilizer particle on the spiral blade is shown in Figure 5, where the location of the fertilizer particle is at a distance from the axis and the spiral lift angle is α.When the spiral blade rotates, the implication velocity of the fertilizer particle is equal to the linear velocity VN at that point, and the implication velocity direction is perpendicular to the wall spread line of the spiral blade.As can be seen from Figure 5, the absolute velocity Va of fertilizer in the conveying process can be decomposed into the circumferential directional partial velocity Vr and the axial directional partial velocity Vz, whose values are ( ) ( ) As can be seen from Figure 5, the absolute velocity V a of fertilizer in the conveying process can be decomposed into the circumferential directional partial velocity V r and the axial directional partial velocity V z , whose values are where: n-spiral auger speed, r/min.According to Equations ( 10)- (12).
In the working process of soil-fertilizer mixing and mulching device, the speed of the auger shaft directly affects the mixing effect.If the speed is larger, the spiral blade produces larger centrifugal force along the circumference, and the manure particles will be thrown to the edge of the spiral blade.While the material axial speed is faster, the material is quickly discharged.It reduces the mixing effect of soil-fertilizer.If the speed is smaller, the manure will cake due to the increasing squeezing pressure between fertilizer and its high water content.It is not conducive to the mixing.Therefore, based on the theoretical analysis and related research, the auger speed range was set from 130 to 150 r/min according to the pre-experimental result of the device.

Soil Separation Plate
The spreading of soil during ditching is often chaotic and irregular, but most soil particle groups will have the same tendency as well as fall within a certain spreading range.So, it is crucial to design a soil-splitting deflector to analyze the landing point of soil particle groups.
The splitting plate not only diverts the back-thrown soil to form the soil mixed with fertilizer and the soil mulched with fertilizer but also controls the amount of both types of soil.Its operation effect directly affects the uniformity of fertilizer mixing and mulching [20].In this paper, the EDEM simulation software and theoretical calculations were applied to obtain the optimal structural parameters of the soil splitting plate by constructing a model of the soil splitting plate [21,22].The force analysis of a soil particle on the splitting plate was analyzed as shown in Figure 6.
Agriculture 2023, 13, x FOR PEER REVIEW 8 of 20 [20].In this paper, the EDEM simulation software and theoretical calculations were applied to obtain the optimal structural parameters of the soil splitting plate by constructing a model of the soil splitting plate [21,22].The force analysis of a soil particle on the splitting plate was analyzed as shown in Figure 6.The combined external force on the soil particles at point D is F.
From Newton's second law.The combined external force on the soil particles at point D is F. Then From Newton's second law.

F = ma
And The derivation of Equation ( 20) yields the absolute velocity V a of the soil particles on the inflow surface, with components in the x and y axes as Different deflection angles of the splitter plates will affect the ratio of the two parts of the soil and the movement of the soil mixed with fertilizer.If the deflection angle of the splitter and the distance between the splitter and the deflector is too small, it will cause the soil to accumulate on top of the splitter and prevent the fertilizer soil from being discharged properly.Based on the above considerations and pre-experiments, the declination angle of the soil divider was finally set at 20~40 • .

EDEM Simulation Model and Simulation Parameters
The model of the soil-fertilizer collision mixing mulching device was created in SolidWorks software and imported into EDEM software.A soil layer of 3000 mm in length, 600 mm in width and 400 mm in height was created, and the soil color was set to red, while the manure was to yellow for better differentiation.The Hertz-mindlin with the JKR model was used between the soil and the manure, and the Hertz-mindlin (no slip) built-in Optimal model was used between the fertilizer particles and the fertilizer discharge device [23].According to scaling theory, the simulation with a large particle size did not affect the results after modifying the discrete element model accordingly [24].Due to the large number of particles, the particle radius was scaled to 5 mm in this paper to shorten the simulation time, and the soil and manure discrete element parameters were referenced as shown in Table 2 [25,26].
According to the simulation pre-experiment, the time step during simulation is 2.5 × 10 −5 s, the data-saving interval is 0.01 s, and the simulation grid is 20 times the particle radius [27].The soil-fertilizer collision mixing mulching device moves in the forward direction with a speed of 0.3 m/s, and the operation time is 7 s.

EDEM Soil Fertilizer Mixing Process Analysis
As shown in Figure 7, during the operation, the ditching and mulching device will produce cutting, throwing and the collision effects on the soil and manure, which affect the soil and fertilizer mixing operation effect.To explain the soil and manure transport mechanism under different stages of the machine operation, the kinetic energy and motion trend analysis of the simulated soil and manure particles are carried out with the EDEM post-processing.In order to observe the flow direction of soil particles, the shape of soil particles is expressed as an arrow shape, and the direction of the arrow indicates the direction of particle movement.The operating characteristics of the machine at a travel speed of 0.3 m/s and a cutter speed of 170 r/min were used as an example.
As shown in Figure 8, when the soil groove was generated, the rotating tillage knife was not in contact with the soil particles, and the soil particles were not subject to the crushing and squeezing effect of the ditching knife on the soil.At this time, the speed of the soil particles was 0; there was no displacement trend.With the compound movement of the ditching knife driven by the knife shaft, the blade-side cutting edge began to contact the soil, and the blade mainly through shear extrusion cut down the soil.At this time, the soil particles were thrown up by the ditching knife.With the ditching knife cutting soil volume increasing, the strength of the soil crushing increased in the ditching knife pushing action.The cut down soil has a tendency along the upward movement, and the soil In order to observe the flow direction of soil particles, the shape of soil particles is expressed as an arrow shape, and the direction of the arrow indicates the direction of particle movement.The operating characteristics of the machine at a travel speed of 0.3 m/s and a cutter speed of 170 r/min were used as an example.
As shown in Figure 8, when the soil groove was generated, the rotating tillage knife was not in contact with the soil particles, and the soil particles were not subject to the crushing and squeezing effect of the ditching knife on the soil.At this time, the speed of the soil particles was 0; there was no displacement trend.With the compound movement of the ditching knife driven by the knife shaft, the blade-side cutting edge began to contact the soil, and the blade mainly through shear extrusion cut down the soil.At this time, the soil particles were thrown up by the ditching knife.With the ditching knife cutting soil volume increasing, the strength of the soil crushing increased in the ditching knife pushing action.The cut down soil has a tendency along the upward movement, and the soil kinetic energy gradually increases [28].
and a cutter speed of 170 r/min were used as an example.
As shown in Figure 8, when the soil groove was generated, the rotating tillage knife was not in contact with the soil particles, and the soil particles were not subject to the crushing and squeezing effect of the ditching knife on the soil.At this time, the speed of the soil particles was 0; there was no displacement trend.With the compound movement of the ditching knife driven by the knife shaft, the blade-side cutting edge began to contact the soil, and the blade mainly through shear extrusion cut down the soil.At this time, the soil particles were thrown up by the ditching knife.With the ditching knife cutting soil volume increasing, the strength of the soil crushing increased in the ditching knife pushing action.The cut down soil has a tendency along the upward movement, and the soil kinetic energy gradually increases [28].As shown in Figure 9, the spreading of soil in the ditching process is often chaotic and irregular, but most of the soil particles will have the same trend and fall within a certain spreading range [29].The deflector cover can deflect the soil in a predetermined direction when the cutter is spreading soil.The cut soil is thrown backward along the inflow surface under the action of the high-speed rotating ditching knife disk axis.When As shown in Figure 9, the spreading of soil in the ditching process is often chaotic and irregular, but most of the soil particles will have the same trend and fall within a certain spreading range [29].The deflector cover can deflect the soil in a predetermined direction when the cutter is spreading soil.The cut soil is thrown backward along the inflow surface under the action of the high-speed rotating ditching knife disk axis.When the kinetic energy is increasing, the particles are thrown to the splitting plate, and they are no longer subject to the action of the knife disk.Due to the inertia and internal force, the soil continues to maintain the movement.Under the action of the inflow plate, the backward thrown soil is divided into mixed fertilizer soil and mulch soil.When the soil collides with the manure, the collision time is short, so the time of the collision process and the displacement occurred can be neglected, and only the kinetic energy change before and after the collision was considered.As shown in Figure 10, the fertilizer soil mixture collides with the manure discharged by the auger.The kinetic energy is instantaneously reduced, and the collision changes the trajectory of movement for some soil particles, thus changing the landing point of the soil-fertilizer mixture.The soilfertilizer mixture falls exactly to the ditch opened by the ditching device, and the soil When the soil collides with the manure, the collision time is short, so the time of the collision process and the displacement occurred can be neglected, and only the kinetic energy change before and after the collision was considered.As shown in Figure 10, the fertilizer soil mixture collides with the manure discharged by the auger.The kinetic energy is instantaneously reduced, and the collision changes the trajectory of movement for some soil particles, thus changing the landing point of the soil-fertilizer mixture.The soil-fertilizer mixture falls exactly to the ditch opened by the ditching device, and the soil mulched with fertilizer continues to do a backward throwing movement downward, covering the soil-fertilizer mixture, thus completing the ditching, fertilizer mixing and mulching operation.When the soil collides with the manure, the collision time is short, so the time of the collision process and the displacement occurred can be neglected, and only the kinetic energy change before and after the collision was considered.As shown in Figure 10, the fertilizer soil mixture collides with the manure discharged by the auger.The kinetic energy is instantaneously reduced, and the collision changes the trajectory of movement for some soil particles, thus changing the landing point of the soil-fertilizer mixture.The soilfertilizer mixture falls exactly to the ditch opened by the ditching device, and the soil mulched with fertilizer continues to do a backward throwing movement downward, covering the soil-fertilizer mixture, thus completing the ditching, fertilizer mixing and mulching operation.

Response Surface Experimental Design
Combined with the above analysis, the Box-Behnken test method was used.The test factors were blade speed, stirrer speed, and soil separating plate deflection angle.The factor levels were set based on the pre-experimental results.Three levels were set for each

Response Surface Experimental Design
Combined with the above analysis, the Box-Behnken test method was used.The test factors were blade speed, stirrer speed, and soil separating plate deflection angle.The factor levels were set based on the pre-experimental results.Three levels were set for each factor, and the test factors and levels are shown in Table 3.The soil-fertilizer collision mixing mulching device was designed to finally improve the fertilizer utilization rate, promote the effective absorption of fruit tree roots and solve the poor mulching quality, so the test indexes were mulching uniformity and soil-fertilizer mixing uniformity [30].In the study of particle mixing, multiple samples are often collected, and the number of target particles in each sample is counted.Then, the coefficient of variation is calculated to characterize the mixing uniformity: the lower the coefficient of variation, the higher the mixing uniformity [31].In this paper, we use the sampling method to calculate the soil fertilizer mixing uniformity.Firstly, we divide the particles in the fertilizer application trench into several samples and count the amount of manure in each sample.Then, we eliminate the samples with the total number of particles less than 20.Finally, we take the average of the results 10 times and calculate the homogeneity of soil-fertilizer mixing.The size of the samples was set as 60 mm × 60 mm × 60 mm cubes, with 20 samples in total.
Mixing uniformity is calculated with reference to the standard "NY/T 1024-2006 Technical Specification for Quality Evaluation of Feed Mixers" for the mixing uniformity M of stable manure and granular fertilizer.
where S-sample standard deviation.X-sample X 1 , X 2 for each measurement . . .mean of X n , kg.X i -mass of tracer in the sample, kg.n-number of samples.

Response Surface Experimental Results Analysis
A Box-Behnken response surface test was carried out, and the test scheme and results are shown in Table 4. ANOVA and significance ANOVA and significance tests were performed on the quadratic regression model, and the results are shown in Table 5.  16 16 Note: p ≤ 0.01 is highly significant, marked as **; 0.01 < p ≤ 0.05 is significant, marked as *; p > 0.05 is not significant, marked as -.
From the analysis of variance in Table 5, it can be seen that the regression model is extremely significant while loss of fit is not significant, indicating that the fitted regression equation of the model constructed for mulch uniformity is consistent with reality, and it can provide a better prediction of the relationship between mulch uniformity and each of the tested factors [32].From the p-values in the table, it can be seen that the interaction of cutter speed and auger speed, the square of the cutter speed, and the square of the splitting plate deflection angle have a highly significant effect on the uniformity of mulching; the splitting plate deflection angle as well as the interaction of the auger speed and splitting plate deflection angle have a significant effect on the uniformity of mulching; and the interaction of cutter speed and the splitting plate deflection angle as well as the square of the auger speed have no significant effect on the mulching uniformity.Combined with the regression equation of uniformity of mulching, the analysis shows that the factors affecting the uniformity of mulching are mainly the square of the splitting plate deflection angle term X 3 2 .Therefore, based on the ANOVA of each factor on the mulching uniformity, after excluding the non-significant, the relationship was Design-Expert 8.0.6 software was used to analyze the influences of turret speed, stirrer speed, and soil dividing plate deflection angle on mulching uniformity and mixing uniformity as well as the relationships among various factors, and to draw response surfaces, as shown in Figures 11 and 12     Figure 11a shows that the mulching uniformity increases with the increase in the a ger speed when the deflection angle of the plate is at the zero level, and the maximu value is achieved at 150 r/min of the auger speed; the mulching uniformity increases w the increase in the auger speed during the change of the blade speed from 130 to 160 r/m    Figure 11a shows that the mulching uniformity increases with the increase in the auger speed when the deflection angle of the plate is at the zero level, and the maximum value is achieved at 150 r/min of the auger speed; the mulching uniformity increases with the increase in the auger speed during the change of the blade speed from 130 to 160 r/min; Figure 11a shows that the mulching uniformity increases with the increase in the auger speed when the deflection angle of the plate is at the zero level, and the maximum value is achieved at 150 r/min of the auger speed; the mulching uniformity increases with the increase in the auger speed during the change of the blade speed from 130 to 160 r/min; the back-throwing speed of the soil increases during the increase in the blade speed, so that it quickly follows the guide hood located on the plate.As the rotating speed of the cutter increases, the back-throwing speed of the soil increases so that it quickly follows the deflector hood located on the splitting plate to carry out the mulching operation, which improves the mulching uniformity.The trend diagram shows that the interaction between the blade speed and the auger speed has a more significant effect.Figure 11b shows that the mulching uniformity decreases and then increases as the deflection angle of the splitting plate increases at zero level; the larger the deflection angle of the splitting plate, the more soil is mulched and the less soil is mixed, which will reduce the soil-fertilizer mixing uniformity.With the increase in auger speed, the uniformity of mulching soil increases, and the trend of change shows that the interaction between the deflection angle of the soil divider and the auger speed has a significant effect.
From the analysis of variance in Table 6, it can be seen that the regression model is extremely significant, and the loss of fit is not significant, which indicates that the fitted regression equation of the model constructed for mixing uniformity is consistent with the reality and can make a better prediction of the relationship between mixing uniformity and each test factor.From the p-values in the table, it can be seen that the influence of the deflection angle of the soil divider and the square of it on the mixing uniformity is highly significant.The square of the deflection angle of the soil divider and the square of the stirrer rotational speed have a significant influence on the mixing uniformity.Lastly, the interaction of the stirrer rotational speed and the deflection angle as well as the square of the stirrer rotational speed do not have a significant influence on the mixing uniformity.Therefore, according to the analysis of variance of each factor on the mixing uniformity, after excluding the insignificant, the relationship is as follows As shown in Figure 12, the uniformity of fertilizer mixing increases from low level to high level with the deflection angle of the soil divider when the rotating speed of the cutter is at the zero level; the uniformity of fertilizer mixing increases with the increase in rotating speed during the change of rotating speed of the stirrer shaft from 130 to 150 r/min.The back-throwing speed of the soil increases during the process of the knife plate rotational speed becoming larger, which makes the collision mixing effect more significant.From the change trend graph, it can be seen that the interaction effect of the soil dividing plate deflection angle and auger speed is more significant.

Optimization of Fertilizer Mixing Unit
On the basis of ensuring that each working parameter of the soil collision mixing and mulching device meets the agronomic requirements and the requirements of fertilization operations, optimization tests were conducted to obtain the optimal range of working parameters of it.To seek the optimal combination of working parameters, the Optimization module in Design-Expert 13 software was used to take the function of multi-objective optimization, with fertilizer mixing and mulching uniformity as the objective functions, and the blade speed, stirrer shaft speed, and splitting plate deflection angle as variables for optimization [33].
The mulch uniformity Y 1 and fertilizer mixing uniformity Y 2 are used as the working performance indexes of soil-fertilizer mixing and the fertilizer discharge device.Among the above response indexes, the mulch uniformity and fertilizer mixing uniformity should be improved as much as possible under the condition of meeting the requirements of fertilization operation.The optimized solution for the maximum fertilizer mixing uniformity is shown in Table 7.

Field Trials
To verify the effectiveness of the soil-fertilizer collision mixing mulching device, a field trial was conducted on 22 June 2022 at a small plant of the School of Mechanical and Electrical Engineering, Shihezi University.The manure used was produced by Xinjiang Shihezi Xinye Agriculture and Animal Husbandry Biological Development Co. Ltd., (Shihezi, China) and a John Deere 454 tractor was used to tow the soil-fertilizer collision mixing mulching test bed.The measurement tools used in the test were a soil moisture meter (TDR150), tape measure (range: 30 m, accuracy: 2 mm), steel ruler (range: 500 mm, index value: 1 mm), and inverter (AMB100-7R5G-T3 produced by Shenzhen Huichuan Technology Co. LTD, Foochow, China).The test site diagram is shown in Figure 13.
Agriculture 2023, 13, x FOR PEER REVIEW 16 of 20 be improved as much as possible under the condition of meeting the requirements of fertilization operation.The optimized solution for the maximum fertilizer mixing uniformity is shown in Table 7.

Field Trials
To verify the effectiveness of the soil-fertilizer collision mixing mulching device, a field trial was conducted on 22 June 2022 at a small plant of the School of Mechanical and Electrical Engineering, Shihezi University.The manure used was produced by Xinjiang Shihezi Xinye Agriculture and Animal Husbandry Biological Development Co. Ltd., (Shihezi, China) and a John Deere 454 tractor was used to tow the soil-fertilizer collision mixing mulching test bed.The measurement tools used in the test were a soil moisture meter (TDR150), tape measure (range: 30 m, accuracy: 2 mm), steel ruler (range: 500 mm, index value: 1 mm), and inverter (AMB100-7R5G-T3 produced by Shenzhen Huichuan Technology Co. LTD, Foochow, China).The test site diagram is shown in Figure 13.In this test, because the physical properties of soil are similar to manure (both are in powder form), it is not possible to accurately calculate the mixing uniformity after mixing.Referring to the standard "NY/T 1024-2006 Technical Specification for Quality Evaluation of Feed Mixers", corn was used as tracer particles; before the soil-fertilizer mixing test, In this test, because the physical properties of soil are similar to manure (both are in powder form), it is not possible to accurately calculate the mixing uniformity after mixing.Referring to the standard "NY/T 1024-2006 Technical Specification for Quality Evaluation of Feed Mixers", corn was used as tracer particles; before the soil-fertilizer mixing test, corn has been mixed with manure.A flat and smooth concrete floor was selected; the length of the measuring area is 10 m, and the width covers one working width of the machine.The device stably goes through the measuring area to complete the fertilizer discharge operation.After the test, 50 cm front and rear of the fertilizer in the test area were removed, and the remaining part was divided into not less than 30 consecutive equal sections according to the length of 10 cm.The mass of corn in each small section was collected and measured to evaluate the uniformity of mixing of corn as tracer particles in the manure.The uniformity of corn mixing in the manure was calculated to be 98.53%.
As there was soil flow during the sampling and measurement process, the actual value could not be measured accurately.So, to ensure that the soil-fertilizer mixture results were closer to the actual value, two acrylic plates were taken and inserted into the bottom of the trench and fixed before sampling, as shown in Figure 14.The two plates were spaced 10 cm apart, and the trench on one side of the plates was cleaned up so that the soil-fertilizer mixture and the overall condition of the mulch could be seen through the transparent plates.The mulching soil between the plates was first taken out with a small shovel and weighed for the mass so as to find out the uniformity of that.Afterwards, the soil-fertilizer mixture in the middle was removed and sieved to find the mass of corn.Each group test was repeated five times to find the mean and standard deviation of the five times for the calculation of soil-fertilizer mixture uniformity [34,35].
As there was soil flow during the sampling and measurement process, the actual value could not be measured accurately.So, to ensure that the soil-fertilizer mixture results were closer to the actual value, two acrylic plates were taken and inserted into the bottom of the trench and fixed before sampling, as shown in Figure 14.The two plates were spaced 10 cm apart, and the trench on one side of the plates was cleaned up so that the soil-fertilizer mixture and the overall condition of the mulch could be seen through the transparent plates.The mulching soil between the plates was first taken out with a small shovel and weighed for the mass so as to find out the uniformity of that.Afterwards, the soil-fertilizer mixture in the middle was removed and sieved to find the mass of corn.Each group test was repeated five times to find the mean and standard deviation of the five times for the calculation of soil-fertilizer mixture uniformity [34,35].The soil-fertilizer mixing uniformity is sensitive to sample size.In order to verify the accuracy of the discrete element model and the reliability of the simulation optimization results, the sample size is set to be consistent with the bench test in the simulation test with the optimal parameters.The test results of soil cover uniformity and soil and fertilizer mixture uniformity are shown in Table 8, and the soil cover uniformity is 94.37%, while the mixing uniformity was 87.02%.The results show that the discrete element model is accurate, the simulation optimization results are reliable, and the mixing performance of the soil-fertilizer collision mixing mulching device processed according to the optimized parameters is good.The soil-fertilizer mixing uniformity is sensitive to sample size.In order to verify the accuracy of the discrete element model and the reliability of the simulation optimization results, the sample size is set to be consistent with the bench test in the simulation test with the optimal parameters.The test results of soil cover uniformity and soil and fertilizer mixture uniformity are shown in Table 8, and the soil cover uniformity is 94.37%, while the mixing uniformity was 87.02%.The results show that the discrete element model is accurate, the simulation optimization results are reliable, and the mixing performance of the soil-fertilizer collision mixing mulching device processed according to the optimized parameters is good.

Discussion
Firstly, the apple-growing areas were researched, and relevant information was collected and organized.Insisting on the principle of combining agricultural machinery and agronomy, a technical scheme of soil-fertilizer mixing and mulching device for deep fertilizer application machine was proposed.Secondly, the conveying of manure particles in the spiral auger discharge parts and of soil in the ditching knife and guide cover were analyzed.Thirdly, the key components of the device were designed, and a simulation test was conducted to analyze the mixing process using EDEM.Finally, the working parameters of the device were optimized according to the test results.

Conclusions
(1) The design of the soil-fertilizer mixing and mulching device for the deep application of manure was carried out, and the key parameters of the ditching mechanism, soildividing plate and auger dragon were determined, which can provide a reference for the research of soil-fertilizer mixing and application technology in a Xinjiang apple orchard.(2) The simulation test of soil and fertilizer mixing and mulching process was conducted by EDEM.When the rotational speed of the cutter, stirrer and deflection angle of the dividing plate were 140 r/min, 146 r/min and 22 • , respectively, the highest uniformity of mulching was 96.86%, and the largest uniformity of mixing was 88.35%.(3) In the soil-fertilizer mixing and mulching bench test, the test results show that the average fertilizer mixing uniformity is 87.02%, with an average error of 1.33% from the optimal value of the simulation test, and the uniformity of mulching is 94.37%, with an average error of 2.49%, which meets the agronomic requirements.It shows that the discrete element model is accurate, the simulation optimization results are reliable, and the soil-fertilizer collision mixing mulching device processed according to the optimized parameters has better mixing performance and can meet the operation requirements.

Figure 3 .
Figure 3. Analysis of ditching tool end point motion.R-Drenching knife disc end point gyration radius, m; t-Time, s; vm-Trencher forward speed, km/h; ω-Knife shaft angular speed, rad/s.

Figure 3 .
Figure 3. Analysis of ditching tool end point motion.R-Drenching knife disc end point gyration radius, m; t-Time, s; v m -Trencher forward speed, km/h; ω-Knife shaft angular speed, rad/s.

Agriculture 2023 , 20 Figure 4 .
Figure 4. Schematic diagram of the structure of screw conveying device.

Figure 4 .
Figure 4. Schematic diagram of the structure of screw conveying device.

Figure 4 .
Figure 4. Schematic diagram of the structure of screw conveying device.

Figure 5 .
Figure 5. Fertilizer speed decomposition diagram.S-pitch, mm; r-distance between the particles on the spiral surface and the spiral axis, mm; VN-implicated velocity, m/s; Va-absolute velocity, m/s; Vr-circumferential directional fractional velocity, m/s; Vz-axial directional fractional velocity, m/s; V0-initial velocity, m/s; α-spiral lift angle of the spiral blade, °; β-external friction angle of the fertilizer particles, °.

Figure 5 .
Figure 5. Fertilizer speed decomposition diagram.S-pitch, mm; r-distance between the particles on the spiral surface and the spiral axis, mm; V N -implicated velocity, m/s; V a -absolute velocity, m/s; Vr-circumferential directional fractional velocity, m/s; V z -axial directional fractional velocity, m/s; V 0 -initial velocity, m/s; α-spiral lift angle of the spiral blade, • ; β-external friction angle of the fertilizer particles, • .

Figure 6 .
Figure 6.Schematic diagram of soil separation plate structure.N1-the force of diversion surface on soil particles, N; f1-the friction force between the diversion surface and soil particles, N; F1-the normal force of trench cutter on soil particles, N; G-the gravity of soil particles, N.

Figure 6 .
Figure 6.Schematic diagram of soil separation plate structure.N 1 -the force of diversion surface on soil particles, N; f 1 -the friction force between the diversion surface and soil particles, N; F 1 -the normal force of trench cutter on soil particles, N; G-the gravity of soil particles, N.

Figure 8 .
Figure 8. Soil cutting process.(a) the rotating tillage knife was not in contact with the soil.(b) the blade-side cutting edge began to contact the soil.(c) The contact area between the ditching knife cutter and the soil becomes larger.

Figure 8 .
Figure 8. Soil cutting process.(a) the rotating tillage knife was not in contact with the soil.(b) the blade-side cutting edge began to contact the soil.(c) The contact area between the ditching knife cutter and the soil becomes larger.

Agriculture 2023 ,
13,  x FOR PEER REVIEW 11 of 20 the kinetic energy is increasing, the particles are thrown to the splitting plate, and they are no longer subject to the action of the knife disk.Due to the inertia and internal force, the soil continues to maintain the movement.Under the action of the inflow plate, the backward thrown soil is divided into mixed fertilizer soil and mulch soil.

Figure 9 .
Figure 9. Soil throwing infusion process.(a) The soil is cut and thrown into the guide cover.(b) Trajectory of soil movement in the guide cover.

Figure 9 .
Figure 9. Soil throwing infusion process.(a) The soil is cut and thrown into the guide cover.(b) Trajectory of soil movement in the guide cover.

Figure 9 .
Figure 9. Soil throwing infusion process.(a) The soil is cut and thrown into the guide cover.(b) Trajectory of soil movement in the guide cover.
Agriculture 2023, 13, x FOR PEER REVIEW 14 of uniformity as well as the relationships among various factors, and to draw response s faces, as shown in Figures 11 and 12 .

Figure 11 .
Figure 11.Effect of interaction factor on the mulching uniformity.(a) Influence of the interact terms between turret speed and stirrer speed on the mulching uniformity.(b) Influence of the in teraction terms of Stirrer speed and soil separator deflection angle on the mulching uniformity.

Figure 12 .
Figure 12.Influence of the interaction terms of Stirrer speed and soil separator deflection angle the mixing uniformity.

Figure 11 .
Figure 11.Effect of interaction factor on the mulching uniformity.(a) Influence of the interaction terms between turret speed and stirrer speed on the mulching uniformity.(b) Influence of the interaction terms of Stirrer speed and soil separator deflection angle on the mulching uniformity. .

Figure 11 .
Figure 11.Effect of interaction factor on the mulching uniformity.(a) Influence of the interaction terms between turret speed and stirrer speed on the mulching uniformity.(b) Influence of the interaction terms of Stirrer speed and soil separator deflection angle on the mulching uniformity.

Figure 12 .
Figure 12.Influence of the interaction terms of Stirrer speed and soil separator deflection angle on the mixing uniformity.

Figure 12 .
Figure 12.Influence of the interaction terms of Stirrer speed and soil separator deflection angle on the mixing uniformity.

Figure 13 .
Figure 13.Soil and fertilizer mixing whole machine test.

Figure 13 .
Figure 13.Soil and fertilizer mixing whole machine test.

Table 1 .
Main technical parameters of the whole machine.

Table 1 .
Main technical parameters of the whole machine.

Table 3 .
Table of factor levels of the interaction test.

Table 4 .
Response surface test scheme and results.

Table 5 .
Analysis variance for the cover uniformity back model.

Table 6 .
Regression model analysis of variance for soil fertilizer mixing uniformity.

Table 8 .
Optimization results test values.

Table 8 .
Optimization results test values.