Next Article in Journal
An Environmental and Economic Analysis of Strawberry Production in Southern Italy
Previous Article in Journal
The Impact of Climate Change on the Urban–Rural Income Gap in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Material and Interaction Properties of Granular Fertilizer Particles Using DEM Simulation and Bench Testing

1
College of Mechanical and Electrical Engineering, Sichuan Agricultural University, Ya’an 625014, China
2
Smart Agriculture Engineering Technology Center of Sichuan Province, Ya’an 625014, China
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(9), 1704; https://doi.org/10.3390/agriculture13091704
Submission received: 25 June 2023 / Revised: 23 August 2023 / Accepted: 25 August 2023 / Published: 29 August 2023
(This article belongs to the Section Agricultural Technology)

Abstract

:
The discrete element method (DEM) is an effective tool for obtaining qualitative and quantitative information on particle motion, which aids in the design and optimization of agricultural equipment structures. The accuracy of the DEM simulation parameters significantly impacts the simulation results. This study employed a combination of high-speed camera measurement, DEM simulation, and validation tests to determine the material and interaction property parameters for fertilizer particles. The basic parameters (triaxial size, bulk density, density, and coefficient of static friction) and coefficients of restitution between fertilizer and material were measured for three fertilizer varieties. There was a significant difference in the angle of repose between various material plates and fertilizer particles. The calibration values of coefficients of restitution and coefficients of rolling friction between fertilizer particles were optimized using the Box–Behnken method. The angle of repose was significantly affected by the coefficient of static friction and the coefficient of rolling friction between the fertilizer particles. The determined values for the coefficient of restitution, coefficient of static friction, and coefficient of rolling friction between the fertilizer particles were 0.323, 0.381, and 0.173, respectively. The error in the angle of the repose test was less than 3.0%, and the variation coefficient for each row consistency was less than 1.68 percentage points under the optimal simulation parameters. DEM simulations of the angle of repose and each row consistency variation coefficient test using the measured parameters can accurately predict the experimental results. The findings of this paper provide a theoretical basis for the DEM study of fertilizer particles.

1. Introduction

Fertilizer has provided balanced nutrition and improved soil chemical and physical properties. Precision fertilizer application is an effective way to increase crop yield and quality [1]. Granular compound fertilizer, which contains nutrients such as N, P2O5, and K2O and is adaptable for mechanized fertilization application, is widely used in agricultural production. Due to the current environmental issues and farmers’ economic needs, precision fertilization relying on fertilizer feeders has become a key component in crop cultivation [2].
With the characteristics of low moisture absorption, high fluidity, and good bonding, it is necessary to have a further understanding of the characteristics of granular compound fertilizer to aid in the design of fertilizer feeders and the understanding of the relationship between fertilizer particles and the feeder. The physical properties of granular compound fertilizer particles mainly include shape, 3-D size, particle friction coefficients, stiffness, and coefficient of restitution. However, due to the complexity of granular compound fertilizer properties, designing a fertilizer feeder still relies heavily on experience and experimental methods, which can be both costly and time-consuming [3]. To reduce the cost and time required for optimal equipment design, the discrete element method (DEM) has been considered an effective simulation method for dealing with granular particles in agriculture [4]. Previous researchers have designed fertilizer spreaders based on the DEM approach [5,6,7,8]. Coetzee et al. [9] investigated the effects of fertilizer feed position, feed rate, and disc structure on the spread pattern of a centrifugal fertilizer using the DEM model. Liedekerke et al. [10] applied DEM simulations to analyze the impacts of fertilizer physical properties on the spread pattern of spinning disc spreaders. Additionally, DEM modeling was proposed to optimize the structure of a dual-band fertilizer applicator and a multi-fertilizer mixing device [11,12]. DEM simulation has proven to be a useful tool for solving agricultural engineering problems and obtaining quantitative information on particles [13].
In the DEM simulation process, the accuracy of simulation parameters determines the accuracy and reliability of results. Some simulation parameters, such as the static friction coefficient, rolling friction coefficient, and coefficient of restitution, are basic and key parameters required in DEM simulations. These coefficients have been determined by methods such as numerical analysis and bench testing. Horabik et al. [14] reported the coefficient of restitution of peas, soybeans, and rapeseed using high-speed camera technology. Huang et al. [15] determined the coefficients of restitution of rapeseed for collisions between rapeseed and harvesting equipment. Wang et al. [16] investigated the coefficient of rolling friction of maize particles and found that the angle of repose was a suitable index to verify the accuracy of the simulation model. Furthermore, Liu et al. [17] analyzed the boundary parameters of the DEM simulation of urea particles and indicated that the static friction coefficient between urea particles and plates significantly affected the simulation results. Lu et al. [18] optimized the DEM simulation parameters of the static friction coefficient between bud seeds and stainless plate and the rolling friction coefficient between bud seeds. Irazábal et al. [19] investigated bounded rolling friction for spherical discrete element particles, and Persson and Frenning [20] analyzed the effects of rolling friction on the shear behavior of granular particles.
Although previous research has analyzed some characteristics of granular fertilizers and fertilizer spreaders that have been optimally designed based on DEM simulations, the parameters of granular fertilizers are mainly empirical, and some boundary parameters remain unclear due to the limitations of experimental conditions. In this work, we investigated the particle static friction coefficient, rolling friction coefficient, and restitution coefficient of granular compound fertilizer by bench testing using a high-speed camera and DEM simulations. The angle of repose and fertilizer performance results were utilized to verify the accuracy of the simulation parameters. This work can improve the accuracy of the simulation model for granular compound fertilizer and aid in the design of fertilizer spreaders.

2. Materials and Methods

