*HEADING
          TWO DIMENSIONAL HEAT TRANSFER INTERFACE  (C3D6T)
          GAP RADIATION
	  DISCRETE RIGID SURFACE
          PENALTY CONTACT
**
** 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
      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.,           1.
     16,          10.,           0.,           0.
     17,          10.,          -.01,           0.
     18,          11.,          -.01,           0.
     19,          10.,          -.01,           1.
     20,          11.,          -.01,           1.
*NSET, NSET=SIDE3
10,11,13,14
*NSET,NSET=QA_TEST
9,12,15,16
*NSET,NSET=NOUT
QA_TEST,
*NODE, NSET=BOT
1,10.5,-.5
*ELEMENT, TYPE=C3D6T
3,15,12,9,14,11,10
4,16,15,9,13,14,10
*ELSET, ELSET=SOLID
3,4
*ELEMENT, TYPE=R3D4, ELSET=RIGID
90002,17,18,20,19
*ELEMENT,TYPE=HEATCAP,ELSET=HEATC
99001, 1
*HEATCAP,ELSET=HEATC
1.0,
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*CONDUCTIVITY
1., 
*SPECIFIC HEAT
46., 
*DENSITY
7850., 
*ELASTIC
1.0E-2,0.0
*SURFACE,TYPE=ELEMENT, NAME=SURF3
4,S1
3,S1
*SURFACE,TYPE=ELEMENT,NAME=SURF4
RIGID, SPOS
*RIGID BODY,REFNODE=1,ELSET=RIGID,ISOTHERMAL=YES
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,1000000.
*CONTACT PAIR, INTERACTION=GAP,WEIGHT=0.0,MECHANICAL CONSTRAINT=PENALTY
SURF3,SURF4
*SURFACE INTERACTION, NAME=GAP
*GAP RADIATION
1.0,1.0
1.0,0.0
1.0,10.0
*BOUNDARY,OP=NEW
1, 1, 6, 0.0
SIDE3,11,,100.
BOT,11,,50.
1, ENCASTRE
******************************
*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=SOLID
**HFL2,
********
*OUTPUT, HISTORY,TIME INTERV=1.0E6
*ENERGY OUTPUT,VAR=ALL
*********
*ENDSTEP