index.for

Go to the documentation of this file.

      SUBROUTINE indexx(STRESS,DDSDDE,SIG,TNG,NTENS,NDI)
C>    INDEXATION: FULL SIMMETRY  IN STRESSES AND ELASTICITY TENSORS
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER II1(6),II2(6),NTENS,NDI,I1,J1
      DOUBLE PRECISION STRESS(ntens),DDSDDE(ntens,ntens)
      DOUBLE PRECISION SIG(ndi,ndi),TNG(ndi,ndi,ndi,ndi)
      DOUBLE PRECISION PP1,PP2
C
      ii1(1)=1
      ii1(2)=2
      ii1(3)=3
      ii1(4)=1
      ii1(5)=1
      ii1(6)=2
C
      ii2(1)=1
      ii2(2)=2
      ii2(3)=3
      ii2(4)=2
      ii2(5)=3
      ii2(6)=3
C
      DO i1=1,ntens
C       STRESS VECTOR
         stress(i1)=sig(ii1(i1),ii2(i1))
         DO j1=1,ntens
C       DDSDDE - FULLY SIMMETRY IMPOSED
            pp1=tng(ii1(i1),ii2(i1),ii1(j1),ii2(j1))
            pp2=tng(ii1(i1),ii2(i1),ii2(j1),ii1(j1))
            ddsdde(i1,j1)=(one/two)*(pp1+pp2)
         END DO
      END DO
C
      RETURN
C
      END SUBROUTINE indexx