*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT 
WITH FRICTIONAL HEAT GENERATION, DEFORMABLE;PENALTY CONTACT
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*NODE,NSET=ALL
1,1.,1.
7,1.,7.
11,2.,1.
17,2.,7.
101,2.,3.
105,2.,5.
111,3.,3.
115,3.,5.
*NGEN, NSET=ALL
1,7,1
11,17,1
101,105,1
111,115,1
*NSET,NSET=QA_TEST
ALL,
*NSET,NSET=BOT,GEN
1,7,1
*NSET,NSET=TOP,GEN
111,115,1
**
*ELEMENT,TYPE=CAX4RT,ELSET=DEAD
1,3,1,11,13
3,7,5,15,17
*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
2,5,3,13,15
101,103,101,111,113
102,105,103,113,115
*ELSET,ELSET=ETOP
102,101
*ELSET,ELSET=EBOT
3,2,1
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.
*DENSITY
1000.
*SPECIFIC HEAT
.001
**
*NSET,NSET=CONST
111,115
*EQUATION
2
CONST,1,1.0,113,1,-1.0
2
CONST,2,1.0,113,2,-1.0
**
*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
4e-3,10.
*CONTACT PAIR,INTERACTION=ROUGH,WEIGHT=0.0,MECHANICAL CONSTRAINT=PENALTY,
CPSET=C1
ASURF,BSURF
*CONTACT CONTROLS,CPSET=C1,SCALE PENALTY=0.05
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION
0.1
*GAP CONDUCTANCE,PRESSURE
20.,0.
20.,100.
*GAP HEAT GENERATION
0.50,
**
*BOUNDARY, TYPE=VELOCITY
BOT,1,1
1,2,2
7,2,2
113,2,2
*DLOAD,AMP=RAMP
ETOP,P3,8.7532
*********************
*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 - GENERATE 120 UNITS OF HEAT, 60 UNITS ABSORBED
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
2.5e-3,100.
*BOUNDARY,AMPLITUDE=SLIDE, OP=MOD
113,2,2,1.5
*END STEP
**
*STEP
STEP 3 - GO TO STEADY STATE - SHOULD GO TO 20.68
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
3e-3,10.0
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E4
*END STEP