Parameter Optimization and Testing of a Conveying and Soil-Removing Device for Tiger Nut ( Cyperus esculentus ) Mechanical Harvesting

: Aimed at solving the large power consumption and high operating cost problems associated with the process of harvesting tiger nuts, the use of a conveying and soil-removing device which removes sandy soil while conveying tiger nuts was proposed. The device was numerically simulated with and without vibration using EDEM software. The results showed that the vibrating force was more conducive to the complete removal of sandy soil and the effective conveyance of tiger nuts. The simulation testing was carried out using spiral speed, vibration amplitude, and vibration frequency as the independent variables and conveyance efﬁciency, sandy soil removal rate, etc., as the dependent variables. The test results showed that the optimal parameter combination was a spiral speed of 107 r/min, a vibration amplitude of 8.5 mm, and a vibration frequency of 10.7 Hz, under which the theoretical value of conveyance efﬁciency was 80.39%, the sandy soil removal rate was 84.61%, and the variation coefﬁcient of sandy soil removal velocity was 3.6%. The ﬁeld test determined that the relative errors of conveyance efﬁciency and sandy soil removal rate with theoretical values were 2.55% and 1.41%, respectively. The tiger nut damage rate was 1.16%, and the leakage rate was 0.52%. The results showed that the conveying and soil-removing device could meet tiger nut harvest performance requirements.


Introduction
The tiger nut, commonly called Cyperus esculentus, is a cyperaceous renascent herb featuring developed root systems, strong adaptability, high production, and ease of management [1,2].The tiger nut is a nutrient-rich raw food material rich in vegetable oil and fat, starch, dietary fiber, and microelements [3,4].The large-scale cultivation of tiger nuts can not only achieve sandy soil improvement, wind prevention, and sand fixation but can also increase people's economic income and improve their living standards.In the tiger nut harvesting process, conventional artificial harvesting has problems, such as high labor intensity, low operating efficiency, and high cost [5,6].Therefore, research on harvesting mechanization for the tiger nut production process is urgently needed.
In terms of research on the mechanized harvesting technology of tiger nuts, an earlier harvester in Spain was mainly comprised of a digging shovel, vibrating screen, walking wheel, and transmission device, and its overall design principle is still used today [7].In China, He et al. [8] optimized the parameters of a caterpillar self-propelled tiger nut harvester hoisting device, studying the influence of vibration frequency, amplitude, and other parameters on the soil-sieving rate.Harbin Dongyu Agricultural Engineering Machinery Co., Ltd. has designed a double-layer drum-type tiger nut harvester.Through the rotating double-layer drum screen and the grid bar screen body swinging in a reciprocating arc, the pellet-shaped straw sundries are released to separate the tiger nut [9].Xinxiang Dilong Pharmaceutical Machinery Co., Ltd.[10] and Beijing Changxin Agricultural Technology Co., Ltd.[11] adopted the structure of post-vibration two-stage separation to separate the mixture of tiger nut, root, and sandy soil.The Shaoyang Agricultural Machinery Research Institute have developed a tiger nut harvester that adopts a spiral screening machine, a reciprocating vibrating screen, and other schemes to complete soil separation [12].The main power consumption in the mechanized harvesting of tiger nuts comes from the process of digging, conveying, and separation.Currently, tiger nuts are widely planted in sandy soil in south Xinjiang, and the harvesting method is to dig a mixture of sand soil and tiger nuts together.There is no ideal technical solution to reduce power consumption.Therefore, the key issue is how to ensure the efficient separation of tiger nuts from sandy soil while reducing the power consumption of the machine as a whole, carrying out a combined conveyance and separation process.Vibration technology is widely applied in the agricultural machinery equipment design and manufacturing process, as it can not only destroy the soil structure to separate crops from soil efficiently, but also reduce the power consumption of the machine as a whole, thus reducing the cost of operation.Regarding research on crop separation techniques, Fu et al. [13], for example, analyzed the law of straw vibration for soil removal by removing the adhering soil and mixing it with yellow corn silage straw.Wei et al. [14] and Lü et al. [15] separated potatoes and soil with the cooperation of a vibration device and a screen, and the influence of vibration parameters on the separation performance was analyzed.Hu et al. [16] and Wang et al. [17] studied the transportation and separation process in peanut harvesting.Peanut fruit, seedlings, and soil are thrown up and the soil removed under the vibration of the screen.Zhang et al. [18] researched the use of vibration for soil removal in the Panax notoginseng harvesting process and analyzed the influence of elevator speed, vibration frequency of the vibrating screen sieve, speed of rototiller parameter on root damage rate, and harvesting net rate.Liu et al. [19] conducted kinematics analysis of a garlic harvester's vibration separation device, achieving the separation of the bulb in the digging mixture from the soil and gravel.With regard to vibration resistance and consumption reduction, Shahgoli et al. [20] discovered that the vibration effect of the soil subsoiling process can destroy the soil structure for chessom, loosening the soil and reducing operation resistance.Rao et al. [21] researched the role of four-bar mechanism vibration in soil cultivation through testing.Niyamapa et al. [22] analyzed soil movement law under conditions with vibration and without vibration.In addition, others researched the effects of eccentric rotation speed and eccentricity on soil loosening in the yam, rice, and radish harvesting processes, and the test results showed that vibratory action could loosen the soil and facilitate artificial harvest [23][24][25].It is concluded that the vibration principle is conducive to the effective separation of rhizome crops, such as potato and peanut, from the soil and the fragmentation of the soil.At the same time, the simulation method can be used to analyze the material transport process and to optimize the mechanism parameters.
To solve the problems existing in the harvesting of tiger nuts planted in sandy soil in south Xinjiang, a tiger-nut-conveying and soil-removing device was proposed by the project team.EDEM discrete element software was used to carry out the numerical simulation of the conveying and soil-removing process of tiger nut harvesting.The aims of the study were to: analyze the influence of the different spiral speeds, vibration amplitudes, and vibration frequencies of the device on different evaluation indexes of tiger nut conveyance efficiency, sandy soil removal rate, and variation coefficient of sandy soil removal velocity; determine the primary and secondary relationship of such influence factors; obtain the optimal parameter combination; conduct experimental verification; and provide a technical and theoretical basis for promoting the operating quality of the tiger nut harvester.

