Télécharger modbet.eso

Retour à la liste

Numérotation des lignes :

  1. C MODBET SOURCE CHAT 05/01/13 01:47:26 5004
  2. SUBROUTINE MODBET(TP1,TP2,SE1,SE2,EFICU,EFICUL,EFIC,CODU,CODL,
  3. & COD,BETJEF,BETFLU)
  4. C
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. CHARACTER*40 FMT,TITRE
  8. DIMENSION BRAN(8),E1(9),E2(9),EFICU(9),EFICUL(8)
  9. DIMENSION CODU(9,9),CODL(8,8),COD(8)
  10. C
  11. *
  12. SEGMENT BETJEF
  13. REAL*8 AA,BETA,FC,PALF,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  14. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  15. INTEGER ICT,ICC,IMOD,IVIS,ITR,
  16. & ISIM,IBB,IGAU,IZON
  17. ENDSEGMENT
  18. SEGMENT BETFLU
  19. REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  20. & TP0,TZER
  21. INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  22. ENDSEGMENT
  23. C
  24. C*******************************************************************
  25. C CALCUL DES COEFFICIENTS DES BRANCHES DU MODELE DE MAXWELL
  26. C*******************************************************************
  27. C
  28. C
  29. MC=NBRC+1
  30. C
  31. C
  32. IF (IFOR.EQ.1) THEN
  33. C Modele visco elasto plastique SOLIDE de MAXWELL
  34. C*******************************************************************
  35. C INITIALISATION
  36. C*******************************************************************
  37. C
  38. EFIC=0.D0
  39. SE1=0.D0
  40. SE2=0.D0
  41. ETPS1=0.D0
  42. ETPS2=0.D0
  43. C
  44. DO 15 K = 1,MC
  45. EFICU(K) = 0.D0
  46. E1(K) = 0.D0
  47. E2(K) = 0.D0
  48. 15 CONTINUE
  49. C
  50. C*******************************************************************
  51. C
  52. DO 20 I=1,MC
  53. C
  54. IF(I.EQ.1) THEN
  55. C
  56. CALL MODMAX(TP1,ETPS1,1,CODU,CODL,COD,BETJEF,BETFLU)
  57. CALL MODMAX(TP2,ETPS2,1,CODU,CODL,COD,BETJEF,BETFLU)
  58. C
  59. BRAN(I) = TAU1
  60. C
  61. EFICU(1)=ETPS1+(ETPS2-ETPS1)/2
  62. E1(I)=ETPS1
  63. E2(I)=ETPS2
  64. ETPS1=0.D0
  65. ETPS2=0.D0
  66. C
  67. SE1=SE1+E1(I)
  68. SE2=SE2+E2(I)
  69. C
  70. ELSE IF (I.GT.1) THEN
  71. C
  72. CALL MODMAX(TP1,ETPS1,I,CODU,CODL,COD,BETJEF,BETFLU)
  73. CALL MODMAX(TP2,ETPS2,I,CODU,CODL,COD,BETJEF,BETFLU)
  74. C
  75. BRAN(I)=10**(I-2)*TAU2
  76. C
  77. E1(I)=ETPS1
  78. E2(I)=ETPS2
  79. C
  80. SE1=SE1+E1(I)
  81. SE2=SE2+E2(I)
  82. C
  83. IF ((TP2-TP1).EQ.0.D0) THEN
  84. EFICU(I)=ETPS1
  85. ELSE
  86. QOEF=BRAN(I-1)/((TP2-TP1)/86400)
  87. AOEF=1.D0-EXP(-1/QOEF)
  88. EFICU(I)=QOEF*(ETPS1*AOEF+(ETPS2-ETPS1)
  89. **(1.D0-QOEF*AOEF))
  90. ENDIF
  91. C
  92. ETPS1=0.D0
  93. ETPS2=0.D0
  94. QOEF=0.D0
  95. AOEF=0.D0
  96. C
  97. ENDIF
  98. C
  99. EFIC=EFIC+EFICU(I)
  100. C
  101. 20 CONTINUE
  102. C
  103. C------------------------------------------------------------------
  104. ELSE IF (IFOR.EQ.2) THEN
  105. C Modele visco elasto plastique LIQUIDE de MAXWELL
  106. C--------------------------------------------------------------------
  107. C
  108. C*******************************************************************
  109. C INITIALISATION
  110. C*******************************************************************
  111. C
  112. EFIC=0.D0
  113. SE1=0.D0
  114. SE2=0.D0
  115. ETPS1=0.D0
  116. ETPS2=0.D0
  117. C
  118. DO 25 K = 1,MC
  119. E1(K) = 0.D0
  120. E2(K) = 0.D0
  121. 25 CONTINUE
  122. DO 26 K = 1,NBRC
  123. EFICUL(K) = 0.D0
  124. 26 CONTINUE
  125. C
  126. C*******************************************************************
  127. C
  128. DO 30 I=1,NBRC
  129. C
  130. CALL MODMAX(TP1,ETPS1,I,CODU,CODL,COD,BETJEF,BETFLU)
  131. CALL MODMAX(TP2,ETPS2,I,CODU,CODL,COD,BETJEF,BETFLU)
  132. C
  133. BRAN(I) = 0.D0
  134. IF (I.EQ.1) THEN
  135. BRAN(I) = TAU1
  136. ELSE
  137. BRAN(I) = 10**(I-2)*TAU2
  138. ENDIF
  139. C
  140. E1(I)=ETPS1
  141. E2(I)=ETPS2
  142. C
  143. SE1=SE1+E1(I)
  144. SE2=SE2+E2(I)
  145. C
  146. IF ((TP2-TP1).EQ.0.D0) THEN
  147. C
  148. EFICUL(I)=ETPS1
  149. ELSE
  150. C
  151. QOEF=BRAN(I)/((TP2-TP1)/86400)
  152. AOEF=1.D0-EXP(-1/QOEF)
  153. EFICUL(I)=QOEF*(ETPS1*AOEF+(ETPS2-ETPS1)
  154. **(1.D0-QOEF*AOEF))
  155. C
  156. ENDIF
  157. C
  158. C
  159. ETPS1=0.D0
  160. ETPS2=0.D0
  161. QOEF=0.D0
  162. AOEF=0.D0
  163. C
  164. C
  165. EFIC=EFIC+EFICUL(I)
  166. C
  167. 30 CONTINUE
  168. C
  169. C--------------------------------------------------------------
  170. ENDIF
  171. C
  172. RETURN
  173. END
  174.  
  175.  
  176.  

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