cmatisomatfic.for

Go to the documentation of this file.

      SUBROUTINE cmatisomatfic(CMISOMATFIC,CBAR,CBARI1,CBARI2,
     1                          diso,unit2,unit4,det,ndi)
C>    ISOTROPIC MATRIX: MATERIAL 'FICTICIOUS' ELASTICITY TENSOR
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      INTEGER NDI,I1,J1,K1,L1
      DOUBLE PRECISION CMISOMATFIC(ndi,ndi,ndi,ndi),UNIT2(ndi,ndi),
     1                 cbar(ndi,ndi),diso(5),
     2                 unit4(ndi,ndi,ndi,ndi)
      DOUBLE PRECISION CBARI1,CBARI2
      DOUBLE PRECISION DUDI1,DUDI2,D2UD2I1,D2UD2I2,D2UDI1I2
      DOUBLE PRECISION AUX,AUX1,AUX2,AUX3,AUX4,DET
      DOUBLE PRECISION UIJ,UKL,CIJ,CKL   
C
      dudi1=diso(1)
      dudi2=diso(2)
      d2ud2i1=diso(3)
      d2ud2i2=diso(4)
      d2udi1i2=diso(5)
C
      aux1=four*(d2ud2i1+two*cbari1*d2udi1i2+
     1           dudi2+cbari1*cbari1*d2ud2i2)
      aux2=-four*(d2udi1i2+cbari1*d2ud2i2)
      aux3=four*d2ud2i2
      aux4=-four*dudi2

      DO i1=1,ndi
        DO j1=1,ndi
           DO k1=1,ndi
              DO l1=1,ndi
                     uij=unit2(i1,j1)
                     ukl=unit2(k1,l1)
                     cij=cbar(i1,j1)
                     ckl=cbar(k1,l1)
                     aux=aux1*uij*ukl+
     1                   aux2*(uij*ckl+cij*ukl)+aux3*cij*ckl+
     3                   aux4*unit4(i1,j1,k1,l1)
                     cmisomatfic(i1,j1,k1,l1)=aux        
              END DO
           END DO
        END DO
      END DO
C
      RETURN
      END SUBROUTINE cmatisomatfic