*HEADING
          TWO DIMENSIONAL HEAT TRANSFER INTERFACE  (CAX4RT)
          GAP CONDUCTANCE
	  ANALYTICAL 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  
**
*NODE,NSET=NALL
1,10.,0.0
2,11.0,0.0
3,11.0,.5
4,10.0,0.5
7,11.,1.0
8,10.,1.0
*NSET, NSET=SIDE3
7,8
*NSET,NSET=QA_TEST
1,2
*NSET,NSET=NOUT
QA_TEST,
*NODE, NSET=BOT
13, 10., -0.5, 0
*ELEMENT,TYPE=CAX4RT,ELSET=MID
1,1,2,3,4
*ELEMENT,TYPE=CAX4RT,ELSET=TOP
2,4,3,7,8
*ELEMENT, TYPE=HEATCAP, ELSET=HEATC
90001, 13
*HEATCAP, ELSET=HEATC
1.0
*ELSET, ELSET=SOLID
1,2
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*CONDUCTIVITY
1.
*SPECIFIC HEAT
46.
*DENSITY
7850.
*ELASTIC
1.0E-2,0.0
*SURFACE, NAME=SURF3
MID,S1
*SURFACE, NAME=SURF4, TYPE=S
START, 10.0,-.1
LINE, 11.0,-.1
*RIGID BODY, ANALYTICAL=SURF4, REFNODE=13, ISOTHERMAL=YES
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,300000.
*CONTACT PAIR, INTERACTION=GAP,WEIGHT=0.0,MECHANICAL CONSTRAINT=PENALTY
SURF3,SURF4
*SURFACE INTERACTION, NAME=GAP
*GAP CONDUCTANCE
10.,0.
10.,10.
*BOUNDARY,OP=NEW
13, 1, 6, 0.0
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=1E4
***NODE OUTPUT, NSET=NOUT
**NT, 
***ELEMENT OUTPUT, ELSET=SOLID
**HFL2,
********
*OUTPUT, HISTORY,TIME INTERV=3.0E5
*ENERGY OUTPUT,VAR=ALL
********
*ENDSTEP