I am trying to apply equal but opposite shear stress to the fixed upper and lower region of the simulation cell. The applied shear will be in the y-direction. Please advice me on the set up and variables definition requirements and the lammps commands to do this job.
I will greatly appreciate any advice or the similar works done before