0 rad 0 0 1 mm 180/pi()*(heaviside((heaviside(z+1E-10+u.z)*atan2(z+1E-10+u.z,y+u.y)+heaviside(-1*(z+1E-10+u.z))*(atan2(z+1E-10+u.z,y+u.y)+2*pi()))-(heaviside(z+1E-10)*atan2(z+1E-10,y)+heaviside(-1*(z+1E-10))*(atan2(z+1E-10,y)+2*pi())))*((heaviside(z+1E-10+u.z)*atan2(z+1E-10+u.z,y+u.y)+heaviside(-1*(z+1E-10+u.z))*(atan2(z+1E-10+u.z,y+u.y)+2*pi()))-(heaviside(z+1E-10)*atan2(z+1E-10,y)+heaviside(-1*(z+1E-10))*(atan2(z+1E-10,y)+2*pi())))+heaviside(-1*((heaviside(z+1E-10+u.z)*atan2(z+1E-10+u.z,y+u.y)+heaviside(-1*(z+1E-10+u.z))*(atan2(z+1E-10+u.z,y+u.y)+2*pi()))-(heaviside(z+1E-10)*atan2(z+1E-10,y)+heaviside(-1*(z+1E-10))*(atan2(z+1E-10,y)+2*pi()))))*((heaviside(z+1E-10+u.z)*atan2(z+1E-10+u.z,y+u.y)+heaviside(-1*(z+1E-10+u.z))*(atan2(z+1E-10+u.z,y+u.y)+2*pi()))-(heaviside(z+1E-10)*atan2(z+1E-10,y)+heaviside(-1*(z+1E-10))*(atan2(z+1E-10,y)+2*pi()))+2*pi())) MPa vm*(s1+s2+s3)/abs(s1+s2+s3)