Optimization and Experiment of Fertilizer-Spreading Device for Wheat Wide-Boundary Sowing Planter under Full Rice Straw Retention

: When sowing with a wide boundary under full rice straw retention in the rice–wheat rotation area of China, conventional fertilization methods have some problems, such as a low fertilizer utilization rate, heap soil around a buried fertilizer device, or blocked fertilizing oriﬁce. Firstly, combined with theoretical analysis, discrete element numerical simulation technology, and central composite test method, the wide-boundary fertilization device for wheat wide-boundary sowing was designed. Then, with the coefﬁcient of variation for particle uniformity (CVPU) as the response value, the central composite experiment was carried out on the key structural parameters (focal length coefﬁcient, lateral span, tilt angle, and ground clearance) of the wide-boundary fertilization device by EDEM software. Finally, the inﬂuential rules of core factors of the device on the CVPU were analyzed by Design-Expert software; then, the optimal parameter combination was determined and veriﬁed by a ﬁeld test. The results showed that all factors had signiﬁcant effects on the CVPU. The primary and secondary factors affecting the CVPU were the tilt angle, lateral span, focal length coefﬁcient, ground clearance, tilt angle × ground clearance, and lateral span × ground clearance, in which there were certain interactions between the tilt angle and ground clearance and lateral span and ground clearance. When the focal length coefﬁcient, lateral span, tilt angle, and ground clearance were 1.5, 60 mm, 30 ◦ , and 192 mm, respectively, the lateral was minimum. In this case, the theoretical value and ﬁeld test value were 14.11% and 17.63%, respectively. The ﬁeld test value is consistent with the theoretical calculation value. This study could provide references for the design of a fertilizer-spreading device with a wide boundary.


Introduction
The rice-wheat rotation system is the main farming system in the Yangtze River Delta region in China. When sowing wheat in this agronomic model, sowing should be completed 2 weeks after the rice is harvested. The sown area and yield of wheat under this agronomic model account for 20% and 22% of China's total wheat sown area and total wheat yield, respectively. It is one of the main modes of grain production in China [1,2]. However, with the background of the ban on straw burning by Chinese law and the change in China's agricultural production model, a large amount of straw has been left in the fields where wheat will be sown after the rice harvest in a rice-wheat rotation system. This kind of land is called primitive land with full rice straw (shown in Figure 1), that is, the field without any straw removal and tillage treatment after the previous rice is harvested, which has become a common phenomenon in a rice-wheat rotation area [3]. When producing wheat in primitive land with full rice straw, wheat sowing cannot be producing wheat in primitive land with full rice straw, wheat sowing cannot be compl on time, which leads to a reduction in yield due to a short sowing time, a large amou straw in the field, and many production procedures [4]. Wide-boundary sowing tech ogy with crushed straw inter-row mulching can effectively solve the above problems. is an integrated operation of technology involving straw returning to the field and w boundary sowing, in which the straw is pulverized and covered on inter-rows while ing wheat [5][6][7]. Meanwhile, the straw-returning mode with crushed straw inter mulching can maintain moisture and hinder weeds' growth, which is conducive to implementation of conservation tillage in rice-wheat rotation areas [8][9][10][11]. International research on the technology is mainly focused on the developme straw returning to the field or sowing in a clayish-wet soil environment, but the tech ogy of wide-boundary fertilization is rarely studied. The current fertilization method the sowing environment are still large-area fertilization with centrifugation or strip lization with ditching.
Przywara et al. [12] clarified the influence of structural and operational param of a centrifugal disc spreader on the spatial distribution of different physical properti fertilizer. Hwang et al. [13] analyzed the distribution pattern of a centrifugal fertilizer tributor with two shutter holes and obtained the optimal optimum structural param using the discrete element method (DEM) for improved distribution uniformity. To prove the operation performance and fertilizer utilization rate of centrifugal fertilizer tributors, different scholars have studied different technologies, such as a centralized tilizer application with an air-driven fertilizer distribution control system, the influ of key parameters of fertilizer spinner on uniformity, and the optimization of structu working parameters [14][15][16][17][18][19]. Coetzee et al. [20] designed a deflector plate to separat flow into two paths along the tree lines; Then, the effects of the key parameters on spread pattern were investigated experimentally and compared to DEM results. Liu [21] optimized the structural parameters of a strip-shunt plate using EDEM software an orthogonal test. Ding et al. [22] developed dual-frequency ditching equipment w can apply starter fertilizer and base fertilizer at the same time. Additionally, they mized the structural parameters with EDEM software. To meet the requirements of fertilization with ditching under different working conditions, Zhu et al. [23], Xiao [24], and Du et al. [25] designed fertilization devices with different structural forms analyzed the impact of key parameters on operation performance.
In a primitive land with full rice straw, there are many problems when using con tional methods to apply fertilizer [26,27]. Large-area fertilization using centrifugal lizer distributors will cause many problems, such as the need for many operations a serious waste of fertilizer. When fertilizing with strip fertilization with ditching, it is for the discharge outlet of submerged components to be blocked by soil due to st stickiness and high moisture of the soil. International research on the technology is mainly focused on the development of straw returning to the field or sowing in a clayish-wet soil environment, but the technology of wide-boundary fertilization is rarely studied. The current fertilization methods in the sowing environment are still large-area fertilization with centrifugation or strip fertilization with ditching.
Przywara et al. [12] clarified the influence of structural and operational parameters of a centrifugal disc spreader on the spatial distribution of different physical properties of fertilizer. Hwang et al. [13] analyzed the distribution pattern of a centrifugal fertilizer distributor with two shutter holes and obtained the optimal optimum structural parameters using the discrete element method (DEM) for improved distribution uniformity. To improve the operation performance and fertilizer utilization rate of centrifugal fertilizer distributors, different scholars have studied different technologies, such as a centralized fertilizer application with an air-driven fertilizer distribution control system, the influence of key parameters of fertilizer spinner on uniformity, and the optimization of structure or working parameters [14][15][16][17][18][19]. Coetzee et al. [20] designed a deflector plate to separate the flow into two paths along the tree lines; Then, the effects of the key parameters on the spread pattern were investigated experimentally and compared to DEM results. Liu et al. [21] optimized the structural parameters of a strip-shunt plate using EDEM software and an orthogonal test. Ding et al. [22] developed dual-frequency ditching equipment which can apply starter fertilizer and base fertilizer at the same time. Additionally, they optimized the structural parameters with EDEM software. To meet the requirements of strip fertilization with ditching under different working conditions, Zhu et al. [23], Xiao et al. [24], and Du et al. [25] designed fertilization devices with different structural forms and analyzed the impact of key parameters on operation performance.
In a primitive land with full rice straw, there are many problems when using conventional methods to apply fertilizer [26,27]. Large-area fertilization using centrifugal fertilizer distributors will cause many problems, such as the need for many operations and a serious waste of fertilizer. When fertilizing with strip fertilization with ditching, it is easy for the discharge outlet of submerged components to be blocked by soil due to strong stickiness and high moisture of the soil.
To meet the demand for wide-boundary fertilization during wheat sowing on primitive land with full rice straw, combined with the working principle of wide-boundary sowing technology for crushed straw inter-row mulching, this study designed a fertilizing device for wide-boundary sowing. Combined with theoretical analysis and discrete element numerical simulation, the structure parameters of the fertilization device were optimized, and a field verification test was carried out. This study can provide a theoretical reference for the design of wheat planters applied to primitive land with full rice straw.

Agronomic Model of Wide-Boundary Sowing of Wheat after Rice
It is common for problems to arise-such as hanging grass on equipment, planting apparatus blockage, and low quality of covering soil-when sowing wheat by strip sowing or hole seeding due to high moisture, strong stickiness of soil, and a large amount of straw in primitive land with full rice straw. Wheat seeding with a wide boundary is a new wheat-seeding technology that combines strip and broadcasting without ridges. It has the advantages of optimizing the population structure and improving the leaf area index. Combining the working conditions of primitive land with full rice straw and the agronomic requirements for wheat seeding with a wide boundary, a sowing model with a sowing width of 240 mm and a spacing of 300 mm was proposed. It relies on a wide-boundary planter with crushed straw inter-row mulching to complete the operation. The agronomic requirements are as follows: three tasks, including crushing straw, laying straw in strips, and sowing in a wide-boundary condition, should be completed in one operation; the width of the crushed straw strip should not exceed 300 mm, and the width of the "no-straw" sowing belt between the adjacent straw belts should not be less than 240 mm. Fertilization and sowing are carried out in the full width of the 240 mm "no-straw" seeding belt, with a depth of 30-50 mm, and the post-sowing seed belt is compacted. The agronomic pattern for wide-boundary sowing of wheat after rice is shown in Figure 2.
To meet the demand for wide-boundary fertilization during wheat sowing on pri itive land with full rice straw, combined with the working principle of wide-bounda sowing technology for crushed straw inter-row mulching, this study designed a fertilizi device for wide-boundary sowing. Combined with theoretical analysis and discrete e ment numerical simulation, the structure parameters of the fertilization device were op mized, and a field verification test was carried out. This study can provide a theoreti reference for the design of wheat planters applied to primitive land with full rice straw

Agronomic Model of Wide-Boundary Sowing of Wheat after Rice
It is common for problems to arise-such as hanging grass on equipment, planti apparatus blockage, and low quality of covering soil-when sowing wheat by strip so ing or hole seeding due to high moisture, strong stickiness of soil, and a large amount straw in primitive land with full rice straw. Wheat seeding with a wide boundary is a ne wheat-seeding technology that combines strip and broadcasting without ridges. It has t advantages of optimizing the population structure and improving the leaf area ind Combining the working conditions of primitive land with full rice straw and the ag nomic requirements for wheat seeding with a wide boundary, a sowing model with sowing width of 240 mm and a spacing of 300 mm was proposed. It relies on a wid boundary planter with crushed straw inter-row mulching to complete the operation. T agronomic requirements are as follows: three tasks, including crushing straw, layi straw in strips, and sowing in a wide-boundary condition, should be completed in o operation; the width of the crushed straw strip should not exceed 300 mm, and the wid of the "no-straw" sowing belt between the adjacent straw belts should not be less than 2 mm. Fertilization and sowing are carried out in the full width of the 240 mm "no-straw seeding belt, with a depth of 30-50 mm, and the post-sowing seed belt is compacted. T agronomic pattern for wide-boundary sowing of wheat after rice is shown in Figure 2.

Machine Structure and Working Principle
The wide-boundary planter with crushed straw inter-row mulching is mainly vided into two key processes: straw returning treatment and wheat sowing with a wi boundary. Its key operating components mainly include a straw crushing device, stra diversion device, wide-boundary fertilization device, rotary tillage device, and seedi device.
The structure of the machine is shown in Figure 3, and the main technical paramete are shown in Table 1. The working principle is as follows: the crushed straw is spray from the rear cover of the straw-crushing device and then divided by the straw deflect The purpose of this process is that a "no-straw" seed strip is formed directly below t wide-boundary fertilization device, and a regular crushed straw belt is formed direc below the interval between adjacent wide-boundary fertilization devices. Meanwhile, f tilization is carried out by the wide-boundary fertilization device located under the stra deflector on the full width of the 240 mm "no-straw" seeding belt. Then, the wide-boun ary seeding device carries out strip rotary plowing and wide-boundary sowing on t seeding belt covered with fertilizer. After the operation, the field showed an alternati

Machine Structure and Working Principle
The wide-boundary planter with crushed straw inter-row mulching is mainly divided into two key processes: straw returning treatment and wheat sowing with a wide boundary.
Its key operating components mainly include a straw crushing device, straw diversion device, wide-boundary fertilization device, rotary tillage device, and seeding device.
The structure of the machine is shown in Figure 3, and the main technical parameters are shown in Table 1. The working principle is as follows: the crushed straw is sprayed from the rear cover of the straw-crushing device and then divided by the straw deflector. The purpose of this process is that a "no-straw" seed strip is formed directly below the wideboundary fertilization device, and a regular crushed straw belt is formed directly below the interval between adjacent wide-boundary fertilization devices. Meanwhile, fertilization is carried out by the wide-boundary fertilization device located under the straw deflector on the full width of the 240 mm "no-straw" seeding belt. Then, the wide-boundary seeding device carries out strip rotary plowing and wide-boundary sowing on the seeding belt covered with fertilizer. After the operation, the field showed an alternating longitudinal arrangement of a 240 mm "no-straw" seeding belt and a 300 mm crushed straw belt. The land surface after sowing is shown in Figure 4.
Agronomy 2022, 12, x FOR PEER REVIEW 4 longitudinal arrangement of a 240 mm "no-straw" seeding belt and a 300 mm cru straw belt. The land surface after sowing is shown in Figure 4.   The structure of a wide-boundary fertilization device is mainly composed of a f lizer pipe, a fertilizer even-distribution plate, and a mounting plate. It is located dir below the crushed straw deflector, fixed on the beam after the crushing device, and responds to four groups of sequential crushed straw deflectors. The working princip as follows: when fertilizers fall on the fertilizer even-distribution plate, disorderly bo ing movements occur under the blocking action of the arc plate so that the fertilizer p cles are uniformly distributed laterally on the "no-straw" seed belt; at the same time crushed straw deflector installed outside the wide-boundary fertilization device prev the fertilizer particles from spreading to crushed straw strips, so that all fertilizer part fall in the "no-straw" seed belt. An assembly drawing of the crushed straw deflector wide-boundary fertilization device is shown in Figure 5.  Agronomy 2022, 12, x FOR PEER REVIEW longitudinal arrangement of a 240 mm "no-straw" seeding belt and a 300 mm c straw belt. The land surface after sowing is shown in Figure 4.   The structure of a wide-boundary fertilization device is mainly composed of lizer pipe, a fertilizer even-distribution plate, and a mounting plate. It is located d The structure of a wide-boundary fertilization device is mainly composed of a fertilizer pipe, a fertilizer even-distribution plate, and a mounting plate. It is located directly below the crushed straw deflector, fixed on the beam after the crushing device, and corresponds to four groups of sequential crushed straw deflectors. The working principle is as follows: when fertilizers fall on the fertilizer even-distribution plate, disorderly bouncing movements occur under the blocking action of the arc plate so that the fertilizer particles are uniformly distributed laterally on the "no-straw" seed belt; at the same time, the crushed straw deflector installed outside the wide-boundary fertilization device prevents the fertilizer particles from spreading to crushed straw strips, so that all fertilizer particles fall in the "no-straw" seed belt. An assembly drawing of the crushed straw deflector and wide-boundary fertilization device is shown in Figure 5.

