I am currently working on indentation using “fix indent” command. I am also able to visualize the indentation using VMD.
But, I am getting the value of force exerted on indenter atoms as zero. I have used the following commands in input script :
fix 6 all indent 1000 sphere 75 75 240 20.0 vel 0.0 0.0 -2.0 units box
fix_modify 6 energy yes
thermo_style custom step temp pe ke f_6 f_6 f_6 f_6
But in thermodynamic output, I am getting only the value of scalar energy but vector forces are zero at every timestep. So, please suggest how do I get the values of forces or if something wrong with the command in input script.
Any suggestion is welcome in this.
If I add this line to lammps/examples/in.indent
thermo_style custom step temp press f_4 f_4 f_4
I get both the energy and force output.
Are you using the most current version of LAMMPS?