Optimal Design and Experiment of Corn-Overlapped Strip Fertilizer Spreader

: As the complex terrain in hilly areas is not conducive to corn mid-tillage precision fertilization, a corn-overlapped strip fertilizer spreader was designed without an external power source. By conﬁguring a passive overlapping spreading method with a three-branch split chamber structure, the uniform spreading of fertilizer in strips was achieved. A horizontal and vertical movement model of fertilizer spreading was developed to determine the angle of the fertilizer extending tube, the width of fattening small plates, and the height of the fertilizer spread as the main factors affecting the fertilizer distribution pattern. The single-factor ternary orthogonal rotational combination response surface simulation test was carried out with pendulum angle, width, and height as test factors and the transversal fertilizer uniformity coefﬁcient and longitudinal fertilizer uniformity coefﬁcient as test indicators. The test results showed that the pendulum angle, height, and width had signiﬁcant effects ( p < 0.05) on the transversal fertilizer uniformity coefﬁcient, and the pendulum angle and width had a considerable impact ( p < 0.05) on the longitudinal fertilizer uniformity coefﬁcient. In the optimal combination of parameters, swing angle 52 ◦ , height 400 mm, and width 50 mm operation, the coefﬁcients of uniformity of both the transversal fertilizer uniformity coefﬁcient and longitudinal fertilizer uniformity coefﬁcient were less than 0.15%. A veriﬁcation test was carried out under the optimal combination of parameters for the simulation tests with the simulation conditions as the standard. The test results were consistent with the simulation results within the error range. The deviation values of the transversal fertilizer uniformity coefﬁcient and longitudinal fertilizer uniformity coefﬁcient were 8.11% and 9.01%, respectively. The corn-overlapped strip fertilizer spreader was able to complete the fertilizer spreading operation smoothly. This study provides evidence for further optimizing the performance of the corn mid-tillage fertilizer applicator.


Introduction
China's geographic environment, where hilly mountainous areas account for 43% of the complex landscape, restricts the development of agricultural mechanization [1].Excellent corn fertilizer chasing machinery can not operate in the field, and due to the lack of a small fertilizer spreading device, fertilization is still dominated by manual spreading [2].Uneven spreading of fertilizer can also lead to low fertilizer utilization, reduce the quality yield of maize, and cause problems such as the environmental pollution of soil slab [3,4].Therefore, it is necessary to study the corn fertilizer fertilization device in hilly areas to enhance the spreading efficiency, increase the emergence rate, and improve the soil environment of the seed bed.
The mechanical fertilizer spreading method is mainly divided into strip application and spreading; the inline gravity type belongs to strip application, applying fertilizer Appl.Sci.2023, 13, 2559 2 of 18 without a power source and low fertilizer utilization, and the centrifugal spreading type belongs to spreading, and fertilizer spreading needs a power source [5,6].In recent years, the discrete element method has been used to simulate and analyze the interaction between agricultural bulk materials and mechanical equipment to good effect [7], providing a new means of digital design for modern agricultural equipment.Especially in the application of fertilizer dischargers, fertilizer application devices and fertilizer spreading devices are more common.Researchers have conducted many studies on fertilizer spreading devices.Hongxin Liu [8] et al. designed an auxiliary roller for the side-throwing of organic fertilizer on the opposite swashplate and studied the influence of the roller speed, spiral angle, and blade number on the uniformity variation coefficient using the discrete element method; Zhengdao Liu [9] et al. designed a pneumatic seed fertilizer cavity application device and used the coupled CFD-DEM method to clarify the flow field distribution within the fertilizer delivery mechanism and its effect on the movement of fertilizer masses under different delivery paths, which ultimately improved the fertilizer yield; Cancan Song [10] et al. designed fertilizer discharge wheels with different groove shapes and number of groove columns, and used EDEM simulation and bench tests to test the discharge range of each wheel, as well as the pulsation and accuracy when discharging fertilizer so as to preferably select a fertilizer discharge wheel that meets the requirements of UAV fertilizer application; Xiaodong Liu [11] et al. designed a spiral cone centrifugal fertilizer application device, based on a discrete element approach to optimize the curved conical discs that affect the uniformity of the fertilizer application device; Shangpeng Ding [12] et al. designed a dual-frequency fertilizer applicator and used the discrete element method (DEM) to model the operating process of the applicator to examine the effect of machine parameters on the ratio of starting fertilizer to base fertilizer discharged and the separation distance of the fertilizer band.The above-mentioned studies on discrete elements in agricultural machinery all provide good technical tools and methods for this study.
Overall, the existing research provides further research for ample fertilizer spreading devices.However, fewer small fertilizer-spreading machines have fully considered nonpowered drives.For the complex terrain of hilly mountainous areas, choosing a minor, lightweight, low-power fertilizer spreading device is more suitable [13,14].Therefore, based on the gravity inline fertilizer spreading technology, the development of new fertilizer spreading machinery structures to increase the spreading width and fertilizer utilization rate under the premise of low power consumption deserves more in-depth research.
Therefore, this article combined the gravity inline fertilizer spreading method and proposed a corn-overlapped strip fertilizer spreader without an external power source.Based on configuring a passive overlapping spreading process with a three-branch split chamber structure, the uniform spreading of fertilizer in strips is achieved.Then, based on the theoretical analysis of the fertilizer particle motion model, the main factors affecting the distribution pattern of fertilizers are determined.A single-factor, ternary orthogonal rotational combined response surface simulation test was conducted using discrete element (EDEM) software, and the analysis of variance was completed for the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient and determined the best parameters for simulation testing.The validation test was verified with 3D printing technology [15] to provide a reference for maize mid-tillage fertilizer chasing machinery design.

