*HEADING
 TEST PRESSURE DEPENDENT GAP CONDUCTANCE (C3D4T)   
**
** 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,           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.5,          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
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
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
RIGHT,WALL
*SURFACE INTERACTION,NAME=GAP
*GAP CONDUCTANCE, PRESSURE
2.,0.
1002.,5000.
*BOUNDARY,TYPE=VELOCITY
BAR1,2
LEFT,1
WALL,1,6,0.0
*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=3.0
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP