<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" dtd-version="3.0" xml:lang="en">
   <front>
      <journal-meta>
         <journal-id journal-id-type="publisher-id">SJAR</journal-id>
         <journal-title-group>
            <journal-title>Spanish Journal of Agricultural Research</journal-title>
            <abbrev-journal-title>SJAR</abbrev-journal-title>
         </journal-title-group>
         <issn pub-type="epub">2171-9292</issn>
         <publisher>
            <publisher-name>Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA)</publisher-name>
         </publisher>
      </journal-meta>
      <article-meta>
         <article-id pub-id-type="publisher-id">14449</article-id>
         <article-id pub-id-type="doi">10.5424/sjar/2019172-14449</article-id>
         <article-categories>
            <subj-group subj-group-type="heading">
               <subject>Research article</subject>
            </subj-group>
         </article-categories>
         <title-group>
            <article-title>
               Parameter optimization of winnowing equipment for machine-harvested
               <italic>Lycium barbarum</italic>
               L.
            </article-title>
         </title-group>
         <contrib-group>
            <contrib contrib-type="author" corresp="no">
               <name>
                  <surname>Zhao</surname>
                  <given-names>Jian</given-names>
                  <aff>
                     <i>Northwest A&amp;F Univ., College of Mechan. &amp; Electron. Eng., Yangling, 712100 Shaanxi, China.</i>
                  </aff>
               </name>
            </contrib>
            <contrib contrib-type="author" corresp="no">
               <name>
                  <surname>Sugirbay</surname>
                  <given-names>Adilet</given-names>
                  <aff>
                     <i>Northwest A&amp;F Univ., College of Mechan. &amp; Electron. Eng., Yangling, 712100 Shaanxi, China.</i>
                     <i>S. Seifullin Kazakh Agro Tech. Univ., Techn. Fac., Astana 010000, Kazakhstan.</i>
                  </aff>
               </name>
            </contrib>
            <contrib contrib-type="author" corresp="no">
               <name>
                  <surname>Liu</surname>
                  <given-names>Fanyi</given-names>
                  <aff>
                     <i>SouthWest Univ., College of Eng. &amp; Technol., Beibei, 400715 Chongqing, China.</i>
                  </aff>
               </name>
            </contrib>
            <contrib contrib-type="author" corresp="no">
               <name>
                  <surname>Chen</surname>
                  <given-names>Yun</given-names>
                  <aff>
                     <i>Northwest A&amp;F Univ., College of Mechan. &amp; Electron. Eng., Yangling, 712100 Shaanxi, China.</i>
                  </aff>
               </name>
            </contrib>
            <contrib contrib-type="author" corresp="no">
               <name>
                  <surname>Hu</surname>
                  <given-names>Guangrui</given-names>
                  <aff>
                     <i>Northwest A&amp;F Univ., College of Mechan. &amp; Electron. Eng., Yangling, 712100 Shaanxi, China.</i>
                  </aff>
               </name>
            </contrib>
            <contrib contrib-type="author" corresp="no">
               <name>
                  <surname>Zhang</surname>
                  <given-names>Enyu</given-names>
                  <aff>
                     <i>Northwest A&amp;F Univ., College of Mechan. &amp; Electron. Eng., Yangling, 712100 Shaanxi, China.</i>
                  </aff>
               </name>
            </contrib>
            <contrib contrib-type="author" corresp="yes">
               <name>
                  <surname>Chen</surname>
                  <given-names>Jun</given-names>
                  <aff>
                     <i>Northwest A&amp;F Univ., College of Mechan. &amp; Electron. Eng., Yangling, 712100 Shaanxi, China.</i>
                  </aff>
               </name>
            </contrib>
         </contrib-group>
         <author-notes>
            <corresp>
               should be addressed to Jun Chen:
               <email xlink:href="chenjun_jdxy@nwsuaf.edu.cn">chenjun_jdxy@nwsuaf.edu.cn</email>
            </corresp>
         </author-notes>
         <pub-date pub-type="epub">
            <day>01</day>
            <month>06</month>
            <year>2019</year>
         </pub-date>
         <pub-date pub-type="collection">
            <year>2019</year>
         </pub-date>
         <volume>17</volume>
         <issue>2</issue>
         <elocation-id content-type="doi">10.5424/sjar/2019172-14449</elocation-id>
         <history>
            <date date-type="recibido">
               <day>27</day>
               <month>12</month>
               <year>2018</year>
            </date>
            <date date-type="aceptado">
               <day>01</day>
               <month>07</month>
               <year>2019</year>
            </date>
         </history>
         <permissions>
            <copyright-statement>© 2019 INIA</copyright-statement>
            <copyright-year>2019</copyright-year>
            <license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by-nc/3.0/">
               <license-p>This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC-by 4.0) License.</license-p>
            </license>
         </permissions>
         <abstract id="abstract01">
            <title>Abstract</title>
            <p>
               To accurately and efficiently remove unripe fruit, flowers, leaves, and other impurities in machine-harvested
               <italic>Lycium barbarum</italic>
               L., winnowing equipment for machine-harvested
               <italic>L. barbarum</italic>
               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
               <italic>L. barbarum</italic>
               .
            </p>
         </abstract>
         <kwd-group>
            <title>Additional key words:</title>
            <kwd>flight coefficient;</kwd>
            <kwd>free flow resistance model;</kwd>
            <kwd>discrete element method;</kwd>
            <kwd>response surface methodology.</kwd>
         </kwd-group>
         <kwd-group>
            <title>Additional key words:</title>
            <kwd>DEM (discrete element method);</kwd>
            <kwd>DF (degree of freedom).</kwd>
         </kwd-group>
         <funding-group>
            <funding-statement>National Key Research and Development Program of China (2017YFD0700402, 2018YFD0701102); Key Research and Development Program of the Ningxia Hui Autonomous Region (nxzdkjxm2016-04-02).</funding-statement>
         </funding-group>
      </article-meta>
      <notes>
         <p>
            <bold>Author's contributions:</bold>
            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.
         </p>
         <p>
            <bold>Citation</bold>
            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
            <italic>Lycium barbarum</italic>
            . Spanish Journal of Agricultural Research, Volume 17, Issue 2, e0203.
            <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.5424/sjar/2019172-14449">https://doi.org/10.5424/sjar/2019172-14449</ext-link>
         </p>
         <p>
            <bold>Competing interests:</bold>
            The authors have declared that no competing interests exist.
         </p>
      </notes>
   </front>
   <body>
      <sec id="S1">
         <title>Introduction</title>
         <p>
            <italic>Lycium barbarum</italic>
            L. is a solanaceous
            <italic>Lycium</italic>
            shrub (
            <xref ref-type="bibr" rid="b21">
               Xu
               <italic>et al.</italic>
               , 2018
            </xref>
            ). Its ripe fruit contains health-promoting bioactive components such as
            <italic>L. barbarum</italic>
            polysaccharides (
            <xref ref-type="bibr" rid="b1">Amagase &amp; Farnsworth, 2011</xref>
            ). The polysaccharides have demonstrated its effects on the improvement of renal function and the alleviation of inflammatory reaction (
            <xref ref-type="bibr" rid="b29">
               Zhao
               <italic>et al.</italic>
               , 2016
            </xref>
            ). For centuries, mechanized harvesting of
            <italic>L. barbarum</italic>
            has been challenging. Therefore, the fruit are harvested manually, which involves low efficiency and high cost (
            <xref ref-type="bibr" rid="b26">
               Zhang
               <italic>et al</italic>
               ., 2015
            </xref>
            ;
            <xref ref-type="bibr" rid="b21">
               Xu
               <italic>et al</italic>
               ., 2018
            </xref>
            ). With the continuous expansion of
            <italic>L. barbarum</italic>
            acreage, labor for fruit harvesting has become increasingly scarce, leading to more
            <italic>L. barbarum</italic>
            harvesting machines to be developed and used (
            <xref ref-type="bibr" rid="b15">So, 2001</xref>
            ,
            <xref ref-type="bibr" rid="b16">2003</xref>
            ;
            <xref ref-type="bibr" rid="b30">Zhou &amp; He, 2010</xref>
            ;
            <xref ref-type="bibr" rid="b26">
               Zhang
               <italic>et al.</italic>
               , 2015
            </xref>
            ;
            <xref ref-type="bibr" rid="b6">
               He
               <italic>et al.</italic>
               , 2017
            </xref>
            ;
            <xref ref-type="bibr" rid="b19">Wang, 2018</xref>
            ;
            <xref ref-type="bibr" rid="b20">
               Wang
               <italic>et al.</italic>
               , 2018
            </xref>
            ;
            <xref ref-type="bibr" rid="b21">
               Xu
               <italic>et al.</italic>
               , 2018
            </xref>
            ;
            <xref ref-type="bibr" rid="b24">
               Zhang
               <italic>et al.</italic>
               , 2018a
            </xref>
            ,
            <xref ref-type="bibr" rid="b25">b</xref>
            ;
            <xref ref-type="bibr" rid="b4">
               Chen
               <italic>et al.</italic>
               , 2019
            </xref>
            ;
            <xref ref-type="bibr" rid="b28">
               Zhao
               <italic>et al.</italic>
               , 2019b
            </xref>
            ). However, many impurities such as unripe fruit, flowers, and leaves in machine-harvested
            <italic>L. barbarum</italic>
            seriously affect subsequent drying, sto­rage, and other processing procedures (
            <xref ref-type="bibr" rid="b26">
               Zhang
               <italic>et al.</italic>
               , 2015
            </xref>
            ;
            <xref ref-type="bibr" rid="b19">Wang, 2018</xref>
            ;
            <xref ref-type="bibr" rid="b21">
               Xu
               <italic>et al.</italic>
               , 2018
            </xref>
            ;
            <xref ref-type="bibr" rid="b25">
               Zhang
               <italic>et al.</italic>
               , 2018b
            </xref>
            ;
            <xref ref-type="bibr" rid="b4">
               Chen
               <italic>et al.</italic>
               , 2019
            </xref>
            ;
            <xref ref-type="bibr" rid="b28">
               Zhao
               <italic>et al.</italic>
               , 2019b
            </xref>
            ). Therefore, it is important to develop accurate and efficient winnowing equipment to remove impurities such as unripe fruit, flowers, and leaves in machine-harvested
            <italic>L. barbarum</italic>
            .
         </p>
         <p>
            Winnowing equipment for machine-harvested
            <italic>L. barbarum</italic>
            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 machine-harvested
            <italic>L. barbarum</italic>
            .
         </p>
      </sec>
      <sec id="S2">
         <title>Material and methods</title>
         <sec id="S2.1">
            <title>Structure and operating principle of winnowing equipment</title>
            <p>
               <italic>Structure of winnowing equipment</italic>
            </p>
            <p>
               Winnowing equipment for machine-harvested
               <italic>L. barbarum</italic>
               is shown in <xref ref-type="fig" rid="F1">Fig. 1</xref>. The equipment consisted of a mechanical transmission system, winnowing sys­tem, 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.
            </p>
            <fig id="F1">
    <label>Figure 1.</label>
    <caption>
    <title>Winnowing equipment for machine-harvested
<italic>Lycium barbarum</italic>. 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.</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f01.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

            <p />
            <p>
               <italic>Operating principle of winnowing equipment</italic>
            </p>
            <p>
               As shown in <xref ref-type="fig" rid="F2">Fig. 2</xref>, ripe fruit in the horizontal airflow was subjected to three forces: gravity
               <italic>G</italic>
               , horizontal airflow force
               <italic>P</italic>
               , and buoyancy
               <italic>F</italic>
               <sub>B</sub>
               . 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 (
               <xref ref-type="bibr" rid="b13">Liu, 2013</xref>
               ). After simplification, the angle
               <italic>&#945;</italic>
               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
               <italic>&#945;</italic>
               (
               <xref ref-type="bibr" rid="b13">Liu, 2013</xref>
               ). The larger the angle
               <italic>&#945;</italic>
               , the farther the distance; the smaller the angle
               <italic>&#945;</italic>
               , the closer the distance. The tangent of angle
               <italic>&#945;</italic>
               can be expressed as the ratio of
               <italic>P</italic>
               to
               <italic>G</italic>
               and is calculated as follows (
               <xref ref-type="bibr" rid="b13">Liu, 2013</xref>
               ).
            </p>
            <graphic id="form1" xlink:href="sjar_e0203_form1.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>P</italic>
               is the horizontal airflow force, N;
               <italic>G</italic>
               is the gravity of ripe fruit, N;
               <italic>K</italic>
               is the resistance coef­ficient, which is related to the shape and surface characteristics of ripe fruit;
               <italic>&#961;</italic>
               is the density of horizontal airflow, kg/m
               <sup>3</sup>
               ;
               <italic>B</italic>
               is the windward area of ripe fruit in the horizontal airflow, m
               <sup>2</sup>
               ;
               <italic>V</italic>
               <sub>a</sub>
               is the velocity of horizontal airflow, m/s; and
               <italic>V</italic>
               <sub>x</sub>
               is the initial velocity of ripe fruit entering the horizontal airflow, m/s.
            </p>
            <fig id="F2">
    <label>Figure 2.</label>
    <caption>
    <title>Force analysis of ripe fruit in the horizontal
airflow. <italic>R </italic>is the resultant force of gravity and horizontal
airflow force, N.</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f02.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

            <p>
               The tangent of angle
               <italic>&#945;</italic>
               is called the flight coeffi­cient and is the characteristic value of the aerodynamic property of ripe fruit in the horizontal airflow (
               <xref ref-type="bibr" rid="b13">Liu, 2013</xref>
               ). As shown in <xref ref-type="table" rid="T1">Table 1</xref>, because flying coefficients of different materials in the same horizontal airflow are different (
               <xref ref-type="bibr" rid="b13">Liu, 2013</xref>
               ), the horizontal airflow can be used to remove impurities such as unripe fruit, flowers, and leaves in machine-harvested
               <italic>L. barbarum</italic>
               .
            </p>
            <table-wrap id="T1">
    <label>Table 1.</label>
    <caption>
    <title>Shape and mechanics parameters of materials. </title>
    </caption>
    <graphic xlink:href="sjar_e0203_t01.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</table-wrap>

         </sec>
         <sec id="S2.2">
            <title>Design of key systems</title>
            <p />
            <p>
               <italic>Mechanical transmission system</italic>
            </p>
            <p>
               Because machine-harvested
               <italic>L. barbarum</italic>
               was trans­­ported 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 &#8804; 20 m/s; thus, the open drive mode can be selected (
               <xref ref-type="bibr" rid="b14">Qiu, 2011</xref>
               ). As the transmission power
               <italic>W</italic>
               was&lt;500 kW, the belt of the input conveyor and output conveyor was made of anti-static flat belt (
               <xref ref-type="bibr" rid="b14">Qiu, 2011</xref>
               ). The design <xref ref-type="disp-formula" rid="form2">equations [2</xref>-<xref ref-type="disp-formula" rid="form8">8</xref>] were referred to this book (
               <xref ref-type="bibr" rid="b14">Qiu, 2011</xref>
               ). The effective tensile force was obtained as follows:
            </p>
           <graphic id="form2" xlink:href="sjar_e0203_form2.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>F</italic>
               is the effective tensile force of the input con­veyor, N;
               <italic>F</italic>
               <sub>1</sub>
               is the tight-side tensile force of the input conveyor, N;
               <italic>F</italic>
               <sub>2</sub>
               is the slack-side tensile force of the input conveyor, N;
               <italic>W</italic>
               is the nominal transmission power of the input conveyor, kW; and
               <italic>X</italic>
               <sub>2</sub>
               is the input conveyor speed, m/s.
            </p>
            <p>The relation equation between the tight-side ten­sile force of the input conveyor and the slack-side tensile force of the input conveyor can be obtained as follows:</p>
           <graphic id="form3" xlink:href="sjar_e0203_form3.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>q</italic>
               is the mass of a meter length belt of the input conveyor, kg/m;
               <italic>e</italic>
               is the Napierian base;
               <italic>&#181;</italic>
               is the friction coefficient between the input conveyor and its belt wheel; and
               <italic>&#955;</italic>
               is the cornerite between the input conveyor and its belt wheel, (&#176;).
            </p>
            <p>
               Because the input conveyor speed was&lt;10 m/s, the centrifugal force can be ignored (
               <xref ref-type="bibr" rid="b14">Qiu, 2011</xref>
               ). Because
               <italic>&#955;</italic>
               =&#960;,
               <italic>&#181;</italic>
               &#8776; 0.3. Therefore, the calculated transmission power
               <italic>W</italic>
               <sub>c</sub>
               can be obtained as follows:
            </p>
          <graphic id="form4" xlink:href="sjar_e0203_form4.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>W</italic>
               <sub>c</sub>
               is the calculated transmission power of the input conveyor, kW; and
               <italic>K</italic>
               <sub>A</sub>
               is the working condition coefficient.
            </p>
            <p>
               Therefore, the equations of the diameter of the driving wheel of the input conveyor
               <italic>D</italic>
               <sub>1</sub>
               , diameter of the driven wheel of the input conveyor
               <italic>D</italic>
               <sub>2</sub>
               , and center distance of the input conveyor
               <italic>s</italic>
               are as follows:
            </p>
           <graphic id="form5" xlink:href="sjar_e0203_form5.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
           <graphic id="form6" xlink:href="sjar_e0203_form6.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
           <graphic id="form7" xlink:href="sjar_e0203_form7.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>D</italic>
               <sub>1</sub>
               is the diameter of the driving wheel of the input conveyor, mm;
               <italic>m</italic>
               <sub>1</sub>
               is the revolving speed of the driving wheel of the input conveyor, r/min;
               <italic>D</italic>
               <sub>2</sub>
               is the diameter of the driven wheel of the input conveyor, mm;
               <italic>i</italic>
               is the transmission ratio of the input conveyor;
               <italic>ε</italic>
               is the sliding rate of the input conveyor, %; and
               <italic>s</italic>
               is the center distance of the input conveyor, mm.
            </p>
            <p>The design of the output conveyor was as outlined above. After calculation, the design parameters of conveyors are listed in <xref ref-type="table" rid="T2">Table 2</xref>.</p>
            <table-wrap id="T2">
    <label>Table 2.</label>
    <caption>
    <title>Design parameters of conveyors. </title>
    </caption>
    <graphic xlink:href="sjar_e0203_t02.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</table-wrap>

            <p>According to the nominal transmission power, the electric motors adopted alternating current electric mo­tors (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:</p>
           <graphic id="form8" xlink:href="sjar_e0203_form8.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>k</italic>
               <sub>z</sub>
               is the tooth number coefficient.
            </p>
            <p />
            <p>
               <italic>Winnowing system</italic>
            </p>
            <p>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; 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.</p>
            <p />
            <p>
               <italic>Control system</italic>
            </p>
            <p>
               The circuit diagram of winnowing equipment is shown in <xref ref-type="fig" rid="F3">Fig. 3</xref>. 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
               <italic>L. bar­barum</italic>
               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 machine-harvested
               <italic>L. barbarum</italic>
               .
            </p>
            <fig id="F3">
    <label>Figure 3.</label>
    <caption>
    <title>Circuit diagram of winnowing equipment.</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f03.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

         </sec>
         <sec id="S2.3">
            <title>DEM simulation of ripe fruit in the horizontal airflow</title>
            <p>
               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 (
               <italic>e.g.</italic>
               , velocity, density, and viscosity) were basically constant, the free flow resistance model was adopted (
               <xref ref-type="bibr" rid="b8">Hu, 2010</xref>
               ). To accurately simulate the movement of ripe fruit in the horizontal airflow, its buoyancy and free flow resistance force should be considered (
               <xref ref-type="bibr" rid="b12">
                  Liu
                  <italic>et al.</italic>
                  , 2015
               </xref>
               ;
               <xref ref-type="bibr" rid="b11">Liu, 2018</xref>
               ;
               <xref ref-type="bibr" rid="b22">
                  Yu
                  <italic>et al.</italic>
                  , 2018
               </xref>
               ). The <xref ref-type="disp-formula" rid="form9">equations [9</xref>-<xref ref-type="disp-formula" rid="form12">12</xref>] were referred to Hu (2010). The equations of the buoyancy
               <italic>F</italic>
               <sub>B</sub>
               and the free flow resistance force
               <italic>F</italic>
               <sub>d</sub>
               are as follows:
            </p>
          <graphic id="form9" xlink:href="sjar_e0203_form9.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
          <graphic id="form10" xlink:href="sjar_e0203_form10.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>F</italic>
               <sub>B</sub>
               is the buoyancy, N;
               <italic>F</italic>
               <sub>d</sub>
               is the free flow resistance force, N;
               <italic>g</italic>
               is the gravitational acceleration, m/s
               <sup>2</sup>
               ;
               <italic>V</italic>
               is the volume of ripe fruit, m
               <sup>3</sup>
               ;
               <italic>A</italic>
               is the projected area of ripe fruit perpendicular to the horizontal airflow direction, m
               <sup>2</sup>
               ;
               <italic>v</italic>
               is the relative velocity between the ripe fruit and horizontal airflow, m/s; and
               <italic>C</italic>
               <sub>D</sub>
               is the resistance coefficient which is calculated as follows:
            </p>
        <graphic id="form11" xlink:href="sjar_e0203_form11.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               Therefore,
               <italic>C</italic>
               <sub>D</sub>
               depends on the Reynolds number
               <italic>Re</italic>
               , and its equation is as follows:
            </p>
           <graphic id="form12" xlink:href="sjar_e0203_form12.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>Re</italic>
               is the Reynolds number;
               <italic>β</italic>
               is the free volume of a grid cell;
               <italic>L</italic>
               is the diameter of ripe fruit, m; and
               <italic>&#951;</italic>
               is the viscosity of horizontal airflow, kg/(m·s).
            </p>
            <p>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 <xref ref-type="table" rid="T1">Table 1</xref>.</p>
            <p>
               Based on previous research and references review (
               <xref ref-type="bibr" rid="b5">
                  González-Montellano
                  <italic>et al.</italic>
                  , 2012
               </xref>
               ;
               <xref ref-type="bibr" rid="b7">Horabik &amp; Molenda, 2016</xref>
               ;
               <xref ref-type="bibr" rid="b3">
                  Castillo-Ruiz
                  <italic>et al.</italic>
                  , 2018
               </xref>
               ;
               <xref ref-type="bibr" rid="b4">
                  Chen
                  <italic>et al.</italic>
                  , 2019
               </xref>
               ;
               <xref ref-type="bibr" rid="b27">
                  Zhao
                  <italic>et al.</italic>
                  , 2019a
               </xref>
               ,
               <xref ref-type="bibr" rid="b28">b</xref>
               ), because the properties of
               <italic>L. barbarum</italic>
               are similar to the olives, some simulation parameters of olive can be cited based on the actual situation. Based on
               <xref ref-type="bibr" rid="b5">
                  González-Montellano
                  <italic>et al.´s</italic>
                  (2012)
               </xref>
               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 (
               <xref ref-type="bibr" rid="b5">
                  González-Montellano
                  <italic>et al.</italic>
                  , 2012
               </xref>
               ), because there are irregular bulges on partial surfaces of ripe fruit and the mechanics parameters of the two conveyors material were listed in <xref ref-type="table" rid="T1">Table 1</xref>, 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 (
               <xref ref-type="bibr" rid="b9">
                  Li
                  <italic>et al.</italic>
                  , 2012
               </xref>
               ;
               <xref ref-type="bibr" rid="b11">Liu, 2018</xref>
               ), 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 (
               <italic>i.e.</italic>
               , 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 <xref ref-type="fig" rid="F4">Fig. 4</xref>.
            </p>
            <fig id="F4">
    <label>Figure 4.</label>
    <caption>
    <title>DEM simulation of ripe fruit in the horizontal
airflow.</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f04.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

            <p>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. <xref ref-type="fig" rid="F5">Fig. 5</xref> shows the trajectory of the ripe fruit which was the most easily removed by mistake.</p>
            <fig id="F5">
    <label>Figure 5.</label>
    <caption>
    <title>Trajectory of the ripe fruit.</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f05.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

         </sec>
         <sec id="S2.4">
            <title>Performance experiment of winnowing equip­ment</title>
            <p>Ningqi 7 was selected as the experiment variety. The following instruments were used: 1) an elec­tronic 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 Ins­trument Co., Ltd., China); 4) a digital anemometer (type: PM6252B; manufactured by Shenzhen New Huayi Instrument Co., Ltd., China).</p>
            <p>
               The winnowing equipment was used to remove as many impurities as possible (
               <italic>e.g.</italic>
               , unripe fruit, flowers, and leaves) in machine-harvested
               <italic>L. barbarum</italic>
               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
               <italic>Y</italic>
               <sub>1</sub>
               and the clearance rate of ripe fruit
               <italic>Y</italic>
               <sub>2</sub>
               were selected as performance indices in the experiment. Corresponding equations to calculate these indices are as follows:
            </p>
       <graphic id="form13" xlink:href="sjar_e0203_form13.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
       <graphic id="form14" xlink:href="sjar_e0203_form14.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               where
               <italic>n</italic>
               <sub>1</sub>
               is the amount of unripe fruit, flowers, and leaves before winnowing;
               <italic>n</italic>
               <sub>2</sub>
               is the amount of ri­pe fruit, unripe fruit, flowers, and leaves before winnowing;
               <italic>n</italic>
               <sub>3</sub>
               is the amount of unripe fruit, flowers, and leaves after winnowing;
               <italic>n</italic>
               <sub>4</sub>
               is the amount of ripe fruit, unripe fruit, flowers, and leaves after winnowing;
               <italic>n</italic>
               <sub>5</sub>
               is the amount of removed ripe fruit after winnowing; and
               <italic>n</italic>
               <sub>6</sub>
               is the amount of ripe fruit before winnowing.
            </p>
            <p>
               The experiment was conducted in Zhongning in the Ningxia Hui Autonomous Region (37&#176;22'56"N, 105&#176;37'21"E) on June 26, 2018; the temperature was 27.2&#176;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 input-output conveyor distance. The previous tests and simulation analysis indicated the following suitable value range for the factors: the airflow speed
               <italic>X</italic>
               <sub>1</sub>
               was 5-6 m/s; the input conveyor speed
               <italic>X</italic>
               <sub>2</sub>
               was 0.4-0.6 m/s; and the input-output conveyor distance
               <italic>X</italic>
               <sub>3</sub>
               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.
            </p>
            <p>
               Three factors and three levels were used in a quadratic orthogonal rotation design; the codes of factors are listed in <xref ref-type="table" rid="T3">Table 3</xref>, and the experimental schemes and results are presented in <xref ref-type="table" rid="T4">Table 4</xref>. 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
               <italic>L. barbarum</italic>
               from a shrub using a handheld
               <italic>L. barbarum</italic>
               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 (
               <xref ref-type="bibr" rid="b23">
                  Yuan
                  <italic>et al.</italic>
                  , 2012
               </xref>
               ;
               <xref ref-type="bibr" rid="b18">
                  Wang
                  <italic>et al.</italic>
                  , 2013
               </xref>
               ).
            </p>
            <table-wrap id="T3">
    <label>Table 3.</label>
    <caption>
    <title>Codes of factors. </title>
    </caption>
    <graphic xlink:href="sjar_e0203_t03.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</table-wrap>

<table-wrap id="T4">
    <label>Table 4.</label>
    <caption>
    <title>Experimental schemes and results. </title>
    </caption>
    <graphic xlink:href="sjar_e0203_t04.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</table-wrap>


         </sec>
      </sec>
      <sec id="S3">
         <title>Results</title>
         <sec id="S3.1">
            <title>Regression analysis</title>
            <p>The regression model of the rate of impurity change (response) using codes of all factors as variables was as follows:</p>
          <graphic id="form15" xlink:href="sjar_e0203_form15.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               ANOVA results of the rate of impurity change are shown in <xref ref-type="table" rid="T5">Table 5</xref>. The model was statistically significant (
               <italic>p</italic>
               &lt;0.0001), and
               <italic>X</italic>
               <sub>1</sub>
               ,
               <italic>X</italic>
               <sub>2</sub>
               ,
               <italic>X</italic>
               <sub>3</sub>
               ,
               <italic>X</italic>
               <sub>1</sub>
               <italic>X</italic>
               <sub>2</sub>
               ,
               <italic>X</italic>
               <sub>1</sub>
               <sup>2</sup>
               ,
               <italic>X</italic>
               <sub>2</sub>
               <sup>2</sup>
               , and
               <italic>X</italic>
               <sub>3</sub>
               <sup>2</sup>
               each had a significant effect on the rate of impurity change (
               <italic>p</italic>
               &lt;0.05). Moreover, the lack-of-fit test indicated that the model had a good fit (
               <italic>p</italic>
               =0.1569).
            </p>
            <table-wrap id="T5">
    <label>Table 5.</label>
    <caption>
    <title>ANOVA results of the rate of impurity change. </title>
    </caption>
    <graphic xlink:href="sjar_e0203_t05.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</table-wrap>

            <p>The regression model of the clearance rate of ripe fruit (response) using codes of all factors as variables was as follows:</p>
           <graphic id="form16" xlink:href="sjar_e0203_form16.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
            <p>
               ANOVA results of the clearance rate of ripe fruit are shown in <xref ref-type="table" rid="T6">Table 6</xref>. The model was statistically significant (
               <italic>p</italic>
               &lt;0.0001), and
               <italic>X</italic>
               <sub>1</sub>
               ,
               <italic>X</italic>
               <sub>2</sub>
               ,
               <italic>X</italic>
               <sub>3</sub>
               ,
               <italic>X</italic>
               <sub>1</sub>
               <italic>X</italic>
               <sub>2</sub>
               ,
               <italic>X</italic>
               <sub>1</sub>
               <sup>2</sup>
               ,
               <italic>X</italic>
               <sub>2</sub>
               <sup>2</sup>
               , and
               <italic>X</italic>
               <sub>3</sub>
               <sup>2</sup>
               each had a significant effect on the clearance rate of ripe fruit (
               <italic>p</italic>
               &lt;0.05). Moreover, the lack-of-fit test indicated that the model had a good fit (
               <italic>p</italic>
               =0.0667).
            </p>
            <table-wrap id="T6">
    <label>Table 6.</label>
    <caption>
    <title>ANOVA results of the clearance rate of ripe fruit. </title>
    </caption>
    <graphic xlink:href="sjar_e0203_t06.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</table-wrap>

         </sec>
         <sec id="S3.2">
            <title>Response surface analysis</title>
            <p>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 (<xref ref-type="disp-formula" rid="form15">Eq. [15]</xref>) are shown in <xref ref-type="fig" rid="F6">Fig. 6</xref>. As shown in <xref ref-type="table" rid="T5">Table 5</xref>, the airflow speed, input conveyor speed, and input-output conveyor distance each had an extremely significant effect on the rate of impurity change. The interaction effect between the airflow speed and input conveyor speed was significant. As shown in <xref ref-type="fig" rid="F6">Fig. 6a</xref>, as the airflow speed increased, the rate of impurity change first increased rapidly and then decreased slowly. As shown in <xref ref-type="fig" rid="F6">Fig. 6b</xref>, as the input-output conveyor distance increased, the rate of impurity change first increased slowly and then decreased slowly. As shown in <xref ref-type="fig" rid="F6">Fig. 6c</xref>, as the input conveyor speed increased, the rate of impurity change first increased slowly and then decreased slowly.</p>
            <fig id="F6">
    <label>Figure 6.</label>
    <caption>
    <title>Response surface and contour plots for effects of all factors on the rate of impurity change: <italic>X</italic><sub>1</sub> and <italic>X</italic><sub>2</sub> (a), <italic>X</italic><sub>1</sub>
and <italic>X</italic><sub>3</sub> (b), and <italic>X</italic><sub>2</sub> and <italic>X</italic><sub>3</sub> (c).</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f06.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

            <p>The response surface methodology was used to ana­lyze the effects of all factors on the clearance rate of ripe fruit; the response surface results of the regression equation (<xref ref-type="disp-formula" rid="form16">Eq. [16]</xref>) are shown in <xref ref-type="fig" rid="F7">Fig. 7</xref>. As shown in <xref ref-type="table" rid="T6">Table 6</xref>, 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 <xref ref-type="fig" rid="F7">Fig. 7a</xref>, as the airflow speed increased, the clearance rate of ripe fruit first decreased lowly and then increased rapidly; in <xref ref-type="fig" rid="F7">Fig. 7b</xref>, as the input-output conveyor distance increased, the clearance rate of ripe fruit first decreased rapidly and then increased rapidly; and in <xref ref-type="fig" rid="F7">Fig. 7c</xref>, as the input conveyor speed increased, the clearance rate of ripe fruit first decreased rapidly and then increased rapidly.</p>
            <fig id="F7">
    <label>Figure 7.</label>
    <caption>
    <title>Response surface and contour plots for effects of all factors on the clearance rate of ripe fruit: <italic>X</italic><sub>1</sub> and <italic>X</italic><sub>2</sub> (a), <italic>X</italic><sub>1</sub>
