1 : C RKFCSI SOURCE AM 12/02/28 00:28:55 7289 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. (dabs(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 :
© Cast3M 2003 - All rights reserved.
Disclaimer