setvol.for

Go to the documentation of this file.

      SUBROUTINE setvol(CVOL,PV,PPV,UNIT2,UNIT4S,NDI)
C>    VOLUMETRIC SPATIAL ELASTICITY TENSOR
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1,K1,L1
      DOUBLE PRECISION UNIT2(ndi,ndi),UNIT4S(ndi,ndi,ndi,ndi),
     1                 cvol(ndi,ndi,ndi,ndi)
      DOUBLE PRECISION PV,PPV
C
      DO i1 = 1, ndi
        DO j1 = 1, ndi
         DO k1 = 1, ndi
           DO l1 = 1, ndi
             cvol(i1,j1,k1,l1)=
     1                 ppv*unit2(i1,j1)*unit2(k1,l1)
     2                 -two*pv*unit4s(i1,j1,k1,l1)
           END DO
         END DO
        END DO
      END DO
C      
      RETURN
      END SUBROUTINE setvol