Hi
I have a doubt in implementation of Hertz model in LIGGGHTS. The following code is implemented in LIGGGHTS to calculate Normal hertz stiffness.
---------------------------------------------------------------------
double kn=4./3.*Yeff[itype][jtype]*sqrtval;
---------------------------------------------------------------------
I have gone through some of Literature and It look like the factor (4/3) should be replaced by 2.
---------------------------------------------------------------------
double kn=2.*Yeff[itype][jtype]*sqrtval;
---------------------------------------------------------------------
Please let me know if I am right ?
richti83 | Tue, 12/01/2020 - 21:24
Years ago I thought about the
Years ago I thought about the same question, here is the solution:
https://pasteboard.co/JCXjteg.png
(answer: depends on definition of effective shear- and youngs modulus (Gstern=G* :-)) )
Rahul Kumar | Wed, 12/02/2020 - 13:52
Thanks
Thanks for making this clear