Parameter Optimization of a Conveying and Separating Device Based on a Five-Stage Screw and Vibrating Screen for Tiger Nut Harvesters

: To tackle problems such as the difficult separation from sand and the high power consumption of tiger nut harvesting in the sandy areas of Xinjiang, a conveying and separating device for tiger nut harvesters was designed. The axial and radial migrations of materials under screw action and the separation process of materials under vibratory action were analyzed dynamically. A simulation analysis was carried out on the conveying and separating process based on EDEM software. The migration trajectories of tiger nuts and sand particles were extracted, the displacement variations of sand particles on the X -axis, Y -axis, and Z -axis were analyzed in the action area of the screen-cleaning spike teeth and the screw action area, respectively, and the conveying and separation law of the tiger nut harvest mixture was clarified. With key parameters such as the screw velocity ratio, amplitude, vibration frequency, and machine operation velocity as test factors, and with the sand removal rate, crushing rate, and power consumption as test evaluation indicators, a four-factor, five-level orthogonal central composite test design was implemented. The test results were analyzed via the regression variance analysis method, and relation models between variable factors and evaluation indicators were constructed. The test results show that under the combined conditions of a screw velocity ratio of 0.88, an amplitude of 4.7 mm, a vibration frequency of 7.5 Hz, and a machine operation velocity of 0.92 km/h, the sand removal rate is 90.40%, the crushing rate is 1.66%, and the power consumption is 2.24 kW in theory. The optimized results were verified by tests. The sand removal rate was 88.92%, the crushing rate was 1.71%, the total power consumption was 2.29 kW, and the errors from the predicted values were 1.6%, 3.0%, and 2.2%, respectively, meeting the requirements for tiger nut harvesting conveyance and separation. This research can provide support for the development of technology and equipment for mechanized harvesting of tiger nuts in the sandy areas of Xinjiang.


Introduction
The tiger nut is a perennial herb of the Cyperaceae family.It is characterized by a developed root system, high adaptability, and strong stress resistance [1][2][3].Large-scale planting of tiger nuts can not only modify sandy land, prevent wind, and fix sand, but also increase farmers' economic income and improve people's living standards [4].In China, tiger nuts are mainly planted in Xinjiang, Inner Mongolia, Jilin, and Henan.Harvesting is an important aspect of the production of tiger nuts [5].At present, the plantation of tiger nuts in Xinjiang mainly faces the problem of high power consumption during harvesting.Therefore, it is particularly important to study how to ensure efficient separation of tiger nuts from sand while reducing the power consumption of the whole machine during the conveying and separation process of the tiger nut excavation mixtures.
In recent years, domestic and foreign scholars have studied the related technologies and devices for the conveyance and soil separation of tiger nut harvest and excavation mixtures, achieving some results.The tiger nut harvester, which was used previously in Spain, mainly uses a vibrating screen structure to separate nuts from soil.Its general design principle is still being used today [6].Zhao Zhuang et al. [7] designed a scraper-type tiger nut conveying device and studied the effects of parameters such as the scraper height, scraper width, and scraper spacing on the conveying performance.He et al. [8] designed a vibrating chain-plate conveying structure and studied the effects of parameters such as the vibration frequency and amplitude on the soil screening efficiency.Liu Xuemin et al. [9] used a vibrating screen structure to separate tiger nuts from clay.Harbin Dongyu Agricultural Engineering Machinery Co., Ltd.[10] separated tiger nuts, straw debris, and soil using a rotating doublelayer drum screen and a grid-bar screen body that oscillates reciprocally along an arc.Xinxiang Dilong Pharmaceutical Machinery Co., Ltd.[11] used a post-vibration two-stage separation structure to separate tiger nuts, roots, and sand mixture.Beijing Xinkechuang Tiger Nut Technology Development Co., Ltd.[12] used a spiral screening machine and a vibrating screen to complete conveying and soil screening during the harvesting of tiger nuts.The above studies show that the vibration principle can achieve effective separation of tiger nuts from soil during the conveying and separation process of tiger nut harvested materials.In addition, the vibration principle is also widely used in the separation of fruit and soil during the harvesting process of root crops such as Panax Notoginseng [13,14], potatoes [15][16][17], and peanuts [18][19][20].
On the basis of the above studies, in order to tackle problems such as difficult separation from sand and the high power consumption of tiger nut harvesting in the sandy areas of Xinjiang, we designed a conveying and separating device combining screw conveying and a vibrating screen, implementing a mixed state of vibration and spiral conveying during the conveying process.Compared with the single-vibration sand removal or spiralconveying threshing and harvesting devices, the sand removal efficiency is higher and the bean damage rate and power consumption are lower.The axial and radial migrations of materials under screw action and the separation process of materials under vibratory action were analyzed dynamically.Simulation analysis was carried out on the conveying and separation process based on EDEM software, and the conveying and separation law of the tiger nut harvest mixture was clarified.A four-factor, five-level orthogonal central composite test was performed to explore the influence law of the screw velocity ratio, amplitude, vibration frequency, machine operation velocity, and other test factors on the conveying and separation performance to determine the optimal operating parameters of the device so as to support the development of the technology and equipment for mechanized harvesting of tiger nuts in the sandy areas of Xinjiang.

Agronomy and Plant Characteristics of Tiger Nuts
Our research was conducted by studying the tiger nuts planted in sand by the 54th Regiment of the Third Division of the Xinjiang Production and Construction Corps.The variety is Zhongyousha No. 1, the distribution of the tiger nuts and their rhizomes is shown in Figure 1, and the planting method is flat planting, with a planting row spacing of 30~35 cm and a plant spacing of 5~7 cm.From the end of September to the beginning of October, tiger nuts were harvested.Before harvesting, the planting soil was cut off from the water supply for 7 to 10 days, and the tiger nut grass was harvested first, followed by the tiger nuts themselves.The tiger nut has a strong tillering regeneration capacity and develops roots in a clump-like shape.Their disorderly distributed fibrous root system is interwoven with sand, forming a composite structure of tiger nuts, roots, and sand (hereinafter referred to as "the complex").The characteristic parameters of tiger nuts were measured by the following measurement tools: a DL91150 digital vernier caliper (maximum measurement length: 150 mm, measurement accuracy: 0.01 mm, manufacturer: Deli Group Co., Ltd., Ningbo City, Zhejiang, China), a Sartorius MA100 electronic moisture analyzer (quality accuracy: 0.001 g, moisture content accuracy: 0.01%, manufacturer: Seydouris, Germany) and so on.The main characteristic parameters of mature tiger nuts are shown in Table 1.
Agriculture 2024, 14, x FOR PEER REVIEW 3 of 27 develops roots in a clump-like shape.Their disorderly distributed fibrous root system is interwoven with sand, forming a composite structure of tiger nuts, roots, and sand (hereinafter referred to as "the complex").The characteristic parameters of tiger nuts were measured by the following measurement tools: a DL91150 digital vernier caliper (maximum measurement length: 150 mm, measurement accuracy: 0.01 mm, manufacturer: Deli Group Co., Ltd., Ningbo City, Zhejiang, China), a Sartorius MA100 electronic moisture analyzer (quality accuracy: 0.001 g, moisture content accuracy: 0.01%, manufacturer: Seydouris, Germany) and so on.The main characteristic parameters of mature tiger nuts are shown in Table 1.The conveying and separating device is mainly composed of five-stage screw conveying device, front screen, back screen, exciter, exciting roller, return spring, collection box and rack, as shown in Figure 2a, and the sectional drawing of the conveying and separating device is shown in Figure 2b.The conveying and separating device is mainly composed of five-stage screw conveying device, front screen, back screen, exciter, exciting roller, return spring, collection box and rack, as shown in Figure 2a, and the sectional drawing of the conveying and separating device is shown in Figure 2b.

Working Principle
During operation, the rotary blade excavator tosses the excavated mixture of tiger nuts, roots, complexes, and a large amount of sand to the first-stage screw conveyor of the conveying and separating device.The screw blade pushes the mixture axially through augers in the same stage, so that the material moves along the peripheral direction of the screw blade while migrating axially.The mixture is conveyed centrally and tossed between adjacent augers through left-handed and right-handed rotation of blades in opposite directions.During rotation, the exciter interacts with the exciting roller, forcing the front screen and back screen to change from a static state to a state of moving away from the exciter.At the same time, the return spring exerts tension on the front screen and

Working Principle
During operation, the rotary blade excavator tosses the excavated mixture of tiger nuts, roots, complexes, and a large amount of sand to the first-stage screw conveyor of the conveying and separating device.The screw blade pushes the mixture axially through augers in the same stage, so that the material moves along the peripheral direction of the screw blade while migrating axially.The mixture is conveyed centrally and tossed between adjacent augers through left-handed and right-handed rotation of blades in opposite directions.During rotation, the exciter interacts with the exciting roller, forcing the front screen and back screen to change from a static state to a state of moving away from the exciter.At the same time, the return spring exerts tension on the front screen and back screen to pull them back to the original positions, thus leading to reciprocating periodic vibrations of the front screen and back screen.Under the combined action of screw conveying and the vibrating screen body, the mixture is in a loose state, which is conducive to the discharge of sand from the screen gaps and also conducive to the destruction of the structural state of the complex, promoting the separation of tiger nuts and roots from sand.In addition, the sand on the surface of the screen body is disturbed by spike teeth to clean the surface of the screen body and avoid blockage of the screen body caused by the accumulation of sand.A schematic diagram of the migration of complexes is shown in Figure 3.The technical parameters of the conveyance and separation test bench are shown in Table 2. Fifth-stage screw conveyor.8. Collection box.9. Back screen.10.Tire.11.Front screen.12. Exciting roller.13.Return Spring.14. Exciter.

Working Principle
During operation, the rotary blade excavator tosses the excavated mixture of tiger nuts, roots, complexes, and a large amount of sand to the first-stage screw conveyor of the conveying and separating device.The screw blade pushes the mixture axially through augers in the same stage, so that the material moves along the peripheral direction of the screw blade while migrating axially.The mixture is conveyed centrally and tossed between adjacent augers through left-handed and right-handed rotation of blades in opposite directions.During rotation, the exciter interacts with the exciting roller, forcing the front screen and back screen to change from a static state to a state of moving away from the exciter.At the same time, the return spring exerts tension on the front screen and back screen to pull them back to the original positions, thus leading to reciprocating periodic vibrations of the front screen and back screen.Under the combined action of screw conveying and the vibrating screen body, the mixture is in a loose state, which is conducive to the discharge of sand from the screen gaps and also conducive to the destruction of the structural state of the complex, promoting the separation of tiger nuts and roots from sand.In addition, the sand on the surface of the screen body is disturbed by spike teeth to clean the surface of the screen body and avoid blockage of the screen body caused by the accumulation of sand.A schematic diagram of the migration of complexes is shown in Figure 3.The technical parameters of the conveyance and separation test bench are shown in Table 2.A screw conveyor is an operating unit in the five-stage screw conveying device.It is mainly composed of a screw shaft, a screw blade, and screen-cleaning spike teeth and performs axial conveyance of the excavation mixture and radial migration among multistage screw conveyors.As the excavation mixture contains tiger nuts, sand, and other granular materials, a solid helical screw conveyor was selected, and its structure is shown in Figure 4.
mainly composed of a screw shaft, a screw blade, and screen-cleaning spike teeth and performs axial conveyance of the excavation mixture and radial migration among multistage screw conveyors.As the excavation mixture contains tiger nuts, sand, and other granular materials, a solid helical screw conveyor was selected, and its structure is shown in Figure 4. .By referring to the literature [21,22], according to Formula (1), it can be obtained that   ≥ 372 mm , 186 mm ≤   ≤ 818 mm , 74 mm ≤  0 ≤ 130 mm , and lead angle   = 33.37° .According to the recommended values in the Transport Machinery Type Selection Design Manual [23], let   = 400 m,   = 400 mm, and  0 = 108 mm.
In order to reduce the overall weight of the device, a hollow structure was designed for the screw shaft and the inner diameter of the screw shaft was set to be 104 mm.-outer diameter of the screw blade, mm;   -pitch of the screw blade, mm;  0 -outer diameter of the screw shaft, mm; K -characteristic coefficient of the material, The screw blade and screw shaft are the main components of the screw conveyor.The key structure and parameters of the screw conveyor were designed.The mass of sand in the excavation mixture was relatively large, so the sand density 2.65 t/m 3 was selected as the bulk density of the material.The feed amount of the conveying and separating device was calculated according to the moving velocity and working width of the machine to be Q = 100 t/h.By referring to the literature [21,22], according to Formula (1), it can be obtained that D S ≥ 372 mm, 186 mm ≤ S p ≤ 818 mm, 74 mm ≤ D 0 ≤ 130 mm, and lead angle φ a = 33.In order to reduce the overall weight of the device, a hollow structure was designed for the screw shaft and the inner diameter of the screw shaft was set to be 104 mm.
where D s -outer diameter of the screw blade, mm; S p -pitch of the screw blade, mm; D 0 -outer diameter of the screw shaft, mm; K-characteristic coefficient of the material, 0.05; Q-conveying capacity, t/h; ψ-filling coefficient, 0.25~0.35;ρ-bulk density, t/m 3 ; C-dip coefficient, 1.0 for a horizontal shaft; ϕ a -lead angle, • ; and ϕ f -frictional angle between sand and screw blade, 23.27 • .The screen-cleaning spike teeth were cylindrical steel teeth with a raised surface and a diameter of 8 mm.They were installed on both sides of the screw blade and uniformly arranged along the helical line, with a spacing of 20 mm.Considering that the particle size of tiger nuts ranges from 8 mm to 12 mm, and in order to avoid blockage of gaps by tiger nut root stubbles, the gaps between screen-cleaning spike teeth and screen body were designed to be 24 mm wide.Given the prerequisite of ensuring the transportation of tiger nuts in mixture and avoiding the phenomenon of damage, the operating line speed of the cleaning nail teeth should be less than 6.75 m/s, referring to the line speed of similar materials such as soybeans in the literature [24].Using Formula (2), it can be determined that the screw speed cannot exceed 280.8 r/min; that is, the maximum value of the screw speed in the later stages of this article is 280 r/min.
where v L -line speed of the screen-cleaning spike teeth, m/min; n-screw speed, r/min; and R q -rotary radius of the screen-cleaning spike teeth, m.

Analysis of the Material Conveying Process under Screw Action
(1) Analysis of the particle velocity during axial conveyance of materials In the process of conveying the complex, the screw conveying shaft rotates around the X-axis at a constant velocity w.The motion state of the complex at a distance r from the center of the screw shaft and on the screw blade surface was analyzed.Under the push and pull of the screw blade, the complex exhibited variable motion.The velocity analysis is shown in Figure 5.
where   -line speed of the screen-cleaning spike teeth, m/min; -screw speed, r/min; and   -rotary radius of the screen-cleaning spike teeth, m.

Analysis of the Material Conveying Process under Screw Action
(1) Analysis of the particle velocity during axial conveyance of materials In the process of conveying the complex, the screw conveying shaft rotates around the X-axis at a constant velocity w.The motion state of the complex at a distance r from the center of the screw shaft and on the screw blade surface was analyzed.Under the push and pull of the screw blade, the complex exhibited variable motion.The velocity analysis is shown in Figure 5.When the complex is not under a frictional action, its velocity is   perpendicular to the screw blade.When it is under a frictional action, its velocity deviates by β to become v [25,26].From the velocity analysis in Figure 3, it can be obtained that: where  0 -carrier velocity of the complex on the screw blade, m/s, and   -velocity of the complex perpendicular to the screw blade, m/s.
Meanwhile, the carrier velocity of the complex on the screw blade is: When the complex is not under a frictional action, its velocity is v n perpendicular to the screw blade.When it is under a frictional action, its velocity deviates by β to become v [25,26].From the velocity analysis in Figure 3, it can be obtained that: where v 0 -carrier velocity of the complex on the screw blade, m/s, and v n -velocity of the complex perpendicular to the screw blade, m/s.Meanwhile, the carrier velocity of the complex on the screw blade is: where r-distance of the complex from the center of the screw shaft, mm; w-angular velocity of the screw conveyor, r/min; and t-conveyance time, s.At this moment, the axial and circumferential components of the velocities v x and v t are: where α-lead angle, ( • ); β-frictional angle between the complex and the helical surface, ( • ); v x -axial velocity of the complex, m/s; and v-resultant velocity of the complex, m/s.Meanwhile, as f = tan β and the linear velocity of the edge of the screw blade is the highest, the criticality principle is applied.Suppose ∆r = r and substitute Formula (4) into Formula (5) to obtain the axial velocity of the complex: It can be found from the calculation that: Agriculture 2024, 14, 682 7 of 25 Substitute Formula (7) into Formula (6) to obtain: Here, f -frictional coefficient, 0.43; r-distance between the complex and the shaft center, 0.2 m; and S-screw blade pitch, 0.4 m.
When the speeds of the screw conveying shaft and the material are known, at the same radius r, the larger the screw blade pitch, the greater the movement velocity of the materials.To ensure that the complex can be conveyed axially, the axial velocity of complex particles must be greater than the circumferential velocity.From Formulae ( 5) and ( 8), the axial velocity and circumferential velocity are obtained, which are 0.81 m/s and 0.74 m/s, respectively, meeting the conditions of axial migration.
(2) Analysis of the particle velocity during radial tossing of materials In the process of screw conveying of materials, when the radial velocity of the materials increases, the materials will exhibit a tossing movement, and a tossing zone will be formed in the migration process of two adjacent screw conveyor stages.Due to the different arc lengths of the U-shaped conveyor channel, the complex will exhibit two kinds of motion: front tossing and rear tossing.The motion analysis is shown in Figure 6.As can be seen from Figure 6, under the action of the screw blades in opposite rotation directions, the axial velocity of the complex in the tossing zone is equal and opposite.Therefore, in the tossing process, the complex only has a radial velocity [27,28], which can be obtained as follows: In this case, when the angle δ of arc length AB of the U-shaped conveyor channel is less than or equal to 90°, a rear tossing motion occurs, and the motion trajectory is ABCDE.
When the angle δ of arc length AB of the U-shaped conveyor channel is greater than 90°, a front tossing motion occurs, and the motion trajectory is ABC′D′E′.As can be seen from Figure 6, under the action of the screw blades in opposite rotation directions, the axial velocity of the complex in the tossing zone is equal and opposite.Therefore, in the tossing process, the complex only has a radial velocity [27,28], which can be obtained as follows: In this case, when the angle δ of arc length AB of the U-shaped conveyor channel is less than or equal to 90 • , a rear tossing motion occurs, and the motion trajectory is ABCDE.
When the angle δ of arc length AB of the U-shaped conveyor channel is greater than 90 • , a front tossing motion occurs, and the motion trajectory is ABC From the above analysis, it can be seen that to ensure the effective conveying of materials between two screw conveyor stages, the angle of arc length AB of the U-shaped conveyor channel should be less than or equal to 90 • , and it was designed to be 65.5 • in this paper.

Structure and Parameter Design of the Vibrating Screen Body
According to the structural form of five-stage screw conveyors, the front screen and back screen of the vibrating screen body were designed as multi-stage curved screens, and the screen body adopted the structural form of a grid screen, as shown in Figure 7.The effective working width was designed to be 1620 mm, and the pores of a single screen were 8 mm × 4 mm.According to the gaps between the screen-cleaning spike teeth and the screen body and the working length of the screen-cleaning spike teeth determined in the previous stage, the radius of the grid screen arc was designed to be 254 mm.Considering the large number, length and load of screen bars, stiffeners were provided to position the screen bars and improve the stiffness of the screen bars.In order to improve the vibration frequency, the exciter was designed as a quad-arc cam structure with a symmetrical distribution of arcs, as shown in Figure 8.To ensure a smooth and stable contact process between the exciter and exciting roller and considering the processing difficulty and cost, the exciting roller was designed to have a circular outline.The parameters of the cam and exciting roller need to satisfy Formula ( 14).Then, based on the characteristics of the sand and the separation requirements, the design cam base circle radius is where   -cam base circle radius, mm;  ℎ1 -cam hub radius, mm;   -radius of the exciting roller, mm;  -radius of the exciting roller shaft, mm;  -center distance between cam and roller center, mm;   -maximum radius of convex profile, mm; and  ℎ2 -roller hub radius, mm.In order to improve the vibration frequency, the exciter was designed as a quad-arc cam structure with a symmetrical distribution of arcs, as shown in Figure 8.To ensure a smooth and stable contact process between the exciter and exciting roller and considering the processing difficulty and cost, the exciting roller was designed to have a circular outline.The parameters of the cam and exciting roller need to satisfy Formula (14).Then, based on the characteristics of the sand and the separation requirements, the design cam base circle radius is R b = 120 mm and the cam hub radius is R h1 = 24 mm.Let R r = 64 mm and let the roller shaft radius r = 21 mm.The roller hub radius was chosen as R h2 = 36 mm, and then let R max = 59 mm.
where R b -cam base circle radius, mm; R h1 -cam hub radius, mm; R r -radius of the exciting roller, mm; r-radius of the exciting roller shaft, mm; L-center distance between cam and roller center, mm; R max -maximum radius of convex profile, mm; and R h2 -roller hub radius, mm.

Analysis of the Material Separation Process under Vibration
When the motion mode of the screen body is vibration along the screen surface, as shown in Figure 9, the displacement, acceleration, and stress of sand particles were analyzed.The stress analysis of the sand particles on the vibrating screen surface shows that the direction of the inertia force of the sand particles was opposite to the vibration direction of the screen surface.At this moment, the screen body vibrated along the screen surface to the left and right, and the sand moved left and right along the screen surface under the action of inertia force.It is assumed that when the inertia force on the sand particles is leftward, the force on the sand particles at point C of the screen surface is as shown in Figure 9.When the screen body vibrates left and right along the screen surface, the vibration of the screen body can be regarded as a cosine acceleration motion of the cam-push rod mechanism (simple harmonic motion) [29].The displacement equation and acceleration equation are shown in Formula (15) and Formula (16).
where   -cam base circle radius, mm;  ℎ1 -cam hub radius, mm;   -radius of the exciting roller, mm;  -radius of the exciting roller shaft, mm;  -center distance between cam and roller center, mm;   -maximum radius of convex profile, mm; and  ℎ2 -roller hub radius, mm.When the motion mode of the screen body is vibration along the screen surface, as shown in Figure 9, the displacement, acceleration, and stress of sand particles were analyzed.The stress analysis of the sand particles on the vibrating screen surface shows that the direction of the inertia force of the sand particles was opposite to the vibration direction of the screen surface.At this moment, the screen body vibrated along the screen surface to the left and right, and the sand moved left and right along the screen surface under the action of inertia force.It is assumed that when the inertia force on the sand particles is leftward, the force on the sand particles at point C of the screen surface is as shown in Figure 9.When the screen body vibrates left and right along the screen surface, the vibration of the screen body can be regarded as a cosine acceleration motion of the cam-push rod mechanism (simple harmonic motion) [29].The displacement equation and acceleration equation are shown in Formula (15) and Formula (16).Through the analysis of Figure 9, it can be obtained that the force balance equation of sand particles on the screen surface [30] is:

Analysis of the Material Separation Process under Vibration
where Fns-normal component at point C of the arc surface for the force applied on the sand on the screen surface, N; Fa-inertia force of sand, N; G-gravity, N; θ-angle between the connecting line of point C on the screen surface, the circle center Os of the arc surface, and the horizontal line, rad; Fts-tangential component at point C on the arc surface for the force applied on the sand on the screen surface, N; Ffs-frictional force of sand on the screen surface, N; fs-friction factor of the screen surface against the sand; ms-mass of the sand, g; as-acceleration of vibrating screen, m/s 2 ; r-arc surface radius of the vibrating screen, m; hs-vertical distance between sand at point C of the screen surface and the bottom of the arc surface of the vibrating screen, m; and Fs-resultant force on the sand, N. Through the analysis of Figure 9, it can be obtained that the force balance equation of sand particles on the screen surface [30] is: where F ns -normal component at point C of the arc surface for the force applied on the sand on the screen surface, N; F a -inertia force of sand, N; G-gravity, N; θ-angle between the connecting line of point C on the screen surface, the circle center Os of the arc surface, and the horizontal line, rad; F ts -tangential component at point C on the arc surface for the force applied on the sand on the screen surface, N; F fs -frictional force of sand on the screen surface, N; f s -friction factor of the screen surface against the sand; m s -mass of the sand, g; a s -acceleration of vibrating screen, m/s 2 ; r-arc surface radius of the vibrating screen, m; h s -vertical distance between sand at point C of the screen surface and the bottom of the arc surface of the vibrating screen, m; and F s -resultant force on the sand, N.
After substituting Formula (16) into Formula (15), it was rearranged to obtain: According to Formulas ( 17) and ( 18), when F ns ≤ 0, the sand particles are on the surface of the vibrating screen.If F ts > 0, the sand will slide along the screen surface to the right of the arc surface of the screen body; if F ts < 0, the sand will slide along the screen surface to the left of the arc surface of the screen body.When F ns > 0, the sand is tossed diagonally upward by the vibrating screen.If F ts > 0, the sand will be tossed up by the screen body to the right of the screen body normal; if F ts < 0, the sand will be tossed up by the screen body to the left of the screen body normal.To ensure the fluidity and separability of sand, it is necessary to appropriately improve the vibration frequency and amplitude of the vibrating screen.When selecting the amplitude and frequency parameters, a large amplitude and a low frequency should be chosen for materials with a larger particle size, and, vice versa, a small amplitude and a high frequency should be chosen for materials with a smaller particle size [31].Since the screening materials mainly consist of loose materials such as sand and soil, paper [32] concluded that the optimal screening efficiency of the vibrating screen is achieved with an amplitude of 4 mm, while paper [33] found that the optimal screening efficiency is achieved with an amplitude of 5 mm.Therefore, the preliminary amplitude value is set at 4.5 mm, and experiments will be conducted with this value as the intermediate level for selecting test variable factors.In the tiger nut and castor bean cleaning process, paper [34,35] used 7 Hz as the intermediate level for vibration frequency testing, while paper [36] found that in the process of separating potato soil, the maximum vibration frequency was 9 Hz.Considering factors such as a smaller mesh size and a finer sand particle size, and in accordance with the sand separation requirements of the tiger nut harvesting process, the maximum vibration frequency parameter is set at 9 Hz, with 7 Hz as the intermediate level for selecting test variable factors in later experiments.

