Hi all,
I am new to LAMMPS. Following the method detailed in a journal article. I am now implementing a jamming simulation of granular materials, where the particle normal interaction is Hooke with damp component (mass_velocity).
Specifically, it is mentioned that “An extra switching function is employed to set the magnitude of normal force to zero” if the damp component is larger than the Hooken component (which will result in attractive force). However, I do not know how I can specify the “switching function” in my script.
I think this modification matters since I failed to reproduce their results (other settings are basically the same), and I am looking forward to your feedback.
Here is my current settings in script:
pair_style granular
pair_coeff * * hooke 1 0.5 damping mass_velocity tangential linear_history 1 1 0
comm_modify vel yes
Best,
Andrew