Calibration of Model Parameters for Soda Saline Soil-Subsoiling Component Interaction Based on DEM

: To apply the discrete element method (DEM) to simulate the interaction process between soda saline‒alkali soil and subsoiling component in Northeast China, establishing the soda saline‒ alkali soil particle model and selecting more accurate simulation parameters are important. In this paper, we studied the soda saline–alkali soil of the Songnen Plain in China. First, we studied the geometric shape of soda saline–alkali soil particles and proposed a modeling method for single soil particles based on the multisphere combination method. Considering the cohesion of soda saline– alkali soil particles, the Hertz‒Mindlin with JKR (JKR) model was used as the contact model between soil particles. Then, the calibration method was used to obtain simulation parameters of soils that are difficult to obtain experimentally. We conducted soil angle of repose (AoR) tests, the Plackett‒Burman (PB) tests, and steepest ascent (SA) tests in turn to perform a sensitivity analysis for microscopic contact parameters and select the parameters that have a significant effect on the response value (static AoR), i.e., soil surface energy, soil‒soil static friction coefficient, and soil‒soil rolling friction coefficient. Then, the optimal combination of simulation parameters was obtained via the Box‒Behnken (BB) tests, using ANOVA to optimize the multiple regression equation. Finally, the optimal parameter combination was verified by the AoR test and the direct shear (DS) test. The results showed that the parameters had good adaptability for the AoR test. However, the simulation results of the DS test were significantly different from the experimental values. Therefore, the contact model needs to be further modified by adding Bonding bonds between soil particles based on the JKR model and further correcting for Rayleigh time step, shear modulus, and surface energy. By comparing the simulation and the experimental results, it was found that the simulation results obtained from both the DS test and AoR test had relatively small errors relative to physical tests, the two trends are the same, and the values are similar. This verified the feasibility and effectiveness of the soda saline–alkali soil particle modeling method and parameter selection proposed in this paper.


Introduction
The Songnen Plain in Northeastern China, is home to one of the world's three major concentrations of soda saline-alkali soil.Recently, soil saline-alkali desertification is on the rise.Subsoiling technology is widely recognized as an effective engineering solution for improving saline lands.To develop an optimal subsoiling component for soda-saline land improvement, we utilized DEM for simulation.Given the inherent contact between subsoiling components and soil particles during the subsoiling process, it is crucial to establish a more accurate model of soda saline-alkali soil particles when using the DEM [1][2][3][4] to analyze these contact effects and optimize the design of related mechanical components [5,6].
Soil exhibits pronounced anisotropic and nonlinear behavior, characterized by a complex constitutive relationship.Its properties are profoundly influenced by factors such as water and organic ma er content, porosity, and texture.A comprehensive review of the literature underscores the abundance of research dedicated to the physicochemical properties of soda saline-alkali soils [7][8][9][10][11][12], yet investigations into their physical and mechanical properties are notably scarce.Given the unique physicochemical a ributes of soda saline-alkali soil, distinct from those of typical soil, further exploration into how these properties influence physical and mechanical characteristics is imperative.The discrete element method (DEM) is a numerical method for analyzing and solving the dynamic problems of complex discrete systems.It is widely used in the simulation of particle structures [13].The discrete element method can integrate the contact mechanics model of soil particles to simulate the microscopic and macroscopic deformation of particles and the interaction between soil particles and machinery [14].Because the microcontact parameters between soil particles cannot be directly measured by experiments, they need to be obtained through repeated simulation experiments.When utilizing the Discrete Element Method (DEM) to simulate and analyze the movement of soil particle assemblies, it is necessary to calibrate the contact parameters of soil particles such as restitution coefficient, static friction coefficient, and rolling friction coefficient [15].
Many scholars have studied the calibration of soil particle contact parameters using different contact models.The JKR model is suitable for the simulation between wet materials and easily bonded particles, and the cohesion in the model can be er characterize the bonding characteristics of soil particles.Zhu, Li, Xing, Hao, et al. used a JKR contact model to calibrate the contact parameters between soil particles with the angle of repose as the response value [16][17][18][19].The Her -Mindlin with Bonding (HMB) model binds two particles together by "bond", which can characterize the cohesion between soils.Du et al. used HMB and JKR contact models to conduct comparative tests, combined with soil bin test results, to study the mixing performance of rotary tillers and spiral horizontal blade rotary tillers in silty clay loam [20].Hang et al. used the HMB contact model to establish a layered soil model to study the soil disturbance during the subsoiling process [21].
Within DEM analyses, soil particles are conventionally defined as spheres due to the computational simplicity of calculating contacts between spherical particles, as it enhances the simulation efficiency [22].Song, Li, et al. defined soil particles using a singlesphere model and established a soil bin model to analyze soil disturbance during tillage [23].However, the real-world size and shape of soil particles are inherently more intricate; the relationship between the macroscopic and microscopic characteristics of soil particle groups has yet to be fully understood [24,25].The equal diameter soil particles cannot fully reflect the influence of the particle size distribution of the actual soil on the interaction between the soil particles and the working parts [26].It is imperative to recognize that spherical particles alone cannot adequately characterize the accurate shape of soil particles.Song et al. calibrated and optimized the bonding parameters between soil particles by taking the non-equal-diameter soil particles in a mulberry field as the research object.[27].Zhang et al. established four soil particle models with different sizes and shapes for simulation and verified the accuracy of the simulation model [28].
At present, there are few studies on the particles and particle aggregates of soda saline-alkali soil in Northeast China.When DEM is used for a simulation, the contact parameters between the soil particles are worth studying; the meticulous verification of the accuracy of the model also demands thorough analysis and consideration.In this study, we used DEM to calibrate soil particle contact parameters in soda saline soils of the Songnen Plain, Northeast China, using the angle of repose as the response value.To verify the accuracy of the calibration parameters, we conducted the DS test and added the bond between soil particles based on the JKR model to accurately simulate the soil bonding characteristics.Finally, the contact model and contact parameter combination of soil particles were determined, which provided a theoretical basis for the subsequent development of subsoiling machinery in soda saline-alkali land in Northeast China.

Materials and Methods
This paper comprehensively examined and analyzed soda saline-alkali soil, encompassing particle shape, size distribution, and physical and mechanical properties.Recognizing the adhesion characteristics of soil, we employed the Her -Mindlin with JKR (JKR) model as the contact model.First, using soil angle of repose (AoR) tests, we conducted the Placke -Burman (PB), steepest ascent (SA), and Box-Behnken (BB) tests sequentially to select and calibrate the sensitive parameters in the model.Second, we obtained the optimal parameter combination by variance analysis.Finally, to improve the reliability of the verification, two different tests were used to verify the calibrated parameter combination, that is, the AoR test and direct shear (DS) test.Considering the adhesion characteristics of saline-alkali soil and the characteristics of high firmness and easy hardening, bonds are added between soil particles based on the JKR model.Combined with the physical test and simulation results, the contact model and sensitivity parameters are further corrected, in order to provide a reference for the discrete element simulation parameters of soda saline-alkali soil and subsoiling machine in Northeast China.
The soda saline-alkali soil samples were collected on 12 May 2022, at an outdoor temperature of 21 °C.The collection site, situated in Sanye Village, Haituo Township, Daan City, Jilin Province (123.96°E; 45.32° N), is centrally located within the Songnen Plain (see Figure 1).This area is characterized by a temperate semi-humid and semi-arid climate, typical of a continental climate.The soil in this region exhibits an average pH value of 9.5 and an average total salt content of 0.7%, which is heavily salinized soil that has not been cultivated and reclaimed, representing the original landform characteristics of soda saline land under native conditions, as shown in Figure 2.  To comprehensively consider the water and salt transport characteristics within soda saline-alkali soil, three sampling plots were randomly selected, and we employed the fivepoint sampling method.Using a shovel, soil samples were obtained at depths of 0-10 cm, 10-20 cm, 20-30 cm, and 30-40 cm below the surface [29].Each collected soil sample was placed in a self-sealing bag, meticulously numbered for identification purposes, and its wet weight was measured on-site.This process was replicated for a total of 20 samples within each land block.

Soda Saline-Alkali Soil Particle Size and pH Value Test
Soil particle size analysis was conducted using the screening method, which determined the percentage of each particle group within the overall mass of soil particles.To ensure the precise determination of soil pH, we entrusted the soil to a specialized testing agency.Soil pH measurements were performed using the PHSJ-6L pH meter, provided by Shanghai Yidian Scientific Instrument Co., Ltd.(Shanghai, China).

Water Content, Density, and Firmness Test of Soda Saline-Alkali Soil
Moisture content analysis was conducted using the drying method.Soil samples were placed in aluminum boxes, and the drying oven was set to a temperature of 105 °C.The soda saline-alkali soil-filled aluminum boxes were dried for a period of 8 h until a constant weight was achieved.The weight lost during this process represented the quantity of water contained within the sample.
Soil samples were obtained using a ring knife with a volume of V = 100 cm 3 .The ring knife, containing the soil sample, was then weighed using an electronic balance.This procedure was repeated three times for each sample, and the average value was calculated as the density of soda saline-alkali soil.
The TYD-2 hardness tester, manufactured by Zhejiang Top Instrument Co., Ltd.(Hangzhou, Zhejiang, China)., was vertically inserted into the soil.Soil firmness measurements for the 0-40 cm soil layer were obtained by recording the readings.

Simulation Contact Model
Considering the distinctive a ributes of soda saline-alkali soil in Northeast China, characterized by the soil's propensity for bonding and elevated water content, we opted for the JKR model, which is particularly well-suited for small particles and materials exhibiting wet bonding properties, in our simulation analysis.Within the JKR model, the calculation of tangential contact force, normal damping force, and tangential damping force closely parallels that of the Her -Mindlin (no-slip) (HM) contact model.Within the JKR model, the tangential contact force, normal damping force, and tangential damping force are calculated similarly to the Her -Mindlin (no-slip) (HM) contact model, and the cohesive force explains the normal contact force in the JKR theory [30,31].In this paper, the macroscopic bonding characteristics between soil particles are explained by surface energy, and the contact force between particles can be expressed as follows: where FJKR is the JKR normal contact force, N; Ft c , Fn d , and Ft d are the tangential contact force, normal damping force, and tangential damping force, N; δn is the normal overlap of particles, m; δt is the tangential overlap of particles, m; α is the contact radius of particles, m; γ is the surface energy, N/m; vn and vt are the normal and tangential components of the relative velocity of particles, m/s; E * is the equivalent elastic modulus, Pa; and R * is the equivalent contact radius, m.The equivalent elastic modulus, E * , and the equivalent contact radius, R * , are defined as shown in Equations ( 6) and ( 7).
In the equation, E1 and E2 are the elastic moduli of particles 1 and 2, Pa; v1 and v2 are Poisson's ratios of particles 1 and 2; and R1 and R2 are the contact radii of particles 1 and 2, m.

Soil Particle Simulation Model
Since the soda saline-alkali soil in the Songnen Plain of Northeast China has high cohesiveness when the water content is 18%, the soil particles will be bonded into agglomerates.Our observation of the agglomerates formed by soil particles revealed that most of them showed shapes similar to spheres, ellipsoids, and cones, as shown in Figure 3. Consequently, the geometric model for soil particles can be constructed based on these three fundamental shapes, as illustrated in Figure 4.

Soil AoR Test
Due to the strong cohesiveness of soda saline-alkali soil in the Songnen Plain of Northeast China, we used the cylinder pulling method to measure the static angle of repose by a 65 Mn steel cylinder with a diameter of 80 mm and a height of 250 mm.During the test, the cylinder was pulled at a constant speed of 0.1 m/s.After conducting the test for 5 repetitions, the static angle of repose of the soil was determined to be 44.32°, with a standard deviation of 2.001°, as shown in Figure 5.

