Next Article in Journal
MSSFNet: A Multiscale Spatial–Spectral Fusion Network for Extracting Offshore Floating Raft Aquaculture Areas in Multispectral Remote Sensing Images
Next Article in Special Issue
Internet of Things-Based Automated Solutions Utilizing Machine Learning for Smart and Real-Time Irrigation Management: A Review
Previous Article in Journal
Designing a Hybrid Energy-Efficient Harvesting System for Head- or Wrist-Worn Healthcare Wearable Devices
Previous Article in Special Issue
YOLOv8 Model for Weed Detection in Wheat Fields Based on a Visual Converter and Multi-Scale Feature Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calibration and Validation of Simulation Parameters for Maize Straw Based on Discrete Element Method and Genetic Algorithm–Backpropagation

1
College of Mechanical and Electronic Engineering, Shandong Agriculture and Engineering University, Jinan 250100, China
2
College of Mechanical and Electrical Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2024, 24(16), 5217; https://doi.org/10.3390/s24165217
Submission received: 26 June 2024 / Revised: 29 July 2024 / Accepted: 9 August 2024 / Published: 12 August 2024
(This article belongs to the Special Issue Sensor and AI Technologies in Intelligent Agriculture: 2nd Edition)

Abstract

:
There is a significant difference between the simulation effect and the actual effect in the design process of maize straw-breaking equipment due to the lack of accurate simulation model parameters in the breaking and processing of maize straw. This article used a combination of physical experiments, virtual simulation, and machine learning to calibrate the simulation parameters of maize straw. A bimodal-distribution discrete element model of maize straw was established based on the intrinsic and contact parameters measured via physical experiments. The significance analysis of the simulation parameters was conducted via the Plackett–Burman experiment. The Poisson ratio, shear modulus, and normal stiffness of the maize straw significantly impacted the peak compression force of the maize straw and steel plate. The steepest-climb test was carried out for the significance parameter, and the relative error between the peak compression force in the simulation test and the peak compression force in the physical test was used as the evaluation index. It was found that the optimal range intervals for the Poisson ratio, shear modulus, and normal stiffness of the maize straw were 0.32–0.36, 1.24 × 108–1.72 × 108 Pa, and 5.9 × 106–6.7 × 106 N/m3, respectively. Using the experimental data of the central composite design as the dataset, a GA–BP neural network prediction model for the peak compression force of maize straw was established, analyzed, and evaluated. The GA–BP prediction model’s accuracy was verified via experiments. It was found that the ideal combination of parameters was a Poisson ratio of 0.357, a shear modulus of 1.511 × 108 Pa, and a normal stiffness of 6.285 × 106 N/m3 for the maize straw. The results provide a basis for analyzing the damage mechanism of maize straw during the grinding process.

1. Introduction

Maize straw is an important renewable resource and feed source for animal husbandry. With the development of the agricultural circular economy, the recycling of agricultural solid waste has become an important research direction. The efficient utilization of crop straw is the key issue in this regard [1,2,3]. The deep processing and utilization of maize straw can alleviate the supply pressure of feed and reduce the residue of agricultural solid waste. It is also significant for promoting the development of biomass energy [4,5,6]. Using the discrete element method to establish a maize straw simulation model and carrying out parameter optimization design of the mechanism can improve straw productivity and utilization rates and promote the rapid development of straw-processing machinery.
Some researchers have carried out multiple studies on the calibration of the simulation parameters of agricultural materials. It is necessary to calibrate the relevant input parameters to establish an accurate discrete element parameter model and obtain a better direct interaction relationship between the material and the contact material [7,8,9]. Ma et al. [10] established a calibration model of mixed parameters for red clover seed and coated powder using the JCR model and verified the accuracy of the model parameters via experiments. Zhang et al. [11] constructed a calibration model for the mixed parameters of a maize root system and soil, providing a reference for the discrete element simulation of maize no-tillage operations. Jung et al. [12] established a simulation model for the flow parameters of soybeans with different pore sizes and water content and predicted the flow performance of soybeans with different water content via experiments. Zhang et al. [13] analyzed the relationship between different structures and bonding parameters of banana straw and established a high-precision discrete element model of banana straw. Zhao et al. [14] established a discrete element model of cotton straw compression and experimentally verified that the established model could accurately simulate the compression process of cotton straw. Tang et al. [15] studied the mechanical properties of rice straws of different lengths under vibration and non-vibration compression conditions, providing basic data for the development of straw compression machinery. Shi et al. [16] established a discrete element model of a wheat straw monomer and conducted physical and simulation experiments on the properties of wheat straw, such as tensile strength, compression, three-point bending, and shear stress. Zheng et al. [17] carried out discrete element parameter calibration for the cutting process of corn stalks. Traditional machine learning models require numerous datasets to be analyzed to obtain satisfactory results. The backpropagation (BP) neural network model can use less data to build the model [18,19]. Many scholars have studied the genetic algorithm (GA)–BP algorithm. Li et al. [20] compared the photovoltaic power generation prediction model constructed using three neural network algorithms, BP, GA–BP, and PSO (particle swarm optimization)–BP (backpropagation), and improved the prediction accuracy of photovoltaic power generation. Wei et al. [21] effectively predicted tool wear based on a GA–BP neural network model. Liu D et al. [22] compared the accuracy of the BP and GA–BP neural network models for soil water prediction. The test results showed that the GA–BP model could be used to predict soil water for ecological protection. However, the simulation parameters for maize straw calibration using a combination of the discrete element method, physical experiments, and machine learning have not been consulted.
This article used a combination of physical experiments, virtual simulation, and machine learning to calibrate the simulation parameters of maize straw. A bimodal-distribution discrete element model of maize straw was established based on the intrinsic and contact parameters via physical experiments. Taking the relative error of the peak compression force from physical tests and simulations as the test index, the Plackett–Burman test, steepest-climb test, and central composite design test were successively carried out. The results of the central composite design experiment were used as the dataset. GA–BP was used for cycle iteration, and the model’s number of cycle iterations was set. After the iteration was completed, the selection was stopped, and the individual with the closest fitness was obtained. The accuracy of the GA–BP prediction model was verified via physical experimentation. The results provide a theoretical basis for analyzing the damage mechanism of maize straw during the grinding process.

2. Materials and Methods

2.1. Physical Test of Radial Compression of Maize Straw

Natural air-dried maize straw harvested in Xujiazhuang Village, Zhoucun District, Zibo City, was used to make the experimental materials. As depicted in Figure 1, the variety was NK 815. During the experiment, 10 stalks were randomly selected, and the leaf sheaths, bracts, and roots of the corn stalks were removed. The third section of maize straw from the ground was selected. As shown in Table 1, the length and diameter of the maize straw were 90 ± 4.92 mm and 23.78 ± 3.21 mm, respectively. The average moisture content of the maize straw was 8.34%. All the maize straw was transferred into a measuring cylinder, and the volume of discharged water indicated the volume of the maize straw. The quality of the maize straw was measured with a balance. The density was 136.4 kg/m3.
A radial compression test of the maize straw was conducted using a universal testing machine, as shown in Figure 2. A corn straw with a length of 90 mm was selected and placed horizontally in the center of the support base. The experimental loading speed was 20 mm/min. The experiment was stopped after the maize straw was visibly crushed. After repeating the process 10 times, the peak compression force of the maize straw was 1798 N.

2.2. Measurement of Contact Mechanics Parameters

2.2.1. Static Friction Coefficient Measurement

The static friction coefficients of the maize straw–maize straw and the maize straw–steel plate were measured with an inclined-plane instrument [23]. Before the experiment, the maize straws were neatly and tightly bonded together to form a straw bottom plate. Each plate was placed on the inclined-plane instrument on a horizontal test bench and stuck to the testing plane of the inclined-plane tester. The angle of the test plane was adjusted until the maize straw was observed to slide. The angle on the digital display inclinometer was recorded, as shown in Figure 3. The static friction coefficient between the maize stalks was calculated using Equation (1). Each group of tests was repeated 10 times, taking the average value as the final value. The static friction coefficients of the maize straw–maize straw and the maize straw–steel plate were 0.28 ± 0.05 and 0.45 ± 0.08, respectively.
μ 1 = tan φ 1
Here, μ 1 is the static friction coefficient, and φ 1 is the critical angle of the static friction coefficient, °.

2.2.2. Rolling Friction Coefficient Measurement

The instruments and testing methods used in the measurement experiment of the rolling friction coefficient were similar to those used for the static friction coefficient [24,25]. The angle of the test plane was adjusted until maize straw rolling was observed. The angle on the digital display inclinometer was recorded. Each group of tests was repeated 10 times, taking the average value as the final value. The rolling friction coefficients of the maize straw–maize straw and the maize straw–steel plate were 0.16 ± 0.15 and 0.24 ± 0.24, respectively.

2.3. Establishment of Discrete Element Simulation Model

2.3.1. Hertz–Mindlin with Bonding Contact Model

The Hertz–Mindlin with bonding contact model [26] was used to bond the particles with a finite size in the maize straw model. A bond could resist tangential and normal displacement until the maximum normal and tangential shear stress was reached, and the bond broke. The Hertz–Mindlin contact model was used to calculate the interaction between particles before bond formation. After the bond was generated at a certain time, the bond force and torque were set to zero, and the superimposed increment in the bond force and torque was found in each time step. The normal and tangential strains were calculated using Equation (2):
F n = v n S n A δ t F t = v t S t A δ t M n = ω n S t J δ t M t = ω t S n J 2 δ t
where A = π R B 2 is the contact area between particles (mm2); F n is the normal force acting on the particles (N); F t is the tangential force acting on the particles (N); M n is the normal moment of the particle (N∙m); M t is the tangential moment of the particles (N∙m); S n is the normal bonding stiffness (N/m3); S t is the tangential bonding stiffness (N/m3); v n is the normal velocity of particle motion (m/s); v t is the tangential velocity of the particles (m/s); ω n is the normal relative velocity of the particles (rad/s); ω t is the tangential relative velocity of the particles (rad/s); δ t is the time step size (s); J = 0.5 π R B 4 is the moment of inertia of the particles (kg∙m2); and R B is the contact radius of the particles (mm).
The bonding force of the maize straw particles mainly depends on five parameters: the normal stiffness coefficient, the tangential stiffness coefficient, the critical normal stress, the critical tangential stress, and the bonding radius. To simplify the simulation calculation, the normal stiffness coefficient of the maize straw particles was equal to the tangential stiffness coefficient, and the critical normal stress was equal to the critical tangential stress. When the normal and tangential stresses on the bond were greater than the set values of the normal or tangential ultimate stress, the bond broke. The calculation formulas for the normal ultimate stress and tangential ultimate stress are shown in Equation (3). The initial range was obtained via the simulation of several maize straw particles’ radial compression, which was the basis for the subsequent simulation test.
σ max < F n A b + 2 M t J R B τ max < F t A + M n J R B

2.3.2. Establishment of Simulation Model for Radial Compression of Maize Straw

In the modeling process, the maize straw was equivalent to an isotropic structure [27,28]. A geometric model of the maize straw particles was created using the Solidworks 2022. The maize straw model was converted into .* stp format and then imported into EDEM 2018. The bimodal-distribution stacking method was used to build the maize straw model, as shown in Figure 4. Spherical particles were selected as the base particles of the maize straw model. The particle radius was set to 1.5 mm. The total number of particles was 3534, among which the large particles occupied the main spatial position, the small particles were closely arranged around the large particles, and the particles had a high coordination number. The adhesive force was stronger, which made the mechanical properties of the particle groups more similar to the actual situation and reduced the computer simulation load [29].
As shown in Figure 5, the Hertz–Mindlin model (no slip) was selected as the contact model between the maize straw and the compression and support base. The simulated loading speed was 20 mm/min for the maize straw compression. The experiment was stopped after the maize straw was visibly crushed. The Rayleigh time step was set to 15%, the data saving interval was 0.01 s, and the grid size was three times the minimum particle radius.

2.4. Calibration of Simulation Parameters for Maize Straw

The Plackett–Burman test was used to analyze the significance of the simulation parameters. The optimal range of the significance parameters was determined using the steepest-climb test. The central composite design was carried out, taking the relative error between the peak compression force of the simulation test and that of the physical test as the experimental indicator. A GA–BP neural network prediction model for the peak compression force of maize straw was established using the experimental data of the central composite design as the dataset. The GA–BP prediction model was analyzed and evaluated. The accuracy of the GA–BP prediction model was verified via physical experiments. Finally, the optimal combination of simulation parameters for the maize straw was obtained.

2.4.1. The Plackett–Burman Test

The parameter range of the Plackett–Burman test was based on the physical experiment results. The other simulation parameters were found in the relevant literature [30,31,32,33,34]. The relative error between the peak compression force of the simulation test and that of the physical test was taken as the test index. The parameters that significantly influenced the response value were selected via the Plackett–Burman test. The minimum and maximum values of the test parameters in Table 2 are coded as −1 and +1, representing the low and high parameter levels, respectively.

2.4.2. Steepest-Climb Test

The Plackett–Burman test was used to analyze the significance of the simulation parameters. Then, the optimal range of the significance parameters was determined by the steepest-climb test. All other non-significant parameters were averaged. The relative error between the peak compression force of the simulation test and that of the physical test is shown in Equation (4):
P = F F f F × 100 %
where P is the relative error between the peak compression force of the simulation test and that of the physical test (%); F is the peak compression force of the simulation test (N); and Ff is the peak compression force of the physical test (N).

2.5. Regression Fitting Modeling Based on Machine Learning Algorithms

2.5.1. Principle of BP Neural Network

The BP (backpropagation) neural network is a multi-layer feedforward neural network trained using the error backpropagation algorithm, mainly composed of three layers: the input layer, hidden layer, and output layer [35,36]. Neurons connect each layer to another, transmitting information between them. The gradient descent method adjusts the weight and threshold of the neural network via error backpropagation to minimize the network error. As shown in Figure 6, the input (xn) of each neuron is multiplied by their respective weights (ωn) and then subtracted from the bias vector (b). The sum of the inputs is passed through an activation function to obtain a specific neuron output (y).

2.5.2. Sample Construction

The optimal range of the significant influencing factors was obtained via the Plackett–Burman test and the steepest-climb test. The central composite design was carried out by taking the relative error between the peak compression force of the simulation test and that of the physical test as the experimental indicator. The horizontal coding of the simulation parameters for maize straw is shown in Table 3.

2.5.3. BP Neural Network Construction and Training

In the construction process of the BP neural network, the number of input layer nodes was determined by the number of input parameters. The number of the output layer nodes was determined by the number of output parameters. Multiple hidden layers are often added between the input and output layers. When the number of hidden layers increases, the accuracy increases, but the network structure is complex, and the learning efficiency decreases. When the input node of the BP artificial neural network is n, the number of hidden layer nodes in the network is selected as 2n + 1 [37]. This article adopted a single hidden layer structure, with 7 hidden layer nodes selected. The BP neural network model adopted a three-layer network, with a neural network structure of 3-7-1. The network topology is shown in Figure 7.
The results of the central composite design experiment were used as the dataset. The experimental design and results are shown in Table 4. The total data (21 groups) were randomly divided into 17 groups (80%) for training, and 6 groups (20%) for testing and verification to avoid over-training and over-parameterization. In the training process, the transfer function from the input layer to the hidden layer was selected as the sigmoid function. The nonlinear damped least-square (LM) optimization algorithm was used in the training algorithm. The mapminmax function was selected to normalize the input and output data to eliminate dimensional effects.

