*HEADING
          TWO DIMENSIONAL HEAT TRANSFER INTERFACE  (C3D4T)
          GAP RADIATION
**
** The temperature at the NOUT nodes, based on a 1D model is:
** NT_out=100+HFL2*R_th; R_th=L/(kA);  
** See ch. 1.7.1 of ABAQUS/Verification Manual  
**
*PHYSICAL CONSTANTS,ABS=-460.,STEFAN=5.E-10
*NODE,NSET=NALL
      1,          10.,          -1.,           0.
      2,          11.,          -1.,           0.
      3,          11.,          -1.,           1.
      4,          10.,          -1.,           1.
      5,          11.,         -0.1,           0.
      6,          11.,         -0.1,           1.
      7,          10.,         -0.1,           0.
      8,          10.,         -0.1,           1.
      9,          11.,           0.,           0.
     10,          11.,           1.,           0.
     11,          11.,           1.,           1.
     12,          11.,           0.,           1.
     13,          10.,           1.,           0.
     14,          10.,           1.,           1.
     15,          10.,           0.,           0.
     16,          10.,           0.,           1.
*NODE
     17,      10.4309,        -0.55,          0.5
     18,      10.5001,          0.5,          0.5
*ELEMENT, TYPE=C3D4T
 1, 17,  4,  3,  8
 2, 17,  1,  3,  4
 3, 17,  8,  3,  6
 4, 17,  5,  8,  6
 5, 17,  7,  8,  5
 6, 17,  5,  1,  7
 7, 17,  5,  6,  2
 8, 17,  2,  1,  5
 9, 17,  2,  6,  3
10, 17,  2,  3,  1
11, 17,  7,  4,  8
12, 17,  1,  4,  7
*ELEMENT, TYPE=C3D4T,ELSET=TOP
13, 18, 13, 14, 10
14, 18, 10, 15, 13
15, 18, 10, 14, 11
16, 18, 14, 12, 11
17, 18, 16, 12, 14
18, 18, 15, 12, 16
19, 18,  9, 11, 12
20, 18,  9, 12, 15
21, 18, 10, 11,  9
22, 18,  9, 15, 10
23, 18, 13, 16, 14
24, 18, 15, 16, 13
*NSET, NSET=SIDE3
10,11,13,14
*NSET, NSET=BOT
1,2,3,4,5,6,7,8
*NSET,NSET=QA_TEST
9,12,15,16
*NSET,NSET=NOUT
QA_TEST,
*ELSET, ELSET=SOLID, GENERATE
1,24,1
*ELSET, ELSET=SOLID
1,2,3,4
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*CONDUCTIVITY
1., 
*SPECIFIC HEAT
46., 
*DENSITY
7850., 
*ELASTIC
1.0,0.0
*SURFACE,TYPE=ELEMENT, NAME=SURF3
18,S3
20,S3
*SURFACE,TYPE=ELEMENT, NAME=SURF4
4,S3
5,S3
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,1000000.
*CONTACT PAIR, INTERACTION=GAP,WEIGHT=1.0
SURF3,SURF4
*SURFACE INTERACTION, NAME=GAP
*GAP RADIATION
1.0,1.0
1.0,0.0
1.0,10.0
*BOUNDARY,OP=NEW
SIDE3,11,,100.
BOT,11,,50.
******************************
*FILE OUTPUT, NUM=1
*NODE FILE,NSET=NOUT
 NT,
******************************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
NT
*NODE OUTPUT,NSET=NALL
NT
********************************
***OUTPUT, HISTORY,TIME INTERV=1E5
***NODE OUTPUT, NSET=NOUT
**NT, 
***ELEMENT OUTPUT, ELSET=TOP
**HFL2,
********
*OUTPUT, HISTORY,TIME INTERV=1.0E6
*ENERGY OUTPUT,VAR=ALL
*********
*ENDSTEP