*HEADING
: THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH FRICTIONAL HEAT GENERATION, DEFORMABLE
*RESTART,WRITE,F=2
*NODE,NSET=ALL
   1, 1., 0.
   7, 7., 0.
  11, 1., 1.
  17, 7., 1.
 101, 3., 1.
 105, 5., 1.
 111, 3., 2.
 115, 5., 2.
1001, 1., 0.5
1007, 7., 0.5
1101, 3., 1.5
1105, 5., 1.5
*NGEN     
1,7,1
11,17,1
101,105,1
111,115,1
1001,1007,1
1101,1105,1
*NSET,NSET=BOT,GEN
1,7,1
*NSET,NSET=TOP,GEN
111,115,1
**
*ELEMENT,TYPE=CPE6MT,ELSET=DEAD
1,1,3,13,2,1003,1002
3,5,7,17,6,1007,1006
4,13,11,1,12,1001,1002
6,17,15,5,16,1005,1006
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELEMENT,TYPE=CPE6MT,ELSET=ACTIVE
2,3,5,15,4,1005,1004
5,15,13,3,14,1003,1004
101,101,103,113,102,1103,1102
102,103,105,115,104,1105,1104
103,113,111,101,112,1101,1102
104,115,113,103,114,1103,1104
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10., 
*DENSITY
1000., 
*SPECIFIC HEAT
.001, 
**
*ELSET,ELSET=ETOP
101,102
*ELSET,ELSET=EBOT
4,5,6
*surface,NAME=ASURF
ETOP,S1
*surface, NAME=BSURF
EBOT,S1
*CONTACT PAIR,INTERACTION=ROUGH,SMOOTH=0.05
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*FRICTION,LAGRANGE
0.1, 
** clearance dependent gap conductance not present
*GAP CONDUCTANCE,PRESSURE
20.,100.
*GAP HEAT GENERATION
0.50, 
**
*NSET,NSET=CONST
111,112,114,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
0.0, 0.0,  0.2, 1.0,  0.4, 0.0,  0.6, -1.0
0.8, 0.0,  1.0, 0.0
**
*INITIAL CONDITIONS,TYPE=TEMP
ALL,0.0
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 1
*COUPLED TEMPERATURE,DELTMX=100.
10.,10.
*CONTROLS,PARAMETERS=FIELD,FIELD=TEMPERATURE
5.0E-3, 
*BOUNDARY
BOT,2,2
1,1,1
7,1,1
113,1,1
*CLOAD
113,2,-200.
**
*CONTACT PRINT,SLAVE=ASURF,SUMMARY=YES,TOTALS=NO
*CONTACT FILE,SLAVE=ASURF,F=10
*NODE PRINT
RF,NT11,RFL
*EL PRINT
S,E,TEMP
HFL,
NFLUX, 
*EL FILE,F=10
S,E
*PRINT,CONTACT=YES
*ENERGY PRINT,F=1
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 2 - Generate 120 units of heat, 60 units absorbed
*COUPLED TEMPERATURE,DELTMX=100.
.05,1.,,.05
*BOUNDARY,AMPLITUDE=SLIDE
113, 1,1, 1.5
*CONTACT FILE,SLAVE=ASURF,F=10
*EL FILE,F=10
S,E,TEMP
NFLUX, 
*NODE FILE,F=10
NT,RFL
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 3 - Go to steady state - TEMP=7.5
*COUPLED TEMPERATURE,DELTMX=100.
100.,500.
*CONTACT FILE,SLAVE=ASURF,F=10
*EL FILE,F=10
S,E,TEMP
*NODE FILE,F=10
NT,RFL
*END STEP