Esope source files

   1 : C MOBFLU    SOURCE    CHAT      05/01/13    01:46:47     5004           
   2 : C MOBFLU     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE MOBFLU(SIGR,STFL,D1,NSTRS,IFLU,RMOY,PHIF,TEMP0,
   4 :      & DTEMP,TEMP1,POAR)
   5 : C
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 :       DIMENSION SIGR(NSTRS),STFL(NSTRS),D1(NSTRS,NSTRS)
   9 : C----------------------------------------------------------------------
  10 : *
  11 :       IF(IFLU.EQ.0) GOTO 10
  12 :       IF(IFLU.EQ.10.OR.IFLU.EQ.11.OR.IFLU.EQ.20.OR.IFLU.EQ.21) GOTO 20
  13 :       IF(IFLU.EQ.30.OR.IFLU.EQ.31) GOTO 30
  14 : *
  15 :       WRITE(*,*) '!! ATTENTION CE MODELE DE FLUAGE N EST PAS DISPONIBLE'
  16 :       WRITE(*,*) '      (DANS MOBFLU)  IFLU= ',IFLU
  17 :       STOP
  18 : c+mdj
  19 : C      GOTO 100
  20 : c+mdj
  21 : *
  22 :    10 CONTINUE
  23 : *
  24 : *    CE MATERIAUX NE FLUE PAS 
  25 : *
  26 :         GOTO 100
  27 : *
  28 :    20 CONTINUE
  29 : *
  30 : *    MODELE DE FLUAGE BETON : LOI BPEL91
  31 : *
  32 : *       IFLU = 10 : FLUAGE (BPEL91_K(t)..) +  BETON ELASTIQUE 
  33 : *       IFLU = 11 : FLUAGE (BPEL91_K(t)..) +  BETON ELASTOPLASTIQUE 
  34 : *       IFLU = 20 : FLUAGE (BPEL91_K(t)=2) +  BETON ELASTIQUE 
  35 : *       IFLU = 21 : FLUAGE (BPEL91_K(t)=2) +  BETON ELASTOPLASTIQUE 
  36 : *
  37 :         TEMPF= TEMP0 + DTEMP
  38 :         TEMPR= TEMPF - TEMP1
  39 :         IF(IFLU.EQ.10.OR.IFLU.EQ.11) THEN
  40 :           DKS  = 1.D0/(1.D0+(20.D0*POAR))  
  41 :           DKE  = 0.4D0
  42 :           DKC  = (120.D0-PHIF)/30.D0+.6667D0*(100.D0-PHIF)/(20.D0+RMOY)
  43 :           DKT1 = 100.D0/(100.D0+TEMP1) 
  44 :           FACT = DKS*(DKE+DKC*DKT1)
  45 :         ELSE 
  46 :           FACT = PHIF
  47 :         ENDIF
  48 : *
  49 :         COEF0 = DSQRT(TEMPR)
  50 :         COEF1 = 5.D0*DSQRT(RMOY)
  51 : *
  52 : *       CALCUL EN DEFORMATION DE FLUAGE TOTALE
  53 : *
  54 :         COEFT = FACT*COEF0/(COEF0+COEF1)
  55 : *
  55 :         CALL BST(D1,SIGR,NSTRS,NSTRS,STFL)
  56 :         DO 2 I=1,NSTRS
  57 :            STFL(I)=COEFT*STFL(I)
  58 :    2    CONTINUE
  59 :         TEMP0 = TEMPF 
  60 :         GOTO 100
  61 : *
  62 :    30 CONTINUE
  63 : *
  64 : *    MODELE DE FLUAGE BETON :LOI1 ----> EPSfl = A1 * T**A2 * Sig/E
  65 : *
  66 : *        IFLU = 31 : FLUAGE (LOI1) +  BETON ELASTOPLASTIQUE
  67 : *        IFLU = 30 : FLUAGE (LOI1) +  BETON ELASTIQUE
  68 : *
  69 :         TEMPF = TEMP0 + DTEMP
  70 :         COEFT = RMOY*(TEMPF**PHIF)
  70 :         CALL BST(D1,SIGR,NSTRS,NSTRS,STFL)
  71 :         DO 3 I=1,NSTRS
  72 :            STFL(I)=COEFT*STFL(I)
  73 :    3    CONTINUE
  74 :         TEMP0 = TEMPF 
  75 :         GOTO 100
  76 : C
  77 :   100 CONTINUE
  78 : C---------------------------------------------------------------------
  79 : 1991  FORMAT(18(1X,E12.5))
  80 :       RETURN
  81 :       END
  82 :  
  83 :