pk2vol.for

Go to the documentation of this file.

      SUBROUTINE pk2vol(PKVOL,PV,C,NDI)
C>    VOLUMETRIC PK2 STRESS
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1
      DOUBLE PRECISION PKVOL(ndi,ndi),C(ndi,ndi),CINV(ndi,ndi)
      DOUBLE PRECISION PV
C
      CALL matinv3d(c,cinv,ndi)
C
      DO i1=1,ndi
        DO j1=1,ndi
          pkvol(i1,j1)=pv*cinv(i1,j1)
        END DO
      END DO
C
      RETURN
      END SUBROUTINE pk2vol