2.5.4. Optimization of the BP Neural Network Model with Genetic Algorithm (GA–BP)

The genetic algorithm was used to optimize the initial weight and threshold of the BP neural network model, improve its computational efficiency and prediction accuracy, and build the collaborative mechanism of the genetic algorithm and BP neural network model. The GA–BP neural network is a hybrid algorithm combining the genetic algorithm and error backpropagation algorithm to train a feedforward artificial neural network, which can accelerate the convergence speed and avoid local minimization. This network converges quickly and easily reaches the optimal solution. The process of the GA–BP optimization algorithm involved using the genetic algorithm to optimize the initial weights and thresholds of the BP network before executing the BP algorithm. After genetic completion, the optimized initial weights and thresholds were assigned to the BP neural network for updating learning so that the optimized BP neural network could better predict the output of the function, as shown in Figure 8.

2.5.5. Data Analysis and Processing

Matlab R2022b software was used as the algorithm-running platform. The predictive performance of the machine learning model was evaluated by the R2, RMSE, and MAE. The larger the R2, the higher the model fit. The lower the RMSE and MAE, the better the model’s accuracy and stability. GA–BP was used for cycle iteration, and the number of cycle iterations of the model was set. After the iteration was completed, the selection was stopped and the individual with the closest fitness was obtained. The accuracy of the model was verified via physical experiments.

3. Results and Analysis

3.1. Analysis of Plackett–Burman Test Results

Since many factors affected the compression test between the maize straw and steel plate, the Plackett–Burman test was needed to determine the significance of each factor’s influence on the radial compression test. The contact parameters between the maize straw and steel plate were screened with the peak compressive force as the response value. A total of 12 groups of tests were carried out. Each group of tests was repeated three times, and the average value was taken. The experimental design and results are shown in Table 5. The significance analysis of the Plackett–Burman test results is shown in Table 6. The significance analysis revealed that X1, X2, and X5 had a significant impact on the peak compression force of the maize straw and steel plate. Therefore, further analysis of the influence law of the X1, X2, and X5 factors on the peak compressive force was needed.
The design scheme and results of the steepest-climb test are shown in Table 7. The peak compression force increased with the increase in the Poisson ratio, shear modulus, and normal stiffness coefficient of the maize straw. The relative error first increased and then decreased. The relative error of the third group was the smallest, and the central point of the central composite design test was selected. The optimal range intervals for X1, X2, and X5 were determined to be 0.32–0.36, 1.24 × 108–1.72 × 108 Pa, and 5.9 × 106–6.7 × 106 N/m3, respectively.

3.2. Regression Model Based on Machine Learning

3.2.1. GA–BP Model Training Results

The change law of the measured and predicted values of the GA–BP is shown in Figure 9. The evaluation index R2 was 0.9069, the RMSE was 0.5524, and the MAE was 0.7763, which showed a good performance in model accuracy, stability, and fitting. This shows that the GA–BP algorithm achieved a better fitting effect in this study and constructed the model with higher precision and less error. This model could be used for further analysis.

3.2.2. Model Evaluation

The MSE performance evaluation of the GA–BP algorithm model was carried out, as shown in Figure 10. The MSE of the model showed a decreasing trend during the training process. The fitting effect of the model to the training data gradually improved as the training progressed. The best performance was obtained at the third step of training, which indicates that the GA–BP model’s training convergence is fast and stable.
The training, validation, testing, and comprehensive performance changes of the GA–BP are shown in Figure 11. The correlation coefficients for the training, validation, testing, and comprehensive data were 0.9930, 0.9994, 0.9989, and 0.9386. This indicates that the model had a strong fitting effect and good generalization ability. The correlation coefficients of the data were very close, indicating no obvious overfitting or underfitting. The GA–BP algorithm was used to obtain a model with high precision and strong generalization ability, which could be used for subsequent experimental research.

3.2.3. GA–BP Optimization Test

The GA–BP algorithm was used for cycle iteration, and the number of cycle iterations of the model was set to 150 times. After the iteration was completed, the selection was stopped, and the individual with the closest fitness was obtained. The results show that the Poisson ratio, shear modulus, and normal stiffness of the maize straw were 0.357, 1.511 × 108 Pa, and 6.285 × 106 N/m3, respectively. The peak compression force simulation experiment of the maize straw was conducted using the GA–BP algorithm optimized parameter combination, and the relative error was 1.14%. The established model can be used for the discrete element simulation of maize straw crushing.

3.3. Discussion

The intrinsic parameters and contact parameters of the maize straw were measured via physical tests. The static friction coefficients of the maize straw–maize straw and the maize straw–steel plate were 0.28 and 0.45, respectively. The rolling friction coefficients of the maize straw–maize straw and the maize straw–steel plate were 0.16 and 0.24, respectively. The other simulation parameters can be found in the relevant literature. The significance analysis of the simulation parameters was conducted using the Plackett–Burman experiment. It was found that the Poisson ratio, shear modulus, and normal stiffness of the maize straw significantly impacted the peak compression force of the maize straw and steel plate. The steepest-climb test was carried out for the significance parameter, and the relative error between the peak compression force in the simulation test and that in the physical test was used as the evaluation index. The optimal range intervals for the Poisson ratio, shear modulus, and normal stiffness of the maize straw were determined to be 0.32–0.36, 1.24 × 108–1.72 × 108 Pa, and 5.9 × 106–6.7 × 106 N/m3, respectively. A GA–BP algorithm neural network prediction model for the peak compression force of the maize straw was established using the experimental data of the central composite design as the dataset and was analyzed and evaluated. The accuracy of the GA–BP algorithm prediction model was verified via experiments. The ideal combination of parameters was found to be 0.357 for the Poisson ratio, 1.511 × 108 Pa for the shear modulus, and 6.285 × 106 N/m3 for the normal stiffness of the maize straw.
In this paper, the Rayleigh time step was set to 15%, which was slightly high for a simulation including breakage. We will change the Rayleigh time step, particle size, and particle ratio to carry out further research.

4. Conclusions

(1) The intrinsic parameters and contact parameters of the maize straw were measured via physical tests. The static friction coefficients of the maize straw–maize straw and the maize straw–steel plate were 0.28 and 0.45, respectively. The rolling friction coefficients of the maize straw–maize straw and the maize straw–steel plate were 0.16 and 0.24, respectively. Other simulation parameters can be found in the relevant literature.
(2) The significance analysis of the simulation parameters was conducted via the Plackett–Burman experiment. It was found that the Poisson ratio, shear modulus, and normal stiffness of the maize straw significantly impacted the peak compression force of the maize straw and steel plate. The steepest-climb test was carried out for the significance parameter, and the relative error between the peak compression force in the simulation test and that in the physical test was used as the evaluation index. The optimal range intervals for the Poisson ratio, shear modulus, and normal stiffness of the maize straw were determined to be 0.32–0.36, 1.24 × 108–1.72 × 108 Pa, and 5.9 × 106–6.7 × 106 N/m3, respectively.
(3) A GA–BP algorithm neural network prediction model for the peak compression force of the maize straw was established using the experimental data of the central composite design as the dataset and was analyzed and evaluated. The accuracy of the GA–BP algorithm prediction model was verified via experiments. The ideal combination of parameters was 0.357 for the Poisson ratio, 1.511 × 108 Pa for the shear modulus, and 6.285 × 106 N/m3 for the normal stiffness of the maize straw.

Author Contributions

Conceptualization, Z.Z.; methodology, H.D. and F.Z.; software, F.Z. and H.D.; data curation, Y.L. and J.C.; formal analysis and writing—original draft preparation, D.J. and H.D.; writing—review and editing, F.Z. and M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by the Program for Innovative Research Team in SDAEU (sgykycxtd2020-03).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

