C RKFCCA SOURCE AM 12/02/28 00:28:54 7289 FUNCTION RKFCCA(t,tmax) c Reduction of compression strength in calcareous concrete c according to Table 3.1 of prEN1992-1-2 IMPLICIT REAL*8 (A-B,D-H,O-Z) implicit integer (I-K,M,N) implicit logical (L) implicit character*10 (C) i1 = 1 r1 = 1.0d0 r0 = 0.0d0 precision = 1.0d-3 if (tmax.le.100.d0) then fc = 1.00d0 else if (tmax.le.200.d0) then fc = 1.00d0-.0003d0*(tmax-100.d0) else if (tmax.le.300.d0) then fc = 0.97d0-.0006d0*(tmax-200.d0) else if (tmax.le.400.d0) then fc = 0.91d0-.0006d0*(tmax-300.d0) else if (tmax.le.500.d0) then fc = 0.85d0-.0011d0*(tmax-400.d0) else if (tmax.le.600.d0) then fc = 0.74d0-.0014d0*(tmax-500.d0) else if (tmax.le.700.d0) then fc = 0.60d0-.0017d0*(tmax-600.d0) else if (tmax.le.800.d0) then fc = 0.43d0-.0016d0*(tmax-700.d0) else if (tmax.le.900.d0) then fc = 0.27d0-.0012d0*(tmax-800.d0) else if (tmax.le.1000.d0) then fc = 0.15d0-.0009d0*(tmax-900.d0) else if (tmax.le.1100.d0) then fc = 0.06d0-.0004d0*(tmax-1000.d0) else fc = max(0.02d0-.0002d0*(tmax-1100.d0),0.0001D0) endif c Attention to the cases when T<20°C and Tmax=20°C if ( (t.lt.tmax) .and. (dabs(tmax-20.).ge.precision) ) then c Cooling down c calculate fc(20°C) = rmul * fc(tmax) if (tmax.lt.300.) then rmul = 1.0-0.1*(t-20.)/280. else rmul = 0.9 endif fc20 = rmul*fc c calculate fc(t) = interpolation between fc(20) and fc(tmax) fc = fc20+(fc-fc20)*(t-20.)/(tmax-20.) endif rkfcca = fc return end