Operating Principle of the Corn-Overlapped Strip Fertilizer Spreader
The structure of the corn-overlapped strip fertilizer spreading device (Figure 1) mainly consists of fattening small plates, upper fertilizer spreading tube, lower fertilizer spreading tube, fertilizer funnel, fertilizer cone, and fattening bulges.The upper fertilizer spreading tube is axially equipped with a fertilizer funnel and cone, the lower fertilizer spreading tube is circumferentially evenly distributed with three fertilizer-type cavities, and the lower fertilizer spreading tube is solidly connected with fattening small plates with fattening bulges.
The structure of the corn-overlapped strip fertilizer spreading device (Figure 1) mainly consists of fattening small plates, upper fertilizer spreading tube, lower fertilizer spreading tube, fertilizer funnel, fertilizer cone, and fattening bulges.The upper fertilizer spreading tube is axially equipped with a fertilizer funnel and cone, the lower fertilizer spreading tube is circumferentially evenly distributed with three fertilizer-type cavities, and the lower fertilizer spreading tube is solidly connected with fattening small plates with fattening bulges.In the process of working, the fertilizer falls by gravity through the fertilizer discharge device into the upper fertilizer spreading tube.The fertilizer funnel collects the fertilizer so that the fertilizer flows to the fertilizer cone after the fertilizer funnel gathers flow.The fertilizer washes into the three cavities of the lower fertilizer tube through the first diversion of the fertilizer cone, then flows to the fattening small plates, achieves the second diversion under the action of the fattening bulge on the fattening small plates, and finally is thrown to the crops on both sides of the monopoly furrow.
When the machine moves forward in operation, three fattening small plates correspond to three single plates of the fertilizer spreading area, which form the left spreading belt, middle spreading belt, and right spreading belt, as shown in Figure 2. To ensure that the fertilizer spreading covers the fertilization monopoly in all directions and avoid leakage, the left applying belt and the right spreading belt should be overlapped with the middle spreading belt to ensure that the fertilizer spreading set becomes a surface.In the process of working, the fertilizer falls by gravity through the fertilizer discharge device into the upper fertilizer spreading tube.The fertilizer funnel collects the fertilizer so that the fertilizer flows to the fertilizer cone after the fertilizer funnel gathers flow.The fertilizer washes into the three cavities of the lower fertilizer tube through the first diversion of the fertilizer cone, then flows to the fattening small plates, achieves the second diversion under the action of the fattening bulge on the fattening small plates, and finally is thrown to the crops on both sides of the monopoly furrow.
When the machine moves forward in operation, three fattening small plates correspond to three single plates of the fertilizer spreading area, which form the left spreading belt, middle spreading belt, and right spreading belt, as shown in Figure 2. To ensure that the fertilizer spreading covers the fertilization monopoly in all directions and avoid leakage, the left applying belt and the right spreading belt should be overlapped with the middle spreading belt to ensure that the fertilizer spreading set becomes a surface.

Analysis of Fertilizer Spreading Motion
The fertilizer particles enter the fertilizer spreading tube with a certain initial velocity, collide randomly with the wall surface, and are finally thrown from the fattening small plates at different angles.The distribution of the fertilizer particles changes when the throwing speed and throwing height change.Referring to the study of Qingjin Lv [16,17] et al. on the vertical organic fertilizer spiral spreading device, the up-throw phase model and the down-throw phase model of the spiral spreading device were determined.This study is based on the down-throw phase model but does not involve the analysis of the specific collision motion of the fertilizer particles within the spreading device.

Analysis of Fertilizer Spreading Motion
The fertilizer particles enter the fertilizer spreading tube with a certain initial velocity, collide randomly with the wall surface, and are finally thrown from the fattening small plates at different angles.The distribution of the fertilizer particles changes when the throwing speed and throwing height change.Referring to the study of Qingjin Lv [16,17] et al. on the vertical organic fertilizer spiral spreading device, the up-throw phase model and the down-throw phase model of the spiral spreading device were determined.This study is based on the down-throw phase model but does not involve the analysis of the specific collision motion of the fertilizer particles within the spreading device.

Analysis of the Spatial Dispersion Motion of Fertilizer Particles
During the fertilizer throwing motion, the swing angle between the fattening small plates and the vertical plane α determines the throwing direction after contact between the particles and the surface.If α is too small, the vertical fall of the particles is not conducive to throwing, and if α is too large, the fertilizer pile in the fertilizer spreader tube is likely to be blocked, and the fertilizer particles tend to become stuck when the friction angle is greater than [18,19].Therefore, reasonable parameters need to be set to achieve the design requirements.
In the process of fertilizer particle throwing, it is necessary to carry out component processing in the X, Y, and Z directions, respectively, as shown in Figure 3. β indicates the combined throwing range angle of the fertilizer particles at the three outlets, and in this paper, the A-A surface is selected within the β throwing range angle for analysis, and main structural parameters in the figure: M indicates the width of the fattening small plates, N indicates the width of the bottom edge of the fattening small plates, α indicates the angle between the fattening small plates and the vertical plane pendulum, V 1 indicates the forward speed of the machine, L indicates the farthest distance the fertilizer can be thrown, L 1 indicates the horizontal distance from the center of the fertilizer spreading tube to the bottom edge of the fertilizer distribution plate, L 2 indicates the length of fertilizer veneer throwing, and L 3 indicates the outer diameter of lower fertilizer spreading tube.According to the geometric relationship in the horizontal throwing motion in Figure 3: Combining the structural parameters of existing surface mini-medium tillage fertilizer application products [20] and set fertilizer spreading tube inner diameter 48 mm, outer diameter 50 mm, and spreading tube width 55 mm, substituting into Equation (1) yields the following relationship: Referring to the agronomic standard for maize mid-tillage monopoly spacing [21] and set corn monopoly spacing of 600 mm, the farthest distance of fertilizer spreading should also satisfy the relationship: According to the geometric relationship in the horizontal throwing motion in Figure 3: Combining the structural parameters of existing surface mini-medium tillage fertilizer application products [20] and set fertilizer spreading tube inner diameter 48 mm, outer diameter 50 mm, and spreading tube width 55 mm, substituting into Equation (1) yields the following relationship: Referring to the agronomic standard for maize mid-tillage monopoly spacing [21] and set corn monopoly spacing of 600 mm, the farthest distance of fertilizer spreading should also satisfy the relationship: Since the inclination angle α of the angle between the fattening small plates and the vertical plane pendulum and the width N of the width of the bottom edge of the fattening small plates are mutually constrained, the range of the width N of the bottom edge of the fattening small plates and L 2, the width of fertilizer veneer throwing, needs to be further investigated in conjunction with the vertical throwing motion of the fertilizer.

Analysis of the Vertical Throwing Motion of Fertilizer Particles
The upper end of this fertilizer spreading device is connected to the fertilizer discharger.The fertilizer blending device based on EDEM software by Chen Haitao [22,23] et al. is similar to this investigation in that its fertilizer allocation device is connected to the fertilizer discharger, and during its analysis of the discharge velocity of the discharger outlet, the velocity at the discharge port of the discharger is derived in the theoretical output velocity range of 0.44 m/s-0.57m/s.As shown in Figure 4, the fertilizer enters the fertilizer spreading tube at speed V in , with the height H 1 of the fertilizer particles from the inlet to the outlet of the spreading device and energy loss E 1 generated during the impact with the inner wall surface of the fertilizer spreading tube, while the process increases the conversion of gravitational potential energy into kinetic energy, and finally, at the outlet of fattening small plates, it is thrown out at speed V out .The kinetic energy of its exit is calculated as: 1  2 mV 2 out = 1 2 mV 2 in + 2mgH 1 − E 1 , and the outlet velocity V out is obtained as Appl.Sci.2023, 13, x FOR PEER REVIEW 6 of 19 calculated as: , and the outlet velocity Vout is obtained as . Ignoring other motions, such as the rotation and collision of fertilizer particles in the fertilizer spreading device, and simplifying the vertical spreading movement of fertilizer, the vertical spreading motion is shown in Figure 4.The direction of fertilizer particles is set to falling in the Z direction, and the direction of movement of the corn along the side of the monopoly is the horizontal X direction.Vx and Vy denote the velocity of fertilizer particles moving in the X and Y directions, respectively; F denotes the sum of air resistance to fertilizer particles; and Fx and Fz signify air resistance to fertilizer particles in the X and Z directions, respectively.
Fertilizer pellets are thrown from the fattening small plates by gravity and air resistance.The air resistance applied is: Ignoring other motions, such as the rotation and collision of fertilizer particles in the fertilizer spreading device, and simplifying the vertical spreading movement of fertilizer, the vertical spreading motion is shown in Figure 4.The direction of fertilizer particles is set to falling in the Z direction, and the direction of movement of the corn along the side of the monopoly is the horizontal X direction.V x and V y denote the velocity of fertilizer particles moving in the X and Y directions, respectively; F denotes the sum of air resistance to fertilizer particles; and F x and F z signify air resistance to fertilizer particles in the X and Z directions, respectively.
Fertilizer pellets are thrown from the fattening small plates by gravity and air resistance.The air resistance applied is: where F is the air resistance of fertilizer particles, N; ρ is the air density, kg/m 3 ; S is the windward area of fertilizer particles, m 2 ; V is the velocity of fertilizer particles, m/s; C is the air resistance coefficient, m 2 ; and H is the throwing height of fertilizer particles, mm.Combining Newton s second law and related studies [24,25], the equation of motion of a fertilizer particle in the vertical plane is: where F z is the air resistance of fertilizer particles in the vertical plane, N; t is the time required to throw to the ground, s; and H z is the displacement of fertilizer particles in the vertical plane, m.It can be obtained as: The three fattening small plates are evenly distributed in a 120 • homogeneous direction; the throwing direction of the two side fattening small plates and the rear fattening small plates is different from the forward direction of the machine.Therefore, the motion of the fertilizer particles in the horizontal direction needs to be solved separately.
When the throwing direction of the fattening small plates is opposite to the forward direction of the machine, the equation of motion in the horizontal direction is: When the throwing direction of fattening small plates is the same as the forward direction of the machine, the equation of motion in the horizontal direction is: During the entire motion, the fertilizer particles are thrown at an initial position velocity V out = V 2 , then, substituting Equation (6) into Equations ( 7) and ( 8), respectively, the horizontal displacement L 2 of the fertilizer particles is found as: Then, the fertilizer granule spreading height H at the farthest spreading distance is: Combining Equations ( 2), ( 3) and ( 12), we can find that α > 30 • .The main structural parameters of the spreading width are the swing angle α of the spreading tube, the width N of the bottom edge of the fattening small plates, and the spreading height H.
The range of structural parameters was determined according to the height of the fuselage and the actual operational requirements; it was determined that 250 mm < height H < 450 mm, 40 mm < width N < 60 mm, and 30 • < pendulum angle α < 60 • according to the sliding friction coefficient of fertilizer.

Simulation Test Model and Parameters
This paper selects Dongping Lake urea fertilizer produced by Shandong Runyin Biological.The fertilizer had a water content of 0.89%, an actual density of 1.471 g/cm 3 , an average triaxial dimension of 2.24 mm × 2.22 mm × 2.24 mm, an equivalent diameter of 2.23 mm, and a sphericity of 0.975.Therefore, fertilizer particles with a similar equivalent diameter were selected to establish a spherical profile model.
As previously expressed, this paper builds a simulation model based on an agronomic model of corn planting with a monopoly spacing of 60 cm.The model was built using Solidworks, as shown in Figure 5, and consists mainly of the ground, corn root stubble, soil particles factory, fertilizer particles factory, and corn-overlapped fertilizer spreader.The model ground had a width of 2000 mm and a width of 1500 mm, corn root stubble evenly spaced at 60 cm, soil particles factory width of 2000 mm and a width of 1500 mm, height of 300 mm, and size from the ground of 100 mm.The fertilizer particles factory combines the outer slotted wheel fertilizer discharger with the dispensing funnel size [26], which is set as a cylinder with a diameter of 40 mm and a height of 60 mm in order to ensure that fertilizer does not accumulate during the operation.The main fertilizer spreading device parameters swing angle α, width N, and height H are changed according to the test requirements.The model was developed with Solidworks and saved as an IGS file, then imported into EDEM software.According to the different simulation objects, a suitable contact model should be s lected, combined with related research [27,28].This paper chooses the Hertz-Mindlin ( slip) model in EDEM as the contact model.During the fertilizer spreading process sim lation test, which simulates a natural operating environment, the soil first covers the lan surface in contact with the corn root stubble, and the fertilizer falls into contact with t soil and corn root stubble, respectively.As a result, contact occurs between fertilizer pa ticles and fertilizer particles, between fertilizer particles and soil particles, between fer lizer particles and corn roots, between fertilizer particles and fertilizer spreading devic between soil particles and fertilizer particles, and between soil particles and fertiliz spreading tubes.Referring to the relevant literature and research on the subject of corn ro stubble collision parameters [29,30], the simulation-related parameters were determined, According to the different simulation objects, a suitable contact model should be selected, combined with related research [27,28].This paper chooses the Hertz-Mindlin (no slip) model in EDEM as the contact model.During the fertilizer spreading process simulation test, which simulates a natural operating environment, the soil first covers the land surface in contact with the corn root stubble, and the fertilizer falls into contact with the soil and corn root stubble, respectively.As a result, contact occurs between fertilizer particles and fertilizer particles, between fertilizer particles and soil particles, between fertilizer particles and corn roots, between fertilizer particles and fertilizer spreading devices, between soil particles and fertilizer particles, and between soil particles and Appl.Sci.2023, 13, 2559 8 of 18 fertilizer spreading tubes.Referring to the relevant literature and research on the subject of corn root stubble collision parameters [29,30], the simulation-related parameters were determined, as shown in Table 1.The material of the fertilizer spreading device was set to PLA to facilitate prototype processing and manufacturing using 3D printing technology at a later stage.

Simulation Test Model and Test Index
A layer of soil needs to be laid on the ground to prevent the fertilizer particles from bouncing before the fertilizer spreading device starts operating.The Hertz-Mindlin (no slip) model in EDEM is chosen as the contact model.A virtual soil particles factory is first established, simplifying the soil particles to spherical particles with a radius of 7 mm and adding material properties to them, as shown in Table 1.Soil particles are generated at a rate of 200,000 particles/s and a total of 20,000 particles at a speed of 3 m/s along the negative direction of the Z-axis.
Then, fertilizer particles factory are created, fertilizer particles radius is set to its 1.165 mm spherical particles model, and material properties are added as shown in Table 1, with a total number of 100,000 and generation rate of 10,000/s.As the fertilizer particles in the fertilizer spreading device are in motion as the machine advances, the fertilizer particles need to be generated instantaneously and synchronously, and the lateral speed of the fertilizer particles needs to be aligned with the fertilizer spreading device during this process, as shown in Figure 6 Therefore, the positive x-axis velocity is set to 0.5 m/s in line with the travel speed of the fertilizer spreading device.At the same time, an initial velocity is given for the falling speed of the fertilizer particles by the above theoretical analysis process.For the simulation test analysis and calculation, the middle value of the above theoretical analysis of the inlet velocity V in velocity range of 0.5 m/s is taken, and the z-axis setting velocity is set to −0.5 m/s.ticles need to be generated instantaneously and synchronously, and the lateral sp the fertilizer particles needs to be aligned with the fertilizer spreading device duri process, as shown in Figure 6 Therefore, the positive x-axis velocity is set to 0.5 m/s with the travel speed of the fertilizer spreading device.At the same time, an initial v is given for the falling speed of the fertilizer particles by the above theoretical a process.For the simulation test analysis and calculation, the middle value of the theoretical analysis of the inlet velocity Vin velocity range of 0.5 m/s is taken, and axis setting velocity is set to −0.5 m/s.When the soil particle factory starts operation and when the soil is spread on the ground, 0 s is set, and 0.72 s is set after the fertilizer particles factory and the fertilizer spreading device then start the spreading operation.The fertilizer spreading device is located in the rightmost starting position of the land model, and the operation speed is synchronized with the X-axis direction of the fertilizer particle factory is 0.5 m/s, the simulation step length is 2.0 × 10 −5 s, and the data recording interval is 0.05 s.The total duration of the simulation is 3 s, and the simulation process is shown in Figure 7.When the soil particle factory starts operation and when the soil is spread on the ground, 0 s is set, and 0.72 s is set after the fertilizer particles factory and the fertilizer spreading device then start the spreading operation.The fertilizer spreading device is located in the rightmost starting position of the land model, and the operation speed is synchronized with the X-axis direction of the fertilizer particle factory is 0.5 m/s, the simulation step length is 2.0 × 10 −5 s, and the data recording interval is 0.05 s.The total duration of the simulation is 3 s, and the simulation process is shown in Figure 7.After the simulation, a fertilizer quality monitoring area of 2000 mm in length, 1500 mm in width, and 300 mm in height was set up on the ground surface, and the monitoring area was divided into fertilizer quality monitoring units of 150 mm in length, 200 mm in After the simulation, a fertilizer quality monitoring area of 2000 mm in length, 1500 mm in width, and 300 mm in height was set up on the ground surface, and the monitoring area was divided into fertilizer quality monitoring units of 150 mm in length, 200 mm in width, and 300 mm in height, as shown in Figure 8.According to the agronomic requirements of corn planting, referring to the relevant literature [31], the evaluation index of fertilizer uniformity is the uniformity coefficient.The exact moment of the monitoring area is selected, and the quality of fertilizer in the monitoring area of unit fertilizer quality is measured.The transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient are calculated and expressed by Y 1 and Y 2 , respectively.The smaller the uniformity coefficient, the more uniformly the fertilizer was spread.The larger the uniformity coefficient [32,33], the more unevenly the fertilizer was distributed, calculating the uniformity coefficient according to Equation (13).
where C V is the uniformity coefficient; m k is the fertilizer mass of the kth monitoring cell row, g; s c is the standard deviation; m is the mean value of the fertilizer mass collected in each column of the grid for the effective operational width, g; and K is the number of grid rows in the effective width region.
pl. Sci.2023, 13, x FOR PEER REVIEW ( ) where CV is the uniformity coefficient; mk is the fertilizer mass of the kth row, g; sc is the standard deviation; m is the mean value of the fertilizer each column of the grid for the effective operational width, g; and K is th rows in the effective width region.

Single-Factor Simulation Test
Single-factor tests were conducted with pendulum angle α, height H test factors and the transversal fertilizer uniformity coefficient (TFUN)Y1 fertilizer uniformity coefficient(LFUC) Y2 as test indexes.Based on the th of the parameter range, five test values were set for every aspect of the fie remaining elements were kept constant in the middle of their respectiv Design-expert 8.0.6 software was applied for data processing and st [34,35], the results of the tests are shown in Table 2, and the trends of the ef  Design-expert 8.0.6 software was applied for data processing and statistical analysis [34,35], the results of the tests are shown in Table 2, and the trends of the effects are shown in Figure 9.  Design-expert 8.0.6 software was applied for data processing and statistical analysis [34,35], the results of the tests are shown in Table 2, and the trends of the effects are shown in Figure 9.

Response Surface Simulation Test
The single-factor test results showed that the angle α, width N, and height H had significant effects on the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient.To further optimize the fertilizer spreading performance, the degree of influence of the above factors and their interaction terms were investigated, and the ternary orthogonal rotational combination response surface simulation test was carried out with angle α, height H, and width N as test factors, denoted by X 1 , X 2 , and X 3 , respectively, and the transversal fertilizer uniformity coefficient Y 1 and longitudinal fertilizer uniformity coefficient Y 2 as test indexes.The test factor code table is shown in Table 3, the simulation results and the test protocol table are shown in Table 4, and the analysis of variance is shown in Tables 5 and 6.The variance of the model the of transversal fertilizer uniformity coefficient Y 1 is shown in Table 5.The model was highly significant, the significance test of the model F 1 = 4.44, p < 0.01, and the misfit term was not significant F 2 = 9.38, p > 0.01, indicating that the regression model is highly significant and fits well within the test.The effects of angle X 1 , height X 2 , and width X 3 were significant (0.01 < p < 0.05), the interaction term X 1 X 3 had an effect (0.05 < p < 0.1), X 3 2 had an extremely significant effect (p < 0.01), and the remaining terms did not have a significant effect on this test index (p > 0.1).
The variance of the model of the longitudinal fertilizer uniformity coefficient Y 2 is shown in Table 6.The model was highly significant, the significance test of the model F 1 = 13.65,p < 0.01, and the misfit term was not significant F 2 = 1.89, p > 0.01, indicating that the regression model is highly significant and fits well within the test.The effects of angle X 1 and width X 3 were significant (0.01 < p < 0.05), the squared terms X 1 2 , X 2 2 , and X 3 2 were extremely significant (p < 0.01), and the remaining terms had no significant effect on this test index (p > 0.1).
Removing the non-significant term, the regression equation of the coefficient of variation of transversal fertilizer uniformity coefficient Y 1 with each factor was: Removing the non-significant term, the regression equation of the coefficient of variation of longitudinal fertilizer uniformity coefficient Y 2 with each factor was:

Response Surface Analysis
The single-factor test analysis showed that the transversal fertilizer uniformity coefficient Y 1 and longitudinal fertilizer uniformity coefficient Y 2 tended to decrease as the height X 2 decreased and reached the lowest value at 400 mm, so the response surface analysis of the interaction term X 1 X 3 was conducted when the height X 2 was selected as 400 mm, as shown in Figure 10.With the angle X 1 increased, the transversal fertilizer uniformity coefficient Y 1 showed a decreasing trend and subsequently became smaller, indicating that the spreading uniformity became better; as the width X 3 increased, the decreasing trend of the transversal fertilizer uniformity coefficient Y 1 slowed down compared with the angle X 1 , but also subsequently became smaller, indicating that the spreading uniformity became better and better with the increase in the width X 3 .

Parameter Optimization
According to the agronomic requirements of corn planting, the coefficient of variation of fertilizer spreading uniformity needs to meet the requirements of NY/T 1003-2006 <Technical Specification for Fertilizer Application Machinery Quality Evaluation> [36,37], and in order to seek the optimal combination between relevant factors, the height X2 is determined to be 400 mm, the transversal fertilizer uniformity coefficient Y1 and longitudinal fertilizer uniformity coefficient Y2 are less than 0.15 as the optimization index, and the solved range of parameters is the constraint for the optimization solution.The optimization results are shown in Figure 11.In order to facilitate the actual processing to reduce the cost and the stability of the fertilizer spreading operation, the median swing angle X1 = 52°, height X2 = 400 mm, and minimum width X3 = 50 mm were determined as the optimal combination of parameters.
The simulation model was re-established with the optimized parameters, and the simulation validation test was conducted.The test was repeated three times, and the test results were averaged.The test results showed that the optimized transversal fertilizer uniformity coefficient Y1 and longitudinal fertilizer uniformity coefficient Y2 are 0.132 and 0.135, which are less than 0.15 and are within the optimization criterion and prove that the optimized simulation parameters are correct.

