*HEADING
 TEST CLEARANCE  DEPENDENT GAP CONDUCTANCE (C3D8RT) 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, 0.,0.,0.
2, .5,0.,0.
3, 1.,0.,0.
4, 0.,0.,-1
5, 0.5,0.,-1
6, 1.,0.,-1
101, 0.,1.,0.
102, 0.5,1.,0.
103, 1.0,1.,0.
104, 0.,1.,-1
105, 0.5,1.,-1
106, 1.0,1.,-1
*NODE
1003, 1.01,1.1,0.1
1004, 1.01,-0.1,0.1
1005, 1.01,-.1,-1.1
1006, 1.01,1.1,-1.1
*NSET,NSET=LEFT
1,101,4,104
*NSET,NSET=QA_TEST
3,6,103,106
*NSET,NSET=NOUT
QA_TEST,
*ELEMENT,TYPE=C3D8RT,ELSET=ROD
1,1,2,5,4,101,102,105,104
2,2,3,6,5,102,103,106,105
*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
*ELEMENT,ELSET=RIGID,TYPE=R3D4
1001,1003,1004,1005,1006
*RIGID BODY, REF NODE=9999,ELSET=RIGID,ISOTHERMAL=YES
*HEATCAP, ELSET=HCAP
1.0,
*SURFACE,TYPE=ELEMENT,NAME=RIGHT
2,S4
*SURFACE,TYPE=ELEMENT,NAME=WALL
RIGID,SNEG
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,5.
*CONTACT PAIR,INTERACTION=GAP,WEIGHT=0.0,MECHANICAL CONSTRAINT=PENALTY,
CPSET=C1
RIGHT,WALL
*CONTACT CONTROLS,CPSET=C1,SCALE PENALTY=0.5
*SURFACE INTERACTION,NAME=GAP
*GAP CONDUCTANCE
2.,0.
0.,.02
*BOUNDARY,TYPE=VELOCITY
BAR,2
LEFT,1
WALL,1,6
*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.25
***NODE OUTPUT, NSET=BAR
**U1,NT
*************************
*OUTPUT, HISTORY, TIME INTERVAL=5.0
*ENERGY OUTPUT,VARIABLE=ALL
*************************
*END STEP