Dear LAMMPS users,
For some reasons, I am trying to move a wall in x direction. I have prepared a simple input file as follows. Harmonic bonds are defined between each wall atom and it’s three neighboring atoms.
When I use setforce command to move the wall, I don’t receive any error and the wall moves without any deflection or disorganization. But if I use SMD or velocity commands, I receive shake atoms error. it seems that some restrictions do not allow some wall atoms to move in x direction and it leads to deterioration of the wall lattice.
Could you please help me to solve the problem?
Best,
Omar
INPUT:
units real
neighbor 2.0 bin
neigh_modify delay 0 every 1 check yes
boundary p p f
#************************* styles **************************
atom_style full
bond_style harmonic
angle_style harmonic
pair_style lj/cut/coul/long 11.0