Analysis of Falling Position of Fertilizer Particles
The structure of the fertilizer even-distribution plate has a crucial influence lateral uniformity of fertilization. When fertilizer particles collide with the fertiliz distribution plate, the motion trajectory is changed by the collision forces so that t lizer flow can be dispersed evenly. The uniform distribution performance of the f even-distribution plate with plane structure is low, which leads to an accumul fertilizer particles in the middle of the seeding belt. To improve the performanc fertilizer even-distribution plate, a new even-distribution plate was designed study, whose cross-section is a quadratic curve and its longitudinal section is a line, as shown in Figure 6. Taking the intersection of the bottom surface, the sym plane, and the end vertical plane of the plate as the origin, the forward directio machine as the positive direction of the y-axis, the vertical upward as the positiv tion of the z-axis, and the right-hand rule to determine the positive direction of th the coordinate system oxyz as shown in Figure 7 was established. The quadratic c the bottom surface and the straight line of the symmetrical surface satisfy Relatio  tan tan  The structure of the fertilizer even-distribution plate has a crucial influence on the lateral uniformity of fertilization. When fertilizer particles collide with the fertilizer evendistribution plate, the motion trajectory is changed by the collision forces so that the fertilizer flow can be dispersed evenly. The uniform distribution performance of the fertilizer even-distribution plate with plane structure is low, which leads to an accumulation of fertilizer particles in the middle of the seeding belt. To improve the performance of the fertilizer even-distribution plate, a new even-distribution plate was designed in this study, whose cross-section is a quadratic curve and its longitudinal section is a straight line, as shown in Figure 6. Taking the intersection of the bottom surface, the symmetrical plane, and the end vertical plane of the plate as the origin, the forward direction of the machine as the positive direction of the y-axis, the vertical upward as the positive direction of the z-axis, and the right-hand rule to determine the positive direction of the x-axis, the coordinate system oxyz as shown in Figure 7 was established. The quadratic curve of the bottom surface and the straight line of the symmetrical surface satisfy Relation (1).
where p is the focal length coefficient; δ is the inclination angle of the device edge line in the symmetry plane, • ; and l is span, mm.
Agronomy 2022, 12, x FOR PEER REVIEW Figure 5. Assembly drawing of crushed straw deflector and wide-boundary fertilization d

Analysis of Falling Position of Fertilizer Particles
The structure of the fertilizer even-distribution plate has a crucial influence lateral uniformity of fertilization. When fertilizer particles collide with the fertilize distribution plate, the motion trajectory is changed by the collision forces so that t lizer flow can be dispersed evenly. The uniform distribution performance of the f even-distribution plate with plane structure is low, which leads to an accumul fertilizer particles in the middle of the seeding belt. To improve the performanc fertilizer even-distribution plate, a new even-distribution plate was designed study, whose cross-section is a quadratic curve and its longitudinal section is a line, as shown in Figure 6. Taking the intersection of the bottom surface, the symm plane, and the end vertical plane of the plate as the origin, the forward directio machine as the positive direction of the y-axis, the vertical upward as the positiv tion of the z-axis, and the right-hand rule to determine the positive direction of the the coordinate system oxyz as shown in Figure 7 was established. The quadratic c the bottom surface and the straight line of the symmetrical surface satisfy Relation  . Schematic diagram of fertilizer particle motion analysis.
A fertilizer particle was selected as the research object, and the contact point between the particle and the fertilizer even-distribution plate is point p 0 (x 0 , y 0 , z 0 ). If the tangent plane passing through the point p 0 intersects the quadratic curve at the point p 1 (x 1 , y 1 , 0), the following relation is established: where  is the angle between the projection of line cp 1 on the xoy plane and the y-axis, °.
AB is the tangent line of the quadratic curve at point p 1 (x 1 , y 1 , 0), and its tangent equation is 1 11 tan ( tan ) 50 Then, the intercepts a and b of line AB on the x and y axes are, respectively: Therefore, the contact point p 0 of the fertilizer particles on the fertilizer even-distribution plate satisfies the following equation: Figure 7. Schematic diagram of fertilizer particle motion analysis.
A fertilizer particle was selected as the research object, and the contact point between the particle and the fertilizer even-distribution plate is point p 0 (x 0 , y 0 , z 0 ). If the tangent plane passing through the point p 0 intersects the quadratic curve at the point p 1 (x 1 , y 1 , 0), the following relation is established: where θ is the angle between the projection of line cp 1 on the xoy plane and the y-axis, • .
AB is the tangent line of the quadratic curve at point p 1 (x 1 , y 1 , 0), and its tangent equation is Then, the intercepts a and b of line AB on the x and y axes are, respectively: a = 50y 1 py 1 tan θ + y 1 tan θ b = py 1 2 tan θ 2 50 + y 1 Therefore, the contact point p 0 of the fertilizer particles on the fertilizer even-distribution plate satisfies the following equation:

Fertilizer Particle Motion Analysis
Because the collision process is very short, and the collision force is much greater than gravity, friction, and air resistance, it is assumed that the displacement change of particles where v 0 is the initial velocity, m · s −1 ; v 1 is the rebound velocity, m · s −1 ; α is the angle between v 0 and line cp 1 , • ; β is the angle between v 1 and line cp 1 , • ; and e is the coefficient of restitution. Because the machine moves at a steady speed during the operation, it is considered that the fertilizer particles are uniformly distributed longitudinally. Therefore, only the lateral homogeneity of the fertilizer granules was investigated. We can decompose the rebound velocity v 1 at point p 0 in the directions of the coordinate axis: where v x1 is the x-axis component of the initial velocity of the particle in oblique throwing motion, m · s −1 ; and v z1 is the z-axis component of the initial velocity of the particle in oblique throwing motion, m · s −1 .
The fertilizer particle performs an oblique throwing motion after bouncing at point p 0 . It is assumed that the time when the fertilizer particle begins to perform oblique throwing motion is time t 0 = 0. Ignoring the spin of the particle, the differential equations of motion along the coordinate axis are established in the xoz plane.
where v x is the x-axis component of the velocity of the particle in oblique throwing motion, m · s −1 ; v z is the z-axis component of the velocity of the particle in oblique throwing motion, m · s −1 ; m is the quality of fertilizer particles, kg; k is the air resistance coefficient, which is related to the shape of the object and air density, etc. When v z > 0, n = 0; when v z < 0, n = 1; f xy is the air resistance of the fertilizer particles on the xoz plane.
The velocity equations of the particle in the directions of the x-axis and z-axis are shown in Equation (10). The kinematic equations of the particle in the directions of the x-axis and z-axis are shown in Equation (11).
Among them: Based on the analysis of Equations (3)-(12), the following can be observed: during the entry movement of fertilizer particles, the location of the collision point is mainly related to the tilt angle and curve shape of the even-distribution plate bottom, and the falling time of particles is mainly determined by the movement path in the z-axis direction, initial velocity, and ground clearance. When the initial velocity is determined, the lateral position of particles on the ground is mainly determined by the position of the collision point and the moving time. Therefore, the focal length coefficient, lateral span, tilt angle, and ground clearance directly affect the lateral uniformity performance of the fertilizer even-distribution plate.
Through a series of single factor tests, it was preliminarily determined that the range of focal length coefficient is 0.5~1.5, the range of lateral span is 35~95 mm, the range of tilt angle is 20~50 • , and the range of ground clearance is 120~260 mm.

Materials and Methods
In the process of wide-boundary fertilization, particles are subjected to the interaction force of even-distribution plate and particle flow, and the complex collision movement of particles occurs in the process. In this study, the simulation platform of fertilizer particles-device-soil tank system was built by EDEM software, and the kinematics simulation of wide-boundary fertilization process was carried out. Taking the coefficient of variation for fertilizer uniformity as the evaluation index, the influence of the core factors on operation performance was analyzed.

Construction of Discrete Element Platform
The simulation model includes fertilizer particle models, device geometry models, and the statistical plate model with soil properties. Since the spherical rate of fertilizer particles is more than 90%, 4 mm spherical particles are directly used as the simulation model of fertilizer particles. The simplified model of the crushed straw deflector and the wide-boundary fertilization device was imported into EDEM software, and the statistical plate model (3000 mm × 300 mm × 20 mm) was established under the devices. The simulation platform of particle-device-statistical plate system is shown in Figure 8. Combined with [16,19], the relevant parameters of particles, devices, and statistical plate are shown in Table 2.

Simulation Test Design
The operating speed of the device is generally 0.8~1.8 m/s, so the forward speed of the device and fertilizer particles-factory during the simulation process are both 1.3 m/s. The fertilizer quality required by agronomy for wheat production in the rice-wheat rotation system is 675~825 kg/hm. Combined with the simulation parameters of machine speed, the generation speed of fertilizer particles was set to 58.5 g/s. Taking focal length coefficient, Agronomy 2022, 12, 2251 9 of 15 lateral span, tilt angle, and ground clearance as the test factors and the coefficient of variation for fertilizer uniformity as the response value, the compound experiment of four factors and five levels was carried out. The test factors and horizontal codes are shown in Table 3. The test schemes and results are shown in Table 4.

Construction of Discrete Element Platform
The simulation model includes fertilizer particle models, device geometry models, and the statistical plate model with soil properties. Since the spherical rate of fertilizer particles is more than 90%, 4 mm spherical particles are directly used as the simulation model of fertilizer particles. The simplified model of the crushed straw deflector and the wide-boundary fertilization device was imported into EDEM software, and the statistical plate model (3000 mm × 300 mm × 20 mm) was established under the devices. The simulation platform of particle-device-statistical plate system is shown in Figure 8. Combined with [16,19], the relevant parameters of particles, devices, and statistical plate are shown in Table 2.

Simulation Test Design
The operating speed of the device is generally 0.8~1.8 m/s, so the forward speed of the device and fertilizer particles-factory during the simulation process are both 1.3 m/s. The fertilizer quality required by agronomy for wheat production in the rice-wheat rotation system is 675~825 kg/hm. Combined with the simulation parameters of machine speed, the generation speed of fertilizer particles was set to 58.5 g/s. Taking focal length coefficient, lateral span, tilt angle, and ground clearance as the test factors and the coefficient of variation for fertilizer uniformity as the response value, the compound experiment of four factors and five levels was carried out. The test factors and horizontal codes are shown in Table 3. The test schemes and results are shown in Table 4.

Quality Evaluation Method for Wide-Boundary Fertilization
The forward direction of the machine is defined as longitudinal, and fertilizer is mainly distributed unevenly in the horizontal direction. Therefore, the CVPU is taken as an index to evaluate the uniformity of fertilizer application. Combined with the static test method of disc fertilizer spreader and the agronomic requirements of wide-boundary fertilization technology, the specific measurement method of CVPU is as follows: a 3000 mm strip is selected in the middle area of the statistical plate, and the Grid Bin Group module of the EDEM post-processing module is applied to divide the statistical plate into grid cells with 10 rows and 10 columns. The size of each grid cell is 300 mm × 30 mm, as shown in Figure 9. Then, the post-processing module of EDEM software was used to measure the quality of fertilizer particles in each cell. The calculation method of the CVPU is shown in Formula (13).
where F is the CVPU, %; w ij is the mass of fertilizer particles in the grid cell in row i and column j, g. w is the average mass of fertilizer particles in each cell, g. m is the number of rows of measurement strips. n is the number of grids for each row of strips. where F is the CVPU, %; wij is the mass of fertilizer particles in the grid cell in row i and column j, g. is the average mass of fertilizer particles in each cell, g. m is the number of rows of measurement strips. n is the number of grids for each row of strips.

Multiple Regression Analysis
The simulation test results based on the scheme are shown in Table 4. The test results are analyzed by Design-Expert software, and the model terms that are not significant are merged into the residual terms. The optimized variance analysis results are shown in Table 5. It can be seen from Table 5 that the term of the model is p < 0.0001, and the term of lack of fit is p = 0.5841 > 0.1, respectively. It is indicated that the term of the model has an extremely significant influence, and the term of lack of fit is less significant; that is, the regression model has a high degree of fit with the simulation results. The quadratic regression model is shown in Formula (14). It can be seen from Table 5 that the terms of 1 X X X X X X X X

Multiple Regression Analysis
The simulation test results based on the scheme are shown in Table 4. The test results are analyzed by Design-Expert software, and the model terms that are not significant are merged into the residual terms. The optimized variance analysis results are shown in Table 5. It can be seen from Table 5 that the term of the model is p < 0.0001, and the term of lack of fit is p = 0.5841 > 0.1, respectively. It is indicated that the term of the model has an extremely significant influence, and the term of lack of fit is less significant; that is, the regression model has a high degree of fit with the simulation results. The quadratic regression model is shown in Formula (14). It can be seen from Table 5 that the terms of X 1 , X 2 , X 3 , X 4 , X 2 1 , X 2 2 , X 2 3 , X 2 4 have an extremely significant influence on the response value (CVPU); the term of X 3 X 4 has a significant influence on the response value; and the term of X 2 X 4 has a generally significant influence on the response value. The primary and secondary order of influence of key factors on CVPU is X 2 > X 1 > X 3 .

Response Surface Analysis
The response surface analysis of interaction terms in Table 5 was carried out by using Design-Expert software, as shown in Figure 10. When the focal length coefficient and tilt angle are taken as the intermediate value (1. 5 and 35 mm), the interaction influence of lateral span and ground clearance on the CVPU is shown in Figure 10a. When ground clearance is determined, with the increase in the lateral span value, CVPU first decreases and then increases; when the lateral span is determined, with the increase in ground clearance, the CVPU first decreases and then increases. When the focal length coefficient and lateral span are taken as the intermediate value (1. 5 and 65 mm), the interaction influence of tilt angle and ground clearance on the CVPU is shown in Figure 10b. When the tilt angle is determined, with the increase in the ground clearance value, the CVPU first decreases and then increases; when ground clearance is determined, with the increase in the tilt angle, the CVPU first decreases and then increases. Note: *** means extremely significant (p < 0.01), ** means significant (0.01 < p < 0.05), * means gener ally significant (0.05 < p < 0.1).

