*HEADING
 THERMAL PROPERTIES: EXPLICIT [C3D8RT]
 GAP CLEARANCE-DEPENDENT CONDUCTIVITY
**
** See ch. 2.3.1 and 1.7.1 of ABAQUS/Verification Manual
**
*PREPRINT,MODEL=YES,HISTORY=YES
*NODE,NSET=PLANEA
1, 
3,,,1.
5,,1.,1.
7,,1.
*NCOPY,CHANGE NUMBER=10,OLD SET=PLANEA,SHIFT,NEW SET=PLANEB
1.,0.,0.

*NCOPY,CHANGE NUMBER=20,OLD SET=PLANEA,SHIFT,NEW SET=PLANEC
1.01,0.,0.

*NCOPY,CHANGE NUMBER=30,OLD SET=PLANEA,SHIFT,NEW SET=PLANED
2.01,0.,0.

*NSET,NSET=ALL
PLANEA,PLANEB,PLANEC,PLANED
*NSET,NSET=CDPLANE
PLANEC,PLANED
*NSET,NSET=BCPLANE
PLANEB,PLANEC
*NSET,NSET=QA_TEST
PLANEB
*ELEMENT,TYPE=c3d8rt,ELSET=ONE
1,11,13,15,17,1,3,5,7
*ELEMENT,TYPE=c3d8rt,ELSET=TWO
2,31,33,35,37,21,23,25,27
*ELSET,ELSET=EALL
ONE,TWO
*SOLID SECTION,ELSET=EALL,MATERIAL=BLOCK
*MATERIAL,NAME=BLOCK
*ELASTIC,TYPE=ISO
1.E5, 
*EXPANSION,TYPE=ISO,ZERO=0.0
1.E-5, 
*CONDUCTIVITY,TYPE=ISO
1., 
*DENSITY
1000.,
*SPECIFIC HEAT
100.,
*BOUNDARY
ALL,2,3
PLANEA,1
PLANEC,1
*SURFACE,TYPE=ELEMENT,NAME=ASURF
ONE,S1
*SURFACE,TYPE=ELEMENT,NAME=BSURF
TWO,S2
*AMPLITUDE,NAME=RAMP, DEFINITION=SMOOTH
0,0,2e5,1
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 , 5.e5
*FIXED MASS SCALING,ELSET=EALL,FACTOR=1.E6
*CONTACT PAIR,INTERACTION=GAPC,weight=0.0
ASURF,BSURF
*SURFACE INTERACTION,NAME=GAPC
*GAP CONDUCTANCE
2.,0.
0.,.02
*BOUNDARY,AMP=RAMP
PLANEA,11,11,400.
CDPLANE,11,11,200.
****************************
*FILE OUTPUT, NUM=1
*NODE FILE, NSET=PLANEB
U, NT, 
****************************
*OUTPUT, FIELD, NUMBER=1
*NODE OUTPUT, NSET=QA_TEST
U,NT,
********
*NODE OUTPUT
NT,U
*ELEMENT OUTPUT
S
*****************************
***OUTPUT, HISTORY, TIME INTERV=2.5E4
***NODE OUTPUT, NSET=PLANEB
**U1,NT, 
*****
*OUTPUT, HISTORY, TIME INTERV=5.E5
*ENERGY OUTPUT,VAR=ALL
*****************************
*END STEP