The authors gratefully acknowledge the financial support provided by Shandong Agriculture and Engineering University. We also appreciate the work of the editors and the reviewers of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Duque-Acevedo, M.; Belmonte-Ureña, L.J.; Yakovleva, N.; Camacho-Ferre, F. Analysis of the circular economic production models and their approach in agriculture and agricultural waste biomass management. Int. J. Environ. Res. Public Health 2020, 17, 9549. [Google Scholar] [CrossRef] [PubMed]
  2. Esposito, B.; Sessa, M.R.; Sica, D.; Malandrino, O. Towards circular economy in the agri-food sector. A systematic literature review. Sustainability 2020, 12, 7401. [Google Scholar] [CrossRef]
  3. Stegmann, P.; Londo, M.; Junginger, M. The circular bioeconomy: Its elements and role in European bioeconomy clusters. Resour. Conserv. Recycl. X 2020, 6, 100029. [Google Scholar] [CrossRef]
  4. Korai, P.K.; Sial, T.A.; Pan, G.; Abdelrahman, H.; Sikdar, A.; Kumbhar, F.; Channa, S.A.; Ali, E.F.; Zhang, J.; Rinklebe, J.; et al. Wheat and maize-derived water-washed and unwashed biochar improved the nutrients phytoavailability and the grain and straw yield of rice and wheat: A field trial for sustainable management of paddy soils. J. Environ. Manag. 2021, 297, 113250. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, W.; Liu, Y.; Liu, G.; Xie, R.; Ming, B.; Yang, Y.; Guo, X.; Wang, K.; Xue, J.; Wang, Y. Estimation of maize straw production and appropriate straw return rate in China. Agric. Ecosyst. Environ. 2022, 328, 107865. [Google Scholar] [CrossRef]
  6. Lu, H.; Chen, Y.; Zhang, P.; Huan, H.; Xie, H.; Hu, H. Impacts of farmland size and benefit expectations on the utilization of straw resources: Evidence from crop straw incorporation in China. Soil Use Manag. 2022, 38, 929–939. [Google Scholar] [CrossRef]
  7. Bu, P.; Li, Y.; Zhang, X.; Wen, L.; Qiu, W. A calibration method of discrete element contact model parameters for bulk materials based on experimental design method. Powder Technol. 2023, 425, 118596. [Google Scholar] [CrossRef]
  8. Coetzee, C.J. Calibration of the discrete element method. Powder Technol. 2017, 310, 104–142. [Google Scholar] [CrossRef]
  9. Li, X.; Liu, Y.; An, F.; Wong, H.; Xu, Y. Flow Field Distribution and Morphology Variation of Particles in Planetary Ball Milling. Acta Armamentarii 2022, 43, 876. [Google Scholar]
  10. Ma, X.; Liu, M.; Hou, Z.; Li, J.; Gao, X.; Bai, Y.; Guo, M. Calibration and Experimental Studies on the Mixing Parameters of Red Clover Seeds and Coated Powders. Processes 2022, 10, 2280. [Google Scholar] [CrossRef]
  11. Zhang, S.; Yang, F.; Dong, J.; Chen, X.; Liu, Y.; Mi, G.; Wang, T.; Jia, X.; Huang, Y.; Wang, X. Calibration of discrete element parameters of maize root and its mixture with soil. Processes 2022, 10, 2433. [Google Scholar] [CrossRef]
  12. Jung, H.; Yoon, W.B. Determination and validation of discrete element model parameters of soybeans with various moisture content for the discharge simulation from a cylindrical model silo. Processes 2022, 10, 2622. [Google Scholar] [CrossRef]
  13. Zhang, S.; Jiang, J.; Wang, Y. A Study on the Physical Properties of Banana Straw Based on the Discrete Element Method. Fluid Dyn. Mater. Process. 2022, 19, 1159–1172. [Google Scholar] [CrossRef]
  14. Zhao, W.; Chen, M.; Xie, J.; Cao, S.; Wu, A.; Wang, Z. Discrete element modeling and physical experiment research on the biomechanical properties of cotton stalk. Comput. Electron. Agric. 2023, 204, 107502. [Google Scholar] [CrossRef]
  15. Tang, H.; Xu, W.; Zhao, J.; Xu, C.; Wang, J. Comparison of rice straw compression characteristics in vibration mode based on discrete element method. Biosyst. Eng. 2023, 230, 191–204. [Google Scholar] [CrossRef]
  16. Shi, Y.; Jiang, Y.; Wang, X.; Thuy, N.T.D.; Yu, H. A mechanical model of single wheat straw with failure characteristics based on discrete element method. Biosyst. Eng. 2023, 230, 1–15. [Google Scholar] [CrossRef]
  17. Zheng, Z.; Zhao, H.; Liu, P.; He, J. Maize straw cutting process modelling and parameter calibration based on discrete element method (DEM). INMATEH-Agric. Eng. 2021, 63, 461–468. [Google Scholar] [CrossRef]
  18. Hammoudi, A.; Moussaceb, K.; Belebchouche, C.; Dahmoune, F. Comparison of artificial neural network (ANN) and response surface methodology (RSM) prediction in compressive strength of recycled concrete aggregates. Constr. Build. Mater. 2019, 209, 425–436. [Google Scholar] [CrossRef]
  19. Veza, I.; Spraggon, M.; Fattah, I.R.; Idris, M. Response surface methodology (RSM) for optimizing engine performance and emissions fueled with biofuel: Review of RSM for sustainability energy transition. Results Eng. 2023, 18, 101213. [Google Scholar] [CrossRef]
  20. Li, Y.; Zhou, L.; Gao, P.; Yang, B.; Han, Y.; Lian, C. Short-term power generation forecasting of a photovoltaic plant based on PSO-BP and GA-BP neural networks. Front. Energy Res. 2022, 9, 824691. [Google Scholar] [CrossRef]
  21. Wei, W.; Cong, R.; Li, Y.; Abraham, A.D.; Yang, C.; Chen, Z. Prediction of tool wear based on GA-BP neural network. Proc. Inst. Mech. Eng. Part B J. Eng. Manuf. 2022, 236, 1564–1573. [Google Scholar] [CrossRef]
  22. Liu, D.; Liu, C.; Tang, Y.; Gong, C. A GA-BP neural network regression model for predicting soil moisture in slope ecological protection. Sustainability 2022, 14, 1386. [Google Scholar] [CrossRef]
  23. Fang, M. Study on the Particle Motion Characteristics and Throwing Mechanism of the Disc Knife Chaff Cutter on Gas-Soild Coupling. Ph.D. Thesis, Inner Mongolia Agricultural University, Hohhot, China, 2022. [Google Scholar]
  24. Xiao, Z.Q. Design and Experimental Study on Dust Removal Device of Straw Feed Harvester. Master’s Thesis, Inner Mongolia Agricultural University, Hohhot, China, 2022. [Google Scholar]
  25. Wang, H.Y. Analysis of Heat and Stress Transfer of Compression with Assisted Vibration for Alfalfa Straw Based on Discrete Element Method. Ph.D. Thesis, Inner Mongolia Agricultural University, Hohhot, China, 2021. [Google Scholar]
  26. Yan, D.X. Particle Modelling of Soybean Seeds and the Simulation Analysis and Experimental Study of the Seed-Throwing and Pressing. Ph.D. Thesis, Jilin University, Changchun, China, 2021. [Google Scholar]
  27. Liu, W.; Su, Q.; Fang, M.; Zhang, J.; Zhang, W.; Yu, Z. Parameters calibration of discrete element model for corn straw cutting based on Hertz-Mindlin with bonding. Appl. Sci. 2023, 13, 1156. [Google Scholar] [CrossRef]
  28. Wang, H.; Fan, Z.; Wulantuya; Wang, C.; Ma, Z. Parameter calibration of discrete element model for simulation of crushed corn stalk screw conveying. J. Agric. Sci. Technol. 2023, 25, 96–106. [Google Scholar]
  29. Zhang, F.W.; Song, X.F.; Zhang, X.K.; Zhang, F.Y.; Wei, W.C.; Dai, F. Simulation and experiment on mechanical characteristics of kneading and crushing process of corn straw. Trans. Chin. Soc. Agric. Eng. 2019, 35, 58–65. [Google Scholar]
  30. Hao, J.J.; Long, S.F.; Li, H.; Jia, Y.L.; Ma, Z.K.; Zhao, J.G. Development of discrete element model and calibration of simulation parameters for mechanically-harvested yam. Trans. Chin. Soc. Agric. Eng. 2019, 35, 34–42. [Google Scholar]
  31. Liu, Y.C.; Zhang, F.W.; Song, X.F.; Wang, F.; Zhang, F.Y.; Li, X.Z.; Cao, X.Q. Study on mechanical properties for corn straw of double-layer bonding model based on discrete element method. J. Northeast Agric. Univ. 2022, 53, 45–54. [Google Scholar]
  32. Pang, X.X. Analysis of positioning error of CNC machine tool table based on GA-BP network. China Met. Equip. Manuf. Technol. 2024, 59, 75–77. [Google Scholar]
  33. Wang, X.; Tian, H.; Xiao, Z.; Zhao, K.; Li, D.; Wang, D. Numerical Simulation and Experimental Study of Corn Straw Grinding Process Based on Computational Fluid Dynamics–Discrete Element Method. Agriculture 2024, 14, 325. [Google Scholar] [CrossRef]
  34. Zhu, H.B.; Qian, C.; Bai, L.Z.; Li, H.; Mou, D.L.; Li, J.J. Optimaztion of discrete element model of corn stalk based on Plackett-Burman design and response surface methodology. J. China Agric. Univ. 2021, 26, 221–231. [Google Scholar]
  35. Chen, G.; Tang, B.; Zeng, X.; Zhou, P.; Kang, P.; Long, H. Short-term wind speed forecasting based on long short-term memory and improved BP neural network. Int. J. Electr. Power Energy Syst. 2022, 134, 107365. [Google Scholar] [CrossRef]
  36. Li, Y.; Li, J.; Huang, J.; Zhou, H. Fitting analysis and research of measured data of SAW micro-pressure sensor based on BP neural network. Measurement 2020, 155, 107533. [Google Scholar] [CrossRef] [PubMed]
  37. Meng, Z.P.; Tian, Y.D.; Lei, Y. Prediction Models of Coal Bed Gas Content Based on BP Neural Networks and Its Applications. J. China Univ. Min. Technol. 2008, 37, 456–461. [Google Scholar]
