Next Article in Journal
Development Results of a Cross-Platform Positioning System for a Robotics Feed System at a Dairy Cattle Complex
Next Article in Special Issue
Research on the Resistance Reduction Law of Self-Excited a Resonant Circular Arc-Surface Bulldozing Plate Based on the Discrete Element Method
Previous Article in Journal
Plant Height and Stem Diameter of Solanum quitoense Lamarck Improved with Applications of AMF and Biostimulants
Previous Article in Special Issue
Self-Excited-Resonance of Soil-Engaging Surface Spectrum: A New Method of Soil Cutting Resistance Reduction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Mechanism of Motion Interaction between Soil and a Bionic Hole-Forming Device

1
College of Mechanical and Electronic Engineering, Tarim University, Alar 843300, China
2
College of Engineering, China Agricultural University, Beijing 100083, China
3
School of Mechanical Engineering, Hunan University of Technology, Zhuzhou 412007, China
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(7), 1421; https://doi.org/10.3390/agriculture13071421
Submission received: 23 June 2023 / Revised: 13 July 2023 / Accepted: 17 July 2023 / Published: 18 July 2023

Abstract

:
Due to the lack of water resources, the main agricultural planting method used in the northwest region of China is plastic film mulching, with precision hole sowing performed on the film after mulching. However, conventional hole-forming devices damage the compactness of the soil hole while moving on the plastic film, causing seed misplacement. Therefore, this study designed a bionic hole-forming device based on the oriental mole cricket. In order to explore the interaction between the hole-forming device and the soil, a typical soil discrete element particle model was established, and its contact parameters were calibrated. An experiment was conducted to compare the performance of the bionic hole-forming device with a conventional device using discrete element method and multi-body dynamics (DEM-MBD) coupled simulations. The results revealed that the bionic hole-forming device caused less soil disturbance during the hole-forming process and could reduce the sowing operation resistance compared to the traditional device. Compared to traditional square and cone-shaped hole-forming devices, the soil resistance of the bionic hole-forming device was the smallest, at 7.51 N. This work provides a reference for the optimization of hole-forming devices for plastic film sowing.

1. Introduction

Due to its unique geographical location and natural environment, the northwest region of China has emerged as a major cotton planting base for the country [1]. The primary method of seeding cotton is through hole seeding on plastic film. However, existing hole-forming devices used in this process create large holes in the plastic film and seriously disturb the soil [2]. Therefore, there is an urgent need to develop hole-forming machinery suitable for the soil conditions of the farmland in Northwest China to reduce soil disturbance and minimize operational resistance. Such equipment would greatly contribute to the sustainable development of the cotton industry.
In order to successfully design soil-engaging components for agricultural machinery, the interaction mechanism between the soil-engaging components and the soil must be clarified [3]. The use of soil trenches and field tests can only provide a macroscopic analysis of soil disturbance and is unable to investigate the soil’s movement patterns and force characteristics at the microscopic scale [4]. Numerous studies have demonstrated the applicability of the discrete element method to simulate and analyze the interaction between soil-engaging components and soil and to optimize the structural parameters of the machinery [5,6,7,8]. However, the type and properties of the soil directly affect the interaction effect between the agricultural machinery components and soil, and the accuracy of the soil model and contact parameters have a significant impact on the simulation results [9,10].
Currently, much research, both domestic and foreign, has been conducted on soil parameter calibration and soil–machinery interactions based on the discrete element method [11,12,13,14]. Aikins et al. developed a high-viscosity soil model by combining the hysteretic spring and linear cohesion models. The static and rolling coefficients of friction were calibrated, and the accuracy of the model was verified through trenching experiments [15]. In a study on the interaction between soil and a plate-type plow using the discrete element method, Ucgul et al. compared the working traction force and the profile of the plowed furrow through simulations and field experiments. The authors revealed the ability of the discrete element method to accurately simulate the actual working processes [16]. Milkevych et al. established a discrete element model for the interaction between the soil and soil-engaging components for weed control operations. The authors analyzed the impact of weed control operations on soil displacement and validated the effectiveness of the discrete element simulation through the consistent results obtained from both simulations and field experiments [17]. Shaikh et al. established a soil–grouser interaction model and simulated the interaction of a single grouser shoe with clay loam terrain at a varying moisture content by the DEM with Hertz–Mindlin contact [18]. Xiang et al. calibrated the contact parameters of a discrete element model for cohesive soil based on soil pile-up experiments and validated the effectiveness of the calibrated parameters through excavation tests [19]. Based on the hysteretic spring and linear cohesion models, Ma Shuai et al. calibrated the contact model parameters between frozen soil particles and validated the accuracy of the model parameters through a comparative analysis of trench and simulation experiments [20]. Zhang et al. constructed a discrete element simulation model of the contact between sticky black soil and agricultural implements based on the Hertz–Mindlin and JKR cohesion models. By combining simulations and trench experiments, the authors evaluated the micro-disturbance mechanism and macro-disturbance state of soil caused by plowing and shoveling [21].
In summary, establishing an accurate soil discrete element particle model is crucial in order to accurately investigate the interaction between agricultural machinery components and soil [22,23,24]. In this study, physical soil characteristic parameters, including moisture content, density, hardness, particle size distribution, and the angle of repose were measured in typical cotton fields within the northwest region of China during the suitable seeding period. A soil discrete element particle model was constructed, and the contact parameters of the model were calibrated based on experiments and discrete element simulations. A bionic hole-forming device for cotton was designed and compared with conventional square and cone-shaped hole-forming devices through simulation experiments. The purpose of this study was to: (1) perform discrete element method and multi-body dynamics (DEM-MBD) modeling to analyze hole-forming processes and performance under various hole-forming devices, and (2) provide optimization approaches for the improvement of hole-forming devices.

2. Materials and Methods

2.1. Experiment Material

Soil samples were collected from a cotton test field located in Tamutugrak Township, Xinhe County, Northwest China, with the geographical coordinates 82°25′26″ E and 41°16′43″ N and an altitude of 952 m. In order to determine the model parameters, the samples were collected using the five-point sampling method from the plowed and harrowed cotton field soil during the suitable sowing period (March to April). The sampling depth was set as 10 cm, as the cotton seeder had an impact within 10 cm of the soil layer.

2.2. Soil Particle Size Distribution

The soil sieving method was employed to determine the particle size distribution of the soil in the cotton fields during the suitable sowing period. A high-frequency vibrating screen (GZS-1, TuoZhan Instrument Equipment Co., Ltd., Taizhou, China) with a 500 times/min vibration frequency and 1.5 mm vibration amplitude was used for the measurements along with different sieve meshes. During the experiment, sieve meshes with pore diameters of 2.0 mm, 1.0 mm, and 0.075 mm were placed on the high-frequency vibrating screen [9]. A certain amount of soil was weighed and placed on the top sieve, the top cover was closed, and the sieve meshes were fixed prior to turning on the high-frequency vibrating screen. The vibrating screen was then run for 10 min and subsequently turned off. The sieve meshes were then removed (from the largest to smallest pore diameter), and the soil samples in each sieve were taken out. The soil samples adhering to the sieve meshes were cleaned with a soft brush, and each soil sample was weighed. The experiment was repeated three times. Table 1 reports the results.
The proportion of gravel (particle size ˃ 1 mm), sand (particle size 0.075–1 mm), and silt (particle size < 0.075 mm) in the soil of the experimental farmland was 34.78%, 52.12%, and 13.1%, respectively (Table 1).

2.3. Soil Moisture Content and Density

The soil moisture content was measured by the constant pressure and temperature drying method using an electric hot air oven (GZX-9140MBE, Bosen Instrument Co., Ltd., Shanghai, China); electronic balance (CL-T500, Chenlong Hengxin Trading Co., Ltd., Beijing, China, 0.001 g accuracy); and vacuum sample drying box. During the experiment, the aluminum boxes of each sample were labeled and weighed. The soil samples were then placed in the aluminum boxes, weighed, and placed in the drying oven with the temperature set at 105 °C. Every two hours, the samples were taken out for weighing until the weight became constant, and the final weight was recorded. The soil moisture content of the farmland was calculated according to Formula (1):
M = m 0 - m 1 m 0 × 100 % ,
where M is the moisture content of the soil (%), m 0 is the mass of the soil before drying (g), and m 1 is the mass of the soil after drying (g).
The soil moisture content test was repeated three times. The results determined that the average moisture content of the soil layer at a 10 cm depth equaled 14.46%.
The soil density of the cotton field during the suitable seeding period was measured with the ring knife method. A soil sampling ring knife (ZLT, Zhonglutong, testing instrument Co., Ltd., Suzhou, China) with a 100 cm3 volume and electronic balance (CL-T500, Chenlong Hengxin Trading Co., Ltd., Beijing, China, 0.001 g accuracy) were used for the measurements. The mass of the aluminum sample box was first weighed using the electronic balance. A 100 cm3 volume of soil was collected using the soil sampling ring knife within the 10 cm soil layer and placed in the aluminum sample box for weighing. The soil density was then calculated according to Formula (2):
ρ s = m s V s ,
where ρ s is the soil density (g/cm3), m s is the weight of the undisturbed soil (g), and V s is the volume of the undisturbed soil (100 cm3).
The soil density measurements were repeated three times. The average soil density in the 10 cm soil layer during the suitable seeding period of the cotton field was determined as 1.37 g/cm3.

2.4. Soil Firmness

A soil firmness tester (LD-TJ, Shandong Lai Ende Intelligent Technology Co., Ltd., 0.1% measurement accuracy) was used to measure the soil firmness of the cotton field during the suitable seeding period. Figure 1 depicts the measurement process. The five-point sampling method was used for the measurements. The testing probe was placed vertically on each of the five testing points. The operation panel data were cleared by selecting a “zero clearing” operation. The probe was pressed into the soil at a depth of 10 cm, and the experimental data were recorded. The average soil firmness was calculated based on the data from the five testing points.
The soil firmness test was repeated three times, and the firmness of the soil in the cotton field during the suitable seeding period was determined as 1547.7 ± 89.6 Pa.

2.5. Soil Stacking Angle

Soil pile experiments were performed to measure the soil pile angle of the cotton field during the suitable seeding period. A funnel (14.0 cm top diameter, 2.6 cm bottom diameter, and 12.5 cm height); stand; material dropping plate; ruler; and inclination angle meter (ZLT, Zhonglutong Testing Instrument Co., Ltd., Suzhou, China, accuracy of 0.2°) were used for the measurements. The height between the funnel mouth and bottom plate was determined by adjusting the position of the bracket. During the experiment, the distance between the funnel mouth and bottom plate was first adjusted to 15 cm, and the soil was evenly poured into the funnel. The soil pile angle was measured with a ruler and an inclination angle meter. Figure 2 presents the measurement process. In order to reduce experimental errors, four test directions were taken for each soil pile angle test, and the average value of the test results was taken as the soil pile angle of that test. The soil pile angle test was repeated three times. Table 2 reports the results.
The minimum and maximum repose static angles of the cotton field soil during the sowing period were determined as 37.45° and 42.60°, respectively, and the soil accumulation as 39.34 ± 1.71°.

2.6. Discrete Element Contact Model of Soil

