push2.for

Go to the documentation of this file.

      SUBROUTINE push2(SIG,PK,F,DET,NDI)
C>        PIOLA TRANSFORMATION
C>      INPUT:
C>       PK - 2ND PIOLA KIRCHOOF STRESS TENSOR
C>       F - DEFORMATION GRADIENT
C>       DET - DEFORMATION DETERMINANT
C>      OUTPUT:
C>       SIG - CAUCHY STRESS TENSOR
       IMPLICIT NONE
       include 'PARAM_UMAT.INC'
C
       INTEGER I1,J1,II1,JJ1,NDI
       DOUBLE PRECISION PK(ndi,ndi),F(ndi,ndi)
       DOUBLE PRECISION SIG(ndi,ndi)
       DOUBLE PRECISION AUX,DET
C
       DO i1=1,ndi
        DO j1=1,ndi
          aux=zero
         DO ii1=1,ndi
          DO jj1=1,ndi
            aux=aux+(det**(-one))*f(i1,ii1)*f(j1,jj1)*pk(ii1,jj1)
         END DO
        END DO
           sig(i1,j1)=aux
        END DO
       END DO
C
       RETURN
      END SUBROUTINE push2