Figure 1. The maize straw.
Figure 1. The maize straw.
Sensors 24 05217 g001
Figure 2. Physical experiment on radial compression of the maize straw.
Figure 2. Physical experiment on radial compression of the maize straw.
Sensors 24 05217 g002
Figure 3. Static friction coefficient measurement device. 1. Digital display inclinometer; 2. tested maize straw; 3. maize straw base plate; 4. test plane.
Figure 3. Static friction coefficient measurement device. 1. Digital display inclinometer; 2. tested maize straw; 3. maize straw base plate; 4. test plane.
Sensors 24 05217 g003
Figure 4. Establishment of the simulation model for maize straw.
Figure 4. Establishment of the simulation model for maize straw.
Sensors 24 05217 g004
Figure 5. Establishment of a simulation model for maize straw compression.
Figure 5. Establishment of a simulation model for maize straw compression.
Sensors 24 05217 g005
Figure 6. Neuron model.
Figure 6. Neuron model.
Sensors 24 05217 g006
Figure 7. Neural network topology.
Figure 7. Neural network topology.
Sensors 24 05217 g007
Figure 8. The training process of GA–BP optimization algorithm.
Figure 8. The training process of GA–BP optimization algorithm.
Sensors 24 05217 g008
Figure 9. The variation pattern between measured and predicted values.
Figure 9. The variation pattern between measured and predicted values.
Sensors 24 05217 g009
Figure 10. The MSE performance evaluation of the model.
Figure 10. The MSE performance evaluation of the model.
Sensors 24 05217 g010
Figure 11. Regression analysis.
Figure 11. Regression analysis.
Sensors 24 05217 g011
Table 1. Size parameters of the maize straw.
Table 1. Size parameters of the maize straw.
Maize Straw12345678910Average Value (mm)
Length of maize straw (mm)8487869186879899938990
Length of maize straw (mm)23.6420.8219.4820.0422.2226.4728.1629.2925.0222.6623.78
Table 2. Plackett–Burman test parameter table.
Table 2. Plackett–Burman test parameter table.
No.Test ParameterCode
−10+1
X1Poisson’s ratio of maize straw0.30.350.4
X2Shear modulus of maize straw (Pa)1 × 1081.6 × 1082.2 × 108
X3Collision recovery coefficient between maize straw0.340.470.60
X4Collision recovery coefficient between maize straw and steel plate0.380.520.66
X5Normal stiffness coefficient (N/m3)5.5 × 1066.5 × 1067.5 × 106
X6Tangential stiffness coefficient (N/m3)5.5 × 1066.5 × 1067.5 × 106
X7Normal critical stress (MPa)4.55.05.5
X8Tangential critical stress (MPa)5.26.06.8
X9Bonding radius (mm)1.82.02.2
Table 3. The horizontal coding table of the simulation parameters.
Table 3. The horizontal coding table of the simulation parameters.
LevelParameter
X1X2X5
−1.6820.2931.413 × 1086.233 × 106
−10.321.44 × 1086.26 × 106
00.341.48 × 1086.30 × 106
+10.361.52 × 1086.34 × 106
+1.6820.3871.547 × 1086.367 × 106
Table 4. Experimental design and results.
Table 4. Experimental design and results.
No.X1X2X5Relative Error P (%)
10001.82
211−11.04
300−1.6824.85
41.682004.4
50002.47
6−1−1−18.76
70002.52
80001.96
9−1−118.93
10001.6827.15
11−1116.41
12−1.682008.24
131117.68
141−1−12.39
150002.88
16−11−15.48
170003.43
180002.15
190002.67
200−1.68207.02
210003.11
221−117.26
2301.68203.54
Table 5. The Plackett–Burman experimental design and results.
Table 5. The Plackett–Burman experimental design and results.
No.ParameterRelative Error P (%)
X1X2X3X4X5X6X7X8X9
11−1−1−11−111−19.14
21−1111−1−1−1122.24
3−1111−1−1−11−14.49
4−1−1−1−1−1−1−1−1−110.29
511−1−1−11−11128.41
6111−1−1−11−118.64
71−111−1111−125.42
8−111−1111−1−115.5
9−11−111−11111.33
10−1−1−11−111−113.52
11−1−11−111−1116.66
1211−1111−1−1−125.29
Table 6. Significance analysis of the Plackett–Burman test.
Table 6. Significance analysis of the Plackett–Burman test.
ParameterEffectMean Square SumImpact RateSignificance Order
X1201.083121,30412.282
X2453.217616,21662.371
X342.91675525.520.569
X499.616729,770.43.016
X5−190.383108,73711.013
X663.7512,192.21.237
X7121.3544,177.54.474
X845.28336151.740.628
X9101.51730,916.93.135
Table 7. Analysis of steepest-climb test results.
Table 7. Analysis of steepest-climb test results.
No.ParameterPeak Compression Force F (N)Relative Error P (%)
X1X2X5
10.31.0 × 1085.5 × 1061437.7125.06%
20.321.24 × 1085.9 × 1061653.798.72%
30.341.48 × 1086.3 × 1061752.782.58%
40.361.72 × 1086.7 × 1061744.613.06%
50.381.96 × 1087.1 × 1062076.4513.41%
60.42.2 × 1087.5 × 1062588.9130.55%
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