Setting of Simulation Parameters and Model Construction
In order to verify the rationality of the design of the conveying and separating device, EDEM software was used for simulation analysis of the conveying and separating process.As the motion process and the forms of the complex and sand mass in the device are complex, and the number of simulated particles included in the material modeling is large; thus, a Hertz-Mindlin contact model was used.Referring to our previous article [37] and based on the determined parameters of tiger nuts, sand, the complex, and other materials, as well as the literature [38][39][40][41][42], the intrinsic parameters (Poisson's ratio, shear modulus, etc.) and the basic contact parameters (collision recovery coefficient, static frictional coefficient and rolling frictional coefficient, etc.) of sand particles, tiger nuts, tiger nut roots, and the geometry of the sand conveying and removing device were set.The intrinsic simulation parameters and contact parameters are shown in Tables 3 and 4.
According to the tiger nut planting pattern on the sandy land and the overall working width of the harvester, a soil tank model was constructed in EDEM 2021.2 software.EDEM software is an advanced engineering simulation software that uses the discrete element method to simulate the particle flow, particle dynamics, and particle material interactions.It has many functions such as particle simulations, multi-physical simulations, visualization, and analysis.EDEM software is a powerful and widely used particle dynamics simulation software in the engineering field, and can help engineers and researchers better understand and optimize the behavior of particle systems.The software version used in this work is EDEM 2021.2.Multiple small squares were generated in the working area according to the tiger nut planting pattern parameters.Supposing each small square is a hole of a tiger nut plant, the location of the tiger nut-root-sand complexes was determined.Then, the tiger nut-root-sand complex and sand particle models established earlier were introduced so that sand particles could fill up the working area and bury and cover the complex.After the settlement of sand particles was complete, the soil tank model was finally established, as shown in Figure 10a.Solidworks three-dimensional software was used to simplify the conveying and separating device.The simplified model was saved in igs format and imported into EDEM software.The total simulation time was 12 s and the step length was 0.02 s.A simulation model was established, as shown in Figure 10b.01 According to the tiger nut planting pattern on the sandy land and the overall working width of the harvester, a soil tank model was constructed in EDEM 2021.2 software.EDEM software is an advanced engineering simulation software that uses the discrete element method to simulate the particle flow, particle dynamics, and particle material interactions.It has many functions such as particle simulations, multi-physical simulations, visualization, and analysis.EDEM software is a powerful and widely used particle dynamics simulation software in the engineering field, and can help engineers and researchers better understand and optimize the behavior of particle systems.The software version used in this work is EDEM 2021.2.Multiple small squares were generated in the working area according to the tiger nut planting pattern parameters.Supposing each small square is a hole of a tiger nut plant, the location of the tiger nut-root-sand complexes was determined.Then, the tiger nut-root-sand complex and sand particle models established earlier were introduced so that sand particles could fill up the working area and bury and cover the complex.After the settlement of sand particles was complete, the soil tank model was finally established, as shown in Figure 10a.Solidworks threedimensional software was used to simplify the conveying and separating device.The simplified model was saved in igs format and imported into EDEM software.The total simulation time was 12 s and the step length was 0.02 s.A simulation model was established, as shown in Figure 10b.

Analysis of Sand Particle Movements and Displacement Variation
The motion process of the sand particles under the double actions of screw conveying and screen body vibration was analyzed.The mixture is under the actions of the screw blade and screen-cleaning spike teeth in the process of screw conveying and is spread on the surface of the screen body.At the same time, the structural state of the mixture changed under the vibratory action of the screen body, accelerating the separation

Analysis of Sand Particle Movements and Displacement Variation
The motion process of the sand particles under the double actions of screw conveying and screen body vibration was analyzed.The mixture is under the actions of the screw blade and screen-cleaning spike teeth in the process of screw conveying and is spread on the surface of the screen body.At the same time, the structural state of the mixture changed under the vibratory action of the screen body, accelerating the separation between tiger nuts and sand in the mixture, and the sand was discharged from the gaps between the screen bars so as to realize the effect of separating sand and discharging it out of the machine while conveying and migrating tiger nuts, as shown in Figure 11.In order to gain a more intuitive understanding of the movement state of sand during the mixture conveying process, sand particles were randomly selected in the action area of the screen-cleaning spike teeth and the screw blade conveying area, respectively, and the variations in the displacement parameters in the X-axis, Y-axis, and Z-axis over time were extracted.
between tiger nuts and sand in the mixture, and the sand was discharged from the gaps between the screen bars so as to realize the effect of separating sand and discharging it out of the machine while conveying and migrating tiger nuts, as shown in Figure 11.In order to gain a more intuitive understanding of the movement state of sand during the mixture conveying process, sand particles were randomly selected in the action area of the screen-cleaning spike teeth and the screw blade conveying area, respectively, and the variations in the displacement parameters in the X-axis, Y-axis, and Z-axis over time were extracted.

Analysis of the Displacement Variation of Sand Particles along the X-axis
The X-axis direction is the direction of advancement along the screw shaft.As can be seen from Figures 12a and 13a, during the movement of sand particles in the action area of the screen-cleaning spike teeth in the reverse direction of the X-axis, the displacement variation amplitude was small at within 0.25~1.70s and increased between 1.70 s and 2.05 s probably because of the indirect contact between the screen-cleaning spike teeth and the sand particles.With the disturbance of the screen-cleaning spike teeth, the displacement of the sand particles under the applied force changed and the sand particles were discharged from the screen body gaps.After 2.55 s, the displacement in the X-axis direction no longer changed.The X-axis direction is the direction of advancement along the screw shaft.As can be seen from Figures 12a and 13a, during the movement of sand particles in the action area of the screen-cleaning spike teeth in the reverse direction of the X-axis, the displacement variation amplitude was small at within 0.25~1.70s and increased between 1.70 s and 2.05 s probably because of the indirect contact between the screen-cleaning spike teeth and the sand particles.With the disturbance of the screen-cleaning spike teeth, the displacement of the sand particles under the applied force changed and the sand particles were discharged from the screen body gaps.After 2.55 s, the displacement in the X-axis direction no longer changed.
out of the machine while conveying and migrating tiger nuts, as shown in Figure 11.In order to gain a more intuitive understanding of the movement state of sand during the mixture conveying process, sand particles were randomly selected in the action area of the screen-cleaning spike teeth and the screw blade conveying area, respectively, and the variations in the displacement parameters in the X-axis, Y-axis, and Z-axis over time were extracted.

Analysis of the Displacement Variation of Sand Particles along the X-axis
The X-axis direction is the direction of advancement along the screw shaft.As can be seen from Figures 12a and 13a, during the movement of sand particles in the action area of the screen-cleaning spike teeth in the reverse direction of the X-axis, the displacement variation amplitude was small at within 0.25~1.70s and increased between 1.70 s and 2.05 s probably because of the indirect contact between the screen-cleaning spike teeth and the sand particles.With the disturbance of the screen-cleaning spike teeth, the displacement of the sand particles under the applied force changed and the sand particles were discharged from the screen body gaps.After 2.55 s, the displacement in the X-axis direction no longer changed.out of the machine while conveying and migrating tiger nuts, as shown in Figure 11.In order to gain a more intuitive understanding of the movement state of sand during the mixture conveying process, sand particles were randomly selected in the action area of the screen-cleaning spike teeth and the screw blade conveying area, respectively, and the variations in the displacement parameters in the X-axis, Y-axis, and Z-axis over time were extracted.

Analysis of the Displacement Variation of Sand Particles along the X-axis
The X-axis direction is the direction of advancement along the screw shaft.As can be seen from Figures 12a and 13a, during the movement of sand particles in the action area of the screen-cleaning spike teeth in the reverse direction of the X-axis, the displacement variation amplitude was small at within 0.25~1.70s and increased between 1.70 s and 2.05 s probably because of the indirect contact between the screen-cleaning spike teeth and the sand particles.With the disturbance of the screen-cleaning spike teeth, the displacement of the sand particles under the applied force changed and the sand particles were discharged from the screen body gaps.After 2.55 s, the displacement in the X-axis direction no longer changed.The displacement variation of sand particles in the screw action area over time was caused by the movement change trends in the forward direction and reverse direction of the X-axis.The displacement of sand particles showed variation in the forward direction of the X-axis within 0~1.65 s, and showed variation in the reverse direction of the X-axis after 1.65 s probably because the sand particles were tossed from the conveyor auger on one side to the conveyor auger on the opposite side and then were conveyed to the ground in the opposite auger, while the displacement in the X-axis direction no longer changed.

Analysis of Displacement Variation of Sand Particles along the Y-axis
The Y-axis direction is along the advancing direction of the machine.It can be seen from Figures 12b and 13b that the sand particles in the action area of the screen-cleaning spike teeth and the screw action area all moved in the forward direction of the Y-axis with the change in time.The displacement variation amplitude of the sand particles in the action area of the screen-cleaning spike teeth first increased gently within 0~1.80 s, then increased sharply from 1.80 s to 2.50 s, and finally became stable and no longer changed after 2.50 s.
Within 0~1.60 s, the displacement of sand particles in the screw action area was consistent with the displacement in the action area of the screen-cleaning spike teeth; that is, there was first a gentle increase and then sharp increase probably because the sand particles were first transferred from the bottom end of the screw blade to the lateral end under the drive of the screw blade, and then were transferred from the lateral end of the blade and to the auger in the next stage under the action of centrifugal force, passing through the gaps of the screen body and finally falling to the ground.However, the displacement decreased after 1.60 s probably because the sand particles rebounded after they collided with and contacted the auger blades, rack, and other components.

Analysis of Displacement Variation of Sand Particles along the Z-axis
In the figure, the Z-axis is perpendicular to the ground, where the forward direction of the Z-axis is away from the ground and the reverse direction of the Z-axis approaches the ground.From Figures 12c and 13c, it can be seen that the sand particles in the action area of the screen-cleaning spike teeth moved in the reverse direction of the Z-axis with the change in time, and the displacement variation gently decreased within 0.25~1.90s, probably because there was still a sand layer between the sand particles and the surface of the screen body or the gaps of the screen body were blocked.With the movement of the screen-cleaning spike teeth, the structure of the sand layer was destroyed or the sand in the gaps of the screen body was cleared, and the sand particles were discharged to the ground through the gaps in the screen body.After 2.60 s, the displacement in the Z-axis direction no longer changed.
The displacement of the sand particles in the screw action area showed change trends of movement in the reverse direction of the Z-axis, forward direction of the Z-axis, and reverse direction of the Z-axis over time.To be specific, the displacement of the sand particles changed in the forward direction of the Z-axis within 1.45~1.60s, probably because the sand particles moved towards the forward direction of the Z-axis under the vibration of the screen body.When they reached a certain position, under the action of gravity, the sand fell back to the screen surface in the reverse direction of the Z-axis and was discharged to the ground through the gaps in the screen body.After 2.60 s, the displacement in the Z-axis direction no longer changed.
In summary, it can be seen that the displacement of the sand particles in the action area of the screen-cleaning spike teeth and the screw action area in the Y-axis direction changed more significantly than in the X-axis direction and the Z-axis direction; that is, the sand could be conveyed stage by stage and out of the machine, finally falling to the ground.At the same time, it can be seen from the displacement variation in the X-axis direction that the sand particles could also be conveyed axially in the augers of the same stage, thus realizing the removal of the sand.

Test Materials and Equipment
In April 2023, experimental research on conveying and separating devices was carried out in the Soil Tank Performance Laboratory, the School of Mechanical and Electrical Engineering, Shihezi University.The tiger nut-root-sand complex material was spread in the soil tank test area according to the planting pattern.The test instruments and equipment included a soil tank test vehicle (TCC-3.0),an electronic scale (500 g, 100 kg), a three-phrase motor (type: YE3-132S-4, power: 5.5 KW, revolution speed: 1469 r/min), an HCNJ-101 dynamic torque sensor (range: 0~500 KN× m, accuracy: ±0.25%, Supplier: Beijing Haibohua Technology Co., Ltd., Beijing, China), a tape measure, a spade, a woven bag, and a Ziplock bag and marker.The device for the test was a self-built tiger nut harvest conveying and separating device.This device was equipped with a control system with the functions of single-stage auger speed adjustment and single-stage auger conveyance torque acquisition and was fitted with a sand collection tank, which can collect the sand discharged by the augers at all stages.The test site is shown in Figure 14.
In April 2023, experimental research on conveying and separating devices was carried out in the Soil Tank Performance Laboratory, the School of Mechanical and Electrical Engineering, Shihezi University.The tiger nut-root-sand complex material was spread in the soil tank test area according to the planting pattern.The test instruments and equipment included a soil tank test vehicle (TCC-3.0),an electronic scale (500 g, 100 kg), a three-phrase motor (type: YE3-132S-4, power: 5.5 KW, revolution speed: 1469 r/min), an HCNJ-101 dynamic torque sensor (range: 0~500 KN× m, accuracy: ±0.25%, Supplier: Beijing Haibohua Technology Co., Ltd., Beijing, China), a tape measure, a spade, a woven bag, and a Ziplock bag and marker.The device for the test was a self-built tiger nut harvest conveying and separating device.This device was equipped with a control system with the functions of single-stage auger speed adjustment and single-stage auger conveyance torque acquisition and was fitted with a sand collection tank, which can collect the sand discharged by the augers at all stages.The test site is shown in Figure 14.

Test Scheme
Considering that the function of the conveyor augers in the first three stages is mainly to remove sand and that the mass of the corresponding mixture shows a decreasing trend, the speed of the screw conveyor in the first three stages was designed to decrease; that is, the speed of the screw conveyor in the first stage is the highest, followed by the speed in the second stage, and it is smallest in the third stage.At the same time, the speeds of the screw conveyor in the third, fourth, and fifth stages were set to be consistent.The speed ratio among the screw conveyor in the first, second, and third stages was set to be 1, 0.95, 0.9, 0.85, and 0.8, as shown in Table 5.

Test Scheme
Considering that the function of the conveyor augers in the first three stages is mainly to remove sand and that the mass of the corresponding mixture shows a decreasing trend, the speed of the screw conveyor in the first three stages was designed to decrease; that is, the speed of the screw conveyor in the first stage is the highest, followed by the speed in the second stage, and it is smallest in the third stage.At the same time, the speeds of the screw conveyor in the third, fourth, and fifth stages were set to be consistent.The speed ratio among the screw conveyor in the first, second, and third stages was set to be 1, 0.95, 0.9, 0.85, and 0.8, as shown in Table 5.With key parameters such as the screw velocity ratio, amplitude, vibration frequency, and machine operation velocity as test factors, and the sand removal rate, crushing rate, and power consumption as test evaluation indicators, an orthogonal test was performed.The test factor levels and coding values are shown in Table 6.The sand removal performance was described based on paper [43] with the sand removal rate as an evaluation indicator.Supposing that the working width of the conveying and separating device is 1.6 m and the working depth is 0.15 m, the sand feed amount was theoretically calculated in m 2 , which is the target sand removal amount.The total mass m 1 of sand discharged by the augers at all stages in the collection tank was weighed, which is the actual sand removal amount.The sand removal rate [44] was calculated according to the formula: In Formula ( 19), R 1 is the sand removal rate, %; m 1 is the actual sand removal amount, kg; and m 2 is the target sand removal amount, kg.
(2) Crushing rate We referred to the national standard GB/T 21962−2020 and the literature [45] and used the tiger nut crushing rate as an evaluation indicator to describe the conveying performance.The nuts at the discharge port were collected and weighed (m 3 ), crushed nuts were removed and their mass was weighed (m 4 ), and the crushing rate was calculated according to the formula: In Formula (20), R 2 is the crushing rate, %; m 3 is the mass of all nuts at the discharge port, g; and m 2 is the mass of crushed nuts at the discharge port, g.
(3) Power consumption With power consumption as an evaluation indicator, the conveying and sand removal performance was described.The torque values of the augers from torque sensors were obtained and the power consumption was calculated according to the following formula in combination with the set rotation speed of the screw shaft.
In Formula ( 21), R 3 is the power consumption of the device in the working process, kW; T is the torque of all conveyor augers, Nm; and n is the speed of the screw shaft, r/min.

Test Design
According to the test factors and evaluation indicators set in the early stage, a fourfactor five-level orthogonal test design was used in the Central Composite Design (CCD) module in the experimental design and analysis software Design Expert [46].A total of 30 groups of tests were conducted, including 24 groups of analysis factor tests and 6 groups of center zero estimation error tests.Each group of tests was repeated three times, and the average value was calculated.During the test, the conveying and separating device was dragged by the soil tank test bench and the operating speed was adjusted by the soil tank test bench.The screw auger rotation speed and vibration frequency parameters were adjusted by the test device control system.The amplitude was adjusted by using vibration exciters of different sizes.The test arrangement and results are shown in Table 7.The results of the regression variance analysis of conveying and sand removal quality are shown in Table 8.It can be seen from the table that in the regression variance analysis of evaluation indicator R 1 , the p values of single factor X 2 and interaction terms X 3 X 4 , X 1 2 , X 2 2 , and X 4 2 are 0.0061, 0.0014, <0.0001, <0.0001, and 0.0014, respectively, much less than 0.01, and are highly significant influence factors of R 1 .The p values of single factors X 3 and X 4 and interaction terms X 2 X 4 and X 3 2 are 0.0171, 0.0338, 0.0375, and 0.0489, respectively, all within the range of 0.01~0.05,and are significant influence factors of R 1 .The p values of interaction terms X 1 X 2 , X 1 X 3 , X 1 X 4 , and X 2 X 3 are all greater than 0.05, and are insignificant influence factors of R 1 .After ignoring the insignificant influence factors of R 1 , the highly significant and significant influence factors of R 1 were sorted according to their p values as: Note: ** = (p < 0.01), * = (0.01 < p < 0.05), not significant = p > 0.05.
As can be seen from Table 8, in the regression variance analysis of evaluation indicator R 2 , the p values of single factors X 1 , X 3 , and X 4 and interaction terms X 1 X 2 , X 1 2 , and X 2 2 are 0.0096, 0.0002, 0.0021, 0.0063, 0.0019, and 0.0094, respectively, all much less than 0.01, and are highly significant influence factors of R 2 .The p value of the interaction term X 1 X 4 is 0.0139, within the range of 0.01~0.05,and is a significant influence factor of R 2 .The p values of the single factor X 2 and interaction terms X 1 X 3 , X 2 X 3 , X 2 X 4 , X 3 X 4 , X 3 2 , and X 4 2 are all greater than 0.05, and are insignificant influence factors of R 2 .After ignoring the insignificant influence factors of R 2 , the highly significant and significant influence factors of R 2 were sorted according to p values as: As can be seen from Table 8, in the regression variance analysis of evaluation indicator R 3 , the p values of single factors X 3 and X 4 and interaction terms X 1 X 2 and X 1 X 4 are < 0.0001, <0.0001, 0.0004, and <0.00010, respectively, all much less than 0.01, and are highly significant influence factors.The p values of interaction terms X 3 X 4 , X 1 2 , and X 2 2 are 0.0419, 0.0273, and 0.045, respectively, within the range of 0.01~0.05,and are significant influence factors of R 3 .The p values of single factors X 1 and X 2 and interaction terms X 1 X 3 , X 2 X 3 , X 2 X 4 , X 3 2 , and X 4 2 are all greater than 0.05, and are insignificant influence factors of R 3 .After ignoring the insignificant influence factors of R 3 , the highly significant and significant influence factors of R 1 were sorted according to p values as: Based on the results of the regression variance analysis, a regression model equation of conveying and sand removal evaluation indicators was established using the sand removal rate R 1 , crushing rate R 2 , and total power consumption R 3 against test independent variables such as the screw velocity ratio X 1 , amplitude X 2 , vibration frequency X 3 , and machine operation velocity X 4 .
The significance of the influence of each variable in the regression equation on the indicators was determined by the F-test, and the smaller the probability p, the higher the significance of the corresponding variable.The p values of regression models R 1 , R 2 , and R 3 are 0.0001, 0.0004, and <0.0001, respectively, all less than 0.01, indicating that the constructed regression models are highly significant, and the p values of lack-of-fit terms in the models are 0.0636, 0.3355, and 0.0944, respectively, all greater than 0.05; that is, the lack-of-fit terms are not significant.This indicates that the predicted values of the regression models have a good fit with the actual values within the test range.At the same time, the coefficients of determination of the models are 0.8852, 0.8596, and 0.911, respectively, and the coefficients of variation are 2.69%, 21.37%, and 8.44%, respectively, indicating that the regression models can explain the test data well, the predicted values are highly correlated with the actual values, and the constructed regression models can be used to analyze and predict the test evaluation indicators of conveying and sand removal quality.

Analysis of the Influence of the Interaction of Key Parameters of the Conveying Device on Conveying and Sand Removal Quality
As can be seen from Table 8: Results of regression variance analysis, the interaction factors with a highly significant or significant influence on evaluation indicator R 1 are X 2 X 4 and X 3 X 4 , respectively; the factors with a highly significant or significant influence on evaluation indicator R 2 are X 1 X 2 and X 1 X 4 , respectively; and the factors with a highly significant or significant influence on evaluation indicator R 3 are X 1 X 2 , X 1 X 4 , and X 3 X 4 , respectively.In the process of analyzing the influence of the interaction of key parameters of the device on the test evaluation indicators, the insignificant interaction factors were not considered.Only the influence law between the interaction factors with a highly significant or significant influence and the evaluation indicators of conveying and sand removal quality was analyzed, and the response surface diagram of each interaction term was extracted.The analysis of the influence law of interaction term X 2 X 4 on the evaluation indicator of the sand removal rate R 1 is shown in Figure 15a.When amplitude X 2 is 1.5 mm and the machine operation velocity increases from 0.6 km/h to 1.4 km/h, the sand removal rate shows a trend of a gradual increase.When the amplitude is 7.5 mm and the machine operation velocity increases from 0.6 km/h to 1.4 km/h, the sand removal rate shows a trend of first increasing then decreasing.It can be seen that when the amplitude is the largest, the sand removal rate can be increased by increasing the machine operation velocity.
increasing and then decreasing, indicating that when the machine operation velocity is the smallest or largest, its influence of the amplitude on the sand removal rate is significant.
The analysis of the influence law of interaction term X3X4 on the evaluation indicator of the sand removal rate R1 is shown in Figure 15b.When the vibration frequency X3 is 5 Hz and the machine operation velocity X4 increases from 0.6 km/h to 1.4 km/h, the sand removal rate shows a trend of a gradual increase.When the vibration frequency is 9 Hz and the machine operation velocity X4 increases from 0.6 km/h to 1.4 km/h, the sand removal rate shows a trend of a gradual decrease.To sum up, when the vibration frequency is small, the sand removal rate can be increased by increasing the machine operation velocity.When the vibration frequency is large, the sand removal rate can be increased by reducing the machine operation velocity.
When the machine operation velocity is 0.6 km/h and the vibration frequency increases from 5 Hz to 9 Hz, the sand removal rate shows a trend of a gradual increase.When the machine operation velocity is 1.4 km/h and the vibration frequency increases from 5 Hz to 9 Hz, the sand removal rate shows a trend of a gradual decrease.Thus, it can be seen that when the machine operation velocity is the largest, the sand removal rate can be increased by reducing the vibration frequency.The analysis of the influence law of interaction term X1X2 on the sand removal rate evaluation indicator R2 is shown in Figure 16a.When the screw velocity ratio X1 is 0.8 and the amplitude X2 increases from 1.5 mm to 7.5 mm, the crushing rate shows a trend of a gradual increase.When the screw velocity ratio is 1 and the amplitude increases from 1.5 mm to 7.5 mm, the crushing rate shows a trend of a gradual decrease.To sum up, when the screw velocity ratio is small, the crushing rate can be reduced by reducing the amplitude.When the screw velocity ratio is large, the crushing rate can be reduced by increasing the amplitude.
When the amplitude is 1.5 mm and the screw velocity ratio increases from 0.8 to 1, the crushing rate shows a trend of a gradual increase.Thus, it can be seen that when the amplitude is small, the crushing rate can be reduced by reducing the amplitude and screw velocity ratio.When the amplitude is 7.5 mm and the screw velocity ratio increases from 0.8 to 1, the crushing rate shows a trend of first decreasing and then increasing.When the machine operation velocity X 4 is 0.6 km/h and the amplitude increases from 1.5 mm to 7.5 mm, the sand removal rate shows a trend of first increasing and then decreasing.When the machine operation velocity X 4 is 1.4 km/h and the amplitude increases from 1.5 mm to 7.5 mm, the sand removal rate also shows a trend of first increasing and then decreasing, indicating that when the machine operation velocity is the smallest or largest, its influence of the amplitude on the sand removal rate is significant.
The analysis of the influence law of interaction term X 3 X 4 on the evaluation indicator of the sand removal rate R 1 is shown in Figure 15b.When the vibration frequency X 3 is 5 Hz and the machine operation velocity X 4 increases from 0.6 km/h to 1.4 km/h, the sand removal rate shows a trend of a gradual increase.When the vibration frequency is 9 Hz and the machine operation velocity X 4 increases from 0.6 km/h to 1.4 km/h, the sand removal rate shows a trend of a gradual decrease.To sum up, when the vibration frequency is small, the sand removal rate can be increased by increasing the machine operation velocity.When the vibration frequency is large, the sand removal rate can be increased by reducing the machine operation velocity.
When the machine operation velocity is 0.6 km/h and the vibration frequency increases from 5 Hz to 9 Hz, the sand removal rate shows a trend of a gradual increase.When the machine operation velocity is 1.4 km/h and the vibration frequency increases from 5 Hz to 9 Hz, the sand removal rate shows a trend of a gradual decrease.Thus, it can be seen that when the machine operation velocity is the largest, the sand removal rate can be increased by reducing the vibration frequency.

Analysis of the Influence Law of Significant Interaction Terms on the Evaluation Indicator Crushing Rate R 2
The analysis of the influence law of interaction term X 1 X 2 on the sand removal rate evaluation indicator R 2 is shown in Figure 16a.When the screw velocity ratio X 1 is 0.8 and the amplitude X 2 increases from 1.5 mm to 7.5 mm, the crushing rate shows a trend of a gradual increase.When the screw velocity ratio is 1 and the amplitude increases from 1.5 mm to 7.5 mm, the crushing rate shows a trend of a gradual decrease.To sum up, when the screw velocity ratio is small, the crushing rate can be reduced by reducing the amplitude.When the screw velocity ratio is large, the crushing rate can be reduced by increasing the amplitude.and the machine operation velocity X4 increases from 0.6 km/h to 1.4 km/h, the crush rate shows a trend of a gradual increase.Thus, it can be seen that when the screw veloc ratio is small, the crushing rate can be reduced by reducing the machine operat velocity.When the screw velocity ratio is 1 and the machine operation velocity increa from 0.6 km/h to 1.4 km/h, the crushing rate shows a trend of a gradual decrea indicating that when the screw velocity ratio is large, the crushing rate can be reduced increasing the machine operation velocity.
When the machine operation velocity is 0.6 km/h and the screw velocity ra increases from 0.8 to 1, the crushing rate shows a trend of a gradual increase.When machine operation velocity is 1.4 km/h and the screw velocity ratio increases from 0.8 1, the crushing rate shows a trend of first decreasing and then increasing.The analysis of the influence law of interaction term X1X2 on the total pow consumption evaluation indicator R3 is shown in Figure 17a.When the screw velocity ra X1 is 0.8 and the amplitude X2 increases from 1.5 mm to 7.5 mm, the total pow consumption shows a trend of a gradual decrease.When the screw velocity ratio is 1 a the amplitude increases from 1.5 mm to 7.5 mm, the total power consumption show trend of a gradual increase.To sum up, when the screw velocity ratio is small, the to power consumption can be reduced by increasing the amplitude.When the screw veloc ratio is large, the total power consumption can be reduced by reducing the amplitude When the amplitude is 1.5 mm and the screw velocity ratio increases from 0.8 to 1, the crushing rate shows a trend of a gradual increase.Thus, it can be seen that when the amplitude is small, the crushing rate can be reduced by reducing the amplitude and screw velocity ratio.When the amplitude is 7.5 mm and the screw velocity ratio increases from 0.8 to 1, the crushing rate shows a trend of first decreasing and then increasing.
The analysis of the influence law of interaction term X 1 X 4 on the evaluation indicator of the sand removal rate R 2 is shown in Figure 16b.When the screw velocity ratio X 1 is 0.8 and the machine operation velocity X 4 increases from 0.6 km/h to 1.4 km/h, the crushing rate shows a trend of a gradual increase.Thus, it can be seen that when the screw velocity ratio is small, the crushing rate can be reduced by reducing the machine operation velocity.When the screw velocity ratio is 1 and the machine operation velocity increases from 0.6 km/h to 1.4 km/h, the crushing rate shows a trend of a gradual decrease, indicating that when the screw velocity ratio is large, the crushing rate can be reduced by increasing the machine operation velocity.
When the machine operation velocity is 0.6 km/h and the screw velocity ratio increases from 0.8 to 1, the crushing rate shows a trend of a gradual increase.When the machine operation velocity is 1.4 km/h and the screw velocity ratio increases from 0.8 to 1, the crushing rate shows a trend of first decreasing and then increasing.

Analysis of the Influence Law of Significant Interaction Terms on the Total Power Consumption Evaluation Indicator, R 3
The analysis of the influence law of interaction term X 1 X 2 on the total power consumption evaluation indicator R 3 is shown in Figure 17a.When the screw velocity ratio X 1 is 0.8 and the amplitude X 2 increases from 1.5 mm to 7.5 mm, the total power consumption shows a trend of a gradual decrease.When the screw velocity ratio is 1 and the amplitude increases from 1.5 mm to 7.5 mm, the total power consumption shows a trend of a gradual increase.To sum up, when the screw velocity ratio is small, the total power consumption can be reduced by increasing the amplitude.When the screw velocity ratio is large, the total power consumption can be reduced by reducing the amplitude.
When the amplitude is 1.5 mm and the screw velocity ratio increases from 0.8 to 1, the total power consumption shows a trend of a gradual decrease.When the amplitude is 7.5 mm and the screw velocity ratio increases from 0.8 to 1, the total power consumption shows a trend of a gradual increase.To sum up, when the amplitude is small, the total power consumption can be reduced by increasing the screw velocity ratio.When the amplitude is large, the total power consumption can be reduced by reducing the screw velocity ratio.

Consumption Evaluation Indicator, R3
The analysis of the influence law of interaction term X1X2 on the total power consumption evaluation indicator R3 is shown in Figure 17a.When the screw velocity ratio X1 is 0.8 and the amplitude X2 increases from 1.5 mm to 7.5 mm, the total power consumption shows a trend of gradual decrease.When the screw velocity ratio is 1 and the amplitude increases from 1.5 mm to 7.5 mm, the total power consumption shows a trend of a gradual increase.To sum up, when the screw velocity ratio is small, the total power consumption can be reduced by increasing the amplitude.When the screw velocity ratio is large, the total power consumption can be reduced by reducing the amplitude.In short, when the screw velocity ratio is small and the amplitude is large, another factor in the interaction term needs to be increased or reduced to obtain a lower power consumption.When the screw velocity ratio increases from 0.8 to 1 and the amplitude increases from 1.5 mm to 7.5 mm, the total power consumption shows a trend of first decreasing and then increasing.
The analysis of the influence law of interaction term X 1 X 4 on the total power consumption evaluation indicator, R 3 , is shown in Figure 17b.When the screw velocity ratio X 1 is 0.8 and the machine operation velocity X 4 increases from 0.6 km/h to 1.4 km/h, the total power consumption shows a trend of a gradual increase.When the screw velocity ratio is 1 and the machine operation velocity increases from 0.6 km/h to 1.4 km/h, the total power consumption shows a trend of a gradual decrease.To sum up, when the screw velocity ratio is small, the operating power consumption can be reduced by reducing the machine operation velocity.When the screw velocity ratio is large, the operating power consumption can be reduced by increasing the machine operation velocity.
When the machine operation velocity is 0.6 km/h and the screw velocity ratio increases from 0.8 to 1, the total power consumption shows a trend of a gradual increase.When the machine operation velocity is 1.4 km/h and the screw velocity ratio increases from 0.8 to 1, the total power consumption shows a trend of a gradual decrease.To sum up, when the machine operation velocity is small, the total power consumption can be reduced by reducing the screw velocity ratio.When the machine operation velocity is large, the total power consumption can be reduced by increasing the screw velocity ratio.
In short, when the screw velocity ratio is small and the machine operation velocity is large, another factor in the interaction term needs to be reduced or increased to obtain a lower power consumption value.When the screw velocity ratio increases from 0.8 to 1 and the machine operation velocity increases from 0.6 km/h to 1.4 km/h, the total power consumption shows a trend of first increasing and then decreasing.
The analysis of the influence law of interaction term X 3 X 4 on the total power consumption evaluation indicator R 3 is shown in Figure 17c.When the vibration frequency X 3 is 5 Hz and the machine operation velocity X 4 increases from 0.6 km/h to 1.4 km/h, the total power consumption shows a trend of a gradual increase.This shows that when the vibration frequency is small, the power consumption can be reduced by reducing the machine operation velocity.When the vibration frequency is 9 Hz and the machine operation velocity increases from 0.6 km/h to 1.4 km/h, the total power consumption shows a trend of first decreasing and then increasing gently.
When the machine operation velocity is 0.6 km/h and the vibration frequency increases from 5 Hz to 9 Hz, the total power consumption shows a trend of a gradual decrease.When the machine operation velocity is 1.4 km/h and the vibration frequency increases from 5 Hz to 9 Hz, the total power consumption shows a trend of a gradual decrease.To sum up, when the machine operation velocity is small or large, the total power consumption can be reduced by increasing the vibration frequency.
3.6.Optimization and Test Verification of Target Parameters 3.6.1.Optimization of Target Parameters By using the optimization function provided in the Design Expert software and taking the value ranges of the screw velocity ratio, amplitude, vibration frequency, and machine operation velocity as constraint conditions, the objective function was set based on the following evaluation indicators: maximum sand removal rate, minimum crushing rate and minimum total power consumption.The optimum combination of key parameters in the conveying and separating device was thus determined.The set objective function and constraint function are as follows.
Objective The optimum combination of influence factor parameters is as follows: screw velocity ratio of 0.88, amplitude of 4.7 mm, vibration frequency of 7.5 Hz, and machine operation velocity of 0.92 km/h.Under these combined conditions, the sand removal rate is 90.40%, the crushing rate is 1.66%, and the power consumption is 2.24 kW in theory.

Test Verification
In order to verify the reliability of the optimization results, the performance indicators of conveying and sand removal were verified in a test, and the optimum combination of operating parameters of the device was adjusted for a soil tank test.The test was repeated three times, and the average value of the test results was taken.The sand removal rate obtained from the test was 88.92%, the crushing rate was 1.71%, and the total power consumption was 2.29 kW.The errors from the predicted values were 1.6%, 3.0%, and 2.2%, respectively.The test verification results show that the regression model has good reliability, and the conveying and separating device can meet the sand removal requirements in the tiger nut harvesting process.
On this basis, with reference to the national standard GB/T5262-2008 [47] Measuring Methods for Agricultural Machinery Testing Conditions-General Rules and DB34/T 534-2022 [48] Technical Specification for Peanut Harvesting Mechanization, a field performance test was carried out in the 54th Regiment, the Third Division of Xinjiang Production and Construction Corps, as shown in Figure 18.The results showed that the conveying and separating device worked normally, and the average total loss rate was 1.8% and the crushing rate was 4.9%, meeting the design requirements.

1.
A conveying and separating device combining screw conveying and a vibrating screen was designed for tiger nut harvesting.The conditions of axial and radial migration under the screw action and of the separation under vibratory action were analyzed.By simulating the separating process, the conveying and separating law of materials was clarified, and the rationality of the design was verified.

2.
A four-factor, five-level orthogonal central composite test was conducted.The test results were analyzed via the regression variance analysis method, and relation models between variable factors and evaluation indicators were constructed.The verification test results show that under the combined conditions of a screw velocity ratio of 0.88, an amplitude of 4.7 mm, a vibration frequency of 7.5 Hz, and a machine operation velocity of 0.92 km/h, the theoretical sand removal rate was 90.40%, the crushing rate was 1.66%, and the power consumption was 2.24 kW.

3.
The optimized results were verified via testing.The sand removal rate was 88.92%, the crushing rate was 1.71%, and the total power consumption was 2.29 kW.The errors from the predicted values were 1.6%, 3.0%, and 2.2%, respectively.A field performance test was carried out, and the conveying and separating device worked normally, meeting the requirements of conveying and separation for tiger nut harvesting.

Figure 1 .
Figure 1.The distribution of tiger nuts and its rhizomes.

Figure 1 .
Figure 1.The distribution of tiger nuts and its rhizomes.

Figure 3 .
Figure 3. Migration process of the complexes.

Figure 3 .
Figure 3. Migration process of the complexes.

Figure 4 .
Figure 4. Screw conveying device.1. Right-handed blade.2. Screen-cleaning spike teeth.3. Lefthanded blade.4. Screw shaft.The screw blade and screw shaft are the main components of the screw conveyor.The key structure and parameters of the screw conveyor were designed.The mass of sand in the excavation mixture was relatively large, so the sand density
37 • .According to the recommended values in the Transport Machinery Type Selection Design Manual [23], let D s = 400 m, S p = 400 mm, and D 0 = 108 mm.

Figure 5 .
Figure 5. Analysis of the velocity of the complex on the front screen surface.

Figure 5 .
Figure 5. Analysis of the velocity of the complex on the front screen surface.

Agriculture 2024 ,
14, x FOR PEER REVIEW 9 of 27 Considering the large number, length and load of screen bars, stiffeners were provided to position the screen bars and improve the stiffness of the screen bars.

Figure 9 .
Figure 9. Force analysis of leftward inertia force of sand particles.

Figure 9 .
Figure 9. Force analysis of leftward inertia force of sand particles.

Figure 10 .
Figure 10.Simulation model.(a) Discrete element model of soil tank.(b) The process of simulation.

Figure 10 .
Figure 10.Simulation model.(a) Discrete element model of soil tank.(b) The process of simulation.

Figure 11 .
Figure 11.Working area of the screw conveying device.

Figure 12 .Figure 11 .
Figure 12.Displacement of the sand particles in the action area of the screen-cleaning spike teeth.(a) X displacement.(b) Y displacement.(c) Z displacement.

Figure 11 .
Figure 11.Working area of the screw conveying device.

Figure 12 .Figure 12 .
Figure 12.Displacement of the sand particles in the action area of the screen-cleaning spike teeth.(a) X displacement.(b) Y displacement.(c) Z displacement.

Figure 11 .
Figure 11.Working area of the screw conveying device.

Figure 12 .Figure 13 .
Figure 12.Displacement of the sand particles in the action area of the screen-cleaning spike teeth.(a) X displacement.(b) Y displacement.(c) Z displacement.

Figure 14 .
Figure 14.The conveying and separating device for tiger nuts.

Figure 14 .
Figure 14.The conveying and separating device for tiger nuts.

3. 5 . 1 .
Analysis of the Influence Law of Significant Interaction Terms on the sand Removal Rate Evaluation Indicator, R 1

Figure 15 .
Figure 15.The influence of the interaction term on sand removal rate.(a) The influence law of interaction term X2X4 on the sand removal rate evaluation indicator, R1.(b) The influence law of interaction term X3X4 on the sand removal rate evaluation indicator, R1.

3. 5 . 2 .
Analysis of the Influence Law of Significant Interaction Terms on the Evaluation Indicator Crushing Rate R2

Figure 15 .
Figure 15.The influence of the interaction term on sand removal rate.(a) The influence law of interaction term X 2 X 4 on the sand removal rate evaluation indicator, R 1 .(b) The influence law of interaction term X 3 X 4 on the sand removal rate evaluation indicator, R 1 .

Figure 16 .
Figure 16.The influence of interaction terms on the crushing rate.(a) The influence of interact term X1X2 on sand removal rate evaluation indicator, R2.(b) The influence of interaction term X on the sand removal rate evaluation indicator, R2.3.5.3.Analysis of the Influence Law of Significant Interaction Terms on Total Powe Consumption Evaluation Indicator, R3

Figure 17 .Figure 16 .
Figure 17.The influence of interaction terms on the total power consumption.(a) The influenc interaction term X1 X2 on the total power consumption evaluation indicator, R3.(b) The influenc

Figure 17 .Figure 17 .
Figure 17.The influence of interaction terms on the total power consumption.(a) The influence of interaction term X1 X2 on the total power consumption evaluation indicator, R3.(b) The influence of Figure 17.The influence of interaction terms on the total power consumption.(a) The influence of interaction term X 1 X 2 on the total power consumption evaluation indicator, R 3 .(b) The influence of interaction term X 1 X 4 on the total power consumption evaluation indicator, R 3 .(c) The influence of interaction term X 3 X 4 on the total power consumption evaluation indicator, R 3 .

Agriculture 2024 , 27 Figure 18 .Figure 18 .
Figure 18.Test prototype.4. Conclusions 1.A conveying and separating device combining screw conveying and a vibrating screen was designed for tiger nut harvesting.The conditions of axial and radial migration under the screw action and of the separation under vibratory action were analyzed.By simulating the separating process, the conveying and separating law of materials was clarified, and the rationality of the design was verified.2. A four-factor, five-level orthogonal central composite test was conducted.The test Figure 18.Test prototype.

Table 1 .
Characteristic parameters of tiger nuts.
2.2.Overall Structure and Working Principle 2.2.1.Overall Structure of Conveying and Separating Device Test Bench

Table 1 .
Characteristic parameters of tiger nuts.
2.2.Overall Structure and Working Principle 2.2.1.Overall Structure of Conveying and Separating Device Test Bench

Table 2 .
Performance parameters of the conveying and separating device.

Table 2 .
Performance parameters of the conveying and separating device.

Table 5 .
The proportioning of the speed ratio.

Table 5 .
The proportioning of the speed ratio.

Table 6 .
Test factor levels and coding values.

Table 7 .
Test arrangement and results.

Table 8 .
Analysis of variance of measurement indexes.