Soil particles adhere to each other due to the influence of water and other chemical substances, resulting in the phenomenon of aggregation [25]. The Hertz–Mindlin (no-slip) contact model cannot accurately simulate this mechanical phenomenon of soil [26]. However, combining the cohesive force of the Johnson–Kendall–Roberts (JKR) model with the Hertz–Mindlin model can accurately represent the influence of various cohesive effects [27]. By fully considering the influence of the Van der Waals forces among particles within the contact range, it is possible to simulate adhesive and cohesive systems such as fine powders, dry powders, and wet materials [28]. Therefore, the combined Hertz–Mindlin and JKR contact model is generally adopted to simulate the mechanical behavior of soil. In this model, the JKR normal elastic contact force is based on the Johnson–Kendall–Roberts theory, which characterizes the overlap between particles, interaction parameters, and surface energy. The other forces are calculated in the same way as the Hertz–Mindlin (no-slip) contact model. The JKR normal elastic force is described in Formula (3):
F JKR = 4 π γ E eq α 3 2 + 4 E eq 3 R eq α 3 ,
δ = α 2 R eq 4 π γ α E eq   ,
where F JKR is the JKR normal elastic force (N), α is the radius of the contact circle of the two contacting particles (m), δ is the normal distance between the particles (m), and γ is the surface energy (J/m2).
When there is a certain gap between two particles, this contact model can also be used to calculate the cohesive force between the particles. The maximum distance between the two particles with a non-zero cohesive force can be calculated using Formula (5):
δ c = α c 2 R eq 4 π γ α c E eq ,
α c = 9 π γ R eq 2 E eq 3 4 1 2 1 3 ,
where δ c is the maximum distance between particles with non-zero cohesive force (m), and α c is the maximum contact circle radius between two contacting particles (m).
When the normal gap between the particles is greater than δ c , the cohesive force between the particles is zero. If the particles are not in contact with each other and the normal gap between them is less than δ c , the maximum value of the cohesive force can be calculated using Formula (7):
F pullout = 3 2 π γ R eq ,
where F pullout is the maximum value of the cohesive force (N).
The combined Hertz–Mindlin and JKR contact model can provide a frictional force for a large cohesive force component in the contact normal and can accurately simulate the mechanical behavior of soil. As the direct target of the hole-forming process was the soil in a typical cotton field in Northwest China during the suitable sowing period, the combined Hertz–Mindlin and JKR model was selected as the contact model between soil particles in the simulation process.

2.7. Stacking Angle Test

In this study, the contact parameters between soil and soil and between soil and steel were calibrated using simulation experiments. The soil pile angle was simulated and measured. During the simulation experiment, the funnel size and distance between the funnel outlet and steel plate were consistent with the soil pile angle measurement experiment in soil stacking angle test. According to the soil particle size distribution experiment, the soil particle size was mainly concentrated around 1 mm, and the soil particle radius was set as 1 mm in the simulation. The soil particles were randomly distributed within 0.8–1.2 times the standard volume and generated in a planar manner on the funnel. The soil particles fell through the funnel in a free fall, and the total mass of soil particle production was 200 g. The simulation ended when the soil particles were piled up and stabilized on the flat plate. The Protractor module in EDEM was used to measure the simulation static angle of the soil. Figure 3 depicts the measurement process.
We determined the Plackett–Burman screening experimental design using Design-Expert 12.0 to screen the contact parameters that had a significant impact on the soil stacking angle. The soil stacking angle was set as the response value. There were a total of seven parameters to be calibrated in the simulation experiment, labeled as X1 to X7, and four virtual parameters labeled as X8 to X11. Based on a large number of preliminary trials and relevant literature, the following parameter values were set: the Poisson’s ratio of the cotton field soil was 0.36; the shear modulus was 1.0 × 106 Pa; the soil-to-soil collision restoration coefficient was 0.2–0.6; the soil-to-soil static friction coefficient was 0.2–0.6; the soil-to-soil dynamic friction coefficient was 0.05–0.35; the soil-to-steel collision restoration coefficient was 0.3–0.7; the soil-to-steel static friction coefficient was 0.2–0.6; the soil-to-steel dynamic friction coefficient was 0.05–0.25; and the JKR surface energy was 0.05–0.35 J/m2. The Plackett–Burman experimental parameter design list is shown in Table 3.
Table 4 reports the design plan and simulation results of the Plackett–Burman screening experiment for the soil simulation parameters. The impacts of the calibrated parameters were obtained through variance analysis (Table 5). Compared with other parameters, the soil-to-soil static friction coefficient X2, soil-to-soil dynamic friction coefficient X3, soil-to-steel plate static friction coefficient X5, and JKR surface energy X7 had a significant impact on the soil simulation repose angle, while the remaining parameters had no significant effect on the soil pile angle. Therefore, in the curved surface response test, the optimization of the soil-to-soil static friction coefficient X2, soil-to-soil dynamic friction coefficient X3, soil-to-steel plate static friction coefficient X5, and JKR surface energy X7 was sufficient.
Based on the results of the screening experiment, with the soil pile angle as the indicator, the significant parameters including the soil-to-soil static friction coefficient X2, soil-to-soil dynamic friction coefficient X3, soil-to-steel plate static friction coefficient X5, and JKR surface energy X7 were gradually increased according to the selected step size, and the remaining parameters were simulated at the intermediate level. Table 6 reports the simulation test factor encoding.
Based on the factor encoding of the simulation test shown in Table 6, a Box–Behnken experiment of the surface response was designed. A total of 30 tests were conducted. Table 7 reports the simulation test plan and results.

2.8. Design of Novel Bionic Hole-Forming Device

The oriental mole cricket has an extremely strong soil digging ability. Its forelegs are nail-shaped (Figure 4). The wedge angle of its claws reduces resistance during the soil digging process [29]. In this study, the outer contour of the foreleg claws of the oriental mole cricket was employed to design a new type of bionic hole-forming device.
The oriental mole cricket has forelegs with four claws, where the outermost claw (claw toe 1 in Figure 4) is well-developed. As such, we used the contour fitting curve of foreleg claw 1 as the bionic prototype [30]. The image analysis tools in Matlab R2018b (MathsWork) were employed to analyze the structure of the oriental mole cricket foreleg claws (Figure 5). In particular, grayscale processing was performed on an image of foreleg claw toe 1, which was captured from Figure 4, followed by binary processing, and the contour of the binary image was then extracted to obtain contour data points. The contour curve was drawn based on these data points.
The outer contour of the lateral claw toe frequently comes into contact with the soil. The design prototype of the hole-forming device was based on the outer contour. The X and Y two-dimensional coordinate data of the obtained outer contour curve were imported into the Curve Fitting Toolbox of Matlab, and polynomial fitting was performed on the data points. By selecting a third-order polynomial model for fitting, the outer contour of the first claw toe was fitted to a curve (Figure 6).
The equation for the fitted curve model is given by Formula (8):
f x = p 1 x 3 + p 2 x 2 + p 3 x + p 4 ,
where the coefficients of p 1 , p 2 , p 3 , and p 4 are 7.144 × 10−7, 2.328 × 10−4, −0.1561, and 34.51, respectively. The sum of the squares of error (SSE) of the fitting equation was 8.032 × 104, the coefficient of determination (R-square) was 0.9992, and the adjusted coefficient of determination (adjusted R-square) was also 0.9992. Furthermore, the root mean square error (RMSE) was 1.077. Thus, the sum of the squares of error of the polynomial fitting model was close to zero, and both the coefficient of determination and the adjusted coefficient of determination were close to 1. This indicated that the model fit the data points well, and the model equation could represent the contour curve of the oriental mole cricket.
The equation for the bionic curve was input into the 3D modeling software Solidworks 2018. Figure 7 presents the resulting design of the three-dimensional model for the bionic hole seeder fixed hole-forming mechanism. The length, width, and thickness of the fixed hole-forming mechanism were 74 mm, 36.3 mm, and 30 mm. The inclination angle was 86°. The structural parameters of the hole-forming mechanism mainly affected the hole size of hole sowing. The length and thickness of the fixed hole-forming mechanism determined the depth and width of the hole.

