I don't know wheter there are other problems in your calculation but your convergence threshold of 1.d-1 looks to me extremely poor... forces and stress are therefore probably completely wrong. stefano