*HEADING
 TEST CLEARANCE DEPENDENT GAP CONDUCTANCE (C3D4T) PENALTY CONTACT 
**
** The target solution is the temp. and displ. at NOUT nodes:
** T_nout=285.4; U1_nout=3.427e-3
** See ch. 1.6.3 of ABAQUS/Benchmarks Manual
**    
*NODE,NSET=BAR
      1,           1.,           0.,           0.
      2,           1.,           1.,           0.
      3,           1.,           1.,           1.
      4,           1.,           0.,           1.
      5,           0.,           1.,           0.
      6,           0.,           1.,           1.
      7,           0.,           0.,           0.
      8,           0.,           0.,           1.
      9,     0.499657,          0.5,          0.5
*NSET,NSET=BAR1,GENERATE
1,8,1
*NSET,NSET=LEFT,GENERATE
5,8,1
*NSET,NSET=QA_TEST
1,2,3,4
*NSET,NSET=NOUT
QA_TEST
*ELEMENT,TYPE=C3D4T,ELSET=ROD
 1, 9, 6, 5, 8
 2, 9, 5, 7, 8
 3, 9, 7, 1, 4
 4, 9, 6, 4, 3
 5, 9, 8, 7, 4
 6, 9, 8, 4, 6
 7, 9, 3, 2, 6
 8, 9, 2, 5, 6
 9, 9, 5, 2, 7
10, 9, 2, 1, 7
11, 9, 4, 1, 3
12, 9, 1, 2, 3
*SOLID SECTION,ELSET=ROD,MATERIAL=A1
1., 
*MATERIAL,NAME=A1
*ELASTIC
1.E5, 
*EXPANSION
1.E-5, 
*CONDUCTIVITY
1., 
*DENSITY
1., 
*SPECIFIC HEAT
1., 
*NODE, NSET=WALL
9999,2.,0.5
*ELEMENT, TYPE=HEATCAP, ELSET=HCAP
9999,9999
*HEATCAP, ELSET=HCAP
1.0,
*SURFACE, NAME=WALL, TYPE=CYLINDER
1.01,-0.1,0,2.01,-0.1,0
1.01,-0.1,-1
START,0,0
LINE,0,1.1
*SURFACE,TYPE=ELEMENT,NAME=RIGHT
12,S3
11,S3
*RIGID BODY, REFNODE=9999, ISOTHERMAL=YES, ANALYTICAL SURFACE =WALL
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,3.
*CONTACT PAIR,INTERACTION=GAP,WEIGHT=0.0,MECHANICAL CONSTRAINT=PENALTY
RIGHT,WALL
*SURFACE INTERACTION,NAME=GAP
*GAP CONDUCTANCE
2.,0.
0.,.02
*BOUNDARY,TYPE=VELOCITY
BAR1,2
LEFT,1
WALL,1,6,0.0
*BOUNDARY
LEFT,11,11,400.
WALL,11,11,200.
*********************
*FILE OUTPUT,NUM=1
*NODE FILE,NSET=NOUT
U,NT 
********************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT
NT
*NODE OUTPUT,NSET=QA_TEST
U,NT
********************
***OUTPUT,HISTORY,TIME INTERVAL=0.05
***NODE OUTPUT, NSET=BAR
**U1,NT
***********
*OUTPUT, HISTORY, TIME INTERVAL=3.0
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP