*HEADING
: THERMALLY COUPLED CONTACT -- FRICTIONAL HEAT GENERATION
*RESTART,WRITE,F=2
*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
  1,  7
 11, 17
 21, 27
 31, 37
101,105
111,115
121,125
131,135
*NSET,NSET=BOT,GEN
1,7
21,27
*NSET,NSET=TOP,GEN
111,115
131,135
**
*ELEMENT,TYPE=C3D8RHT,ELSET=DEAD
1,1,3,13,11,21,23,33,31
2,3,5,15,13,23,25,35,33
3,5,7,17,15,25,27,37,35
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELEMENT,TYPE=C3D8RHT,ELSET=ACTIVE
101,101,103,113,111,121,123,133,131
102,103,105,115,113,123,125,135,133
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
*surface, NAME=SLAVE
ACTIVE, S3
*surface, NAME=MASTER
DEAD, S5
*CONTACT PAIR, INTERACTION=SLIDE
SLAVE, MASTER
*SURFACE INTERACTION, NAME=SLIDE
*FRICTION,LAGRANGE
0.1, 
** clearance dependent gap conductance not present
*GAP CONDUCTANCE,PRESSURE
20.,100.
*GAP HEAT GENERATION
0.5, 
**
*NSET,NSET=CONST
111,115,131,133,135
*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
*AMPLITUDE, NAME=REMOVE, VALUE=RELATIVE, TIME=STEP
0.0, 1.0,  30.0, 0.0,  100.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,3
TOP,3,3
1,1,1
7,1,1
21,1,1
27,1,1
113,1,1
*CLOAD
113,2,-200.
*OUTPUT, FIELD, FREQ=10
*NODE OUTPUT
 RF,NT,RFL
*ELEMENT OUTPUT
 S,E,TEMP
 NFLUX,
*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
*CONTROLS,PARAMETER=FIELD,FIELD=DISPLACEMENT
,,,,,,,1.E-7
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 3 - Go to steady state - should go to Temp=7.5
*COUPLED TEMPERATURE,DELTMX=100.
100.,500.
*CONTROLS,PARAMETER=FIELD,FIELD=DISPLACEMENT
,,,,,,,1.E-7
*END STEP