Simulation and Experiment of the Spiral Digging End-Effector for Hole Digging in Plug Tray Seedling Substrate

: Usually the planters used in agricultural machinery face two problems: The planters have to spread the soil twice using the tray seed planter, and it is difﬁcult for them to dig the holes before planting. This study has designed a spiral digging end-effector to dig holes in parallel effectively and quickly so that the seeds can be planted in the holes easily. A spiral digging end-effector was designed with ﬁve spiral bits, a synchronous belt, a cylinder, a gear motor, and a connecting plate based on the tray size and the pumpkin seed characteristics. Before the optimization of the end-effector’s parameters, the substrate discrete element model parameters were calibrated using the method of “material funneling” applied for the EDEM discrete element model simulation. The contact parameters and model parameters (substrate–substrate static friction coefﬁcient, substrate–substrate rolling friction coefﬁcient, and surface energy) that have signiﬁcant impacts on the substrate AOR (angle of repose) were selected by applying the Plackett–Burman Design and the path of steepest ascent method, and their respective optimal value range were determined. The optimal parameter values were obtained through the Central Composite Design response surface analysis test, and the parameters of the spiral digging end-effector were optimized combining the substrate particle simulation. The veriﬁcation test results indicate that the substrate discrete element model parameters were accurate and reliable if the substrate static friction coefﬁcient was 0.427; the substrate rolling friction coefﬁcient was 0.039, and the surface energy was 0.228. When the cone angle is 30 ◦ , the spiral angle is 80 ◦ , and the rotational speed is 240 r/min, the section width at the hole depth of 13 mm is 23.1 mm, and the particle overﬂow proportion is about 2.05%. Compared with the protrusions penetrating end-effector, soil porosity is increased, and soil aeration is improved. Therefore, the spiral digging end-effector can effectively and quickly dig holes in the seedling substrate and can sow pumpkin seeds, which provides a research basis for the design and improvement of automatic pumpkin seed sowing equipment.


Introduction
The development of agriculture has solved the great problem of insufficient vegetable supply in China from a long time ago, and the green leaf product industry is able to provide sufficient stock for the people all year round. Seedling cultivation is the first and most important link in agricultural production. It is also an important guarantee for the high yield of vegetables and their quality. The precision seed planter with tray is the main equipment used for seedling and production. They can lower the labor intensity of

The Working Principle of the Spiral Digging End-Effector
The pumpkin seeds are usually planted 10 mm~15 mm below the surface of the soil, and a hole needs to be dug in the center of the plug tray substrate for sowing. To avoid the compaction of the substrate during digging and affecting the germination and growth of seeds, the spiral digging end-effector was designed, as shown in Figure 2. The endeffector is composed of the following parts: five spiral bits, a synchronous belt, a cylinder, a gear motor, and a connecting plate. An MGPM20-50 dual-axial cylinder (SMC Corporation, Tokyo, Japan) was connected to the connecting plate; the five spiral bits were fixed on the connecting plate through the bearing with 50 mm of space between each other. Since the electric motor will drive the five spiral bits to rotate simultaneously, and the driving ratio of electric motor gears and the spiral blade gears is 1:6, a high rotation torque of the electric motor is required. It was found during the trial test that the faster the spiral bits rotated, the easier the substrate would lose and fall apart. That's why the ZGA37RG electric motor (Zhejiang Zhengke Electric Motor Co., Ltd. Yueqing, China) equipped with a ZYTD-32-C deceleration box was selected. It had a loading rotational speed of 60 r/min, and the loading moment was roughly 1 N·m. The electric gear motor was fixed on the connecting plate, and the driving gear would drive one spiral blade to rotate, and the remaining four spiral bits would be driven by the synchronous gear belt.

