[lammps-users] Regarding force exerted by indenter in "fix indent"

Hello Sir/friends,

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[1] f_6[2] f_6[3] 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.

Sarang Muley.


If I add this line to lammps/examples/in.indent

thermo_style custom step temp press f_4 f_4[1] f_4[2]

I get both the energy and force output.

Are you using the most current version of LAMMPS?