Calibration of the DEM Parameters of the Soda Saline-Alkali Soil
Numerous parameters exert influence on the static angle of repose of soil.In this study, we initially conducted the Placke -Burman (PB) test to identify the most significant parameters impacting the static angle of repose of soil.Subsequently, the steepest ascent (SA) test was employed to narrow down the range of these influential parameters.Following this, a quadratic regression orthogonal test was executed to formulate a multiple regression equation.The optimal combination of simulation parameters, tailored to the soda saline-alkali soil of the Songnen Plain, was determined through the solution optimization of this equation.Finally, the obtained parameter combinations were used to simulate the soil rest angle accumulation test and compared with the experimental values.Table 1 shows the high (+1) and low (−1) levels of parameters A~I in the PB test.[32][33][34], which cannot reflect the magnitude of the soil water content.The bonding between particles needed to be characterized by the JKR surface energy, which was subsequently calibrated by the simulation test in this paper, and the Rayleigh time step of the simulation test was 1.5 × 10 −4 s.The simulation screenshot of the soil AoR test is shown in Figure 6 below.

DS Test
The accuracy and effectiveness of the simulation parameters were further verified by the direct shear test.As shown in Figure 7, a ZJ-type automatic quadruple direct shear tester was used to determine the shear strength of soil under four vertical loads (50, 100, 200, and 300 kPa) in soda saline soils in the northeast.Three soil specimens were taken for each vertical pressure, and three tests were conducted.The test data were recorded to obtain the soil shear strength versus vertical load curves under different vertical pressures.

Soda Saline-Alkali Soil Particle Size and pH Value
The test results are shown in Figures 8 and 9.In Figure 8, the particle size distribution under the tillage layer is roughly the same, except for the surface layer of soda salinealkali soil at 0-10 cm, where the particle size distribution is special.The soda saline-alkali soil sample texture is sandy loam.Figure 9 shows the soda saline-alkali soil pH value test results under different tillage layers.It can be seen in the Figure 9 that the soil particles are alkaline under the 0-10 cm tillage layer.Figure 11 shows the density of soda saline-alkali soil under different plow layers.It can be seen in the figure that the density increases with the deepening of the plow layer.

Water Content, Density, and Firmness Test Results of Soda Saline-Alkali Soil
Figure 12 shows the firmness of soda saline-alkali soil under different plow layers.It can be seen in the figure that the difference in firmness is not significant with the deepening of the plow layer.

Calibration Results of the DEM Parameters of the Soda Saline-Alkali Soil
Table 2 shows the PB design and results.The ANOVA (shown in Appendix A, Table A1) showed that the model was significant (p-value < 0.05), and the parameter I (soil surface energy) had the most significant effect on the static angle of repose.After that, the parameters that had a great influence on the static angle of repose were D (soil-soil static friction coefficient), E (soil-soil rolling friction coefficient), A (Poisson's ratio), F (soil-65 Mn steel restitution coefficient), B (shear modulus), and C (soil-soil restitution coefficient).After that, the SA test was carried out on I (soil surface energy), D (soil-soil static friction coefficient), and E (soil-soil rolling friction coefficient) to narrow the range of parameters with significant influence.The remaining parameters with less influence were taken as the intermediate value of the PB test.

No
The results of the SA test are shown in Table 3.At the level of the No. 3 test, when D was 0.6, E was 0.2, and I was 0.55, the static angle of repose was 44.44° ± 0.82, the static angle of repose of the physical test was 44.32° ± 2.001, and the minimum relative error was 0.28%.The No. 3 test level was used as the central point in the subsequent experiments, and No. 2 and No. 3 were used as low and high levels for the experimental design.From the BB test ANOVA (shown in Appendix A, Table A2), we can see that D (soilsoil static friction coefficient), E (soil-soil rolling friction coefficient), and I (soil surface energy) have significant effects on the static angle of repose.The parameters D (soil-soil static friction coefficient) and I (soil surface energy) have more significant effects on the static angle of repose than E (soil-soil rolling friction coefficient).
Additionally, from the relationship curve between the experimental values and the predicted values of the fi ed model (Figure 13), it can be seen that the data points are distributed around the fi ed line with small errors, indicating that the model fits well and can be used for a predictive analysis of the static rest angle.The response surface analysis of the regression model can obtain the influence of the interaction of various factors on the angle of repose, where the regression model of the angle of repose and each significant parameter can be described as Equation (8).Taking the physical repose angle (44.32° ± 2.001) of soda saline-alkali soil in Northeast China as the target value and then solving the regression model, the parameter range can be determined as follows: Using Design-Expert 11.0 software and taking the angle of repose of the physical experiment as the target value, the second-order regression equation was obtained by analyzing the BB experimental data.The equation was optimized to obtain a set of parameters: the soil-soil static friction coefficient was 0.734, the soil-soil rolling friction coefficient was 0.151, the JKR surface energy was 0.56, and the other nonsignificant parameters were the same as the steepest climbing test.
The best parameter combination was used as the EDEM simulation parameter to verify the accuracy and reliability of soil parameter calibration.The static AoR simulation test was conducted many times.The average static AoR of the soil is 44.034°, the standard deviation is 0.911°, the physical static AoR is 44.32°, the standard deviation is 2.001°, and the relative error is 0.645%, as shown in Figure 14; there is no significant difference between the simulated and the physical static AoR of the soil.

DS Test Simulation Results
The DS simulation test is shown in Figure 15, using the combination of parameters obtained from the calibration.The plate was added to apply pressure to the soil sample to simulate different vertical loads in the DS test, and the total shear duration was set to 7 s.When the simulation is 1.5 s, the soil particles have been completely generated.At 1.5 s in the simulation, the plate begins to load in the +Z direction.

Contact Model and Parameter Correction
To enhance the validation's reliability, we employed two distinct validation methods: the angle of repose (AoR) test and the direct shear (DS) test.During the verification process, it became evident that employing the same set of simulation parameters yielded different outcomes for these two tests.After that, we modified the relevant model parameters.Equation ( 11) is the calculation formula of model-related parameters.
In Equation (11), Ab is the contact area between soil particles, m 2 ; Rb is the bond radius between soil particles, m; J is the cross-sectional polar moment of inertia of soil particles, m 4 ; Sn is the normal stiffness of bonded particles, N/m 3 ; Sτ is the tangential stiffness of bonded particles, N/m 3 ; vn and vτ are the normal and tangential components of particle velocity, m/s; ωn and ωτ are the normal and tangential components of particle angular velocity, rad/s; and Δt is the time step, s.
When the external force exceeds a preset value, the bond is destroyed.The normal critical stress, σmax, and the tangential critical stress, τmax, are defined as follows: max The bond parameters in the HMB contact model mainly include the normal stiffness, Sn; tangential stiffness, Sτ; critical normal stress, σmax; critical tangential stress, τmax; and bond radius, Rb.Previous studies have shown that the particle behavior is not sensitive to the change of bond stiffness parameters under this model [6].After multiple pre-tests, the normal stiffness per unit area is taken as 1 × 10 5 N/m 3 , the tangential stiffness per unit area is taken as 1 × 10 5 N/m 3 , and σmax = τmax is 8 × 10 3 Pa.The particle contact radius is 1.1 times the particle radius [12].
Meanwhile, the shear modulus affects the time step variation of the DEM solution [36]; the larger the shear modulus is, the smaller the Rayleigh time step, according to Equation (13).In general, the smaller the time step is, the more time the calculation will take, and the higher the calculation accuracy will be [37].To obtain an accurate combination of simulation parameters, we adjusted the shear modulus and Rayleigh time step to obtain more accurate calculation results for subsequent simulations.In the direct shear simulation test, the Rayleigh time step is 4.76 × 10 −5 s, and the shear modulus is 3 × 10 7 Pa.
In Equation ( 14), Tr is the Rayleigh time step, R is the particle radius, ρ is the particle density, v is Poisson's ratio, and G is the shear modulus.
After modifying the model parameters and conducting the DS simulation again (Figure 17), was is obvious that the bonds between soil particles were constantly broken during the DS test process.Because the bonds between soil particles have a certain strength, the simulated shear strength is improved compared with the single JKR model.The simulation results with a small relative error with the physical test were obtained, as shown in Figure 18.This simulation parameter combination is used to simulate the AoR test again, as shown in Figure 19.After multiple tests, the static angle of repose is 44.2°, and the standard deviation is 1.406°, which is not significantly different from the physical test angle of repose of 44.32° and the standard deviation of 2.001°, which verifies the reliability of the obtained parameter combination.Table 6 shows the calibrated soil simulation parameters.

Discussion
In this study, the contact parameters between soda saline-alkali and soil-subsoiling components in Northeast China were calibrated.Based on the DEM, the JKR model and HMB model were used as contact models, and the accuracy of the calibrated contact parameter combination was verified by the AoR test and DS test.
Through the PB test results (Table A1), it can be found that both the soil intrinsic parameters and the contact parameters will affect the AoR.The contact parameters between the soils (soil-soil rolling friction coefficient, soil-soil static friction coefficient, and JKR surface energy of the soil) have a more significant effect on the AoR, and the JKR surface energy of the soil has a greater effect than the static friction coefficient and the rolling friction coefficient; the results are similar to those of Xing et al. [19].After that, the SA test and BB test were conducted on the parameters with significant influence, and the BB tests were analyzed by ANOVO (Table A2).The results showed a model p-value < 0.001; a coefficient of determination of R 2 = 0.9661; a coefficient of determination adjusted of R 2 = 0.9225, close to 1; and a coefficient of variation of C.V. = 1.02%, thus indicating a significant relationship between the response values and the parameters.Through the ANOVO analysis, it can be seen that the interaction of two significant parameters will also have a significant impact on the angle of repose.The interaction term containing soil JKR surface energy is more significant than other interaction terms.The results of this study are similar to those of Hao, Li, et al. [16,18].We speculate that because the soil moisture content cannot be characterized when the EDEM software (EDEM version: EDEM 2018, DEM Solutions, Edinburgh, UK, 2002) is used to simulate the soil particles, when the JKR model is used as the simulation contact model, the soil moisture content can only be simulated by changing the surface energy of the JKR, which also indicates that the soil water content has a large effect on its AoR and friction coefficient.
When using EDEM to calibrate soil contact parameters, most scholars use the angle of repose as the response value and use a single contact model as the contact model [6,14,17,21,23].However, during the subsoiling process, the soil is sheared by the subsoiling components.Taking this into account, we conducted direct shear tests to verify the calibrated contact parameters.The results show that there is a certain error between the simulation results of the direct shear test and the physical test results when only using the JKR model as the contact parameter obtained by the contact model (Figure 16).This discrepancy can be a ributed to the fact that the AoR test and DS test characterize the soil's flow and shear characteristics, respectively.A single model cannot precisely simulate both aspects of soil behavior.Through repeated experiments, it was found that the addition of bonds between soil particles to simulate the characteristics of high hardness and easy hardening of saline-alkali soil in Northeast China can obtain results closer to physical experiments.Notably, these bonds between soil particles can be disrupted during the DS test.However, even in the event of bond disruption during the shear process, the soil retains cohesion due to the surface energy.In order to faithfully reflect the shear and flow characteristics of the soil particle population, we conducted the DS test by amalgamating the JKR contact model and the HMB contact model, thereby modifying the discrete element simulation parameters.
We used the modified parameter combination to carry out the repose angle simulation test and the direct shear simulation test, and the results are very close to those of the physical test.Therefore, future research could consider combining multiple contact models to simulate different mechanical properties of soil to obtain more accurate simulation results.

Conclusions
This study utilized DEM to model saline soil particles with a multisphere combination approach and validated the model's accuracy using soda saline soil from Northeast China's Songnen Plain.We introduced a modeling method for individual soil particles based on their geometric shape.The effectiveness of this method was confirmed through comparisons between simulation results and experimental data from the angle of repose (AoR) and direct shear tests.
1.The physical properties of the soil, including the texture, density, water content, and particle size distribution, were tested.According to the test analysis, the soil studied in this paper is sandy loam.2. Most soil particles can be approximated as spheres, ellipsoids, and cones.Based on the above three shapes, the geometric model of soil particles was established using the method of multisphere combination.Considering the adhesion between soil particles, the JKR model was used to calculate the force between particles.3. The simulation parameters were calibrated by the soil AoR test.The results of the PB test showed that the static friction coefficient, rolling friction coefficient, and JKR surface energy between saline-alkali soil particles in Northeast China significantly affected the static angle of repose.Then, SA and BB tests were carried out to further determine the optimal parameter combination for soda saline-alkali soil simulation in Northeast China.ANOVA was used to analyze the BB test results, a quadratic regression model was established, and the model fit was good.The results showed that the soil-soil static friction coefficient, rolling friction coefficient, and JKR surface energy and parameters significantly affected the resting angle by a two-two interaction.

Figure 1 .
Figure 1.Administrative map of Daan City.

Figure 3 .
Figure 3. Three different shapes of soil particles.

Figure 4 .
Figure 4. Simulation model of three different shapes of soil particles.

Figure 6 .
Figure 6.Simulation screenshot of the soil AoR test.
(a) ZJ-type automatic quadruple direct shear tester.(b)Soil samples after direct shear.

Figure 8 .
Figure 8. Particle size of soda saline-alkali soil under different tillage layers.

Figure 9 .
Figure 9. pH value of soda saline-alkali soil under different tillage layers.

Figure 10
Figure 10 shows the soda saline-alkali soil moisture content under different tillage layers.It can be seen in the figure that the moisture content increases slowly with the deepening of the tillage layer.Figure11shows the density of soda saline-alkali soil under different plow layers.It can be seen in the figure that the density increases with the deepening of the plow layer.Figure12shows the firmness of soda saline-alkali soil under different plow layers.It can be seen in the figure that the difference in firmness is not significant with the deepening of the plow layer.

Figure 10 .
Figure 10.Moisture content of soda saline-alkali soil under different plow layers.

Figure 11 .
Figure 11.Density of soda saline-alkali soil under different plow layers.

Figure 12 .
Figure 12.Variation of soda saline-alkali soil firmness under different tillage depths.

Figure 13 .
Figure 13.Relationship curve between the actual value and predicted value of soil AOR.

Figure 14 .
Figure 14.Comparison of simulation static AoR test and physical static AoR test results.
At 2.5 s, the down box begins to move along the +X direction at a speed of 0.002 m/s until the end of the simulation.The soil shear strength was obtained by exporting the data through the postprocessing module of EDEM software (EDEM version: EDEM 2018, DEM Solutions, Edinburgh, UK, 2002) [35].(a) t = 2 s.(b) t = 3.5 s.(c) t = 5 s.

Figure 15 .
Figure 15.Screenshot of the DS test simulation.

Figure 16 .
Figure 16.Comparison of simulation results and physical results.

Figure 17 .
Figure 17.Simulation screenshot of DS test after correction.

Figure 18 .
Figure 18.Comparison of simulation results and physical results after correction.

Figure 19 .
Figure 19.Comparison of static AoR between simulation test and physical test after correction.

Table 1 .
High and low levels of the PB test.

Table 2 .
PB test design and results.

Table 3 .
Design and results of the SA test.Finally, we conducted the BB test.The test parameter level is shown in Table4, the test design and results are shown in Table5, and the nonsignificant parameters are the same as the SA test.

Table 4 .
Table of parameter levels of the BB test.

Table 5 .
Design and results of the BB test.

Table 6 .
Soil simulation parameters table.

Table A1 .
Tangential stiffness of bonded particles, N/m 3 ANOVA of PB test.

Table A2 .
ANOVA of BB test.