Next Article in Journal
Effect of Cymbopogon citratus on Enteric Methane Emission, Nutrients Digestibility, and Energy Partition in Growing Beef Cattle
Previous Article in Journal
A Review on Soil Nitrogen Sensing Technologies: Challenges, Progress and Perspectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Experimental and Numerical Study for Discrete Element Model Parameters Calibration: Gluten Pellets

1
College of Engineering, Nanjing Agricultural University, Nanjing 210031, China
2
College of Biotechnology and Food Engineering, Chuzhou University, Chuzhou 239000, China
3
Anhui Bi Lv Chun Biotechnology Co., Ltd., Chuzhou 239200, China
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(4), 744; https://doi.org/10.3390/agriculture13040744
Submission received: 24 February 2023 / Revised: 21 March 2023 / Accepted: 22 March 2023 / Published: 23 March 2023
(This article belongs to the Section Agricultural Technology)

Abstract

:
Discrete element method (DEM) simulation is widely used to calculate the flow characteristics of particles under certain conditions. DEM input parameters are the prerequisite for the accurate modeling and simulation of particles. In order to explore the mechanical properties and breaking behavior of gluten pellets, the pellet material property, the interaction parameters of pellet–stainless steel and pellet–pellet (multi-spheres autofill model), and the bonding parameters (bonded particle model) were calibrated by experiments and simulations. The relative error of the angle of repose, the breaking displacement, and the breaking force between simulated and experimental values were 0.28%, 0.66%, and 1.09%, respectively. Based on the regression analysis in the Design-Expert 12.0 software, the relationships among evaluating indicators (angle of repose, breaking displacement, and breaking force) and their corresponding influencing factors were established, respectively. Meanwhile, the feasibility of applying the interaction parameters of the multi-spheres autofill model to the bonded particle model was verified through the free fall test, the inclined plane sliding test, and the inclined plane tumbling time test. This work can provide a reference for the design of pellet feed processing and transportation machinery.

1. Introduction

Gluten contains fifteen amino acids that the human body needs and it is a natural protein isolated from wheat flour that meets the nutritional needs of modern diets [1,2]. It is widely used in bread [3], noodles [4], instant noodles [5], meat product additives [6], beverage additives [7], and processing high-grade aquatic feed [8]. After extrusion molding, gluten pellets undergo conveying, drying, and packaging processes. When the gluten pellets move in the equipment mentioned above, they collide and squeeze with the inner wall of the equipment and other pellets, resulting in serious problems such as pellet breakage and pulverization, thus reducing the economic value of the product [9,10,11]. Therefore, studying the mechanical interaction parameters of gluten pellets is of great significance.
Cundall proposed the discrete element method (DEM) in 1979 to simulate and analyze particle behavior by establishing a parametric model of a solid particle system [12]. In recent years, commercial software based on the DEM has been widely used in agricultural engineering. The properties of gluten pellets were affected by factors such as hygroscopicity, composition, particle size, etc., so the parameter calibration must be conducted according to specific conditions [13,14,15]. Scholars at home and abroad have conducted much work on the parameter calibration of agricultural materials. The fundamental physical property and mechanical interaction parameters were usually obtained through experiments and simulations. Niu and Peng [10,16] calibrated the discrete element parameters of pellet feed by measuring the angle of repose and used EDEM simulation software to research the calibration of discrete element simulation parameters of pellet feed uniaxial compression damage. Zhao et al. [17] calibrated the discrete element parameters of coconut bran particles and verified the coconut bran’s stacking angle and fluidity by experiments. Kong et al. [18] combined computational fluid dynamics (CFDs) and the DEM to study the effects of the inlet velocity and bending radius on pellet breakage rates and energy loss during the pneumatic conveying of aquafeed pellets. Wang et al. [19] measured the angle of repose of alfalfa with different moisture contents and obtained the interaction parameters by combining it with a discrete element simulation. Li et al. [20] proposed a corn kernel model-building method for Rocky DEM simulations and a calibration method for calibrating the interaction parameters. There are few reports on the research on the discrete element model of gluten pellets both at home and abroad [21].
Pelletization is a polymerization process of raw powder materials by mechanical extrusion, which is most widely used in the feed and biomass energy industry [22,23]. In the process of pelletization, friction between raw materials and mold generates a large amount of heat, which causes the gelatinization of starch and the denaturation of protein to become a binder and bond raw material particles together [24]. At the same time, high temperature and high pressure in the pelletization process can eliminate salmonella and reduce microbial activity, thus improving feed safety [25]. Extrusion molding can significantly improve the bulk density of raw materials, thereby improving their transportation and storage performance and reducing transportation and storage costs [26].
This study takes gluten pellets as the research object, and their physical parameters (size, density, and shear modulus), interaction parameters between pellets and stainless steel plates (static friction coefficient, rolling friction coefficient, and restitution coefficient), and angle of repose (image method) were measured. Based on the Hertz–Mindlin contact model of EDEM software, the multi-spheres autofill model and bonded particle model were established. The contact and bonded parameters between pellets were calibrated by the angle of repose and compression simulation, combined with the response surface optimization test. Finally, the feasibility of using the interaction parameters of the multi-spheres autofill model to the bonded particle model was verified.

2. Materials and Methods

2.1. Material Properties and DEM Model of Gluten Pellets

2.1.1. Material Properties Determination

The gluten pellets were provided by Bilvchun Biotechnology Co., Ltd. And produced in Lai’an County, Chuzhou City, Anhui Province, China. The gluten pellets are pale yellow with a protein content of 85%, and the initial moisture content (dry-basis) is 10.18 ± 0.10% [27].
  • Size and density
The electronic digital vernier caliper (DL91150, Shanghai Yuanmai Trading Co., Ltd., Shanghai, China) was used to measure the length and height. The electronic balance (ME55/02, METTLER TOLEDO Instruments (Shanghai) Co., Ltd., Shanghai, China) was used to measure the mass of pellets. The density was calculated by Equation (1).
ρ ¯ = m ¯ π 4 d ¯ 2 L ¯ × 1000
where ρ ¯ is the average density of pellets, 103 kg·m−3; m ¯ is the average mass of pellets, kg; d ¯ is the average height of pellets, m; and L ¯ is the average length of pellets, mm.
Fifty pellets were randomly selected, and the final result was expressed as means ± standard deviation. The mass, height, length, and density of gluten pellet are 0.178 × 10−3 ± 0.041 × 10−3 kg, 4.24 × 10−3 ± 0.13 × 10−3 m, 11.43 × 10−3 ± 2.43 × 10−3 m, and 1099 ± 67 kg·m−3, respectively.
  • Shear modulus
Gluten pellets are brittle materials, and Poisson’s ratio can be obtained from reference [16] as 0.4. A texture analyzer (TA.XT.plus, Stable Micro System, Surrey, UK) with a P/36R probe was used to carry out a compression test on the gluten pellets until the pellets were broken. Gluten pellets were placed horizontally on a rigid base. The test speed and the trigger force were 0.5 × 10−3 m·s−1 and 0.049 N, respectively. The texture compression experiment device is shown in Figure 1. The experiment was repeated five times, and the average value was taken. The force–displacement curve is shown in Figure 2. The average compression displacement of gluten pellets was 1.209 × 10−3 m, and the average ultimate breaking force was 94.656 N. The slope of the stress–strain curve, which was automatically drawn by the computer, was elastic modulus. The shear modulus was calculated by Equation (2).
G = E 2 ( 1 + ν )
where G is shear modulus, Pa; ν is Poisson’s ratio; and E is elastic modulus, Pa.
The gluten pellets’ elastic modulus and shear modulus can be obtained as 2.207 × 108 ± 4.771 × 107 Pa and 7.881 × 107 ± 1.704 × 107 Pa, respectively.

2.1.2. DEM Model

(1)
Gluten pellet model
In this study, a multi-spheres autofill model and a bonded particle model were established and used for the angle of repose simulation and bonding parameters calibration, respectively. Considering factors such as shape fitting, number of particles, or computation time, the two modeling methods are as follows.
  • Multi-spheres autofill model: According to the average geometric size of gluten pellets measured by the test above, a three-dimensional model was established in SolidWorks software, and the file (“.stl” format) was exported. Then, the file was imported into EDEM as an autofill template. The smoothing value and the minimum sphere radius were set to 2 and 0.5 × 10−3 m, respectively. The generated pellet model with a total of 177 particles is shown in Figure 3a. In order to be as close to the actual size of the pellets as possible in the simulation, a random size distribution with a minimum and maximum value was 0.85 and 1.15, respectively. The pellet’s size was scaled by volume.
  • Bonded particle model: The pellet model file (“.stl” format) was imported into EDEM. Then, a box with a length, width, and height of 6 × 10−3 m, 6 × 10−3 m, and 13 × 10−3 m was created to contain the pellet model. The box and pellet model types were set to physical and virtual, respectively. The box was filled with particles, and the box and pellet model types were changed to virtual and physical, respectively. When the particle motion inside the pellet model was stable, and the external particles disappeared under gravity, the contact model of the particles was set to Hertz–Mindlin (no slip), and the connection parameters of the particles were set. Then, we started to generate bonding keys to connect the particles and deleted the box and the pellet model shell to achieve a gluten pellet model. The particle filling effect can be verified by Equation (3).
    α · r 2 · L = N · 4 3 · R 3
    where α is fill rate; r is the radius of actual gluten pellet, m; L is the length of actual gluten pellet, m; N is the number of particles filled; and R is the radius of particles filled, m.
The bonded particle model is shown in Figure 3b. Particles with a radius of 0.3×10−3 m were used for filling, and the number of particles was 839. The filling rate is 0.59, and the filling effect is ideal [28].
(2)
Compression model
After the bonded particle model was established, the indenter and the support plate models were established, respectively. In order to ensure that the quality of the bonding model was equal to the actual pellet, the density amplification method was used to enlarge the density of the pellets from 1099 kg·m−3 to 1869 kg·m−3, and the other parameters remained unchanged. The compression speed was set to 0.5 × 10−3 m·s−1 to simulate the compression process of the pellets, as shown in Figure 3c.

2.2. Test Design and Indicators Determination

2.2.1. Interaction Parameters Determination of Pellet–Stainless Steel

(1)
Static friction coefficient
The static friction coefficient measurement principle is shown in Figure 4b. Five pellets were glued together and placed on a stainless steel inclinometer to prevent the gluten pellets from rolling. The inclination angle gradually increased and was read by a digital protractor (PT180, Shenzhen Red Dragon Instruments Co., Ltd., Shenzhen, China), as shown in Figure 4a. When the pellets just began to slide on the inclined surface, the inclination angle was the static friction angle. Each group of experiments was repeated five times, and the mean value was calculated. The static friction coefficient was calculated by Equation (4) [29].
μ s = tan α
where μ is the static friction coefficient and α is the static friction angle.
The static friction coefficient between the gluten pellets and the stainless steel plate was measured to be 0.51 ± 0.02.
(2)
Rolling friction coefficient
The rolling friction coefficient between the gluten pellets and the stainless steel plate adopted the inclined surface rolling method. The measurement principle of the rolling friction coefficient is shown in Figure 4c. The gluten pellets were placed along the radial directions on the stainless steel plate. The initial velocity of the pellets was zero. When the rolling was stationary, the rolling distance was measured, and the rolling friction coefficient was calculated by Equation (5) [30]. Experiments were repeated five times, and the rolling friction coefficient between the gluten pellets and the stainless steel plate was measured to be 0.27 ± 0.004. Furthermore, the average time for each group of pellets to roll off the stainless steel plate was recorded as 0.45 s.
μ r = h 1 H 2 4 c o s 2 β ( h 2 H t a n β ) L c o s β
(3)
Restitution coefficient
The collision test method of pellet to the fixed surface was used to determine the restitution coefficient between gluten pellets and stainless steel plate. The measurement principle of the restitution coefficient is shown in Figure 4d. A stainless steel plate with a horizontal plane was set at 45°. Equations (6) and (7) were used to calculate the restitution coefficient [31,32]. Each group was repeated five times, and the average value was obtained.
v x = g s 1 s 2 ( s 1 s 2 ) 2 ( h 1 s 2 h 2 s 1 ) v y = h 1 v x s 1 g s 1 2 v x
e = v x 2 + v y 2 · c o s 45 ° + a r c t a n ( v y v x ) 2 g H · s i n 45 °
where v x is horizontal speed, m·s−1; g is gravitational acceleration, m·s−2; s 1 and s 2 are horizontal displacement for two measurements, respectively, m; h 1 and h 2 are vertical displacement for two measurements, respectively, m; v y is vertical velocity; e is restitution coefficient; and H is the drop height, m.
The restitution coefficient between the gluten pellets and the stainless steel plate was measured to be 0.50 ± 0.04.
(4)
Angle of repose
The angle of repose of gluten pellets was determined by injection [10]. The angle of repose reflects the internal friction and scattering properties of gluten pellets. It is related to pellet shape, size, moisture content, stacking conditions, etc. The larger the angle of repose, the greater the internal friction force and the smaller the scatter.
The angle of the repose measuring device is shown in Figure 5. The outlet diameter of the 304 stainless steel (shear modulus is 2.7 × 1010 Pa; Poisson’s ratio is 0.3; and density is 7930 kg·m−3) funnel is 3.05 × 10−2 m, the diameter of the plexiglass disc is 0.09 m, the height is 0.0175 m, and the distance between the funnel outlet and the pellets in the plexiglass disc is 0.05 m. During the test, first, a baffle was used to close the outlet of the funnel, and the funnel was filled with gluten pellets. Then, the baffle was quickly removed until the pellets no longer fell and a pellet pile formed in the glass disc. In order to achieve the exact angle of repose, a mobile phone was used to take a vertical photo of the front of the pellets heap, and then Python processed the image. Specific treatment methods were as follows: Firstly, the image was enhanced by Gaussian filtering to remove noise points and histogram homogenization. An appropriate threshold was set for binary segmentation, and then an expansion was carried out. The hole-filling algorithm was used three times to fill the internal holes, remove excess edges and corners, and smooth the contour. Finally, the canny operator was used to detect the edge, and the edge contour was obtained. The edge contour was imported into Origin 2021 and transformed into coordinate data by an image digitization tool, and then the linear fitting was performed. The slope obtained by linear fitting was converted into the angle of repose. The image processing is shown in Figure 6. Approximately 0.1 kg pellets were selected for each group of experiments and repeated five times. The mean value was obtained. The actual measured gluten pellets angle of repose was 28.61 ± 2.65°.

2.2.2. Interaction Parameters Determination of Pellet–Pellet

The Generic EDEM material model (GEMM) database contains thousands of material models representing various materials such as corn, coal powder, sand, etc. It addresses one of the critical challenges of DEM simulation to estimate simulation parameters by inputting material size, density, and angle of repose. The ranges of the static friction coefficient ( μ s ), rolling friction coefficient ( μ r ), and restitution coefficient ( e ) between pellets were determined by inputting the above three parameters of the gluten pellets studied in this paper and combining them with the related research in works from the literature [10,16,33], as shown in Table 1.
The steepest ascent test was carried out to further narrow the factor values range. The simulated angle of repose θ and its relative error δ with the actual angle of repose θ were indicators. The relative error was calculated by Equation (8).
δ = θ θ θ
Based on the above results, the three factors quadratic orthogonal rotation combination test was carried out. The Design-Expert 12.0 software was used to analyze test results by the design principle of CCD (Central Composite Design). Nine center points were set to acquire 23 groups of tests in total. Each group of tests was repeated five times, and the average value was taken as the final simulation result. The factor codes are shown in Table 2.

2.2.3. Bonding Parameters Calibration Test Design

Normal stiffness per unit area X1, shear stiffness per unit area X2, critical normal stress X3, critical shear stress X4, and bond disk radius X5 were taken as test factors, and breaking displacement Y1 and ultimate breaking force Y2 were taken as indicators. First, the Plackett–Burman test at 2-level was used, as shown in Table 3, and ANOVA and t-test were performed to determine the significance of the factors. Then, the steepest ascent test was carried out further to reduce the range of parameters of each factor. Finally, the Box–Behnken test was carried out to determine the optimal parameter combination of each factor.
The simulation process of extrusion rupture is shown in Figure 7. With the downward movement of the indenter, when the displacement was 1.23 × 10−3 m, the force on the gluten pellet reached the maximum value, and when the displacement was 1.74 × 10−3 m, the gluten pellet was broken entirely.

3. Results

3.1. Parameters Calibration of Pellet–Pellet Test

3.1.1. Steepest Ascent Test

The design and results of the steepest ascent tests are shown in Table 4. With the increase in μ s , μ r , and e , the simulation value of the angle of repose rose first and then it decreased. Among them, the relative error of the fourth group of tests was the smallest. Therefore, the third, fourth, and fifth groups of the parameter values were used as the low, center point, and high levels of the quadratic orthogonal rotation combination design test.

3.1.2. Quadratic Orthogonal Rotation Combination Test

The test results of the quadratic orthogonal rotation combination design are shown in Table 5. The simulation results compared different fitting models, showing that the quadratic full model fitting was the best.
The results of the significance analysis and variance analysis of the regression model are shown in Table 6. The fitting degree of the regression model was very significant (p < 0.0001). The interaction terms of the static friction coefficient with the collision restitution coefficient and the rolling friction coefficient had no significant effect on the angle of repose (p > 0.05). All the other effects were significant. The order of influence of each factor was μ s > e > μ r . The lack of fit of the regression model was p = 0.0527 > 0.05, which was not significant, indicating that there were no other main factors affecting the index in the model. The R2 and adjusted R2 of the regression model were 0.936 and 0.891, respectively, indicating that the predicted value of the regression model fitted well with the actual value. In order to improve the accuracy of the model, the regression model of the interaction parameters between the pellets and angle of repose was Equation (9), after the non-significant terms were removed.
θ = 29.16 + 2.19 μ s 1.42 e + 0.93 μ r 0.94 e μ r + 1.30 μ s 2 + 2.09 e 2 + 1.29 μ r 2
The angle of repose θ = 28.61° obtained in the actual test was input into the Optimization-Numerical module of Design-Expert 12.0 software. The parameters were further optimized in the range of 0.44 ≤ μ s ≤ 0.68, 0.39 ≤ e ≤ 0.63, and 0.11 ≤ μ r ≤ 0.17. The optimal calibration parameter combinations of the angle of repose were μ s = 0.52, e = 0.54, and μ r = 0.12. The simulation angle of repose was 28.69° using the optimal parameter combination for simulation verification, and the relative error between the simulation and the actual test was 0.28%. The comparison between the actual and the simulation test is shown in Figure 8. The results showed that the shape and angle of the simulation test of the angle of repose were consistent with the actual test, which indicated that the simulation parameters obtained by optimization were reasonable.

3.2. Bonding Parameters Calibration Test

3.2.1. Plackett–Burman Test

The 11-factor Plackett–Burman test protocol was selected, and the test results are shown in Table 7. Design-Expert 12.0 was used to conduct ANOVA and t-test on the test data, and the results are shown in Table 8 and Figure 9. The influence of each factor on breaking displacement Y1 was X5 > X4 > X2 > X1 > X3, among which X4, X5, X2, and X1 had significant effects (p < 0.05), and the effective value of X4 and X1 on the target was positive. In contrast, the effect value of X5 and X2 on the target was negative. Similarly, the influence of each factor on the breaking force Y2 was X5 > X4 > X2 > X3 > X1, among which X5 and X4 had significant effects (p < 0.05), and the effective value of X4 on the target was positive, and the effective value of X5 on the target was negative. A comprehensive analysis showed that X3 had no significant effect on the results (p > 0.05), so 8 × 106 Pa was used in the subsequent simulation test.

3.2.2. Steepest Ascent Test

According to the Plackett–Burman test results, the appropriate step size was chosen to carry out the steepest ascent test further, as shown in Table 9. The breaking displacement Y1 and the breaking force Y2 decreased first and then increased, and the errors between the simulation results and the actual values of the fourth and fifth groups were the smallest. Moreover, the actual values Y1 and Y2 were between the third and fifth simulated values.

3.2.3. Box–Behnken Test

According to the steepest ascent test results, the third and fifth sets of parameters were chosen as the high and low levels of the Box–Behnken test, and the test results and ANOVA results are shown in Table 10 and Table 11. It can be seen from Table 11 that X1, X2, X4, and X5 had significant influences on Y1 and Y2 (p < 0.05), and the influence sizes were X5 > X2 > X4 > X1 and X4 > X1 > X5> X2, respectively. The lack of fit of the regression model p = 0.1597 > 0.05 and p = 0.1013 > 0.05 were not significant. The functional relation between Y1 and the significant factors was Equation (10). The R2 and adjusted R2 of the regression model were 0.978 and 0.956, respectively. The functional relation between Y2 and the significant factors was Equation (11). The R2 and adjusted R2 of the regression model were 0.939 and 0.878, respectively. It indicated that the fitting model was reliable.
Y 1 = 48.33 12.62 X 1 37.38 X 2 1.24 X 4 + 25.77 X 5 4.08 X 1 X 2 + 3.43 X 1 X 4 14.63 X 4 X 5 + 11.95 X 2 2 + 78.53 X 5 2
Y 2 = 18485.10 4606.06 X 1 12411.53 X 2 + 519.33 X 4 17162.46 X 5 + 1370.83 X 1 X 2 + 182.26 X 1 X 4 + 8804.46 X 2 X 5 + 208.93 X 1 2 + 1646.45 X 2 2 75.63 X 4 2
In order to further determine the optimal parameter combination, the Numerical module of Design-Expert 12.0 software was used to optimize the solution. The solution boundary conditions were set as 1.88 × 1010 N·m−3X1 ≤ 2.12 × 1010 N·m−3, 1.86 × 1010 N·m−3X2 ≤ 2.00 × 1010 N·m−3, 5.8 × 106 Pa ≤ X4 ≤ 6.2 × 106 Pa, 0.36 × 10−3 m ≤ X5 ≤ 0.40 × 10−3 m, Y1 = 1.209 × 10−3 m, and Y2 = 94.656 N. It can be obtained that X1, X2, X4, and X5 were 1.988 × 1010 N·m−3, 1.931 × 1010 N·m−3, 6.010 × 106 Pa, and 0.381 × 10−3 m, respectively. Simulation tests on the optimized parameters show that Y1 and Y2 are 1.217 × 10−3 m and 93.623 N, respectively, and the relative errors with the real values are 0.66% and 1.09%, respectively.

3.3. Adaptability Verification of Model Parameters

In this study, the multi-spheres autofill model and the bonded particle model were, respectively, used to carry out the free fall test, inclined plane sliding test, and inclined plane tumbling time test, in order to verify the feasibility of using the calibration parameters of the multi-spheres autofill model on the bonded particle model.
(1) Free fall test: The pellet fell freely from a fixed height and bounced after collision with the stainless steel bottom plate. The ratio between the maximum height of the first rebound and the initial height was calculated.
(2) Inclined plane sliding test: A single pellet was statically placed on the stainless steel plate, and the stainless steel plate rotated at an angular speed of 5°·s−1 until the pellet began sliding. The rotation angle of the stainless steel plate was measured and recorded.
(3) Inclined plane rolling time test: The stainless steel plate was fixed at 45°, a single pellet was placed statically at the end of the stainless steel plate and released, and the total time was recorded until it rolled out of the stainless steel plate.
The simulation test results of the two models are shown in Table 12. It can be seen from the test results that the relative errors of the three test results of the two models were 0.578%, 0.610%, and 0.937%, respectively. Furthermore, the motion states of the two models were consistent, indicating that it is feasible to apply the interaction parameters of the multi-spheres autofill model to the bonded particle model.

4. Discussion

The variance analysis results (Table 6) show a strong correlation between the static friction coefficient, the restitution coefficient, and the rolling friction coefficient of pellet–pellet with the angle of repose. It is consistent with the research results of Peng et al. [16] but inconsistent with those of Liao et al. [34,35] and Guo et al. [36]. This is explained by the fact that agricultural materials have complex properties. The same materials have different properties, such as hygroscopic, moisture content, constituent, dimension, etc. The properties of different materials vary greatly. Generally, the greater the friction coefficient between pellets, the more excellent the resistance to downward movement during accumulation. The larger the restitution coefficient, the smaller the energy dissipation of pellets after a collision, which increases the stability time in pellet accumulation. The material upon which the pellets fall has an important influence on the angle of repose. Therefore, under the joint action of the three factors, the angle of repose will be affected to varying degrees.
The breaking displacement and ultimate breaking force were both selected as indicators, which was not similar to that reported by other authors in the bonding parameter’s calibration simulations [10,37]. At the same time, the placement state (along the diameter or axial direction) of the pellet during compression will also affect the calibration results of its bond parameters. However, considering the actual pellet end shape and the pellet stacking state in the stacking experiment (Figure 8a), it is more practical to calibrate the bonding parameters along the radial compression of the pellets. It is assumed that the material properties of the pellet are isotropic in the simulations. The actual situation is that there may be defects in the process of pellet manufacturing, such as cracks, uneven moisture content, etc., resulting in anisotropy of material properties. Therefore, standards for pellet parameter calibration shall be formulated, and uniform samples should be adopted, including size, moisture content, etc.
After extruding the gluten pellets, they must be cooled and dried in the drum dryer. The accurate calibration of the discrete element simulation parameters of the pellet feed is of great significance for optimizing the mechanical parts, contacting the pellets, and determining the drying process. Based on previous research methods and practical tests, this paper calibrates the parameters by a discrete element simulation, but there are still some limitations. In this test, only pellets with a moisture content of 10.18% (d.b.) are studied, and the change in the pellet’s properties with the moisture content is not considered. At the same time, the simulation pellet size is the average value without considering the impact of the pellet’s size. In future research, the above factors should be comprehensively considered to obtain the discrete element parameters of gluten pellets under different conditions.

5. Conclusions

(1) For gluten pellets, the parameters required for constructing the discrete element model and the bonded particle model are yet to be calibrated, such as the pellets’ material parameters, the contact parameters (pellet–pellet and pellet–stainless steel), and the bonding parameters.
(2) It has been confirmed that the interaction parameters (static friction coefficient, restitution coefficient, and rolling friction coefficient) between pellets significantly impact the angle of repose and have a very close correlation (R2 = 0.936) with the angle of repose.
(3) It has also been found that the bonding parameters (normal stiffness per unit area, shear stiffness per unit area, critical shear stress, and bond disk radius) significantly influence the breaking displacement and ultimate breaking force. Moreover, the correlation coefficients are R2 = 0.978 and R2 = 0.939, respectively.
(4) The simulation results show that the motion states of the multi-spheres autofill model and the bonded particle model with the same interaction parameters are the same. The simulations can reference pellet breakage in production, processing, and transportation.

Author Contributions

Conceptualization, Z.B.; methodology, K.C.; resources, Z.B.; software, Z.B.; validation, D.Y.; visualization, X.Z.; writing—original draft, Z.B.; writing—review and editing, K.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Chuzhou eight industrial chain strong chain reinforcement project, grant number: 2022GJ011, and by the Chuzhou science and technology project, grant number: 2022XJZD24.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data collected in this research are available when be required.

Acknowledgments

