*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT 
WITH INTERFACE CONDUCTIVITY, DEFORMABLE; PENALTY CONTACT
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*NODE,NSET=ALL
  1, 1., 0.,-0.01
  7, 7., 0.,-0.01
 11, 1., 1.,-0.01
 17, 7., 1.,-0.01
 21, 1., 0., 1.01
 27, 7., 0., 1.01
 31, 1., 1., 1.01
 37, 7., 1., 1.01
101, 3., 1.
105, 5., 1.
111, 3., 2.
115, 5., 2.
121, 3., 1., 1.0
125, 5., 1., 1.0
131, 3., 2., 1.0
135, 5., 2., 1.0
*NGEN, NSET=ALL
  1,  7
 11, 17
 21, 27
 31, 37
101,105
111,115
121,125
131,135
*NSET,NSET=QA_TEST
ALL,
*NSET,NSET=BOT,GEN
1,7
21,27
*NSET,NSET=TOP,GEN
111,115
131,135
**
*ELEMENT, TYPE=C3D4T
1,33,13,23,31
2,3,23,13,1
3,11,31,1,13
4,21,1,31,23
101,133,113,123,131
102,103,123,113,101
103,111,131,101,113
104,121,101,131,123
**
*ELEMENT,TYPE=C3D4T,ELSET=DEAD5
25,5,27,17,35
35,3,25,15,33
45,1,23,13,31
*ELEMENT,TYPE=C3D4T,ELSET=ACTIVE5
125,101,123,113,131
135,103,125,115,133
**
*ELGEN, ELSET=DEAD1
1,3,2,4
*ELGEN, ELSET=DEAD2
2,3,2,4
*ELGEN, ELSET=DEAD3
3,3,2,4
*ELGEN, ELSET=DEAD4
4,3,2,4
*ELSET, ELSET=DEAD
DEAD1,DEAD2,DEAD3,DEAD4,DEAD5
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELGEN, ELSET=ACTIVE1
101,2,2,4
*ELGEN, ELSET=ACTIVE2
102,2,2,4
*ELGEN, ELSET=ACTIVE3
103,2,2,4
*ELGEN, ELSET=ACTIVE4
104,2,2,4
*ELSET, ELSET=ACTIVE
ACTIVE1,ACTIVE2,ACTIVE3,ACTIVE4,ACTIVE5
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELSET,ELSET=ETOP, GENERATE
102,108,2
*ELSET,ELSET=EBOT, GENERATE
1,11,2
*ELSET, ELSET=ELOAD, GENERATE
101,107,2
*ELSET, ELSET=EFILM
6,8
**
*AMPLITUDE, NAME=SLIDE, VALUE=RELATIVE, TIME=STEP, DEFINITION=SMOOTH STEP
0.0, 0.0,  20.0, 1.0,  40., 0.0,  60., -1.0
80., 0.0,  100.0, 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,S2
*SURFACE, NAME=BSURF
EBOT,S2
**
*STEP
STEP 1
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
1.75e-3,10.
*CONTACT PAIR,INTERACTION=ROUGH,WEIGHT=0.0,MECHANICAL CONSTRAINT=PENALTY
ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*GAP CONDUCTANCE,PRESSURE
20.,0.
20.,100.
20.,110.
 5.,190.
 5.,200.
*GAP HEAT GENERATION
0.0, 
**
*BOUNDARY, TYPE=VELOCITY
BOT,2,3
TOP,3,3
1,1,1
7,1,1
21,1,1
27,1,1
113,1,1
133,1,1
*DLOAD, AMP=RAMP
ELOAD,P2,100.
*DFLUX, AMP=RAMP
ELOAD,S2, 100.
*FILM, AMP=RAMP
EFILM,F2,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
1.75e-3,100.
*BOUNDARY,TYPE=VELOCITY
TOP,2,2
*BOUNDARY,AMPLITUDE=SLIDE
113,1,1, 1.5
133,1,1, 1.5
*END STEP
**
*STEP
STEP 3 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
3e-3,5.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*BOUNDARY,OP=NEW,TYPE=VELOCITY
BOT,2,3
TOP,3,3
1,1,1
7,1,1
21,1,1
27,1,1
113,1,1
133,1,1
*END STEP
**
*STEP
STEP 4 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
2.0e-3,10.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE,NAME=RAMP2
0.0,0.5,10.0,1.0
*DLOAD, AMP=RAMP2
ELOAD, P2,200.
*END STEP
**
*STEP
STEP 5
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
2.0e-3,40.0
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE, NAME=REMOVE, VALUE=RELATIVE, TIME=STEP, 
DEFINITION=SMOOTH STEP
0.0, 1.0,  30.0, 0.0,  40.0, 0.0
*DFLUX, AMP=REMOVE
ELOAD, S2, 100.
*BOUNDARY, TYPE=VELOCITY
113, 1,1, 0.1
133, 1,1, 0.1
*BOUNDARY,FIXED
TOP,2,2
*END STEP