*HEADING
          TWO DIMENSIONAL HEAT TRANSFER INTERFACE  (CPE6MT)
          GAP RADIATION
	  RIGID BODY-ANALYTICAL RIGID SURFACE
**
** 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.11.6 of ABAQUS/
** Verification Manual 
**
*PHYSICAL CONSTANTS,ABS=-460.,STEFAN=5.E-10
*NODE,NSET=NALL
1,10.,0.0
101,10.5,0.0
201,11.0,0.0
5,10.,1.0
105,10.5,1.0
205,11.0,1.0
1001,10.5,-.5
*NGEN,NSET=NALL
1,5
101,105
201,205
*NSET, NSET=SIDE3
5,205,105
*NSET,NSET=QA_TEST
1,201,101
*NSET,NSET=NOUT
QA_TEST,
*ELEMENT,TYPE=CPE6MT
1,1,201,203,101,202,102
2,1,203,3,102,103,2
3,3,203,205,103,204,104
4,3,205,5,104,105,4
*ELSET,ELSET=BOT
1
*ELSET,ELSET=TOP
4
*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.0e-2,0.0
*RIGIDBODY,ISOTHERMAL=NO,REFNODE=1001,ELSET=SOLID
*NODE,NSET=BOT
11,10.5,-.5
*ELEMENT,TYPE=HEATCAP,ELSET=HEATC
99001, 11
*HEATCAP,ELSET=HEATC
1.0,
*SURFACE,TYPE=ELEMENT, NAME=SURF3
BOT,S1
*SURFACE,NAME=SURF4,TYPE=S
START, 10.0,-.01
LINE, 11.0,-.01
*RIGIDBODY,ISOTHERMAL=YES,REFNODE=11, ANALYTICAL SURFACE =SURF4
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT,DIRECT USER CONTROL
1000,1000000.
*CONTACT PAIR, INTERACTION=GAP,PENALTY, WEIGHT=0.0
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
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=5E4
*NODE OUTPUT, NSET=NOUT
NT, 
*ELEMENT OUTPUT, ELSET=SOLID
HFL2,
********
*OUTPUT, HISTORY,TIME INTERV=1E6
*ENERGY OUTPUT,VAR=ALL
********
*ENDSTEP