Response Surface Analysis
The response surface analysis of interaction terms in Table 5 was carried out by using Design-Expert software, as shown in Figure 10. When the focal length coefficient and til angle are taken as the intermediate value (1.5 and 35 mm), the interaction influence o lateral span and ground clearance on the CVPU is shown in Figure 10a. When ground clearance is determined, with the increase in the lateral span value, CVPU first decreases and then increases; when the lateral span is determined, with the increase in ground clear ance, the CVPU first decreases and then increases. When the focal length coefficient and lateral span are taken as the intermediate value (1.5 and 65 mm), the interaction influence of tilt angle and ground clearance on the CVPU is shown in Figure 10b. When the tilt angle is determined, with the increase in the ground clearance value, the CVPU first decreases and then increases; when ground clearance is determined, with the increase in the tilt an gle, the CVPU first decreases and then increases.

Parameter Optimization Analysis
To optimize the operational performance of the wide-boundary fertilization device the parameter combinations of focal length coefficient, lateral span, tilt angle, and ground clearance were optimized. The objective function and bundle function were established as shown in Formula 15, respectively, and the mathematical model was analyzed and solved by Design-Expert software. Among a series of parameter combination schemes the one with the highest confidence level was selected as the best parameter-combination scheme. The optimization results of parameter combination are as follows: the focal length coefficient was 1.35, the lateral span was 59.37 mm, the tilt angle was 29.21°, and the ground clearance was 191.56 mm. In this case, the CVPU was 14.11%, and the operation

Parameter Optimization Analysis
To optimize the operational performance of the wide-boundary fertilization device, the parameter combinations of focal length coefficient, lateral span, tilt angle, and ground clearance were optimized. The objective function and bundle function were established, as shown in Formula 15, respectively, and the mathematical model was analyzed and solved by Design-Expert software. Among a series of parameter combination schemes, the one with the highest confidence level was selected as the best parameter-combination scheme. The optimization results of parameter combination are as follows: the focal length coefficient was 1.35, the lateral span was 59.37 mm, the tilt angle was 29.21 • , and the ground clearance was 191.56 mm. In this case, the CVPU was 14.11%, and the operation effect met the relevant standards and requirements.

Field Verification Test
To further verify the performance of the wide-boundary fertilization device under complex field conditions, a field test was conducted at the rice farm of Jiangsu Academy of Agricultural Sciences in October 2021. The experimental field is primitive land with full rice straw; the average height of rice straw is ≥400 mm. The type of soil is heavy, cohesive soil, and the type of fertilizer is compound fertilizer. Using the wide-boundary planter with crushed straw inter-row mulching as the test platform, field verification tests were carried out, as shown in Figure 11. The parameters of the wide-boundary fertilization device are as follows: the focal length coefficient was 1.4, the lateral span was 60 mm, the tilt angle was 30, and the ground clearance was 192 mm. During the test, the height of the rotary tillage device was adjusted upward to keep the rotary plowing knife at a certain height from the ground to avoid the interference of the rotary plowing knife on the fertilizer particles. The forward speed of the machine was set to 1.4 m/s, and the other working parameters of the machine were matched at the same time; that is, the rotational speed of the crushing device was 2000 r/min, and the rotational speed of the rotary tillage device was 300 r/min. A total of 5 repeated experiments were carried out, and the average value was taken.

Field Verification Test
To further verify the performance of the wide-boundary fertilization device u complex field conditions, a field test was conducted at the rice farm of Jiangsu Acad of Agricultural Sciences in October 2021. The experimental field is primitive land with rice straw; the average height of rice straw is ≥400 mm. The type of soil is heavy, cohe soil, and the type of fertilizer is compound fertilizer. Using the wide-boundary pla with crushed straw inter-row mulching as the test platform, field verification tests w carried out, as shown in Figure 11. The parameters of the wide-boundary fertilization vice are as follows: the focal length coefficient was 1.4, the lateral span was 60 mm, th angle was 30, and the ground clearance was 192 mm. During the test, the height o rotary tillage device was adjusted upward to keep the rotary plowing knife at a ce height from the ground to avoid the interference of the rotary plowing knife on the f lizer particles. The forward speed of the machine was set to 1.4 m/s, and the other wor parameters of the machine were matched at the same time; that is, the rotational spee the crushing device was 2000 r/min, and the rotational speed of the rotary tillage de was 300 r/min. A total of 5 repeated experiments were carried out, and the average v was taken. The field test results showed that the CVPU was 17.63%, and the deviation from theoretical value and the field test value was less than 5%. The deviation was ma caused by the poor flatness of the field ground, so the test results were basically consis which indicated the accuracy of the parameter's combination.

Discussions
To solve poor fertilization uniformity at low rotational speeds, Zha et al. [28 signed a blocking wheel-type screw fertilizer distributor. Simulation tests based on ED software were carried out to optimize the distributor variables, and then bench veri tion tests were built under the same conditions as the simulation tests. The results of ulation tests revealed that the minimum coefficient of variation of fertilization unifor was 19.27%. Hwang et al. [13] analyzed the behavior of fertilizer particles on the spi using the discrete element method (DEM). When the optimal combination of param The field test results showed that the CVPU was 17.63%, and the deviation from the theoretical value and the field test value was less than 5%. The deviation was mainly caused by the poor flatness of the field ground, so the test results were basically consistent, which indicated the accuracy of the parameter's combination.

Discussions
To solve poor fertilization uniformity at low rotational speeds, Zha et al. [28] designed a blocking wheel-type screw fertilizer distributor. Simulation tests based on EDEM software were carried out to optimize the distributor variables, and then bench verification tests were built under the same conditions as the simulation tests. The results of simulation tests revealed that the minimum coefficient of variation of fertilization uniformity was 19.27%. Hwang et al. [13] analyzed the behavior of fertilizer particles on the spinner using the discrete element method (DEM). When the optimal combination of parameters was used for the operation, the CV values were <20%. To maximize the fertilizer-utilization rate, Zhao et al. [29] designed side deep fertilizing components of a no-tillage planter. The prototype performance experiment results showed that the sowing width variation coefficient was 28.4%, which could satisfy the side deep fertilizing quality requirements. Fertilizer-application devices were extensively studied. These devices have shown good operational performance in their respective fields, with the CVPU generally being less than 30%. When the device designed in this study was used for fertilization, the field test results showed that the CVPU was 17.63%, which was consistent with the results of other related studies. The results show that it can work well under full rice straw retention.
In the case of fixed values of other factors, the CVPU first decreases and then increases with the increase in the lateral span value. The reasons for this are as follows: when the lateral span value is too small, the size of the fertilizer even-distribution plate is very small, which results in a small probability of the particles bouncing twice; so, fertilizer particles are mainly concentrated in the middle area of the seed belt. When the lateral span value is too large, the particles bounce on the fertilizer even-distribution plate many times, which leads to a larger lateral displacement of the particles; so, the distribution of the particles on the seed belt is less in the middle and more on the sides.
In the case of fixed values of other factors, the CVPU first decreases and then increases with the increase in the ground clearance value. The reasons for this are as follows: when the ground clearance value is too small, the particles' movement time in the air is too short, which leads to a small lateral displacement of particles, so the particles are concentrated in the middle area. When the ground clearance value is too large, the particles' movement time in the air is longer, which leads to a larger width of spreading fertilizer. However, under the action of the crushed straw deflector, the particles are blocked from flying out of the boundary of the seed belt. Therefore, there were more particles in the boundary and fewer particles in the middle of the seed belt.
In the case of fixed values of other factors, the CVPU first decreases and then increases with the increase in the tilt angle value. The reasons for this are as follows: when the tilt angle is too small, the number of collisions between the particles and the plate surface increases, which leads to larger velocity loss and smaller vertical velocity; therefore, the particle movement time increases, and the lateral displacement increases. When the tilt angle is too large, the incidence angle of the collision is too small, which leads to a smaller velocity change and a larger vertical velocity; so the particle motion time is shorter, which leads to insufficient lateral displacement of the particles and aggregation in the middle of the seed belt.
In addition, this study did not consider the influence of factors such as the machine operating speed, variability of fertilizer types, and particle incident speed on the CVPU, which will be further analyzed and tested in subsequent experiments. Different soil characteristics can have some influence on the field test results. On the one hand, the device designed in this study was mainly applied in the middle and lower areas of the Yangtze River plain in China, where the soil properties are more consistent. On the other hand, only representative flat plots were selected in this study to conduct multiple sets of field experiments due to the limitation of research conditions. So, the influence of soils with different characteristics on the CVPU will also be the focus of our next research.

Conclusions
Wide-boundary fertilization technology with crushed straw inter-row mulching can effectively solve the problems, such as a low fertilizer utilization rate, heap soil around a buried fertilizer device, or blocked fertilizing orifice. However, the existing scientific and technological reports are seriously deficient in the basic research and structural parameters of the wide-boundary fertilization equipment, which leads to the lack of theoretical basis in the process of production machineries. Therefore, combined with theoretical analysis and discrete element numerical simulation technology, the wide-boundary fertilization device for wheat wide-boundary sowing was designed and systematically studied. The specific conclusions are as follows: (1) Aiming at the agronomic mode of wide-boundary sowing with crushed straw interrow mulching under full straw croplands in the rice-wheat rotation system, the parameters of the wide-boundary fertilization device were designed. Through the theoretical analysis of the fertilizer movement process, it was concluded that the key structural parameters affecting the performance of the device were the focal length coefficient, lateral span, tilt angle, and ground clearance. (2) The central composite simulation test was conducted on the key structural parameters (focal length coefficient, span, tilt angle, and ground clearance height) with the CVPU as the response value, and the effective laws of factors on the CVPU were analyzed. The results showed that each factor has a significant effect on the CVPU, and its primary and secondary order is tilt angle, lateral span, focal length coefficient, ground clearance, tilt angle × ground clearance, and lateral span × ground clearance. There were certain interactions between tilt angle and ground clearance, lateral span, and ground clearance.