*HEADING
NAFEMS TEST T4, 2-Dimensional Heat Transfer w/ Convection  [EC3D8RT]
**
** The target solution is the temperature at NOUT nodes, T_nout=18.3;
** See ch. 4.3.4 of ABAQUS/Benchmarks Manual
**
** Multimaterial Eulerian HT is tested by using 3 identical material
**  segments with different volume fractions that add up to 1.0
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE
 1,0.0,0.0
 301,0.6,0.0
 11,0.0, 1.0
 311,0.6,1.0
*NGEN,NSET=EDGE1
 1,11
*NGEN,NSET=EDGE2
 301,311
*NFILL,NSET=BACK
 EDGE1,EDGE2,6,50
*NODE
 2001,0.0,0.0,1.0
 2301,0.6,0.0,1.0
 2011,0.0, 1.0,1.0
 2311,0.6,1.0,1.0
*NGEN,NSET=EDGE3
 2001,2011
*NGEN,NSET=EDGE4
 2301,2311
*NFILL,NSET=FRONT
 EDGE3,EDGE4,6,50
*NFILL,NSET=NALL
BACK,FRONT,2,1000
*NSET,NSET=NOUT
303,2303
*NSET,NSET=QA_TEST
303,2303
*ELEMENT,TYPE=EC3D8RT, ELSET=EALL
 1, 1,51, 52, 2, 2001,2051,2052,2002
*ELGEN, ELSET=EALL
 1,6,50,1,10,1,10
*EULERIAN SECTION, ELSET=EALL
ELMAT1, ELMAT_PART1
ELMAT2, ELMAT_PART2
ELMAT3, ELMAT_PART3
*MATERIAL,NAME=ELMAT1
*EXPANSION
0.0, 
*ELASTIC
1.,0.
*DENSITY
7200., 
*SPECIFIC HEAT
440.5, 
*CONDUCTIVITY
 52.0, 
*MATERIAL,NAME=ELMAT2
*EXPANSION
0.0, 
*ELASTIC
1.,0.
*DENSITY
7200., 
*SPECIFIC HEAT
440.5, 
*CONDUCTIVITY
 52.0, 
*MATERIAL,NAME=ELMAT3
*EXPANSION
0.0, 
*ELASTIC
1.,0.
*DENSITY
7200., 
*SPECIFIC HEAT
440.5, 
*CONDUCTIVITY
 52.0, 
*ELSET,ELSET=SIDE
 6,16,26,36,46,56,66,76,86,96
*ELSET,ELSET=TOP
 91,92,93,94,95,96
*Initial Conditions, type=volume fraction
EALL, ELMAT_PART1, 0.1
EALL, ELMAT_PART2, 0.4
EALL, ELMAT_PART3, 0.5
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
 ,20000.,
*BOUNDARY
   1,11,,100.0
  51,11,,100.0
 101,11,,100.0
 151,11,,100.0
 201,11,,100.0
 251,11,,100.0
 301,11,,100.0
2001,11,,100.0
2051,11,,100.0
2101,11,,100.0
2151,11,,100.0
2201,11,,100.0
2251,11,,100.0
2301,11,,100.0
*FILM
SIDE,F4,0.0,750.0
TOP ,F5,0.0,750.0
***************************
*OUTPUT,FIELD,NUMBER INTERVAL=1
*NODE OUTPUT,NSET=QA_TEST
NT, 
*OUTPUT,FIELD,NUMBER INTERVAL=20
*NODE OUTPUT
NT, 
*ELEMENT OUTPUT
TEMP,HFL 
***************************
*OUTPUT,HISTORY,time interval=1000
*NODE OUTPUT,NSET=NOUT
NT,
*ELEMENT OUTPUT, ELSET=TOP
TEMP,HFL
*ENERGY OUTPUT, VAR=ALL
***************************
*END STEP