*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
*NODE
1003, 1.0,1.1,1.1
1004, 1.0,-0.1,1.1
1005, 1.0,-.1,-0.1
1006, 1.0,1.1,-0.1
*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,
*ELEMENT,ELSET=RIGID,TYPE=R3D4
1001,1003,1004,1005,1006
*RIGID BODY, REF NODE=9999,ELSET=RIGID,ISOTHERMAL=YES
*SURFACE,TYPE=ELEMENT,NAME=RIGHT
12,s3
11,s3
*SURFACE,TYPE=ELEMENT,NAME=WALL
RIGID,SNEG
*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
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=5.0
*ENERGY OUTPUT,VARIABLE=ALL
****************************
*END STEP