Validation and Calibration of Maize Seed–Soil Inter-Parameters Based on the Discrete Element Method

An appropriate contact mechanics model and parameters are key to achieving accurate results in discrete element analyzis. This is necessary to predict the process of contact collision between the soil and maize seed during deposition. In this paper, the contact process between maize seed and soil is analyzed using the maize seed variety (Liangyu 99) and maize-sowing field soil (with three different moisture contents) as research objects. Based on this, the contact process between maize seeds and soil has been analyzed, on the basis of which a mechanical model suitable for simulating the contact process between maize seeds and soil has been explored, and the selection of parameters between heterogeneous particles (maize seed and soil particles) has been investigated. The results showed that adhesion forces have a significant effect on the collision process between seed and soil particles. While the presence of tangential adhesion force can be replaced by increasing the static and rolling friction coefficients, the normal adhesion force cannot be compensated in this way. The Edinburgh Elasto-Plastic Adhesive (EEPA) model is selected in this paper to describe the contact between seed and soil particles. The significance of the input parameters in the EEPA model is investigated using the Plackett–Burman test. The parameters between soil and seed particles are optimized using the central composite design, and the optimal parameter combinations are obtained. The relative error between the simulation and test result of the slope test for the three soil moisture contents is within 5.4%, validating the accuracy of the calibrated parameters.


Introduction
During the process of sowing maize, the three stages of seeding, dropping, and bedding are all interrelated and affect the quality of sowing [1]. The bedding stage, which is the final stage in the seed-sowing process, plays an important role in ensuring that the seeds are placed into the soil in a uniformly and orderly way, ultimately affecting the uniformity of the final seed distribution in the field [2]. However, during the bedding stage, the maize seeds and soil are prone to bouncing and rolling, especially under highspeed operating conditions [2]. This bouncing and rolling seriously affects the quality of sowing and restricts further increases in the yield of maize per unit area. Therefore, there is a need for in-depth and systematic research on the mechanism of maize seedsoil contact and collision bouncing and its influencing factors under high-speed sowing operations. This research will enable improvements in the quality of sowing, which has important theoretical significance and practical application value. Most of the experimental methods used to study the contact collision process between seed particles and soil are time-consuming, labor-intensive, and restricted by time and season. It is also difficult to obtain information on the force between seed particles and soil, particle displacement, and For the contact between maize seed and soil, further research is needed to determine whether soil adhesion to the seed needs to be taken into account, and whether the presence of adhesion can be mitigated by using methods that increase the static or rolling coefficient of friction.
Agronomy 2023, 13, 2115 3 of 16 In this section, the adhesion of soil to maize seeds during vertical contact of maize seed particles with soil is investigated using a texture instrument (TA.XTC-18, BosinTech, shanghai, China), as shown in Figure 1. The adhesion force is then compared to the gravity of the seed particles. The mass of the seed particles used in the test is weighed using a balance with an accuracy of 0.01 g. Three seeds of similar size are selected and their dimensions are 8.22 mm wide, 12.21 mm high and 4.83 mm thick. If the adhesion force is close to or greater than the gravity of the seed particles, it can be concluded that the adhesion force of the soil to the seed particles will significantly affect the contact collision between the seed particles and the soil.

Testing of Contact Processes between Soil and Maize Seeds
For the contact between maize seed and soil, further research is needed to determine whether soil adhesion to the seed needs to be taken into account, and whether the presence of adhesion can be mitigated by using methods that increase the static or rolling coefficient of friction.
In this section, the adhesion of soil to maize seeds during vertical contact of maize seed particles with soil is investigated using a texture instrument (TA.XTC-18, BosinTech, shanghai, China), as shown in Figure 1. The adhesion force is then compared to the gravity of the seed particles. The mass of the seed particles used in the test is weighed using a balance with an accuracy of 0.01 g. Three seeds of similar size are selected and their dimensions are 8.22 mm wide, 12.21 mm high and 4.83 mm thick. If the adhesion force is close to or greater than the gravity of the seed particles, it can be concluded that the adhesion force of the soil to the seed particles will significantly affect the contact collision between the seed particles and the soil. The test procedure is as follows: First, the seeds are fixed onto the probe of the texture instrument so that the flat side of the seeds is in contact with the soil, and the probe of the texture instrument is adjusted to the appropriate height. Then, the probe of the texture instrument is moved downward at a speed of 1 mm/s until the probe compress the soil by 30% of its deformation, and the probe is moved upward. Finally, the curves of force versus displacement and the data are obtained. Three replicate tests are performed for each set of tests.

Calibration of Parameters between Maize Seed Particles and Soil Particles
The particle model of maize seeds is adopted from the seed model established in previous studies (shown in Figure 2) [10,11]. The mechanical modeling between maize seeds is adopted from the EEPA model, where the constant pull-off force (f0) is 0, slope exponent is 1.5, tensile exponent is 0, surface energy is 0, contact plasticity ratio is 0, and tangential stiff multiplier is 1. The remaining parameters are adopted from the literature. The soil particle model is used as established in previous studies (shown in Figure 3) [20], and the mechanical modeling between soil particles is based on the EEPA model, with the parameters in the model using data from the literature. In this paper, simulations are performed using EDEM 2018 (4.0.0) software. The test procedure is as follows: First, the seeds are fixed onto the probe of the texture instrument so that the flat side of the seeds is in contact with the soil, and the probe of the texture instrument is adjusted to the appropriate height. Then, the probe of the texture instrument is moved downward at a speed of 1 mm/s until the probe compress the soil by 30% of its deformation, and the probe is moved upward. Finally, the curves of force versus displacement and the data are obtained. Three replicate tests are performed for each set of tests.

Calibration of Parameters between Maize Seed Particles and Soil Particles
The particle model of maize seeds is adopted from the seed model established in previous studies (shown in Figure 2) [10,11]. The mechanical modeling between maize seeds is adopted from the EEPA model, where the constant pull-off force (f 0 ) is 0, slope exponent is 1.5, tensile exponent is 0, surface energy is 0, contact plasticity ratio is 0, and tangential stiff multiplier is 1. The remaining parameters are adopted from the literature. The soil particle model is used as established in previous studies (shown in Figure 3) [20], and the mechanical modeling between soil particles is based on the EEPA model, with the parameters in the model using data from the literature. In this paper, simulations are performed using EDEM 2018 (4.0.0) software.
In this paper, the parameters between maize seed particles and soil particles are calibrated by means of angle of repose tests. Firstly, a sensitivity study is carried out on the contact parameters to be artificially input into the EEPA model using the Plackett-Burman test. Then, the sensitivity parameters are calibrated and optimized using the Central Composite Design test.  Particle models of maize seeds of (a) horse-tooth shape, (b) truncated triangular pyramid shape, (c) ellipsoid cone shape, (d) spheroid shape. In this paper, the parameters between maize seed particles and soil particles are calibrated by means of angle of repose tests. Firstly, a sensitivity study is carried out on the contact parameters to be artificially input into the EEPA model using the Plackett-Burman test. Then, the sensitivity parameters are calibrated and optimized using the Central Composite Design test.

EEPA Model
The Edinburgh Elasto-Plastic Adhesive (EEPA) model [21] can capture the shear stresses associated with the stress history and the cohesion of the bulk material. Figure 4 shows a schematic diagram of the EEPA contact spring for the normal direction of f-δ (force-overlap). The normal force fn is equal to the sum of the hysteresis spring force fhys and the normal damping force fnd: (1) where u is the unit normal vector, directed from the point of contact to the centre of the particle; fhys can be calculated from the following equation: Figure 2. Particle models of maize seeds of (a) horse-tooth shape, (b) truncated triangular pyramid shape, (c) ellipsoid cone shape, (d) spheroid shape. Particle models of maize seeds of (a) horse-tooth shape, (b) truncated triangular pyramid shape, (c) ellipsoid cone shape, (d) spheroid shape. In this paper, the parameters between maize seed particles and soil particles are calibrated by means of angle of repose tests. Firstly, a sensitivity study is carried out on the contact parameters to be artificially input into the EEPA model using the Plackett-Burman test. Then, the sensitivity parameters are calibrated and optimized using the Central Composite Design test.

EEPA Model
The Edinburgh Elasto-Plastic Adhesive (EEPA) model [21] can capture the shear stresses associated with the stress history and the cohesion of the bulk material. Figure 4 shows a schematic diagram of the EEPA contact spring for the normal direction of f-δ (force-overlap). The normal force fn is equal to the sum of the hysteresis spring force fhys and the normal damping force fnd: where u is the unit normal vector, directed from the point of contact to the centre of the particle; fhys can be calculated from the following equation:

EEPA Model
The Edinburgh Elasto-Plastic Adhesive (EEPA) model [21] can capture the shear stresses associated with the stress history and the cohesion of the bulk material. Figure 4 shows a schematic diagram of the EEPA contact spring for the normal direction of f-δ (force-overlap).

Figure 2.
Particle models of maize seeds of (a) horse-tooth shape, (b) truncated triangular p shape, (c) ellipsoid cone shape, (d) spheroid shape. In this paper, the parameters between maize seed particles and soil particles a ibrated by means of angle of repose tests. Firstly, a sensitivity study is carried out contact parameters to be artificially input into the EEPA model using the Plackett-B test. Then, the sensitivity parameters are calibrated and optimized using the Centra posite Design test.

EEPA Model
The Edinburgh Elasto-Plastic Adhesive (EEPA) model [21] can capture the stresses associated with the stress history and the cohesion of the bulk material. Fi shows a schematic diagram of the EEPA contact spring for the normal direction (force-overlap). The normal force fn is equal to the sum of the hysteresis spring force fhys and th mal damping force fnd: ( ) where u is the unit normal vector, directed from the point of contact to the centre particle; fhys can be calculated from the following equation: The normal force f n is equal to the sum of the hysteresis spring force f hys and the normal damping force f nd : where u is the unit normal vector, directed from the point of contact to the centre of the particle; f hys can be calculated from the following equation: Agronomy 2023, 13, 2115

of 16
The normal damping force f nd is calculated by the following equation: where υ n is the normal relative velocity; β n is the normal damping factor, which can be calculated from the following equation: where m * is the equivalent mass, m* = (mimj/mi + mj), where mi and mj are the masses of the individual particles; e is the coefficient of restitution.
The tangential force f t is the sum of the tangential elastic force f ts and the tangential damping force f td .
The tangential elasticity fts can be expressed in incremental terms as: f ts = f ts(n−1) + ∆ f ts (6) where f ts(n−1) is the tangential elastic force at the previous step; ∆ f ts is the increment in the tangential elastic force, calculated from the following equation.
where k t is the tangential stiffness coefficient, the ratio of the tangential stiffness coefficient to the normal stiffness coefficient is 1 in EDEM, i.e., the two values are equal; in the luding model, the ratio is 0.2; in LAMMPS and PFC, this value is 2/7; however, the ratio of the tangential stiffness to the normal stiffness of a real elastic material is between 2/3 and 1, subject to the Poisson's ratio; δ t is the tangential displacement increment. The tangential damping force ftd is equal to the product of the tangential damping factor β t and the tangential relative velocity vt, calculated by the following formula.
The tangential damping factor β t is calculated by the following formula.
The ultimate tangential friction is calculated using the Coulomb friction criterion with the normal force modified by the adhesion force: where f ct is the ultimate tangential friction force and µ is the coefficient of friction. In this paper, the EDEM default rolling friction model is used and the total applied torque τ i is calculated by the following equation.
where µ r is the coefficient of rolling friction; R i is the distance from the point of contact to the centre of mass of the particle and w i is the unit angular velocity at the point of contact.

Experimental Setup for Angle of Repose Tests
The angle of repose test for maize seed particles and soil particles was conducted using an electronic universal testing machine. The plexiglass cylinder was lifted upward at a certain speed to complete the accumulation process. To ensure full contact between the seed particles and soil particles and calibrate the contact parameters, maize seed particles were homogeneously mixed with soil particles at a ratio of 4:3. The angle of repose tests were conducted as shown in Figure 5.
where r  is the coefficient of rolling friction; i R is the distance from the poin tact to the centre of mass of the particle and i w is the unit angular velocity at t of contact.

Experimental Setup for Angle of Repose Tests
The angle of repose test for maize seed particles and soil particles was co using an electronic universal testing machine. The plexiglass cylinder was lifted at a certain speed to complete the accumulation process. To ensure full contact the seed particles and soil particles and calibrate the contact parameters, maize s ticles were homogeneously mixed with soil particles at a ratio of 4:3. The angle o tests were conducted as shown in Figure 5. The test procedure is as follows: First, fill the bottom container with soil of a water content and make the surface flat. Next, mix 0.2 kg of soil and 0.15 kg of mai with the same water content in the plexiglass cylinder. Put the cylinder on the so bottom container, and connect the cylinder and the electronic universal testing to the fixed position. Then, the electronic universal testing machine lifts the cylin wards at a speed of 300 mm/min, and the mixed seed particles and soil flow ou cylinder, eventually piling up on the subsoil to form an angle of repose. Finally, us recognition to obtain the size of the angle of repose. Three replicate trials are pe for each set of trials.

Simulation Setup for Angle of Repose Tests
In the simulation of the angle of repose test, the cylinder size and diamete bottom cylinder are the same as in the actual test. The simulation steps are as follo soil particles are generated in the bottom container and the soil surface is levele the bottom soil particles stabilize, a certain mass of soil particles and maize seed are generated in the cylinder, which are stabilized in the cylinder with the suppo bottom soil particles (see Figure 6a). Then, the cylinder moves upwards at a spee mm/min, and soil particles and maize seed particles flow out of the cylinder and late on the underlying soil particles, forming an angle of repose (see Figure 6b). the angle of repose is measured using an image recognition method. Each set of repeated three times. The test procedure is as follows: First, fill the bottom container with soil of a certain water content and make the surface flat. Next, mix 0.2 kg of soil and 0.15 kg of maize seeds with the same water content in the plexiglass cylinder. Put the cylinder on the soil in the bottom container, and connect the cylinder and the electronic universal testing machine to the fixed position. Then, the electronic universal testing machine lifts the cylinder upwards at a speed of 300 mm/min, and the mixed seed particles and soil flow out of the cylinder, eventually piling up on the subsoil to form an angle of repose. Finally, use image recognition to obtain the size of the angle of repose. Three replicate trials are performed for each set of trials.

