gamman and gammat specified in gran model hooke/stiffness

Submitted by zhou_yu on Wed, 12/23/2015 - 15:35

when using gran model hooke/stiffness, I noticed that the gamman and gammat specified are weird for the case of absolute_damping = 'off'.
In another dem code PFC, the parameter specified for viscous damping is damping ratio, however liggghts gives the expression:
gamma(n)=m(eff)*gamma(n,specified)
I wonder if this gamma(n,specified) has a specific meaning that I ignore.

In PFC ,the expression is:
gamma(n)=2sqrt[m(eff)*kn]*gamma(n,specified)
and the gamma(n,specified) is called damping ratio in the normal direction.

thank you.

ckloss's picture

ckloss | Tue, 12/29/2015 - 15:39

Hi zhou_yu,

LIGGGHTS is following the implementation as it is in LAMMPS. However, if you prefer the PFC version, it is pretty easy to implement a modified versioN!

Best wishes
Christoph