The Working Principle of the Spiral Digging End-Effector
The pumpkin seeds are usually planted 10 mm~15 mm below the surface of the soil, and a hole needs to be dug in the center of the plug tray substrate for sowing. To avoid the compaction of the substrate during digging and affecting the germination and growth of seeds, the spiral digging end-effector was designed, as shown in Figure 2. The end-effector is composed of the following parts: five spiral bits, a synchronous belt, a cylinder, a gear motor, and a connecting plate. An MGPM20-50 dual-axial cylinder (SMC Corporation, Tokyo, Japan) was connected to the connecting plate; the five spiral bits were fixed on the connecting plate through the bearing with 50 mm of space between each other. Since the electric motor will drive the five spiral bits to rotate simultaneously, and the driving ratio of electric motor gears and the spiral blade gears is 1:6, a high rotation torque of the electric motor is required. It was found during the trial test that the faster the spiral bits rotated, the easier the substrate would lose and fall apart. That's why the ZGA37RG electric motor (Zhejiang Zhengke Electric Motor Co., Ltd. Yueqing, China) equipped with a ZYTD-32-C deceleration box was selected. It had a loading rotational speed of 60 r/min, and the loading moment was roughly 1 N·m. The electric gear motor was fixed on the connecting plate, and the driving gear would drive one spiral blade to rotate, and the remaining four spiral bits would be driven by the synchronous gear belt. During digging, the conveyor belt would transport the tray right und digging end-effector. At this moment, the gear motor would drive the five rotate ( Figure 2I), and at the same time, the end-effector would move do pushing the cylinder ( Figure 2Ⅱ). The spiral bits would dig into the substr rotating for several seconds. Finally, the cylinder was retracted, and a hole w

Spiral Bit Design
The spiral bit includes rotating shaft, drill tip, and spiral blade. Durin operation, the spiral bit rotated around the rotating shaft. The substrate pa the rotating shaft were transported to the spiral blade by the drill tip. After strate particles slid upward along the spiral blade under the cutting and ro Thereby, the substrate was lifted, and the hole was formed.

Drill Tip Design
In order to prevent the rotating shaft from compacting the substrate, an time centering and transporting the substrate particles laterally, the drill tip end of the rotating shaft. Drill tips include four types: flat triangle, twisted w small conical spiral, and fork. Among them, the small conical spiral drill tip h cutting resistance and working stability and has a lifting effect on the partic other three exert downward pressure on the substrate particles. Therefore, t small conical spiral was adopted in this paper.
According to the diameter of the rotating shaft, the bottom radius of the designed to be 2.5 mm, and the height was 5.0 mm. In order to improv transport capacity of the drill tip on the substrate particles, 4 spiral blades During digging, the conveyor belt would transport the tray right under the spiral digging end-effector. At this moment, the gear motor would drive the five spiral bits to rotate ( Figure 2I), and at the same time, the end-effector would move downwards by pushing the cylinder ( Figure 2II). The spiral bits would dig into the substrate and keep rotating for several seconds. Finally, the cylinder was retracted, and a hole was formed.

Spiral Bit Design
The spiral bit includes rotating shaft, drill tip, and spiral blade. During the digging operation, the spiral bit rotated around the rotating shaft. The substrate particles below the rotating shaft were transported to the spiral blade by the drill tip. After that, the substrate particles slid upward along the spiral blade under the cutting and rotating action. Thereby, the substrate was lifted, and the hole was formed.

Drill Tip Design
In order to prevent the rotating shaft from compacting the substrate, and at the same time centering and transporting the substrate particles laterally, the drill tip was set at the end of the rotating shaft. Drill tips include four types: flat triangle, twisted wing triangle, small conical spiral, and fork. Among them, the small conical spiral drill tip has better soil cutting resistance and working stability and has a lifting effect on the particles, while the other three exert downward pressure on the substrate particles. Therefore, the drill tip of small conical spiral was adopted in this paper.
According to the diameter of the rotating shaft, the bottom radius of the drill tip was designed to be 2.5 mm, and the height was 5.0 mm. In order to improve the lateral transport capacity of the drill tip on the substrate particles, 4 spiral blades were evenly distributed and rotated 1/4 circle, and the blades protruded 2.5 mm from the drill tip. The main dimensions of the drill tip are shown in Figure 3.

