Calibration of Simulation Parameters for Fresh Tea Leaves Based on the Discrete Element Method

: To address the problem of a lack of accurate parameters in the discrete element simulation study of the machine-picked fresh tea leaf mechanized-sorting process, this study used machine-picked fresh tea leaves as the research object, established discrete element models of different fresh tea leaf components in EDEM software version 7.0.0. based on the bonded particle model using three-dimensional scanning inverse-modeling technology, and calibrated the simulation parameters through physical tests and virtual simulation tests. Firstly, the intrinsic parameters of machine-picked tea leaves were measured using physical tests; the physical-stacking tea leaf test was conducted using the cylinder lifting method, the tea leaf repose angle being 32.62 ◦ as measured from the stacking images using CAD. With the physical repose angle as the target value, the Plackeet–Burman test, the steepest-ascent test and the Box–Behnken optimization test were conducted in turn, and the results showed that the static friction coefficient between tea leaves, the rolling friction coefficient between tea leaves and the static friction coefficient between tea leaves and PVC have a major effect on the repose angle, and the optimal combination of the three significant parameters was determined. Finally, five simulations were conducted using the optimal combination of parameters, the relative error between the repose angle measured by the simulation test and the physical repose angle being just 0.28%. Moreover, the t -test obtained p > 0.05, indicating that there was no significant difference between the simulation test results and the physical test results. The results showed that the calibrated discrete element simulation parameters obtained could provide a reference for the discrete element simulation study of fresh tea leaves.


