Hi
Thanks for your reply. My confusion was because after using the subroutine, the stresses are computed correctly at the integration points but the strains are not. I thought that in FEM the order of solution is (i) nodal disp> strains> stresses. So, I was wondering how the stresses are correct...