Spiral Blade Design
Since the shape of sowing holes was the frustum of a cone, the lapsed easily after the cylindrical spiral blade dug the holes, and the sta shape was poor. The screw pitch of the spiral blades increased continu to top, which was beneficial to avoid the extrusion and blockage of the the blades [13,14]. The single-line spiral blade was not easily blocked fore, the single-line variable pitch conical spiral blade was designed [1 the spiral blade was a conical helix with constant inclination and varia and its parameter equations were Equation (1) and the polar coordina where n = a sinγ, b = a cosγ, m = sinγ cosα, k = a tanα, a is a constant, γ i t is an independent variable, and α is the spiral angle (°).

Spiral Blade Design
Since the shape of sowing holes was the frustum of a cone, the loose substrate collapsed easily after the cylindrical spiral blade dug the holes, and the stability of the groove shape was poor. The screw pitch of the spiral blades increased continuously from bottom to top, which was beneficial to avoid the extrusion and blockage of the substrate between the blades [13,14]. The single-line spiral blade was not easily blocked by the soil. Therefore, the single-line variable pitch conical spiral blade was designed [15]. The helix line of the spiral blade was a conical helix with constant inclination and variable pitch (Figure 4), and its parameter equations were Equation (1) and the polar coordinate Equation (2).
where n = a sinγ, b = a cosγ, m = sinγ cosα, k = a tanα, a is a constant, γ is the cone angle ( • ), t is an independent variable, and α is the spiral angle ( • ).

Spiral Blade Design
Since the shape of sowing holes was the frustum of a cone lapsed easily after the cylindrical spiral blade dug the holes, and shape was poor. The screw pitch of the spiral blades increased co to top, which was beneficial to avoid the extrusion and blockage the blades [13,14]. The single-line spiral blade was not easily bl fore, the single-line variable pitch conical spiral blade was design the spiral blade was a conical helix with constant inclination and and its parameter equations were Equation (1) where n = a sinγ, b = a cosγ, m = sinγ cosα, k = a tanα, a is a consta t is an independent variable, and α is the spiral angle (°). Since the sowing holes were the frustum of a cone, the upp Since the sowing holes were the frustum of a cone, the upper diameter D should be smaller than the hole diameter of the plug, and the lower diameter d should be larger than the long axis of the seed. The sowing depth H 0 was about 10~15 mm. So, D = 45 mm, d = 22 mm, and H = 15 mm. According to the geometric relationship of the envelope of the conic curve (Equation (3)), the maximum γ was about 37.5 • .

Seedling Substrate Particle Discrete Element Model Parameter Calibration
The seedling substrate particle discrete element model parameters were needed for the simulation analysis and optimization of the spiral digging end-effector before the discrete element simulation. However, no particle model parameters available for the test were found in the literature. Hence, a combinational method of simulation test and the physical test was adopted for the first time to calibrate the parameters of the substrate simulation model.

Test Method
The "Material funneling" method was adopted in the physical test, and the physical AOR was obtained by measuring the substrate particle pile. Then the EDEM software was used for the simulation test. In this procedure, the substrate discrete model parameters were filtered using the Plackett-Burman Design to obtain the parameters that had a significant impact on the AOR. The path of steepest ascent method was carried out to determine the optimal value range of the parameters. The Central Composite Design (CCD) response surface method (RSM) was used to establish the regression model involving the substrate AOR and the significant parameters. The optimal AOR simulation parameters were obtained by solving the regression equation, the target value of which was the physical AOR.

Physical Test on the Substrate AOR
The substrate used in this study was produced by Shaanxi Qili Agriculture Co., Ltd., and the moisture content was detected at (30.55 ± 0.22%) using the halogen moisture tester (HE53, Mettler-Toledo Company, Zurich, Switzerland). A mixture of substrate and water was created by putting 100 g substrate in the flask and adding V mL distilled water to it. The substrate density of (1.38 ± 0.04 g/cm 3 ) was calculated according to Equation (4). The moisture content and bulk density were measured three times.
A total of 500 substrate particles were measured with vernier caliper with digital display; the particle diameters conform with positive distribution, as shown in Figure 5. Figure 6 illustrates the device components and the size of each component. About 50 g substrate was poured through the funnel at a constant speed for about 5 s to form the particle pile on the steel plate. The actual measurement of AOR was conducted using the image processing method. The open-source software OpenCV 3.4 based on Python was used for the particle pile image processing. First, the redundant background of the image was eliminated to obtain the original image ( Figure 7a). After that, grayscale ( Figure 7b) and binarization methods ( Figure 7c) were used to process the original image (Figure 7a), followed by the Canny edge detection algorithm to extract the boundary profile ( Figure 7d) and then save the file. Last, Digitizer of Origin 2019b software was used to import the processed boundary profile image of the particle pile and obtain its pixel spot coordinates after setting the image pixel and coordinate axis ( Figure 7e). The non-linear fitting was used to obtain the Gaussian distribution function. Please refer to the Gaussian function deduction process to obtain the AOR of the substrate particles [16,17]. The test was repeated 20 times for the average value. The ultimate actual substrate AOR was (39.03 ± 3.29 • ).
The measurement error was ±3.29 • , which fitted the purpose of this research. The same method was applied for the measurement of the substrate simulation test AOR.
A total of 500 substrate particles were measured w play; the particle diameters conform with positive dist Figure 5. Particle diameter distribution. Figure 6 illustrates the device components and the g substrate was poured through the funnel at a consta particle pile on the steel plate. The actual measuremen image processing method. The open-source software used for the particle pile image processing. First, the re was eliminated to obtain the original image (Figure 7a  Figure 7c) were used to process the original followed by the Canny edge detection algorithm to extract the bound 7d) and then save the file. Last, Digitizer of Origin 2019b software was processed boundary profile image of the particle pile and obtain its pix after setting the image pixel and coordinate axis ( Figure 7e). The non used to obtain the Gaussian distribution function. Please refer to the deduction process to obtain the AOR of the substrate particles [16,17 peated 20 times for the average value. The ultimate actual substrate A 3.29°). The measurement error was ±3.29°, which fitted the purpose of same method was applied for the measurement of the substrate simula

