*HEADING
 TEST PRESSURE DEPENDENT GAP CONDUCTANCE (C3D6T)
**
** 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.,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
*NSET,NSET=LEFT
1,101,4,104
*NSET,NSET=QA_TEST
3,6,103,106
*NSET,NSET=NOUT
QA_TEST,
*ELEMENT,TYPE=C3D6T,ELSET=ROD
1,1,5,4,101,105,104
2,1,2,5,101,102,105
3,2,6,5,102,106,105
4,2,3,6,102,103,106
*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
*HEATCAP, ELSET=HCAP
1.0,
*SURFACE, NAME=WALL, TYPE=CYLINDER
1.0,-0.1,0.,2.0,-0.1,0.
1.0,-0.1,-1.
START,0.,0.
LINE,0.,1.1
*SURFACE,TYPE=ELEMENT,NAME=RIGHT
4,S4
*RIGID BODY, REFNODE=9999, ISOTHERMAL=YES, ANALYTICAL SURFACE =WALL
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,2.
*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,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.05
***NODE OUTPUT, NSET=BAR
**NT
***ELEMENT OUTPUT,ELSET=ROD
**S11
**********
*OUTPUT, HISTORY, TIME INTERVAL=2.0
*ENERGY OUTPUT,VARIABLE=ALL
****************************
*END STEP