2.1. Test Method of Material and Interaction Properties of Granular Fertilizer Particle

2.1.1. Shape and Size Analysis of Granular Fertilizer

In this work, we tested three types of granular compound fertilizers that are widely used in crop production. These types are Zhushang (Summit Fertilizer Company Limited, Qingdao, China), Zhonghua (Sinochem Fertilizer Company Limited, Linyi, Shandong, China), and Sanning (Sanning Fertilizer Company Limited, Yichang, Hubei, China), as shown in Table 1. We measured the triaxial dimensions, namely length (L), width (W), and thickness (T), using a digital vernier caliper with an accuracy of 0.01 mm, and the size distribution is shown in Figure 1 with samples of 100 particles. Table 1 represents the triaxial size for three types granular compound fertilizer.
Figure 1 shows that the triaxial dimensions of granular compound fertilizer follow a normal distribution. There were significant differences in the triaxial dimensions among the three types of granular compound fertilizer. The length and width of the Sanning type were lower than those of the Zhonghua and Zhushang types. The length, width, and thickness mainly varied between 3.5–5.1 mm, 3.3–4.3 mm, and 2.9–4.1 mm, respectively. The sphericity of these types was over 89%, indicating that the granular compound fertilizer was elliptical.

2.1.2. Density

The density of fertilizer particles was measured by the weighting method. An appropriate amount of fertilizer particles with a known mass was poured into a measuring cylinder with 500 mL of distilled water. The liquid level reading of the measuring cylinder was recorded with five replications. The average density of fertilizer particles was obtained as 1542 kg/m3.

2.1.3. Coefficient of Restitution

The test principle of the coefficient of restitution is presented in Figure 2. The coefficient of restitution is defined as the absolute ratio of the instantaneous velocity after collision to that before collision. When a particle was dropped from a height (h), the coefficient of restitution was calculated using Equation (1).
e = v v 0 = v x 2 + v y 2 + v z 2 2 g h
where e is the coefficient of restitution; v is an instantaneous velocity after a collision, m s−1; v0 is an instantaneous velocity before a collision, m s−1; vx, vy, and vz are the velocity in x-direction, y-direction, and z-direction, respectively; and h is the dropping height, m.
Figure 3 shows the process for measuring the coefficient of restitution. The experimental setup consists of a mirror (1), a first coordinate plane (2), a second coordinate plane (3), a base plate (4), a vacuum nozzle (5), a test plate (7), and a high-speed camera (8). In the experimental setup, the angle between the first and second coordinate planes was 90°, and the angle between the first coordinate plane and the mirror was 135°. The first and second coordinate planes were used to measure the particle coordinates in the three-axis directions. The velocities of the two planes were obtained through a mirror (1) and a high-speed camera (FASTCAM Mini UX100; Photron Limited; Tokyo, Japan). The particle, held by a vacuum nozzle (5), was dropped freely from a height (h). The test plate (7) materials included ABS (Acrylonitrile-butadiene-styrene copolymer), 304 stainless steel, and aluminum. The particle movement trajectories were extracted from images. The programming process for extracting the trajectory is shown in Figure 4. The velocities and coefficients of restitution were obtained by extracting the particle trajectory.

2.1.4. Coefficient of Static Friction

The multifunctional comprehensive test platform (Figure 5) was utilized to measure the coefficient of static friction. The platform consisted of a 304 stainless steel (or acrylonitrile-butadiene-styrene copolymer, aluminum) test plate, which was bonded to the inclined plate. Initially, the fertilizer particles were placed at the upper end of the inclined test plate, and then the drive rod was rotated to gradually raise the inclined plate. When the fertilizer particle began to slide, the rotation of the drive rod was halted, and the inclination angle (θ) was recorded. The coefficient of static friction was calculated as μ = tanθ. The measured coefficients of static friction between fertilizer particles and 304 stainless steel, acrylonitrile-butadiene-styrene copolymer, and aluminum were found to be in the range of 0.47–0.51, 0.45–0.50, and 0.47–0.50, respectively.

2.2. DEM Simulation and Validation Experiment Method

2.2.1. Particle-Contact Model

EDEM 2018 (provided by DEM Solutions Limited, Edinburgh, UK) was used to verify the accuracy of the simulation parameters. The behavior of granular particles is governed by interactions between particles and between particles and surfaces. The surface of granular fertilizer particles was sooth and the hardness of fertilizer particles was high. Thus, the Hertz–Mindlin (no-slip) model, as shown in Figure 6, was chosen as the particle contact model, and the granular compound fertilizer was modeled as a hard-sphere model in the simulation. The particle-contact model was based on the triaxial dimensions.
In the DEM formulation, the particle dynamics equations defined the translational and rotational motion of particles. These equations could be expressed as follows:
m i υ ¨ = m i g + j = 1 n F n + F n d + F t + F t d
where mi is the mass of particle i, υ ¨ is the particle’s centroid acceleration, g is the vector of gravitational acceleration, F n and F t is the normal force and tangential force, respectively, and F n d and F t d is the normal damping force and tangential damping force, respectively.
The normal force F n was defined as follows:
F n = 4 3 E * R * 1 2 δ n 3 2
where E * was the equivalent elastic modulus, R * was the equivalent radius, and δ n was the normal overlap amount.
F n d = 2 5 6 ln γ ln γ 2 + π 2 S n m * υ n
where γ was the coefficient of restitution, Sn was the normal stiffness, m * was the equivalent mass, and υ n was the relative normal velocity.
F t = S n δ t
where, δ t was the tangential overlap amount.
F t d = 2 5 6 ln γ ln γ 2 + π 2 S t m * υ t
where St was the tangential stiffness, and υ t was the relative tangential velocity.
I i ω ˙ = j = 1 n T t + T r
T t = R i F t + F t d
where, I i was the moment of the inertia of particle I, ω ˙ was the angular acceleration, T t was the tangential torque, T r was the rolling friction torque, R i and was the distance from the barycenter of particle i to the contact point.
T r = β R r F n ω ω
where β was the rolling friction coefficient, R r was the effective rolling contact radius, and ω was the relative angular velocity of two particles in contact.

2.2.2. DEM Simulation Model

The 3D model of the fertilizer is depicted in Figure 7, where the fertilizer particle is constructed by filling overlapping spheres to satisfy the particle parameters specified in Table 1. According to the particle size distribution, the particle size distribution coefficient ranged from 0.9 to 1.1. The main simulation parameters include particle characteristics and particle-material parameters, as listed in Table 2.

2.2.3. DEM Simulation Design

As some parameters were measured or from previous research [14,22], the values of uncertain parameters included the coefficient of restitution, the coefficient of static friction, the coefficient of rolling friction between fertilizer particles, and the coefficient of static friction of fertilizer-plate material. The four factors with three levels of experimental design were carried out using the Box-Behnken method. The coding of factors is shown in Table 3. The Design-Expert version 13.0 was used to analyze the regression model and variance.

2.2.4. Simulation Method and Validation Experiment of Angle of Repose

The angle of repose is an important parameter for verifying the accuracy of the parameters. The angle of repose was tested using the hollow cylinder method, which consists of a hollow cylinder and a bottom plate. The bottom plate material included 304 stainless steel, acrylonitrile-butadiene-styrene copolymer (ABS), and aluminum. The diameter and height of the hollow cylinder were 42 mm and 126 mm, respectively. The procedure for the hollow cylinder test was as follows in the simulations: 5000 fertilizer particles were poured into the cylinder, and the particle surface was flattened first (Figure 8a). Second, the cylinder was lifted at a speed of 50 mm per second, causing the particles to flow out of the cylinder and contact the bottom plate. Finally, when the particles were stable, a cone-shaped pile was formed, as shown in Figure 8b. Each group of tests was repeated four times.
In the validation experiment, the cylinder (5) was lifted by the lifting frame of a universal material testing machine (3). The lifting velocity was controlled by the hydraulic system in the universal material testing machine (3), which was operated through the controller (4) and PC (2). In front of the universal material testing machine (3), a high-speed camera system (1) was positioned to capture images of the cone-shaped pile (Figure 9). Each treatment was repeated five times.
Figure 10 presents the image processing method for measuring the repose angle. The repose angle in both DEM simulations and bench tests was measured using MATLAB. The original image was imported into MATLAB software (Version 2012R), and then the image was subjected to gray processing and binary processing. The outlines of the particles were extracted, and a fitting line for the repose angle was obtained. The repose angle was calculated according to Wang et al. [16].

2.2.5. Validation Experiment of Fertilization Uniformity

To further validate the accuracy of simulation parameters, the fertilization process of a fertilizer metering device was investigated using DEM simulation (Figure 11b) and bench tests (Figure 11a). The Zhonghua fertilizer was tested in the bench test. The fertilizer metering device was designed by our research group and consisted of a fertilizer box (1), fertilizer metering mechanism (2), and fertilizer tube (3). The fertilizer metering mechanism (2) was a key component of the fertilizer metering device and consisted of three spiral groove-type fertilizer feeding wheels. Three rows of fertilizer were discharged by the three spiral groove-type fertilizer feeding wheels through fertilizer tubes. The fertilizer mass was measured at the 10–50 rpm rotational speed of the fertilizer metering mechanism for both the DEM simulation and bench test. The fertilization mass in each fertilizer tube was collected. The coefficient of variation for each row was calculated to evaluate the uniformity of fertilization across the three rows. The coefficient of variation for each row (cv) was calculated as follows:
c v = S D X ¯ × 100 % S D = i = 1 n ( X i X ¯ ) 2 / ( n 1 ) 1 2 X ¯ = 1 n i = 1 n X i ( i = 1 , 2 , 3 )
where, cv was the coefficient of variation for each row, %; SD was the standard deviation of fertilization mass, g; X i was the fertilization mass of the ith fertilizer tube, g; X ¯ was the average fertilization mass for each fertilizer tube, g; and n was the total number of samples.

3. Results and Analysis

3.1. Coefficient of Restitution and Angle of Repose of Experimental Study

Table 4 presents the coefficients of restitution (COR) for granular compound fertilizer against different types of materials. The CORs between the fertilizer particle and various plate materials showed significant differences. The COR between the fertilizer particle and ABS was markedly larger than the CORs between the fertilizer particle and 304 stainless steel and aluminum. The coefficient of restitution decreased with increasing dropping height, and the deviation was within 0.03 for the same material. No significant difference in COR was observed for fertilizer types or dropping height. The COR between the fertilizer particle and 304 stainless steel, ABS, and aluminum was 0.35, 0.52, and 0.32, respectively.
The angle of repose between fertilizer particles and 304 stainless steel, ABS, and aluminum is shown in Table 5. The average angle of repose between fertilizer particles and 304 stainless steel, ABS, and aluminum was 33.78°, 31.34°, and 30.70°, respectively.

3.2. Parameters Determination in DEM Simulation

Table 6 represents the effects of the coefficient of restitution (x1), the coefficient of static friction (x2), the coefficient of rolling friction (x3) between fertilizer particles, and the coefficient of static friction of fertilizer-plate material (x4) on the angle of repose in the simulation results between fertilizer particles and 304 stainless steel. The reduced quadratic model was established between the angle of repose and parameters as Equation (11).
y = 30.63 + 0.0079 x 1 + 1.62 x 2 0.5606 x 3 + 5.40 x 4 0.4856 x 1 x 3 0.6025 x 2 x 4 0.6338 x 3 x 4 1.26 x 2 2 0.0717 x 3 2 0.8186 x 4 2
Variance analysis for the reduced quadratic model is shown in Table 7. The p-value of the quadratic model was less than 0.001, and the adequacy precision was 21.578, indicating that the quadratic model reflected the relationship between the angle of repose and simulation parameters. Furthermore, the angle of repose was significantly affected by the coefficient of static friction and the coefficient of rolling friction between fertilizer particles.
The simulation parameters were optimized to achieve an objective angle of repose of 33.78°, which was obtained from the hollow cylinder bench test. Using numerical optimization analysis, the coefficients of restitution, static friction, and rolling friction between fertilizer particles were determined to be 0.323, 0.381, and 0.173, respectively. The coefficient of static friction between fertilizer particles and 304 stainless steel was 0.453.

3.3. Parameters Verification between DEM Simulation and Bench Test

To verify the accuracy of EDM simulation parameters, the optimal simulation parameters of fertilizer particles and fertilizer particles with 304 stainless steel were imported into the EDEM software (Version 2018). The average angle of repose was found to be 34.02° in the DEM simulation, indicating a relative error of only 0.72% between the DEM simulation and experiment (Figure 12a). When the coefficient of static friction between fertilizer particles and both aluminum and ABS copolymers was set to 0.5, the relative errors of the coefficient of static friction between fertilizer particles and aluminum and acrylonitrile-butadiene-styrene copolymers were 2.44% and 3.22%, respectively (Figure 12b,c). Qualitative comparisons were also made between the simulation and bench test to analyze the coefficient of variation for each row of the fertilizer metering device. Figure 12d shows that the relative error of the coefficient of variation for each row between the simulation and bench test was less than 1.68 percentage points, indicating a good agreement between the experimental and simulation results.

4. Discussion

4.1. Fertilizer Particle and Interaction Properties Determination Method

Accurate simulation results rely on the calibration of appropriate parameters for the DEM simulation model. The coefficients of restitution, static friction, and rolling friction for both particle-particle and particle-material interactions were significant DEM simulation parameters. The coefficient of restitution is a measure of the energy transferred between particles and the surface of impact [23]. It is expressed as the ratio of the rebound velocity to the initial velocity. The coefficient of restitution between particle and plate is tested between fertilizer particle and plate using two coordinate planes and a high-speed camera in this research. Therefore, the coefficient of restitution between fertilizer particles and 304 stainless steel, ABS, and aluminum was obtained to be 0.35, 0.52, and 0.32 using bench test. However, determining the coefficient of restitution between fertilizer particles is challenging due to the uncertain movement direction and unpredictable collision. We used simulation optimization to determine the coefficient of restitution between fertilizer particles, and the coefficient of restitution between fertilizer particles was determined to be 0.323 in this simulation experiment.
The coefficients of static friction and rolling friction are important parameters that influence particle behavior. The coefficient of static friction between particles and a plate can be obtained using a platform, as demonstrated by Wang et al. [24]. However, determining the static friction between particle-particle interactions and the coefficient of rolling friction through experimental methods is challenging. A common approach to improving simulation accuracy and reducing computational requirements is to select some significant parameters. Su Y. [25] found that combinations of factors were important parameters using a Plackette–Burman experimental design. We used the Box–Behnken method to determine static friction between particle-particle interactions and the coefficient of rolling friction. We found that the coefficients of static friction and rolling friction between fertilizer particles significantly affected the angle of repose. An increase in the coefficient of rolling friction led to an increase in the angle of repose. Our findings were consistent with those from other studies (Wang et al. [16]). Song et al. [26] optimized the DEM models for fertilizer particles by training a radial basis function neural network, and the approximate model’s prediction accuracy was less than 5%. In our research, the coefficient of restitution between particles and other materials, and static friction between particles and other materials were determined by bench test. The coefficient of restitution, coefficient of static friction, and coefficient of rolling friction between fertilizer particles were determined by simulation optimization. By conducting experimental measurements in combination with simulation optimization, the interactive parameters between fertilizer particles and between fertilizer particles and other materials were determined. A DEM simulation model of fertilizer particles was constructed to accurately represent their real characteristics.

4.2. Accuracy Analysis of Simulation Parameters

The angle of repose is a commonly used parameter to evaluate the accuracy of simulation parameters [23]. The accuracy of the simulation parameters is assessed by comparing the angle of repose obtained between the simulation and experimental data, typically by measuring the accumulation of particles.
Our study indicated that the revolving cylinder method was applied to test the angle of repose in both simulation and bench tests, and the relative error between the simulation and experiment was found to be less than 3.0%. Zhou et al. [27] and Yan et al. [28] utilized the fixed funnel method to test the angle of repose for maize and soybeans. These studies confirmed the effectiveness of the angle of repose as a parameter for ensuring higher accuracy in determining particle properties. This implied that the angle of repose was an effective evaluation index to test the accuracy of DEM simulation and helped to construct the DEM model.
To further verify the accuracy and applicability of the DEM simulation parameters for fertilizer particles, the coefficient of variation for each row was measured using a fertilizer metering device in this research. The fertilizer feeding process was consistent between the simulation and bench test. The relative error of each row consistency variation coefficient was found to be less than 1.68 percentage points, and both values were below 5.0%. This verification indicated that the fertilizer DEM modeling method and the simulation parameters obtained above are accurate and can be applied in simulations of fertilizer spreaders. These results demonstrated that accurately determining simulation parameters enabled the creation of precise DEM simulation models, which were beneficial in understanding the fertilizer feeding process, fertilizer particle movement characteristics, and optimizing the structural parameters of the fertilizer feeder.

