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