metvol.for

Go to the documentation of this file.

      SUBROUTINE metvol(CVOL,C,PV,PPV,DET,NDI)
C>    VOLUMETRIC MATERIAL ELASTICITY TENSOR
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1,K1,L1
      DOUBLE PRECISION C(ndi,ndi),CINV(ndi,ndi),
     1                 cvol(ndi,ndi,ndi,ndi)
      DOUBLE PRECISION PV,PPV,DET
C
      CALL matinv3d(c,cinv,ndi)
C
      DO i1 = 1, ndi
        DO j1 = 1, ndi
         DO k1 = 1, ndi
           DO l1 = 1, ndi
             cvol(i1,j1,k1,l1)=
     1                 det*ppv*cinv(i1,j1)*cinv(k1,l1)
     2           -det*pv*(cinv(i1,k1)*cinv(j1,l1)
     3                      +cinv(i1,l1)*cinv(j1,k1))
           END DO
         END DO
        END DO
      END DO
C
      RETURN
      END SUBROUTINE metvol