Parameter optimization of winnowing equipment for machine- harvested Lycium barbarum L

To accurately and efficiently remove unripe fruit, flowers, leaves, and other impurities in machine-harvested Lycium barbarum L., winnowing equipment for machine-harvested L. barbarum based on the principle that different materials have different flight coefficients was designed. To optimize the structure and working parameters of winnowing equipment, this study adopted the free flow resistance model to establish a horizontal airflow model based on C++ in Microsoft Visual Studio. A discrete element method (DEM) simulation of ripe fruit in the horizontal airflow was performed using EDEM software. Results showed that the optimal parameters included an airflow speed of 5-6 m/s, input conveyor speed of 0.4-0.6 m/s, and input-output conveyor distance of 260-270 mm. We used three factors and three levels in a quadratic orthogonal rotation design to establish mathematical models regarding the rate of impurity change and the clearance rate of ripe fruit based on the airflow speed, input conveyor speed, and input-output conveyor distance. We also analyzed the effects of all factors on the rate of impurity change and the clearance rate of ripe fruit. The optimal parameter combination was an airflow speed of 5.52 m/s, input conveyor speed of 0.5 m/s, and input-output conveyor distance of 265.04 mm. The field experiment showed that the rate of impurity change and the clearance rate of ripe fruit were 89.74% and 8.71%, respectively. Findings provide a design basis for future research on winnowing equipment for machine-harvested L. barbarum. Additional keywords: flight coefficient; free flow resistance model; discrete element method; response surface methodology. Abbreviations used: DEM (discrete element method); DF (degree of freedom). Authors’ contributions: Conceived and designed the experiments: JZ, AS and JC. Performed the experiments: JZ, YC, GH, and EZ. Analyzed the data: JZ and FL. Wrote the paper: JZ, AS, and JC. All authors read and approved the final manuscript. Citation: Zhao, J.; Sugirbay, A.; Liu, F. Y.; Chen, Y.; Hu, G. G.; Zhang, E.; Chen, J. (2019). Parameter optimization of winnowing equipment for machine-harvested Lycium barbarum. Spanish Journal of Agricultural Research, Volume 17, Issue 2, e0203. https://doi. org/10.5424/sjar/2019172-14449 Received: 27 Dec 2018. Accepted: 01 Jul 2019. Copyright © 2019 INIA. This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC-by 4.0) License. Funding: National Key Research and Development Program of China (2017YFD0700402, 2018YFD0701102); Key Research and Development Program of the Ningxia Hui Autonomous Region (nxzdkjxm2016-04-02). Competing interests: The authors have declared that no competing interests exist. Correspondence should be addressed to Jun Chen: chenjun_jdxy@nwsuaf.edu.cn


Introduction
Lycium barbarum L. is a solanaceous Lycium shrub . Its ripe fruit contains healthpromoting bioactive components such as L. barbarum polysaccharides (Amagase & Farnsworth, 2011). The polysaccharides have demonstrated its effects on the improvement of renal function and the alleviation of inflammatory reaction (Zhao et al., 2016). For centuries, mechanized harvesting of L. barbarum has been challenging. Therefore, the fruit are harvested manually, which involves low efficiency and high cost Xu et al., 2018). With the continuous expansion of L. barbarum acreage, labor for fruit harvesting has become increasingly scarce, leading to more L. barbarum harvesting machines to be developed and used (So, 2001(So, , 2003Zhou & He, 2010;Zhang et al., 2015;He et al., 2017;Wang, 2018;Wang et al., 2018;Xu et al., 2018;Zhang et al., 2018a,b;Chen et al., 2019;Zhao et al., 2019b). However, many impurities such as unripe fruit, flowers, and leaves in machine-harvested L. barbarum seriously affect subsequent drying, storage, and other processing procedures Wang, 2018;Xu et al., 2018;Zhang et al., 2018b;Chen et al., 2019;Zhao et al., 2019b). Therefore, it is important to develop accurate and efficient winnowing equipment to remove impurities such as unripe fruit, flowers, and leaves in machine-harvested L. barbarum.
Winnowing equipment for machine-harvested L. barbarum was designed based on the principle that different materials have different flight coefficients. Using a discrete element method (DEM) simulation and the field experiment, the structure and working parameters of winnowing equipment were optimized. The results provide a design basis for future research and development of winnowing equipment for machineharvested L. barbarum.

