Esope source files

   1 : C RKFTCO    SOURCE    AM        12/02/28    00:28:56     7289           
   2 :         FUNCTION rkftco(t,tmax)
   3 : 
   4 : c       Reduction of tension strength in siliceous concrete 
   5 : c     according to 3.2.2.2 (3) of prEN1992-1-2
   6 : 
   7 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
   8 :       implicit integer (I-K,M,N)
   9 :       implicit logical (L)
  10 :       implicit character*10 (C)
  11 :       
  12 :       i1 = 1
  13 :       r1 = 1.0d0
  14 :       r0 = 0.0d0
  15 :       precision = 1.0d-3
  16 : 
  17 :       if (tmax.le.100.d0) then
  18 :         ft = 1.0d0
  19 :       else if (tmax.le.600.d0) then
  20 :         ft = 1.00d0-(tmax-100.d0)/500.d0
  21 :       else
  22 :         ft = 0.0d0
  23 :       endif               
  24 : 
  25 : c     Attention to the cases where T<20 C et Tmax=20 C
  26 :       if ( (t.lt.tmax) .and. (dabs(tmax-20).ge.precision) ) then
  27 : c       Cooling down
  28 : c       calculate  ft(20øc) = rmul * ft(tmax)
  29 :         if (tmax.lt.300.) then
  30 :           rmul = r1-0.1*(tmax-20.)/280.
  31 :         else
  32 :           rmul = 0.9
  33 :         endif
  34 :         ft20 = rmul*ft
  35 : c       calculate fc(t) = interpolation between fc(20) and fc(tmax)
  36 :         ft = ft20+(ft-ft20)*(t-20.)/(tmax-20)
  37 :       endif
  38 : 
  39 :         rkftco = ft
  40 : 
  41 :         return
  42 :         end
  43 :         
  44 :