Simulation Setup for Angle of Repose Tests
In the simulation of the angle of repose test, the cylinder size and diameter of the bottom cylinder are the same as in the actual test. The simulation steps are as follows: first, soil particles are generated in the bottom container and the soil surface is leveled. After the bottom soil particles stabilize, a certain mass of soil particles and maize seed particles are generated in the cylinder, which are stabilized in the cylinder with the support of the bottom soil particles (see Figure 6a). Then, the cylinder moves upwards at a speed of 300 mm/min, and soil particles and maize seed particles flow out of the cylinder and accumulate on the underlying soil particles, forming an angle of repose (see Figure 6b). Finally, the angle of repose is measured using an image recognition method. Each set of trials is repeated three times.

Verification Tests
In this section, the parameters between the calibrated soil and seed are validated the inclined slide test, which is a good basis for the following simulation test of maize se bouncing by touching the soil.

Inclined Slide Test Setup
A sketch of the setup for the inclined slide test of soil and seed particles is presen in Figure 7. The test setup consists of a high-speed video camera, a computer, an in nometer, and an angle meter. The steps of the inclined slide test are as follows: first, the square slot of the inclinometer with soil of a certain water content and make the surf soil level. Then, place the maize seed particles statically on the surface of the soil and pl the inclinometer horizontally, as shown in Figure 7. Next, slowly lift the inclinometer un the maize seed grain slides off the soil surface. Finally, use a high-speed camera to obt an indication  of the angle meter at the time when the sliding of the maize seed partic on the inclined plane first occurs. Repeat each experiment three times.

Verification Tests
In this section, the parameters between the calibrated soil and seed are validated by the inclined slide test, which is a good basis for the following simulation test of maize seed bouncing by touching the soil.

Inclined Slide Test Setup
A sketch of the setup for the inclined slide test of soil and seed particles is presented in Figure 7. The test setup consists of a high-speed video camera, a computer, an inclinometer, and an angle meter. The steps of the inclined slide test are as follows: first, fill the square slot of the inclinometer with soil of a certain water content and make the surface soil level. Then, place the maize seed particles statically on the surface of the soil and place the inclinometer horizontally, as shown in Figure 7. Next, slowly lift the inclinometer until the maize seed grain slides off the soil surface. Finally, use a high-speed camera to obtain an indication ϕ of the angle meter at the time when the sliding of the maize seed particles on the inclined plane first occurs. Repeat each experiment three times.
Agronomy 2023, 13, x FOR PEER REVIEW 7 Figure 6. Simulation setup for the angle of repose test between maize seed particles and soil cles: (a) cylinder erected on soil, (b) angle of repose formation.

Verification Tests
In this section, the parameters between the calibrated soil and seed are validate the inclined slide test, which is a good basis for the following simulation test of maize bouncing by touching the soil.

Inclined Slide Test Setup
A sketch of the setup for the inclined slide test of soil and seed particles is prese in Figure 7. The test setup consists of a high-speed video camera, a computer, an nometer, and an angle meter. The steps of the inclined slide test are as follows: firs the square slot of the inclinometer with soil of a certain water content and make the su soil level. Then, place the maize seed particles statically on the surface of the soil and p the inclinometer horizontally, as shown in Figure 7. Next, slowly lift the inclinometer the maize seed grain slides off the soil surface. Finally, use a high-speed camera to o an indication  of the angle meter at the time when the sliding of the maize seed par on the inclined plane first occurs. Repeat each experiment three times. The simulation of inclined sliding uses the calibrated parameters, while the othe rameters remain consistent with those used in the actual test. The simulation setu

Simulation Setup for Inclined Slide Test
The simulation of inclined sliding uses the calibrated parameters, while the other parameters remain consistent with those used in the actual test. The simulation setup for inclined sliding is shown in Figure 8. The simulation steps are as follows: first, a seed particle is generated on the soil surface with a frontal downward attitude; then, the soil layer is rotated anticlockwise at 10 deg/s until the seed particle slip off the soil. Finally, the angle of inclination of the seed particles as they slip off the soil is calculated. Each set of simulations is repeated three times. inclined sliding is shown in Figure 8. The simulation steps are as follows: first, a seed particle is generated on the soil surface with a frontal downward attitude; then, the soi layer is rotated anticlockwise at 10 deg/s until the seed particle slip off the soil. Finally, th angle of inclination of the seed particles as they slip off the soil is calculated. Each set o simulations is repeated three times.

