Esope source files

   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
  55 :       CALL MODMAX(TP1,ETPS1,1,CODU,CODL,COD,BETJEF,BETFLU)
  55 :       CALL MODMAX(TP2,ETPS2,1,CODU,CODL,COD,BETJEF,BETFLU)
  56 : C
  57 :       BRAN(I) = TAU1
  58 : C
  59 :       EFICU(1)=ETPS1+(ETPS2-ETPS1)/2
  60 :       E1(I)=ETPS1
  61 :       E2(I)=ETPS2
  62 :       ETPS1=0.D0
  63 :       ETPS2=0.D0
  64 : C
  65 :       SE1=SE1+E1(I)
  66 :       SE2=SE2+E2(I)
  67 : C
  68 :       ELSE IF (I.GT.1) THEN
  69 : C
  69 :       CALL MODMAX(TP1,ETPS1,I,CODU,CODL,COD,BETJEF,BETFLU)
  69 :       CALL MODMAX(TP2,ETPS2,I,CODU,CODL,COD,BETJEF,BETFLU)
  70 : C
  71 :       BRAN(I)=10**(I-2)*TAU2
  72 : C
  73 :       E1(I)=ETPS1
  74 :       E2(I)=ETPS2
  75 : C
  76 :       SE1=SE1+E1(I)
  77 :       SE2=SE2+E2(I)
  78 : C
  79 :       IF ((TP2-TP1).EQ.0.D0) THEN
  80 :       EFICU(I)=ETPS1
  81 :       ELSE
  82 :       QOEF=BRAN(I-1)/((TP2-TP1)/86400)
  83 :       AOEF=1.D0-EXP(-1/QOEF)
  84 :       EFICU(I)=QOEF*(ETPS1*AOEF+(ETPS2-ETPS1)
  85 :      **(1.D0-QOEF*AOEF))
  86 :       ENDIF
  87 : C
  88 :       ETPS1=0.D0
  89 :       ETPS2=0.D0
  90 :       QOEF=0.D0
  91 :       AOEF=0.D0
  92 : C      
  93 :       ENDIF
  94 : C
  95 :       EFIC=EFIC+EFICU(I)        
  96 : C
  97 :    20 CONTINUE
  98 : C
  99 : C------------------------------------------------------------------
 100 :       ELSE IF (IFOR.EQ.2) THEN
 101 : C     Modele visco elasto plastique LIQUIDE de MAXWELL   
 102 : C--------------------------------------------------------------------
 103 : C
 104 : C*******************************************************************
 105 : C     INITIALISATION
 106 : C*******************************************************************
 107 : C
 108 :       EFIC=0.D0         
 109 :       SE1=0.D0
 110 :       SE2=0.D0
 111 :       ETPS1=0.D0
 112 :       ETPS2=0.D0
 113 : C
 114 :       DO 25 K = 1,MC
 115 :       E1(K) = 0.D0
 116 :       E2(K) = 0.D0
 117 :    25 CONTINUE
 118 :       DO 26 K = 1,NBRC
 119 :       EFICUL(K) = 0.D0
 120 :    26 CONTINUE   
 121 : C   
 122 : C*******************************************************************  
 123 : C   
 124 :       DO 30 I=1,NBRC
 125 : C
 125 :       CALL MODMAX(TP1,ETPS1,I,CODU,CODL,COD,BETJEF,BETFLU)
 125 :       CALL MODMAX(TP2,ETPS2,I,CODU,CODL,COD,BETJEF,BETFLU)
 126 : C
 127 :       BRAN(I) = 0.D0
 128 :       IF (I.EQ.1) THEN
 129 :       BRAN(I) = TAU1
 130 :       ELSE
 131 :       BRAN(I) = 10**(I-2)*TAU2
 132 :       ENDIF
 133 : C
 134 :       E1(I)=ETPS1
 135 :       E2(I)=ETPS2
 136 : C
 137 :       SE1=SE1+E1(I)
 138 :       SE2=SE2+E2(I)
 139 : C
 140 :       IF ((TP2-TP1).EQ.0.D0) THEN
 141 : C      
 142 :       EFICUL(I)=ETPS1
 143 :       ELSE
 144 : C      
 145 :       QOEF=BRAN(I)/((TP2-TP1)/86400)
 146 :       AOEF=1.D0-EXP(-1/QOEF)
 147 :       EFICUL(I)=QOEF*(ETPS1*AOEF+(ETPS2-ETPS1)
 148 :      **(1.D0-QOEF*AOEF))
 149 : C     
 150 :       ENDIF
 151 : C
 152 : C
 153 :       ETPS1=0.D0
 154 :       ETPS2=0.D0
 155 :       QOEF=0.D0
 156 :       AOEF=0.D0
 157 : C      
 158 : C
 159 :       EFIC=EFIC+EFICUL(I)       
 160 : C
 161 :    30 CONTINUE
 162 : C
 163 : C--------------------------------------------------------------
 164 :       ENDIF
 165 : C
 166 :       RETURN
 167 :       END
 168 :  
 169 :