[lammps-users] creep simulation

Dear Developers and Users of LAMMPS
I am trying to do MD simulation for the estimation of creep properties of metals. But after many attempts, I am unable to find satisfactory results as creep strain vs time at constant applied stress as reported similarly published praper . I have used fix deform, addforce and setforce command during LAMMPS simulations but not getting actual results and also simulation output showing some arbitrary energy value. Following is my script. Kindly suggest me any modification. I will be too much grateful.

13631.full.pdf (1.55 MB)

nanomaterials-10-01693.pdf (4.48 MB)