Dear LAMMPS Users;
I am currently using fix "indent command" in my simulations. The
simulations work fine and I think the force exerted on system should
be evaluated as: F = -k(r-R), where k is the force constant of
indenter, r is the final position of indenter, and R is the initial
position. However, when I see the lammps documentation webpage
(http://lammps.sandia.gov/doc/fix_indent.html : UNDER DESCRIPTION), it
says lammps calculates the force as: F = -k(r-R)^2. I am wondering why
there is power 2 for r-R term. Is it a typo on the webpage or lammps
calculates force in s different way or I am making some mistake.
Please let me know the exact procedure of evaluating the force in
lammps for fix indenter command.