Esope source files

   1 : C FLUA      SOURCE    AM        09/12/08    21:15:14     6582           
   2 :       SUBROUTINE FLUA(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 PROPRE 
   9 : C   D APRES LE MODELE AMORTISSEUR
  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 :       VIDV = CMAT(23) 
  21 :       TDT  = CMAT(NMAT)
  22 :       XNU  = CMAT(24)
  23 : C
  24 : C RECUPERATION DES VARIABLES INTERNES INITIALES
  25 : C
  26 : C    REMPLISSAGE DES VECTEURS DE DEFORMATION DE FLUAGE
  27 : C
  28 :       EPSFLU(1) = VARF(45)
  29 :       EPSFLU(2) = VARF(46)
  30 :       EPSFLU(3) = VARF(47)
  31 :       EPSFLU(4) = VARF(48)
  32 :       EPSFLU(5) = VARF(49)
  33 :       EPSFLU(6) = VARF(50)
  34 : C
  35 : C
  36 : C CONSTRUCTION DE LA MATRICE 
  37 : C
  38 :       ADR = 0.D0
  39 :       BDR = TDT/(2.D0*VIDV)
  40 :       CDR = TDT/(2.D0*VIDV)
  41 : C
  42 :       DO 10 I=1,NSTRS
  43 :           AN(I) = ADR * EPSFLU(I)
  44 :   10  CONTINUE
  45 : C
  46 :       IF(IFOUR.EQ.-2) THEN
  47 : C
  48 : C  - CONTRAINTES PLANES 
  49 : C
  50 :          BN(1,1)=1.D0
  51 :          BN(1,2)=-XNU
  52 :          BN(2,1)=BN(1,2)
  53 :          BN(2,2)=BN(1,1)
  54 :          BN(3,3)=1.D0+XNU
  55 :       ELSE IF((IFOUR.EQ.-1).OR.(IFOUR.EQ.0)) THEN
  56 : C
  57 : C  - DEFORMATION PLANE OU AXISYMETRIQUE 
  58 : C
  59 :          BN(1,1)=1.D0
  60 :          BN(1,2)=-XNU
  61 :          BN(1,3)=BN(1,2)
  62 :          BN(2,1)=BN(1,2)
  63 :          BN(2,2)=BN(1,1)
  64 :          BN(2,3)=BN(1,2)
  65 :          BN(3,1)=BN(1,2)
  66 :          BN(3,2)=BN(1,2)
  67 :          BN(3,3)=BN(1,1)
  68 :          BN(4,4)=1.D0+XNU
  69 :       ELSE
  70 : C
  71 : C  - CALCUL TRIDIMENSIONEL 
  72 : C
  73 :          BN(1,1)=1.D0
  74 :          BN(1,2)=-XNU
  75 :          BN(1,3)=BN(1,2)
  76 :          BN(2,1)=BN(1,2)
  77 :          BN(2,2)=BN(1,1)
  78 :          BN(2,3)=BN(1,2)
  79 :          BN(3,1)=BN(1,2)
  80 :          BN(3,2)=BN(1,2)
  81 :          BN(3,3)=BN(1,1)
  82 :          BN(4,4)=1.D0+XNU
  83 :          BN(5,5)=BN(4,4)
  84 :          BN(6,6)=BN(4,4)
  85 :       ENDIF
  86 : C
  87 :       DO 20 I=1,NSTRS
  88 :         DO 20 J=1,NSTRS
  89 :           CN(I,J) = CDR * BN(I,J)
  90 :           BN(I,J) = BDR * BN(I,J)
  91 :   20  CONTINUE
  92 : C
  93 :       END
  94 :  
  95 :