Discrete element simulation is an effective method to reveal the interaction between tillage components and work objects. However, due to the lack of discrete element modelling parameters of maize root and its mixture with soil, existing tillage models cannot accurately simulate the farmland environment under a no-tillage system. This study developed single maize root (SMR) with different diameters and maize root-soil mixture (MRSM) DEM models based on calibrated parameters through the angle of repose (AOR) tests. First, the Plackett–Burman and the steepest climb tests were performed to identify the range of essential parameters for the AOR of the SMR. Then, the optimal parameters for the SMR and MRSR models were obtained by Box–Behnken design (BBD) testing. The results showed that the static friction coefficient of SMR-SMR and the rolling friction coefficient of SMR-SMR and SMR-steel significantly affected the AOR. In addition, the AOR of MRSM was extremely sensitive to the restitution coefficient and surface energy coefficient of root soil. Based on optimal parameters, the relative errors between the simulated and measured AOR and pixel peak values of the piles’ contour curve were less than 5% for SMR and MRSM. The error of the dynamic AOR of the measured and simulated MSRM was less than 10%. These results indicate that the parameter calibration method and the developed models can be valuable references for DEM simulation for maize stubble and tillage.

Mechanical subsoiling is an effective practice to promote better water infiltration and crop root development. The bending subsoiling tool (BST) is a primary subsoiling tool and is used to remove soil compaction and restore soil productivity. In this study, a discrete element model was developed and validated using laboratory soil bin tests to investigate the effects of the mounting angle of the BST (5°–33°) on soil disturbance behaviors and draft forces. The results show that the upheaval, failure and fragmentation of soil was achieved by successive shearing, uplifting, extrusion, tension and turning actions from the cutting share and cambered shank of the BST. Increasing soil depths gave smaller soil disturbance ranges in lateral, forward and upward directions. With an increase in mounting angle, both the draft force and soil rupture distance ratio initially decreased and then increased, whereas the soil loosening efficiency initially increased and then decreased. Overall, increasing the mounting angle of the BST from 5° to 33° gave a greater soil surface flatness that increased rapidly when the mounting angle increased from 26° to 33°. Appropriately increasing mounting angle of the BST from 5° to 26° could lift more moist soil from the deep seed and middle layers (5.0–15.5% increase) into the shallow seed layer (depth of <50 mm) without seriously affecting the mixing of the deep layer and other layers. Considering both the soil disturbance characteristics and draft forces, a mounting angle of 26° appeared to outperform the other angles.

Discrete element modelling (DEM) parameters are of great importance for the accurate prediction of soil properties and disturbance. This study aimed to provide an efficient method for accurately determining the DEM parameters of a paddy soil with a high moisture content (45.66%). The DEM parameters of the paddy soil modelled using the Hertz–Mindlin with JKR contact model were determined by using the Plackett–Burman (PB), steepest ascent, and central composite tests. The accuracies of the developed DEM models were evaluated using actual slumping tests. Based on the PB test, the surface energy of soil, coefficients of soil–soil rolling friction, and coefficients of soil–steel static friction exerted larger influences on the overall relative error (δ_ZH). The optimization results showed that to achieve a minimum δ_ZH (5.96%), the surface energy of soil, the coefficients of soil–soil rolling friction, and the coefficients of soil–steel static friction should be 0.808 J m⁻², 0.11, and 0.6, respectively. The optimized DEM model had an overall relative error of 7.27% with a coefficient of variation of 1.32%, indicating that the DEM parameters of the calibrated paddy soil had good accuracy.

Oil peony is an important oil crop that is primarily sown by using artificial methods at present. Its seeder encounters the problems of low efficiency of seeding that significantly limits the highly efficient mechanized production of high-quality peony oil. In this study, Fengdan white oil peony seeds were used as the research object and repose angle as the response value to establish a discrete element model (DEM) and parameter calibration. The range of parameters of oil peony seeds was first obtained through an experiment, and their repose angle was obtained by an inclinometer. A three-dimensional DEM of oil peony seeds was then established. The Plackett–Burman (PB) test was utilized to screen the parameters that had a significant influence on the repose angle, and the steepest ascent (SA) test was applied to determine their optimum range of testing. Following this, based on Box–Behnken (BBD) test results, a second-order regression model between the important parameters and the repose angle was constructed. Finally, the absolute minimum difference between simulated and measured repose angles was utilized as the objective of optimization to obtain the following optimum combination of parameters: The values of the seed–steel collision recovery coefficient (CRC), seed–seed static friction coefficient (SFC), seed–steel SFC, and seed–seed rolling friction coefficient (RFC) were 0.704, 0.324, 0.335, and 0.045, respectively. This optimal combination of parameters was confirmed through simulations, and the error between simulated and measured repose angles was only 0.67%, indicating that the calibrated DEM of oil peony seeds was reliable.

Subsoiling practice is an essential tillage practice in modern agriculture. Tillage forces and energy consumption during subsoiling are extremely high, which reduces the economic benefits of subsoiling technology. In this paper, a cicada-inspired biomimetic subsoiling tool (CIST) was designed to reduce the draught force during subsoiling. A soil–tool interaction model was developed using EDEM and validated using lab soil bin tests with sandy loam soil. The validated model was used to optimize the CIST and evaluate its performance by comparing it with a conventional chisel subsoiling tool (CCST) at various working depths (250–350 mm) and speeds (0.5–2.5 ms⁻¹). Results showed that both simulated draught force and soil disturbance behaviors agreed well with those from lab soil bin tests, as indicated by relative errors of <6.1%. Compared with the CCST, the draught forces of the CIST can be reduced by 17.7% at various working depths and speeds; the design of the CIST obviously outperforms some previous biomimetic designs with largest draught force reduction of 7.29–12.8%. Soil surface flatness after subsoiling using the CIST was smoother at various depths than using the CCST. Soil loosening efficiencies of the CIST can be raised by 17.37% at various working speeds. Results from this study implied that the developed cicada-inspired subsoiling tool outperforms the conventional chisel subsoiling tool on aspects of soil disturbance behaviors, draught forces, and soil loosening efficiencies. This study can have implications for designing high-performance subsoiling tools with reduced draught forces and energy requirements, especially for the subsoiling tools working under sandy loam soil.

Discrete element parameters of the black bean (BLB) are key to developing high-performance BLB machineries (e.g., seeders and shellers), which are still lacking in previous literature. In this study, the effects of the radius and lifting speed of cylinder-in-cylinder lifting method (CLM) simulations were investigated to efficiently and accurately obtain the repose angle. Discrete element method (DEM) parameters of the BLB were determined by combining the Plackett–Burman Design test, the steepest ascent design test, and the central composite design test. The results show that the measurement moment (i.e., 12 s) of repose angles should be determined when kinetic energy reaches the minimal threshold (1 × 10⁻⁶ J) to efficiently and accurately obtain repose angles; too early or too late a measurement can result in inaccurate repose angles or excessive computation time of the computer, respectively. The lifting speed and cylinder radius affected the lateral displacements of BLBs and came at the cost of higher computation time and memory usage. A lifting speed of 0.015 m·s⁻¹ and a radius of 40 mm of the cylinder were determined in CLM simulations. The static friction coefficient and rolling friction coefficient between BLBs significantly affected the repose angles. A static friction coefficient of 0.202 and rolling friction coefficient of 0.0104 between BLBs were obtained based on the optimization results. A low relative error (0.74%) and insignificant difference (p > 0.05) between the simulated and measured repose angles were found. The suggested method can be potentially used to calibrate the DEM parameters of BLBs with good accuracy. The results from this study can provide implications for investigating interactions of BLBs and various BLB processing machines and for the efficient and accurate determination of DEM parameters of crop grains.