Télécharger la plaquette Cast3M


fr en

Télécharger densgf.eso

Retour à la liste

Numérotation des lignes : 

densgf

  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