Esope source files

   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 :