*HEADING
          TWO DIMENSIONAL HEAT TRANSFER INTERFACE (CAX3T)
          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
      1,          10.,        -0.01
      4,          11.,        -0.01
      5,          10.,           0.
      6,          11.,           0.
      7,          11.,           1.
      8,          10.,           1.
*NSET, NSET=SIDE3
7,8
*NSET,NSET=QA_TEST
5,6
*NSET,NSET=NOUT
QA_TEST,
*NODE,NSET=GAP1
11,10.5, -0.5
*ELEMENT,TYPE=CAX3T
3, 5, 6, 8
4, 8, 6, 7
*ELSET,ELSET=SOLID
3,4
*ELEMENT, TYPE=RAX2, ELSET=RIGID
90002, 1, 4
*ELEMENT, TYPE=HEATCAP, ELSET=HEATC
90001, 11
*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
3,S1
*SURFACE,TYPE=ELEMENT, NAME=SURF4
90002, SPOS
*RIGID BODY, ELSET=RIGID, REFNODE=11, 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
11, 1, 6, 0.0
SIDE3,11,,100.
GAP1,11,,50.
11,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
*********
*END STEP