Structure of winnowing equipment
Winnowing equipment for machine-harvested L. barbarum is shown in Fig. 1. The equipment consisted of a mechanical transmission system, winnowing system, and control system. The mechanical transmission system consisted of an input conveyor, output conveyor, electric motors, and drive systems; the winnowing system consisted of a fan, electronic speed regulator, and lithium battery; and the control system consisted of a controller, electronic speed governors, and relays.

Operating principle of winnowing equipment
As shown in Fig. 2, ripe fruit in the horizontal airflow was subjected to three forces: gravity G, horizontal airflow force P, and buoyancy F B . Because the buoyancy has a slight influence on the movement of ripe fruit in the horizontal airflow compared with the gravity and horizontal airflow force, it can be ignored (Liu, 2013). After simplification, the angle α between the resultant force of gravity and horizontal airflow force and the gravity can be obtained. When ripe fruit entered the horizontal airflow, its movement under the resultant force was parabolic and the distance of the movement of ripe fruit was related to the angle α (Liu, 2013). The larger the angle α, the farther the distance; the smaller the angle α, the closer the distance. The tangent of angle α can be expressed as the ratio of P to G and is calculated as follows (Liu, 2013).
where P is the horizontal airflow force, N; G is the gravity of ripe fruit, N; K is the resistance coeffi cient, which is related to the shape and surface characteristics of ripe fruit; ρ is the density of horizontal airflow, kg/m 3 ; B is the windward area of ripe fruit in the horizontal airflow, m 2 ; V a is the velocity of horizontal airflow, m/s; and V x is the initial velocity of ripe fruit entering the horizontal airflow, m/s.
The tangent of angle α is called the flight coeffi cient and is the characteristic value of the aerodynamic property of ripe fruit in the horizontal airflow (Liu, 2013). As shown in Table 1, because flying coefficients of different materials in the same horizontal airflow are different (Liu, 2013), the horizontal airflow can be used to remove impurities such as unripe fruit, flowers, and leaves in machineharvested L. barbarum.

Figure 1.
Winnowing equipment for machine-harvested Lycium barbarum. 1, frame; 2, electronic speed governor of the output conveyor; 3, relay of the output conveyor; 4, electric motor of the output conveyor; 5, drive system of the output conveyor; 6, output conveyor; 7, lithium battery; 8, electronic speed regulator; 9, fan; 10, input conveyor; 11, relay of the input conveyor; 12, electronic speed governor of the input conveyor; 13, drive system of the input conveyor; 14, electric motor of the input conveyor; 15, controller. [1] Parameter optimization of winnowing equipment for machine-harvested Lycium barbarum L. where W c is the calculated transmission power of the input conveyor, kW; and K A is the working condition coefficient. Therefore, the equations of the diameter of the driving wheel of the input conveyor D 1 , diameter of the driven wheel of the input conveyor D 2 , and center distance of the input conveyor s are as follows: where D 1 is the diameter of the driving wheel of the input conveyor, mm; m 1 is the revolving speed of the driving wheel of the input conveyor, r/min; D 2 is the diameter of the driven wheel of the input conveyor, mm; i is the transmission ratio of the input conveyor; ε is the sliding rate of the input conveyor, %; and s is the center distance of the input conveyor, mm.
The design of the output conveyor was as outlined above. After calculation, the design parameters of conveyors are listed in Table 2.
According to the nominal transmission power, the electric motors adopted alternating current electric motors (type: 5GN30KB; rated power: 60W; reduction ratio: 1:30; manufactured by Guangdong Shengkai Motor Manufacturing Co., Ltd., China). The input conveyor and output conveyor can be driven to rotate by the toothed chain. The toothed chain was designed according to the calculated transmission power as follows: where k z is the tooth number coefficient.

