Télécharger la plaquette Cast3M


fr en

Télécharger adens.eso

Retour à la liste

Numérotation des lignes : 

adens

  1. C ADENS     SOURCE    CB215821  16/04/21    21:15:05     8920
  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)*(LOG(BUMAX/BUMI)))
  45. C
  46.         RETURN
  47.         END
  48.  
  49.  
  50.  
  51.  
  52.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales