Discrete Meta-Modeling Method of Breakable Corn Kernels with Multi-Particle Sub-Area Combinations

Simulation is an important technical tool in corn threshing operations, and the establishment of the corn kernel model is the core part of the simulation process. The existing modeling method is to treat the whole kernel as a rigid body, which cannot be crushed during the simulation process, and the calculation of the crushing rate needs to be considered through multiple criteria such as the contact force, the number of collisions, and so on. Aiming at the issue that kernel crushing during maize threshing cannot be accurately modeled in discrete element simulations, in this study, a sub-area crushing model was constructed; representative samples with 26%, 30% and 34% moisture content were selected from a double-season maturing region in China; based on the physical dimensions and biological structure of the maize kernel, three stress regions were defined; and mechanical property tests were conducted on each of the three stress regions using a texturometer as a way to determine the different crushing forces due to the heterogeneity of the maize structure. The correctness of the model was verified by stacking angle and mechanical property experiments. A discrete element model of corn kernels was established using the Bonding V2 method and sub-area modeling. Bonding parameters were calculated by combining stacking angle tests and mechanical property tests. The flattened corn kernel was used as a prototype, and the bonding parameters were determined through size and mechanical property tests. A 22-ball bonding model was developed using dimensional parameters, and the kernel density was recalculated. Results showed that the relative error between the stacking angle test and the measured mean value was 0.31%. The maximum deviation of axial compression simulation results from the measured mean value was 22.8 N, and the minimum deviation was 3.67 N. The errors between simulated and actual rupture forces at the three force areas were 5%, 10%, and 0.6%, respectively. The decreasing trend of the maximum rupture force for the three moisture levels in the simulation matched that of the actual rupture force. The discrete element model can accurately reflect the rupture force, energy relationship, and rupture process on both sides, top, and bottom of the grain, and it can solve the error problem caused by the contact between the threshing element and the grain line in the actual threshing process to achieve the design optimization of the threshing drum. The modeling method provided in this study can also be applied to breakable discrete element models for wheat and soybean, and it provides a reference for optimizing the design of subsequent threshing devices.