5. Conclusions

This work aimed to determine the parameters of fertilizer particles and fertilizer interactions for DEM simulation in order to ensure the accuracy of the fertilizer particle model. The simulation model has been validated by comparing it with experimental and numerical results. The conclusions were as follows:
(1) Experimentally, the triaxial size, bulk density, and density of three fertilizer varieties were measured, along with the coefficients of restitution and static friction between the fertilizers and different materials. The coefficients of restitution between the fertilizer particle and 304 stainless steel, ABS, and aluminum were 0.35, 0.52, and 0.32, respectively. The angles of repose between fertilizer particles and 304 stainless steel, ABS, and aluminum were 33.78°, 31.34°, and 30.70°, respectively.
(2) Numerically, the DEM modeling of fertilizer particles based on the triaxial size was established and the fertilizer particle model was closer to the actual fertilizer size. The angle of repose was an effective parameter to assess the accuracy of DEM simulation parameters and was found to be significantly influenced by the coefficient of static friction and the coefficient of rolling friction between fertilizer particles. The determined values of the coefficient of restitution, coefficient of static friction, and coefficient of rolling friction between fertilizer particles were 0.323, 0.381, and 0.173, respectively. The coefficient of static friction between fertilizer particles and 304 stainless steel, aluminum, and ABS was 0.453, 0.5, and 0.5, respectively.
(3) The accuracy of simulation parameters was verified by comparing the angle of repose and the coefficient of variation for each row between the DEM simulation and bench test. The relative error of the angle of repose and the coefficient of variation for each row between the simulation and experiment was less than 3.0% and 1.68 percentage points, respectively. This verified the feasibility and effectiveness of the DEM simulation parameters for fertilizer particles proposed in this research. This research can provide the accurate parameters, including fertilizer particles and particle plates, needed to construct a DEM model. The interaction between fertilizer particles and the component of fertilizer spreader, the analysis of fertilizer particle motion characteristics, and the optimization of the fertilizer spreader can be further studied.

Author Contributions

