Discussion: Implementation of tangential spring in pair_style gran/hooke
Hi all,
I think there are some problems with the implementation of the tangential spring in the pair_style gran/hooke, especially for soft particles where the overlap between the particles gets high.
So far the relative rotational velocities between 2 particles are calculated by using the radius of the particles:
wr1 = (radi*omega[i][0] + radj*omega[j][0]) * rinv;
The force application point for the tangential forces also lies on the radius of the particles:
torque[i][0] -= radi*tor1;