Structure and Working Principle of Tiger Nut Harvester
The project team proposed that the tiger nut harvester be mainly composed of rotary tillage digging device, conveying and soil-removing device, scraper lifting device, threshing device, cleaning device, stock bin, crawler walking device, etc., as shown in Figure 1.In the harvester's operation process, the crawler traveling device effectively ensures traveling performance in the sand.The digging device digs a mixture of tiger nuts, roots, and a large amount of sandy soil and throws it to the conveying and soil-removing device, removing most of the sandy soil in the mixture.The scraper lifting device then transports the mixture of tiger nuts and roots containing a small amount of sandy soil to the threshing device.After threshing, whole tiger nuts and a small part of the roots fall into the cleaning device, with most of the roots being discharged from the straw outlet.The mixture of tiger nuts and fibrous roots is discharged from the rear end of the vibrating screen through the combined action of the fan and vibrating screen, and the tiger nuts are transported by the collection and transportation device to the stock bin.

Structure and Working Principle of Conveying and Soil-Removing Device
The conveying and soil-removing device is mainly composed of a feed scraper, feed scraper shaft, spiral blade, spiral shaft, bar screen, vibration exciter, vibratory element, vibration plate, regulating port, return spring, discharge scraper, frame, etc.Four vibration exciters are averagely arranged in the radial direction of the feed scraper shaft.The vibratory element is installed on the bar screen, which is connected with the frame through the return spring.The structure is shown in Figure 2.During the operation of the tiger nut harvester, the combination of the rotary tillage digging device and the fixed digging shovel will throw the excavated tiger nuts, roots, and a large amount of sandy soil mixture to the feed scraper, which will then throw the mixture to the screw conveyor.In this process, the feed scraper shaft drives the vibration exciter to rotate, and the vibration exciter contacts with the vibratory element, forcing the bar screen to move away from the vibration exciter.The return spring then reacts to the movement of the bar screen to complete the reciprocating movement.The mixture of tiger nuts, roots, and sandy soil is in a loose state under the vibration of the bar screen, and the sandy soil falls to the ground in the gaps in the grid sieve to complete the removal operation.The mixture moves along the circumferential direction of the spiral blade and in the axial direction to the end under the thrust of the spiral blade.A discharge scraper is arranged at the end of the spiral shaft to throw the tiger nuts, the roots, and a small amount of soil to the scraper lifting device, thereby completing the tiger nut conveyance operation.

