Sources Esope | Cast3M

bpel
C BPEL      SOURCE    PV        07/11/23    21:15:32     5978
       SUBROUTINE BPEL(RM,RH,FC28,T0,T,FLU)
*
*  Calcul de la fonction de fluage selon l'eurocode 2
*
*
* entree
*        RM   rayon moyen de la piece [cm]
*        RH  pourcentage d'humidite (0.&lt;roh&lt;100.)
*        FC28  resistance moyenne en compression  a 28 jours [MPa]            type de ciment
*        T0   temps de chargement [j]
*        T    temps actuel [j]
*
* sortie
*        FLU  J(to,t) [Pa-1]
*
*
      IMPLICIT INTEGER(I-N)
       IMPLICIT REAL*8(A-H,O-Z)
       real*8 j
*
*
       J=T0
       FCJ=J/(4.76D0+0.83D0*J)*FC28
       IF (FCJ.GT.(1.1D0*FC28)) THEN
         FCJ=1.1D0*FC28
       ENDIF
       IF (FCJ.GT.40.D0) THEN
         FCJ=J/(1.4D0+0.95D0*J)*FC28
       ENDIF
       IF (FCJ.GT.(1.1D0*FC28)) THEN
         FCJ=1.1D0*FC28
       ENDIF
       ECT0=11000.D0*FCJ**(1.D0/3.D0)
       EC28=11000.D0*FC28**(1.D0/3.D0)
       AKT0=100.D0/(100.D0+T0)
       AKC=(120.D0-RH)/30.D0+2.D0*(100.D0-RH)/
     &                   (3.D0*(20.D0+RM))
       AKFL=0.4D0+AKT0*AKC
       FTT0=SQRT(T-T0)/(SQRT(T-T0)+5.D0*SQRT(RM))
       PHI=AKFL*FTT0
       FLU=(1.D0/ECT0+PHI/EC28)/1.D6
*       write(6,*) ' rm rh fc28 t t0 ', rm,th,fc28,t ,t0
*       write(6,*) ' flu ' , flu
       END