2.9. Interaction Model between Hole-Forming Device and Soil

Currently, square and cone-shaped hole-forming devices are commonly used in cotton hole seeders due to their simple structure and manufacturing process and stress condition analysis. In this study, a comparative experiment was conducted between the bionic hole seeder and traditional square and cone-shaped hole-forming devices to eliminate the influence of dimensional parameters on the experimental results. The design dimensions of the three types of hole seeders were essentially consistent. Figure 8 depicts the corresponding three-dimensional models.
We adopted the discrete element method and multi-body dynamics (DEM-MBD) method to analyze the interaction mechanism between the main soil-contacting component of the seeder, the hole-forming device, and the soil. First, the hole seeder was modeled in Solidworks, and the three-dimensional model was then imported into the multibody dynamics software Recurdyn V9R2 (FunctionBay). The material properties, connections, forces, motions, and contact parameters were set in the software. The dynamic ring, dynamic disc, and fixed hole-forming device in the model were fixed to each other and rotated along the central axis. The central axis moved forward at a certain speed, and the dynamic hole-forming device rotated about the axis of the fixed hole-forming device with a spring connection between them (Figure 9). After all parameters were set, motion simulations were performed to ensure the correct motion before importing the model into EDEM 2020 (DEM-Solutions).
In order to facilitate the analysis of the soil disturbance induced by the hole-forming device, the soil in the simulation soil trough was layered according to the main soil layer subjected to the hole-forming device. In particular, the soil was divided into three layers, namely, the upper, middle, and lower layers. The length, width, and height of the simulated soil trough were 600 mm, 150 mm, and 50 mm, respectively, and the contact parameters and size distribution of soil particles in each layer were consistent with the soil parameters calibrated in Section 2.2. The mechanism did not move during the generation of soil particles in the soil trough. The bottom layer of soil was generated first, followed by the middle and upper layers. Approximately 460,000 particles were generated. Figure 10 presents the simulation operating state of the hole-forming device.

3. Results and Discussion

3.1. Stacking Angle Test

Table 8 reports the significance test results of the regression model. The fitting of the regression model was highly significant (p < 0.01). p-values less than 0.01 were determined for the first-order terms (X2, X5, and X7) of the soil-to-soil static friction coefficient X2, soil-to-steel static friction coefficient X5, and JKR surface energy X7, as well as the second-order term (X52) of the soil-to-steel static friction coefficient X5. This indicated that they had a very significant effect on the soil simulation repose angle. Moreover, the first-order term (X3) of the soil-to-soil dynamic friction coefficient X3 had a p-value less than 0.05, indicating a significant effect on the soil simulation pile angle. The p-value of the lack-of-fit term was 0.1569, which was not significant, indicating that the regression model did not have any other major factors affecting the index. The goodness of fit R2 of the regression equation was 0.9373, revealing a strong fit between the predicted and actual values, and a high degree of explanation between the independent and dependent variable.
Based on the variance analysis results of the regression model, the insignificant factors affecting the soil simulation repose angle in Table 8 were removed. The regression equation for the soil simulation repose angle could be obtained using Formula (9):
Y = 41.86 − 0.93X2 + 0.84X3 + 1.49X5 + 3.65X7 − 0.78X52
With the target of the measured soil repose angle of 39.37°, the soil simulation repose angle regression model was optimized and solved. The optimal solution was obtained when the soil-to-soil static friction coefficient X2 was 0.49, the soil-to-soil dynamic friction coefficient X3 was 0.26, the soil-to-steel static friction coefficient X5 was 0.46, and the JKR surface energy X7 was 0.14 J/m2. Figure 11 depicts the simulation results based on this parameter combination. Comparing the simulation soil pile angle and measured values revealed that they essentially had the same shape.

3.2. Hole-Forming Test

