setjr.for

Go to the documentation of this file.

      SUBROUTINE setjr(CJR,SIGMA,UNIT2,NDI)
C>    JAUMAN RATE CONTRIBUTION FOR THE SPATIAL ELASTICITY TENSOR
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1,K1,L1
      DOUBLE PRECISION UNIT2(ndi,ndi),
     1                 cjr(ndi,ndi,ndi,ndi),sigma(ndi,ndi)
C
      DO i1 = 1, ndi
        DO j1 = 1, ndi
         DO k1 = 1, ndi
           DO l1 = 1, ndi
              cjr(i1,j1,k1,l1)=
     1             (one/two)*(unit2(i1,k1)*sigma(j1,l1)
     2             +sigma(i1,k1)*unit2(j1,l1)+unit2(i1,l1)*sigma(j1,k1)
     3             +sigma(i1,l1)*unit2(j1,k1))
           END DO
         END DO
        END DO
      END DO
C
      RETURN
      END SUBROUTINE setjr