Winnowing system
A brushless motor (type: Sunnysky; KV value: 1400; manufactured by Zhongshan Langyu Model Co., Ltd., China) was used as the fan, powered by the lithium battery (type: Geshi ACE 5C; capacity: 3300 mAh;

Mechanical transmission system
Because machine-harvested L. barbarum was transported on the conveyors and we needed to adjust the input conveyor speed, the design of the input conveyor and output conveyor was highly important. The two transmission shafts were parallel, the rotation direction of the two transmission shafts was the same, and the input conveyor speed was ≤ 20 m/s; thus, the open drive mode can be selected (Qiu, 2011). As the transmission power W was<500 kW, the belt of the input conveyor and output conveyor was made of antistatic flat belt (Qiu, 2011). The design equations [2][3][4][5][6][7][8] were referred to this book (Qiu, 2011). The effective tensile force was obtained as follows: where F is the effective tensile force of the input conveyor, N; F 1 is the tight-side tensile force of the input conveyor, N; F 2 is the slack-side tensile force of the input conveyor, N; W is the nominal transmission power of the input conveyor, kW; and X 2 is the input conveyor speed, m/s.
The relation equation between the tight-side tensile force of the input conveyor and the slack-side tensile force of the input conveyor can be obtained as follows: where q is the mass of a meter length belt of the input conveyor, kg/m; e is the Napierian base; μ is the friction coefficient between the input conveyor and its belt wheel; and λ is the cornerite between the input conveyor and its belt wheel, (°).
Because the input conveyor speed was<10 m/s, the centrifugal force can be ignored (Qiu, 2011). Because λ=π, μ ≈ 0.3. Therefore, the calculated transmission power W c can be obtained as follows:  (2015). [2] [4] [5] [6] [7] manufactured by Shenzhen Grepow Battery Co., Ltd., China). The controller generated pulse width modulation waves to control the run and stop of the fan by controlling the electronic speed regulator (type: Skywalker 40A; manufactured by Shenzhen Hobbywing Technology Co., Ltd., China) and to adjust the airflow speed.

Control system
The circuit diagram of winnowing equipment is shown in Fig. 3. The type of central processing unit of controller (type: STM32-F103; manufactured by Guangzhou Hard Rock Technology Co., Ltd., China) was STM32-F103 ZET6, and its development environment was Keil uVision5. When machine-harvested L. bar barum entered the input conveyor, the fan started. The controller controlled the run and stop of electric motors by controlling the relays (type: SRD-05VDC-SL-C; manufactured by Ningbo Songle Relay Co., Ltd., China) and also adjusted the revolving speed of electric motors by controlling the electronic speed governors, so as to control the run and stop of the two conveyors and to adjust the input conveyor speed. The two conveyors rotated continuously to remove impurities such as unripe fruit, flowers, and leaves in machineharvested L. barbarum.

DEM simulation of ripe fruit in the horizontal airflow
The horizontal airflow model was established based on C++ in Microsoft Visual Studio. Because the size of ripe fruit was approximately consistent and the horizontal airflow parameters (e.g., velocity, density, and viscosity) were basically constant, the free flow resistance model was adopted (Hu, 2010). To accurately simulate the movement of ripe fruit in the horizontal airflow, its buoyancy and free flow resistance force should be considered (Liu et al., 2015;Liu, 2018;Yu et al., 2018). The equations [9-12] were referred to Hu (2010). The equations of the buoyancy F B and the free flow resistance force F d are as follows: where F B is the buoyancy, N; F d is the free flow resistance force, N; g is the gravitational acceleration, m/s 2 ; V is the volume of ripe fruit, m 3 ; A is the projected area of ripe fruit perpendicular to the horizontal airflow Table 2. Design parameters of conveyors.

Width of the belt of the input conveyor 200
Width of the belt of the output conveyor 280 Diameter of the driving wheel of the input conveyor

