A modelling and verification approach for soybean seed particles using the discrete element method

Nonlinear effect of grain elongation on the flow rate in silo discharge

By means of two-dimensional numerical simulations based on contact dynamics, we present a systematic analysis of the joint effects of grain shape (i.e., grain elongation) and system size on silo discharge for increasing orifice sizes D. Grains are rounded-cap rectangles whose aspect ratio are varied from 1 (disks) to 7. In order to clearly isolate the effect of grain shape, the mass of the grains is keeping constant as well as the condition of the discharge by reintroducing the exiting grains at the top of the silo. In order to quantify the possible size effects, the thickness W of the silos is varied from 7 to 70 grains diameter, while keeping the silos aspect ratio always equal to 2. We find that, as long as size effects are negligible, the flow rate Q increases as a Beverloo-like function with D, also for the most elongated grains. In contrast, the effects of grain elongation on the flow rate depend on orifice size. For small normalized orifice sizes, the flow rate is nearly independent with grain elongation. For intermediate normalized orifice sizes the flow rate first increases with grain elongation up to a maximum value that depends on the normalized size of the orifice and saturates as the grains become more elongated. For larger normalized orifice size, the flow rate is an increasing function of grains' aspect ratio. Velocity profiles and packing fraction profiles close to the orifice turn out to be self-similar for all grain shapes and for the whole range of orifice and system sizes studied. Following the methodology introduced by Janda et al. [Phys. Rev. Lett. 108, 248001 (2012)PRLTAO0031-900710.1103/PhysRevLett.108.248001], we explain the nonlinear variation of Q with grain elongation, and for all orifice sizes, from compensation mechanisms between the velocity and packing fraction measured at the center of the orifice. Finally, an equation to predict the evolution of Q as a function of the aspect ratio of the grains is deduced.

Study on Verification Approach and Multicontact Points Issue When Modeling Cyperus esculentus Seeds Based on DEM

In this paper, the Multisphere (MS) models of three varieties of Cyperus esculentus seeds are modeled based on DEM. In addition, for comparison, other particle models based on automatic filing in EDEM software are also introduced. Then, the direct shear test, piling test, bulk density test, and rotating hub test are used to verify the feasibility of particle models of Cyperus esculentus seeds that we proposed. By comparing the simulated results and experimental results, combined with the CPU computation time, the proposed particle models achieved better simulation accuracy with fewer filing spheres. According to simulation results, some limitation was present when using one single verification test; varieties of verification tests used could improve the verification reliability, and a more appropriate particle model could be selected. Additionally, the issue of multicontact points in the MS model was studied. The Hertz Mindlin (no slip) (HM) model and Hertz Mindlin new restitution (HMNR) model were both considered in simulations for comparison. The rotating hub test and particle–wall impact test were used, and the influences of multiple contact points on the motion behavior of individual particles and particle assemblies were analyzed. Simulation results showed that the multiple contact points affected the motion behavior of individual particles; in contrast, the influence of multiple contact points on the motion behavior of the particle assembly was insignificant. Moreover, the relationships between moisture content of seeds and Young’s modulus, Young’s modulus, and the number of contact points were also considered. Young’s modulus decreased with increasing moisture content. The number of contact points increased with a decreasing Young’s modulus.

Determination and Validation of Discrete Element Model Parameters of Soybeans with Various Moisture Content for the Discharge Simulation from a Cylindrical Model Silo

This study investigates the physical parameters that affect the flow patterns of soybeans with various moisture content (12% to 60%) at varying orifice sizes (20, 40, and 60 mm) in a cylindrical silo. The flow conditions required to obtain a steady mass flow during discharge were evaluated via experiments and three-dimensional discrete element method (DEM) simulation. The discharged mass flow rates at different flow conditions provided the critical size of the orifice. If the reduced diameter (Dred) of an orifice is >5.59, the flow showed a steady state. Based on the mass flow index (MFI), the flow patterns at 40% and 60% moisture content at 40 and 60 mm orifice sizes, respectively, showed funnel flows. although these flow conditions were satisfied to maintain a steady flow. The maximum wall pressure for the funnel flow showed the location of the interlocking phenomenon where the stagnant zone began during discharging. DEM simulation was validated through the mass profiles using the parameters obtained by the experiments. This study demonstrates that the experimental and analytical results with DEM simulation predict the flow behaviors of soybeans well at various moisture contents. These results are useful for designing silos for continuous food processing.

Test and Simulation Analysis of Soybean Seed Throwing Process

In order to analyze the effect of different factors on the bouncing and rolling distance of soybeans at the time of seed throwing, tests and discrete element method (DEM) are employed to analyze test soil and three representative soybean varieties. The parameters between soybean seed particles and soil particles are calibrated by means of a piling test and simulation. A seed throwing test apparatus is improved to analyze the effects of seed throwing height, soil plane inclination angle and collision orientation on the bouncing and rolling distance of soybean seeds. The effect of relative seed throwing speed on the bouncing and rolling distance of soybean seeds is analyzed using a computer vision seeding test bench. On this basis, the above-mentioned test procedure is simulated and compared with the test results. The results showed that the bouncing distance of the soybean seed particles was not significant. The rolling distance had a certain randomness when the seed throwing height was different. When the inclination of the soil plane became larger, the rolling distance increased and when the sphericity of the soybean seed particles was high, the effect of different collision orientations was not obvious if sphericity was low. The rolling distance was the shortest when colliding in the horizontal orientation and the longest when colliding in the vertical orientation. The larger the relative seed throwing speed, the larger the rolling distance of the soybean seed particles.

Study on the Discharge Process and Mechanism of Anti-Corrosion Pill Particles in the Oil and Gas Field Wellbore Casing Annulus Based on the Discrete Element Method

This research studies the discharge process and its mechanism using the discrete element method (DEM) with self-developed annular corrosion pill particles and the discharge device as an example in order to optimize the oil and gas field wellbore casing annular corrosion process. The object of study was chosen from four different grid numbers and four different grid widths, and EDEM software was utilized to simulate and assess the pill particle discharge process based on preliminary experimental research. Under five different pill wheel rotation speeds, the effects of the grid number and grid width on the filling amount, filling density, discharge variation coefficient, and compressive force of pill particles were investigated from macroscopic and microscopic viewpoints. The findings reveal that the grid number, grid width, and rotation speed all have a significant impact on pill filling and discharge performance. As a result, the discharge wheel’s structure and operating characteristics were optimized. The discharge wheel performs best when the grid number is 8, the grid width is 75 mm, and the rotation speed is 15 rpm; the pill filling density is 692.26 kg/m3, the discharge variation coefficient is 0.022, and the maximum compressive force is 188 N. This study establishes the groundwork for enhancing wellbore integrity management in oil and gas fields by providing a reference for the optimal design of wellbore casing annular corrosion prevention devices in oil and gas fields.

A DEM-Based Modeling Method and Simulation Parameter Selection for Cyperus esculentus Seeds

To build a DEM model of Cyperus esculentus seed particles, the shape and size of the Cyperus esculentus seed particles were measured and analyzed. The results showed that the dispersity in size had a normal distribution. Additionally, a certain functional relationship between the primary dimension and secondary dimensions was determined. The width of the seed was the primary dimension, and the other secondary dimensions (length and thickness) were calculated based on their relationships with the primary dimension. On this basis, an approach for modeling Cyperus esculentus seed particles based on the multi-sphere (MS) method was proposed. The discrete element analysis models of three varieties of Cyperus esculentus seeds were established with different numbers of filing spheres. Moreover, to obtain more accurate simulation parameters, first, a range of values of the simulation parameters was obtained by the experimental method. Second, the Plackett–Burman (PB) test and the path of steepest ascent method were both adopted to correct and calibrate the simulation parameters, which were difficult to obtain through experiments, and simulation of the direct shear test was used for calibration. All of the methods guaranteed that the selected parameters were reasonable. The test results showed that the static friction coefficient of seed–seed had a significant effect on the simulation results. Finally, piling tests and the bulk density test were used for modeling verification. By comparing the simulated results and experimental results in the piling tests and bulk density test, when the number of filing spheres increased, the simulated results were close to those obtained experimentally. Therefore, the feasibility and validity of the modeling method for Cyperus esculentus seed particles that we proposed and the simulation parameters that were obtained were verified.

A Review of the Application of Discrete Element Method in Agricultural Engineering: A Case Study of Soybean

The discrete element method has become a common method for analyzing the contact interaction between particulate materials and between particles and mechanical components. It has been widely used in agricultural engineering and other fields. Taking soybean as an example, soybean seed particles always have contact effects between particles and mechanical components in the process of planting, harvesting, threshing, separation, cleaning, and processing. The discrete element method can be used to obtain information on the contact forces between seed particles and mechanical parts, as well as the velocity and displacement of seed particle motion from a microscopic perspective. This paper summarizes the application of the discrete element method in soybean cultivation and production processes in recent years. This will help future researchers to conduct relevant test studies, develop and improve existing research methods. It can also serve as a guide and reference for the production and processing of other granular materials and the optimization of agricultural machinery components.

Simulation and Analysis of the Working Process of Soil Covering and Compacting of Precision Seeding Units Based on the Coupling Model of DEM with MBD