In order to analyze the disturbance of different types of hole seeders on the soil, the soil cross-sections at the bottom of the hole following the seeder operation were clipped using the Clipping tool in the EDEM Analyst template. The cross-sections were taken along the X-axis for the top view and along the Z-axis for the side view (Figure 12 and Figure 13, respectively). The movement range of soil particles in both views caused by the square hole-forming device was larger than that of the other devices, and soil particles from the upper and middle layers were brought into the lower layer. The cone-shaped and bionic hole-forming devices disturbed the soil less, bringing fewer particles from the upper and middle layers into the lower layer. From the top view, it can be seen that the soil disturbance caused by the square and cone-shaped hole-forming devices formed a trapezoidal shape. Moreover, the disturbance range caused by the cone-shaped hole-forming device was smaller than that of the square hole-forming device. The disturbance caused by the bionic hole-forming device formed a triangular shape and created the smallest disturbance range. From the side view, it can be seen that the hole-forming device of the seeder created holes in the soil, and the square and cone-shaped hole-forming devices formed soil ridges behind the holes. However, the bionic hole-forming device hardly formed any soil ridges.
In order to investigate the effect of different types of hole-forming devices on the movement of soil particles, the velocity vector of soil particle groups under the action of different types of hole-forming device was analyzed. Figure 14 presents the velocity vector map of soil particles for hole seeders with different hole-forming devices at 0.44 s. In the figure, the dynamic hole-forming device is about to open, as the soil pressure exerted on it is greater than the pre-tension force of the spring. The maximum velocity of soil particles caused by the square hole-forming device was the highest, reaching 0.386 m/s, followed by the cone-shaped hole-forming device at 0.301 m/s, and the bionic hole-forming device had the lowest velocity at 0.298 m/s. As the contact surface of the square hole-forming device with the soil was a plane, the velocity direction of the soil particles was essentially perpendicular to the contact surface. The lower layers of the soil had a larger linear velocity caused by contact with the hole-forming device, resulting in a larger disturbance to the soil. A similar phenomenon also existed for the cone-shaped hole-forming device, but the disturbance range on the soil was smaller than that of the square hole-forming device. As for the bionic hole-forming device, the contact surface with the soil was a curved contour line, and the velocity direction of the soil particles was along the normal direction of the contour line. This resulted in a more uniform distribution of velocity compared to the other devices. Thus, the simulation results showed that the bionic hole-forming device caused the least disturbance on the soil.
Figure 15 presents the velocity vector plot of soil particles for hole seeders with different hole-forming devices at 0.6 s. In the figure, the hole-forming device is beginning to leave the soil as the seeder rotates, and the dynamic hole-forming device starts to close. The maximum velocity of the soil particles caused by the square hole-forming device was the highest, reaching 0.402 m/s, followed by the cone-shaped hole-forming device at 0.349 m/s, and the bionic hole-forming device had the lowest velocity at 0.306 m/s. The movement of the fixed square hole-forming device followed a parabolic trend as the seeder rotated, causing the soil particles to also move backwards in a parabolic trend. As a result, ridges were formed in the soil behind the hole after the operation. The cone-shaped hole-forming device also caused ridges in the soil, yet the ridges induced by the square hole-forming device were larger, as its contact surface with the soil was larger. Although the fixed bionic hole-forming device moved in a parabolic trend, its contact surface with the soil was formed by a curved contour line and did not cause the soil particles to be propelled backwards. As a result, there were barely any soil ridges created behind the hole in the soil. Thus, the simulation results revealed that the bionic hole-forming device caused less disturbance on the soil and did not create any soil ridges.
Figure 16 presents the maximum force on different hole openers during the hole-forming operation. The force exerted on the hole-forming device gradually increased with the depth of insertion into the soil. When the top of the hole-forming device reached the lowest point in the soil, the force exerted on the hole opener reached its maximum value. Following this, the force exerted on the hole-forming device began to decrease until it completely left the soil. The square hole-forming device experienced more force than the cone-shaped and bionic hole-forming devices, with a maximum force of 9.50 N throughout the entire operation. This could be attributed to the similar upper and lower widths of the square hole-forming device, resulting in a larger contact surface with the soil compared to the cone-shaped and bionic hole-forming devices. This led to a larger force being exerted. Furthermore, the cone-shaped and bionic hole-forming devices had a larger upper width and smaller lower width, resulting in a smaller contact surface area, and hence less force was exerted. The maximum force exerted on the cone-shaped hole-forming device during the entire operation was 7.96 N, and that on the bionic hole-forming device was 7.51 N. The analysis showed that although the bionic hole-forming device had a larger contact surface with the soil than the cone-shaped hole-forming device, its maximum force during the hole-forming device process was smaller. This indicated a disturbance and resistance reduction effect on the soil during the operation. The soil moisture content is a critical factor that influences the performance of hole-forming devices. An increase in moisture content weakens the cohesion between soil particles, resulting in a reduction in the shear and compressive strength of the soil. Moreover, a high moisture content increases the fluidity of the soil and reduces hole formation performance. On the other hand, a low moisture content raises the resistance faced during the sowing process. Therefore, planting under optimal soil conditions is advantageous for enhancing the quality of sowing.

4. Conclusions

A typical soil discrete element particle model was established, and a bionic hole-forming device was designed. The hole-forming processes and performance of various hole-forming devices were investigated by DEM-MBD simulation tests. The mechanism of motion interaction between the soil and hole-forming device was analyzed, and the following conclusions could be drawn:
(1)
The moisture content of the soil in the cotton field during the appropriate sowing period was 14.46%, the density was 1.37 g/cm3, and the bulk density was 1547.7 Pa. Gravel (particle diameter ˃ 1–2 mm); sand (particle diameter 0.075–1 mm); and silt (particle diameter < 0.075 mm) accounted for 34.78%, 52.12%, and 13.1%, respectively, of the particle size distribution of the soil. The natural rest angle of the soil was 39.34°.
(2)
Based on physical experiments and simulation experiments, the combined Hertz–Mindlin and JKR contact model was used to establish a second-order regression model for the soil accumulation angle and soil contact parameters through optimized experimental methods. The target optimization was based on the measured soil accumulation angle of 39.37°. The optimal solution was obtained when the static friction coefficient between soil and soil was 0.49, the dynamic friction coefficient between soil and soil was 0.26, the static friction coefficient between soil and steel was 0.46, and the surface energy of JKR was 0.14 J/m2.
(3)
A bionic hole-forming device was designed based on bionic technology using the oriental mole cricket as a prototype. A discrete element method (DEM) and multi-body dynamics (MBD) coupled algorithm was used to establish a discrete element simulation model of the hole-forming device and soil. A hole-formation experiment was carried out using an EDEM–Recurdyn joint simulation. The interaction between different structural forms of the hole-forming device and soil during the hole-formation process was analyzed. The simulation results showed that the disturbance of the soil caused by the bionic hole-forming device was small, and it did not cause soil ridges to form. Moreover, the soil resistance of the hole-forming device was the smallest, at 7.51 N.
This work should be helpful in improving hole-forming performance and reducing soil resistance during the cotton planting process. However, there was still some discrepancy between the results of the simulation experiment and the actual field operation. Future work will explore the performance of hole-forming devices and the interaction mechanism between hole-forming devices and soil under different soil conditions in field experiments.

Author Contributions

Conceptualization, L.W. and J.X.; methodology, L.W. and J.X.; software, L.W.; validation, L.W., X.W. and J.X.; formal analysis, X.L.; investigation, W.G. and X.H.; data curation, S.H.; writing—original draft preparation, L.W.; writing—review and editing, L.W. and J.X.; project administration, L.W.; funding acquisition, L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Bintuan Science and Technology Program (2022CB001-06 and 2023AB005) and the President Fund from Tarim University (TDZKYS202302).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

