Hertz-Mindlin (no slip)

This contact model is the default model used in EDEM due to its accurate and efficient force calculation. In this model the normal force component is based on Hertzian contact theory (Hertz 1882). The tangential force model is based on Mindlin-Deresiewicz work (Mindlin 1949) (Mindlin and Deresiewicz 1953). Both normal and tangential forces have damping components where the damping coefficient is related to the coefficient of restitution as described in (Tsuji, Tanaka and Ishida 1992). The tangential friction force follows the Coulomb law of friction model as in, for example, (Cundall and Strack 1979). The rolling friction is implemented as the contact independent directional constant torque model, see, for example, (Sakaguchi, Ozaki and Igarashi 1993).

 

 

In particular, the normal force, F_n, is a function of normal overlap δ_n and is given by:

(1)

 

 

 

Where the equivalent Young’s Modulus E*, the equivalent radius R* are defined as:

		(2)
		(3)

 

 

 

 

with Eⁱ, vⁱ, Rⁱ, and E^j, v^j, R^j, being the Young’s Modulus, Poisson ratio and Radius of each sphere in contact. Additionally there is a damping force, F_n^d, given by:

(4)

 

 

Where is the equivalent mass, v_n^rel is the normal component of the relative velocity and β and S_n(the normal stiffness) are given by:

		(5)
		(6)

 

 

 

 

With e the coefficient of restitution. The tangential force, F_t, depends on the tangential overlap δ_t and the tangential stiffness S_t.

(7)

 

with

 

(8)

 

 

Here G* is the equivalent shear modulus. Additionally, tangential damping is given by:

(9)

 

 

 

where v_trel is the relative tangential velocity. The tangential force is limited by Coulomb friction μ_sF_n where μ_s is the coefficient of static friction.