Télécharger epsthsi.eso

Retour à la liste

Numérotation des lignes :

epsthsi
  1. C EPSTHSI SOURCE CB215821 16/04/21 21:16:39 8920
  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. (ABS(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.  
  55.  
  56.  
  57.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales