Dear LAMMPS community,
I am trying to apply a constant stress by using ‘fix aveforce’ command.
I define two regions ‘upper’ and ‘lower’, where I apply aveforce. The shear stress is derived as stress=force_total/Area.
The following inscript seems to lead a shear stress of about 160 MPa.
However, the shear stress is actually too large, so that it catastrophically destroys my model.
Any idea about applying a stress with reasonable magnitude?