I have been trying to simulate the cyclic triaxial test using the SANISAND04 subroutine. When I plotting the stress path (p-q plot), I am getting only the positive value of q. Maybe because the definition of q i.e; sqrt(J2) always gives a positive value. How to get the negative values of q?. Am I missing something in modeling in Abaqus?. Please help me to figure things out. I have attached the exact figure in which I am getting the values below the p axis.