ETL 1110-2-365
31 Aug 94
C 99 FORMAT(3G15.6,4I8)
RETURN
END
*** USER SUBROUTINE HETVAL ***
SUBROUTINE HETVAL (CMNAME,TEMP,TIME,DTIME,SVAR,FLUX,PREDEF)
IMPLICIT REAL (A-H,O-Z)
DIMENSION SVAR(1),PREDEF(1)
CHARACTER*8 CMNAME
C******************************************************************
C
C Calculate and return Volumetric Heating Rate at each integration
C point for each concrete element by material name R James 10/27/92
C This method is independent of element numbers - Each lift is
C given a material name using *HEAT GENERATION material option
C
C NQ1 is no. of points in T1 & Q1 arrays
C Start time for each lift based on 5 day placement increments
C NOTE - FLUXES IN Q1 ARRAY OBTAINED FROM WES
C
- ASSUMED UNITS ARE BTU/(hr-in**3)
C
- BASED ON SPEC. WT = 0.08449, SPEC.HEAT = .22 & HOURS
C
C******************************************************************
C
PARAMETER (NQ1=21)
DIMENSION Q1(NQ1),T1(NQ1)
SAVE Q1,T1
C
"T1 contains the time (in hours) associated with each value of
flux given in Q1."
DATA T1 /
6.00, 12.00, 18.00, 24.00, 30.00,
$
36.00, 48.00, 60.00, 84.00, 108.00,
$
120.00, 144.00, 168.00, 192.00, 216.00,
$
240.00, 264.00, 288.00, 336.00, 360.00,
$
648.00 /
DATA Q1 /0.00846090,0.01001042,0.01284686,0.00521171,0.00539755,
$
0.00567338,0.00520116,0.00414147,0.00338928,0.00279045,
$
0.00224723,0.00194696,0.00160559,0.00138042,0.00117012,
$
0.00091299,0.00071519,0.00058468,0.00037927,0.00030453,
$
0.00015165 /
C
FLUX = 0.0
C
C
Find lift number of current material
C Assume name is of form xxLIFTii where ii is lift number
C
L = INDEX(CMNAME,'LIFT')
IF (L .EQ. 0) RETURN
READ (CMNAME(L+4:L+5),100) LIFT
A3-3
Previous Page
