Esope source files

   1 : C HUJFLO    SOURCE    CHAT      05/01/13    00:26:48     5004           
   2 :       SUBROUTINE HUJFLO(IFOUR,SIG,XINV,PI,XKSI,XMAT,DFDS,
   3 :      .                  DFDPI,DFDKSI,DGDS,HV,HQ)
   4 : *----------------------------------------------------------------
   5 : *                      MODELE HUJEUX
   6 : *            QUANTITES NECESSAIRES A L'ECOULEMENT
   7 : *----------------------------------------------------------------
   8 : *   ENTREE :
   9 : *          IFOUR     CF CCOPTIO
  10 : *          SIG       CONTRAINTES
  11 : *          PI,XKSI   VARIABLES INTERNES
  12 : *          XMAT      PROPRIETES MATERIELLES
  13 : *
  14 : *   SORTIE :
  15 : *          DFDS      DF / DSIG
  16 : *          DGDS      DG / DSIG
  17 : *          DFDPI     DF / DPI
  18 : *          HV,HQ     LOIS D'ECROUISSAGE
  19 : *----------------------------------------------------------------
  20 : *
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  23 :       DIMENSION SIG(*),XINV(*),XMAT(*),DFDS(*),DGDS(*)
  24 : *
  25 :       COHE = XMAT(6)
  26 :       GM   = XMAT(7)
  27 :       BETA = XMAT(8)
  28 :       HUA  = XMAT(9)
  29 :       HUB  = XMAT(10)
  30 :       HUR  = XMAT(11)
  31 : *
  32 :       UNSP=1.D0/(COHE-XINV(1))
  33 :       IF(XINV(2).EQ.0.D0) THEN
  34 :          UNIT=0.D0
  35 :          UNSQ=0.D0
  36 :       ELSE
  37 :          UNIT=1.D0
  38 :          UNSQ=1.D0/(XINV(2)*GM)
  39 :       ENDIF
  40 :       DENOM=(PI+COHE)*EXP(-1.D0/HUB)
  41 :       RAP=(COHE-XINV(1))/DENOM
  42 :       XLRAP=LOG(RAP)
  43 : *
  44 :       FAC=(UNSP*XINV(2)/GM - HUB*XKSI)/3.D0
  45 :       DFDS(1)=UNSP*(UNSQ*(SIG(1)-0.5D0*(SIG(2)+SIG(3))) + FAC)
  46 :       DFDS(2)=UNSP*(UNSQ*(SIG(2)-0.5D0*(SIG(1)+SIG(3))) + FAC)
  47 :       DFDS(3)=UNSP*(UNSQ*(SIG(3)-0.5D0*(SIG(1)+SIG(2))) + FAC)
  48 :       DFDS(4)=UNIT*UNSP*UNSQ*1.5D0*SIG(4)
  49 :       IF(IFOUR.EQ.2) THEN
  50 :          DFDS(5)=UNIT*UNSP*UNSQ*1.5D0*SIG(5)
  51 :          DFDS(6)=UNIT*UNSP*UNSQ*1.5D0*SIG(6)
  52 :       ENDIF
  53 : *
  54 :       FAC=(UNSP*XINV(2)/GM - 1.D0)/3.D0
  55 :       DGDS(1)=UNSP*(UNSQ*(SIG(1)-0.5D0*(SIG(2)+SIG(3))) + FAC)
  56 :       DGDS(2)=UNSP*(UNSQ*(SIG(2)-0.5D0*(SIG(1)+SIG(3))) + FAC)
  57 :       DGDS(3)=UNSP*(UNSQ*(SIG(3)-0.5D0*(SIG(1)+SIG(2))) + FAC)
  58 :       DGDS(4)=DFDS(4)
  59 :       IF(IFOUR.EQ.2) THEN
  60 :          DGDS(5)=DFDS(5)
  61 :          DGDS(6)=DFDS(6)
  62 :       ENDIF
  63 : *
  64 :       DFDPI =-HUB*XKSI/(PI+COHE)
  65 :       DFDKSI= HUB*XLRAP - 1.D0
  66 : *
  67 :       HV=-BETA*UNSP*(PI+COHE)*(1.D0-UNSP*XINV(2)/GM)
  68 :       IF(HUA.EQ.0.D0) THEN
  69 :          HQ=0.D0
  70 :       ELSE
  71 :          HQ= UNSP*((HUR+1.D0-XKSI)**2)/(GM*HUA)
  72 :       ENDIF
  73 : *
  74 :       RETURN                                                                    
  75 :       END                                                                       
  76 :  
  77 :  
  78 :