Substrate Particle Contact Model
At present, the traditional interparticle contact models are mainly JKR Cohesio ear Cohesion, Bradley, DMT, and Hertz-Mindlin with bonding, etc. These contact cannot describe the stress change process of bonded elastic-plastic particles effe The main components of the seedling substrate include peat, grass charcoal, wood etc. The substrate contains some moisture, and there is a complicated interaction b the particles, and they are adhesive. Elastic and plastic deformation will happen wet particles during compression. A certain amount of overlap will remain betw particles after the pressure is unloaded, and the contact area between the partic

Substrate Particle Contact Model
At present, the traditional interparticle contact models are mainly JKR Cohesion, Linear Cohesion, Bradley, DMT, and Hertz-Mindlin with bonding, etc. These contact models cannot describe the stress change process of bonded elastic-plastic particles effectively. The main components of the seedling substrate include peat, grass charcoal, wood chips, etc. The substrate contains some moisture, and there is a complicated interaction between the particles, and they are adhesive. Elastic and plastic deformation will happen to the wet particles during compression. A certain amount of overlap will remain between the particles after the pressure is unloaded, and the contact area between the particles will impact the adhesion intensity directly [18,19]. Hence, the Edinburgh Elasto-Plastic Adhesion (EEPA) model was selected as the contact model of substrate particles and the Hertz-Mindlin no slip model for the contact model of the substrate and the steel cylinder [20][21][22].

Substrate and Funnel Model
The substrate and the funnel shape model must be created before the simulation. The shape of the substrate particle is close to a sphere. The particle size distribution range is normally small, and the Rayleigh time step of the numerical simulation is limited by the smallest particle size in the model. Hence, according to the substrate particle diameter distribution, the particle diameter was set at 1.38 mm and the standard deviation at 0.04 mm. According to the test, the simulation time was reduced by over 90% if the positive distribution simulation was applied compared with multiple distributions of the particle diameter ( Figure 5). According to Han et al. [23] and COSMOS material base [24], the substrate particle shearing modulus was in the range of 2.9~10.0 MPa. The Poisson's ratio of the substrate particles was 0.464, referring to the method of Wang et al. [25].
A 3D model based on the same funnel and steel plate used in the test was created in Solidworks and saved as an STP file and imported into EDEM. The Poisson's ratio was 0.3; the density was 7.865 g/cm 3 , and shear modulus was 79,700 MPa after referring to the relevant literature [26].

Setting of Simulation Parameters
In the simulation test, 50 g particles were generated at the rate of 10 g/s. The Rayleigh time step of the simulation test was about 20%; the total simulation time was 5.5 s; the data-saving interval was 0.1 s, and the grid size was 2.5 times the minimum particle radius.

Plackett-Burman Design Test for Filtering Significant Parameters
Not all contact parameters and contact model parameters had a significant impact on the AOR. The insignificant parameters cannot be calibrated based on the AOR, or the calibration result would be inaccurate [27]. The parameters to be input in the EDEM of the EEPA model include constant pull-off force (f 0 ), surface energy (∆γ), contact plasticity ratio (λ p ), slope exp. (n), tensile exp. (X), and tangential stiff multiplier (k tm ). In order to reduce the factors that will influence the test, we defined f 0 = 0 and n = 1.5, according to the reference [22].
The Design-Expert software was employed for the PBD test design and analysis. The substrate particle shear modulus, contact parameter (substrate-substrate restitution coefficient, substrate-substrate static friction coefficient, substrate-substrate rolling friction coefficient, substrate-steel restitution coefficient, and substrate-steel static friction coefficient), and contact model parameters (surface energy, contact plasticity ratio, tensile exp., and tangential stiff multiplier) were filtered, and the parameters that are influential on the AOR were selected. According to the reference literature [28][29][30][31], the substrate-substrate restitution coefficient range was 0.10~0.70; the static friction coefficient range was 0.30~1.16, and the rolling friction coefficient range was 0.00~0.56. The restitution coefficient range between substrate and steel plate was 0.1~0.6, and the static friction coefficient range was 0.20~0.60. According to the reference of Xie et al. [22], Ma et al. [32], Wang et al. [33], and Thakur et al. [18], the surface energy range was 0.00~1.00 J·m −2 ; the contact plasticity ratio range was 0.30~0.70; the tensile exp. range was 1.50~5.00, and the tangential stiff multiplier range was 0.05~1.00.
The factors levels of the PBD test are shown in Table 1, and 13 test schemes and results of 10 factors are shown in Table 2. No. After the variation analysis on the test result using Design-Expert software, the influential effect and reliability analysis of 10 parameters on AOR are shown in Table 3. The results show that the static friction coefficient and the rolling friction coefficient between the substrate and the surface energy had a significant effect on the AOR (p < 0.05), while the rest of the parameters demonstrated no significant effect on it. The variation coefficient, CV = 3.03%, was less than 15%, indicating that the test results were all normal and worthy of analysis. The calibration determination coefficient R adj 2 = 0.9956 was greater than 0.8, indicating a good fitting of the regression model. Precision is the effective signal to noise ratio. This value is normally the bigger the better and should be greater than 4, generally speaking. The precision of this test was 46.4244, indicating good reliability of the test.

