[lammps-users] crack growth simulation

Dear All,

I am using this script to simulate crack propagation in FCC Ni. The code is running and I am able to get the snapshots from XMOVIE.

I expect the crack to open between white and green region but the crack doesnot seem to open between rather the atoms burst out. I think the velocity command or the fix/setforce command I have specified is not right.
Can anyone suggest where am I going wrong. ?


stage 2.png

Stage 1.png

Could be a lot of things. I suggest you run it slowly with a lot
of output and see if it is doing what you expect. By "slowly"
I mean pull it apart with a slow deformation rate.