Parameter Optimization
According to the agronomic requirements of corn planting, the coefficient of variation of fertilizer spreading uniformity needs to meet the requirements of NY/T 1003-2006 <Technical Specification for Fertilizer Application Machinery Quality Evaluation> [36,37], and in order to seek the optimal combination between relevant factors, the height X 2 is determined to be 400 mm, the transversal fertilizer uniformity coefficient Y 1 and longitudinal fertilizer uniformity coefficient Y 2 are less than 0.15 as the optimization index, and the solved range of parameters is the constraint for the optimization solution.The optimization results are shown in Figure 11.

Parameter Optimization
According to the agronomic requirements of corn planting, the coefficient of variation of fertilizer spreading uniformity needs to meet the requirements of NY/T 1003-2006 <Technical Specification for Fertilizer Application Machinery Quality Evaluation> [36,37], and in order to seek the optimal combination between relevant factors, the height X2 is determined to be 400 mm, the transversal fertilizer uniformity coefficient Y1 and longitudinal fertilizer uniformity coefficient Y2 are less than 0.15 as the optimization index, and the solved range of parameters is the constraint for the optimization solution.The optimization results are shown in Figure 11.In order to facilitate the actual processing to reduce the cost and the stability of the fertilizer spreading operation, the median swing angle X1 = 52°, height X2 = 400 mm, and minimum width X3 = 50 mm were determined as the optimal combination of parameters.
The simulation model was re-established with the optimized parameters, and the simulation validation test was conducted.The test was repeated three times, and the test results were averaged.The test results showed that the optimized transversal fertilizer uniformity coefficient Y1 and longitudinal fertilizer uniformity coefficient Y2 are 0.132 and 0.135, which are less than 0.15 and are within the optimization criterion and prove that the optimized simulation parameters are correct.In order to facilitate the actual processing to reduce the cost and the stability of the fertilizer spreading operation, the median swing angle X 1 = 52 • , height X 2 = 400 mm, and minimum width X 3 = 50 mm were determined as the optimal combination of parameters.
The simulation model was re-established with the optimized parameters, and the simulation validation test was conducted.The test was repeated three times, and the test results were averaged.The test results showed that the optimized transversal fertilizer uniformity coefficient Y 1 and longitudinal fertilizer uniformity coefficient Y 2 are 0.132 and 0.135, which are less than 0.15 and are within the optimization criterion and prove that the optimized simulation parameters are correct.

Validation Test
A validation test of the corn-overlapped strip spreader was carried out to verify the feasibility of the simulation method mentioned above and the correctness of the optimization results of the corn-overlapped strip spreader.The experiment was conducted in September 2021 at the Intelligent Agricultural Equipment Engineering Laboratory of Harbin Cambridge University.The experiment was conducted realistically by replicating the grid division of the statistical simulation area, with a total of 160 fertilizer collection points in 8 columns and 20 rows in a test area of 2 m × 1.2 m, with a row spacing of 15 cm and a column spacing of 20 cm.
The experimental procedure is shown in Figure 12 The test device mainly consisted of a fertilizer tank, a fertilizer discharger, a corn-overlapped strip fertilizer spreader, and a walking device, and the corn-overlapped strip spreader was obtained through a 3D printing process [38,39].A chain drove the walking device and the fertilizer spreader, and the machine s forward speed was controlled by adjusting the transmission ratio to 0.5 m/s before the test.Before the test, five groups of the same mass of fertilizer were weighed and poured into the fertilizer tank.The front wheel of the whole machine was located at the middle starting line of the test grid along the center line of the grid to advance the operation.At the end of the test, the fertilizer in each grid was collected with a brush and then weighed and counted, and each fertilizer collection bag was labeled with the corresponding grid position to facilitate the tallying of test results.The test was repeated five times, the data from each test were weighed with an electronic balance, and a statistical interval was selected to calculate its horizontal and vertical uniformity coefficients, as shown in Table 7.