Introduction
Tea is rich in polyphenols, proteins, amino acids, vitamins, and other nutrients, offers good health benefits and is one of the most important cash crops in China.By the end of 2022, the area under tea plantations in China had increased to nearly 50 million mu, and tea production has continued to increase steadily [1].With the increasing scale of the tea industry and its increasing labor costs, the picking of fresh tea leaves has become an important factor limiting the industry's development.Consequently, mechanized picking has become an inevitable trend, with the current mechanized picking process being primarily a rigid, non-selective picking method, resulting in problems such as the uneven length, uneven age, and low uniformity of machine-picked fresh tea leaves [2]; this affects both the quality of finished tea and its economic benefits, making it necessary to sort the tea leaves to obtain different grades of tea leaves.In recent years, the discrete element method (DEM) has been widely used to study the kinematic behavior of bulk materials in agriculture in terms of the interaction between the bulk material and the relevant machinery [3,4].
To improve the accuracy of discrete element simulation, it is necessary to accurately establish the discrete element model of the material and accurately define its intrinsic and contact parameters [5,6].The intrinsic parameters are usually measured directly, using physical bench tests; the contact parameters of materials may be difficult to measure, resulting in inconsistent simulation and physical test results [7].Consequently, the calibration of the intrinsic and contact parameters is required for the simulation tests of agricultural material models.
To date, many studies on the calibration of discrete element parameters have been conducted for different agricultural materials, including soil [8,9], grain seeds [10][11][12], crop straws [11,[13][14][15][16], and fruit and vegetables [17,18], amongst others.Qiu et al. [9] used the Plackett-Burman design and response surface method to calibrate the discrete element simulation parameters of cinnamon soil based on the Hertz-Mindlin model with the JKR contact model.Hou et al. [10] established a discrete element model for Agropyron seed and calibrated the contact parameters using a combination of physical and simulation tests.Bart et al. [16] developed a bendable straw-stalk model to study the grain-straw separation process and conducted a sensitivity study on the mechanism influencing the crop characteristics and separation rate.Du et al. [17] used reverse engineering technology to build a discrete element model of pod pepper and calibrated its contact parameters using the response surface method.Ren et al. [19] constructed a discrete element model of sugarcane leaves using the multi-sphere aggregation method, and the contact parameters of sugarcane leaves were optimized and calibrated using the response surface method.Zhu et al. [20] used the DEM to calibrate the lunar soil simulant parameters to accurately simulate the interaction between the lunar rover wheels and the lunar soil simulant.Zhang et al. [21] established a single-root model of maize of different diameters and a maize-rootsoil mixture model, the discrete element parameters of the maize-root-soil mixture being calibrated using the response surface method.Yu et al. [22] established a discrete element model of fresh Goji berries using the multi-spherical particle method, and calibrated the discrete element contact parameters using a combination of physical and simulation tests, the physical repose angle being the target value.
In this study, the intrinsic parameters of different components of fresh tea leaves are measured using physical tests, combined with three-dimensional scanning technology.A discrete element model of different components of fresh tea leaves based on the bonded particle model was developed; using a combination of physical and simulation tests, the physical repose angle of fresh tea leaves was used as the response value, and the DEM contact parameters of fresh tea leaves were calibrated using the design of the experimental method, to conduct simulation tests for verification purposes, and to provide reliable discrete element simulation parameters for the simulation of tea mechanization operations and equipment development-such as fresh tea leaf sorting.

Geometric Measurement of Fresh Tea Leaves
To accurately develop the fresh tea leaf discrete element model, this study selected "tea 108" from the Hangzhou Institute of Agricultural Science of Tea, with one hundred single buds, one bud and one leaf, and one bud and two leaves without insect damage or pathological characteristics being randomly selected from the tea garden.The leaf length (L) and width (W) were measured using digital display vernier calipers (accuracy of 0.1 mm), separately (as shown in Figure 1), and the leaf thickness was measured using a micrometer (as shown in Figure 2); the test was repeated 3 times, the measurements of which are shown in Table 1.
According to the measured data, the average leaf length and leaf spreading of a single bud were 19.6 and 3.8 mm, respectively; the average leaf length and leaf spreading of one bud and one leaf were 29.5 and 13.8 mm, respectively; and the average leaf length and leaf spreading of one bud and two leaves were 45.2 and 29.3 mm, respectively.Among them, the thickness of the tea leaves is generally in the range of 0.15~0.3mm.

Density Measurement of Fresh Tea Leaves
Fresh machine-picked tea leaves belong to bulk materials, and the density measured in this paper is the true density, the measurement standard being based on the measurement method of solid density in the GB/T4472-2011 standard [23].First, an electronic balance with an accuracy of 0.01 g was used to weigh the fresh tea leaf mass (M) of 2~5 g, after which a measuring cylinder with an accuracy of 0.1 mL was used to measure the According to the measured data, the average leaf length and leaf spreading of a single bud were 19.6 and 3.8 mm, respectively; the average leaf length and leaf spreading of one bud and one leaf were 29.5 and 13.8 mm, respectively; and the average leaf length and leaf spreading of one bud and two leaves were 45.2 and 29.3 mm, respectively.Among them, the thickness of the tea leaves is generally in the range of 0.15~0.3mm.

Density Measurement of Fresh Tea Leaves
Fresh machine-picked tea leaves belong to bulk materials, and the density measured in this paper is the true density, the measurement standard being based on the measurement method of solid density in the GB/T4472-2011 standard [23].First, an electronic balance with an accuracy of 0.01 g was used to weigh the fresh tea leaf mass (M) of 2~5 g, after which a measuring cylinder with an accuracy of 0.1 mL was used to measure the  According to the measured data, the average leaf length and leaf spreading of a single bud were 19.6 and 3.8 mm, respectively; the average leaf length and leaf spreading of one bud and one leaf were 29.5 and 13.8 mm, respectively; and the average leaf length and leaf spreading of one bud and two leaves were 45.2 and 29.3 mm, respectively.Among them, the thickness of the tea leaves is generally in the range of 0.15~0.3mm.

Density Measurement of Fresh Tea Leaves
Fresh machine-picked tea leaves belong to bulk materials, and the density measured in this paper is the true density, the measurement standard being based on the measurement method of solid density in the GB/T4472-2011 standard [23].First, an electronic balance with an accuracy of 0.01 g was used to weigh the fresh tea leaf mass (M) of 2~5 g, after which a measuring cylinder with an accuracy of 0.1 mL was used to measure the pure water volume (V 0 ).The fresh tea leaves were completely immersed in the measuring cylinder using a glass rod, and the total volume (V 1 ) of water and fresh tea leaves in the measuring cylinder was recorded (as shown in Figure 2).The difference between the two values was the volume of fresh tea leaves.Based on the principle of solid density measurement, we could calculate the fresh tea leaf density using Equation (1), repeated 10 times to obtain an average value, the density of fresh tea leaves being 851.4 kg m −3 .The calculation as shown in Equation ( 1) is the following: where M denotes the mass of fresh tea leaves (kg), V 0 denotes the volume of water in the measuring cylinder without the addition of fresh tea leaves (m 3 ), V 1 denotes the total volume in the measuring cylinder after adding the fresh tea leaves, and ρ denotes the density of the fresh tea leaves (kg m −3 ).

Modulus of Elasticity and Shear Modulus
The modulus of elasticity is one of the important parameters in discrete element simulation [24].Since tea leaves are thin and soft materials, it is impossible to conduct compression tests on them, so their stems were cut into small circular segments of length 3 mm using a utility knife, and uniaxial compression tests were conducted using a TMS-PRO mass spectrometer (Zhuohao Laboratory Equipment Co., Ltd, Shanghai, China), as shown in Figure 3 [25].During the test, the diameter of the stem segment was measured using a vernier caliper, after which the stem segment was placed vertically on the loading table of the mass spectrometer, and a 38.1 mm-diameter circular probe was used to set the compression speed to 10 mm/min and the loading displacement to 1 mm, the probe being automatically returned to its original position at the end of the test.The test was repeated 10 times, the average value of the elastic modulus of tea leaves being 9.24 MPa, the shear modulus being 3.3 MPa, and Poisson's ratio being 0.4 using Equation (2), Equation (3), and Equation (4), respectively.The calculation as shown in Equations ( 2)-( 4) is the following: where E is the modulus of elasticity (MPa), F is the axial load applied to the fresh leaf stem (N), A is the contact area (mm 2 ), ε is the strain, G is the shear modulus (MPa), µ is Poisson's ratio, l 0 is the pre-test length of the fresh leaf stem (mm), d 0 is the pre-test diameter of the fresh leaf stem (mm), ∆l is the change in the length of the fresh leaf stem after the test, ∆d is the change in diameter of the fresh leaf stem after the test (mm), ε 2 is the transverse strain, and ε 1 is the longitudinal strain.

Measurement of Physical Angle of Repose
The repose angle of fresh tea leaves is one of the more important indicators for characterizing the macroscopic properties of tea materials-including the flow properties and internal friction characteristics [26].In this study, a bulk material stacking-angle measurement method was used to determine the actual repose angle of fresh tea leaves based on the JB/T9014.7-1999standard [23], using the cylinder lifting method (Figure 4).The test sample was machine-picked fresh tea leaves; using a PVC hollow cylinder of inner diameter 100 mm and height 200 mm, after filling the cylinder with fresh tea leaves and then lifting the cylinder at a constant speed of 0.2 m/s, so that the tea leaves fall out of the cylinder naturally, a pile of tea fresh leaf materials forms on the horizontal bottom plate, and the angle between the bus bar and the horizontal plane of the accumulation is measured.We repeated the above test 10 times to obtain an average value of the repose angle of fresh tea leaves, which was 32.62

Measurement of Physical Angle of Repose
The repose angle of fresh tea leaves is one of the more important indicators for characterizing the macroscopic properties of tea materials-including the flow properties and internal friction characteristics [26].In this study, a bulk material stacking-angle measurement method was used to determine the actual repose angle of fresh tea leaves based on the JB/T9014.7-1999standard [23], using the cylinder lifting method (Figure 4).The test sample was machine-picked fresh tea leaves; using a PVC hollow cylinder of inner diameter 100 mm and height 200 mm, after filling the cylinder with fresh tea leaves and then lifting the cylinder at a constant speed of 0.2 m/s, so that the tea leaves fall out of the cylinder naturally, a pile of tea fresh leaf materials forms on the horizontal bottom plate, and the angle between the bus bar and the horizontal plane of the accumulation is measured.We repeated the above test 10 times to obtain an average value of the repose angle of fresh tea leaves, which was 32.62°.

Measurement of Physical Angle of Repose
The repose angle of fresh tea leaves is one of the more important indicators for characterizing the macroscopic properties of tea materials-including the flow properties and internal friction characteristics [26].In this study, a bulk material stacking-angle measurement method was used to determine the actual repose angle of fresh tea leaves based on the JB/T9014.7-1999standard [23], using the cylinder lifting method (Figure 4).The test sample was machine-picked fresh tea leaves; using a PVC hollow cylinder of inner diameter 100 mm and height 200 mm, after filling the cylinder with fresh tea leaves and then lifting the cylinder at a constant speed of 0.2 m/s, so that the tea leaves fall out of the cylinder naturally, a pile of tea fresh leaf materials forms on the horizontal bottom plate, and the angle between the bus bar and the horizontal plane of the accumulation is measured.We repeated the above test 10 times to obtain an average value of the repose angle of fresh tea leaves, which was 32.62°.Owing to the irregular appearance of fresh tea leaves, conventional modeling methods cannot accurately restore their real characteristics.The aim is to establish more accurate 3D models of fresh tea leaves and improve the accuracy of simulation tests while also considering the limitations of non-spherical particle modeling using EDEM software version 7.0.0., to reduce the simulation time and calculation of overheads, as shown in Figure 5. Fresh tea leaves-namely, a single bud, one bud and one leaf, and one bud and two leaveswhose leaf length and leaf spreading were close to the average values were selected as the research object.rate 3D models of fresh tea leaves and improve the accuracy of simulation tests wh considering the limitations of non-spherical particle modeling using EDEM software 7.0.0., to reduce the simulation time and calculation of overheads, as shown in Fi Fresh tea leaves-namely, a single bud, one bud and one leaf, and one bud and two le whose leaf length and leaf spreading were close to the average values were selected research object.As is evident from Figure 5, the Tianyuan 3D OKIO scanner is used to scan th contour of the fresh tea leaves using blue light non-contact photography, to accu obtain the three-dimensional coordinates of the outer surface of the tea leaves, and their point cloud data.The point cloud data were then imported into Geomagic software version 2017 for inverse modeling, and the operations of coloring, noise tion, point cloud triangulation, merging, smoothing and model correction were con in turn to obtain a 3D solid model.Finally, the solid model was imported into Hype software version 2021 for meshing to obtain the tea mesh model.

Discrete Element Model
The tea leaf discrete element model adopts a multi-sphere bonded particle comprising several spherical particles of equal diameter bonded by "Bond" to si the characteristics similar to those of real tea leaves.The smaller the radius of the sp particles used, the larger the number of bonded spherical particles, and the smal radius of the spherical particles used, the closer they are to the tea leaves' profile however, the time cost of simulation greatly increases.The fresh tea leaf 3D contour is imported into EDEM 2021 software as a geometry, and the spherical particle mat added, a spherical particle of radius 0.25 mm being used for filling, using the pre test.After the filling is completed, the fresh tea leaf discrete element model As is evident from Figure 5, the Tianyuan 3D OKIO scanner is used to scan the outer contour of the fresh tea leaves using blue light non-contact photography, to accurately obtain the three-dimensional coordinates of the outer surface of the tea leaves, and obtain their point cloud data.The point cloud data were then imported into Geomagic Wrap software version 2017 for inverse modeling, and the operations of coloring, noise reduction, point cloud triangulation, merging, smoothing and model correction were conducted in turn to obtain a 3D solid model.Finally, the solid model was imported into Hypermesh software version 2021 for meshing to obtain the tea mesh model.

Discrete Element Model
The tea leaf discrete element model adopts a multi-sphere bonded particle model comprising several spherical particles of equal diameter bonded by "Bond" to simulate the characteristics similar to those of real tea leaves.The smaller the radius of the spherical particles used, the larger the number of bonded spherical particles, and the smaller the radius of the spherical particles used, the closer they are to the tea leaves' profile model; however, the time cost of simulation greatly increases.The fresh tea leaf 3D contour model is imported into EDEM 2021 software as a geometry, and the spherical particle material is added, a spherical particle of radius 0.25 mm being used for filling, using the pre-filling test.After the filling is completed, the fresh tea leaf discrete element model can be obtained, as shown in Figure 6.In the EDEM post-processing interface, the radius size, ID number, and spherical center coordinate parameters of spherical particles are exported.According to statis, the number of filled particles of discrete element models of different components of fresh tea leaves-that is, the single bud, one bud and one leaf, and one bud and two leaves-are 284, 2286 and 4608, respectively.
Since tea leaves belong to flexible materials, they easily produce a certain degree of bending after collision in the process of movement.To improve the realism of fresh tea leaf simulations, it is necessary to build a flexible discrete element model and adopt the meta-particle model of Bonding V2.The contact radius of the discrete element multi-sphere model can be set to detect whether the particles are bonded, and when the contact radius of the spheres detects contact a "Bond" key is generated to bond the two spheres together.If the contact radius is too small, the bonding model is brittle, and if the contact radius is too large, a bond will be generated between non-contacting particles in the bonding model.To reduce the effect of the contact radius on the simulation results, the contact radius should be 20-30% larger than the physical radius [23].
lture 2024, 14, x FOR PEER REVIEW 7 of 15 obtained, as shown in Figure 6.In the EDEM post-processing interface, the radius size, ID number, and spherical center coordinate parameters of spherical particles are exported.According to statis, the number of filled particles of discrete element models of different components of fresh tea leaves-that is, the single bud, one bud and one leaf, and one bud and two leaves-are 284, 2286 and 4608, respectively.
(a) (b) (c) Since tea leaves belong to flexible materials, they easily produce a certain degree of bending after collision in the process of movement.To improve the realism of fresh tea leaf simulations, it is necessary to build a flexible discrete element model and adopt the meta-particle model of V2.The contact radius of the discrete element multisphere model can be set to detect whether the particles are bonded, and when the contact radius of the spheres detects contact a "Bond" key is generated to bond the two spheres together.If the contact radius is too small, the bonding model is brittle, and if the contact radius is too large, a bond will be generated between non-contacting particles in the bonding model.To reduce the effect of the contact radius on the simulation results, the contact radius should be 20-30% larger than the physical radius [23].
Based on the coordinate information of spherical particles filled with fresh tea leaves mentioned above, the physical radius of particles was set to 0.25 mm and the contact radius was set to 0.3 mm.Combined with the discrete element simulation parameters of fresh tea leaves and PVC in the literature [27,28], the range of discrete element simulation parameters in this study are as shown in Table 2. Based on the data obtained from the compression test, after calculation and multiple simulation adjustments, the Bonding parameters of the discrete elements of fresh tea leaves as set in Table 3 could be obtained, and a flexible bond model of discrete elements established (as shown in Figure 7).Based on the coordinate information of spherical particles filled with fresh tea leaves mentioned above, the physical radius of particles was set to 0.25 mm and the contact radius was set to 0.3 mm.Combined with the discrete element simulation parameters of fresh tea leaves and PVC in the literature [27,28], the range of discrete element simulation parameters in this study are as shown in Table 2. Based on the data obtained from the compression test, after calculation and multiple simulation adjustments, the Bonding parameters of the discrete elements of fresh tea leaves as set in Table 3 could be obtained, and a flexible bond model of discrete elements established (as shown in Figure 7).dius was set to 0.3 mm.Combined with the discrete element simulation parameters of fresh tea leaves and PVC in the literature [27,28], the range of discrete element simulation parameters in this study are as shown in Table 2. Based on the data obtained from the compression test, after calculation and multiple simulation adjustments, the Bonding parameters of the discrete elements of fresh tea leaves as set in Table 3 could be obtained, and a flexible bond model of discrete elements established (as shown in Figure 7).

DEM Simulation Test
In the EDEM 2021 software, a hollow cylinder (of inner diameter 100 mm and height 200 mm) was added, and a virtual plane created above the cylinder as a particle plant.The fresh tea leaf simulation stacking test is shown in Figure 8.

DEM Simulation Test
In the EDEM 2021 software, a hollow cylinder (of inner diameter 100 mm and height 200 mm) was added, and a virtual plane created above the cylinder as a particle plant.The fresh tea leaf simulation stacking test is shown in Figure 8.The particle generation method is dynamic; the generation rates of the single bud, one bud and one leaf, and one bud and two leaves are set to 0.002, 0.018, and 0.035 kg s −1 , respectively, and the total generated fresh tea leaf mass is 55 g.The total simulation time is 3 s, and when the fresh tea leaf particles are stabilized, the hollow cylinder is lifted at a speed of 0.2 m s −1 , and the particles form a material pile when they are stationary on the plane.The particle generation method is dynamic; the generation rates of the single bud, one bud and one leaf, and one bud and two leaves are set to 0.002, 0.018, and 0.035 kg s −1 , respectively, and the total generated fresh tea leaf mass is 55 g.The total simulation time is 3 s, and when the fresh tea leaf particles are stabilized, the hollow cylinder is lifted at a speed of 0.2 m s −1 , and the particles form a material pile when they are stationary on the plane.
Using Protractor, a built-in angle measurement tool in the software, the repose angle is measured in both the +X and +Y directions using the center of the stacked body as the measurement origin, and the results are averaged.The simulated repose-angle measurement diagram is shown in Figure 9.
Using Protractor, a built-in angle measurement tool in the software, the repose angle is measured in both the +X and +Y directions using the center of the stacked body as the measurement origin, and the results are averaged.The simulated repose-angle measurement diagram is shown in Figure 9.

Analysis of the Simulation Results of the Plackett-Burman Test
The Plackett-Burman test was used to determine the significance of each factor by comparing the differences between the two levels of each factor based on the relationship between the target response and each factor, to quickly screen out the factors that had a significant effect on the response values [29].In this study, the Plackett-Burman test design was conducted using Design-Expert software 10.0, and the physical repose angle of fresh tea leaves was used as the response value to screen out the factors that had significant effects on the repose angle of machine-harvested fresh tea leaves.Each test parameter was set at two levels-that is, high (+1) and low (−1), denoted by X1 to X6.Based on the existing literature and a large number of simulation pre-tests, the value range of each simulation parameter was determined, as shown in Table 4, and a group of central point tests was added, totaling 13 test groups.The Plackett-Burman test protocol and results are shown in Table 5, and the analysis of variance (ANOVA) of the test results is shown in Table 6.As is evident from Table 6, the model shows p < 0.05 and the coefficient of determination R 2 = 0.8842, indicating that the regression model is significant and can predict the trends of each parameter well.Moreover, p < 0.01 for X2 (tea-tea static friction coefficient) and X3 (tea-tea rolling friction coefficient), indicating that X2 and X3 have an extremely significant effect on the formation of the repose angle; p < 0.05 for X5 (tea-PVC static friction coefficient), indicating that X5 has a significant effect on the formation of the repose angle.By comparing the magnitude of the contribution of each parameter, the order of the effect of each parameter on the repose angle was X2, X3, X5, X1, X4, and X6.
In the subsequent steepest-ascent test and Box-Behnken test, only the three most significant parameters were considered.The standardized effects of X2, X3, and X5 were all greater than 0, so their effects on the repose angle were positive.Consequently, in the

Analysis of the Simulation Results of the Plackett-Burman Test
The Plackett-Burman test was used to determine the significance of each factor by comparing the differences between the two levels of each factor based on the relationship between the target response and each factor, to quickly screen out the factors that had a significant effect on the response values [29].In this study, the Plackett-Burman test design was conducted using Design-Expert software 10.0, and the physical repose angle of fresh tea leaves was used as the response value to screen out the factors that had significant effects on the repose angle of machine-harvested fresh tea leaves.Each test parameter was set at two levels-that is, high (+1) and low (−1), denoted by X 1 to X 6 .Based on the existing literature and a large number of simulation pre-tests, the value range of each simulation parameter was determined, as shown in Table 4, and a group of central point tests was added, totaling 13 test groups.The Plackett-Burman test protocol and results are shown in Table 5, and the analysis of variance (ANOVA) of the test results is shown in Table 6.Repose Angle ( • ) As is evident from Table 6, the model shows p < 0.05 and the coefficient of determination R 2 = 0.8842, indicating that the regression model is significant and can predict the trends of each parameter well.Moreover, p < 0.01 for X 2 (tea-tea static friction coefficient) and X 3 (tea-tea rolling friction coefficient), indicating that X 2 and X 3 have an extremely significant effect on the formation of the repose angle; p < 0.05 for X 5 (tea-PVC static friction coefficient), indicating that X 5 has a significant effect on the formation of the repose angle.By comparing the magnitude of the contribution of each parameter, the order of the effect of each parameter on the repose angle was X 2 , X 3 , X 5 , X 1 , X 4 , and X 6 .
In the subsequent steepest-ascent test and Box-Behnken test, only the three most significant parameters were considered.The standardized effects of X 2 , X 3 , and X 5 were all greater than 0, so their effects on the repose angle were positive.Consequently, in the subsequent steepest-ascent test, the factors that showed positive effects were gradually increased in fixed increments.

Analysis of the Steepest-Ascent Test Simulation Results
According to the significant parameters screened using the Plackett-Burman test, the parameters are increased or decreased in certain steps based on the degree of parameter influence: the parameters with less influence are taken up to the middle level of the values in Table 4 for the steepest-ascent test, and the relative error Y between the repose angle θ of the physical test and the repose angle θ ′ of the steepest-ascent test is calculated using Equation ( 5), so that the nearby area of the optimal value can be determined quickly.The calculation is as shown in Equation ( 5): The design and results of the steepest-ascent test program are shown in Table 7, with the repose-angle relative error tending to decrease first before increasing, and the reposeangle relative error being the smallest under Test No. 3. Consequently, the parameters of Test No. 3 were used as the center point in the subsequent tests, and the parameters of Tests No. 2 and No. 4 were used as the low and high levels for response surface design, respectively.The optimal values of the three parameters are X 2 = 0.865, X 3 = 0.067, and X 3 = 0.737.To verify the accuracy of the optimal parameter combination, five repose-angle simulation tests were conducted with the above optimized parameter combinations as EDEM simulation parameters, the simulated repose angle measurements being 32.21 • , 32.74 • , 32.83 • , 32.56 • , and 32.33 • , with an average value of 32.53 • and a relative error of 0.28% with respect to the physical repose angle.The t-test was applied to analyze the simulated repose angle and the physical repose angle, and p = 0.506 > 0.05 was obtained, indicating that there was no significant difference between the simulated repose angle and the physical repose angle.

Conclusions
Using a 3D scanner and reverse engineering technology, a discrete element model of representative fresh tea leaves with different components was developed, which could provide a reference for DEM modeling of irregularly shaped crops.The intrinsic parameters of fresh tea leaves were measured using physical tests, and the Plackett-Burman test was conducted based on the simulation parameters obtained from physical tests and the related literature.The results were analyzed using ANOVA, and the parameters with a significant effect on the repose angle were obtained as the tea-tea static friction coefficient, the tea-tea rolling friction coefficient, and the static friction coefficient between tea and PVC, with the steepest-climb test being conducted to determine the range of significant parameters.Based on the Box-Behnken test results, a quadratic polynomial regression model of three significance parameters and a repose angle was determined, the repose angle of fresh tea leaves being used as the target value for the regression equation to find the optimal combination of the significance parameters; the best combination of the significance parameters was 0.865 for the static friction coefficient of tea leaves, 0.067 for the rolling friction coefficient of tea leaves, and 0.737 for the static friction coefficient between tea leaves and PVC.A t-test was then conducted and p > 0.05 was obtained, indicating that the simulation results were not significantly different from the physical test results, further verifying the reliability of the simulated parameter combinations.The results showed that the contact parameter calibration method for fresh tea leaves was feasible and could provide a basis for the discrete element study of machine-picked fresh tea leaf sorting and other work.

Figure 2 .
Figure 2. Density measurement of fresh tea leaves: (a) quality measurement of tea leaves; (b) volume measurement of tea leaves.

Figure 2 .
Figure 2. Density measurement of fresh tea leaves: (a) quality measurement of tea leaves; (b) volume measurement of tea leaves.

Figure 4 .
Figure 4. Physical repose angle test of fresh tea leaves.

Figure 4 .
Figure 4. Physical repose angle test of fresh tea leaves.Figure 4. Physical repose angle test of fresh tea leaves.

Figure 4 .
Figure 4. Physical repose angle test of fresh tea leaves.Figure 4. Physical repose angle test of fresh tea leaves.

Figure 5 .
Figure 5.The acquisition process of a fresh tea leaf 3D model.

Figure 5 .
Figure 5.The acquisition process of a fresh tea leaf 3D model.

Figure 6 .
Figure 6.Discrete element model of fresh tea leaves: (a) single bud; (b) one bud and one leaf; (c) one bud and two leaves.

Figure 7 .
Figure 7. BondingV2 model of fresh tea leaves: (a) single bud (b) one bud and one leaf (c) one bud and two leaves.

Figure 6 .
Figure 6.Discrete element model of fresh tea leaves: (a) single bud; (b) one bud and one leaf; (c) one bud and two leaves.

Figure 7 .
Figure 7. BondingV2 model of fresh tea leaves: (a) single bud (b) one bud and one leaf (c) one bud and two leaves.Figure 7. BondingV2 model of fresh tea leaves: (a) single bud (b) one bud and one leaf (c) one bud and two leaves.

Figure 7 .
Figure 7. BondingV2 model of fresh tea leaves: (a) single bud (b) one bud and one leaf (c) one bud and two leaves.Figure 7. BondingV2 model of fresh tea leaves: (a) single bud (b) one bud and one leaf (c) one bud and two leaves.

Figure 8 .
Figure 8. Simulation test of the repose angle.

Figure 8 .
Figure 8. Simulation test of the repose angle.

Figure 9 .
Figure 9. Measurement of repose angle in simulation test: (a) tea leaf stacking +X direction, (b) tea leaf stacking +Y direction.

Figure 9 .
Figure 9. Measurement of repose angle in simulation test: (a) tea leaf stacking +X direction, (b) tea leaf stacking +Y direction.

Table 5 .
Design and results of Plackett-Burman test.

Table 6 .
Analysis of significance of parameters in Plackett-Burman test.

Table 7 .
Design and results of steepest-ascent test.