proj_lagrangian.for

Go to the documentation of this file.

      SUBROUTINE projlag(C,AA,PL,NDI)
C>    LAGRANGIAN PROJECTION TENSOR
C      INPUTS:
C          IDENTITY TENSORS - A, AA
C          ISOCHORIC LEFT CAUCHY GREEN TENSOR - C
C          INVERSE OF C - CINV
C      OUTPUTS:
C          4TH ORDER SYMMETRIC LAGRANGIAN PROJECTION TENSOR - PL
C
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER I,J,K,L,NDI
C
      DOUBLE PRECISION CINV(ndi,ndi),AA(ndi,ndi,ndi,ndi),
     1                 pl(ndi,ndi,ndi,ndi),c(ndi,ndi)
C
      CALL matinv3d(c,cinv,ndi)
C
      DO i=1,ndi
         DO j=1,ndi
          DO k=1,ndi
             DO l=1,ndi
              pl(i,j,k,l)=aa(i,j,k,l)-(one/three)*(cinv(i,j)*c(k,l))
           END DO
          END DO
        END DO
      END DO
C
      RETURN
      END SUBROUTINE projlag