Esope source files

   1 : C DENSGF    SOURCE    CHAT      05/01/12    22:41:57     5004           
   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=DLOG(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 :