Analysis of Conveyance Velocity of Tiger Nut
During the operation process of the screw conveyor, the tiger nuts, roots, and sandy soil mixture move upwards along the circumferential direction of the spiral blade under a pushing action.When the mixture reaches a certain position, it collapses downward and flows in the axial direction at the same time.The movement process is required to improve the tiger nut conveyance efficiency; therefore, the tiger nut conveyance conditions are analyzed.The changes in radial velocity and axial velocity of the tiger nuts in the spiral conveyance process are examined.During the operation, the spiral conveying shaft rotates around the x-axis at a constant speed.Taking any tiger nut particle, the movement state of the particle at a distance from the center of the spiral shaft and on the surface of the spiral blade is analyzed and a variable velocity movement is performed under the push of the spiral blade.The velocity analysis is shown in Figure 3.According to the velocity analysis in the figure, the component v x of the tiger nut particle velocity in the axial direction is: where α is the helix angle ( • ), β denotes the friction angle between the particles and the helicoid ( • ), v x represents the axial velocity of particle (m/s), and v stands for the resultant velocity of the particle (m/s).
where v 0 is the implicated velocity of the particle (m/s) and v n represents the vertical velocity of the particle and the spiral blade (m/s).
From the above analysis, when there is no friction, the particle velocity is v n .When there is friction, the particle velocity will become v after the deflection of β in the direction: where r is the distance from the particle to the center of the spiral axis (mm), w refers to the spiral speed (r/min), and t denotes the conveyance time (s).At the same time, because f = tanβ, and considering the maximum linear velocity at the edge of the spiral blade, the critical principle is applied ∆r = r.Substitute Equations ( 2) and (3) into Equation ( 1) to obtain: During the calculation, it can be seen that: According to Equations ( 4) and ( 5), it can be found that: According to the same calculation method, the radial velocity of the particle can be obtained: It can be seen from Equations ( 6) and ( 7) that the axial velocity and radial velocity of tiger nut particle are both related to the rotating speed of the spiral shaft.Therefore, based on the determination of the structural parameters of the spiral auger, further research on the rotating speed is required, with the critical condition being that the axial velocity is greater than the radial velocity to ensure the tiger nut conveyance efficiency.

Analysis of Bar Screen Motion Process
The machine travels at speed v, the vibration exciter interacts with the vibratory element when rotating, and the vibratory element and the bar screen are welded and fixed together, forcing the bar screen to move, as shown in the Figure 4.The adjusting hole can adjust the position of the vibration exciter between point O 1 and point O 2 , and the adjustment range is h, which can adjust the values of different amplitudes.During the contact between the vibration exciter and the vibratory element, the vibration exciter makes a circular motion with point O as the center, and the contact process with the vibratory element is arc length ab.The contact starts at point a, the bar screen starts to move towards position A 1 B 1 , and the amplitude begins to increase gradually.When reaching the highest point b, the bar screen reaches the maximum limit position A 2 B 2 and the amplitude reaches the maximum value.Therefore, the distance between AB and A 2 B 2 is the amplitude variation range of the bar screen.When the vibration exciter leaves point b, the vibratory element loses force, the bar screen returns from position A 2 B 2 to AB under the tension of the return spring, and the next vibration exciter repeats the previous process to realize the reciprocating motion of the bar screen.The sandy soil, tiger nuts, and roots used in the experiment were sampled from the 54th Regiment of the 3rd Division, Xinjiang.The planting pattern is shown in Figure 5.The planting method was: flat ground, planting spacing M of 35-40 cm, row spacing L of 35-40 cm, stubble height H of 8-10 cm, and growth depth h of 10-12 cm.The sandy soil particle with a radius of 1.5 mm was modeled as a single ball.According to a previous determination, the average size of three axes of tiger nut particle is 10.89 mm; the average length and diameter of the main root of tiger nuts are 53.75 mm and 6.38 mm, respectively; and the average diameter and length of fibrous roots are 0.36 mm and 67.43 mm, respectively.SolidWorks software was used for the three-dimensional modeling of a tiger nut particle and root saved in IGS format.The model was imported into EDEM software, and a multi-ball model was used for filling.A Hertz-Mindlin non-sliding contact model was used between balls [26][27][28], and the simulation models are shown in Figure 6.In addition, to facilitate the simulation calculation, the parts that were irrelevant to the conveying and soil-removing device were removed [29].

