Abstract
In order to ensure the reliability of the simulation test on the shelling process of ginkgo nuts, the physical parameters (characteristic dimensions, thousand grain mass, moisture content, density, Poisson’s ratio, elastic modulus, and shear modulus), contact parameters (restitution coefficient, static friction coefficient, and dynamic friction coefficient) and actual stacking angle of the ginkgo nuts were measured by physical test methods. A stacking angle test model was established, and simulation tests were performed. Parameter intervals for the simulation were selected according to the measured parameters. Three factors affecting the stacking angle significantly were determined using the Plackett–Burman test: static and dynamic friction coefficients of ginkgo nut–ginkgo nut and the restitution coefficient of the ginkgo nut–steel plate. Each factor’s range was selected using the steepest climbing test. A quadratic polynomial regression model of the stacking angle considering these factors was obtained using the Box–Behnken test. Taking the actual stacking angle of 39.59° as the optimization objective, the optimal combination was formed: the static and dynamic friction coefficients were 0.5 and 0.36, respectively, and the restitution coefficient was 0.42. This combination was used for the verification test. The stacking angle was found to be 41.12°. The relative error was 3.86%, indicating that the model is reliable and can provide a useful reference for relevant research.
Download PDF
Full Article
Experimental Calibration of Parameters of Discrete Element Model for Ginkgo Nut
Minji Liu,a Jiannan Wang,a ZhaoyanYou,a Yanghua Wu,a Ni Wang,b,* and Huanxiong Xie a,*
In order to ensure the reliability of the simulation test on the shelling process of ginkgo nuts, the physical parameters (characteristic dimensions, thousand grain mass, moisture content, density, Poisson’s ratio, elastic modulus, and shear modulus), contact parameters (restitution coefficient, static friction coefficient, and dynamic friction coefficient) and actual stacking angle of the ginkgo nuts were measured by physical test methods. A stacking angle test model was established, and simulation tests were performed. Parameter intervals for the simulation were selected according to the measured parameters. Three factors affecting the stacking angle significantly were determined using the Plackett–Burman test: static and dynamic friction coefficients of ginkgo nut–ginkgo nut and the restitution coefficient of the ginkgo nut–steel plate. Each factor’s range was selected using the steepest climbing test. A quadratic polynomial regression model of the stacking angle considering these factors was obtained using the Box–Behnken test. Taking the actual stacking angle of 39.59° as the optimization objective, the optimal combination was formed: the static and dynamic friction coefficients were 0.5 and 0.36, respectively, and the restitution coefficient was 0.42. This combination was used for the verification test. The stacking angle was found to be 41.12°. The relative error was 3.86%, indicating that the model is reliable and can provide a useful reference for relevant research.
DOI: 10.15376/biores.19.4.8154-8172
Keywords: Ginkgo nut; Stacking angle; Discrete element; Parameter calibration; Experiment
Contact information: a: Nanjing Institute of Agricultural Mechanization, Ministry of Agriculture and Rural Affairs, Nanjing, 210014, China; b: Jiangsu Academy of Agricultural Sciences, Nanjing 210014, China;
* Corresponding author: wangni0209@163.com; xhxnq@163.com
INTRODUCTION
Ginkgo biloba (gingko) is planted widely in China (Liu et al. 2022). The number of cultivated plants exceeds 2.5 billion and is still increasing rapidly every year (Zhang et al. 2022a). Ginkgo biloba is an important biomass resource. Many studies have been carried out on the extraction of biomass resources such as cellulose and lignin from gingko trees, leaves, and shells (Usuki et al. 2011; Andersson et al. 2015; Jiang et al. 2018). The ginkgo nutshell contains about 44.3% cellulose, 46.3% lignin, and 7.0% hemicelluloses (Ni et al. 2020), which has been successfully used as a raw material for the preparation of activated carbon and biosorbents (Crouvisier-Urio et al. 2017; Qin et al. 2019). Moreover, ginkgo nutshell is widely used in medicine (Nowak et al. 2021; Shen et al. 2021) and chemistry (Jiang et al. 2018; Hong et al. 2021).
Shelling is the routine method for obtaining ginkgo nutshell. The shell is subjected to mechanical effects such as extrusion and friction during the shelling process. Therefore, the study of the extrusion process is beneficial to obtain a reasonable combination of shelling parameters and thereby improve the shelling quality. Discrete element method (DEM) is an effective method to study particulate materials, which has been widely used in different fields. Discrete element simulation can simulate the extrusion process of ginkgo nut, which provides a useful reference for force analysis (Yu et al. 2005).
For discrete element simulation, the intrinsic parameters and contact parameters of the simulation model need to be defined. The intrinsic parameters include density, shear modulus, etc., and the values are consistent with the real values. Contact parameters include restitution coefficient, as well as the static friction coefficients between particles and between particles and contact materials. Due to the difference in geometry between the particle simulation model and the real particles, the simulation contact parameters are in error from the real values, and the contact parameters need to be calibrated. Parameter calibration can effectively ensure the reliability of simulation tests (Horabik et al. 2016; Coetzee. 2017; Hou et al. 2020). Li et al. (2019) increased the average particle size of the wheat particle simulation model and calibrated the relevant parameters of wheat particles. Wen et al. (2020) proposed a simulation parameter calibration method by comparing the measurement results of the repose angle using different particle characteristic test methods. Liu et al. (2018a) conducted uniaxial compression and impact tests by combining a nonlinear viscoelastic contact model with a hysteretic model to calibrate the normal damping parameters between wheat grains, and between wheat grains and steel plates. Liu et al. (2018b) optimized and calibrated the simulation parameters of micro-potatoes by combining experimental measurements and simulation tests. Santos et al. (2015) used the DEM to measure the repose angle of a cherry using a central combination test design and experiment; they then optimized the relevant parameters of the DEM simulation. Coetzee et al. (2010) determined the particle stiffness through compression tests and obtained relevant parameters by trial. The DEM modeling process can be simplified by calibrated simulation parameters.
In summary, several researchers have conducted parameter calibration for wheat, potato, cherry, and other foods. However, most of the existing studies have been conducted on materials having relatively regular shapes and uniform compositions. The composition of ginkgo kernels and shells varies greatly. Meanwhile, due to the gap between the kernel and shell, they can move relative to each other. Therefore, compared to the existing research objects, the properties of ginkgo nut are more complex. There have been few studies on ginkgo nuts. Ch’ng et al. (2013) tested some physical characteristics data of ginkgo nuts at a moisture content of 45.53% (±2.07) (wet basis) and ginkgo kernels at 60.13% (±2.00) (wet basis). At the same time, the coefficient of static friction for nuts and kernels was determined by using plywood, glass, rubber, and galvanized steel sheet. The study provided useful inputs for this study.
In this study, the ginkgo nut discrete element simulation parameters were calibrated by comparing the simulated value and the physical value of the stacking angle in order to improve the reliability of the simulation test of ginkgo nut shelling process. It is important for the design and optimization of the key components of the sheller, which is conducive to improving the quality and efficiency of shelling. There are two main contributions: (1) Compared with materials with uniform composition, the properties of ginkgo nut are more complex, and the parameter calibration study has been rarely reported. This paper provides a useful reference for the study of materials having such complex aspects. (2) Due to the special properties of ginkgo nut, there is a large error in measuring the restitution coefficient using the traditional oblique collision method. In this paper, with the help of a high-speed camera, the velocities of ginkgo nuts before and after collision with the horizontal contact surface were measured, which could significantly reduce the calculation error.
EXPERIMENTAL
Materials
The ginkgo nut breed selected for this study was “Da Fo Zhi”. It was obtained from Liuchen Town, Taixing City, Jiangsu Province, China. Taixing is a famous ginkgo-producing area and plays a vital role in China. The growth period of ginkgo nuts is about 90 days. In this study, ginkgo nuts were selected for the experiment after 85 days of growth. When the fruit had matured, it was harvested manually, peeled, and dried. The selected nuts are spindle-shaped, and their coated shells are compact. The chemical composition of the shell of ginkgo nut used in this study is shown in Table 1. The results differed from other studies due to differences in varieties, growth periods, etc.
Table 1. Chemical Composition of the Ginkgo Nutshell
Measurement of Physical Parameters
Measurement of basic physical parameters
A thousand ginkgo nuts were selected randomly. Then, the characteristic dimensions were measured by a vernier caliper (Harbin Measuring & Cutting Tool Group Company Ltd., Harbin, China). The test method drew on the relevant studies (Hou et al. 2020). The characteristic dimensions of the ginkgo nut are shown in Fig. 1.
Fig. 1. Characteristic dimensions of the ginkgo nut
The thousand grain mass of the ginkgo nut was measured according to GB/T 5519 (2018) standard. Five hundred intact ginkgo nuts were randomly selected and weighed with an electronic balance (Puchun Measure Instrument Co., Ltd., Shanghai, China). The moisture content was measured using an electronic balance and an electric convection oven (Nanjing Experimental Instrument Factory, Nanjing, China) according to GB/T 3543.6 (1995) standard. The particle density was obtained using an electronic balance and a measuring cylinder (Jiangsu Sanaisi scientific Co., Ltd., Yancheng, China). Each parameter was tested ten times and averaged. The results are presented in Table 2.
Table 2. Basic Physical Parameters of the Ginkgo Nut
Measurement of Poisson’s Ratio
Poisson’s ratio was measured according to GB/T 22315 (2008). Ten ginkgo nuts were randomly selected, and their original dimensions in the Z and X directions were recorded. An extrusion test of the ginkgo nuts was carried out using an electronic universal testing machine (Shenzhen Suns Technology Stock Co., Ltd., Shenzhen, China). The Poisson’s ratio was obtained by testing the dimensional change of a ginkgo nut at the cracking limit in the Z and X directions before and after pressure. In the test, the plane indenter was used to pressure along the X-axis direction at a constant speed (0.2 mm/s), and the pressure was stopped when the ginkgo shell broke. The vernier caliper recorded the deformation at the load crack limit along the Z-axis direction, and the universal testing machine recorded the deformation along the X-axis direction. The changes in X-axis and Z-axis directions were measured individually with a testing accuracy of 0.02 mm. The Poisson’s ratio of ginkgo nuts was obtained using Eq. 1,
where μ is the Poisson’s ratio, εx is the X-axis strain, εz is the Z-axis strain, ΔX is the X-axis deformation (mm), X is the X-axis dimension (mm), ΔZ is the Z-axis deformation (mm), and Z is the Z-axis dimension (mm). The Poisson’s ratio was tested ten times and averaged. The Poisson’s ratio was 0.34 and the standard deviation was 0.018. Poisson’s ratio is an important property of agricultural materials, usually within a certain range that varies with the type and condition of the material. Therefore, the Poisson’s ratio obtained in this study differs from other agricultural materials (Liu et al. 2018b; Hou et al. 2020).
Measurement of elastic modulus and shear modulus
The Hertz contact stress method was used to determine the elastic modulus of ginkgo nuts according to ASAE S368.4 (2017). Ten ginkgo nuts were randomly selected, and their characteristic dimensions were tested before the experiment. During the test, the chosen sample was put on the working plane of the electronic tester and compressed along the X-axis direction. The pressing speed was maintained at 0.1 mm/s. The loading time was set 10 s. The elastic modulus of the ginkgo nuts was calculated using the Hertz formula,
where E is the elastic modulus (MPa), F is the compression force (N), D is the dimensional change of the ginkgo nut (mm), μ is the Poisson’s ratio, RU is the minimum radius of curvature of the convex surface of the sample at the point of contact with the upper plate (mm), RU’ is the maximum radius of curvature of the convex surface of the sample at the point of contact with the upper plate (mm), RL is the minimum radius of curvature of the convex surface of the sample at the point of contact with the lower plate (mm), RL’ is the maximum radius of curvature of the convex surface of the sample at the point of contact with the lower plate (mm), KU is a constant related to principal curvature radius, KL is a constant related to principal curvature radius, and θ is the angle between the contact surface of the ginkgo nut and the contact surface of the probe (°).
KU and KL can be retrieved from the ASAE S368.4 (2017) standard by calculating cos θ. As the radius of curvature of the ginkgo nut is almost the same as that of the upper and lower contact surfaces, Eq. 2 can be simplified as
The mean elastic modulus of the ginkgo nut was 18.92 MPa and the standard deviation was 1.91. Therefore, the shear modulus was obtained as 7.09 MPa using Eq. 5,
(5)
where G is the shear modulus (MPa).
Measurement of Contact Parameters
Measurement of restitution coefficient
The restitution coefficient is the ratio of the post-collision rebound normal velocity to the pre-collision normal velocity. The oblique collision method is usually used to measure the restitution coefficient, which is calculated by manually measuring the drop and rebound distances of particles (Wang et al. 2021). The conventional method has problems such as large measurement errors (Wang et al. 2021). The method seems to be mainly suitable for materials with relatively regular shape and uniform composition. Ginkgo kernels and shells are highly variable in composition and can move relative to each other. In order to minimize the error, a camera was used to record the collision process of a vertically falling ginkgo nut and a horizontal contact material, and the normal velocities before and after the collision were calculated. Then, the restitution coefficient was calculated. This method can effectively reduce the measurement error.
The test principle of the restitution coefficient is shown in Fig. 2. The ginkgo nut falls freely from the placement hole of a drop frame. The vertical displacement and corresponding instant time before the collision with the contact material are recorded with the aid of a background plate and a high-speed camera. Then, the pre-collision normal velocity (v0) is calculated. The same method is used to record the vertical displacement and corresponding time of the ginkgo nut before it reaches the highest point after collision with the contact material. The next step is to calculate the post-collision normal velocity (vz), which is the velocity of the post-collision velocity (v) in the z-direction. The restitution coefficient can be calculated using Eq. 6,
(6)
where e is the restitution coefficient; vz is the post-collision normal velocity of ginkgo nut (m/s); and v0 is the pre-collision normal velocity of ginkgo nut (m/s).
Fig. 2. Test principle of the restitution coefficient.1: Ruler; 2: Drop frame; 3: Contact material; 4: Ginkgo nut; 5:Background plate; 6: Placement hole
Note: v is the post-collision velocity; vx, vy, and vz are the post-collision velocity in x-direction, y-direction, and z-direction, respectively.
According to the above principles, the restitution coefficients of the ginkgo nut–steel plate and ginkgo nut–ginkgo nut were tested. The measurement setup of the test is shown in Fig. 3, which mainly consists of a display, a drop frame, a background plate, a light source, a high-speed camera, and a base plate, etc. During the test, the contact material to be tested was placed on the base plate. The ginkgo nut was released from the placement hole with an initial velocity of 0 m/s and rebounded after colliding with the contact material. Videos and photographs of the ginkgo nut collision were captured using a high-speed camera (Fastec Imaging, San Diego, CA, USA). The camera was set at 800 f/s. Therefore, the time interval between each photo was 1.25 ms. Coordinate paper was taped to the background plate, and the falling distance could be obtained through the background plate. The pre-collision normal velocity (v0) was calculated by measuring the falling distance of the particle before collision within 5 ms. Similarly, the post-collision normal velocity (vz) was calculated by measuring the rise distance of the particle after collision within 5 ms. González-Montellano et al. (2012) utilized a similar method to measure restitution coefficients, and in order to avoid the particles from being subjected to other forces or rotating, the particles were released by changing the vacuum level. In this study, the ginkgo nut fell through the placement hole and was subjected only to gravity, avoiding other forces. Photographs of the bottom surface of the particle colliding with the contact material were selected. When the particles rotated so that the protruding prongs collided with the contact material, the direction of motion changed significantly, so these photographs were excluded. The restitution coefficient is an intrinsic property of granular materials and does not change due to the change in the drop height of the material. In this study, the height of the drop frame was set to 200 mm.
Fig. 3. Drop impact test.1: Display; 2: Drop frame; 3: Background plate; 4: Light source; 5: Base plate; 6: High-speed camera
The restitution coefficient was tested 20 times, and the calculation results of Eq. 6 were averaged. The restitution coefficients of the ginkgo nut–steel plate was 0.32 and the standard deviation was 0.022. The restitution coefficients of the ginkgo nut–ginkgo nut was 0.38 and the standard deviation was 0.034.
Measurement of static friction coefficient
In this study, the slope method was used to determine the coefficient between a ginkgo nut and different materials according to ASTM D4918-97(2002). The principle of measurement is presented in Fig. 4.
Fig. 4. Measurement principle of the static friction coefficient
When the ginkgo nut is in the critical state of downward sliding in the X direction,
(7)
where F is the component force of gravity on the ginkgo nut in the X direction, N; and f is the static friction between the two materials (N).
F and f were determined from Eqs. 8 and 9, respectively,
(8)
(9)
where G is the force of gravity (N); α is the angle of static friction, °; γ is the coefficient of static friction; and N is the component force of gravity (N).
In summary, Eq. 10 can be obtained.
When the ginkgo nut is in the critical state of downward sliding, the slope angle is the angle of static friction, and the coefficient of static friction can be obtained using Eq. 10. The dynamic friction coefficient is usually smaller than the static friction coefficient. When the particles are placed on an inclined plane and the angle of inclination is slowly raised, the particles roll before they begin to slide, making the static friction coefficient impossible to measure. Therefore, sticking ginkgo nuts together can avoid rolling and facilitate the testing of static friction coefficient. When testing the static friction coefficient of the ginkgo nut–steel plate, four ginkgo nuts were stuck and placed on the self-made inclinometer to prevent rolling. The inclinometer was slowly increased by lifting the handle until the ginkgo nuts exhibited a downward sliding trend. The angle of static friction was tested by an angle ruler (Deqing Shengtaixin Electronic Technology Co., Ltd, Huzhou, China). When measuring the static friction coefficient of ginkgo nut–ginkgo nut, the ginkgo nuts were arranged to form a material plate and pasted on the inclinometer. Then, four ginkgo nuts were stuck and placed on the material plate of ginkgo nuts. The test is presented in Fig. 5. The test was repeated 20 times, and the calculated results of Eq. 10 were averaged. The static friction coefficient of the ginkgo nut–steel plate was 0.58 and the standard deviation was 0.028. The static friction coefficient of the ginkgo nut–ginkgo nut was 0.73 and the standard deviation was 0.048.
Fig. 5. Static friction coefficient test.1: Inclinometer; 2: Ginkgo nut; 3: Angle ruler
Measurement of dynamic friction coefficient
The dynamic friction coefficient, also known as the rolling friction coefficient, is tested in different ways. Wójcik et al. (2020) conducted tests on three plant particulate materials (peas, maize, triticale) to obtain a direct relationship between maximum static friction and mean dynamic friction. Khazaei et al. (2010) proposed a new method for simultaneous measurement of angle of repose and coefficient of friction of wheat grains based on the motion of the grains in a partially-filled horizontal rotating drum at the slumping and rolling modes of motion. Kweon et al. (2007) developed a device to measure the dynamic friction coefficient of particles accelerated in a tube rotating in a horizontal plane. Combarros et al. (2014) calibrated the dynamic friction coefficient by comparing the macroscopic behavior of granular matter in standard bulk experiments, which have the advantage of being highly repeatable and reliable. Liu et al. (2018a) calculated the dynamic friction coefficient by testing the rolling distance of the particles on inclined and horizontal planes. In this study, calculating the dynamic friction coefficient by measuring the angle of inclination of the inclined plane when the particles start rolling on the inclined plane is another common method, which can also be referred to relevant research (Hou et al. 2020). The method used to measure the dynamic friction coefficient can refer to that used to measure the static friction coefficient. When measuring the dynamic friction coefficient of ginkgo nut–steel plate, a ginkgo nut was placed on the steel inclinometer in the thickness direction. The inclined plane was raised slowly until the ginkgo nut rolled downward; the inclination angle of the inclined plane was the dynamic friction angle. The test is presented in Fig. 6. When measuring the dynamic friction coefficient of ginkgo nut–ginkgo nut, the ginkgo nuts were arranged to form a material plate and pasted on the inclinometer. Then, a ginkgo nut was placed on the material plate of ginkgo nuts. The results were averaged after 20 repetitions of the test. The dynamic friction coefficients of the ginkgo nut–steel plate was 0.46, and the standard deviation was 0.045. The mean dynamic friction coefficient of the ginkgo nut–ginkgo nut was 0.53, and the standard deviation was 0.061.
Fig. 6. Dynamic friction coefficient test
Calibration of Simulation Parameters
Many scholars have carried out research on stacking angle testing methods. Du et al. (2023) measured the stacking angle of granular materials in DEM software based on Python’s linear regression analysis and particle coordinate acquisition algorithm, and the consistency and stability of the measurement results were good. Madrid et al. (2022) calculated the stacking angle by measuring the height of the granular pile directly. In this study, the cylinder lifting method was chosen to measure the stacking angle. The height of the selected cylinder was 150 mm, and the diameter was 80 mm. During the test, the cylinder was placed on a steel plate chosen and filled with ginkgo nuts to be tested. It was lifted at a uniform speed of 0.05 m/s so that the ginkgo nuts formed a pile. A camera was used to take vertical photos of the ginkgo nut pile, and the photos were imported into computer aided design (CAD) software (Autodesk Inc., San Rafael, CA, USA) to measure the angle (Dai et al. 2019; Liao et al. 2020). The results were averaged after 20 repetitions of the test. The stacking angle measured was 39.59° and the standard deviation was 0.535.
Fig. 7. Stacking angle test of ginkgo nuts
Simulation model establishment and simulation test
In this study, EDEM software (DEM Solutions Ltd., Edinburgh, UK) was used to carry out the simulation tests. Since the ginkgo nuts are spindle-shaped, if a single regular spherical particle model is used to build the simulation model, the characteristics of the material and its interaction will not be accurately reflected. Therefore, this study used the method of aggregating multiple spherical particles to build a model of the ginkgo nut to approximate its actual shape (Kruggel-Emden et al. 2008; Markauskas et al. 2010). The simulation model of the ginkgo nut is shown in Fig. 8. The coordinates and radius of each sphere are shown in Table 3.
Fig. 8. Simulation model of a ginkgo nut
Table 3. Coordinates and Radius of Each Sphere of the Ginkgo Nut Model
The liquid bridge force and adhesion force between ginkgo nuts can be ignored, and the ginkgo nut can be regarded as a non-viscous body. Considering the calculation accuracy and simulation speed, the Hertz–Mindlin non-slip model was chosen. Based on the establishment of the model of the ginkgo nut, the stacking angle test model was established according to the stacking angle measuring device, and the lifting process of the cylinder was simulated. Then, a virtual face was built above the cylinder as a particle factory. The ginkgo nuts were dynamically generated in the particle factory, and the initial velocity was 0.5 m/s downward along the cylinder center. The size distribution of ginkgo nuts is normal, with a generation rate of 200 grains/s and a generation time of 2.5 s. A total of 500 ginkgo nut granules were generated. After the particles were stable, the cylinder moved uniformly upward along the center of the cylinder, so that the ginkgo nuts naturally formed a pile. The determination model of stacking angle is shown in Fig. 9.
Fig. 9. Determination model of stacking angle
RESULTS AND DISCUSSION
Plackett–Burman test
The Plackett–Burman test design was developed using the Design-Expert software (Stas-Ease Inc., Minneapolis, MN, USA). The experimental parameters and ranges were selected based on the physical parameters measured in the early stages and related literature (Zhang et al. 2022b). The stacking angle of the ginkgo nuts was considered as the research object, and the factors that had an essential impact on the stacking angle were selected. The test factors and results are presented in the Tables below.
Table 4. Parameters of Plackett-Burman Test
Table 5. Factor Levels and Plackett-Burman Test Results
The variance of the obtained results was analyzed using the Design-Expert software. The analysis results are presented in Table 6. The p-value corresponding to each factor can reflect the impact of the factor on the stacking angle. When p < 0.01, it indicates that the effect of the factor on the stacking angle is very significant. When 0.01 < p < 0.05, it indicates that the effect of the factor on the stacking angle is significant. When p > 0.05, it indicates that the effect of the factor on stacking angle is not significant.
Table 6. Analysis of Factors in Plackett-Burman Test
According to Table 6, the static friction coefficient of ginkgo nut–ginkgo nut X4 and the dynamic friction coefficient of ginkgo nut–ginkgo nut X5 significantly affected the stacking angle of the ginkgo nuts (P < 0.01). The restitution coefficient of ginkgo nut–steel plate X6 had an important effect on the stacking angle of the ginkgo nuts (P < 0.05). The effects of the remaining parameters on the response value were not significant (P > 0.05).
Steepest Climbing Test
The factors that have a significant impact on the stacking angle were determined using the Plackett–Burman test. These factors were studied using the steepest climbing test to obtain a range of values for each factor. The steepest climbing test is an optimization algorithm experiment that simulates the process of finding the steepest climbing in a curvilinear terrain, and gradually approximates the optimal solution by continuously adjusting the values of the parameters. In the simulation test, the average value was selected for other non-significant factors, that is, the Poisson’s ratio of ginkgo nut is 0.3, the restitution coefficient of ginkgo nut–ginkgo nut is 0.35, the shear modulus of ginkgo nut is 7.5 MPa, the static friction coefficient of ginkgo nut–steel plate is 0.55, and the dynamic friction coefficient of ginkgo nut–steel plate is 0.5. The steepest climbing test scheme and results are presented in Table 7.
Table 7. Steepest Climbing Experimental Scheme and Results
The results indicate that the relative error of the stacking angle was the smallest when the No. 2 test combination was used. It means that the simulated stacking angle obtained using No. 2 test combination was closest to the actual value, which ensures the accuracy of the simulation test. The relative errors of combinations No. 1 and No. 3 were closest to those of No. 2. Therefore, No. 2 was used as the intermediate value, and Nos. 1 and 3 were used as the low level and high level of the subsequent response surface test, respectively.
Box–Behnken Test
The Plackett-Burman test determined that the number of factors that significantly affect the stacking angle was three. They were the static friction coefficient of ginkgo nut–ginkgo nut X4, the dynamic friction coefficient of ginkgo nut–ginkgo nut X5, and the restitution coefficient of ginkgo nut–steel plate X6. The steepest climbing test determined the intermediate value, low level, and high level of the significant influencing factors, so the level of each factor was three. Based on the previous test, the static friction coefficient of ginkgo nut–ginkgo nut X4, the dynamic friction coefficient of ginkgo nut–ginkgo nut X5, and the restitution coefficient of ginkgo nut–steel plate X6. were considered as the factors of test. The stacking angle was used as an assessment index. A response surface test analysis was developed using the three-factor three-level quadratic regression test design scheme. The factors and levels of the response surface tests are shown in Table 8. The remaining factors were consistent with the previously performed climbing test values.
Table 8. Factors and Levels of Response Surface Test
The Box–Behnken test was carried out using the Design-Expert software. The total number of tests was 17, including 12 analysis tests and 5 estimation error tests. The design and results of the experiment are presented in Table 9.
Table 9. Experimental Design and Response Values
According to the experimental results, a multiple regression fitting analysis was developed using the Design-Expert software. A quadratic polynomial regression model of the stacking angle θ on three independent variables (X4, X5, X6) was obtained, as shown in Eq. 11. Then, a variance analysis of the obtained regression equation was carried out. The results obtained are presented in Table 10.
(11)
As shown in Table 10, the P value of the stacking angle response surface model was less than 0.01, indicating that the regression model was very significant. The P value of the misfit term was greater than 0.05, indicating a good fit of the regression equation. The model coefficient of determination R2 was 0.9648, indicating that this model could explain more than 96% of the evaluation indicators, and less than 4% of the total variation could not be explained by the model, indicating a high correlation between the predicted value of the model and the actual value and a small experimental error. Therefore, it is feasible to predict the stacking angle θ of ginkgo nuts using the model.
Table 10. Variance Analysis of the Regression Equation
Table 10 shows that the P value of X4 was less than 0.01, indicating that the effect of this factor on the model was highly significant. The P values of X5, X6, X42, and E2 were all less than 0.05, and these factors significantly influenced on the model. The effects of X4X5, X4X6, X5X6, and X62 were not important. Different studies have shown that dynamic and static friction coefficients between the materials have extremely significant effect on the stacking angle (Fu et al. 2023; Lei et al. 2023), which is consistent with the results of this study. In addition, the friction coefficient and restitution coefficient between the particles and contact materials also have a significant effect on the stacking angle in different studies, which is related to the different particles as well as contact materials. By analyzing the coefficients of each regression term in the regression model, it was found that the strength of the influence of different factors on the stacking angle of ginkgo nuts was: X4>X5>X6. The significant terms of the above model were retained, and the insignificant terms were eliminated to ensure that the model was highly significant and the lack-of-fit terms were not significant. The model was re-fitted to obtain the optimized regression model:
(12)
Verification Test
The effects of the various factors are different. Regression Eq. 12 needs to be optimized so that the simulated value of the stacking angle is as close as possible to the experimental value. The optimized parameter combination was captured: the static friction coefficient of ginkgo nut–ginkgo nut, the dynamic friction coefficient of ginkgo nut–ginkgo nut, and the restitution coefficient of ginkgo nut–steel plate were 0.50, 0.36, and 0.42, respectively. Because the optimized optimal combination was not included in the response surface test, a DEM simulation test was developed using the combination determined above to test the reliability of the predicted value and the rationality of the test results. The results were averaged after five repetitions of the test. The simulated stacking angle was 41.12°. The standard deviation was 1.79°, and the relative error between the simulated and actual values was 3.86%. The comparison shows that the calibrated parameters are accurate and can be used for subsequent discrete element simulations test of ginkgo nut shelling.
CONCLUSIONS
In this study, three factors that had an essential impact on the stacking angle of the ginkgo nut were captured: the static friction coefficient of ginkgo nut–ginkgo nut, the dynamic friction coefficient of ginkgo nut–ginkgo nut, and the restitution coefficient of ginkgo nut–steel plate. A Box–Behnken experiment was conducted to obtain a quadratic polynomial regression model of the three factors on the stacking angle. The model was optimized, and the experimental value of 39.59° was used as the target. The optimal parameter combination was obtained as follows: the static friction coefficient of ginkgo nut–ginkgo nut, the dynamic friction coefficient of ginkgo nut–ginkgo nut, and the restitution coefficient of ginkgo nut–steel plate were 0.50, 0.36, and 0.42, respectively. The experiment was conducted using an optimal parameter combination. The simulated stacking angle obtained was 41.12°. The relative error between the simulated and actual values was 3.86%. This indicates that the model is reliable, and that the calibrated parameters can provide a useful reference for relevant research.
ACKNOWLEDGMENTS
The authors are grateful for the support of primary processing equipment for major grain cash crops of the innovation project of the Chinese Academy of Agricultural Sciences, Grant No. 31-NIAM-09.
REFERENCES CITED
Andersson, S., Wang, Y., Pönni, R., Hänninen, T., Mononen, M Ren, H., Serimaa, R., and Saranpää, P. (2015). “Cellulose structure and lignin distribution in normal and compression wood of the Maidenhair tree (Ginkgo biloba L.),” Journal of Integrative Plant Biology 57(4), 388-395. DOI: 10.1111/jipb.12349
ASAE S368.4 DEC2000 (R2017). “Compression test of food materials of convex shape,” American Society of Agricultural and Biological Engineers, Michigan, USA.
ASTM D4918-97(2002). “Test method for coefficient of static friction of uncoated writing and printing paper by use of the inclined plane method,” ASTM International, West Conshohocken, PA, USA.
Ch’ng, P. E., Abdullah, M. H. R. O., Mathai, E. J., and Yunus, N. A. (2013). “Some physical properties of gingko nuts and kernels,” International Agrophysics 27(4), 485-489. DOI: 10.2478/intag-2013-0019
Coetzee, C. J., Els, D. N. J., and Dymond, G. F. (2010). “Discrete element parameter calibration and the modelling of dragline bucket filling,” Journal of Terramechanics 47, 33-44. DOI: 10.1016/j.jterra.2009.03.003.
Coetzee, C. J. (2017). “Calibration of the discrete element method,” Powder Technology 310, 104-142.
Combarros, M., Feise, H. J., Zetzener, H., and Kwade, A. (2014). “Segregation of particulate solids: Experiments and DEM simulations,” Particuology 12, 25-32. DOI: 10.1016/j.partic.2013.04.005
Crouvisier-Urion, K., Lagorce-Tachon, A., Lauquin, C., Winckler, P., Tongdeesoontorn, W., Domenek, S., Debeaufort, F., and Karbowiak, T. (2017). “Impact of the homogenization process on the structure and antioxidant properties of chitosan-lignin composite films,” Food Chemistry 236, 120-126. DOI: 10.1016/j.foodchem.2017.03.094
Dai, F., Zhao, W. Y., Liu, G. C., Zhang, S. L., Shi, R. J., and Wei, B. (2019). “Design and experiment of separating and cleaning machine for flax threshing material,” Transactions of the Chinese Society for Agricultural Machinery 50, 140-147. DOI: 10.6041/j.issn.1000-1298.2019.08.016.
Du, X., Liu, C. L., Liu, C. Q., Sun, Q. X., and Chen, S. F. (2023). “A novel method for measurement of the angle of repose of granular seeds in discrete element methods,” Journal of Agricultural Engineering 54(2). DOI: 10.4081/jae.2023.1504
Fu, J. J., Cui, Z. C., Chen, Y. S., Guan, C. S., Chen, M. J., and Ma, B. (2023). “Simulation and experiment of compression molding behavior of substate block suitable for mechanical transplanting based on discrete element method (DEM),” Agriculture 13, article 883. DOI: 10.3390/agriculture13040883.
GB/T 22315 (2008). “Metallic materials – Determination of modulus of elasticity and Poisson’s ratio,” Standardization Administration of China, Beijing, China.
GB/T 3543.6 (1995). “Rules for agricultural seed testing – Determination of moisture content,” Standardization Administration of China, Beijing, China.
GB/T 5519 (2018). “Cereals and pulses – Determination of the mass of 1000 grains,” Standardization Administration of China, Beijing, China.
González-Montellano, C., Fuentes, J. M., Ayuga-Téllez, E., and Ayuga, F. (2012). “Determination of the mechanical properties of maize grains and olives required for use in DEM simulations,” Journal of Food Engineering 111(4), 553-562. DOI: 10.1016/j.jfoodeng.2012.03.017
Hong, W., Wang, X., Zheng, H.Y.; Li, R., Wu, R., and Chen, J.S. (2021). “Molten-salt-assisted synthesis of nitrogen-doped carbon nanosheets derived from biomass waste of ginkgo shells as efficient catalyst for oxygen reduction reaction,” Processes 9. DOI: 10.3390/PR9122124
Horabik, J., and Molenda, M. (2016). “Parameters and contact models for DEM simulations of agricultural granular materials: A review,” Biosystems Engineering, 147, 206-225.
Hou, Z. F., Dai, N. Z., Chen, Z., Qiu, Y., and Zhang, X. W. (2020). “Measurement and calibration of physical property parameters for agropyron seeds in a discrete element simulation,” Transactions of the Chinese Society of Agricultural Engineering 36, 46–54. DOI: 10.11975/j.issn.1002-6819.2020.24.006
Jiang, Bo., Zhang, Y., Guo, T. Y., Zhao, H. F., and Jin, Y. C. (2018). “Structural characterization of lignin and lignin-carbohydrate complex (LCC) from ginkgo shells (Ginkgo biloba L.) by comprehensive NMR spectroscopy,” Polymers 10(7). DOI: 10.3390/polym10070736.
Khazaei, J., and Ghanbari, S. (2010). “New method for simultaneously measuring the angles of repose and frictional properties of wheat grains,” International Agrophysics 24(3), 275-286.
Kruggel-Emden, H., Rickelt, S., Wirtz, S., and Scherer, V. (2008). “A study on the validity of the multi-sphere discrete element method,” Powder Technology 188, 153-165.
Kweon, G., Grift, T. E., and Miclet, D. (2007). “A spinning-tube device for dynamic friction coefficient measurement of granular fertiliser particles,” Biosystems Engineering 97(2), 145-152. DOI: 10.1016/j.biosystemseng.2007.02.013
Lei, X.L., Wu, W. C., Deng, X., Li, T., Liu, H. N., Guo, J. Y., Li, J., Zhu, P. X., and Yang, K. (2023). “Determination of material and interaction properties of granular fertilizer particles using DEM simulation and bench testing,” Agriculture 13(9), article 1704. DOI: 10.3390/agriculture13091704
Li, Y. X., Li, F. X., Xu, X. M., Shen, C. P., Meng, K. P., Chen, J., and Chang, D. T. (2019). “Parameter calibration of wheat flour for discrete element method simulation based on particle scaling,” Transactions of the Chinese Society of Agricultural Engineering 35, 320-327. DOI: 10.11975/j.issn.1002-6819.2019.16.035
Liao, Y. T., Wang, Z. T., Liao, Q. X., Wang, X. Y., Zhou, Y., and Liang, F. (2020). “Calibration of discrete element model parameters of forage rape stalk at early pod stage,” Transactions of the Chinese Society for Agricultural Machinery 51(s1), 236-243. DOI: 10.6041/j.issn.1000-1298.2020.S1.027
Liu, F. Y., Zhang, J., and Chen, J. (2018a). “Construction of visco-elasto-plasticity contact model of vibratory screening and its parameters calibration for wheat,” Transactions of the Chinese Society of Agricultural Engineering 34, 37-43. DOI: 10.11975/j.issn.1002-6819.2018.15.005
Liu, M. J., Wang, J. N., Wang, N., Xie, H. X., Zhang, H. J., Shen, H. Y., and Gao, X. M. (2022). “Design and optimization of sheller for ginkgo nut: A study about multifunctional ginkgo nut,” PLoS ONE 17(10), article e0276139. DOI: 10.1371/journal.pone.0276139.
Liu, W. Z., He, J., Li, H. W., Li, X .Q., Zheng, K., and Wei, Z. C. (2018b). “Calibration of simulation parameters for potato minituber based on EDEM,” Transactions of the Chinese Society for Agricultural Machinery 49, 125-135, 142. DOI: 10.6041/j.issn.1000-1298.2018.05.014
Madrid, M., Fuentes, J., Ayuga, F., and Gallego, E. (2022). “Determination of the angle of repose and coefficient of rolling friction for wood pellets,” Agronomy 12, article 424. DOI: 10.3390/agronomy12020424
Markauskas, D., Kačianauskas, R., Džiugys, A., and Navakas, R. (2010). “Investigation of adequacy of multi-sphere approximation of elliptical particles for DEM simulations,” Granular matter 12, 107-123. DOI: 10.1007/s10035-009-0158-y
Ni, Y., Li, J., and Fan, L. (2020). “Production of nanocellulose with different length from ginkgo seed shells and applications for oil in water Pickering emulsions,” International Journal of Biological Macromolecules 149, 617-626. DOI: 10.1016/j.ijbiomac.2020.01.263
Nowak, A., Kojder, K., Zielonka-Brzezicka, J., Wróbel, J., Bosiacki, M., Fabiańska, M., Wróbel, M., Sołek-Pastuszka, J., and Klimowicz, A. (2021). “The use of Ginkgo biloba L. as a neuroprotective agent in Alzheimer’s disease,” Frontiers in Pharmacology 12, article 775034. DOI: 10.3389/fphar.2021.775034
Qin, L., Feng, L. H., Li, C., Fan, Z., Zhang, G. L., Shen, C., and Meng, Q. (2019). “Amination/oxidization dual-modification of waste ginkgo shells as bio-adsorbents for copper ion removal,” Journal of Cleaner Production 228, 112-123. DOI:10.1016/j.jclepro. 2019.4.249.
Santos, K. G., Campos, A. V. P., Oliveira, O. S., Ferreira, L. V., Francisquetti, M. C., and Barrozo, M. A. S. (2015). “DEM simulations of dynamic angle of repose of acerola residue: A parametric study using a response surface technique,” Blucher Chemical Engineering Proceedings 1, 11326-11333. DOI: 10.5151/chemeng-cobeq2014-0187-26615-184001
Shen, N., Zeng, W., Leng, F., Lu, J. K., Lu, Z. G., Cui, J., Wang, L., and Jin, B. (2021). “Ginkgo seed extract promotes longevity and stress resistance of Caenorhabditis elegans,” Food & Function 12(24), 12395-12406. DOI: 10.1039/d1fo02823e
Usuki, T., Yasuda, N., Yoshizawa-Fujita, M., and Rikukawa, M. (2011). “Extraction and isolation of shikimic acid from Ginkgo biloba leaves utilizing an ionic liquid that dissolves cellulose,” Chemical Communications (Cambridge, England) 47(38), 10560-10562. DOI: 10.1039/c1cc13306c
Wang, L. J., Liu, T. H., Feng, X., Gao, Y. P., Wang, B., and Zhang, S. (2021). “Research progress of the restitution coefficients of collision of particles in agricultural and food fields,” Transactions of the Chinese Society of Agricultural Engineering 37, 313-322. DOI: 10.11975/j.issn.1002-6819.2021.20.035
Wen, X. Y., Yuan, H. F., Wang, G., and Jia, H. L. (2020). “Calibration method of friction coefficient of granular fertilizer by discrete element simulation,” Transactions of the Chinese Society for Agricultural Machinery 51, 115-122, article 142. DOI: 10.6041/j.issn.1000-1298.2020.02.013
Wójcik, A., Frączek, J., and Niemczewska-Wójcik, M. (2020). “The relationship between static and kinetic friction for plant granular materials,” Powder Technology 361, 739-747.
Yu, J. Q., Fu, H., Li, H., and Shen, Y. F. (2005). “Application of discrete element method to research and design of working parts of agricultural machines,” Transactions of the Chinese Society of Agricultural Engineering 21, 1-6. DOI: 10.3321/j.issn:1002-6819.2005.05.001
Zhang, Q. L, Geng, T., and Zhang, L. (2022a). “Development situation and prospects of ginkgo food industry,” Food Industry 43, 208-212.
Zhang, S. W., Zhang, R. Y., Chen, T. Y., Fu, J., and Yuan, H. F. (2022b). “Calibration of simulation parameters of mung bean seeds using discrete element method and verification of seed-metering test,” Transactions of the Chinese Society for Agricultural Machinery 53, 71-79. DOI: 10.6041/j.issn.1000-1298.2022.03.007
Article submitted: June 3, 2024; Peer review completed: July 31, 2024; Revised version received: August 20, 2024; Accepted: August 25, 2024; Published: September 10, 2024.
DOI: 10.15376/biores.19.4.8154-8172