Télécharger hujflo.eso

Retour à la liste

Numérotation des lignes :

  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.  
  79.  

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