Esope source files

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