Esope source files

   1 : C ADENS     SOURCE    CHAT      05/01/12    21:17:54     5004           
   2 :        SUBROUTINE ADENS(XMAT,AA,KERRE)
   3 : C
   4 : C ==================================================================
   5 : C ENTREES:
   6 : C -------
   7 : C XMAT(NCOMAT) = CARACTERISTIQUES THERMOMECANIQUES DU MATERIAU
   8 : C
   9 : C SORTIES:
  10 : C -------
  11 : C AA           = COEFFICIENT INTERVENANT DANS LE CALCUL DE L INCREMENT
  12 : C                DE LA DEFORMATION DE DENSIFICATION
  13 : C KERRE        = GESTION DES ERREURS
  14 : C ==================================================================
  15 :       IMPLICIT INTEGER(I-N)
  16 :        IMPLICIT REAL*8(A-H,O-Z)                                                 
  17 :        DIMENSION XMAT(*)                               
  18 : C
  19 :         ADEN  = XMAT(23)
  20 :         XKGON = XMAT(24)
  21 :         RHO0  = 1.D0-XMAT(25)
  22 :         BUMI  = XMAT(26)
  23 : C
  24 :         KERRE=0
  25 : C
  26 : C---------------------------------------------------------------------
  27 : C Calcul du coefficient A intervenant dans la vitesse de densification 
  28 : C---------------------------------------------------------------------
  29 :         BUMAX0=60.D0*BUMI
  30 :         CRIT=1.D-10
  31 :         DO 50 I=1,150          
  32 :           BUMAX = BUMI*EXP(1.D0-(ADEN/(XKGON*BUMAX0)))
  33 :           BUMAX = (0.2D0*BUMAX) + (0.8D0*BUMAX0)
  34 :           TEST=ABS((BUMAX-BUMAX0)/BUMAX0)
  35 :           IF (TEST.LT.CRIT) GOTO 60 
  36 :           BUMAX0=BUMAX
  37 :  50     CONTINUE
  38 :         KERRE=2
  39 :         PRINT *,'&&&& ADENS: IMPOSSIBLE DE CALCULER LE COEFFICIENT AA'
  40 :         PRINT *,'&&&& INTERVENANT DANS LA VITESSE DE DENSIFICATION'
  41 :         RETURN
  42 :  60     CONTINUE
  43 :         AA = (-1.D0*RHO0*(ADEN-(XKGON*BUMAX))) /
  44 :      &         ((1.D0+ADEN)*(DLOG(BUMAX/BUMI)))
  45 : C
  46 :         RETURN
  47 :         END
  48 :  
  49 :