The authors are grateful to the Bilvchun Biotechnology Co., Ltd. for providing experimental materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Biesiekierski, J.R. What is gluten? J. Gastroenterol. Hepatol. 2017, 32, 78–81. [Google Scholar] [CrossRef] [Green Version]
  2. Shewry, P. What is gluten—Why is it special? Front. Nutr. 2019, 6, 101. [Google Scholar] [CrossRef]
  3. Verbauwhede, A.E.; Lambrecht, M.A.; Jekle, M.; Lucas, I.; Fierens, E.; Shegay, O.; Brijs, K.; Delcour, J.A. Microscopic investigation of the formation of a thermoset wheat gluten network in a model system relevant for bread making. Int. J. Food Sci. Technol. 2020, 55, 891–898. [Google Scholar] [CrossRef]
  4. Wang, Y.H.; Zhang, Y.R.; Yang, Y.Y.; Shen, J.Q.; Zhang, Q.M.; Zhang, G.Z. Effect of wheat gluten addition on the texture, surface tackiness, protein structure, and sensory properties of frozen cooked noodles. LWT 2022, 161, 113348. [Google Scholar] [CrossRef]
  5. Khatkar, A.B.; Kaur, A. Effect of protein incorporation on functional, thermal, textural and overall quality characteristics of instant noodles. J. Food Meas. Charact. 2018, 12, 2218–2229. [Google Scholar] [CrossRef]
  6. Samard, S.; Gu, B.Y.; Ryu, G.H. Effects of extrusion types, screw speed and addition of wheat gluten on physicochemical characteristics and cooking stability of meat analogues. J. Sci. Food Agric. 2019, 99, 4922–4931. [Google Scholar] [CrossRef] [PubMed]
  7. Speranza, B.; Bevilacqua, A.; Campaniello, D.; Sinigaglia, M.; Musaico, D.; Corbo, M.R.; Lamacchia, C. The impact of gluten friendly flour on the functionality of an active drink: Viability of Lactobacillus acidophilus in a fermented milk. Front. Microbiol. 2018, 9, 2042. [Google Scholar] [CrossRef]
  8. Basto-Silva, C.; Guerreiro, I.; Oliva-Teles, A.; Neto, B. Life cycle assessment of diets for gilthead seabream (Sparus aurata) with different protein/carbohydrate ratios and fishmeal or plant feedstuffs as main protein sources. Int. J. Life Cycle Assess. 2019, 24, 2023–2034. [Google Scholar] [CrossRef]
  9. Jang, H.Q. Study on the Mechanism and Key Technology of Rotary Extrusion Forming of Feed Powder. Ph.D. Thesis, Nanjing University of Science & Technology, Nanjing, China, 2018. [Google Scholar]
  10. Niu, Z.Y.; Kong, X.R.; Sheng, B.S.; Li, H.C.; Geng, J.; Liu, J. Parameters Calibration of Discrete Element Simulation for Pellet Feed Attrition. Trans. Chin. Soc. Agric. Mach. 2022, 53, 132–140+207. [Google Scholar]
  11. Liu, X.P.; Li, J.W.; Liu, Y.Z.; Li, B.; Yang, L.; Zhang, Y.L. Contact Parameters Calibration and Flowing Property Analysis of Granular Feed. Feed Ind. 2021, 42, 6–17. [Google Scholar] [CrossRef]
  12. Cundall, P.A. A computer model for simulating progressive, large-scale movement in blocky rock system. In Proceedings of the International Symposium on Rock Mechanics, Nancy, France, 4–6 October 1971. [Google Scholar]
  13. Xie, C.; Yang, J.; Wang, B.; Zhuo, P.; Li, C.; Wang, L. Parameter calibration for the discrete element simulation model of commercial organic fertilizer. Int. Agrophys. 2020, 35, 107–117. [Google Scholar] [CrossRef]
  14. Ghodki, B.M.; Patel, M.; Namdeo, R.; Carpenter, G. Calibration of discrete element model parameters: Soybeans. Comput. Part Mech. 2019, 6, 3–10. [Google Scholar] [CrossRef]
  15. Zhang, R.; Han, D.L.; Qiao, L.; Ji, Q.L.; He, Y.; Li, J.Q. Calibration methods of sandy soil parameters in simulation of discrete element method. Trans. Chin. Soc. Agric. Mach. 2017, 48, 49–56. [Google Scholar] [CrossRef]
  16. Peng, F.; Wang, H.Y.; Fang, F.; Liu, Y.D. Calibration of Discrete Element Model Parameters for Pellet Feed Based on Injected Section Method. Trans. Chin. Soc. Agric. Mach. 2018, 49, 140–147. [Google Scholar] [CrossRef]
  17. Zhao, L.; Zhou, H.; Xu, L.; Song, S.; Zhang, C.; Yu, Q. Parameter calibration of coconut bran substrate simulation model based on discrete element and response surface methodology. Powder Technol. 2022, 395, 183–194. [Google Scholar] [CrossRef]
  18. Kong, X.; Liu, J.; Yang, T.; Su, Y.; Geng, J.; Niu, Z. Numerical simulation of feed pellet breakage in pneumatic conveying. Biosyst. Eng. 2022, 218, 31–42. [Google Scholar] [CrossRef]
  19. Wang, H.; Wu, P.; He, H.; Ma, Y.; Bu, K.; Xue, J. Calibration of Parameters for Discrete Element Simulation Model for Alfalfa with Different Moisture Contents Based on Angle of Repose Test. BioResources 2022, 17, 1467–1484. [Google Scholar] [CrossRef]
  20. Li, H.; Zeng, R.; Niu, Z.; Zhang, J. A calibration method for contact parameters of maize kernels based on the discrete element method. Agriculture 2022, 12, 664. [Google Scholar] [CrossRef]
  21. Horabik, J.; Wiącek, J.; Parafiniuk, P.; Stasiak, M.; Bańda, M.; Kobyłka, R.; Molenda, M. Discrete element method modelling of the diametral compression of starch agglomerates. Materials 2020, 13, 932. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kanmani, N.; Romano, N.; Ebrahimi, M.; Amin, S.N.; Kamarudin, M.S.; Karami, A.; Kumar, V. Improvement of feed pellet characteristics by dietary pre-gelatinized starch and their subsequent effects on growth and physiology in tilapia. Food Chem. 2018, 239, 1037–1046. [Google Scholar] [CrossRef]
  23. Styks, J.; Wróbel, M.; Frączek, J.; Knapczyk, A. Effect of Compaction Pressure and Moisture Content on Quality Parameters of Perennial Biomass Pellets. Energies 2020, 13, 1859. [Google Scholar] [CrossRef] [Green Version]
  24. Gageanu, I.; Cujbescu, D.; Persu, C.; Tudor, P.; Cardei, P.; Matache, M.; Vladut, V.; Biris, S.; Voicea, I.; Ungureanu, N. Influence of Input and Control Parameters on the Process of Pelleting Powdered Biomass. Energies 2021, 14, 4104. [Google Scholar] [CrossRef]
  25. Behnke, K.C. Factors influencing pellet quality. Feed Tech. 2001, 5, 19–22. [Google Scholar]
  26. Said, N.; Daiem, M.; García-Maraver, A.; Zamorano, M. Influence of densification parameters on quality properties of rice straw pellets. Fuel Process. Technol. 2015, 138, 56–64. [Google Scholar] [CrossRef]
  27. American Society of Agricultural and Biological Engineers. S352. 2: Moisture Measurement—Unground Grain and Seeds; American Society of Agricultural and Biological Engineers: Saint Joseph, MI, USA, 2006. [Google Scholar]
  28. Wang, M.M.; Wang, W.Z.; Yang, L.G.; Hou, M.T. Research of discrete element modeling method of maizekernel based on EDEM. J. Henan Agric. Univ. 2018, 52, 80–84. [Google Scholar]
  29. Sun, J.X.; Guo, Y.M.; Yang, Z.M.; Cui, Q.L.; Wu, X.H.; Zhang, Y.Q. Experimental study on biomechanical properties of buckwheat grain and viscoelastic properties of buckwheat powder. Trans. Chin. Soc. Agric. Eng. 2018, 34, 287–298. [Google Scholar] [CrossRef]
  30. Li, B.; Cheng, Y.; Sun, P.; Li, B.; Liu, W.F. Engineering Measurement and Verification of Rolling Friction Coefficients. Constr. Mach. Equip. 2017, 48, 29–32. [Google Scholar] [CrossRef]
  31. Xing, J.J.; Zhang, R.; Wu, P.; Zhang, X.R.; Dong, X.H.; Chen, Y.; Ru, S.F. Parameter calibration of discrete element simulation model for latosol particles in hot areas of Hainan. Trans. Chin. Soc. Agric. Eng. 2020, 36, 158–166. [Google Scholar]
  32. Feng, B.; Sun, W.; Shi, L.R.; Sun, B.G.; Zhang, T.; Wu, J.M. Determination of restitution coefficient of potato tubers collision in harvest and analysis of its influence factors. Trans. Chin. Soc. Agric. Eng. 2017, 33, 50–57. [Google Scholar] [CrossRef]
  33. Zhang, X.H. Computer Numerical Simulation of Pellet Feed in Cooling Process. Ph.D. Thesis, Huazhong Agricultural University, Wuhan, China, 2004. [Google Scholar]
  34. Liao, Y.T.; Wang, Z.T.; Liao, Q.X.; Wan, X.Y.; Zhou, Y.; Liang, F. Calibration of discrete element model parameters of forage rape stalk at early pod stage. Trans. Chin. Soc. Agric. Mach. 2020, 51, 236–243. [Google Scholar]
  35. Liao, Y.T.; Liao, Q.X.; Zhou, Y.; Wang, Z.T.; Jang, Y.J.; Liang, F. Parameters calibration of discrete element model of fodder rape crop harvest in bolting stage. Trans. Chin. Soc. Agric. Mach. 2020, 51, 73–82. [Google Scholar] [CrossRef]
  36. Guo, S.Q.; Wan, P.; Yang, J.; Deng, Z.Y.; Wang, S.S.; Zhu, M.; Tan, H.Q. Discrete element based calibration and test of simulation parameters for fish feed. J. Huazhong Agric. Univ. 2022, 41, 286–295. [Google Scholar] [CrossRef]
  37. Peng, F.; Huang, Z.; Fang, F. Modeling and experiments of chewing mechanical properties of pellet feed using discrete element method. Int. J. Agric. Biol. Eng. 2020, 13, 37–44. [Google Scholar] [CrossRef]
