*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH
INTERFACE CONDUCTIVITY, DEFORMABLE
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*NODE,NSET=ALL
*NGEN, NSET=ALL
*NSET,NSET=QA_TEST
ALL,
*NSET,NSET=BOT,GEN
*NSET,NSET=TOP,GEN
**
*ELEMENT,TYPE=CAX4RT,ELSET=DEAD
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000,
*SPECIFIC HEAT
0.001,
*DENSITY
1000.,
**
*ELEMENT,TYPE=CAX4RT,ELSET=ACTIVE
*ELSET,ELSET=ETOP
*ELSET,ELSET=EBOT
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.,
*DENSITY
1000.,
*SPECIFIC HEAT
.001,
**
*AMPLITUDE, NAME=SLIDE, VALUE=RELATIVE, TIME=STEP, DEFINITION=SMOOTH STEP
0.0, 0.0,  20., 1.0,  40., 0.0,  60., -1.0
80., 0.0,  100., 0.0
*AMPLITUDE, NAME=RAMP, DEFINITION=SMOOTH STEP
0.0,0.0,10.,1.0
**
*INITIAL CONDITIONS,TYPE=TEMP
ALL,20.0
**
*ELSET, ELSET=EALL
ACTIVE, DEAD
**
*SURFACE, NAME=ASURF
ETOP,S1
*SURFACE, NAME=BSURF
EBOT,S3
**
*STEP
STEP 1
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
*CONTACT PAIR,INTERACTION=ROUGH,WEIGHT=0.0
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*GAP CONDUCTANCE,PRESSURE
*GAP HEAT GENERATION
**
*BOUNDARY, TYPE=VELOCITY
BOT,1,1
1,2,2
7,2,2
113,2,2
*DLOAD, AMPLITUDE=RAMP
101,P3,100.
102,P3,100.
*DFLUX, AMPLITUDE=RAMP
101,S3,100.
102,S3,100.
*FILM, AMPLITUDE=RAMP
2,F1,20.,10.
*********************
*FILE OUTPUT,N=1
*NODE FILE,NSET=ALL
NT,
********************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
NT,
*NODE OUTPUT,NSET=ALL
U,
************
*CONTACT OUTPUT
CSTRESS,
*ELEMENT OUTPUT
S,
********************
***OUTPUT,HISTORY, FREQUENCY=500
***NODE OUTPUT,NSET=ALL
**U,NT
*******************
*OUTPUT,HISTORY, FREQUENCY=99999
*ENERGY OUTPUT,VAR=ALL
******************
*END STEP
**
*STEP
STEP 2
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
*BOUNDARY,TYPE=VELOCITY
TOP,1,1
*BOUNDARY,AMPLITUDE=SLIDE
113,2,2,1.5
*END STEP
**
*STEP
STEP 3 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*BOUNDARY,OP=NEW, TYPE=VELOCITY
BOT,1,1
1,2,2
7,2,2
113,2,2
*END STEP
**
*STEP
STEP 4 - RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE,NAME=RAMP2, DEFINITION=SMOOTH STEP
0.0,0.5,10.0,1.0
*DLOAD, AMPLITUDE = RAMP2
101,P3,200.
102,P3,200.
*END STEP
**
*STEP
STEP 5
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE, NAME=REMOVE, VALUE=RELATIVE, TIME=STEP
0.0, 1.0,  30.0, 0.0,  40.0, 0.0
*DFLUX, AMP=REMOVE
101, S3, 100.0
102, S3, 100.0
*BOUNDARY, TYPE=VELOCITY
113, 2,2,0.1
111,1,1
113,1,1
115,1,1
*END STEP