In precision seeding, the final displacements of the seeds are determined as a working result of the profiling mechanism, opener, seed-metering device, covering apparatus and compacting machine. For a better understanding of the disturbance of seed displacement during soil covering and compaction in the actual working process, experiments and simulations have been performed. In this paper, a type of soybean seeding monomer was taken as the research object, and a soil bin test of soil covering and compacting was executed. The experimental results showed that the traction velocity and the open angle of the covering discs had a significant influence on the changes in the horizontal and vertical displacements of seeds during the soil covering processing. With an increasing traction velocity, the vertical displacements of seeds increased after soil covering; in contrast, the horizontal displacements decreased. When the covering apparatus had a larger open angle it had a smaller disturbance influence on the soil. Therefore, with an increase in the opening angle, the changes in the vertical and horizontal displacements of seeds showed a decreasing trend. Inversely, in the process of compacting, the forward velocity had little effect on the three-dimensional displacement change in the seeds after compacting. The analysis model of the precision seeding unit was established based on the coupling model of the DEM (discrete element method) with MBDs (multi-body dynamics). The process of soil covering and compacting was simulated and analyzed. The comparison between the experimental results and the simulated results showed that the trend was similar, and the two results were close. Thus, the feasibility and applicability of the coupling method were verified. It also provided a new method for the design and optimization of covering and compacting components of a precise seeding monomer.

High-Throughput Legume Seed Phenotyping Using a Handheld 3D Laser Scanner

High-throughput phenotyping involves many samples and diverse trait types. For the goal of automatic measurement and batch data processing, a novel method for high-throughput legume seed phenotyping is proposed. A pipeline of automatic data acquisition and processing, including point cloud acquisition, single-seed extraction, pose normalization, three-dimensional (3D) reconstruction, and trait estimation, is proposed. First, a handheld laser scanner is used to obtain the legume seed point clouds in batches. Second, a combined segmentation method using the RANSAC method, the Euclidean segmentation method, and the dimensionality of the features is proposed to conduct single-seed extraction. Third, a coordinate rotation method based on PCA and the table normal is proposed to conduct pose normalization. Fourth, a fast symmetry-based 3D reconstruction method is built to reconstruct a 3D model of the single seed, and the Poisson surface reconstruction method is used for surface reconstruction. Finally, 34 traits, including 11 morphological traits, 11 scale factors, and 12 shape factors, are automatically calculated. A total of 2500 samples of five kinds of legume seeds are measured. Experimental results show that the average accuracies of scanning and segmentation are 99.52% and 100%, respectively. The overall average reconstruction error is 0.014 mm. The average morphological trait measurement accuracy is submillimeter, and the average relative percentage error is within 3%. The proposed method provides a feasible method of batch data acquisition and processing, which will facilitate the automation in high-throughput legume seed phenotyping.

Measurement and Calibration of the Parameters for Discrete Element Method Modeling of Rapeseed

The discrete element method (DEM) for modeling the behavior of particulate material is highly dependent on the use of appropriate and accurate parameters. In this study, a seed metering DEM simulation was used to measure, calibrate, and verify the physical and interactional parameters of rapeseed. The coefficients of restitution and static friction between rapeseeds and three common materials (aluminum alloy, acrylic, and high-density polyethylene) were measured using free drop and sliding ramp tests, respectively. The angle of repose was determined using a hollow cylinder experiment, which was duplicated using a DEM simulation, to examine the effects of static and rolling friction coefficients on the angle of repose. Response surface optimization was performed to determine the optimized model parameters using a Box–Behnken design test. A metering device was made with three materials, and rapeseed seeding was simulated at different working speeds to verify the calibrated parameters. The validation results showed that the relative errors between the seed metering model and experiments for the single qualified seeding, missed seeding, and multiple seeding rates were −0.15%, 3.29%, and 5.37%, respectively. The results suggest that the determined physical and interactional parameters of rapeseed can be used as references for future DEM simulations.

A Comparative Study on the Modelling of Soybean Particles Based on the Discrete Element Method

To solve the poor universality in the existing modelling approaches of soybean particles, we proposed a soybean particle modelling approach by combining five, nine, and 13 balls. The soybean seeds from three varieties (Suinong42, Jidou17, and Zhongdou39 with a sphericity of 94.78%, 86.86%, and 80.6%, respectively) are chosen as the study objects. By the comparisons between the simulation results and the test results in the “self-flow screening” and “piling angle” tests, it is concluded that the soybean particle modelling approach we presented in this paper is a universal modelling approach appropriate for soybean particles with different sphericities. The five-ball model is appropriate for the soybean particles with high sphericity, and the nine- or 13-ball models are applicable to those with low sphericity. The soybean particle modelling approach we presented is also compared with the ellipsoidal equation modelling approach for soybean particles and with the modelling approaches presented by other researchers. From an overall perspective, the soybean particle modelling approach we presented is better than the ellipsoidal equation modelling approach and those modelling approaches presented by other researchers. Additionally, it is shown that the multiple contacts issue in the multi-ball model has a little influence on the simulation results of soybean particle models. The study in this paper provides a new modelling approach for soybean particles in the DEM simulation of the contacts between soybean particles and the related machines.