fslip.for

Go to the documentation of this file.

      SUBROUTINE fslip(F,FBAR,DET,NDI)
C>     DISTORTION GRADIENT
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1
      DOUBLE PRECISION F(ndi,ndi),FBAR(ndi,ndi)
      DOUBLE PRECISION DET,SCALE1
C     
C     JACOBIAN DETERMINANT
      det = f(1,1) * f(2,2) * f(3,3)
     1    - f(1,2) * f(2,1) * f(3,3)
C
      IF (ndi .EQ. 3) THEN
          det = det + f(1,2) * f(2,3) * f(3,1)
     1              + f(1,3) * f(3,2) * f(2,1)
     2              - f(1,3) * f(3,1) * f(2,2)
     3              - f(2,3) * f(3,2) * f(1,1)
      END IF 
C
      scale1=det**(-one /three)
C      
      DO i1=1,ndi
        DO j1=1,ndi
          fbar(i1,j1)=scale1*f(i1,j1)
        END DO
      END DO
C
      RETURN      
      END SUBROUTINE fslip