160
Diameter of the driving wheel of the output conveyor 224 Figure 3. Circuit diagram of winnowing equipment. [9] [10] Parameter optimization of winnowing equipment for machine-harvested Lycium barbarum L.
Spanish Journal of Agricultural Research June 2019 • Volume 17 • Issue 2 • e0203 5 direction, m 2 ; v is the relative velocity between the ripe fruit and horizontal airflow, m/s; and C D is the resistance coefficient which is calculated as follows: Therefore, C D depends on the Reynolds number Re, and its equation is as follows: where Re is the Reynolds number; β is the free volume of a grid cell; L is the diameter of ripe fruit, m; and η is the viscosity of horizontal airflow, kg/ (m·s).
The movement of ripe fruit in the horizontal airflow was simulated based on DEM using EDEM software. In this study, the Moving Plane model was selected; the parameters of materials are listed in Table 1.
Based on previous research and references review (González-Montellano et al., 2012;Horabik & Molenda, 2016;Castillo-Ruiz et al., 2018;Chen et al., 2019;Zhao et al., 2019a,b), because the properties of L. barbarum are similar to the olives, some simulation parameters of olive can be cited based on the actual situation. Based on González-Montellano et al.´s (2012) work, because ripe fruit were distributed dispersedly with less collision compared with the olives, the restitution coefficient between ripe fruit and ripe fruit and the same coefficient between ripe fruit and the input conveyor were 0.1 and 0.1, respectively. Based on the reference (González-Montellano et al., 2012), because there are irregular bulges on partial surfaces of ripe fruit and the mechanics parameters of the two conveyors material were listed in Table 1, the static friction coefficient between ripe fruit and the input conveyor was 0.5. Besides, the friction force between ripe fruit and ripe fruit was less than the friction force between ripe fruit and the input conveyor, the static friction coefficient between ripe fruit and ripe fruit was 0.3. Based on references review (Li et al., 2012;Liu, 2018), the rolling friction coefficient between ripe fruit and ripe fruit and the same coefficient between ripe fruit and the input conveyor were select as the default value (i.e., 0.01). Because the mass of ripe fruit is not very light and ripe fruit were distributed dispersedly, the movement of ripe fruit in the horizontal airflow was mainly affected by the gravity, horizontal airflow force, buoyancy, and free flow resistance force. Based on the pre-simulation, we came to a conclusion that the six coefficients have little effect on the simulation results. Therefore, the six coefficients were appropriate for this simulation. Furthermore, the ripe fruit model used a sphere with a radius of 4.5 mm. A 4-edged polygon served as the particle factory surface, and the type of particle factory was unlimited quantity. The DEM simulation of ripe fruit in the horizontal airflow is shown in Fig. 4.
According to the simulation, the airflow speed, input conveyor speed, and input-output conveyor distance greatly affected the rate of impurity change and the clearance rate of ripe fruit. Simulation results showed that winnowing equipment met these conditions, with most of ripe fruit falling on the output conveyor and most of unripe fruit, flowers, and leaves removed at an airflow speed of 5-6 m/s, input conveyor speed of 0.4-0.6 m/s, and input-output conveyor distance of 260-270 mm. Fig. 5 shows the trajectory of the ripe fruit which was the most easily removed by mistake.

Performance experiment of winnowing equipment
Ningqi 7 was selected as the experiment variety. The following instruments were used: 1) an electronic vernier caliper (type: AIRAJ second-generation product; range: 0-300 mm; precision: 0.01 mm; manufactured by Qingdao Yigou Hardware Tools Co., Ltd., China); 2) a digital illuminometer (type: PM6612L; manufactured by Shenzhen New Huayi Instrument Co., Ltd., China); 3) a tachometer (type: DT2236B; manufactured by Shenzhen Sanpo Instrument Co., Ltd., China); 4) a digital anemometer (type: PM6252B; manufactured by Shenzhen New Huayi Instrument Co., Ltd., China).  The winnowing equipment was used to remove as many impurities as possible (e.g., unripe fruit, flowers, and leaves) in machine-harvested L. barbarum while removing as little ripe fruit as possible. The purpose of the experiment was to investigate the performance of the winnowing equipment; therefore, the rate of impurity change Y 1 and the clearance rate of ripe fruit Y 2 were selected as performance indices in the experiment. Corresponding equations to calculate these indices are as follows: [13] [14] where n 1 is the amount of unripe fruit, flowers, and leaves before winnowing; n 2 is the amount of ripe fruit, unripe fruit, flowers, and leaves before winnowing; n 3 is the amount of unripe fruit, flowers, and leaves after winnowing; n 4 is the amount of ripe fruit, unripe fruit, flowers, and leaves after winnowing; n 5 is the amount of removed ripe fruit after winnowing; and n 6 is the amount of ripe fruit before winnowing.
The experiment was conducted in Zhongning in the Ningxia Hui Autonomous Region (37°22′56″N, 105°37′21″E) on June 26, 2018; the temperature was 27.2°C, the humidity was 30.1%, and the illuminance was 520.6 Lux. Based on the theoretical analysis, the main factors affecting the performance of the winnowing equipment were determined to be the airflow speed, input conveyor speed, and inputoutput conveyor distance. The previous tests and simulation analysis indicated the following suitable value range for the factors: the airflow speed X 1 was 5-6 m/s; the input conveyor speed X 2 was 0.4-0.6 m/s; and the input-output conveyor distance X 3 was 260-270 mm. The airflow speed was adjusted by the electronic speed regulator and measured with the digital anemometer. The input conveyor speed was adjusted by the electronic speed governor and measured with the tachometer. The input-output conveyor distance was adjusted by lifting the input conveyor and measured with the electronic vernier caliper.
Three factors and three levels were used in a quadratic orthogonal rotation design; the codes of factors are listed in Table 3, and the experimental schemes and results are presented in Table 4. Seventeen groups of tests were carried out in this experiment, and each group was repeated 5 times. The mean value of 5 results was used as the group result. The machine-harvested L. barbarum from a shrub using a handheld L. barbarum harvester was used as a group sample. The amounts of ripe fruit, unripe fruit, flowers and leaves were manually counted before a test and the amounts of that were also manually counted after a test. Experiment design and analysis were performed in Design-Expert (Yuan et al., 2012;Wang et al., 2013).

Regression analysis
The regression model of the rate of impurity change (response) using codes of all factors as variables was as follows: ANOVA results of the rate of impurity change are shown in Table 5. The model was statistically significant (p<0.0001), and X 1 , X 2 , X 3 , X 1 X 2 , X 1 2 , X 2 2 , and X 3 2 each had a significant effect on the rate of impurity change (p<0.05). Moreover, the lack-of-fit test indicated that the model had a good fit (p=0.1569).
The regression model of the clearance rate of ripe fruit (response) using codes of all factors as variables was as follows:  Y 1 = 90.04 + 0.61X 1 + 0.27X 2 + 0.24X 3 -0.095X 1 X 2 + + 0.04X 1 X 3 + 0.025X 2 X 3 -1.24X 1 2 -0.7X 2 2 -0.68X 3 ANOVA results of the clearance rate of ripe fruit are shown in Table 6. The model was statistically significant (p<0.0001), and X 1 , X 2 , X 3 , X 1 X 2 , X 1 2 , X 2 2 , and X 3 2 each had a significant effect on the clearance rate of ripe fruit (p<0.05). Moreover, the lack-of-fit test indicated that the model had a good fit (p=0.0667).

