Hello all
I am carrying out nanoindentation simulation. I want to calculate forces on the indenter. For this I used variable fcm command, however i got zero forces
on indenter during both loading and unloading.
variable l_rate equal -0.5 # loading rate
variable ul_rate equal 0.1 # unloading rate
velocity indenter set 0.0 ${l_rate} 0.0 sum no units box
variable fx equal fcm(indenter,x)
variable fy equal fcm(indenter,y)
variable fz equal fcm(indenter,z)
thermo_style custom step elapsed elaplong temp pe ke etotal press time v_fx v_fy v_fz
How can I get the exact forces?
A-indentation.txt (3.11 KB)
log.lammps (417 KB)
If you dump the forces on the indenter atoms at the
end of the timestep, what do you get? If should
be the same as your fx,fy,fz variables if you sum
them. If that is zero, then you probably did
something else before the end of the timestep
to zero the forces?
Steve
In dump files, atoms have some values of forces.
Is the scalar value of forces zero due to fix setforce command set to zero in all three direction
In dump files, atoms have some values of forces.
Is the scalar value of forces zero due to fix setforce command set to zero
in all three direction
please see the documentation of fix setforce for an explanation.
you have to obtain the forces *before* setting them to zero from the fix.
axel.