Télécharger modbet.eso

Retour à la liste

modbet
C MODBET    SOURCE    CB215821  17/11/30    21:16:50     9639           
      SUBROUTINE MODBET(TP1,TP2,SE1,SE2,EFICU,EFICUL,EFIC,CODU,CODL,
     &                  COD,BETJEF,BETFLU)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      CHARACTER*40 FMT,TITRE
      DIMENSION BRAN(8),E1(9),E2(9),EFICU(9),EFICUL(8)
      DIMENSION CODU(9,9),CODL(8,8),COD(8)
C
*
      SEGMENT BETJEF
       REAL*8 AA,BETA,FC,PALF,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
     &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
       INTEGER ICT,ICC,IMOD,IVIS,ITR,
     &                ISIM,IBB,IGAU,IZON
      ENDSEGMENT
      SEGMENT BETFLU
       REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
     &         TP0,TZER
       INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
      ENDSEGMENT
C
C*******************************************************************
C     CALCUL DES COEFFICIENTS DES BRANCHES DU MODELE DE MAXWELL
C*******************************************************************
C
C
      MC=NBRC+1
C
C
      IF (IFOR.EQ.1) THEN
C     Modele visco elasto plastique SOLIDE de MAXWELL
C*******************************************************************
C     INITIALISATION
C*******************************************************************
C
      EFIC=0.D0
      SE1=0.D0
      SE2=0.D0
      ETPS1=0.D0
      ETPS2=0.D0
C
      DO 15 K = 1,MC
      EFICU(K) = 0.D0
      E1(K) = 0.D0
      E2(K) = 0.D0
   15 CONTINUE
C
C*******************************************************************
C
      DO 20 I=1,MC
C
      IF(I.EQ.1) THEN
C
      CALL MODMAX(TP1,ETPS1,1,CODU,CODL,COD,BETJEF,BETFLU)
      CALL MODMAX(TP2,ETPS2,1,CODU,CODL,COD,BETJEF,BETFLU)
C
      BRAN(I) = TAU1
C
      EFICU(1)=ETPS1+(ETPS2-ETPS1)/2
      E1(I)=ETPS1
      E2(I)=ETPS2
      ETPS1=0.D0
      ETPS2=0.D0
C
      SE1=SE1+E1(I)
      SE2=SE2+E2(I)
C
      ELSE IF (I.GT.1) THEN
C
      CALL MODMAX(TP1,ETPS1,I,CODU,CODL,COD,BETJEF,BETFLU)
      CALL MODMAX(TP2,ETPS2,I,CODU,CODL,COD,BETJEF,BETFLU)
C
      BRAN(I)=10**(I-2)*TAU2
C
      E1(I)=ETPS1
      E2(I)=ETPS2
C
      SE1=SE1+E1(I)
      SE2=SE2+E2(I)
C
      IF ((TP2-TP1).EQ.0.D0) THEN
      EFICU(I)=ETPS1
      ELSE
      QOEF=BRAN(I-1)/((TP2-TP1)/86400)
      AOEF=1.D0-EXP(-1/QOEF)
      EFICU(I)=QOEF*(ETPS1*AOEF+(ETPS2-ETPS1)
     **(1.D0-QOEF*AOEF))
      ENDIF
C
      ETPS1=0.D0
      ETPS2=0.D0
      QOEF=0.D0
      AOEF=0.D0
C
      ENDIF
C
      EFIC=EFIC+EFICU(I)
C
   20 CONTINUE
C
C------------------------------------------------------------------
      ELSE IF (IFOR.EQ.2) THEN
C     Modele visco elasto plastique LIQUIDE de MAXWELL
C--------------------------------------------------------------------
C
C*******************************************************************
C     INITIALISATION
C*******************************************************************
C
      EFIC=0.D0
      SE1=0.D0
      SE2=0.D0
      ETPS1=0.D0
      ETPS2=0.D0
C
      DO 25 K = 1,MC
      E1(K) = 0.D0
      E2(K) = 0.D0
   25 CONTINUE
      DO 26 K = 1,NBRC
      EFICUL(K) = 0.D0
   26 CONTINUE
C
C*******************************************************************
C
      DO 30 I=1,NBRC
C
      CALL MODMAX(TP1,ETPS1,I,CODU,CODL,COD,BETJEF,BETFLU)
      CALL MODMAX(TP2,ETPS2,I,CODU,CODL,COD,BETJEF,BETFLU)
C
      BRAN(I) = 0.D0
      IF (I.EQ.1) THEN
      BRAN(I) = TAU1
      ELSE
      BRAN(I) = 10**(I-2)*TAU2
      ENDIF
C
      E1(I)=ETPS1
      E2(I)=ETPS2
C
      SE1=SE1+E1(I)
      SE2=SE2+E2(I)
C
      IF ((TP2-TP1).EQ.0.D0) THEN
C
      EFICUL(I)=ETPS1
      ELSE
C
      QOEF=BRAN(I)/((TP2-TP1)/86400)
      AOEF=1.D0-EXP(-1/QOEF)
      EFICUL(I)=QOEF*(ETPS1*AOEF+(ETPS2-ETPS1)
     **(1.D0-QOEF*AOEF))
C
      ENDIF
C
C
      ETPS1=0.D0
      ETPS2=0.D0
      QOEF=0.D0
      AOEF=0.D0
C
C
      EFIC=EFIC+EFICUL(I)
C
   30 CONTINUE
C
C--------------------------------------------------------------
      ENDIF
C
      RETURN
      END
 
 
 
 

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