ETL 1110-2-365
31 Aug 94
ANNEX 3: HEAT GENERATION SUBROUTINES
A3-1. User Subroutines
As discussed in the text, there are two user subroutines available for modeling heat generation in the concrete
due to the heat of hydration. These two subroutines are DFLUX and HETVAL. A listing of each subroutine is
provided with comments discussing various portions of the program. Comments which are for the purpose of
this appendix only will be placed in double quotations. Please note that these subroutines were developed for
use with version 4.9 of ABAQUS and may need to be modified for later ABAQUS versions.
*** USER SUBROUTINE DFLUX ***
SUBROUTINE DFLUX(FLUX,TEMP,KSTEP,KINC,TIME,NOEL,NPT,COORDS,
and JLTYP)
IMPLICIT REAL*8 (A-H,O-Z)
C**************************************************************************************
C
C VERSION 2.0
C THE ADIABATIC CURVE IN THIS VERSION OF DFLUX IS BASED ON THE
C ORIGINAL CURVE USED FOR L&D26. UNITS IN THE T ARRAY ARE HOURS.
C UNITS IN THE HEAT ARRAY ARE BTU/(LB-IN**3)
C
C NQ IS THE NUMBER OF POINTS IN ARRAYS T & Q. ENTIME IS THE ENDTIME
C FOR DFLUX. STTIME GIVES THE START TIMES FOR ELEMENTS IN HOURS.
C THE DIMENSION OF STTIME MUST BE AS LARGE AS THE NUMBER OF ELEMENTS.
C YOU MUST CHANGE THE VALUES IN STTIME TO CONFORM TO YOUR PROBLEM.
C FOR INSTANCE, IF THE FIRST POUR IS MODELED USING 50 ELEMENTS,
C 50*0.0 WOULD START DFLUX AT TIME 0 FOR THE FIRST 50 ELEMENTS.
C
C**************************************************************************************
"The array COORDS is simply for the coordinates, Q and T are for the arrays given
below and PROP is an array for the density and specific heat as given on the
DATA PROP card"
DIMENSION COORDS(3),Q(20),T(20),PROP(2)
"STTIME is defined below."
COMMON /ELDEF/ STTIME(736)
DATA PROP/.08681,.21/
DATA ENTIME/648.1/
DATA NQ/20/
"Array T is the time associated with each heat flux given
in array Q"
DATA T/.25,.5,.75,1.,1.25,1.5,1.75,2.,
$
2.5,3.,3.5,4.,6.,7.,8.,
$
9.,10.0,13.,15.,27./
A3-1
Previous Page