Conceptualization, X.L.; methodology, X.L., W.W. and X.D.; investigation, X.L., W.W., H.L., T.L., X.D., J.G., P.Z. and K.Y.; data curation, X.L. and W.W.; writing—original draft preparation, X.L.; writing—review and editing, X.L. and J.L.; project administration, X.L.; funding acquisition, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (31901413), Rape Industry Cluster Project of Sichuan, and the Leading the Charge with Open Competition Project of Tianfu New District (XZY1-11). Thanks for all your support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khalid, A.A.; ElKamil, T.; Ahmed, A.A.; Rangaswamy, M.; Samy, A.M.; Ahmed, M.Z.; Mohamed, K.E. Control and monitoring systems used in variable rate application of solid fertilizers: A review. J. King Saud Univ. Sci. 2023, 35, 102574. [Google Scholar] [CrossRef]
  2. Dintwa, E.; Tijskens, E.; Olieslagers, R.; Baerdenaeker, J.D.; Ramon, H. Calibration of a spinning disc spreader simulation model for accurate site-specific fertiliser application. Biosys. Eng. 2004, 88, 49–62. [Google Scholar] [CrossRef]
  3. Lv, H.; Yu, J.Q.; Fu, H. Simulation of the operation of a fertilizer spreader based on an outer groove wheel using a discrete element method. Math. Comput. Model. 2013, 58, 842–851. [Google Scholar] [CrossRef]
  4. Pasha, M.; Hare, C.; Ghadiri, M.; Gunadi, A.; Piccione, P.M. Effect of particle shape on flow in discrete element method simulation of a rotary batch seed coater. Powder Technol. 2016, 296, 29–36. [Google Scholar] [CrossRef]
  5. Yuan, J.; Liu, Q.H.; Liu, X.M.; Zhang, T.; Zhang, X.H. Granular multi-flows fertilization process simulation and tube structure optimization in nutrient proportion of variable rate fertilization. Transact. Chin. Soc. Agric. Mach. 2014, 45, 81–87. [Google Scholar]
  6. Liu, C.L.; Li, Y.N.; Song, J.N.; Ma, T.; Wang, M.M.; Wang, X.J.; Zhang, C. Performance analysis and experiment on fertilizer spreader with centrifugal swing disk based on EDEM. Transact. Chin. Soc. Agric. Eng. 2017, 33, 32–39. [Google Scholar]
  7. Wang, Y.X.; Liang, Z.J.; Cui, T.; Zhang, D.X.; Qu, Z.; Yang, L. Design and experiment of layered fertilization device for corn. Transact. Chin. Soc. Agric. Mach. 2016, 47, 163–169. [Google Scholar]
  8. Chen, G.B.; Wang, Q.G.; Xu, D.J.; Li, H.W.; He, J.; Lu, C.Y. Design and experimental research on the counter roll differential speed solid organic fertilizer crusher based on DEM. Comput. Elect. Agric. 2023, 207, 107748. [Google Scholar] [CrossRef]
  9. Coetzee, C.J.; Lombard, S.G. Discrete element method modelling of a centrifugal fertilizer spreader. Biosys. Eng. 2011, 109, 308–325. [Google Scholar] [CrossRef]
  10. Liedekerke, P.V.; Tijskens, E.; Ramon, H. Discrete element simulations of the influence of fertilizer physical properties on the spread pattern from spinning disc spreaders. Biosys. Eng. 2009, 102, 392–405. [Google Scholar] [CrossRef]
  11. Ding, S.P.; Bai, L.; Yao, Y.X.; Yue, B.; Fu, Z.L.; Zheng, Z.Q.; Huang, Y.X. Discrete element modelling (DEM) of fertilizer dual-banding with adjustable rates. Comput. Elect. Agric. 2018, 152, 32–39. [Google Scholar] [CrossRef]
  12. Xi, X.B.; Wang, R.Y.; Wang, X.T.; Shi, Y.J.; Zhao, Y.; Zhang, B.F.; Qu, J.W.; Gan, H.; Zhang, R.H. Parametric optimization and experimental verification of multi-fertilizer mixing by air blowing and blade stirring based on discrete element simulations. Comput. Elect. Agric. 2023, 210, 107895. [Google Scholar] [CrossRef]
  13. Kim, Y.S.; Siddique, M.A.A.; Kim, W.S.; Kim, Y.J.; Lee, S.D.; Lee, D.K.; Hwang, S.J.; Nam, J.S.; Park, S.U.; Lim, R.G. DEM simulation for draft force prediction of moldboard plow according to the tillage depth in cohesive soil. Comput. Elect. Agric. 2021, 189, 106368. [Google Scholar] [CrossRef]
  14. Horabik, J.; Beczek, M.; Mazur, R.; Parafiniuk, P.; Ryzak, M.; Molenda, M. Determination of the restitution coefficient of seeds and coefficients of visco-elastic Hertz contact models for DEM simulations. Biosys. Eng. 2017, 161, 106–119. [Google Scholar] [CrossRef]
  15. Huang, X.M.; Zha, X.T.; Pan, H.B.; Zong, W.Y.; Chen, H. Measurement and analysis of rapeseeds’ restitution coefficient in point-to-plate collision model. Transact. Chin. Soc. Agric. Eng. 2014, 30, 22–29. [Google Scholar]
  16. Wang, L.J.; Li, R.; Wu, B.X.; Wu, Z.C.; Ding, Z.J. Determination of the coefficient of rolling friction of an irregularly shaped maize particle group using physical experiment and simulations. Particuology 2018, 38, 185–195. [Google Scholar] [CrossRef]
  17. Liu, C.L.; Wei, D.; Song, J.N.; Li, Y.N.; Du, X.; Zhang, F.Y. Systematic study on boundary parameters of discrete element simulation of granular fertilizer. Transact. Chin. Soc. Agric. Mach. 2018, 49, 82–89. [Google Scholar]
  18. Lu, F.Y.; Ma, X.; Tan, S.Y.; Chen, L.T.; Zeng, L.C.; An, P. Simulative calibration and experiment on main contact parameters of discrete elements for rice bud seeds. Transact. Chin. Soc. Agric. Mach. 2018, 49, 93–99. [Google Scholar]
  19. Irazábal, J.; Salazar, F.; Oñate, E. Numerical modelling of granular materials with spherical discrete particles and the bounded rolling friction model. Application to railway ballast. Comput. Geotech. 2017, 85, 220–229. [Google Scholar] [CrossRef]
  20. Persson, A.S.; Frenning, G. The influence of rolling friction on the shear behaviour of non-cohesive pharmaceutical granules-An experimental and numerical investigation. Eur. J. Pharm. Sci. 2013, 49, 241–250. [Google Scholar] [CrossRef]
  21. Lei, X.L.; Yang, W.H.; Yang, L.J.; Liu, L.Y.; Liao, Q.X.; Ren, W.J. Design and experiment of seed hill-seeding centralized metering device for rapeseed. Transact. Chin. Soc. Agric. Mach. 2020, 51, 54–64. [Google Scholar]
  22. Zhu, Q.Z.; Wu, G.W.; Chen, L.P.; Zhao, C.J.; Meng, Z.J. Influences of structure parameters of straight flute wheel on fertilizing performance of fertilizer apparatus. Transact. Chin. Soc. Agric. Eng. 2018, 34, 12–20. [Google Scholar]
  23. Beakawi Al-Hashemi, H.M.; Baghabra Al-Amoudi, O.S. A review on the angle of repose of granular materials. Powder Technol. 2018, 330, 397–417. [Google Scholar] [CrossRef]
  24. Wang, S.; Yu, Z.H.; Aorigele; Zhang, W.J. Study on the modeling method of sunflower seed particles based on the discrete element method. Comput. Elect. Agric. 2022, 198, 107012. [Google Scholar] [CrossRef]
  25. Su, Y.; Xu, Y.; Cui, T.; Gao, X.J.; Xia, G.Y.; Li, Y.B.; Qiao, M.M. Determination and interpretation of bonded particle model parameters for simulation of maize kernels. Biosys. Eng. 2021, 210, 193–205. [Google Scholar] [CrossRef]
  26. Song, X.F.; Dai, F.; Zhang, F.W.; Wang, D.M.; Liu, Y.C. Calibration of DEM models for fertilizer particles based on numerical simulations and granular experiments. Comput. Elect. Agric. 2023, 204, 107507. [Google Scholar] [CrossRef]
  27. Zhou, L.; Yu, J.Q.; Wang, Y.; Yan, D.X.; Yu, Y.J. A study on the modelling method of maize-seed particles based on the discrete element method. Powder Technol. 2020, 174, 353–376. [Google Scholar] [CrossRef]
  28. Yan, D.X.; Yu, J.Q.; Wang, Y.; Zhou, L.; Yu, Y.J. A general modelling method for soybean seeds based on the discrete element method. Powder Technol. 2020, 372, 212–226. [Google Scholar] [CrossRef]
Figure 1. Triaxial size of various granular fertilizers: (a) length of fertilizer, (b) width of fertilizer, and (c) thickness of fertilizer.
Figure 1. Triaxial size of various granular fertilizers: (a) length of fertilizer, (b) width of fertilizer, and (c) thickness of fertilizer.
Agriculture 13 01704 g001
Figure 2. Test principle of collision recovery coefficient.
Figure 2. Test principle of collision recovery coefficient.
Agriculture 13 01704 g002
Figure 3. Measurement method of particle velocity: (a) particle rebound test platform, (b) test system, (c) particle movement images, (d) particle movement trajectory, (e) tracking particle velocity, (1) mirror, (2) first coordinate plane, (3) second coordinate plane, (4) base plate, (5) vacuum nozzle, (6) particle rebound platform, (7) test plate, and (8) high-speed camera.
Figure 3. Measurement method of particle velocity: (a) particle rebound test platform, (b) test system, (c) particle movement images, (d) particle movement trajectory, (e) tracking particle velocity, (1) mirror, (2) first coordinate plane, (3) second coordinate plane, (4) base plate, (5) vacuum nozzle, (6) particle rebound platform, (7) test plate, and (8) high-speed camera.
Agriculture 13 01704 g003
Figure 4. Programming process of extracting trajectory.
Figure 4. Programming process of extracting trajectory.
Agriculture 13 01704 g004
Figure 5. Multifunctional comprehensive test platform.
Figure 5. Multifunctional comprehensive test platform.
Agriculture 13 01704 g005
Figure 6. Particle-contact model.
Figure 6. Particle-contact model.
Agriculture 13 01704 g006
Figure 7. 3D model parameters of fertilizer: (a) 3D model; (b) front view; and (c) lateral view.
Figure 7. 3D model parameters of fertilizer: (a) 3D model; (b) front view; and (c) lateral view.
Agriculture 13 01704 g007
Figure 8. DEM simulation of the hollow cylinder test method.
Figure 8. DEM simulation of the hollow cylinder test method.
Agriculture 13 01704 g008
Figure 9. Angle of repose setup: (1) high-speed camera, (2) PC, (3) universal material testing machinery, (4) controller, (5) cylinder, and (6) lifting frame.
Figure 9. Angle of repose setup: (1) high-speed camera, (2) PC, (3) universal material testing machinery, (4) controller, (5) cylinder, and (6) lifting frame.
Agriculture 13 01704 g009
Figure 10. Image processing of repose angle: (a) original image, (b) gray processing, (c) binary processing, (d) fitting line of repose angle.
Figure 10. Image processing of repose angle: (a) original image, (b) gray processing, (c) binary processing, (d) fitting line of repose angle.
Agriculture 13 01704 g010
Figure 11. Fertilization performance setup. (a) experiment setup and (b) simulation setup, (1) fertilizer box, (2) fertilizer metering mechanism, (3) fertilizer tube.
Figure 11. Fertilization performance setup. (a) experiment setup and (b) simulation setup, (1) fertilizer box, (2) fertilizer metering mechanism, (3) fertilizer tube.
Agriculture 13 01704 g011
Figure 12. Simulation parameters verification results. (a) Angle of repose verification between simulation and experiment, (b) effects of coefficient of static friction between fertilizer particle and aluminum on angle of repose, (c) effects of coefficient of static friction between fertilizer particle and ABS on the angle of repose, and (d) coefficient of variation for each row verification between simulation and experiment. Note: The light red zone represents the error bar, and the green zone represents the experimental value.
Figure 12. Simulation parameters verification results. (a) Angle of repose verification between simulation and experiment, (b) effects of coefficient of static friction between fertilizer particle and aluminum on angle of repose, (c) effects of coefficient of static friction between fertilizer particle and ABS on the angle of repose, and (d) coefficient of variation for each row verification between simulation and experiment. Note: The light red zone represents the error bar, and the green zone represents the experimental value.
Agriculture 13 01704 g012
Table 1. Triaxial size of various types of granular compound fertilizer.
Table 1. Triaxial size of various types of granular compound fertilizer.
Fertilizer TypesNutrient Ratio (%)Length (mm)Width (mm)Thickness (mm)Equivalent Diameter (mm)Sphericity (%)Bulk Density (kg m−3)
ZhushangN:P2O5:K2O = 20:9:114.24 ± 0.403.77 ± 0.243.35 ± 0.293.76 ± 0.2189.21 ± 5.67941.90
ZhonghuaN:P2O5:K2O = 15:15:154.52 ± 0.464.04 ± 0.363.63 ± 0.354.04 ± 0.3189.68 ± 5.42799.00
SanningN:P2O5:K2O = 15:15:153.55 ± 0.253.43 ± 0.223.32 ± 0.223.43 ± 0.2296.80 ± 2.04826.00
Average-4.103.753.433.7491.90865.97
Table 2. Parameters used in the DEM simulations.
Table 2. Parameters used in the DEM simulations.
ParametersFertilizer304 Stainless SteelABSAluminumPlasticSource/Method of Measurement
Poisson’s ratio0.250.30.3940.360.47[14]
Shear modulus (Pa)2.8 × 1077.9 × 10108.9 × 1086.85 × 10102.9 × 109[21]
Density (kg m−3)15427 9301 06027071350Measured for fertilizer and [21] for other materials
Coefficient of restitution
Particle-particle0.3–0.4 *----
Particle-material-0.350.520.320.45Measured
Coefficient of static friction
Particle-particle0.30–0.40 *----
Particle-material-0.45–0.55 *0.45–0.55 *0.45–0.55 *0.3Measured
Coefficient of rolling friction
Particle-particle0.10–0.20----
Particle-material-0.10.10.10.1[22]
Gravitational acceleration (m s−2)9.81
Note: * represents the uncertain DEM simulation values.
Table 3. Coding of factors.
Table 3. Coding of factors.
Parameter−101
Coefficient of restitution of fertilizer-fertilizer (x1)0.300.350.40
Coefficient of static friction of fertilizer-fertilizer (x2)0.300.350.40
Coefficient of static friction of fertilizer-plate material (x3)0.450.500.55
Coefficient of rolling friction of fertilizer- fertilizer (x4)0.100.150.20
Table 4. Coefficient of restitution for granular compound fertilizer against different types of materials.
Table 4. Coefficient of restitution for granular compound fertilizer against different types of materials.
Fertilizer TypesType of Contact SurfaceDropping Height (mm)
100200300
ZhushangAcrylonitrile-butadiene-styrene copolymer0.5240.5200.518
304 stainless steel0.3810.3790.364
Aluminum0.3800.3300.326
ZhonghuaAcrylonitrile-butadiene-styrene copolymer0.5330.5230.522
304 stainless steel0.3500.3460.330
Aluminum0.3310.3220.311
SanningAcrylonitrile-butadiene-styrene copolymer0.5230.5170.490
304 stainless steel0.3760.3750.357
Aluminum0.3460.3210.312
AverageAcrylonitrile-butadiene-styrene copolymer0.5270.5200.510
304 stainless steel0.3690.3670.350
Aluminum0.3520.3240.316
F valueFertilizer types0.088
Material type305.62 **
Dropping height0.613
Note: ** represent significantly different at the 0.01 probability level.
Table 5. Angle of repose between fertilizer particle and plate material.
Table 5. Angle of repose between fertilizer particle and plate material.
Type of Contact SurfaceAngle of Repose (°)
Test 1Test 2Test 3Test 4Test 5Average
304 stainless steel32.8033.2432.0434.1636.6833.78
Acrylonitrile-butadiene-styrene copolymer31.4031.1531.0332.9330.1831.34
Aluminum31.0228.4733.7329.530.4430.70
Table 6. Simulation scheme and results between fertilizer particle and 304 stainless steel.
Table 6. Simulation scheme and results between fertilizer particle and 304 stainless steel.
Test NumberTest FactorsAngle of Repose (°)
x1x2x3x4
1−110031.75 ± 3.09
20−10−121.73 ± 3.08
3−100133.69 ± 1.49
400−1−124.74 ± 3.16
5011031.92 ± 0.98
61−10027.25 ± 2.37
70−10134.21 ± 0.96
8000031.00 ± 3.11
9−1−10028.00 ± 1.36
100−11027.14 ± 2.06
1101−1032.52 ± 1.97
12010−123.45 ± 2.18
13001−124.15 ± 2.76
1410−1030.91 ± 0.87
15−10−1030.58 ± 1.00
1600−1137.40 ± 4.37
17000030.69 ± 1.20
18−101029.99 ± 1.16
19000031.54 ± 2.08
200−1−1026.46 ± 2.99
21010133.52 ± 0.81
22001134.28 ± 1.31
23100135.53 ± 1.75
24101028.38 ± 2.04
25110031.03 ± 2.10
26000031.03 ± 0.27
27100−125.40 ± 0.43
28−100−124.39 ± 1.87
29000029.99 ± 1.49
Table 7. ANOVA of quadratic model.
Table 7. ANOVA of quadratic model.
SourceSum of SquaresdfMean SquareF-Valuep-Value
Model402.231040.2229.98<0.0001 **
x10.000810.00080.00060.9814
x231.41131.4123.420.0001 **
x33.7713.772.810.1109
x4349.461349.46260.5<0.0001 **
x1x30.943310.94330.70320.4127
x2x41.4511.451.080.3119
x3x41.6111.611.20.2882
x2x210.62110.627.920.0115 *
x3x30.034610.03460.02580.8742
x4x44.5114.513.360.0834
Residual24.15181.34
Lack of fit22.85141.635.050.0648
Pure error1.2940.3234
Total426.3828
Note: * and ** represent significantly different at the 0.05 and 0.01 probability level, respectively.
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

