Télécharger densgf.eso

Retour à la liste

Numérotation des lignes :

  1. C DENSGF SOURCE CB215821 16/04/21 21:16:23 8920
  2. SUBROUTINE DENSGF(BU,BUPT,XMAT,NCOMAT,AA,DEPSD,DEPSG)
  3. C
  4. C ==================================================================
  5. C ENTREES:
  6. C -------
  7. C BU = TAUX DE COMBUSTION INSTANCIE
  8. C BUPT = VITESSE DU TAUX DE COMBUSTION
  9. C XMAT(NCOMAT) = CARACTERISTIQUES THERMOMECANIQUES DU MATERIAU
  10. C AA = COEFFICIENT INTERVENANT DANS LE CALCUL DE L INCREMENT
  11. C DE LA DEFORMATION DE DENSIFICATION
  12. C
  13. C SORTIES:
  14. C -------
  15. C DEPSD(3) = VITESSE DE LA DEFORMATION DE DENSIFICATION
  16. C DEPSG(3) = VITESSE DE LA DEFORMATION DE GONFLEMENT
  17. C ==================================================================
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. DIMENSION XMAT(*),DEPSD(*),DEPSG(*)
  21. C
  22. XKGON=XMAT(24)
  23. RHO0=1.D0-XMAT(25)
  24. BUMI=XMAT(26)
  25. C
  26. C------------------------------------------------------------
  27. C VITESSES DES DEFORMATIONS DE DENSIFICATION ET DE GONFLEMENT
  28. C------------------------------------------------------------
  29. C
  30. IF (BU.LT.BUMI) THEN
  31. EPSDPT=0.D0
  32. ELSE
  33. XLN=LOG(BU/BUMI)
  34. DEN=RHO0+(AA*XLN)
  35. XN1=DEN*XKGON*AA*XLN
  36. XN2=AA*RHO0*(1.D0+(XKGON*BU))/BU
  37. XND=(XN1+XN2)/(DEN*DEN)
  38. EPSDPT=XND*BUPT
  39. EPSDPT=-1.D0*(EPSDPT/3.D0)
  40. ENDIF
  41. EPSGPT=XKGON*BUPT/3.D0
  42. DO 70 I=1,3
  43. DEPSD(I)=EPSDPT
  44. DEPSG(I)=EPSGPT
  45. 70 CONTINUE
  46. C
  47. RETURN
  48. END
  49.  
  50.  
  51.  
  52.  
  53.  

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