2.2.2. Contact Model Settings
Edem software 2018 integrates many options to apply to different material types of contact. We selected Hertz-Mindlin with Bonding along with the characteristics that the soil particles are difficult to reconnect after fracture in this paper. At the same time, the most commonly used contact theoretical model Hertz-Mindlin (no slip), is selected for ordinary contact calculation.
- (a)
Hertz-Mindlin (no slip)
In this model, the normal force component model is based on the Hertzian contact theory. The tangential force component model is based on the work of Mindlin-Deresiewicz. Both normal and tangential forces have damping components, where the damping coefficient is related to the coefficient of restitution. The tangential friction force follows the Coulomb friction law, and the rolling friction adopts the research results of Sakaguchi. The relevant calculation principles are described in the literature [
23,
24,
25,
26,
27,
28].
The normal force
Fn is expressed as follows:
where
Fn is the normal force, with the unit of N,
E* is equivalent to Young’s modulus, with the unit of Pa,
R* is the equivalent radius, with the unit of m,
δn is the amount of normal overlap, with the unit of m,
vi,j is Poisson’s ratio,
Ei,j is Particle Young’s modulus, with the unit of Pa,
Ri,j is Particle radius, with the unit of m.
The normal damping force
Fn d is expressed as follows:
where
Fn d is Normal damping force, with the unit of N,
β is viscous damping coefficient,
e is coefficient of restitution,
Sn _is normal stiffness, with the unit of N/m,
m* is equivalent mass, with the unit of g,
mi,j is Particle mass, with the unit of g,
is Normal relative velocity, with the unit of m/s.
The expression for the tangential force
Ft is as follows:
where
Ft is the normal force, N,
St is tangential stiffness, with the unit of N/m,
δt is tangential overlap, with the unit of m,
G* is equivalent shear modulus, with the unit of Pa.
The tangential damping force
Ftd is expressed as follows:
where
Ftd is angential damping force, with the unit of N,
St is tangential stiffness, with the unit of N/m,
is Tangential relative velocity, with the unit of m/s.
The total tangential force is limited by the Coulomb friction force. When the tangential force reaches the Coulomb friction force, that is,
Ft >
μsFn, relative sliding occurs between the particles, where
s is the static friction coefficient. To simplify the contact model, The model does not take into account the sliding friction of the contact model. The rolling friction is expressed in terms of the contact torque of the particle model, and the expression is as follows:
where
τi is torque, with the unit of N·m,
μr is coefficient of rolling friction,
Ri is the distance from the contact point to the centroid, with the unit of m,
ωi—angular velocity of particles at the contact point, with the unit of r/min.
- (b)
Hertz-Mindlin with Bonding
Hertz-Mindlin with Bonding contact model bonds particles together. Since this model includes Hertz-Mindlin (no slip), the default Hertz-Mindlin (no slip) contact model needs to be deleted from the software Creator Tree list when generating the bond.
Hertz-Mindlin with Bonding contact model uses finite-sized “bonding” bonds to bond particles. The bond can withstand the normal and tangential motion between the contacting particles to generate normal and tangential shear stress until the bond breaks at the critical maximum stress. The particles interact as rigid spheres, and no new bonds are generated, which is in line with the soil stress mechanical behavior characteristics of remaining loose after crushing [
29,
30,
31].
In this model, the first generated particles interact with each other according to the Hertz-Mindlin contact model, and the connection starts after the set bond generation time point. After the particles are connected, the force and torque between the particles are both set to 0, according to the following Stepwise adjustment at each time step:
where
Fn,t is normal and tangential force, with the unit of N,
vn,t is normal and tangential velocity, with the unit of m/s,
Sn,t is normal and tangential stiffness, with the unit of N/m,
A is the cross-sectional area of bonding bond, with the unit of m
2,
δt is time step, with the unit of s,
Mn,t is normal and tangential moments, with the unit of N·m,
ωn,t is normal and tangential angular velocities, with the unit of r/min,
J is the moment of inertia, with the unit of m
4,
RB is the bond-bond radius.
When the normal and tangential stresses exceed preset critical values, the bonded bonds will break, the particles will interact as rigid spheres, and no new bonds will be formed. Spawn bonds in this model can create cohesive forces when particles are generated, even if the particles are not in contact, so the contact radius should be set larger than the radius of the particle body.
Table 3 lists the parameters adopted in this study. In this article, we use Edem’s dynamic generation method to generate soil particle beds layer by layer, and the model is shown in
Figure 5.
The particle contact model can quantify the interaction between discrete elements within a given contact range. Adding a contact model to the particle model is a key step in the test. Appropriate contact model to accelerate simulation test.
2.2.3. Initial Conditions for Simulation
In this simulation test, in order to ensure that the tillage depth was equal to 200 mm, the rotation center was set to the corresponding position into the soil according to the radius of rotation of each blade, as shown in
Figure 6.
According to the research objectives of this article, combined with the existing re-search experience and actual working conditions, the main indicators were determined as the power (P), the rate of crushed soil (RC), and the rate of straw burial (RB).
Figure 1b shows that the value of B is fixed when R and β are certain, so the factor B is removed. The initial conditions of the experiment are shown in
Table 4, and the factor level design is shown in
Table 5.
In this study, for convenience of calculation, the ratio of the number of bond-bond breaks per unit volume to the number of bond bonds at the beginning is used as the index of soil breakage. The sampling position is set 500 mm from the center of rotation after rotary tiller operation, and 1000 mm is set by grid-grid function after EDEM treatment.1000 mm × 500 mm × 250 mm grid solver calculates the bonding bond between particles as a problem to be solved.
The value of power (P) is selected from the data within 2 s of the operation stabilization stage. Each blade is tested three times, and the average value is obtained, which is directly exported after EDEM software post-processing.
The crushed soil rate (RC) is measured by the ratio of bond breakage between particles in the tillage area. In this study, for the convenience of calculation, the ratio of the number of soil bond breakage per unit volume to the number of bonded bonds at the initial time was used as the indicator of crushed soil rate. The sampling position was set at 500 mm from the center of rotation after the rotary tillage knife operation, and the grid solver of 1000 mm × 500 mm × 250 mm was set using the grid function of EDEM post-processing, and the particle-to-particle bond was calculated as the problem to be solved.
Straw burial rate (RB) is the percentage of straw residual mass on the ground post-tillage to the pre-tillage. In the post-processing of the EDEM software, L0 is the straw layer. The grid solver was divided into four layers (L1, L2, L3, L4) in the vertical direction, and the percentage of the total straw mass from L1 to L4 to the initial L0 mass was calculated as the straw burial rate.