Analyzsis of Contact Modelling between Soil and Maize Seeds
The force-versus-displacement curves for three moisture contents of soil in norma contact with seeds are shown in Figure 9 (where gf is the unit of force equal to th magnitude of the force of gravity on an object with a mass of 1 g). It should be noted tha the horizontal axis of the graph represents the change in probe displacement, while th vertical axis represents the magnitude of the contact force between the maize seed and th soil. The maximum negative force in the graph represents the maximum adhesion of th soil to the seed.
The experimental data are tallied, and the mass of the seed particles used in the tes is weighed using a balance with an accuracy of 0.01 g, as shown in Figure 10. The mass o the seed particles is 0.31 g, and the weight of the seed particles is calculated to be 0.003 N. The maximum adhesion force of the soil in normal contact with the seed at thre moisture contents is statistically analyzed, as shown in Table 1 below. From the table, i can be concluded that the maximum adhesion of soil to maize seed particles is 0.099993 N when the soil moisture content is 15%. This can be compared to the gravitational force o seed particles, which is 32.26 times the gravitational force of the seed particles. When th soil moisture content is 20%, the maximum adhesion of soil to maize seed particles i

Analyzsis of Contact Modelling between Soil and Maize Seeds
The force-versus-displacement curves for three moisture contents of soil in normal contact with seeds are shown in Figure 9 (where gf is the unit of force equal to the magnitude of the force of gravity on an object with a mass of 1 g). It should be noted that the horizontal axis of the graph represents the change in probe displacement, while the vertical axis represents the magnitude of the contact force between the maize seed and the soil. The maximum negative force in the graph represents the maximum adhesion of the soil to the seed.
inclined sliding is shown in Figure 8. The simulation steps are as follows: first, a se particle is generated on the soil surface with a frontal downward attitude; then, the s layer is rotated anticlockwise at 10 deg/s until the seed particle slip off the soil. Finally, t angle of inclination of the seed particles as they slip off the soil is calculated. Each set simulations is repeated three times.

Analyzsis of Contact Modelling between Soil and Maize Seeds
The force-versus-displacement curves for three moisture contents of soil in norm contact with seeds are shown in Figure 9 (where gf is the unit of force equal to the mag tude of the force of gravity on an object with a mass of 1 g). It should be noted that t horizontal axis of the graph represents the change in probe displacement, while the ver cal axis represents the magnitude of the contact force between the maize seed and the so The maximum negative force in the graph represents the maximum adhesion of the s to the seed.
The experimental data are tallied, and the mass of the seed particles used in the t is weighed using a balance with an accuracy of 0.01 g, as shown in Figure 10. The mass the seed particles is 0.31 g, and the weight of the seed particles is calculated to be 0.00 N. The maximum adhesion force of the soil in normal contact with the seed at thr moisture contents is statistically analyzed, as shown in Table 1 below. From the table can be concluded that the maximum adhesion of soil to maize seed particles is 0.099993 when the soil moisture content is 15%. This can be compared to the gravitational force seed particles, which is 32.26 times the gravitational force of the seed particles. When t soil moisture content is 20%, the maximum adhesion of soil to maize seed particles 0.158897 N, which is compared to the gravity of seed particles, and the adhesion force 51.26 times the gravity of seed particles. When the soil moisture content is 25%, the ma The experimental data are tallied, and the mass of the seed particles used in the test is weighed using a balance with an accuracy of 0.01 g, as shown in Figure 10. The mass of the seed particles is 0.31 g, and the weight of the seed particles is calculated to be 0.0031 N. In addition, for example, the curve of the normal contact force and displacemen tween the soil and seed with 25% water content, as well as the loading process, unloa process, and adhesion process of the force of this curve, are similar to the curve i nonlinear model of normal contact of the EEPA model (f0 = 0), as shown in Figure 10  In summary, the EEPA model was chosen as the contact model between seed par and soil particles in this paper.

