pk2iso.for

Go to the documentation of this file.

      SUBROUTINE pk2iso(PKISO,PKFIC,PL,DET,NDI)
C>    ISOCHORIC PK2 STRESS TENSOR
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1
      DOUBLE PRECISION PKISO(ndi,ndi),
     1                 pl(ndi,ndi,ndi,ndi),pkfic(ndi,ndi)
      DOUBLE PRECISION DET,SCALE2      
C
      CALL contraction42(pkiso,pl,pkfic,ndi)
C
      scale2=det**(-two/three)
      DO i1=1,ndi
        DO j1=1,ndi
          pkiso(i1,j1)=scale2*pkiso(i1,j1)
        END DO
      END DO
C
      RETURN
      END SUBROUTINE pk2iso