Hello everyone,
I am running a simulation where instead of walls I use particles. These particles are inserted at the desired locations and their torque and force are set to 0 with the command freeze. By doing this, these particles (let's call them type 1) are fixed in the space and they work as a walls for the rest of the particles (type 2) that are not frozen. However, I would like to know the contact forces that the particles of type 2 exert on every wall particle (type 1). If I would have a wall I would use the command store_force but I cannot apply it on the wall particles. If I try to save the particle data using the compute stress/atom, the forces on every particle is 0 because it's already truncated by the command freeze. Then, is there a way of having access to the forces on the particles that are frozen? or, is there another command that fixes the particle positions and allow us to have the forces on every particle?
Many thanks for your help.
Ivan