Télécharger epsc1co.eso

Retour à la liste

Numérotation des lignes :

epsc1co
  1. C EPSC1CO SOURCE AM 12/02/28 00:28:15 7289
  2. FUNCTION EPSC1CO(tmax,rloc)
  3.  
  4. c Strain to peak stress in concrete
  5. c according to Table A.2 of ENV1992-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. dimension R(10,4)
  13.  
  14. i1 = 1
  15. i2 = 2
  16. i3 = 3
  17. i4 = 4
  18. i5 = 5
  19. i6 = 6
  20. r0 = 0.0d0
  21. r1 = 1.0d0
  22. precision = 1.0d-3
  23.  
  24. data R
  25. c temperatures
  26. * /20., 100., 200., 300., 400., 500., 600., 700., 800., 900.,
  27. c minimum values
  28. * 2.5, 2.5, 3.0, 4.0, 4.5, 5.5, 6.5, 7.5, 8.5, 10.,
  29. c recommended values
  30. * 2.5, 3.5, 4.5, 6.0, 7.5, 9.5, 12.5, 14.0, 14.5, 15.,
  31. c maximum values
  32. * 2.5, 4.0, 5.5, 7.0, 10.0, 15.0, 25.0, 25.0, 25.0, 25./
  33.  
  34. if (rloc.lt.r0) then
  35. jloc = i2
  36. else if (rloc.eq.r0) then
  37. jloc = i3
  38. else
  39. jloc = i4
  40. endif
  41.  
  42. if (tmax.le.20.) then
  43. epsc1 = R(i1,jloc)
  44. else if (tmax.ge.900.) then
  45. epsc1 = R(10,jloc)
  46. else
  47. iloc = i1
  48. do while (tmax.gt.R(iloc,i1))
  49. iloc = iloc+i1
  50. enddo
  51. epsc1 = R(iloc-i1,jloc) + (R(iloc,jloc)-R(iloc-1,jloc))
  52. * *(Tmax-R(iloc-i1,i1))/(R(iloc,i1)-R(iloc-i1,i1))
  53. endif
  54.  
  55. epsc1 = epsc1*1.d-3
  56. EPSC1CO = epsc1
  57.  
  58. return
  59. end
  60.  
  61.  
  62.  

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