Calibration of Parameters between Maize Seed Particles and Soil Particles
The results of the angle of repose test between maize seed particles and soil par are shown in Table 2. From the data in the table, it is worth noting that the angle of re does not increase with increasing soil moisture content. For example, the angle of re of soil mixed with seed particles at 25% moisture content is less than the angle of re of soil mixed with seed particles at 20% moisture content. This is due to the fact tha soil with 25% moisture content mixed with the seed particles has a higher adhesio tween the soil and the seed particles. This adhesion prevents the seed particles from tering on the subsoil at the beginning stage of the ascent of the cylinder; instead, con dation occurs, as shown in Figure 11. As the cylinder continues to rise, the mixed soi seed particles collapse, resulting in the formation of a smaller angle of repose. The maximum adhesion force of the soil in normal contact with the seed at three moisture contents is statistically analyzed, as shown in Table 1 below. From the table, it can be concluded that the maximum adhesion of soil to maize seed particles is 0.099993 N when the soil moisture content is 15%. This can be compared to the gravitational force of seed particles, which is 32.26 times the gravitational force of the seed particles. When the soil moisture content is 20%, the maximum adhesion of soil to maize seed particles is 0.158897 N, which is compared to the gravity of seed particles, and the adhesion force is 51.26 times the gravity of seed particles. When the soil moisture content is 25%, the maximum adhesion of soil to maize seed particles is 0.56029 N, which is compared to the gravity of seed particles, and the adhesion force is 180.76 times the gravity of seed particles. Therefore, it is necessary to consider the presence of adhesion forces for the contact of seed particles with soil, which can significantly affect the process of seed-soil contact collisions. In addition, the presence of tangential adhesion can be replaced by increasing the coefficient of static friction and the coefficient of rolling friction, but the adhesion in the normal direction cannot be compensated in this way. In addition, for example, the curve of the normal contact force and displacement between the soil and seed with 25% water content, as well as the loading process, unloading process, and adhesion process of the force of this curve, are similar to the curve in the nonlinear model of normal contact of the EEPA model (f 0 = 0), as shown in Figure 10.
In summary, the EEPA model was chosen as the contact model between seed particles and soil particles in this paper.

Calibration of Parameters between Maize Seed Particles and Soil Particles
The results of the angle of repose test between maize seed particles and soil particles are shown in Table 2. From the data in the table, it is worth noting that the angle of repose does not increase with increasing soil moisture content. For example, the angle of repose of soil mixed with seed particles at 25% moisture content is less than the angle of repose of soil mixed with seed particles at 20% moisture content. This is due to the fact that the soil with 25% moisture content mixed with the seed particles has a higher adhesion between the soil and the seed particles. This adhesion prevents the seed particles from scattering on the subsoil at the beginning stage of the ascent of the cylinder; instead, consolidation occurs, as shown in Figure 11. As the cylinder continues to rise, the mixed soil and seed particles collapse, resulting in the formation of a smaller angle of repose.   Figure 11. Consolidation between seeds and soil with 25% water content.

Plackett-Burman Test
In the EEPA model, the contact parameters to be determined between ma particles and soil particles include restitution coefficient, coefficient of static fricti ficient of rolling friction, constant pull-off force, surface energy, contact plastici tensile exp, slope exp and tangential stiff multiplier. The constant pull-off force is 0, the slope exp is taken as 1.5, and the ranges of the other seven parameters are s Table 3. Table 3. Parameter-level settings for the Plackett-Burman test.

Symbolic
Parameter Low Level (−1) High Lev Figure 11. Consolidation between seeds and soil with 25% water content.

Plackett-Burman Test
In the EEPA model, the contact parameters to be determined between maize seed particles and soil particles include restitution coefficient, coefficient of static friction, coefficient of rolling friction, constant pull-off force, surface energy, contact plasticity ratio, tensile exp, slope exp and tangential stiff multiplier. The constant pull-off force is taken as 0, the slope exp is taken as 1.5, and the ranges of the other seven parameters are shown in Table 3. Contact plasticity ratio 0.4 0.7 X 6 Tensile exp 1 4 X 7 Tangential stiff multiplier 0.5 0.9 X 8 -X 11 Virtual parameters -- The design scheme of the Plackett-Burman test and the results of the test are shown in Table 4. Table 4. Plackett-Burman test results for angle of repose tests.

Sequences
X 1 X 2 X 3 X 4 X 5 X 6 X 7 X 8 X 9 X 10 X 11 Angle of Repose, deg The test results are then analyzsed by ANOVA to determine the degree of influence of each parameter on the angle of repose, as shown in Table 5. Note: ** indicates that the item is highly significant (p < 0.01) and * indicates that the item is significant (p < 0.05).
From the results of the Plackett-Burman test, it can be seen that the coefficient of static friction and surface energy (p < 0.01) between soil and seed particles has a highly significant effect on the angle of repose, while the contact plasticity ratio between soil and seed particles has a significant effect on the angle of repose (p < 0.05). For the other parameters, the restitution coefficient is assumed to be 0.7, the coefficient of rolling friction to be 0.05, the tensile exp to be 3, and the tangential stiff multiplier to be 0.67.

