*HEADING
NAFEMS TEST T4, 2-Dimensional Heat Transfer w/ Convection  [C3D8RHT]
**
** The target solution is the temperature at NOUT nodes, T_nout=18.3;
** See ch. 4.3.4 of ABAQUS/Benchmarks Manual
**
*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=C3D8RHT, ELSET=EALL
 1, 1,101,103, 3, 2001,2101,2103,2003
*ELGEN, ELSET=EALL
 1,3,100,1,5,2,10
*SOLID SECTION,MATERIAL=ELMAT, ELSET=EALL
*MATERIAL,NAME=ELMAT
*EXPANSION
0.0, 
*ELASTIC
1.,0.
*DENSITY
7200., 
*SPECIFIC HEAT
440.5, 
*CONDUCTIVITY
 52.0, 
*ELSET,ELSET=SIDE
 3,13,23,33,43
*ELSET,ELSET=TOP
 41,42,43
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
*BOUNDARY
 NALL,1,3
   1,11,,100.0
 101,11,,100.0
 201,11,,100.0
 301,11,,100.0
2001,11,,100.0
2101,11,,100.0
2201,11,,100.0
2301,11,,100.0
*FILM
SIDE,F4,0.0,750.0
TOP ,F5,0.0,750.0
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=QA_TEST
NT, 
*ELEMENT OUTPUT
TEMP, 
*NODE PRINT,NSET=QA_TEST
NT,
*EL PRINT
TEMP,
*NODE FILE,NSET=QA_TEST
NT,
*EL FILE
TEMP,
*END STEP