Hi all,
I had a look at the gran/hook/history/simple pair_style today. I thought that this pair_style should make you able to directly set the values for kn, kt, gamman and gammat depending on the contact-partners. So I was wondering why the values for the damping coefficients (gamman and gammat) are multiplied by meff. The same modification is done for wall-contacts.
I think that should be changed otherwise you have to take the multiplications into account when you set the values in the input-script.
Regards,
Yansan
ckloss | Tue, 08/30/2011 - 09:10
Hi Yansan, this is the
Hi Yansan,
this is the implementation as it is in LAMMPS - see http://lammps.sandia.gov/doc/pair_gran.html. But you're right and option to directly set it would be good - I will add a switch to turn this on/off in the next release
Christoph