Zeng, F.; Diao, H.; Liu, Y.; Ji, D.; Dou, M.; Cui, J.; Zhao, Z. Calibration and Validation of Simulation Parameters for Maize Straw Based on Discrete Element Method and Genetic Algorithm–Backpropagation. Sensors 2024, 24, 5217. https://doi.org/10.3390/s24165217

AMA Style

Zeng F, Diao H, Liu Y, Ji D, Dou M, Cui J, Zhao Z. Calibration and Validation of Simulation Parameters for Maize Straw Based on Discrete Element Method and Genetic Algorithm–Backpropagation. Sensors. 2024; 24(16):5217. https://doi.org/10.3390/s24165217

Chicago/Turabian Style

Zeng, Fandi, Hongwei Diao, Yinzeng Liu, Dong Ji, Meiling Dou, Ji Cui, and Zhihuan Zhao. 2024. "Calibration and Validation of Simulation Parameters for Maize Straw Based on Discrete Element Method and Genetic Algorithm–Backpropagation" Sensors 24, no. 16: 5217. https://doi.org/10.3390/s24165217

APA Style

Zeng, F., Diao, H., Liu, Y., Ji, D., Dou, M., Cui, J., & Zhao, Z. (2024). Calibration and Validation of Simulation Parameters for Maize Straw Based on Discrete Element Method and Genetic Algorithm–Backpropagation. Sensors, 24(16), 5217. https://doi.org/10.3390/s24165217

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