Path of Steepest Ascent Method for Determining the Optimal Value Range of the Significant Parameter
Before applying the RSM method to establish the regression model to solve the optimal value, it has to be verified that the optimal value of the factors is within the selected range of the highest and lowest values, and the path of steepest ascent method can be used to determine the optimal value range of the factors. According to the PBD test results, the middle level of the factors that had an insignificant impact on the AOR, namely, the substrate particle shear modulus of 6.45 MPa, substrate-substrate restitution coefficient 0.4, substrate-steel restitution coefficient of 0.35, substrate-steel static friction coefficient of 0.4, contact plasticity ratio of 0.5, tensile exp. of 3.25, and tangential stiff multiplier of 0.525. The path of steepest ascent method on 3 significant parameters (C-substrate-substrate static friction coefficient, D-substrate-substrate rolling friction coefficient, and G-surface energy) was conducted with reference to the PBD test factor level table. Meanwhile, the relevant error of simulative substrate AOR and the actual value were calculated. The test scheme and results are shown in Table 4.  The results showed that as C, D, and G increased; the substrate simulation AOR increased, and the relative error of the simulative and the actual AOR decreased first and then increased. Since the test result of the second group demonstrated the smallest relative error, the optimal range of the 3 significant parameters was selected in the proximity of that selected for the second group. Hence, the levels selected for the 1st, 2nd, and 3rd groups were selected for the RSM test, and a regression model was established to solve the optimal values of the parameters.

Establish the Regression Model Based on the RSM to Solve Significant Parameters
The optimal parameters were determined by adopting the RSM test during the CCD test referring to the path of steepest ascent method results. The middle level was selected for the nonsignificant parameters as in Section 3.4.2., and the simulation factor code for the significant parameters was shown in Table 5. A total of 23 testing spots were included in the CCD test scheme, which included 14 analysis factors and 9 zero-point estimation errors. The test design scheme and the response values are shown in Table 6.  The quadratic regression model was obtained by analyzing the CCD test results using Design-Expert software. Variance analysis results of the quadratic regression model are shown in Table 7. The p value of the regression model was <0.01, indicating a very significant relationship between the AOR and regression equation. The p value of lack of fit was 0.6548 (>0.05), indicating a small percentage of abnormal errors in the actual fitting and the regression equation and a good fitting. The variation coefficient CV = 4.57% indicated good test reliability. The determination coefficient R 2 = 0.9227 and the calibration determination coefficient R adj 2 = 0.9584 indicated a good regression equation reliability. The precision = 27.1262 indicated a good precision of the regression model. In addition, C, D, G, CD, and C 2 had an extremely significant effect on the AOR of the substrate (p < 0.01). The interaction term coefficients CG had a significant effect on the AOR (p < 0.05). DG, D 2 , and G 2 had no significant influence on the AOR (p > 0.05).
The quadratic regression model was optimized by eliminating the items with insignificant impact. The variance analysis results of the optimized regression model are shown in Table 8

