*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH 
FRICTIONAL HEAT GENERATION, DEFORMABLE [CPS6MT]
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*NODE,NSET=ALL
  1, 1., 0.
  7, 7., 0.
***
 11, 1., 0.5
 17, 7., 0.5
**
 21, 1., 1.
 27, 7., 1.
**
101, 3., 1.
105, 5., 1.
**
111, 3., 1.5
115, 5., 1.5
**
121, 3., 2.
125, 5., 2.
*NGEN,NSET=ALL
  1,  7,1
 11, 17,1
 21, 27,1
101,105,1
111,115,1
121,125,1
*NSET,NSET=QA_TEST
ALL
*NSET,NSET=BOT,GEN
1,7,1
*NSET,NSET=TOP,GEN
121,125,1
**
*ELEMENT,TYPE=CPS6MT,ELSET=DEAD
1,1,3,21,2,12,11
2,3,23,21,13,22,12
5,5,7,25,6,16,15
6,7,27,25,17,26,16 
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELEMENT,TYPE=CPS6MT,ELSET=ACTIVE
3,3,5,23,4,14,13
4,5,25,23,15,24,14
***
101,101,103,121,102,112,111
102,103,123,121,113,122,112
103,103,105,123,104,114,113
104,105,125,123,115,124,114
*ELSET,ELSET=ETOP
101,103
*ELSET,ELSET=EBOT, GENERATE
2,6,2
*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
*NSET,NSET=CONST
121,122,124,125
*EQUATION
2, 
CONST,1,1.0,123,1,-1.0
2
CONST,2,1.0,123,2,-1.0
**
**
*SURFACE, NAME=ASURF
ETOP,S1
*SURFACE, NAME=BSURF
EBOT,S2
**
*STEP
STEP 1
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
1.25e-3,10.
**
*CONTACT PAIR,INTERACTION=ROUGH,WEIGHT=0.0
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION
0.1,
*GAP CONDUCTANCE,PRESSURE
20.,  0.
20.,100.
20.,110.
 5.,190.
 5.,200.
*GAP HEAT GENERATION
0.5, 
**
*BOUNDARY, TYPE=VELOCITY,OP=NEW
BOT,2,2
1,1,1
7,1,1
123,1,1
*CLOAD, AMP=RAMP
123,2,-200.
*********************
*FILE OUTPUT,N=1
*NODE FILE,NSET=ALL
NT,
********************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
NT,
*NODE OUTPUT,NSET=ALL
U,
*ELEMENT OUTPUT
TEMP,
************
*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
1.25e-3,100.
*BOUNDARY,AMPLITUDE=SLIDE,OP=MOD
123,1,1,1.5
*END STEP
**
*STEP
STEADY STATE --TEMP SHOULD GO TO 27.5
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
2.5e-3,10.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E4
*END STEP