C FONRE1    SOURCE    CHAT      05/01/13    00:06:17     5004           
      SUBROUTINE FONRE1(RELAX,TZO,BETJEF,BETFLU)
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
C*******************************************************************
C
C TZO(MAX)= DATE D'APPLICATION DE LA CHARGE
C DATCOU= DATE DE COULAGE DU BETON
C DMT0(NAX)=TEM(NAX,MAX)-TZERO(MAX)
C E28=4734*SQR(FC)
C
C*******************************************************************
C

*
      SEGMENT BETJEF
       REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
     &          TCON,DPSTF1,DPSTF2,TETA,PDT
       INTEGER ICT,ICC,IMOD,IVIS,ITER,
     &                ISIM,IBB,IGAU,IZON
      ENDSEGMENT
      SEGMENT BETFLU
       REAL*8 DATCOU,TP0,E28,TAU1
       INTEGER NBRC,NCOE,NTPS,NTZERO,IFLU,MC,NC
      ENDSEGMENT
*
*
*     COMMON /DBETJEF/AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
*    &          TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVIS,ITER,
*    &                ISIM,IBB,IGAU,IZON
C
*     COMMON /DBETFLU/DATCOU,TP0,E28,TAU1,NBRC,NCOE,NTPS,NTZERO,IFLU,
*    &                MC,NC
C
C
C
C*******************************************************************
C              RELAX(NAX,MAX)=FONCTION DE RELAXATION
C*******************************************************************
C    
      DIMENSION RELAX(21,21),DLTRAX(21,21),TZO(21),
     &          DMT0(21),CR(21,42)           
C  
      CALL FONFLU(CR,TZO,BETJEF,BETFLU)
C
      MAX=21
      DO 100 J=1,MAX
C
      NAX=21
      DO 200 I=1,NAX      
C       
      IF (I.EQ.1) THEN
      DLTRAX(I,J)=E28*SQRT(TZO(J)/(TZO(J)*0.85+4))
      RELAX(I,J)=DLTRAX(I,J)
C      
      ELSE IF (I.GT.1) THEN
C       
      NAT=I-1
C
      SOM=0.D0
      SOMI=0.D0
C      
      DO 300 K=1,NAT
C
      JACK1=J+K-1
      JAC=J+K
      SOMI=(CR(I,JAC)+CR(I,JACK1)-CR(I-1,JAC)-CR(I-1,JACK1)
     *)*DLTRAX(K,J)
      SOM=SOM+SOMI
C
  300 CONTINUE
C
      RELAX(I,J)=-SOM/(CR(I,J+I)+CR(I,J+I-1))+RELAX(I-1,J)
      DLTRAX(I,J)=RELAX(I,J)-RELAX(I-1,J)
C
      ENDIF
C      
  200 CONTINUE
C
  100 CONTINUE
C
C      DO 500 I=1,NAX
C
C      DO 600 J=1,MAX
C      
C      IF (RELAX(I,J).LT.0.D0) THEN
C      RELAX(I,J)=0.D0
C      ENDIF
C      
C  600 CONTINUE
C
C  500 CONTINUE
C      
      RETURN
      END