Verification Test of Parameter Calibration
The Design-Expert software was used to solve the optimal value of the optimized regression equation with the angle of 39.03 • as the target. Several groups of AOR were simulated and verified, and a group of optimal solutions closest in the shape of the physical test was obtained, namely, the substrate-substrate static friction coefficient was 0.427; the substrate-substrate rolling friction coefficient was 0.039, and the surface energy was 0.228. EDEM software was used for the substrate AOR simulation test, and an intermediate level was selected for the nonsignificant parameters. The other settings were the same as the previous ones. The simulated substrate AOR was 38.51 • as tested. The relative error with the actual substrate AOR 39.03 • was 1.33%, indicating the reliability of the optimal values of 3 significant parameters. The comparison result between the simulation test and the physical test is shown in Figure 8, from which we can find that the respective substrate particle pile contours are similar.

Verification Test of Parameter Calibration
The Design-Expert software was used to solve the optimal value of the optimized regression equation with the angle of 39.03° as the target. Several groups of AOR were simulated and verified, and a group of optimal solutions closest in the shape of the physical test was obtained, namely, the substrate-substrate static friction coefficient was 0.427; the substrate-substrate rolling friction coefficient was 0.039, and the surface energy was 0.228. EDEM software was used for the substrate AOR simulation test, and an intermediate level was selected for the nonsignificant parameters. The other settings were the same as the previous ones. The simulated substrate AOR was 38.51° as tested. The relative error with the actual substrate AOR 39.03° was 1.33%, indicating the reliability of the optimal values of 3 significant parameters. The comparison result between the simulation test and the physical test is shown in Figure 8, from which we can find that the respective substrate particle pile contours are similar.

End-Effector Parameter Optimization
Due to the fact that the planting depth of pumpkin seeds is between 10~15 mm, the digging performance of the spiral digging end-effector means the cross-sectional width W was defined when the planting depth reached 13 mm (as shown in Figure 9) and by particle overflow percentage and dig time. The particle overflow percentage is the percentage of spilling substrate particles in the original particles after spiral digging. The dig time is

