Linear Spring Model

Also set the yield strength in units of Pa for each particle or Geometry that has an interaction. EDEM offers a reasonable default value for this parameter estimated from the material's Shear Modulus. You can, however, overwrite this default value. A linear spring with stiffness k is in parallel with a dashpot with coefficient, c. The magnitude of the normal force between two particles, F_N, is defined as:

Where k is the linear spring stiffness, c is the dashpot coefficient, δ is the overlap, and δ_dot is the overlap velocity. You can apply a similar force to the tangential direction.

The spring stiffness and the dashpot coefficient are the parameters in this model and it is common practice to estimate the spring stiffness and calculate the dashpot coefficient based on this stiffness. The simulation Time Step is then estimated based on the spring stiffness.

The spring constant and dashpot coefficient can be calculated based on a combination of material properties and kinematic constraints. One common method is obtained by equating the maximum strain energy in a purely Hertzian contact (E_hertzian) with the maximum strain energy of the existing contact (E_max) as follows:

Where the equivalent mass m^*, the equivalent radius R^*, and the equivalent Young’s modulus E* were defined earlier. V is the typical impact velocity.

For two identical spherical particles with masses of 7.63e-03 kg, radius of 9mm and Young’s modulus of 2.6e+08 Pa, colliding at a velocity of 3 m/s, k ≈ 2.0e+05 N/m.

The impact velocity in an EDEM simulation can usually be taken as a characteristic velocity in the simulation. You can base this velocity as the maximum velocity in the simulation, for example, for a blending operation with the blender operating at Ω rad/s, the characteristic velocity is equal to r Ω m/s, where r is the radius of the blender. The dashpot coefficient is related to the Coefficient of Restitution as follows:

Where e is the Coefficient of Restitution.

The tangential stiffness is usually estimated as a ratio to the normal spring stiffness (Cundall and Strack 1979). v has the tangential stiffness equal to the normal stiffness. The dashpot coefficient is calculated using the tangential stiffness in the previous equation. The tangential force is defined as:

Where k_t and c_t are the tangential spring and dashpot coefficient, μ is the Coefficient of Friction.

The simulation Time Step is usually a small percentage of the contact duration of the particles. The contact duration for the Linear Spring model is obtained using the normal stiffness is defined as:

Where β = π /ln(e) and e is the Coefficient of Restitution.

For e = 0.5, the contact time is 0.00043 sec. The simulation Time Step must be less than this value for better integration of the particle states. It is recommended to have a value of about 5-10 % of this contact time for accurate results. The details of the soft particle contact model are relatively unimportant due to the fact that a lumped parameter approach which neglects the details of the contact force (Coefficient of Restitution) is sufficient to describe the collision dynamics.