and <italic>X</italic><sub>3</sub> (b), and <italic>X</italic><sub>2</sub> and <italic>X</italic><sub>3</sub> (c).</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f07.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

            <p>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.</p>
         </sec>
         <sec id="S3.3">
            <title>Field experiment verification</title>
            <p />
            <p>
               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 input-output conveyor distance of 265.04 mm. As depicted in <xref ref-type="fig" rid="F8">Fig. 8</xref>, a high-speed camera (type: OLYMPUS i-speed TR; manufactured by Keymed (Medical &amp; 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 (
               <xref ref-type="bibr" rid="b17">
                  Torregrosa
                  <italic>et al.</italic>
                  , 2014
               </xref>
               ;
               <xref ref-type="bibr" rid="b10">
                  Liang
                  <italic>et al.</italic>
                  , 2018
               </xref>
               ). 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.
            </p>
            <fig id="F8">
    <label>Figure 8.</label>
    <caption>
    <title>Winnowing process: 0 ms (a), 40 ms (b), 80 ms (c), 120 ms (d), 160 ms (e), and 200 ms (f).</title>
    </caption>
    <graphic xlink:href="sjar_e0203_f08.jpg" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</fig>

         </sec>
      </sec>
      <sec id="S4">
         <title>Discussion</title>
         <p>
            The machine-harvested
            <italic>L. barbarum</italic>
            contained ma­ny impurities such as unripe fruit, flowers, and leaves (
            <xref ref-type="bibr" rid="b26">
               Zhang
               <italic>et al.</italic>
               , 2015
            </xref>
            ;
            <xref ref-type="bibr" rid="b19">Wang, 2018</xref>
            ;
            <xref ref-type="bibr" rid="b21">
               Xu
               <italic>et al.</italic>
               , 2018
            </xref>
            ;
            <xref ref-type="bibr" rid="b25">
               Zhang
               <italic>et al.</italic>
               , 2018b
            </xref>
            ;
            <xref ref-type="bibr" rid="b4">
               Chen
               <italic>et al.</italic>
               , 2019
            </xref>
            ;
            <xref ref-type="bibr" rid="b28">
               Zhao
               <italic>et al.</italic>
               , 2019b
            </xref>
            ). 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
            <italic>L. barbarum</italic>
            . In the present study, winnowing equipment for machine-harvested
            <italic>L. barbarum</italic>
            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 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. These results provide a design basis for future research on winnowing equipment for machine-harvested
            <italic>L. barbarum</italic>
            . But there are still some problems to be solved. For example, it is important to make miniaturized winnowing equip­ment, so as to real-time winnowing during harvesting. Furthermore, a more accurate and efficient winnowing method is also worth further studying.
         </p>
      </sec>
   </body>
   <back>
      <ref-list id="S5">
         <title>References</title>
         <ref id="b1">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Amagase</surname>
                     <given-names>H</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Farnsworth</surname>
                     <given-names>NR</given-names>
                  </name>
                  ,
               </person-group>
               <year>2011</year>
               .
               <article-title>A review of botanical characteristics, phytochemistry, clinical relevance in efficacy and safety of Lycium barbarum fruit (Goji).</article-title>
               <source>Food Res Int</source>
               <volume>44</volume>
               (
               <issue>7</issue>
               ):
               <fpage>1702</fpage>
               -
               <lpage>1717</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.foodres.2011.03.027">https://doi.org/10.1016/j.foodres.2011.03.027</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b2">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Baidu</surname>
                     <given-names>WK</given-names>
                  </name>
                  ,
               </person-group>
               <year>2015</year>
               .
               <article-title>Elastic Modulus, Poisson's Ratio, and Density of Common Materials.</article-title>
               <comment>[19 May 2015].</comment>
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://wenku.baidu.com/view/6d0a67940b4e767f5acfcecf.html">https://wenku.baidu.com/view/6d0a67940b4e767f5acfcecf.html</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b3">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Castillo-Ruiz</surname>
                     <given-names>FJ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Tombesi</surname>
                     <given-names>S</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Farinelli</surname>
                     <given-names>D</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Olive fruit detachment force against pulling and torsional stress.</article-title>
               <source>Span J Agric Res</source>
               <volume>16</volume>
               <issue>1</issue>
               <comment>e0202</comment>
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.5424/sjar/2018161-12269">https://doi.org/10.5424/sjar/2018161-12269</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b4">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Chen</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhao</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>Y</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Bu</surname>
                     <given-names>LX</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Hu</surname>
                     <given-names>GR</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>EY</given-names>
                  </name>
                  ,
               </person-group>
               <year>2019</year>
               .
               <article-title>Design and experiment on vibrating and comb brushing harvester for Lycium barbarum.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>50</volume>
               (
               <issue>1</issue>
               ):
               <fpage>152</fpage>
               -
               <lpage>161</lpage>
               ,
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.6041/j.issn.1000-1298.2019.01.016">http://dx.doi.org/10.6041/j.issn.1000-1298.2019.01.016</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b5">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>González-Montellano</surname>
                     <given-names>C</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Fuentes</surname>
                     <given-names>JM</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ayuga-Téllez</surname>
                     <given-names>E</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ayuga</surname>
                     <given-names>F</given-names>
                  </name>
                  ,
               </person-group>
               <year>2012</year>
               .
               <article-title>Determination of the mechanical properties of maize grains and olives required for use in DEM simulations.</article-title>
               <source>J Food Eng</source>
               <volume>111</volume>
               (
               <issue>4</issue>
               ):
               <fpage>553</fpage>
               -
               <lpage>562</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.jfoodeng.2012.03.017">https://doi.org/10.1016/j.jfoodeng.2012.03.017</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b6">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>He</surname>
                     <given-names>M</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Kan</surname>
                     <given-names>Z</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>CS</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Wang</surname>
                     <given-names>LH</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Yang</surname>
                     <given-names>LT</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Wang</surname>
                     <given-names>Z</given-names>
                  </name>
                  ,
               </person-group>
               <year>2017</year>
               .
               <article-title>Mechanism analysis and experiment on vibration harves­ting of wolfberry.</article-title>
               <source>T Chin Soc Agr Eng</source>
               <volume>33</volume>
               (
               <issue>11</issue>
               ):
               <fpage>47</fpage>
               -
               <lpage>53</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b7">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Horabik</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Molenda</surname>
                     <given-names>M</given-names>
                  </name>
                  ,
               </person-group>
               <year>2016</year>
               .
               <article-title>Parameters and contact models for DEM simulations of agricultural granular materials: A review.</article-title>
               <source>Biosyst Eng</source>
               <volume>147</volume>
               :
               <fpage>206</fpage>
               -
               <lpage>225</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.biosystemseng.2016.02.017">https://doi.org/10.1016/j.biosystemseng.2016.02.017</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b8">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Hu</surname>
                     <given-names>GM</given-names>
                  </name>
                  ,
               </person-group>
               <year>2010</year>
               .
               <article-title>Analysis and simulation of granular system by discrete element method using EDEM: industrial applica­tion of discrete element method and EDEM software.</article-title>
               <source>Wuhan Univ Technol Press, Wuhan, China.</source>
               <lpage>144</lpage>
            </element-citation>
         </ref>
         <ref id="b9">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Li</surname>
                     <given-names>HC</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>YM</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Gao</surname>
                     <given-names>F</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhao</surname>
                     <given-names>Z</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Xu</surname>
                     <given-names>LZ</given-names>
                  </name>
                  ,
               </person-group>
               <year>2012</year>
               .
               <article-title>CFD-DEM simulation of material motion in air-and-screen cleaning device.</article-title>
               <source>Comput Electr Agr</source>
               <volume>88</volume>
               :
               <fpage>111</fpage>
               -
               <lpage>119</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.compag.2012.07.006">https://doi.org/10.1016/j.compag.2012.07.006</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b10">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Liang</surname>
                     <given-names>N</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ni</surname>
                     <given-names>FP</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>KL</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Tang</surname>
                     <given-names>Y</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Hu</surname>
                     <given-names>YH</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Optimized installation angle and distance of a grading channel for dried jujube fruit with a push-pull actuating mechanism.</article-title>
               <source>Comput Electr Agr</source>
               <volume>150</volume>
               :
               <fpage>134</fpage>
               -
               <lpage>142</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.compag.2018.04.006">https://doi.org/10.1016/j.compag.2018.04.006</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b11">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Liu</surname>
                     <given-names>FY</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Discrete element modelling of the wheat particles and short straw in cleaning devices.</article-title>
               <source>Doctoral thesis. Northwest A&amp;F Univ., Yangling, China.</source>
            </element-citation>
         </ref>
         <ref id="b12">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Liu</surname>
                     <given-names>LY</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Hao</surname>
                     <given-names>SY</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>M</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Liu</surname>
                     <given-names>DM</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Jia</surname>
                     <given-names>FG</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Quan</surname>
                     <given-names>LZ</given-names>
                  </name>
                  ,
               </person-group>
               <year>2015</year>
               .
               <article-title>Numerical simulation and experiment on paddy ventilation resistance based on CFD-DEM.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>46</volume>
               (
               <issue>8</issue>
               ):
               <fpage>27</fpage>
               -
               <lpage>32</lpage>
               +158.
            </element-citation>
         </ref>
         <ref id="b13">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Liu</surname>
                     <given-names>P</given-names>
                  </name>
                  ,
               </person-group>
               <year>2013</year>
               .
               <article-title>Research and design of large air classifier for high impurity raw grain.</article-title>
               <source>Master's thesis. Henan Univ. Technol., Zhengzhou, China.</source>
            </element-citation>
         </ref>
         <ref id="b14">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Qiu</surname>
                     <given-names>XH</given-names>
                  </name>
                  ,
               </person-group>
               <year>2011</year>
               .
               <article-title>Machine design.</article-title>
               <source>Higher Education Press, Beijing, China.</source>
               <fpage>176</fpage>
               <lpage>292</lpage>
            </element-citation>
         </ref>
         <ref id="b15">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>So</surname>
                     <given-names>JD</given-names>
                  </name>
                  ,
               </person-group>
               <year>2001</year>
               .
               <article-title>Vibration characteristics of boxthorn (Lycium chinense Mill) branches.</article-title>
               <source>Appl Eng Agric</source>
               <volume>17</volume>
               (
               <issue>6</issue>
               ):
               <fpage>755</fpage>
               -
               <lpage>760</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.13031/2013.6834">https://doi.org/10.13031/2013.6834</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b16">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>So</surname>
                     <given-names>JD</given-names>
                  </name>
                  ,
               </person-group>
               <year>2003</year>
               .
               <article-title>Vibratory harvesting machine for boxthorn (Lycium chinense Mill) berries.</article-title>
               <source>T ASAE</source>
               <volume>46</volume>
               (
               <issue>2</issue>
               ):
               <fpage>211</fpage>
               -
               <lpage>221</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.13031/2013.12963">https://doi.org/10.13031/2013.12963</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b17">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Torregrosa</surname>
                     <given-names>A</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Albert</surname>
                     <given-names>F</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Aleixos</surname>
                     <given-names>N</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ortiz</surname>
                     <given-names>C</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Blasco</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
               </person-group>
               <year>2014</year>
               .
               <article-title>Analysis of the detachment of citrus fruits by vibration using artificial vision.</article-title>
               <source>Biosyst Eng</source>
               <volume>119</volume>
               :
               <fpage>1</fpage>
               -
               <lpage>12</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.biosystemseng.2013.12.010">https://doi.org/10.1016/j.biosystemseng.2013.12.010</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b18">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Wang</surname>
                     <given-names>CL</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Peng</surname>
                     <given-names>HB</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ding</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhao</surname>
                     <given-names>BJ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Jia</surname>
                     <given-names>F</given-names>
                  </name>
                  ,
               </person-group>
               <year>2013</year>
               .
               <article-title>Optimization for vortex pump based on response surface method.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>44</volume>
               (
               <issue>5</issue>
               ):
               <fpage>59</fpage>
               -
               <lpage>65</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b19">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Wang</surname>
                     <given-names>YL</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Research on key technology of wolfberry vibration harvest.</article-title>
               <source>Master's thesis. Northwest A&amp;F Univ., Yangling, China.</source>
            </element-citation>
         </ref>
         <ref id="b20">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Wang</surname>
                     <given-names>YL</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>Y</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Han</surname>
                     <given-names>B</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Research on laws of wolfberry dropping based on high-speed camera.</article-title>
               <source>J Agric Mechan Res</source>
               <volume>40</volume>
               (
               <issue>11</issue>
               ):
               <fpage>166</fpage>
               -
               <lpage>170</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b21">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Xu</surname>
                     <given-names>LM</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>JW</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Wu</surname>
                     <given-names>G</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Yuan</surname>
                     <given-names>QC</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ma</surname>
                     <given-names>S</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Yu</surname>
                     <given-names>CC</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Duan</surname>
                     <given-names>ZZ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Xing</surname>
                     <given-names>JJ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Liu</surname>
                     <given-names>XD</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Design and operating parameter optimization of comb brush vibratory harvesting device for wolfberry.</article-title>
               <source>T Chin Soc Agr Eng</source>
               <volume>34</volume>
               (
               <issue>9</issue>
               ):
               <fpage>75</fpage>
               -
               <lpage>82</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b22">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Yu</surname>
                     <given-names>LM</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Xu</surname>
                     <given-names>Z</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Yang</surname>
                     <given-names>JR</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Fan</surname>
                     <given-names>WB</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>N</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Long</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018</year>
               .
               <article-title>Numerical simulation of water and sediment movement in screen filter based on coupled CFD-DEM.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>49</volume>
               (
               <issue>3</issue>
               ):
               <fpage>303</fpage>
               -
               <lpage>308</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b23">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Yuan</surname>
                     <given-names>X</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Qi</surname>
                     <given-names>LJ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Wang</surname>
                     <given-names>H</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Huang</surname>
                     <given-names>SK</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ji</surname>
                     <given-names>RH</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>JH</given-names>
                  </name>
                  ,
               </person-group>
               <year>2012</year>
               .
               <article-title>Spraying parameters optimization of swing, automatic variables and greenhouse mist sprayer with response surface method.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>43</volume>
               (
               <issue>4</issue>
               ):
               <fpage>45</fpage>
               -
               <lpage>54</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b24">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhang</surname>
                     <given-names>WQ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>ZZ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Tan</surname>
                     <given-names>YZ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>W</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018a</year>
               <article-title>. Optimal design and experiment on variable pacing combing brush picking device for Lycium barbarum.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>49</volume>
               (
               <issue>8</issue>
               ):
               <fpage>83</fpage>
               -
               <lpage>90</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b25">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhang</surname>
                     <given-names>WQ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>MM</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>JX</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>W</given-names>
                  </name>
                  ,
               </person-group>
               <year>2018b</year>
               <article-title>. Design and experiment of vibrating wolfberry harvester.</article-title>
               <source>T Chin Soc Agr Machin</source>
               <volume>49</volume>
               (
               <issue>7</issue>
               ):
               <fpage>97</fpage>
               -
               <lpage>102</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b26">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhang</surname>
                     <given-names>Z</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Xiao</surname>
                     <given-names>HR</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Ding</surname>
                     <given-names>WQ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Mei</surname>
                     <given-names>S</given-names>
                  </name>
                  ,
               </person-group>
               <year>2015</year>
               .
               <article-title>Mechanism simulation analysis and prototype experiment of Lycium barbarum harvest by vibration mode.</article-title>
               <source>T Chin Soc Agr Eng</source>
               <volume>31</volume>
               (
               <issue>10</issue>
               ):
               <fpage>20</fpage>
               -
               <lpage>28</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b27">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhao</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Sugirbay</surname>
                     <given-names>A</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>Y</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Zhang</surname>
                     <given-names>S</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Liu</surname>
                     <given-names>FY</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Bu</surname>
                     <given-names>LX</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>Y</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Wang</surname>
                     <given-names>ZW</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
               </person-group>
               <year>2019a</year>
               <article-title>. FEM explicit dynamics simulation and NIR hyperspectral reflectance imaging for determination of impact bruises of Lycium barbarum L..</article-title>
               <source>Postharvest Biol Technol</source>
               <volume>155</volume>
               :
               <fpage>102</fpage>
               -
               <lpage>110</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.postharvbio.2019.05.024">https://doi.org/10.1016/j.postharvbio.2019.05.024</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b28">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhao</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>Y</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Wang</surname>
                     <given-names>YL</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
               </person-group>
               <year>2019b</year>
               <article-title>. Experimental research on parameter optimization of portable vibrating and harvesting device of Chinese wolfberry.</article-title>
               <source>J Agric Mechan Res</source>
               <volume>41</volume>
               (
               <issue>3</issue>
               ):
               <fpage>176</fpage>
               -
               <lpage>182</lpage>
               .
            </element-citation>
         </ref>
         <ref id="b29">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhao</surname>
                     <given-names>QH</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Li</surname>
                     <given-names>JJ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Yan</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Liu</surname>
                     <given-names>S</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Guo</surname>
                     <given-names>YL</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Chen</surname>
                     <given-names>DJ</given-names>
                  </name>
                  ,
                  <name>
                     <surname>Luo</surname>
                     <given-names>Q</given-names>
                  </name>
                  ,
               </person-group>
               <year>2016</year>
               .
               <article-title>Lycium barbarum polysaccharides ameliorates renal injury and inflammatory reaction in alloxan-induced diabetic nephropathy rabbits.</article-title>
               <source>Life Sci</source>
               <volume>157</volume>
               :
               <fpage>82</fpage>
               -
               <lpage>90</lpage>
               .
               <comment>
                  <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.lfs.2016.05.045">https://doi.org/10.1016/j.lfs.2016.05.045</ext-link>
               </comment>
            </element-citation>
         </ref>
         <ref id="b30">
            <element-citation publication-type="journal">
               <person-group person-group-type="author">
                  <name>
                     <surname>Zhou</surname>
                     <given-names>B</given-names>
                  </name>
                  ,
                  <name>
                     <surname>He</surname>
                     <given-names>J</given-names>
                  </name>
                  ,
               </person-group>
               <year>2010</year>
               .
               <article-title>Design of simulate hand wolfberry picking machine.</article-title>
               <source>T Chin Soc Agr Eng</source>
               <volume>26</volume>
               <issue>S1</issue>
               <fpage>13</fpage>
               <lpage>17</lpage>
            </element-citation>
         </ref>
      </ref-list>
   </back>
</article>