Validation Test
A validation test of the corn-overlapped strip spreader was carried out to verify the feasibility of the simulation method mentioned above and the correctness of the optimization results of the corn-overlapped strip spreader.The experiment was conducted in September 2021 at the Intelligent Agricultural Equipment Engineering Laboratory of Harbin Cambridge University.The experiment was conducted realistically by replicating the grid division of the statistical simulation area, with a total of 160 fertilizer collection points in 8 columns and 20 rows in a test area of 2 m × 1.2 m, with a row spacing of 15 cm and a column spacing of 20 cm.
The experimental procedure is shown in Figure 12 The test device mainly consisted of a fertilizer tank, a fertilizer discharger, a corn-overlapped strip fertilizer spreader, and a walking device, and the corn-overlapped strip spreader was obtained through a 3D printing process [38,39].A chain drove the walking device and the fertilizer spreader, and the machine′s forward speed was controlled by adjusting the transmission ratio to 0.5 m/s before the test.Before the test, five groups of the same mass of fertilizer were weighed and poured into the fertilizer tank.The front wheel of the whole machine was located at the middle starting line of the test grid along the center line of the grid to advance the operation.At the end of the test, the fertilizer in each grid was collected with a brush and then weighed and counted, and each fertilizer collection bag was labeled with the corresponding grid position to facilitate the tallying of test results.The test was repeated five times, the data from each test were weighed with an electronic balance, and a statistical interval was selected to calculate its horizontal and vertical uniformity coefficients, as shown in Table 7.The results showed that the corn-overlapped strip fertilizer spreading device could complete the spreading operation smoothly, and the average values of the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient were 0.1427 and 0.1471, with deviation values of 8.11% and 9.01%, respectively, from the simulation test, which was in line with the ±10% deviation range.A possible reason for the deviation might be the deviation of the urea particles from the ideal particles in the simulation, and the bouncing the urea particles on the cement ground during the actual operation, resulting in a slight deviation of the fertilizer spreading process on the cement ground from the simulation.The departure was not significant, indicating that the whole machine worked well.

Conclusions
(1) In this study, we designed a corn-overlapped strip fertilizer spreading device in without an external power source.By configuring a passive overlapping spreading method with a three-branch split chamber structure, uniform spreading of fertilizer in strips was achieved.Based on the theoretical analysis of the fertilizer particle motion model, the main factors affecting the distribution pattern of fertilizers were angle α, width N, and height H. (2) The single-factor ternary orthogonal rotational combination response surface simulation test was carried out with angle α, width N, and height H as test factors and the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient as test indicators.The regression model was established using Design-expert8.0.6 software to derive the variation relationship of the test factors on the test indexes.The test results showed that the optimized transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient were less than 0.15% when the pendulum angle = 52 • , height = 400 mm, and width = 50 mm, which were under the optimization criterion.(3) A verification test was carried out under the optimal combination of parameters for the simulation tests with the simulation conditions as the standard.The test results were consistent with the simulation results within the error range, and the deviation values of the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient were 8.11% and 9.01%, respectively, which was in line with the ± 10% deviation range.The corn-overlapped strip fertilizer spreader was able to complete the fertilizer spreading operation smoothly.

Patents
Two patents have been applied in China in this manuscript (Patent No. ZL202121578688.3No. ZL202021054986.8).
Appl.Sci.2023,13,  x FOR PEER REVIEW 8 of spreader.The model ground had a width of 2000 mm and a width of 1500 mm, corn ro stubble evenly spaced at 60 cm, soil particles factory width of 2000 mm and a width 1500 mm, height of 300 mm, and size from the ground of 100 mm.The fertilizer particl factory combines the outer slotted wheel fertilizer discharger with the dispensing funn size[26], which is set as a cylinder with a diameter of 40 mm and a height of 60 mm order to ensure that fertilizer does not accumulate during the operation.The main fer lizer spreading device parameters swing angle α, width N, and height H are changed a cording to the test requirements.The model was developed with Solidworks and sav as an IGS file, then imported into EDEM software.

Figure 8 .
Figure 8. Schematic diagram of fertilizer monitoring area.

Figure 8 .
Figure 8. Schematic diagram of fertilizer monitoring area.

3. 3 .
Parameter Optimization Test 3.3.1.Single-Factor Simulation Test Single-factor tests were conducted with pendulum angle α, height H, and width N as test factors and the transversal fertilizer uniformity coefficient (TFUN)Y 1 and longitudinal fertilizer uniformity coefficient(LFUC) Y 2 as test indexes.Based on the theoretical analysis of the parameter range, five test values were set for every aspect of the field of values.The remaining elements were kept constant in the middle of their respective values degrees.

Figure 8 .
Figure 8. Schematic diagram of fertilizer monitoring area.

Figure 9 .Figure 9 .
Figure 9.The influence curve of the single factor test.(a) Effect of angle on the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient, (b) effect of width on theFigure 9.The influence curve of the single factor test.(a) Effect of angle on the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient, (b) effect of width on the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient, (c) effect of height on the transversal fertilizer uniformity coefficient and longitudinal fertilizer uniformity coefficient.

Table 1 .
Global variable parameter setting.

Table 2 .
Single factor test results of swing angle, height, and width.

Table 3 .
Factor level coding table.

Table 4 .
Simulation test scheme and results.

Table 5 .
Analysis of variance for lateral uniformity coefficient.

Table 6 .
Analysis of variance of average velocity of fertilizer particles.