push4.for

Go to the documentation of this file.

      SUBROUTINE push4(SPATIAL,MAT,F,DET,NDI)
C>        PIOLA TRANSFORMATION
C>      INPUT:
C>       MAT - MATERIAL ELASTICITY TENSOR
C>       F - DEFORMATION GRADIENT
C>       DET - DEFORMATION DETERMINANT
C>      OUTPUT:
C>       SPATIAL - SPATIAL ELASTICITY TENSOR
       IMPLICIT NONE
       include 'PARAM_UMAT.INC'
C
       INTEGER I1,J1,K1,L1,II1,JJ1,KK1,LL1,NDI
C
       DOUBLE PRECISION MAT(ndi,ndi,ndi,ndi),F(ndi,ndi)
       DOUBLE PRECISION SPATIAL(ndi,ndi,ndi,ndi)
       DOUBLE PRECISION AUX,DET
C
       DO i1=1,ndi
        DO j1=1,ndi
         DO k1=1,ndi
          DO l1=1,ndi
           aux=zero
           DO ii1=1,ndi
            DO jj1=1,ndi
             DO kk1=1,ndi
              DO ll1=1,ndi
              aux=aux+(det**(-one))*
     +        f(i1,ii1)*f(j1,jj1)*
     +        f(k1,kk1)*f(l1,ll1)*mat(ii1,jj1,kk1,ll1)
              END DO
             END DO
            END DO
           END DO
           spatial(i1,j1,k1,l1)=aux
          END DO
         END DO
        END DO
       END DO
C
       RETURN
      END SUBROUTINE push4