Télécharger rkfcsi.eso

Retour à la liste

Numérotation des lignes :

  1. C RKFCSI SOURCE CB215821 16/04/21 21:18:19 8920
  2. FUNCTION RKFCSI(t,tmax)
  3.  
  4. c Reduction of compression strength in siliceous concrete
  5. c according to Table 3.1 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. fc = r1
  19. else if (tmax.le.200.d0) then
  20. fc = 1.00d0-.0005d0*(tmax-100.d0)
  21. else if (tmax.le.300.d0) then
  22. fc = 0.95d0-.0010d0*(tmax-200.d0)
  23. else if (tmax.le.400.d0) then
  24. fc = 0.85d0-.0010d0*(tmax-300.d0)
  25. else if (tmax.le.500.d0) then
  26. fc = 0.75d0-.0015d0*(tmax-400.d0)
  27. else if (tmax.le.600.d0) then
  28. fc = 0.60d0-.0015d0*(tmax-500.d0)
  29. else if (tmax.le.700.d0) then
  30. fc = 0.45d0-.0015d0*(tmax-600.d0)
  31. else if (tmax.le.800.d0) then
  32. fc = 0.30d0-.0015d0*(tmax-700.d0)
  33. else if (tmax.le.900.d0) then
  34. fc = 0.15d0-.0007d0*(tmax-800.d0)
  35. else if (tmax.le.1000.d0) then
  36. fc = 0.08d0-.0004d0*(tmax-900.d0)
  37. else if (tmax.le.1100.d0) then
  38. fc = 0.04d0-.0003d0*(tmax-1000.d0)
  39. else
  40. fc = max(0.01d0-.0001d0*(tmax-1100.d0),0.0001D0)
  41. endif
  42.  
  43. c Attention to the cases where T<20 C et Tmax=20 C
  44. if ( (t.lt.tmax) .and. (ABS(tmax-20).ge.precision) ) then
  45. c Cooling down
  46. c calculate fc(20øc) = rmul * fc(tmax)
  47. if (tmax.lt.300.) then
  48. rmul = 1.-0.1*(t-20.)/280.
  49. else
  50. rmul = 0.9
  51. endif
  52. fc20 = rmul*fc
  53. c calculate fc(t) = interpolation between fc(20) and fc(tmax)
  54. fc = fc20+(fc-fc20)*(t-20.)/(tmax-20.)
  55. endif
  56.  
  57. rkfcsi = fc
  58.  
  59. return
  60. end
  61.  
  62.  
  63.  
  64.  
  65.  

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