Esope source files

   1 : C FONRE1    SOURCE    CHAT      05/01/13    00:06:17     5004           
   2 :       SUBROUTINE FONRE1(RELAX,TZO,BETJEF,BETFLU)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 : C
   7 : C*******************************************************************
   8 : C
   9 : C TZO(MAX)= DATE D'APPLICATION DE LA CHARGE
  10 : C DATCOU= DATE DE COULAGE DU BETON
  11 : C DMT0(NAX)=TEM(NAX,MAX)-TZERO(MAX)
  12 : C E28=4734*SQR(FC)
  13 : C
  14 : C*******************************************************************
  15 : C
  16 : 
  17 : *
  18 :       SEGMENT BETJEF
  19 :        REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  20 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT
  21 :        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  22 :      &                ISIM,IBB,IGAU,IZON
  23 :       ENDSEGMENT
  24 :       SEGMENT BETFLU
  25 :        REAL*8 DATCOU,TP0,E28,TAU1
  26 :        INTEGER NBRC,NCOE,NTPS,NTZERO,IFLU,MC,NC
  27 :       ENDSEGMENT
  28 : *
  29 : *
  30 : *     COMMON /DBETJEF/AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  31 : *    &          TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVIS,ITER,
  32 : *    &                ISIM,IBB,IGAU,IZON
  33 : C
  34 : *     COMMON /DBETFLU/DATCOU,TP0,E28,TAU1,NBRC,NCOE,NTPS,NTZERO,IFLU,
  35 : *    &                MC,NC
  36 : C
  37 : C
  38 : C
  39 : C*******************************************************************
  40 : C              RELAX(NAX,MAX)=FONCTION DE RELAXATION
  41 : C*******************************************************************
  42 : C    
  43 :       DIMENSION RELAX(21,21),DLTRAX(21,21),TZO(21),
  44 :      &          DMT0(21),CR(21,42)           
  45 : C  
  45 :       CALL FONFLU(CR,TZO,BETJEF,BETFLU)
  46 : C
  47 :       MAX=21
  48 :       DO 100 J=1,MAX
  49 : C
  50 :       NAX=21
  51 :       DO 200 I=1,NAX      
  52 : C       
  53 :       IF (I.EQ.1) THEN
  54 :       DLTRAX(I,J)=E28*SQRT(TZO(J)/(TZO(J)*0.85+4))
  55 :       RELAX(I,J)=DLTRAX(I,J)
  56 : C      
  57 :       ELSE IF (I.GT.1) THEN
  58 : C       
  59 :       NAT=I-1
  60 : C
  61 :       SOM=0.D0
  62 :       SOMI=0.D0
  63 : C      
  64 :       DO 300 K=1,NAT
  65 : C
  66 :       JACK1=J+K-1
  67 :       JAC=J+K
  68 :       SOMI=(CR(I,JAC)+CR(I,JACK1)-CR(I-1,JAC)-CR(I-1,JACK1)
  69 :      *)*DLTRAX(K,J)
  70 :       SOM=SOM+SOMI
  71 : C
  72 :   300 CONTINUE
  73 : C
  74 :       RELAX(I,J)=-SOM/(CR(I,J+I)+CR(I,J+I-1))+RELAX(I-1,J)
  75 :       DLTRAX(I,J)=RELAX(I,J)-RELAX(I-1,J)
  76 : C
  77 :       ENDIF
  78 : C      
  79 :   200 CONTINUE
  80 : C
  81 :   100 CONTINUE
  82 : C
  83 : C      DO 500 I=1,NAX
  84 : C
  85 : C      DO 600 J=1,MAX
  86 : C      
  87 : C      IF (RELAX(I,J).LT.0.D0) THEN
  88 : C      RELAX(I,J)=0.D0
  89 : C      ENDIF
  90 : C      
  91 : C  600 CONTINUE
  92 : C
  93 : C  500 CONTINUE
  94 : C      
  95 :       RETURN
  96 :       END
  97 : 
  98 : 
  99 :  
 100 :