Figure 1. Texture compression experiment.
Figure 1. Texture compression experiment.
Agriculture 13 00744 g001
Figure 2. Force–displacement curve.
Figure 2. Force–displacement curve.
Agriculture 13 00744 g002
Figure 3. Discrete element model of gluten pellet. (a) Multi-spheres autofill model; (b) bonded particle model; and (c) compression simulation model.
Figure 3. Discrete element model of gluten pellet. (a) Multi-spheres autofill model; (b) bonded particle model; and (c) compression simulation model.
Agriculture 13 00744 g003
Figure 4. The measurement detection device and principle. (a) The measurement detection device; (b) static friction coefficient measurement principle; (c) rolling friction coefficient measurement principle; and (d) restitution coefficient measurement principle.
Figure 4. The measurement detection device and principle. (a) The measurement detection device; (b) static friction coefficient measurement principle; (c) rolling friction coefficient measurement principle; and (d) restitution coefficient measurement principle.
Agriculture 13 00744 g004
Figure 5. The angle of repose measuring device.
Figure 5. The angle of repose measuring device.
Agriculture 13 00744 g005
Figure 6. The angle of repose image processing. (a) Original image; (b) binary image; (c) contour extraction; and (d) contour line fitting.
Figure 6. The angle of repose image processing. (a) Original image; (b) binary image; (c) contour extraction; and (d) contour line fitting.
Agriculture 13 00744 g006
Figure 7. The simulation process of extrusion rupture. (a) Compression 0 mm; (b) compression 1.23 × 10−3 m; and (c) compression 1.74 × 10−3 m.
Figure 7. The simulation process of extrusion rupture. (a) Compression 0 mm; (b) compression 1.23 × 10−3 m; and (c) compression 1.74 × 10−3 m.
Agriculture 13 00744 g007
Figure 8. Comparison of actual test and simulation test. (a) Actual test. (b) Simulation test.
Figure 8. Comparison of actual test and simulation test. (a) Actual test. (b) Simulation test.
Agriculture 13 00744 g008
Figure 9. Pareto charts. (a) Breaking displacement. (b) Ultimate breaking force.
Figure 9. Pareto charts. (a) Breaking displacement. (b) Ultimate breaking force.
Agriculture 13 00744 g009
Table 1. The ranges of the simulation parameters of pellet–pellet.
Table 1. The ranges of the simulation parameters of pellet–pellet.
ParametersRange
μ s 0.20–0.80
μ r 0.05–0.20
e 0.15–0.75
Table 2. Factors and codes of quadratic regression orthogonal rotating test.
Table 2. Factors and codes of quadratic regression orthogonal rotating test.
FactorsCodes
−1.682−1011.682
μ s 0.440.490.560.630.68
e 0.390.440.510.580.63
μ r 0.110.120.140.160.17
Table 3. Plackett–Burman test factors and levels.
Table 3. Plackett–Burman test factors and levels.
FactorsValue
Low LevelHigh Level
X1/1010 N·m−31.642.24
X2/1010 N·m−31.82.12
X3/106 Pa7.258.5
X4/106 Pa5.46.4
X5/10−3 m0.340.44
Table 4. Results of the steepest ascent test.
Table 4. Results of the steepest ascent test.
Test No.FactorsResults
μ s μ r e θ δ /%
10.200.050.1536.0626.03
20.320.080.2733.4917.04
30.440.110.3926.646.88
40.560.140.5127.683.25
50.680.170.6330.295.89
60.800.200.7530.928.07
Table 5. Results of quadratic regression orthogonal rotating test.
Table 5. Results of quadratic regression orthogonal rotating test.
Test No. μ s e μ r θ δ
10.490.440.1230.566.82
20.630.440.1235.2823.31
30.490.580.1228.11.78
40.630.580.1234.5320.69
50.490.440.1635.3823.66
60.630.440.1639.1836.95
70.490.580.1629.161.92
80.630.580.1634.6521.11
90.440.510.1430.717.34
100.680.510.1436.3627.09
110.560.390.1437.3830.65
120.560.630.1434.1719.43
130.560.510.1132.6514.12
140.560.510.1734.3420.03
150.560.510.1429.733.91
160.560.510.1429.563.32
170.560.510.1430.075.10
180.560.510.1428.21.43
190.560.510.1429.091.68
200.560.510.1428.931.12
210.560.510.1429.171.96
220.560.510.1429.794.12
230.560.510.1427.653.36
Table 6. Variance analysis of angle of repose regression model.
Table 6. Variance analysis of angle of repose regression model.
Source of VariationSum of SquaresDegree of FreedomF-Valuep-Value
Model235.47921.05<0.0001 **
μ s 65.65152.82<0.0001 **
e 27.44122.080.0004 **
μ r 11.8919.570.0086 **
μ s e 1.4511.160.3005
μ s μ r 0.432510.34790.5654
e μ r 7.1115.720.0326 *
μ s 2 26.96121.690.0004 **
e 2 69.71156.09<0.0001 **
μ r 2 26.38121.220.0005 **
Residual16.1613
Lack of fit11.2053.610.0527
Pure error4.968
Cor Total251.6322
Note: p < 0.05 (significant, *), p < 0.01 (highly significant, **).
Table 7. Plackett–Burman Test Results.
Table 7. Plackett–Burman Test Results.
Test No.X1/1010 N·m−3X2/1010 N·m−3X3/106 PaX4/106 PaX5/10−3 mY1/10−3 mY2/N
12.242.127.256.40.441.17182.171
21.642.128.55.40.441.05269.787
32.241.88.56.40.341.336102.652
41.642.127.256.40.441.18385.844
51.641.88.55.40.441.10774.813
61.641.87.256.40.341.30294.112
72.241.87.255.40.441.15168.353
82.242.127.255.40.341.192103.427
92.242.128.55.40.341.19295.357
101.642.128.56.40.341.262110.69
112.241.88.56.40.441.25375.765
121.641.87.255.40.341.22490.113
Table 8. Analysis of significance of parameters in Plackett–Burman test.
Table 8. Analysis of significance of parameters in Plackett–Burman test.
Source of VariationY1Y2
Sum of SquaresDegree of FreedomMean SquareF-Valuep-ValueSum of SquaresDegree of FreedomMean SquareF-Valuep-Value
Model0.068950.013859.93<0.0001 **1973.555394.7114.000.0029 **
X10.002310.00239.870.0200 *0.466510.46650.01650.9018
X20.008610.008637.340.0009 **143.301143.305.080.0650
X30.000010.00000.15980.70322.1212.120.07520.7931
X40.028910.0289125.71<0.0001 **203.231203.237.210.0363 *
X50.029110.0291126.57<0.0001 **1624.4311624.4357.630.0003 **
Residual0.001460.0002 169.13628.19
Cor Total0.070311 2142.6811
Note: p < 0.05 (significant, *), p < 0.01 (highly significant, **).
Table 9. Results of the steepest ascent test.
Table 9. Results of the steepest ascent test.
Test No.X1/1010 N·m−3X2/1010 N·m−3X4/106 PaX5/10−3 mY1/10−3 m δ 1 Y2/N δ 2
11.642.125.40.441.1207.43865.00731.326
21.762.065.60.421.1435.53766.57029.675
31.882.005.80.401.1802.47987.0378.053
42.001.9460.381.2331.92890.9433.926
52.121.866.20.361.2805.78595.7471.148
62.241.806.40.341.3077.989101.9437.694
Table 10. Results of the Box–Behnken test.
Table 10. Results of the Box–Behnken test.
Test No.X1/1010 N·m−3X2/1010 N·m−3X4/106 PaX5/10−3 mY1/10−3 mY2/N
11.881.8660.381.177119.803
22.121.8660.381.28688.246
31.88260.381.322100.241
42.12260.381.294114.744
521.935.80.361.22899.062
621.936.20.361.36592.517
721.935.80.41.23191.846
821.936.20.41.13483.291
91.881.9360.361.283102.673
102.121.9360.361.29394.935
111.881.9360.41.22492.116
122.121.9360.41.21196.124
1321.865.80.381.23699.242
14225.80.381.297111.627
1521.866.20.381.23788.913
16226.20.381.33794.114
171.881.935.80.381.274109.722
182.121.935.80.381.12589.003
191.881.936.20.381.14294.014
202.121.936.20.381.32290.792
2121.8660.361.301119.643
222260.361.39695.811
2321.8660.41.21887.514
242260.41.297112.987
2521.9360.381.22697.913
2621.9360.381.20893.764
2721.9360.381.21394.407
2821.9360.381.21294.463
2921.9360.381.20192.757
Table 11. Box–Behnken test program and results.
Table 11. Box–Behnken test program and results.
Source of Variation Y1 Y2
Sum of SquaresDegree of FreedomMean SquareF-Valuep-ValueSum of SquaresDegree of FreedomMean SquareF-Valuep-Value
Model0.1215140.008744.33<0.0001 **2498.4414178.4615.42<0.0001 **
A-X10.001010.00105.060.0411 *166.691166.6914.410.0020 **
B-X20.019810.0198101.37<0.0001 **57.04157.044.930.0434 *
C-X40.001810.00189.070.0093 *269.431269.4323.290.0003 **
D-X50.025310.0253129.24<0.0001 **138.471138.4711.970.0038 **
AB0.004710.004723.970.0002 **530.381530.3845.84<0.0001 **
AC0.027110.0271138.23<0.0001 **76.54176.546.610.0222 *
AD0.000110.00010.67550.424934.49134.492.980.1062
BC0.000410.00041.940.185112.90112.901.120.3089
BD0.000110.00010.32690.5765607.751607.7552.53<0.0001 **
CD0.013710.013769.92<0.0001 **1.0111.010.08730.7720
A20.000110.00010.41560.529658.91158.915.090.0406 *
B20.022710.0227115.99<0.0001 **422.721422.7236.53<0.0001 **
C21.126 × 10−611.126 × 10−60.00580.940659.16159.165.110.0402 *
D20.006710.006734.02<0.0001 **0.004710.00470.00040.9842
Residual0.0027140.0002 161.991411.57
Lack of Fit0.0024100.00022.880.1597146.881014.693.890.1013
Pure Error0.000340.0001 15.1143.78
Cor Total0.124228 2660.4328
Note: p < 0.05 (significant, *), p < 0.01 (highly significant, **).
Table 12. Simulation results of two pellet models.
Table 12. Simulation results of two pellet models.
TsetInvestigation ParameterModelRelative Error/%
Multi-Spheres Autofill ModelBonded Particle Model
Collision testInitial height/10−3 m128.647128.5740.578
Post collision height/10−3 m64.24263.873
Sliding testInitial sliding angle/°26.73326.5710.610
Rolling testRolling off time/s0.4310.4270.937
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ben, Z.; Zhang, X.; Yang, D.; Chen, K. An Experimental and Numerical Study for Discrete Element Model Parameters Calibration: Gluten Pellets. Agriculture 2023, 13, 744. https://doi.org/10.3390/agriculture13040744

AMA Style

Ben Z, Zhang X, Yang D, Chen K. An Experimental and Numerical Study for Discrete Element Model Parameters Calibration: Gluten Pellets. Agriculture. 2023; 13(4):744. https://doi.org/10.3390/agriculture13040744

Chicago/Turabian Style

Ben, Zongyou, Xubo Zhang, Duoxing Yang, and Kunjie Chen. 2023. "An Experimental and Numerical Study for Discrete Element Model Parameters Calibration: Gluten Pellets" Agriculture 13, no. 4: 744. https://doi.org/10.3390/agriculture13040744

APA Style

Ben, Z., Zhang, X., Yang, D., & Chen, K. (2023). An Experimental and Numerical Study for Discrete Element Model Parameters Calibration: Gluten Pellets. Agriculture, 13(4), 744. https://doi.org/10.3390/agriculture13040744

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop