The Spinning Friction contact model is used to account for the friction that would
occur if a particle face is rotating against another particle or Geometry.
Friction involves calculating the contact area and then using this value to obtain the
torque on the particles as follows:
Where Fn is the normal force, μ is the coefficient of friction, and
Rdisk is the effective disk radius given by approximating the contact as
a disk with area equal to the normal area Acontact of the overlap
region:
The torque acts in the opposite direction to the normal part of the relative angular
velocity.
A Limit is applied to this torque to avoid oscillating behavior when the angular velocity
is small. This is done by finding the torque which would completely eliminate the
angular velocity in one Time Step and is described as:
The ωrel,n is the normal component of the relative angular velocity, Δt is the
Time Step, and Ix,y,z is the minimum component of the moment of inertia of
either of the two contacting objects excluding geometries. The extra Csafety
factor is to try to avoid possible instability in cases with multiple contacts and at
the moment has a value of 0.125.
The damping to damp rocking is defined as:
Where γn is the linear damping coefficient in the normal direction from the
base model, and ωrel,t is the tangential part of the relative angular
velocity.