fibdir.for

Go to the documentation of this file.

      SUBROUTINE fibdir(FIB,ST0,ST,NE,NOEL,NDI,VORIF,VD,DISTGR,DFGRD1)
C
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI, NE, NOEL,INOEL,I,J,I1,J1
      DOUBLE PRECISION SUM1, DFGRD1(3,3), DNORM
      DOUBLE PRECISION VORIF(3),ST(3,3),VD(3),ST0(3,3),DISTGR(3,3)
      DOUBLE PRECISION FIB(ne,4)
C
        inoel=0
        i=0
        DO i=1,ne
C               ELEMENT IDENTIFICATION
            IF(noel.EQ.int(fib(i,1))) THEN
                inoel=i
            ENDIF
        ENDDO
C
C     FIB - FIBER ORIENTATION
             dnorm=dsqrt(fib(inoel,2)*fib(inoel,2)+
     1                   fib(inoel,3)*fib(inoel,3)+
     2                   fib(inoel,4)*fib(inoel,4))
C
C       UNDERFORMED FIBER ORIENTATION TENSOR
C
        DO i1=1,ndi
        j1=i1+1
C       FIBER ORIENTATION NORMALIZED VECTOR - FAMILY 1
        vorif(i1)=fib(inoel,j1)/dnorm
        END DO
C

      DO i=1,ndi
         sum1=zero
         DO j=1,ndi
          sum1=sum1+dfgrd1(i,j)*vorif(j)
         ENDDO
C     FIBER DIRECTIONS IN THE DEFORMED CONFIGURATION
C               -FAMILY 1
         vd(i)=sum1
      ENDDO
      dnorm=dsqrt(vd(1)*vd(1)+
     1             vd(2)*vd(2)+
     2             vd(3)*vd(3))
C           COSINE OF THE ANGLE BETWEEN FIBERS
C
C
C--------------------------------------------------------------------------
      DO i=1,ndi
       DO j=1,ndi
C       STRUCTURAL TENSOR - FAMILY 1
       st0(i,j)=vorif(i)*vorif(j)
       END DO
      END DO
C
C       STRUCTURE TENSOR IN THE DEFORMED CONFIGURATION - FAMILY 1
      st=matmul(st0,transpose(distgr))
      st=matmul(distgr,st)
C
C
      RETURN
      END SUBROUTINE fibdir