Lei, X.; Wu, W.; Deng, X.; Li, T.; Liu, H.; Guo, J.; Li, J.; Zhu, P.; Yang, K. Determination of Material and Interaction Properties of Granular Fertilizer Particles Using DEM Simulation and Bench Testing. Agriculture 2023, 13, 1704. https://doi.org/10.3390/agriculture13091704

AMA Style

Lei X, Wu W, Deng X, Li T, Liu H, Guo J, Li J, Zhu P, Yang K. Determination of Material and Interaction Properties of Granular Fertilizer Particles Using DEM Simulation and Bench Testing. Agriculture. 2023; 13(9):1704. https://doi.org/10.3390/agriculture13091704

Chicago/Turabian Style

Lei, Xiaolong, Wencheng Wu, Xuan Deng, Tao Li, Hongnan Liu, Jinyue Guo, Ju Li, Peixu Zhu, and Ke Yang. 2023. "Determination of Material and Interaction Properties of Granular Fertilizer Particles Using DEM Simulation and Bench Testing" Agriculture 13, no. 9: 1704. https://doi.org/10.3390/agriculture13091704

APA Style

Lei, X., Wu, W., Deng, X., Li, T., Liu, H., Guo, J., Li, J., Zhu, P., & Yang, K. (2023). Determination of Material and Interaction Properties of Granular Fertilizer Particles Using DEM Simulation and Bench Testing. Agriculture, 13(9), 1704. https://doi.org/10.3390/agriculture13091704

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