Hi LAMMPS users,
I’m writing to get some help from you about the problem of domain size-dependent energy minimization.
I’m trying to relax a subdomain of a deformed structure using the energy minimizer. I want to treat some parts of my domain as a rigid body and only allow relaxation of a small part of my domain. To achieve this, I used fix setforce 0 0 null to fix the x and y position of my ‘rigid body’ and allow it to relax along the z direction. The mobile part is allowed to relax along all directions. The part of setting is following: