stretch.for

Go to the documentation of this file.

      SUBROUTINE stretch(C,B,U,V,NDI)
C>    STRETCH TENSORS
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI
      DOUBLE PRECISION C(ndi,ndi),B(ndi,ndi),U(ndi,ndi),V(ndi,ndi)
      DOUBLE PRECISION EIGVAL(ndi),OMEGA(ndi),EIGVEC(ndi,ndi)
C
      CALL spectral(c,omega,eigvec)
C
      eigval(1) = dsqrt(omega(1))
      eigval(2) = dsqrt(omega(2))
      eigval(3) = dsqrt(omega(3))
C
      u(1,1) = eigval(1)
      u(2,2) = eigval(2)
      u(3,3) = eigval(3)
C
      u = matmul(matmul(eigvec,u),transpose(eigvec))
C
      CALL spectral(b,omega,eigvec)
C
      eigval(1) = dsqrt(omega(1))
      eigval(2) = dsqrt(omega(2))
      eigval(3) = dsqrt(omega(3))      
C      write(*,*) eigvec(1,1),eigvec(2,1),eigvec(3,1)
C
      v(1,1) = eigval(1)
      v(2,2) = eigval(2)
      v(3,3) = eigval(3)
C
      v = matmul(matmul(eigvec,v),transpose(eigvec))
      RETURN
      END SUBROUTINE stretch