Central Composite Design
The sensitive parameters (coefficient of static friction, surface energy, and contact plasticity ratio) are calibrated and optimized based on Plackett-Burman tests. The high and low levels of the coefficients of static friction between soil and seed particles are 0.05 and 0.2, respectively. The high and low levels of surface energy between soil and seed particles are 1 and 5, respectively. The high and low levels of contact plasticity ratio between soil and seed particles are 0.4 and 0.7, respectively. Table 6 shows a coded table of factor levels for the Central Composite Design trial. The design solutions and simulation results of Central Composite Design test are shown in Table 7. From the experimental results obtained in the table, the following binary multiple regression equations are derived using the coefficient of static friction (X 2 ), surface energy (X 4 ), and contact plasticity ratio (X 5 ) as independent variables, and the angle of repose (Y) as the response value: The results of the tests in the table are also analyzed by ANOVA, and the results are shown in Table 8. Note: ** indicates that the item is highly significant (p < 0.01) and * indicates that the item is significant (p < 0.05).
As can be seen from the data in the table, the regression model has a p value of <0.0001, indicating that the difference reaches a highly significant level. Additionally, the level of the misfit term is not significant (p > 0.05), indicating a good fit of the equation. These findings suggest that it is feasible to apply a modified mathematical model to characterize the extent of influence of the factors on the response values.
In addition, as shown in the table, X 2 (coefficient of static friction) and X 4 (surface energy) have a highly significant effect (p < 0.01) on the angle of repose, while X 5 (contact plasticity ratio) has no significant effect (p > 0.05) on the angle of repose. Among the secondary terms, X 2 2 has a significant effect (p < 0.05) on the angle of repose. None of the interaction terms have a significant interaction effect (p > 0.05).
The response surface plots among the factors are obtained based on the quadratic regression equation, which reveals the interaction between the factors and the optimal parameter combinations. From Figure 12, it can be observed that the optimal parameters exist within the designed range of factor levels. Note: ** indicates that the item is highly significant (p < 0.01) and * indicates that the item is significant (p < 0.05).
As can be seen from the data in the table, the regression model has a p value of < 0.0001, indicating that the difference reaches a highly significant level. Additionally, the level of the misfit term is not significant (p > 0.05), indicating a good fit of the equation. These findings suggest that it is feasible to apply a modified mathematical model to characterize the extent of influence of the factors on the response values.
In addition, as shown in the table, X2 (coefficient of static friction) and X4 (surface energy) have a highly significant effect (p < 0.01) on the angle of repose, while X5 (contact plasticity ratio) has no significant effect (p > 0.05) on the angle of repose. Among the secondary terms, X22 has a significant effect (p < 0.05) on the angle of repose. None of the interaction terms have a significant interaction effect (p > 0.05).
The response surface plots among the factors are obtained based on the quadratic regression equation, which reveals the interaction between the factors and the optimal parameter combinations. From Figure 12, it can be observed that the optimal parameters exist within the designed range of factor levels. Optimization of the established mathematical model by means of statistical analyzis software resulted in a set of optimal parameters for soil with a moisture content of 15%. These parameters include a static coefficient of friction between soil and seed particles of 0.17, a surface energy of 0.97, and a contact plasticity ratio of 0.48; the optimal parameters for a soil with 20% moisture content are: a static friction coefficient of 0.19 between soil and seed particles, a surface energy of 1.12, and a contact plasticity ratio of 0.63; for a soil Optimization of the established mathematical model by means of statistical analyzis software resulted in a set of optimal parameters for soil with a moisture content of 15%. These parameters include a static coefficient of friction between soil and seed particles of 0.17, a surface energy of 0.97, and a contact plasticity ratio of 0.48; the optimal parameters for a soil with 20% moisture content are: a static friction coefficient of 0.19 between soil and seed particles, a surface energy of 1.12, and a contact plasticity ratio of 0.63; for a soil with 25% moisture content, the optimal parameters are: a static friction coefficient of 0.1 between soil and seed particles, a surface energy of 2.7, and a contact plasticity ratio of 0.41.

Results of the Validation Test
A comparison of the test and simulation results for the inclined slip test is shown in Figure 13. For soils with 15%, 20%, and 25% moisture content, the relative error between the simulation results and the actual test results of the slope test are 4.3%, 0.8%, and 5.4%, respectively. This proves the accuracy of the parameters obtained from the calibration.
with 25% moisture content, the optimal parameters are: a static friction coefficien between soil and seed particles, a surface energy of 2.7, and a contact plasticity 0.41.

Results of the Validation Test
A comparison of the test and simulation results for the inclined slip test is sh Figure 13. For soils with 15%, 20%, and 25% moisture content, the relative error b the simulation results and the actual test results of the slope test are 4.3%, 0.8%, an respectively. This proves the accuracy of the parameters obtained from the calibra Figure 13. A comparison of the test and simulation results for the inclined slip test.

Conclusions
This paper analyzes the contact process between maize seeds and soil. Based analysis, it explores a mechanical model that is suitable for simulating the contact between maize seeds and soil. Additionally, it investigates the selection method of eters between heterogeneous particles. The following conclusions are drawn: (1) In this paper, a contact mechanics model suitable for modeling the contact b heterogeneous particles (maize seeds and soil particles) is investigated. achieved by studying the magnitude of the adhesion force of the soil on th seeds when the maize seed particles are in vertical contact with the soil. The a force is then compared with the gravity force of the seed particles. The res follows: when the moisture content of soil particles is between 15 and 25%, th cal adhesion between the seed and soil is 32-180 times the weight of the seed fore, when considering the contact between seed particles and soil, it is nece take into account the presence of adhesion, which can significantly affect the c process between the seed and soil. In addition, the presence of tangential a can be reduced by increasing the coefficients of static and rolling friction, but hesion in the normal direction cannot be compensated for in this way. Moreo loading process, unloading process, and adhesion process of the normal conta and displacement curves between the soil and seed are similar to the curve nonlinear model of normal contact (f0 = 0) of the EEPA model. In summary, th model was chosen as the contact model between seed particles and soil par this paper. (2) For the test on the angle of repose between maize seed particles and soil par is found that the angle of repose does not increase with an increase in soil m content, but instead increases and then decreases. For example, the angle o of soil mixed with seed particles with 25% moisture content is smaller than th of repose of soil mixed with seed particles with 20% moisture content. This i the fact that the soil with 25% moisture content mixed with seed particles has

Conclusions
This paper analyzes the contact process between maize seeds and soil. Based on this analysis, it explores a mechanical model that is suitable for simulating the contact process between maize seeds and soil. Additionally, it investigates the selection method of parameters between heterogeneous particles. The following conclusions are drawn: (1) In this paper, a contact mechanics model suitable for modeling the contact between heterogeneous particles (maize seeds and soil particles) is investigated. This is achieved by studying the magnitude of the adhesion force of the soil on the maize seeds when the maize seed particles are in vertical contact with the soil. The adhesion force is then compared with the gravity force of the seed particles. The result is as follows: when the moisture content of soil particles is between 15 and 25%, the vertical adhesion between the seed and soil is 32-180 times the weight of the seed. Therefore, when considering the contact between seed particles and soil, it is necessary to take into account the presence of adhesion, which can significantly affect the collision process between the seed and soil. In addition, the presence of tangential adhesion can be reduced by increasing the coefficients of static and rolling friction, but the adhesion in the normal direction cannot be compensated for in this way. Moreover, the loading process, unloading process, and adhesion process of the normal contact force and displacement curves between the soil and seed are similar to the curves of the nonlinear model of normal contact (f 0 = 0) of the EEPA model. In summary, the EEPA model was chosen as the contact model between seed particles and soil particles in this paper. (2) For the test on the angle of repose between maize seed particles and soil particles, it is found that the angle of repose does not increase with an increase in soil moisture content, but instead increases and then decreases. For example, the angle of repose of soil mixed with seed particles with 25% moisture content is smaller than the angle of repose of soil mixed with seed particles with 20% moisture content. This is due to the fact that the soil with 25% moisture content mixed with seed particles has a higher adhesion between the soil and the seed particles. During the ascent of the cylinder, it does not scatter on the subsoil at the beginning stage but instead is cemented together. As the cylinder continues to rise, the mixed soil and the seed particles collapse, resulting in the formation of a smaller angle of repose. (3) A significance study of the input parameters in the EEPA model using the Plackett-Burman test shows that the coefficient of static friction and surface energy between the maize seed and the soil particles has a highly significant effect on the angle of repose test. Additionally, the contact plasticity ratio has a significant effect on the angle of repose test. The parameters between soil particles and seed particles are optimized using Central Composite Design, and the optimum combination of parameters is obtained. For a soil with 15% moisture content, the optimal parameters are as follows: a static coefficient of friction between soil and seed particles of 0.17, a surface energy of 0.97, and a contact plasticity ratio of 0.48. The optimal parameters for a soil with 20% moisture content are: a static friction coefficient of 0.19 between soil and seed particles, a surface energy of 1.12, and a contact plasticity ratio of 0.63. For a soil with 25% moisture content, the optimal parameters are: a static friction coefficient of 0.1 between soil and seed particles, a surface energy of 2.7, and a contact plasticity ratio of 0.41.