*HEADING
 TEST PRESSURE DEPENDENT GAP CONDUCTANCE (CPS4RT) 
**
** The target solution is the temp and pressure at the right side
** of the rod: T_nout=122.6; p=s11=161.3
** See ch. 1.6.3 of ABAQUS/Benchmarks Manual
**   
*NODE,NSET=BAR
1,0.
2,.5
3,1.
101,0.,1.
102,0.5,1.
103,1.0,1.
*NODE
1001,1.0,-0.1
1002,1.0,1.1
*NSET,NSET=LEFT
1,101
*NSET,NSET=QA_TEST
3,103
*NSET,NSET=NOUT
QA_TEST,
*ELEMENT,TYPE=CPS4RT,ELSET=ROD
1,1,2,102,101
2,2,3,103,102
*SOLID SECTION,ELSET=ROD,MATERIAL=A1
1., 
*MATERIAL,NAME=A1
*ELASTIC
1.E5, 
*EXPANSION
1.E-5, 
*CONDUCTIVITY
10., 
*DENSITY
1., 
*SPECIFIC HEAT
1., 
*NODE, NSET=WALL
9999,2.,0.5
*ELEMENT, TYPE=HEATCAP, ELSET=HCAP
9999,9999
*ELEMENT,ELSET=RIGID,TYPE=R2D2
1001,1001,1002
*RIGID BODY, REFNODE=9999,ELSET=RIGID,ISOTHERMAL=YES
*HEATCAP, ELSET=HCAP
1.0,
*SURFACE,TYPE=ELEMENT,NAME=RIGHT
2,S2
*SURFACE,TYPE=ELEMENT,NAME=WALL
RIGID,SPOS
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 , 5.
*CONTACT PAIR,INTERACTION=GAP,WEIGHT=0.0
RIGHT,WALL
*SURFACE INTERACTION,NAME=GAP
*GAP CONDUCTANCE, PRESSURE
2.,0.
1002.,5000.
*BOUNDARY,TYPE=VELOCITY
BAR,2
LEFT,1
WALL,1,2
WALL,6,6
*BOUNDARY
LEFT,11,11,200.
WALL,11,11,100.
***************************
*FILE OUTPUT,NUM=1
*NODE FILE,NSET=NOUT
NT 
***************************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT
NT
*ELEMENT OUTPUT
S
*NODE OUTPUT,NSET=QA_TEST
NT
****************************
***OUTPUT,HISTORY,TIME INTERVAL=0.1
***NODE OUTPUT, NSET=BAR
**NT
***ELEMENT OUTPUT,ELSET=ROD
**S11
**********
*OUTPUT, HISTORY, TIME INTERVAL=5.0
*ENERGY OUTPUT,VARIABLE=ALL
****************************
*END STEP