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, Fn, is a function of normal overlap δn and is given by:





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






with Ei, vi, Ri, and Ej, vj, Rj, being the Young’s Modulus, Poisson ratio and Radius of each sphere in contact. Additionally there is a damping force, Fnd, given by:




Where is the equivalent mass, vnrel is the normal component of the relative velocity and β and Sn(the normal stiffness) are given by:






With e the coefficient of restitution. The tangential force, Ft, depends on the tangential overlap δt and the tangential stiffness St.








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





where vtrel is the relative tangential velocity. The tangential force is limited by Coulomb friction μsFn where μs is the coefficient of static friction.





(c) 2022 Altair Engineering Inc. All Rights Reserved.

Intellectual Property Rights Notice | Technical Support