Esope source files

   1 : C CAMECA    SOURCE    AM        09/12/08    21:15:01     6582           
   2 :       SUBROUTINE CAMECA(NSTRS,DEPST,VAR0,NVARI,CMAT,NMAT,
   3 :      *            IFOUR,IB,IGAU,SIGF,VARF,
   4 :      * SIGI,DEP,EPSILO,EPSEQU,EPSFLU,DSIGT,SIGEF,AN,BN,CN,
   5 :      * KERRE)
   6 : 
   7 : C
   8 : C ROUTINE APPELE DANS FLUTRA 
   9 : C___________________________________________________________________
  10 : C
  11 : C             MODELE BETON
  12 : C                   FLUAGE 
  13 : C
  14 : C___________________________________________________________________
  15 : C
  16 :       IMPLICIT INTEGER(I-N)
  17 :       IMPLICIT REAL*8(A-H,O-Z)
  18 :       DIMENSION CMAT(NMAT)
  19 :       DIMENSION VAR0(NVARI),VARF(NVARI)
  20 : *
  21 :       DIMENSION SIGI(NSTRS),DEP(NSTRS,NSTRS)
  22 :       DIMENSION DEPST(NSTRS),EPSILO(NSTRS)
  23 :       DIMENSION EPSEQU(NSTRS),EPSFLU(NSTRS)
  24 :       DIMENSION DSIGT(NSTRS),SIGF(NSTRS),SIGEF(NSTRS)
  25 :       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  26 : C
  27 : C______________________________________________________
  28 : C Recuperation variable & intialisation
  29 : C______________________________________________________
  30 : C
  30 :       CALL ZERO(AN,NSTRS,1)
  30 :       CALL ZERO(BN,NSTRS,NSTRS)
  30 :       CALL ZERO(CN,NSTRS,NSTRS)
  30 :       CALL ZERO(EPSFLU,NSTRS,1)
  31 : C
  32 :       TDT  = CMAT(NMAT)
  33 :       BETA = CMAT(22)
  34 : C
  35 : C RECUPERATION DES VARIABLES INTERNES INITIALES
  36 : C
  37 :       SIGEF(1) = VAR0(3)
  38 :       SIGEF(2) = VAR0(4)
  39 : C
  40 : C RECUPERATION DES CONTRAINTES INITIALES EFFECTIVES
  41 : C
  42 :       IF (IFOUR.EQ.-2) THEN
  43 :         SIGEF(3) = VAR0(5)
  44 :       ELSE
  45 :         SIGEF(3) = VAR0(6)
  46 :         SIGEF(4) = VAR0(5)
  47 :         IF (IFOUR.EQ.2) THEN
  48 :           SIGEF(5) = VAR0(7)
  49 :           SIGEF(6) = VAR0(8)
  50 :         ENDIF
  51 :       ENDIF
  52 : C
  53 : C_______________________________________________________
  54 : C
  55 : C CONSTRUCTION DES MATRICES DE COMPLAISANCE ET ELASTIQUE
  56 : C_______________________________________________________
  57 : C
  58 : C Matrice elastique
  59 : C
  59 :       CALL MATLAS(CMAT,NMAT,NSTRS,IFOUR,DEP,KERRE)
  60 : C
  61 : C Matrice de fluage
  62 : C
  63 : C   DFLUT = AN + BN * SIGMA_N + CN * SIGMA_N+1
  64 : C
  65 :       IF (TDT.NE.0.D0) THEN
  65 :         CALL FLUPR1(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
  66 :      $              AN,BN,CN,KERRE)
  67 :       ENDIF
  68 : C
  69 : C_______________________________________________
  70 : C
  71 : C  CALCUL DE L INCREMENT DE CONTRAINTE EFFECTIVE
  72 : C_______________________________________________
  73 : C
  73 :       CALL ELASF1(DEP,AN,BN,CN,DEPST,SIGEF,
  74 :      $             NSTRS,DSIGT,KERRE)
  75 : C
  76 : C___________________________________________
  77 : C
  78 : C  CALCUL DES CONTRAINTES EFFECTIVES FINALES
  79 : C___________________________________________
  80 : C
  81 :       DO 10 I=1,NSTRS
  82 :         SIGF(I) = SIGEF(I) + DSIGT(I)
  83 :    10 CONTINUE 
  84 : C__________________________________________________
  85 : C
  86 : C SAUVEGARDE DES DEFORMATIONS DE FLUAGE
  87 : C  ET CALCUL DE LA DEFORMATION DE FLUAGE
  88 : C   A LA FIN DU PAS DE TEMPS EPSFLU
  89 : C__________________________________________________
  90 : C
  91 :       IF (TDT.NE.0.D0) THEN
  91 :         CALL SAUFL1(VARF,NVARI,CMAT,NMAT,SIGF,
  92 :      $     SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,KERRE)
  93 :       ENDIF
  94 : C__________________________________________________
  95 : C
  96 : C SAUVEGARDE DES CONTRAINTES EFFECTIVES
  97 : C__________________________________________________
  98 : C
  99 :       VARF(3) = SIGF(1)
 100 :       VARF(4) = SIGF(2)
 101 :       IF (IFOUR.EQ.-2) THEN
 102 :         VARF(5) = SIGF(3)
 103 :       ELSE
 104 :         VARF(6) = SIGF(3)
 105 :         VARF(5) = SIGF(4)
 106 :         IF (IFOUR.EQ.2) THEN
 107 :           VARF(7) = SIGF(5) 
 108 :           VARF(8) = SIGF(6)
 109 :         ENDIF
 110 :       ENDIF
 111 : C__________________________________________________
 112 : C
 113 : C  CALCUL DE LA DEFORMATION ELASTIQUE EN FIN DE PAS
 114 : C__________________________________________________
 115 : C
 115 :       CALL DEFELA(DEP,EPSILO,NSTRS,SIGF,IFOUR,KERRE)
 116 : C
 117 : C__________________________________________________
 118 : C
 119 : C  CALCUL DE LA DEFORMATION EQUIVALENTE
 120 : C__________________________________________________
 121 : C
 122 :       DO 20 I=1,NSTRS
 123 :         EPSEQU(I) = EPSILO(I) + ( BETA * EPSFLU(I) )
 124 :    20 CONTINUE 
 125 : C
 126 : C____________________________________________________________
 127 : C
 128 : C  CALCUL ENDOMMAGEMENT ET DES CONTRAINTES APPARENTES FINALES
 129 : C____________________________________________________________
 130 : C
 131 : C ENDOMMAGEMENT ET CONTRAINTES APPARENTES
 132 : C
 132 :       CALL ENDO(EPSEQU,NSTRS,NVARI,CMAT,DEP,SIGF,VAR0,
 133 :      $          VARF,NMAT,IFOUR)
 134 : C
 135 :       END
 136 :  
 137 :