Simulation Parameter Setting
The intrinsic factors (Poisson's ratio, shear modulus, etc.) and basic contact parameters (collision coefficient of restitution, static friction coefficient, and rolling friction coefficient) of sandy soil particles, tiger nuts, and the conveying and soil-removing device geometry were set.In the process of simulation parameter setting, the intrinsic factors of tiger nuts and basic contact parameters between tiger nuts and between tiger nuts and geometry were determined via calibration [30].The intrinsic parameters of sandy soil and the contact parameters between sandy soil and tiger nuts and between sandy soil and geometry were obtained through reference [31][32][33].The geometry material was steel, and the intrinsic parameters of geometry were obtained through reference [34,35].See Tables 1 and 2 for the simulation parameter setting.A Hertz-Mindlin non-sliding contact model was used as the contact model between tiger nut, root, and sandy soil particle and between them and geometry [36,37].The generation planes of tiger nut particle, tiger nut root, and sandy soil particle were set at the feeding inlet of the conveying and soil-removing device, and the property was set as virtual.At the outlet of the device, a hexahedral box with the top surface removed was set for collecting sandy soil, tiger nuts, and root stock from the discharge scraper.According to the planting mode of tiger nuts and the working width and depth, the mass ratio of the tiger nuts, root, and sandy soil fed into the conveying device was measured as 1:1.66:156.25.Considering the simulation speed, the mass of sandy soil particles was set as 25 kg, the quantity of tiger nut particles was set as 20 (160 g), the quantity of tiger nut roots was set as 5 (265.6 g) with random generation location, the Rayleigh time step was set as 29%, and a simulation process duration of 12 s was determined.
After simulation was completed, the vibratory sandy soil removal process was divided according to net region into six regions: Region 1 was the sandy soil discharge volume from the discharge scraper transportation to the hexahedral box; Region 2 was the sandy soil discharge volume of the discharge scraper and grid screen; Region 3 was the radial sandy soil discharge volume of the spiral blade; Region 4 was the sandy soil discharge volume below in the spiral conveyer conveyance process; Region 5 was the sandy soil discharge volume of the feed scraper; and Region 6 was the sandy soil discharge volume at the feeding inlet.See Figure 7 for the schematic diagram.

Contrastive Analysis of Sandy Soil Removal Process Simulation with and without Vibratory Action
To verify the reliability of the simulation analysis using EDEM software, a contrastive analysis of the conditions with and without vibration forces during the operation of the conveying and soil-removing device was carried out.The sandy soil feed quantity, spiral speed, and other factors had a fixed value; the vibration amplitude, vibration frequency, and other factors were variables in the simulation.To make clear the motion state of the sandy soil with and without vibratory action, velocity simulation pictures at 1.4 s and 2.8 s in the operating process under the two conditions were extracted, among which Figure 8a,b show the sandy soil removal process without vibratory action and Figure 8c,d display the sandy soil removal process with vibratory action.Sandy soil velocity variation is represented by color bars, among which the color changing from blue to green and then to red means that the sandy soil velocity increases.As shown in the figures, sandy soil is principally blue at 1.4 s without vibratory action, indicating a low velocity, with the discharged sandy soil concentrated in Regions 4 and 5. Sandy soil is principally green with vibratory action, indicating a low velocity, with the discharged sandy soil concentrated in Regions 3 and 4. At 2.8 s, the discharge distance of the sandy soil in Region 3 with vibratory action is longer than that without vibratory action, and the discharged sandy soil is loose in the former conditions and in a stacked state in the latter conditions.Under such conditions, the subsequent feeding of the tiger nut, root, and sandy soil mixture will result in a late reset of the reset spring and an easy blockage in the geometry, which are both unfavorable for stable and continual operation.Meanwhile, the sandy soil is in a stacked state without vibratory action when delivered in geometry, which will not effectively destroy the mixture structure during delivery and is detrimental to the separation of the tiger nut root from the sandy soil, subsequently causing an increase in sandy soil quantity after the tiger nut enters the threshing device and severely influencing threshing quality.In summary, the device can be used to finish the sandy soil removal operation under two such conditions, but for constant feed quantity, the device with vibratory action can accomplish a greater sandy soil removal effect, achieve superior sandy soil removal distribution quality, and ensure better machine stability and continuity than that without vibratory action.This conclusion is consistent with the results of Xie's study [38].Combined with the sandy soil particle velocity analysis results, the vibratory action ensures that the sandy soil particles maintain an active flow process in the geometry, extract kinetic energy from the sandy soil particle group, and further analyze the change in sand particle flow velocity.Figure 9 shows that the kinetic energy change trend is generally the same-increasing first, then maintaining in a certain numerical range, and finally decreasing.In the case of vibration, the kinetic energy value remains stable in the range of 1 s to 4 s, gradually decreasing after 4 s.In the case of no vibration, the kinetic energy value remains stable in the range of 1.5 s to 3.5 s, gradually decreasing and approaching 0 after 3.5 s.The kinetic energy of the sandy soil particles with vibration action is higher than that without vibration action.To observe the kinetic energy change more clearly, the total kinetic energy of the sandy soil particles at 3 s under the two conditions was extracted, as shown in Figure 10.Setting the kinetic energy range to be the same, the flow energy variation is represented by color bars.Figure 10a shows the sandy soil removal process without vibratory action.Most of the sandy soil particles in the geometry are blue in color, and a few are green.Figure 10b displays the sandy soil removal process with vibratory action.Most of the sandy soil particles in the geometry are green, and a few are red.The red particles are distributed between the bar screen and the spiral conveying blade.This may be due to the reciprocating vibration force exerted by the bar screen on the particles, with the particles being pushed and centrifuged by the spiral blade at the moment of being bounced up, resulting in an increase in particle flow kinetic energy.At the same time, the color of the sandy soil particles at the discharge outlet is mostly green under the two conditions due to the same rotation speed of the discharge scraper and the same conveyance speed of the sandy soil particles, so the particle flow energy tends to be consistent.The results show that the total kinetic energy of the particles with vibration is significantly higher than that without vibration, which is consistent with the velocity analysis results.To understand the movement of the sandy soil particles in the geometric body more intuitively, a position change diagram of the sandy soil particles under the double action of grid screen and spiral conveying blade at 3 s was extracted, as shown in Figure 11, and in which Figure 11a displays the sandy soil removal process without vibratory action.It can be seen from Figure 11a that the sandy particle flow moves downward along the screen surface.This may be because the sandy soil particle flow first moves upward along the screen under the action of only the thrust and centrifugal force of the spiral blade.At the moment the force disappears, the sandy soil particles fall to the ground through the bar screen under their own gravity, and most of the sand is discharged below the bar screen.Figure 11b shows the sandy soil removal process with vibratory action.The sandy soil particle flow moves upward along the screen surface, which is probably because the sandy soil particle flow is forced by the reciprocating vibration of the screen.The particles are pushed and centrifuged by the spiral blade at the moment of being bounced up, making the particles pass through the screen quickly with an initial velocity.Due to the initial velocity of the particles, most of the sandy soil is discharged behind the screen.It can be seen from the comparison chart that the active state of the particles with vibration is significantly higher than that without vibration, and the sand-blocking phenomenon can be avoided under the bar screen, which is conducive to the sand removal operation.This conclusion is consistent with the above velocity and kinetic energy analysis results.Through the above analysis, it can be seen that vibratory sandy soil removal is the better operation mode in the tiger nut harvesting process.Therefore, the simulation of vibratory sand removal was analyzed.The sandy soil discharge volume from large to small is as follows: Region 3, which is the radial sandy soil discharge volume of the spiral blade; Region 4, which is the sandy soil discharge volume below in the spiral conveyer conveyance process; Region 2, which is the sandy soil discharge volume of the discharge scraper and grid screen; Region 5, which is the sandy soil discharge volume of the feed scraper; and Region 6, which is the sandy soil discharge volume at the feeding inlet.The sandy soil discharge volume of the former three regions is approximately 71% of the total discharge volume, indicating that sandy soil can be separated and discharged when the spiral conveyer is used for conveying tiger nuts, roots, and sandy soil.The weight of the sandy soil discharged is greatest in Region 3, which is good for the device discharge volume and prevents contact of the bar screen with the ground and bar screen blockage caused by too much sandy soil being discharged.
We extracted any movement locus of tiger nuts and roots with vibratory action at the end of the simulation.Figure 12 shows that tiger nuts and roots will appear at the feeding inlet and go through the feed scraper, spiral blade, and discharge scraper over time, ultimately ending in the discharge port, indicating that vibratory action can achieve the separation of the sandy soil and effectively convey the tiger nuts and roots.

Simulation Test Design
The test was designed in accordance with Box-Behnken central composite design theory to analyze the influence of different spiral speeds, vibration amplitudes, and vibration frequencies on the sandy soil removal rate and the variation coefficient of the sandy soil removal velocity of the vibratory sandy-soil-removing device.The level range of factors was determined according to the preliminary experiment.Sandy soil cannot be separated adequately if the spiral speed is too high, as this results in an increased sandy soil volume in the downstream lifting and conveying device and is detrimental to the threshing and separation operation.If the spiral speed is too low, the conveyance efficiency is lowered, although the sandy soil separation effect is achieved.According to the characteristics of low water content and good fluidity of sandy soil planted with tiger nuts, the high-frequency, low-amplitude vibration scheme is selected.When the vibration frequency and amplitude are too high, it can provide a large vibration force, destroying the structural form of the tiger nut root and sandy soil mixture, and maximizing the quality of sand removal.However, it is easy to cause the phenomenon whereby the vibrating screen is not reset in time, which leads to the decline in the stability of sand removal.Based on the analysis of the single factor pre-test results, the spiral speed range is 100-120 r/min, the amplitude range is 5-15 mm, and the vibration frequency is 10-12 Hz.The test factors and levels are shown in Table 3.

Evaluation Index
The conveyance efficiency is taken as the evaluation index describing the tiger nut conveyance performance.The ratio of the average velocity of tiger nuts in the spiral conveyance to the axial velocity of the spiral blade is taken as the tiger nut conveyance efficiency.It is calculated according to Equation (8): where Y 1 is the tiger nut conveyance efficiency (%), v 0 represents the average velocity of the tiger nuts (m/s), and v x is the axial velocity of the spiral blade (m/s).
The performance of the vibratory sandy-soil-removing device is described using the sandy soil removal rate as the evaluation index.It is calculated according to Equation ( 9) when the target and actual sandy soil discharge volume are determined: where Y 2 is the sandy soil removal rate (%); M 0 stands for the actual sandy soil discharge volume, which is the sum of the sandy soil discharge volumes of Regions 2, 3, 4, 5, and 6 (kg); and M represents the target sandy soil discharge volume (25 kg).The stability of the sandy soil removal velocity is described by the coefficient of variation.The smaller the coefficient of variation, the better the stability of the sandy soil removal velocity.This is calculated according to Equation (10): where Y 3 is the variable coefficient, S represents the standard deviation of samples, and V denotes the mean value of samples.

Field Test Method
The tiger nut, root, and sandy soil mixtures were fed into the screw conveyor, and the continuous passing time of the tiger nut through the screw conveyor was recorded.The actual axial velocity of the tiger nut was calculated by the length of the screw auger, and the conveyance efficiency of the tiger nut was calculated according to Equation (8).
In addition, the total mass of sandy soil was calculated within a length of 5 m of the test area as the target sandy soil discharge volume, and the total mass of sandy soil was collected and weighed in the conveying device and the threshing device.The difference between the target sandy soil discharge volume and the total mass weighed was used to obtain the actual sandy soil discharge volume, with the sandy soil removal rate calculated according to Equation (9).

Analysis of the Simulation Test Results
The three-factor and three-level test was designed using spiral speed, vibration amplitude, and vibration frequency as the factors and the conveyance efficiency, sandy soil removal rate, and variation coefficient of sandy soil removal velocity as the evaluation index.The test programs included 17 test points, among which 12 test points were analysis factors, with five zero-estimation errors.See Table 4 for the design scheme and results.The variance analysis of the test results in Table 4 was conducted using Design-Expert 8.0.6 software.See Table 5 for the analysis results.The regression model equation was calculated between the evaluation indexes of conveyance efficiency, sandy soil removal rate, and variation coefficient of sandy soil removal velocity and the independent variables of spiral speed, vibration amplitude, and vibration frequency.The significance of the index influence of various variables in the regression equation was judged by the F value, the lower probability p, and the higher significance of the corresponding variables [39].The regression model p value is less than 0.0001, indicating that the model is extremely significant.The lack-of-fit p value is 0.0525 and more than 0.05; that is, the lack of fit is not significant, indicating that it is within the experiment range and that the prediction value of the regression model has a good degree of fit with the actual value.The determination coefficient of conveyance efficiency, sandy soil removal rate, and variation coefficient of sandy soil removal velocity are 0.9810, 0.9843, and 0.8893, respectively, indicating that the prediction value is highly correlated to the actual value and that the model can be used for the analysis and prediction of the experiment index.

Influence of Various Factors on Conveyance Efficiency
In the regression equation, the absolute value of the factor coefficient determined the influence degree of the test index.According to regression Equation (11) and Table 5, factors X 1 , X 3 , X 1 X 2 , X 2 2 , and X 3 2 are extremely significant, factor X 2 is significant, and the other factors are not significant in the conveyance efficiency model.The influence and importance of various factors on conveyance efficiency from high to low are: spiral speed X 1 , vibration frequency X 3 , and vibration amplitude X 2 .To observe the influence law more clearly, the three-dimensional response surface was extracted, as shown in Figure 13.It shows that when the amplitude is constant, with the increase in spiral speed, the conveyance efficiency will increase slowly and then decrease slowly.This is probably due to the fact that as the spiral speed increases, the tiger nut axial velocity also increases.When the spiral speed is too great, the tiger nut will move in a circle with the spiral blade, leading to a decrease in the conveyance efficiency.When the spiral speed is constant, with the increase in amplitude, the conveyance efficiency will increase and then decrease.This is possibly due to the fact that when the amplitude is too large, the tiger nut will move back and forth with the bar screen, reducing contact with the spiral auger.

Influence of Various Factors on Sandy Soil Removal Rate
According to regression Equation ( 12) and Table 5, factors X 1 , X 2 , X 3 , X 1 2 , X 2 2 , and X 3 2 are extremely significant; factor X 1 X 2 is significant; and the other factors are not significant in the sandy soil removal rate model.The influence and importance of various factors on the sandy soil removal rate from high to low are: vibration amplitude X 2 , vibration frequency X 3 , and spiral speed X 1 .To observe the influence law more clearly, a threedimensional response surface was extracted, as shown in Figure 14.It shows that when the spiral speed is at the 0 level, the sandy soil removal rate will increase slowly and then decrease slowly as the vibration amplitude increases.This is probably due to the facts that the vibration amplitude determines the vibratory action and that the sandy soil removal rate will increase with increasing vibration amplitude.However, if the vibration amplitude is too high, the reset spring will not reset the bar screen in a timely manner and its vibratory action will decrease for the sandy soil, resulting in stacked sandy soil on the bar screen and a decrease in the sandy soil removal rate.When the vibration amplitude is at the 0 level, the sandy soil removal rate will increase slowly and then decrease rapidly as the vibration amplitude increases.This is possibly because the spiral speed correlates to the axial and tangential speed in the sandy soil conveyance process.If the tangential speed increases, sandy soil will be thrown out along the spiral blade and be discharged out of the machine via the bar screen, so the sandy soil removal rate will rise.However, if the tangential speed is too high, sandy soil cannot be thrown out in a timely manner and can only move in the axial direction, which will result in low conveyance efficiency, stacked sandy soil in the machine, and a decrease in the sandy soil removal rate.

Influence of Various Factors on Variation Coefficient of Sandy Soil Removal Velocity
According to regression Equation ( 13) and Table 5, factor X 3 is extremely significant; factors X 1 , X 1 X 2 , and X 1 2 are significant; and the other factors are not significant in the variation coefficient of the sandy-soil-removing velocity model.The influence and importance of various factors on the sandy soil removal rate from high to low are: vibration frequency X 3 , spiral speed X 1 , and vibration amplitude X 2 .To observe the influence law more clearly, a four-dimensional slice and three-dimensional response surface were extracted, as shown in Figure 15.It can be seen that with the increase in spiral speed, the variation coefficient first decreases and then increases and that the surface is steep.When the spiral speed is approximately 110 r/min, the variation coefficient reaches the minimum.When the spiral speed takes a certain value, with the increase in amplitude, the coefficient of variation first decreases and then increases, and the surface changes gently, indicating that the spiral speed has a significant effect on the coefficient of variation compared with the amplitude.

Operating Parameter Optimization
To find the optimal combination of various factors, the maximum values of conveyance efficiency and sandy soil removal rate, and the minimum value of the variation coefficient of sandy soil removal velocity were calculated using the optimization function provided in the Design-Expert software: The optimal influence factor parameter combination is a spiral speed of 107 r/min, a vibration amplitude of 8.5 mm, and a vibration frequency of 10.7 Hz.In this combination, the conveyance efficiency is 80.39%, the sandy soil removal rate is 84.61%, and the variation coefficient of sandy soil removal velocity is 3.6%.

Field Test Result
To verify the reliability of the optimization results, experimental verification was conducted in October 2020 in the 54th Regiment of the 3rd Division, Xinjiang, as shown in Figure 16.The parameters of tiger nut harvesting are shown in Table 6.During the tiger nut harvesting period, the average moisture content of the sandy soil 15 cm below the ground surface was 4.91%, as measured using a TDR300 soil moisture meter.The working width of the machines was 1.3 m, and the working depth was 15 cm.It was found that the optimal combination of operating parameters for the conveying and soil-removing device is a spiral speed of 107 r/min, a vibration amplitude of 8.5 mm, and a vibration frequency of 10.7 Hz.The experiment was carried out three times, and the average value of the test results was taken.The mean value of the conveyance efficiency was 78.34%, the sandy soil removal rate was 83.42%, and the relative errors with theoretical values were 2.55% and 1.41%, respectively.Under this parameter combination, the tiger nut damage rate was 1.16% and the tiger nut leakage rate was 0.52%.The results showed that the conveying and soil-removing device could meet the tiger nut harvest performance requirements.In addition, in the field test process, the water moisture content of the sandy soil directly affected the sandy soil removal rate index.If the water moisture content of sandy soil is too high, the bar screen gap will be blocked.In future, it will be necessary to further study the water moisture content of sandy soil particle simulation modeling to ensure the accuracy of the simulation test and the field test to determine the best harvest date for tiger nuts.

Conclusions
A conveying and soil-removing device that removes sandy soil while conveying tiger nuts was proposed.The factors affecting conveying and soil-removing performance were analyzed, and the two operational processes of the device with and without vibration were numerically simulated using EDEM software.The sandy soil movement velocity changes under the two conditions were compared and analyzed, and the movement tracks of tiger nut particles and tiger nut roots under the conditions of vibration were extracted.The showed that the vibrating force was more conducive to the completion of the sandy soil removal operation and the effective conveyance of tiger nuts and roots.
The simulation test was carried out using spiral speed, vibration amplitude, and vibration frequency as the independent variables and with tiger nut conveyance efficiency, sandy soil removal rate, and the variation coefficient of sandy soil removal velocity as the dependent variables.A Box-Behnken central composite design was applied for the parameter optimization of the three factors affecting the dependent variables.The test results showed that optimal parameter combination was a spiral speed of 107 r/min, a vibration amplitude of 8.5 mm, and a vibration frequency of 10.7 Hz, under which the theoretical value of the conveyance efficiency was 80.39%, the sandy soil removal rate was 84.61%, and the variation coefficient of the sandy soil removal velocity was 3.6%.
The field experiment conducted to determine the optimal parameter combination showed the conveyance efficiency was 78.34%, the sandy soil removal rate was 83.42%, and the relative errors with theoretical values were 2.55% and 1.41%, respectively.Under this parameter combination, the tiger nut damage rate was 1.16% and the tiger nut leakage rate was 0.52%.The test results showed that the conveying and soil-removing device could meet the tiger nut harvest performance requirements.This study could provide the basis for the design and optimization of a tiger nut harvester.

Figure 8 .
Figure 8. Velocity simulation comparison diagram: (a) without vibratory action at 1.4 s; (b) without vibratory action at 2.8 s; (c) with vibratory action at 1.4 s; (d) with vibratory action at 2.8 s.

Figure 10 .
Figure 10.Kinetic energy simulation comparison diagram: (a) sandy soil removal process without vibratory action; (b) sandy soil removal process with vibratory action.

Figure 11 .
Figure 11.Position simulation comparison diagram: (a) sandy soil removal process without vibratory action; (b) sandy soil removal process with vibratory action.

Figure 12 .
Figure 12.Movement locus: (a) movement locus of tiger nut; (b) movement locus of tiger nut root.

Figure 13 .
Figure 13.Influence of factor interaction on conveyance efficiency.

Figure 14 .
Figure 14.Influence of factor interaction on sandy soil removal rate.

Figure 15 .
Figure 15.Influence of factor interaction on the variation coefficient of sandy soil removal velocity.
Note: a indicates coefficient of restitution; b indicates coefficient of static friction; and c indicates coefficient of rolling friction.

Table 3 .
Test factor and level.

Table 4 .
Experimental design scheme and results.

Table 5 .
Variance analysis results of response surface model.