Response surface analysis
The response surface methodology was used to analyze the effects of all factors on the rate of impurity change; the response surface results of the regression equation (Eq. [15]) are shown in Fig. 6. As shown in Table 5, the airflow speed, input conveyor speed, and input-output conveyor distance each had an extremely significant effect on the rate of impurity change. The  [16] Y 2 = 8.51 + 0.5X 1 + 0.27X 2 + 0.23X 3 + 0.082X 1 X 2 + + 0.01X 1 X 3 + 0.012X 2 X 3 + 1.02X 1 2 + 1.28X 2 2 + 1.24X 3 interaction effect between the airflow speed and input conveyor speed was significant. As shown in Fig. 6a, as the airflow speed increased, the rate of impurity change first increased rapidly and then decreased slowly. As shown in Fig. 6b, as the input-output conveyor distance increased, the rate of impurity change first increased slowly and then decreased slowly. As shown in Fig.  6c, as the input conveyor speed increased, the rate of impurity change first increased slowly and then decreased slowly. The response surface methodology was used to analyze the effects of all factors on the clearance rate of ripe fruit; the response surface results of the regression equation (Eq. [16]) are shown in Fig. 7. As shown in Table 6, the airflow speed, input conveyor speed, and input-output conveyor distance each had an extremely significant effect on the clearance rate of ripe fruit. The interaction effect between the airflow speed and input conveyor speed was significant. As shown in Fig. 7a, as the airflow speed increased, the clearance rate of ripe fruit first decreased lowly and then increased rapidly; in Fig. 7b, as the input-output conveyor distance increased, the clearance rate of ripe fruit first decreased rapidly and then increased rapidly; and in Fig. 7c, as the input conveyor speed increased, the clearance rate of ripe fruit first decreased rapidly and then increased rapidly.
After prediction, the optimal parameter combination was determined to be an airflow speed of 5.52 m/s, input conveyor speed of 0.5 m/s, and input-output conveyor distance of 265.04 mm.

Field experiment verification
The field experiment was completed on June 28, 2018 and 15 groups were conducted in this experiment to eliminate random errors. The chosen experiment parameter combination was an airflow speed of 5.52 m/s, input conveyor speed of 0.5 m/s, and inputoutput conveyor distance of 265.04 mm. As depicted in  . Response surface and contour plots for effects of all factors on the rate of impurity change: X 1 and X 2 (a), X 1 and X 3 (b), and X 2 and X 3 (c). Fig. 8, a high-speed camera (type: OLYMPUS i-speed TR; manufactured by Keymed (Medical & Industrial Equipment) Co., Ltd., UK) was used to record the winnowing process; the duration was 200 ms, the interval was 40 ms, and we used 1000 f/s (Torregrosa et al., 2014;Liang et al., 2018). Red, green, purple, and blue points were used to indicate the respective positions of the ripe fruit, unripe fruit, flower, and leaf every 40 ms. The field experiment showed that the rate of impurity change and the clearance rate of ripe fruit were 89.74% and 8.71%, respectively.

Discussion
The machine-harvested L. barbarum contained many impurities such as unripe fruit, flowers, and leaves Wang, 2018;Xu et al., 2018;Zhang et al., 2018b;Chen et al., 2019;Zhao et al., 2019b). It seriously affected subsequent drying, storage, and other processing procedures. Therefore, it is important to develop winnowing equipment which can accurately and efficiently remove unripe fruit, flowers, leaves, and other impurities in machine-harvested L. barbarum. In the present study, winnowing equipment for machineharvested L. barbarum based on the principle that different materials have different flight coefficients was designed. The structure and working parameters of winnowing equipment were optimized based on a DEM simulation and the field experiment. Moreover, we established mathematical models regarding the rate of impurity change and the clearance rate of ripe fruit based on the airflow speed, input conveyor speed, and input-output conveyor distance and analyzed the effects of all factors on the rate of impurity change and the clearance rate of ripe fruit. We came to a conclusion that the optimal parameter combination was an airflow speed of 5.52 m/s, input conveyor speed of 0.5 m/s, and inputoutput conveyor distance of 265.04 mm. The field experiment showed that the rate of impurity change Figure 7. Response surface and contour plots for effects of all factors on the clearance rate of ripe fruit: X 1 and X 2 (a), X 1 and X 3 (b), and X 2 and X 3 (c). and the clearance rate of ripe fruit were 89.74% and 8.71%, respectively. These results provide a design basis for future research on winnowing equipment for machine-harvested L. barbarum. But there are still some problems to be solved. For example, it is important to make miniaturized winnowing equipment, so as to real-time winnowing during harvesting. Furthermore, a more accurate and efficient winnowing method is also worth further studying.