*HEADING
          TWO DIMENSIONAL HEAT TRANSFER INTERFACE  (C3D4T)
          GAP CONDUCTANCE
	  DISCRETE 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.7.1 of ABAQUS/Verification Manual    
**
*NODE,NSET=NALL
      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.
     18,      10.5001,          0.5,          0.5
*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
*ELEMENT, TYPE=R3D4, ELSET=RIGID
90002, 7,5,6,8
*NSET, NSET=SIDE3
10,11,13,14
*NODE, NSET=BOT
1, 10.5,-0.5, 0.
*NSET,NSET=QA_TEST
9,12,15,16
*NSET,NSET=NOUT
QA_TEST,
*ELEMENT,TYPE=HEATCAP,ELSET=HEATC
99001, 1
*HEATCAP,ELSET=HEATC
1.0,
*ELSET, ELSET=SOLID, GENERATE
13,24,1
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*CONDUCTIVITY
1.
*SPECIFIC HEAT
46.
*DENSITY
7850.
*ELASTIC
1.0E-2,0.0
*RIGID BODY, ELSET=RIGID, REFNODE=1, ISOTHERMAL=YES
*SURFACE, NAME=SURF3
18,S3
20,S3
*SURFACE, NAME=SURF4
90002, SPOS
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,200000.
*CONTACT PAIR, INTERACTION=GAP,WEIGHT=0.0
SURF3,SURF4
*SURFACE INTERACTION, NAME=GAP
*GAP CONDUCTANCE
10.,0.
10.,10.
*BOUNDARY,OP=NEW
1, 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=TOP
**HFL2,
********
*OUTPUT, HISTORY,TIME INTERV=2.0E5
*ENERGY OUTPUT,VAR=ALL
********
*ENDSTEP