End-Effector Parameter Optimization
Due to the fact that the planting depth of pumpkin seeds is between 10~15 mm, the digging performance of the spiral digging end-effector means the cross-sectional width W was defined when the planting depth reached 13 mm (as shown in Figure 9) and by particle overflow percentage and dig time. The particle overflow percentage is the percentage of spilling substrate particles in the original particles after spiral digging. The dig time is the time from the spiral bit down to the stable shape of the substrate hole. The smaller the dig time, the higher the digging efficiency.
Due to the fact that the planting depth of pumpkin seeds is b digging performance of the spiral digging end-effector means the c was defined when the planting depth reached 13 mm (as shown in ticle overflow percentage and dig time. The particle overflow perce of spilling substrate particles in the original particles after spiral d the time from the spiral bit down to the stable shape of the substrat dig time, the higher the digging efficiency. According to the pre-experimental results, the digging performance of the spiral digging end-effector was mainly influenced by the spiral bit structure parameters and rotational speed of the end effector. The spiral bit structure parameters mainly included the cone angle and the spiral angle. Combined with the pre-experiment results, a single factor test was carried out. Since the maximum γ was 37.5 • , the cone angles were 25 • , 27.5 • , 30 • , 32.5 • , and 35 • , and the spiral angles were 65 • , 70 • , 75 • , 80 • , and 85 • , and the rotational speeds were 60, 120, 180, 240, 300, and 360 r/min. On the basis of the single factor test, the value range of each factor in the three-factor and three-level orthogonal test was obtained, as shown in Table 9. Each test was repeated three times, and the average value was taken as the test result.

Creating Simulation Model and the Determination of Parameters
According to the seedling substrate calibration results, the EEPA model was selected as the substrate particle contact model and the Hertz-Mindlin no slip model as the contact model between substrate and tray and the spiral bits. Please refer to the above-mentioned values for the particle size, Poisson's ratio of the model, shear modulus, contact parameters, and contact model parameters.
The spiral bits and the tray were modeled and assembled in SolidWorks, and the file was saved as an STP file and imported into EDEM, and the boundary model was obtained. In EDEM software, the tray was used as a particle generation vessel, and 25,000 particles were generated at a speed of 35,000 piece/s. The Rayleigh time step was roughly 20%, and the grid size was 2.5 times the minimal particle radius.

Method of Measuring the Bulk Density of the Substrate
The greater the soil porosity, the better its air permeability, which is beneficial to the growth of plant roots [34]. At the same time, soil bulk density is negatively correlated with soil porosity [35]. Therefore, the spiral digging end-effector and the traditional protrusions penetrating end-effector (as shown in Figure 10) were used for the digging test. The upper and lower surface area (S 1 and S 2 ) and the height (H s ) of the substrate were measured, and the substrate volume of the truncated pyramid (V 0 ) was calculated according to Equation (6). (6) omy 2022, 12, x FOR PEER REVIEW Figure 10. The traditional protrusions penetrating end-effector.
The volume of the hole (the shape of the hole of the spiral and the shape of the traditional protrusions penetrating end-effec amid) was recorded as V1, and the weight of the substrate (m) density (ρb) of the substrate was calculated according to Equation

Simulation Results and Analysis
(1) Single factor test Figure 11 shows the effects of cone angle, spiral angle, and ro section width, particle overflow percentage, and dig time under The volume of the hole (the shape of the hole of the spiral end-effector was a cone, and the shape of the traditional protrusions penetrating end-effector was a truncated pyramid) was recorded as V 1 , and the weight of the substrate (m) was weighed. The bulk density (ρ b ) of the substrate was calculated according to Equation (7).

Simulation Results and Analysis
(1) Single factor test Figure 11 shows the effects of cone angle, spiral angle, and rotational speed on crosssection width, particle overflow percentage, and dig time under the simulation test.

Simulation Results and Analysis
(1) Single factor test Figure 11 shows the effects of cone angle, spiral angle, and rotational speed on crosssection width, particle overflow percentage, and dig time under the simulation test.  Figure 11a shows that with the increase of the cone angle, both the cross-section width and the particle overflow percentage increased, while the dig time had no obvious change. When the cone angle was less than 30 • , the particle overflow percentage was less than 1.5% due to the small hole. At the same time, the smaller cone angle caused the substrate particles to slip and collapse, which cannot meet the requirement of 22 mm cross-section width. On the contrary, when the cone angle was 35 • , the particle overflow percentage was higher than 7.2%. Therefore, considering all the test indicators, the cone angle was 27.5~32.5 • . Figure 11b shows that as the spiral angle increased, the particle overflow percentage decreased. The results show that when the spiral angle was small, the pitch was large, and at the same time, more substrates were cut laterally, and the longitudinal lifting ability of the substrates was poor, resulting in the accumulation and overflow of substrate particles. As the spiral angle increased from 65 • to 80 • , the cross-section width increased from an average of 21.4 mm to 22.7 mm, and the dig time decreased from 3.00 s to 2.54 s. But when the spiral angle continued to increase to 85 • , the section width decreased slightly (21.9 mm), and the dig time increased slightly (2.68 s). Therefore, considering all the test indicators, the spiral angle was 75~85 • . Figure 11c shows that with the increase of the rotational speed of the spiral bit, the cross-section width and particle overflow percentage showed an obvious upward trend, while the dig time gradually decreased. The results show that the greater the rotational speed, the greater the longitudinal lifting force of the particles, the higher the hole digging efficiency, and the greater the kinetic energy of the substrate particles, which easily caused the particles to splash and increase the overflow. However, when the rotational speed was less than 180 r/min, although the particle overflow percentage was very low, the section width could not meet the seeding requirements, and the efficiency was low. Therefore, considering all the test indicators, the rotational speed was 180~300 r/min.
(2) Analysis of variance Table 10 shows the orthogonal test results of the spiral bits to dig holes. γ, α, and r represent the levels of the cone angle, the spiral angle, and the rotational speed. Table 11 shows the results of variance analysis. The primary and secondary orders of the experimental factors affecting the cross-section width, particle overflow percentage, and dig time were: γ > r > α, γ > α > r, and r > α > γ. When the cone angle was too large, the particle overflow percentage was large. However, if the cone angle was too small, the cross-section width could not meet the requirements of sowing. The spiral angle had no significant effect on the cross-section width. In addition, when the spiral angle was small, the particles overflowed more, and when the intermediate level was taken, the dig time was short. When the rotation speed was small, the dig time was longer, which affects the efficiency of digging. However, the larger the rotational speed, the more particles overflow. Therefore, the best combination of parameters for digging holes was γ 2 α 2 r 2 ; the cone angle was 30 • , the spiral angle was 80 • , and the rotational speed was 240 r/min. Since γ = 30 • , α = 80 • , and a = 22, the parameter equations of the helix line were Equation (8)    x = 11e 0.087t cos t y = 11e 0.087t sin t z = 19.05e 0.09t (8) where t ∈ [0, 450 • ].

Simulation Test Verification
First of all, the spiral bit model (Figure 12(1)) was designed and imported into EDEM. The substrate particles (Figure 12(2)) were put into the plug tray (Figure 12(3)) through the particles factory (Figure 12a), and the natural settlement of those particles was simulated (Figure 12b). Then the spiral bit rotated at a speed of 240 r/min, spun into the substrate particles (Figure 12c), and rotated for 2.09 s (Figure 12d). Finally, the spiral blade was lifted up (Figure 12e), and a hole was formed (Figure 12f). The width of cross-section of 23.1 mm at the depth of 13 mm was obtained. And then, an optimized spiral digging end-effector was used, as shown in Figure 13, to carry out the digging test. In order to obtain an accurate hole profile, when the digging was finished, the hole profile was measured by the self-made, small-sized soil profile gauge, as per test method [36,37] and as shown in Figure 14a. According to Table 12, the relative error of the cross-sectional width was less than 13.85%, and the profile of the substrate hole particles from the physical test was almost the same as that of the simulation test, as shown in Figure 14b. The width of cross-section of 23.1 mm at the depth of 13 mm was obtained. And then, an optimized spiral digging end-effector was used, as shown in Figure 13, to carry out the digging test. In order to obtain an accurate hole profile, when the digging was finished, the hole profile was measured by the self-made, small-sized soil profile gauge, as per test method [36,37] and as shown in Figure 14a. According to Table 12, the relative error of the cross-sectional width was less than 13.85%, and the profile of the substrate hole particles from the physical test was almost the same as that of the simulation test, as shown in Figure 14b. method [36,37] and as shown in Figure 14a. According to Table 12 cross-sectional width was less than 13.85%, and the profile of the from the physical test was almost the same as that of the simu Figure 14b.   Through experiments, the bulk density of the substrate excavated by the spiral digging end-effector was 0.622 ± 0.019 g/cm 3 , and the bulk density of the protrusions penetrating end-effector was 0.688 ± 0.005 g/cm 3 . The results show that the bulk density of the substrate excavated by the spiral type was reduced by about 9.6% compared with the penetrating type; soil porosity was increased, and soil aeration was improved. Therefore, the spiral digging end-effector helps to reduce the bulk density of the substrate, increase soil porosity, and promote root growth.

Conclusions
(1) A spiral digging end-effector for soil digging before planting pumpkin seeds was designed. It was mainly composed of five spiral bits, a synchronous belt, a dual-axial cylinder, a gear motor, and a connecting plate. (2) The EDEM simulation method was applied to calibrate the discrete element model parameters of the substrate, through the method of substrate falling from the funnel. The EEPA model was adopted as the particle contact model. The contact parameters and the model parameters (substrate-substrate static friction coefficient, substratesubstrate rolling friction coefficient, and surface energy) that had a significant impact

Comparison of the Substrate Bulk Density of Spiral Digging End-Effector with Traditional End-Effector
Through experiments, the bulk density of the substrate excavated by the spiral digging end-effector was 0.622 ± 0.019 g/cm 3 , and the bulk density of the protrusions penetrating end-effector was 0.688 ± 0.005 g/cm 3 . The results show that the bulk density of the substrate excavated by the spiral type was reduced by about 9.6% compared with the penetrating type; soil porosity was increased, and soil aeration was improved. Therefore, the spiral digging end-effector helps to reduce the bulk density of the substrate, increase soil porosity, and promote root growth.

Conclusions
(1) A spiral digging end-effector for soil digging before planting pumpkin seeds was designed. It was mainly composed of five spiral bits, a synchronous belt, a dual-axial cylinder, a gear motor, and a connecting plate. (2) The EDEM simulation method was applied to calibrate the discrete element model parameters of the substrate, through the method of substrate falling from the funnel. The EEPA model was adopted as the particle contact model. The contact parameters and the model parameters (substrate-substrate static friction coefficient, substratesubstrate rolling friction coefficient, and surface energy) that had a significant impact on the substrate AOR were sorted out by PBD and the path of steepest ascent method, and their optimal value ranges were determined. The optimal values for the parameters were obtained through the Central Composite Design response surface analysis test: 0.427 for the substrate-substrate static friction coefficient, 0.039 substratesubstrate rolling friction coefficient, and 0.228 for surface energy. According to the test verification results, there was no significant difference between the simulative AOR of the substrate and the actual AOR, and the calibrated substrate discrete element model parameters were reliable. (3) EDEM software was used to simulate and optimize the parameters of the spiral bit.
The analysis showed that when the cone angle was 30 • , the spiral angle was 80 • , and the rotational speed was 240 r/min; the section width at the hole depth of 13 mm was 23.1 mm, and the particle overflow percentage was about 2.05%. Compared with the traditional end-effector, soil porosity was increased, and soil aeration was improved. This result laid the academic foundation for the design and improvement of automatic pumpkin seed planting equipment.