The authors would like to thank their schools and colleges, as well as the funding providers of the project. All support and assistance are sincerely appreciated.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, L.; He, X.W.; Hu, C.; Guo, W.S.; Wang, X.F.; Xing, J.F. Measurement of physical parameters of coated cotton seed and parameter calibration of discrete element. Chin. Agric. Univ. 2022, 27, 71–82. [Google Scholar]
  2. Wang, J.K.; Guo, K.Q.; Lu, X.M.; Fu, W.; Qi, C.S. Experiment and improvement on clamping cotton precision seeding dibbler. Chin. Soc. Agric. Mach. 2011, 42, 43–47. [Google Scholar]
  3. Aikins, K.A.; Ucgul, M.; Barr, J.B.; Awuah, E.; Antille, D.L.; Jensen, T.A.; Desbiolles, J.M.A. Review of discrete element method simulations of soil tillage and furrow opening. Agriculture 2023, 13, 541. [Google Scholar] [CrossRef]
  4. Sun, J.B.; Liu, Q.; Yang, F.Z.; Liu, Z.J.; Wang, Z. Calibration of discrete element simulation parameters of sloping soilon loess plateau and its interaction with rotary tillage components. Chin. Soc. Agric. Mach. 2022, 53, 63–73. [Google Scholar]
  5. Song, Z.H.; Li, H.; Yan, Y.F.; Tian, F.Y.; Li, Y.D.; Li, F.D. Calibration method of contact characteristic parameters of soil in mulberry field based on unequal-diameter particles dem theory. Chin. Soc. Agric. Mach. 2022, 53, 21–33. [Google Scholar]
  6. Zeng, Z.W.; Ma, X.; Cao, X.L.; Li, Z.H.; Wang, X.C. Critical review of applications of discrete element method in agricultural engineering. Chin. Soc. Agric. Mach. 2021, 52, 1–20. [Google Scholar]
  7. Dong, X.Q.; Su, C.; Zheng, H.N.; Han, R.Q.; Li, Y.L.; Wan, L.P.C.; Song, J.N.; Wang, J.C. Analysis of soil disturbance process by vibrating subsoiling based on DEM-MBD coupling algorithm. Trans. Chin. Soc. Agric. Eng. 2022, 38, 34–43. [Google Scholar]
  8. Makange, N.R.; Ji, C.Y.; Torotwa, I. Prediction of cutting forces and soil behavior with discrete element simulation. Comput. Electron. Agric. 2020, 179, 105848. [Google Scholar] [CrossRef]
  9. Song, S.L.; Tang, Z.H.; Zheng, X.; Liu, J.B.; Meng, X.J.; Liang, Y.C. Calibration of the discrete element parameters for the soil model of cotton field after plowing in Xinjiang of China. Trans. Chin. Soc. Agric. Eng. 2021, 37, 63–70. [Google Scholar]
  10. Wang, X.L.; Zhong, X.K.; Geng, Y.L.; Wei, Z.C.; Hu, H.; Geng, D.Y.; Zhang, X.C. Construction and parameter calibration of the nonlinear elastoplastic discrete element model for no-tillage soil compaction. Trans. Chin. Soc. Agric. Eng. 2021, 37, 100–107. [Google Scholar]
  11. Zhao, H.B.; Huang, Y.X.; Liu, Z.D.; Liu, W.Z.; Zheng, Z.D. Applications of discrete element method in the research of agricultural machinery: A review. Agriculture 2021, 11, 425. [Google Scholar] [CrossRef]
  12. Ucgul, M. Simulating soil-disc plough interaction using discrete element method-multi-body dynamic coupling. Agriculture 2023, 13, 305. [Google Scholar] [CrossRef]
  13. Sadek, M.A.; Chen, Y.; Zeng, Z. Draft force prediction for a high-speed disc implement using discrete element modelling. Biosyst. Eng. 2021, 202, 133–141. [Google Scholar] [CrossRef]
  14. Asaf, Z.; Rubinstein, D.; Shmulevich, I. Determination of discrete element model parameters required for soil tillage. Soil Tillage Res. 2007, 92, 227–242. [Google Scholar] [CrossRef]
  15. Aikins, K.A.; Ucgul, M.; Barr, J.B.; Jensen, T.A.; Antille, D.L.; Desbiolles, J.M.A. Determination of discrete element model parameters for a cohesive soil and validation through narrow point opener performance analysis. Soil Tillage Res. 2021, 213, 105123. [Google Scholar] [CrossRef]
  16. Ucgul, M.; Saunders, C. Simulation of tillage forces and furrow profile during soil-mouldboard plough interaction using discrete element modelling. Biosyst. Eng. 2020, 190, 58–70. [Google Scholar] [CrossRef]
  17. Milkevych, V.; Munkholm, L.J.; Chen, Y.; Nyord, T. Modelling approach for soil displacement in tillage using discrete element method. Soil Tillage Res. 2018, 183, 60–71. [Google Scholar] [CrossRef]
  18. Shaikh, S.A.; Li, Y.M.; Ma, Z.; Chandio, F.A.; Tunio, M.H.; Liang, Z.W.; Solangi, K.A. Discrete element method (DEM) simulation of single grouser shoe-soil interaction at varied moisture contents. Comput. Electron. Agric. 2021, 191, 106538. [Google Scholar] [CrossRef]
  19. Xiang, W.; Wu, M.L.; Lü, J.N.; Quan, W.; Ma, L.; Liu, J.J. Calibration of simulation physical parameters of clay loam based on soil accumulation test. Trans. Chin. Soc. Agric. Eng. 2019, 35, 116–123. [Google Scholar]
  20. Ma, S.; Xu, L.M.; Yuan, Q.C.; Niu, C.; Zeng, J.; Chen, C.; Wang, S.S.; Yuan, X.T. Calibration of discrete element simulation parameters of grapevine antifreezing soil and its interaction with soil-cleaning components. Trans. Chin. Soc. Agric. Eng. 2020, 36, 40–49. [Google Scholar]
  21. Zhang, C.; Fan, X.H.; Li, M.S.; Li, G.; Zhao, C.K.; Sun, W.L. Simulation analysis and experiment of soil disturbance by chiesel plow based on edem. Chin. Soc. Agric. Mach. 2022, 53, 52–59. [Google Scholar]
  22. Ucgul, M.; Fielke, J.M.; Saunders, C. Three-dimensional discrete element modelling of tillage: Determination of a suitable contact model and parameters for a cohesionless soil. Biosyst. Eng. 2014, 121, 105–117. [Google Scholar] [CrossRef]
  23. Hossein, R.K.; Mahmood, M.E.; Hadi, A.N. Simulation of Draught Force During Chisel Ploughing Using Discrete Element Method. J. Biosyst. Eng. 2022, 47, 152–166. [Google Scholar]
  24. Zheng, K.; He, J.; Li, H.W.; Diao, P.S.; Wang, Q.J.; Zhao, H.B. Research on polyline soil-breaking blade subsoiler based on subsoiling soil model using discrete element method. Chin. Soc. Agric. Mach. 2016, 47, 62–72. [Google Scholar]
  25. Wang, X.L.; Hu, H.; Wang, Q.J.; Li, H.W.; He, J.; Chen, W.Z. Calibration method of soil contact characteristic parameters based on DEM theory. Chin. Soc. Agric. Mach. 2017, 48, 78–85. [Google Scholar]
  26. Wang, X.; Zhang, S.; Pan, H.; Zheng, Z.; Huang, Y.; Zhu, R. Effect of soil particle size on soil-subsoiler interactions using the discrete element method simulations. Biosyst. Eng. 2019, 182, 138–150. [Google Scholar] [CrossRef]
  27. Du, J.; Heng, Y.; Zheng, K.; Luo, C.; Zhu, Y.; Zhang, J.; Xia, J. Investigation of the burial and mixing performance of a rotary tiller using discrete element method. Soil Tillage Res. 2022, 220, 105349. [Google Scholar] [CrossRef]
  28. Qiu, Y.Q.; Guo, Z.J.; Jin, X.; Zhang, P.G.; Si, S.J.; Guo, F.G. Calibration and verification test of cinnamon soil simulation parameters based on discrete element method. Agriculture 2022, 12, 1082. [Google Scholar] [CrossRef]
  29. Xiao, M.H.; Wang, K.X.; Yang, W.; Wang, W.C.; Jiang, F. Design and experiment of bionic rotary blade based on claw toe of cryllotalpa orientalis burmeister. Chin. Soc. Agric. Mach. 2021, 52, 55–63. [Google Scholar]
  30. Zou, X.X.; Liao, Y.L.; Ma, Q.F.; Liu, S.H.; Zhu, D.Y. Bionic design and finite element analysis for digging shovel of cassava harvester. Sci. Technol. Eng. 2013, 13, 122–128. [Google Scholar]
Figure 1. Firmness determination of soil.
Figure 1. Firmness determination of soil.
Agriculture 13 01421 g001
Figure 2. Stacking angle determination of soil.
Figure 2. Stacking angle determination of soil.
Agriculture 13 01421 g002
Figure 3. Determination of simulated soil stacking angle.
Figure 3. Determination of simulated soil stacking angle.
Agriculture 13 01421 g003
Figure 4. Forefoot outline of oriental mole cricket.
Figure 4. Forefoot outline of oriental mole cricket.
Agriculture 13 01421 g004
Figure 5. Process of extracting contours: (a) grayscale image; (b) binary image; (c) contour extraction; (d) contour curve.
Figure 5. Process of extracting contours: (a) grayscale image; (b) binary image; (c) contour extraction; (d) contour curve.
Agriculture 13 01421 g005
Figure 6. External contour fitting curve of claw toe.
Figure 6. External contour fitting curve of claw toe.
Agriculture 13 01421 g006
Figure 7. Soil contact components of bionic fixed hole-forming device: (a) front view; (b) lateral view; (c) isometric view.
Figure 7. Soil contact components of bionic fixed hole-forming device: (a) front view; (b) lateral view; (c) isometric view.
Agriculture 13 01421 g007
Figure 8. Soil contact components of bionic fixed hole-forming device: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device.
Figure 8. Soil contact components of bionic fixed hole-forming device: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device.
Agriculture 13 01421 g008
Figure 9. RecurDyn simulation model.
Figure 9. RecurDyn simulation model.
Agriculture 13 01421 g009
Figure 10. Simulation operation status of the hole-forming device: 1—upper soil; 2—middle soil; 3—subsoil.
Figure 10. Simulation operation status of the hole-forming device: 1—upper soil; 2—middle soil; 3—subsoil.
Agriculture 13 01421 g010
Figure 11. Comparison of simulation and measured values of soil stacking angle: (a) simulation test; (b) actual test.
Figure 11. Comparison of simulation and measured values of soil stacking angle: (a) simulation test; (b) actual test.
Agriculture 13 01421 g011
Figure 12. Front view of soil particle state after hole-forming operation: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device. Blue particles are upper soil, Pink particles are middle soil, yellow particles are subsoil. The regions between dashed lines represent the disturbed soil particles.
Figure 12. Front view of soil particle state after hole-forming operation: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device. Blue particles are upper soil, Pink particles are middle soil, yellow particles are subsoil. The regions between dashed lines represent the disturbed soil particles.
Agriculture 13 01421 g012
Figure 13. Lateral view of soil particle state after hole-forming operation: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device. Blue particles are upper soil, pink particles are middle soil, yellow particles are subsoil. The regions between dashed lines represent the disturbed soil particles.
Figure 13. Lateral view of soil particle state after hole-forming operation: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device. Blue particles are upper soil, pink particles are middle soil, yellow particles are subsoil. The regions between dashed lines represent the disturbed soil particles.
Agriculture 13 01421 g013
Figure 14. Vector plot of soil particle velocity at 0.44 s: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device.
Figure 14. Vector plot of soil particle velocity at 0.44 s: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device.
Agriculture 13 01421 g014
Figure 15. Vector plot of soil particle velocity at 0.6 s: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device.
Figure 15. Vector plot of soil particle velocity at 0.6 s: (a) square hole-forming device; (b) cone-shaped hole-forming device; and (c) bionic hole-forming device.
Agriculture 13 01421 g015
Figure 16. Maximum stress of different hole-forming devices.
Figure 16. Maximum stress of different hole-forming devices.
Agriculture 13 01421 g016
Table 1. Particle size distribution of soil.
Table 1. Particle size distribution of soil.
SamplesParticle Size Ratio of Soil (%)
Gravel (>1 mm)Sand (0.075–1 mm)Silt (<0.075 mm)
134.1252.2613.62
235.2451.6713.09
334.9852.4312.59
Average34.7852.1213.10
Table 2. Soil stacking angle measurement results.
Table 2. Soil stacking angle measurement results.
No.Stacking Angle of Soil (°)
Orientation 1Orientation 2Orientation 3Orientation 4AverageTotal Average
139.3038.2537.6040.5538.9339.34
239.2542.6040.5037.4539.95
338.6039.4541.8536.7039.15
Table 3. List of Plackett–Burman test parameters.
Table 3. List of Plackett–Burman test parameters.
SymbolParameterParameter Level
−101
X1Soil-to-soil collision restoration coefficient 0.20.40.6
X2Soil-to-soil static friction coefficient 0.20.40.6
X3Soil-to-soil dynamic friction coefficient 0.050.20.35
X4Soil-to-steel collision restoration coefficient 0.30.50.7
X5Soil-to-steel static friction coefficient 0.20.40.6
X6Soil-to-steel dynamic friction coefficient 0.050.150.25
X7JKR surface energy (J/m2)0.050.250.35
X8X11Virtual parameters
Table 4. Schemes and results of Plackett–Burman test.
Table 4. Schemes and results of Plackett–Burman test.
No.X1X2X3X4X5X6X7X8X9X10X11Stacking Angle of Soil (°)
111−1111−1−1−11−128.57
2−111−1111−1−1−1151.90
31−111−1111−1−1−150.16
4−11−111−1111−1−145.32
5−1−11−111−1111−137.19
6−1−1−11−111−111142.74
71−1−1−11−111−11149.88
811−1−1−11−111−1123.18
9111−1−1−11−111−145.17
10−1111−1−1−11−11126.42
111−1111−1−1−11−1140.51
12−1−1−1−1−1−1−1−1−1−1−124.43
Table 5. Parameter significance analysis of Plackett–Burman test.
Table 5. Parameter significance analysis of Plackett–Burman test.
ParameterDegrees of FreedomMean Square SumF-Valuep-ValueSignificance
X117.476.660.0613
X2149.4144.020.0027**
X31115.51102.910.0005**
X410.32340.28810.6199
X51141.93126.450.0004**
X610.33670.30000.6130
X71916.48816.51<0.0001**
Note: ** denotes an extremely significant impact from the parameter (p < 0.01).
Table 6. Coding of simulation test factors.
Table 6. Coding of simulation test factors.
CodeX2X3X5X7
−20.20.050.20.05
−10.30.1250.30.125
00.40.200.40.20
10.50.2750.50.275
20.60.350.60.35
Table 7. Box–Behnken experimental design and results.
Table 7. Box–Behnken experimental design and results.
No.X2X3X5X7Stacking Angle of Soil (°)
1−1−11145.09
2002043.30
3111144.73
4020043.75
5000040.99
6111−140.47
7000041.91
8000041.88
9−1−1−1−136.09
100−20038.60
11−111147.94
12000042.41
13000041.66
14−11−1−138.59
15200038.00
16000−234.61
1711−1143.64
1811−1−135.24
19000249.71
20−200044.52
211−11−136.48
2200−2034.04
23−111−140.18
241−11144.57
251−1−1−135.17
26−1−1−1144.56
27−11−1143.04
281−1−1144.04
29000043.89
30−1−11−138.05
Table 8. Analysis of variance for regression model.
Table 8. Analysis of variance for regression model.
Source of VarianceSum of SquaresDegrees of FreedomMean Square SumF-Valuep-ValueSignificance
Model438.28171431.305816.0279<0.0001**
X220.6091120.609110.55140.0054**
X316.8002116.80038.60140.0103*
X552.9848152.984827.12710.0001**
X7319.30221319.3022163.4758<0.0001**
X2X30.286210.28620.14650.7072
X2X50.042010.04200.02150.8853
X2X70.225610.22560.11550.7387
X3X54.494414.49442.30100.1501
X3X73.6113.611.84820.1941
X5X70.577610.57760.29570.5946
X220.752410.75240.38520.5441
X320.957910.95790.49040.4945
X5218.1350118.13509.28470.0082**
X720.096710.09670.04950.8269
Residual29.2981151.9532
Lack of fit24.4910102.44912.54740.1569
Pure error4.807150.9614
Cor total467.579816729
Note: ** denotes an extremely significant impact from the parameter (p < 0.01), and * represents a significant impact from the parameter (p < 0.05).
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

Wang, L.; Xing, J.; He, X.; Li, X.; Guo, W.; Wang, X.; Hou, S. Study on the Mechanism of Motion Interaction between Soil and a Bionic Hole-Forming Device. Agriculture 2023, 13, 1421. https://doi.org/10.3390/agriculture13071421

AMA Style

Wang L, Xing J, He X, Li X, Guo W, Wang X, Hou S. Study on the Mechanism of Motion Interaction between Soil and a Bionic Hole-Forming Device. Agriculture. 2023; 13(7):1421. https://doi.org/10.3390/agriculture13071421

Chicago/Turabian Style

Wang, Long, Jianfei Xing, Xiaowei He, Xin Li, Wensong Guo, Xufeng Wang, and Shulin Hou. 2023. "Study on the Mechanism of Motion Interaction between Soil and a Bionic Hole-Forming Device" Agriculture 13, no. 7: 1421. https://doi.org/10.3390/agriculture13071421

APA Style

Wang, L., Xing, J., He, X., Li, X., Guo, W., Wang, X., & Hou, S. (2023). Study on the Mechanism of Motion Interaction between Soil and a Bionic Hole-Forming Device. Agriculture, 13(7), 1421. https://doi.org/10.3390/agriculture13071421

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