Esope source files

   1 : C FLUDE1    SOURCE    AM        09/12/08    21:15:16     6582           
   2 :       SUBROUTINE FLUDE1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
   3 :      $                       AN,BN,CN,NDI)
   4 : C
   5 : C_______________________________________________________________________
   6 : C
   7 : C SOUS PROGRAMME QUI CALCUL LES MATRICES DE DEFORMATION 
   8 : C  DE FLUAGE DESSICCATION 
   9 : C   D APRES LE MODELE BAZANT
  10 : C
  11 : C
  12 : C_______________________________________________________________________
  13 : C
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8(A-H,O-Z)
  16 :       DIMENSION VARF(NVARI),CMAT(NMAT)
  17 :       DIMENSION AN(*),BN(NDI,*),CN(NDI,*)
  18 : C
  18 :       CALL ZERO(AN,NSTRS,1)
  18 :       CALL ZERO(BN,NDI,NSTRS)
  18 :       CALL ZERO(CN,NDI,NSTRS)
  19 : C
  20 : C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  21 : C
  22 :       XFDE = CMAT(25) 
  23 :       XNU  = CMAT(24)
  24 : C
  25 : C RECUPERATION DES TEMPERATURES
  26 : C
  27 :       TINI = CMAT(29)
  28 :       TFIN = CMAT(30)
  29 :       HINI = DEXP(0.018*TINI/(1000.*8.314*293.))
  30 :       HFIN = DEXP(0.018*TFIN/(1000.*8.314*293.))
  31 : C
  32 : C CONSTRUCTION DE LA MATRICE 
  33 : C
  34 :       ADR = 0.D0
  35 :       BDR = XFDE*( DABS(HFIN - HINI) )/2.D0
  36 :       CDR = BDR
  37 : C
  38 :       DO 10 I=1,NSTRS
  39 :           AN(I) = 0.D0
  40 :   10  CONTINUE
  41 : C
  42 :       IF(IFOUR.EQ.-2) THEN
  43 : C
  44 : C  - CONTRAINTES PLANES 
  45 : C
  46 :          BN(1,1)=1.D0
  47 :          BN(1,2)=-XNU
  48 :          BN(2,1)=BN(1,2)
  49 :          BN(2,2)=BN(1,1)
  50 :          BN(3,3)=1.D0+XNU
  51 :       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  52 : C
  53 : C  - DEFORMATION PLANE OU AXISYMETRIQUE 
  54 : C
  55 :          BN(1,1)=1.D0
  56 :          BN(1,2)=-XNU
  57 :          BN(1,3)=BN(1,2)
  58 :          BN(2,1)=BN(1,2)
  59 :          BN(2,2)=BN(1,1)
  60 :          BN(2,3)=BN(1,2)
  61 :          BN(3,1)=BN(1,2)
  62 :          BN(3,2)=BN(1,2)
  63 :          BN(3,3)=BN(1,1)
  64 :          BN(4,4)=1.D0+XNU
  65 :       ELSE
  66 : C
  67 : C  - CALCUL TRIDIMENSIONEL 
  68 : C
  69 :          BN(1,1)=1.D0
  70 :          BN(1,2)=-XNU
  71 :          BN(1,3)=BN(1,2)
  72 :          BN(2,1)=BN(1,2)
  73 :          BN(2,2)=BN(1,1)
  74 :          BN(2,3)=BN(1,2)
  75 :          BN(3,1)=BN(1,2)
  76 :          BN(3,2)=BN(1,2)
  77 :          BN(3,3)=BN(1,1)
  78 :          BN(4,4)=1.D0+XNU
  79 :          BN(5,5)=BN(4,4)
  80 :          BN(6,6)=BN(4,4)
  81 :       ENDIF
  82 : C
  83 :       DO 20 I=1,NSTRS
  84 :         DO 20 J=1,NSTRS
  85 :           CN(I,J) = CDR * BN(I,J)
  86 :           BN(I,J) = BDR * BN(I,J)
  87 :   20  CONTINUE
  88 : C
  89 :       END
  90 :  
  91 :