Esope source files

   1 : C EPSTHSI   SOURCE    AM        12/02/28    00:28:16     7289           
   2 :         FUNCTION EPSTHSI(t,tmax)
   3 : 
   4 : c       Thermal elongation of siliceous concrete 
   5 : c     according to equations A.1 and A.2 of ENV1992-1-2
   6 : c                same as     3.3.1 (1) of prEN 1992-1-2
   7 : 
   8 :       IMPLICIT REAL*8 (A-B,D-H,O-Z)
   9 :       implicit integer (I-K,M,N)
  10 :       implicit logical (L)
  11 :       implicit character*10 (C)
  12 :       
  13 :       precision = 1.0d-3
  14 :       r1 = 1.0d0
  15 :       r0 = 0.d0
  16 : 
  17 :       if (tmax.le.700.) then
  18 :         epsth = -1.80184d-4+tmax*(9.d-6+2.3d-11*tmax*tmax)
  19 :       else
  20 :         epsth = 14.d-3
  21 :       endif
  22 : 
  23 : c     Attention to the cases where T<20°C and Tmax=20°C
  24 :       if ( (t.lt.tmax) .and. (dabs(tmax-20.).ge.precision) ) then
  25 :         if (tmax.le.300.) then
  26 :           epsres = -.58d-3*(tmax-20.)/280.
  27 :         else if (tmax.le.400.) then
  28 :           epsres = -.58d-3+0.29d-3*(tmax-300.)/100.
  29 :         else if (tmax.le.600.) then
  30 :           epsres = -.29d-3+2.00d-3*(tmax-400.)/200.
  31 :         else if (tmax.le.800.) then
  32 :           epsres = 1.71d-3+1.58d-3*(tmax-600.)/200.
  33 :         else if (tmax.le.900.) then
  34 :           epsres = 3.29d-3+1.71d-3*(tmax-800.)/100.
  35 :         else
  36 :           epsres = 5.00d-3
  37 :         endif
  38 :         epsrevers = epsth-epsres
  39 :         dtmax = tmax-20.
  40 :         dt = t-20.
  41 :         dtsdtmax = dt/dtmax
  42 :         epsth = epsres+epsrevers*dtsdtmax*(.6+.4*dtsdtmax)
  43 :       endif
  44 : 
  45 :       epsthsi = epsth
  46 : 
  47 :       return
  48 :       end
  49 : 
  50 : 
  51 : 
  52 : 
  53 : 
  54 :