A Calibration Method for Contact Parameters of Maize Kernels Based on the Discrete Element Method

: Clarifying the maize kernel movement during the crushing process is critical for improving the design and optimization of the impact mill. Rather than through experiments, maize kernel movement can be quantitatively analyzed through the discrete element method (DEM)


Introduction 1.Maize Kernel Movement
Maize kernels are one of the most commonly used cereals in poultry diets worldwide [1].Crushing is an important process in the journey from maize kernel to feed product.The maize kernel movement plays an important role in the crushing process [2].Understanding the maize kernel flow behavior during the crushing process is necessary for the design and optimization of the impact mill [3].Due to the complex working environment and limited available sensors in the crushing processes, it is difficult to quantify the particle movement using traditional experimental methods [4].Simultaneously, the experimental method for particle motion analysis necessitates the use of expensive experimental equipment, such as a high-speed camera, a laser, and divergent light sheet optics [5], which increases the experimental expense.

Discrete Element Method
The discrete element method (DEM) has become a common method used in the study of bulk material crushing systems.It can be employed instead of an experiment to quantitatively analyze the movement and breakage of particles [6,7].The Ab-T10 breakage model in Rocky DEM [8] can better predict the size distribution and shape of progeny particles in crushing simulations, which is helpful for in-depth research on the crushing mechanism and optimization of the mechanical structure [8].However, to assure the accuracy of the DEM simulation, an accurate particle model and input parameters, including intrinsic and contact parameters, are required.Unfortunately, there are no standard methods for the construction of the particle model and the calibration of contact parameters [9].
In recent years, the construction methods of the particle model have been studied.The approximate representation of real granular material as regular bodies is one of the most common methods [10][11][12].This is suitable for granular materials with simple structures, such as potatoes and rapeseed grains.Another method is to obtain the contours of the real granular material by 3D scanning, CT, and other technologies, and then fill the contours with spherical particles to complete the particle establishment [6,13,14].These particles still maintain the structural characteristics of the curved surface, such as in maize kernels and barley grains [15].However, the breakage models of Rocky DEM only work with convex polyhedral particles for the crushing process simulation, and no spherical, rounded, or concave shape particles are able to be used [16].Therefore, the above methods are not suitable for particle construction in Rocky DEM, especially for agricultural granular materials with complex structures, such as maize kernels.
One critical step when developing a DEM model is determining the model input parameters, including the intrinsic parameters of granular materials (such as bulk density, Young's modulus, and Poisson's ratio) and the contact parameters (mainly the restitution coefficients and friction coefficients of p-w and p-p).The intrinsic parameters can be measured directly or values reported in the literature can be used [6].However, it is necessary to calibrate the contact parameters in the simulation because of differences in geometry between particles and real granular materials [9].The method of combining an experiment with a simulation is commonly used to calibrate the contact parameters [9].
The restitution coefficient affects the movement of particles after they collide with the contact surface in the simulation [17,18].Restitution coefficient calibration methods have been investigated for agricultural granular materials, such as Panax notoginseng seeds, rapeseed, and wheat.One of the common restitution coefficient calibration methods is the trajectory method.High-speed cameras are widely used with the trajectory method to record the trajectories of grains during collisions and calibrate the restitution coefficients of p-w and p-p by changing the parameter values many times until the simulated particle trajectory is consistent with the measured grain trajectory [17,[19][20][21].The high-speed camera is an indispensable part of this method, which increases the cost of the experiment.Additionally, the acquisition and analysis process of the grains' trajectory is complex, which reduces the calibration efficiency.Therefore, in order to improve the efficiency of the restitution coefficient calibration, the use of a self-made experiment device has become a more effective approach [6,22].For example, the inclined surface drop method was employed to calibrate the restitution coefficients of potato-contact surface [22], maize kernel-contact surface, and wheat-contact surface [6].The inclined surface drop method avoids the use of a high-speed camera, which reduces the cost, simplifies the process, and effectively improves the calibration efficiency.Nevertheless, at present, the inclined surface drop method is basically only used to calibrate the restitution coefficient between particles and the contact surface, while the calibration of the p-p restitution coefficient is rarely studied.Friction coefficients (static friction and dynamic friction coefficients of p-w and p-p) are the key parameters affecting particle flow characteristics [21].At present, the funnel method is commonly used to calibrate friction coefficients, but only a single index like the angle of repose (AoR) or discharging time is used to obtain and calibrate the friction coefficient [6,14,17,20,[23][24][25].In fact, calibrating the friction coefficients with a single index only will reduce the accuracy of parameter calibration [26].
In conclusion, to improve the accuracy of maize kernel crushing process simulation, exploring a maize kernel model construction method suitable for the Ab-T10 breakage model and optimizing the current contact parameter calibration method to improve the accuracy of maize kernel model parameter calibration is critical.

Aim of the Paper
In order to improve the accuracy of the maize kernel crushing simulation, in this study, we propose a method for maize kernel model construction based on the section method and a method for the calibration of contact parameters involving experiments and simulations using Rocky DEM.The inclined surface drop method device was optimized, and the restitution coefficient for both p-w and p-p were calibrated at the same time.In addition, both the AoR and the discharging time were used as indexes to calibrate the dynamic friction coefficients of the maize kernel by the funnel method [17].By comparing the simulation and experimental results, the feasibility and accuracy of the maize kernel model and contact parameter values were verified.

Experiment Material
Maize kernels are an important raw feed material.According to different feed formulas, the wet basis moisture content of maize kernels in the crushing process ranges from 10% to 17%.Generally, it is best to crush maize kernels with a moisture content of 13% [27].The maize kernel variety "Xianyu-335," widely planted in the north and northeast of China, was chosen as the experiment material.Three samples of maize kernels were prepared by following the standard method [28].The wet basis moisture content of three samples maize kernels was 10.9% ± 0.8 (sample 1), 13.9% ± 0.7 (sample 2), and 17.2% ± 1.2 (sample 3), respectively.The bulk density of the maize kernels in the three samples was 782.57± 3.78 g/L, 769.06 ± 1.48 g/L, and 739.00 ± 6.69 g/L, respectively.The wet basis moisture content and bulk density were analyzed by following the standard experiment/formulae [29,30] and are presented as the mean ± standard deviation.

Physical Experiment 2.2.1. Inclined Surface Drop Method
In order to improve the calibration efficiency and accuracy of the restitution coefficients of p-w and p-p, this study optimized and improved the experiment device of the inclined surface drop method [6].As shown in Figure 1a, the improved experiment device consists of a particle releasing pipe, a PVC division plate, a particle collection box, and two impact plates, one of which is a low-carbon steel plate (impact plate 1), and the other is a lowcarbon steel plate covered with maize kernels (impact plate 2).The geometric dimensions of the experiment device were determined by pre-experiment measures and are shown in Figure 1b.In particular, the spacing of the PVC division plates was set to 30 mm, as suggested by [6].
The particle releasing pipe, the impact plate 2, and quartz sand were added to the improved experiment device.The particle releasing pipe can guarantee that the contact positions of each particle are close to the impact plate, avoiding the influence of the contact position on the particle landing positions and preventing particles from falling out of the particle collection box after contacting the edge of the impact plate.When the particles fall down from the releasing pipe, the result of the collision on impact plate 1 can be used to calibrate the restitution coefficient of p-w.Accordingly, the result of the collision on impact plate 2 can be used to calibrate the restitution coefficient of p-p.As shown in Figure 1a, impact plate 2 is covered with maize kernels.Impact plate 2 was prepared as follows: firstly, a number of maize kernels of uniform particle size and thickness were selected, and then their tip caps were cut off.Secondly, the treated maize kernels were uniformly pasted to the surface of a steel plate with double-sided adhesive.Finally, the steel plate surface with maize kernels was placed downwards on a flat table, and a 5 kg weight was pressed onto the surface of the steel plate and checked every 20 min until the maize kernels were completely fixed and the surface with maize kernels was flat.Subsequently, impact plate 2 was completed.The bottom plate of the particle collection box was uniformly paved with quartz sand to prevent the particles rebounding during the experiment and affecting the accuracy of the experimental results.For each sample of maize kernels, 300 maize kernels [6] in random orientations were put into the release position and released one-by-one, and then the particles collided with the impact plate.After bouncing off of the impact plate, the maize kernels were collected by the collection box and fell into separate bins, depending on their moving distance after impact.The experiment process is shown in Figure 1b.The number of maize kernels in each bin was recorded.The average bin number (i) was defined as shown in Equation (1).Each experiment was repeated three times to obtain the average bin number in the inclined surface drop experiment.The statistical results are shown in Table 1.
Here, i is the average bin number, i is the bin number, n is the total number of bins, N i denotes the number of kernels that fall into each bin i, and N total represents the total number of kernels used in the experiment.

Funnel Method
The funnel method device consists of a funnel, a baffle, a base plate, and a fixed support, as shown in Figure 2a.The geometric dimensions of the experimental device were determined by the pre-experiment measures and are shown in Figure 2b.The pre-experiment results show that when the mass of maize kernels was 0.45 kg, a stable AoR of maize kernels could be formed, and the AoR would not change as the mass increased.Therefore, 0.45 kg of maize kernels were released into the funnel evenly in the funnel experiment.After standing for 5-10 s, maize kernels in the funnel were basically stationary, and the baffle was quickly pulled out, and the maize kernels to fall evenly on the base plate.Subsequently, the maize kernels gradually formed an AoR, as shown in Figure 2c.In addition, the video of the falling process of the maize kernels was analyzed in a frame-by-frame manner by PotPlayer, and then the discharging time was determined with an accuracy of 0.001 s.
Each group of experiments was repeated 3 times to obtain the mean AoR value and the discharging time, and the results are listed in Table 2.

Simulation Set-Up
Fundamentally, the DEM solves Newton's equations of motion to resolve particle motion and uses a contact law to describe the inter particle contact forces.Forces are typically integrated explicitly in time to predict the time history response of the material using an appropriate numerical integration method.The hysteretic linear spring model for normal force calculation and the linear spring Coulomb limit model for tangential force calculation were used to simulate the inclined surface drop test and the funnel test in Rocky DEM [16].

Construction of the Maize Kernel Model
Maize kernels were divided into wedge, ellipsoid, and irregular shapes, among which wedge-shaped maize kernels accounted for the highest proportion of kernels, about 94.72% [31].The effect of the moisture content on the geometric structure of maize kernels is negligible [32].Therefore, this study took the wedge-shaped maize kernel from sample 2 as the object to establish the particle model for the DEM simulation.For granular materials with complex geometry, the geometry can be characterized by the sizes of three axes, length (L), width (W), and thickness (T) [33], as shown in Figure 3.In the crushing simulation, the resulting fragments preserve both the mass and volume when a particle breaks [16].Hence, the sizes of the three axes and particle volume (V) should be indicators of the accuracy of the particle model.In this study, 20 maize kernels [6] were selected from sample 2, and then the sizes of the three axes and the volume information were obtained by 3D scanning.L was 12.56 ± 0.59 mm, W was 8.52 ± 0.32 mm, T was 4.42 ± 0.19 mm, and V was 293.52 ± 25.24 mm 3 .
The 3D-scanned maize kernel (Figure 4a) was imported into 3D software.The central section plane was established based on the centroid of the maize kernel and was arranged so that it was equally spaced along the positive and negative directions of the Z axis.Combining the thickness of the maize kernel and the particle construction method [33], the section plane array distance was set to 0.30 mm, 0.40 mm, 0.50 mm, 0.60 mm, 0.70 mm, 0.80 mm, 0.90 mm, and 1.00 mm, respectively, to explore the influence of the array distance on the accuracy of the maize kernel model and to obtain the optimal array distance for use in the maize kernel model.The maize kernel model construction process was as follows.Firstly, a Boolean operation was applied between the 3D scanned maize kernel and each section plane to obtain the contour lines of the 3D scanned maize kernel (Figure 4b).Secondly, multi-segment straight lines were used to fit the contour line (Figure 4c).This step was done to ensure that the maize kernel model was a convex polyhedral particle.Finally, the multi-segment straight lines were filled to form a surface, and the adjacent surfaces were blended to obtain the maize kernel model (Figure 4d).As shown in Figure 5, the relative errors of the sizes of the three axes and volume and the number of maize kernel model faces were all affected by the array distance.The relative errors of the sizes of the three axes and volume between the maize kernel model Firstly, a Boolean operation was applied between the 3D scanned maize kernel and each section plane to obtain the contour lines of the 3D scanned maize kernel (Figure 4b).Secondly, multi-segment straight lines were used to fit the contour line (Figure 4c).This step was done to ensure that the maize kernel model was a convex polyhedral particle.Finally, the multi-segment straight lines were filled to form a surface, and the adjacent surfaces were blended to obtain the maize kernel model (Figure 4d).
As shown in Figure 5, the relative errors of the sizes of the three axes and volume and the number of maize kernel model faces were all affected by the array distance.The relative errors of the sizes of the three axes and volume between the maize kernel model constructed with an array distance of 0.8 mm and the real maize kernel were maximal, because the section plane could not capture the contours of the maize kernel edge structure under this array distance, resulting in the geometric structure of the constructed maize kernel model being quite different from that of the real maize kernel.When the array distance was less than 0.5 mm, the relative errors of the sizes of the three axes and the volume between the maize kernel model and the real maize kernel were lower and tended to be stable.However, the number of maize kernel model faces increased significantly as the array distance decreased.The simulation efficiency decreased as the number of particle faces increased when the number of particles was the same in the DEM simulation [16].Therefore, the maize kernel model with an array distance of 0.5 mm was used.This model is similar to the geometric structure of the maize kernel, and the number of particle faces is relatively low.In other words, this maize kernel model can guarantee both the accuracy and efficiency of the simulation.As shown in Figure 5, the relative errors of the sizes of the three axes and volume between the maize kernel model and the real maize kernel were 3.81% and 0.48%, respectively.In a discrete element simulation, Poisson's ratio and Young's modulus were used to define the physical properties of particles.The results of a study performed by Ma Wenpeng et al. showed that the particle moisture content has no significant effect on Poisson's ratio [34], so the value of Poisson's ratio in the three samples was set to 0.4 [35,36].The Young's modulus of maize kernels with different moisture contents can be calculated with Equation (2) [37]: where Y is Young's modulus (Mp), m is the moisture content (as a percentage, wet basis), and Z is a dummy variable (Z = 0 for dent kernels and Z = 1 for flint kernels).Since the Xianyu-335 maize kernel is a partial flint kernel [38], Z is equal to 1.
The experimental device models and maize kernel models were imported into Rocky DEM.Table 3 summarizes the parameters adopted to define the material properties of the particles and the low-carbon steel plate.The geometry of the experimental device duplicated that used in the simulation of the inclined surface drop, and the schematic diagram is shown in Figure 6.A virtual circular surface consistent with the particle releasing pipe in terms of diameter (20 mm) was established to generate particles in Rocky DEM.The particles were generated continuously, and the generation rate was 0.005 kg/s, the generation time was 23 s, and the total number of particles was 300.A flat plate was used to replace the bottom plate of the quartz sand to improve the simulation efficiency, and the restitution coefficient from particle to bottom plate was set to the minimum value of 0.1.The simulation duration and output frequency were 30 s and 0.01 s, respectively.The pre-simulation results showed that the particle would not rebound after falling on the bottom plate, indicating that the parameters were reasonable.After the simulation finished, the number of maize kernels in each bin was recorded, and the average bin number was calculated using Equation (1).The relationship between the restitution coefficient and the average bin number was explored through the single factor test, and the regression equation between the average bin number and the restitution coefficient (C) was established using Matlab.The average bin number of measured values was substituted into the regression equation to determine the restitution coefficient.In the calibration process, the particle to particle and particle to low carbon steel restitution coefficients were both set to 0.10-0.90, with an increment of 0.05 [1].

Calibration of the Friction Coefficients
As shown in Figure 7a, the geometry and testing conditions used in the experimental test setup were duplicated in the DEM simulations.A virtual circular surface consistent with the upper end funnel in terms of diameter (100 mm) was established to generate particles in the Rocky DEM.The particles were generated continuously, and the generation rate was 0.45 kg/s, the generation time was 1 s, and the total mass of particles was 0.45 kg.The simulation duration and output frequency were 3 s and 0.001 s, respectively.After finishing the simulation, the AoR of the maize kernel models shown in Figure 7b and the discharging process were recorded in Rocky DEM.The contact parameters were demonstrated to affect the AoR and discharging time.Specifically, these were the particle to particle static friction coefficient (pps), the dynamic friction coefficient (ppd), the restitution coefficient (C P ), the particle to low carbon steel static friction coefficient (pbs), the dynamic friction coefficient (pbd), and the restitution coefficient (C B ).The two-level factorial design, the response surface methodology, and the numerical optimization in Design Expert 12.0 were used to calibrate the friction coefficients.The two-level factorial design was used to screen the major parameters that significantly influence the AoR and discharging time.In this test design, the above mentioned factors were set to two levels.Referring to the relevant literature [1,9,31], the threshold and levels of the contact parameters were designed as shown in Table 4.The response surface methodology based on central composite design (CCD) was used to analyze the influence law of major factors on indexes, and the numerical optimization was used to determine the optimal values of major factors to bring the indexes closer to the target value.

Restitution Coefficient Calibration
Figures 8 and 9 show the average bin number of simulation results varied with the restitution coefficients of p-p and p-w.The average bin number increased as the restitution coefficient increased, which corresponds with Kingsly Ambrose's conclusion [6], indicating that the simulation was reliable.Furthermore, the average bin number of each sample of maize kernels differed, although the restitution coefficient was the same.In order to ensure the accuracy of the restitution coefficient when applying the regression equation, the simulation results of each sample of maize kernels were used to perform a regression analysis by Matlab to determine the regression equation of the average bin number and the restitution coefficient, as shown in Figures 8 and 9.The adjustment R 2 values of the regression equations were all higher than 0.99, indicating an excellent fitting effect.In the regression equations, A P-B-n represents the average bin number when the maize kernels of the n th sample collided with the low carbon steel plate, where n = 1, 2, 3, with the units being pieces.C B-n is the restitution coefficient of maize kernels of the nth sample contacting the low-carbon steel plate.A P-P-n is the average bin number of collisions between maize kernels in the nth sample.C P-n is the restitution coefficient between maize kernels of the nth sample.
The average bin numbers obtained from the experiments in Section 1.1 were substituted into the regression equations to determine the restitution coefficients, as shown in Table 5.The restitution coefficient decreased as the moisture content of maize kernels increased, and the particle to low-carbon steel plate restitution coefficient was higher than that of the particle to particle restitution coefficient.The results show that the energy loss during the collision between maize kernels and the low-carbon steel plate was lower than that during the collision between maize kernels, which corresponds to Yan Hui's conclusion [39], indicating that the simulation was reliable.The two-level factorial test scheme and simulation results for Sample 2 are shown in Tables 6 and 7 shows the results of the variance analysis on the simulation results as well as the significance results for each simulation parameter.It can be seen from Table 7 that the ppd and pbd (p < 0.01) have very significant influences on the AoR and discharging time of the simulation, and the effects of the other parameters (p > 0.05) on the AoR and discharging time of the simulation results were found to be very small.Therefore, only the significant parameters of the ppd and pbd were considered in the subsequent response surface experiments.The restitution coefficients of non-significant factors were set in accordance with Table 5.The contact parameters were calibrated for a future crushing simulation of the hammer mill.Since the static friction coefficients have little effect on the simulation accuracy of the crushing process, the pps and pbs were set to 0.31 [40] and 0.40 [39], respectively, by referring to the literature.Since the results of the two-level factor analysis for all samples of maize kernels were consistent, only the simulation results and analysis of variance results of sample 2 are listed in this paper.

Response Surface Methodology Analysis
According to the results of the two-level factorial test, the ppd and pbd were taken as factors for the response surface methodology test.A two factor and three level response surface methodology test was designed.The test scheme and simulation results are shown in Tables 8 and 9 shows the variance of AoR and the analysis of the discharging time of the response surface linear model.It can be seen from Table 9 that the linear model p-value was less than 0.01, which proves that the model is reliable.ppd and pbd were shown to have significant impacts on the AoR and discharging time, and a subsequent analysis of the ppd and pbd was carried out.The contours of the AoR and discharging time were determined as functions of ppd and pbd for the three samples of maize kernels based on the results of the response surface methodology analysis and are shown in Figure 10.As with increases the values of ppd and pbd, the AoR and discharging time gradually increased.Because larger ppd and pbd values are associated with greater relative movement resistance between particles and the contact surface, the AoR and the discharging time will increase accordingly.
According to Table 2, contour lines were drawn using the AoR and discharging time from the experiment results and are shown as red lines in Figure 10 when a single index (AoR or discharging time) was used for parameter calibration (for example, when calibrating the parameters of ppd and pbd, only the AoR was selected as the index).It can be seen from Figure 10 that the combinations of ppd and pbd consistent with the points on the contour lines of the AoR will be the results of parameter calibration, and it is obvious that there can be many groups of calibrated ppd and pbd combinations.However, the ppd and pbd of each granular material should be unique both in simulation and in reality.Therefore, the accuracy of parameter calibration will be reduced when only a single response value is used to calibrate ppd and pbd.This conclusion is consistent with Wensrich C.M's research results [26].Based on the above analysis, the AoR and discharging time were considered together to calibrate the ppd and pbd.The numerical optimization was employed to solve ppd and pbd.The solution process involved taking the experiment results for the AoR and discharging time of each sample of maize kernels and the thresholds of ppd and pbd as the constraints and then solving the values of ppd and pbd using the numerical optimization.Table 10 shows the solutions to ppd and pbd, and a unique solution can be obtained for ppd and pbd for each sample of maize kernels.Meanwhile, the desirability of these solutions was calculated to be 1, indicating that the calibrating method is reliable.The results show that it is feasible to calibrate the ppd and pbd with the AoR and discharging time as indexes using the funnel method, and this can improve the accuracy and reliability of parameter calibration.

Test Validation
The accuracy of the DEM simulation, measured as the closeness of DEM results to experimental values, should be evaluated by an accredited criterion.The relative error was used as a criterion for the assessment of a simulation's accuracy [26].The relative error measures the accuracy of a DEM simulation using the difference between the results of a simulation (S) and an experiment (E).The relative error was calculated with Equation (3): where P error is the relative error, and E and S are the measured and simulated values, respectively.

Verification of Restitution Coefficients
The accuracy and reliability of the calibrated restitution coefficients should be verified with a new series of experimental results.For this purpose, the experimental results of the inclined surface drop test for 300, 400, and 500 maize kernels were compared with the results of the calibrated DEM simulation [26].Table 11 illustrates the comparisons between the experimental and simulation results of the average bin number for different amounts of maize kernels.The maximal relative errors of the average bin number of the particle to low-carbon steel plate and particle to particle were, respectively, 2.63% and 4.38%, which verifies the accuracy of the DEM simulation and the reliability of the calibrated method [6].

Verification of Dynamic Friction Coefficients
The accuracy and reliability of the calibrated dynamic friction coefficients should be verified by a new series of experimental results.For this purpose, the experimental results of the funnel test for maize kernel quantities of 0.45 kg, 0.60 kg, and 0.90 kg were compared with the results of the calibrated DEM simulation [26].Table 12 shows comparisons between the experimental and simulation results for the AoR and the discharging time for different amounts of maize kernels.The maximal relative errors of the AoR and discharging time were, respectively, 5.49% and 6.98%, which also verifies the accuracy of the DEM results and the reliability of the calibrated method.

Limitations of the Study
In the inclined surface drop experiment and the funnel experiments, maize kernels with similar structures and sizes were selected for the experiments, but maize kernel models of unified structure and size were used to complete the simulation in this study.Because the flow characteristics of the particles are impacted by their size and shape, as well as the percentage of particles of varied sizes, this could have compromised the accuracy of the parameter calibration.Therefore, in future scientific research, the particle size and shape proportion of maize kernels in the experiment should be counted, followed by the determination of the particle size distributions and the conduction of maize kernel models of various sizes.The number of particle with various sizes will be defined in simulation according to actual particle size distributions, bringing the simulation closer to actual experiments and increasing the practical significance of the simulation.

Conclusions
(1) The maize kernel model was constructed by combining 3D scanning and the section method.The relative errors of the sizes of the three axes and the volume between the particle model and the real maize kernel were 3.81% and 0.48%, respectively.(2) The restitution coefficients between the maize kernel and low-carbon steel plate and between maize kernels were calibrated by improving the inclined surface drop method.Moreover, the dynamic friction coefficients between the maize kernel and low-carbon steel plate and between maize kernels were calibrated using the funnel method by considering the AoR with the discharging time.
(3) Validation tests were carried out, and the maximum relative errors between the simulated and measured values obtained with the inclined surface drop method and funnel method were 4.38% and 6.98%, respectively, which verified the feasibility and accuracy of the particle model construction and parameter calibration method.
In this study, we proposed a method for constructing a convex polyhedral particle model suitable for Rocky DEM and a calibration method to calibrate contact parameters, and then the research methods were verified which are feasible and accurate.The research method used in this study could provide support for the construction of the particle model and calibration of the contact parameters of agricultural granular materials with complex geometric structures, as well as the maize kernel model and the values of the contact parameters in this paper can provide a reference for the DEM simulation of the maize kernel crushing process to design and optimize of the impact mill.

Figure 3 .
Figure 3. Sizes of the three axes of the maize kernel.
Agriculture 2022, 12, x FOR PEER REVIEW 7 of 18 for use in the maize kernel model.The maize kernel model construction process was as follows.

Figure 5 .
Figure 5. Maize kernel models with different array distances.

Figure 6 .
Figure 6.Simulation setup used for the inclined surface drop method.

Figure 7 .
Figure 7. Simulation setup used for the funnel method.(a).Simulation device; (b).AoR of the maize kernel model.

Figure 8 .
Figure 8. Relationship between the average bin number and the restitution coefficient of particle to low-carbon steel restitution coefficient.

Figure 9 .
Figure 9. Relationship between the average bin number and the particle to particle restitution coefficient.

Figure 10 .
Figure 10.Effects of dynamic friction coefficients on the AoR and discharging time of maize kernels.(a) Sample 1; (b) Sample 2; and (c) Sample 3 of the maize kernel models.

Table 1 .
Average bin number of maize kernels.

Table 2 .
AoR and discharging time of the maize kernels.

Table 3 .
Parameters used in the simulation model for maize kernel particles.

Table 4 .
Factors and levels of the two-level factor design.

Table 5 .
The results for the maize kernel restitution coefficients.

Table 6 .
The scheme and results of the two-level factor test.

Table 7 .
Parameter significance determined by the analysis of variance.

Table 8 .
The scheme and results related to the response surface design.

Table 9 .
The analysis of variance of the AoR and discharging time for the response surface linear model.

Table 10 .
Constraints and solutions for dynamic friction coefficients.

Table 11 .
Relative error between simulation and experiment values of the average bin number for maize kernels.

Table 12 .
Verification of the calibration of the AoR and discharging time.