*HEADING
EIA3TSSF: THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH FRICTIONAL HEAT GENERATION, DEFORMABLE WITH SHELL ELEMENTS AND SOLID ELEMENTS
*RESTART,WRITE,F=2
*NODE,NSET=ALL
   1, 1.0, 1.0
   7, 1.0, 7.0
  11, 2.0, 1.0
  17, 2.0, 7.0
 101, 2.0, 3.001
 105, 2.0, 4.999
 111, 3.0, 3.001
 115, 3.0, 4.999
1001, 1.5, 1.0
1007, 1.5, 7.0
1101, 2.5, 3.0
1105, 2.5, 5.0
*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=CAX8T,ELSET=DEAD
1,3,1,11,13,2,1001,12,1003
3,7,5,15,17,6,1005,16,1007
*ELEMENT,TYPE=SAX2T,ELSET=SDEAD
1001,13,12,11
1003,17,16,15
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*SHELL SECTION,ELSET=SDEAD,MATERIAL=DEAD
0.2,
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELEMENT,TYPE=CAX8T,ELSET=ACTIVE
2,5,3,13,15,4,1003,14,1005
101,103,101,111,113,102,1101,112,1103
102,105,103,113,115,104,1103,114,1105
*ELEMENT,TYPE=SAX2T,ELSET=SACTIVE
1002,15,14,13
2001,103,102,101
2002,105,104,103
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*SHELL SECTION,ELSET=SACTIVE,MATERIAL=ACTIVE
0.2,
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10., 
*DENSITY
1000., 
*SPECIFIC HEAT
.001, 
**
*ELSET,ELSET=ETOP
2002,2001
*ELSET,ELSET=EBOT
1003,1002,1001
*surface,NAME=ASURF
ETOP,SNEG
*surface,NAME=BSURF
EBOT,SPOS
*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.5, 
**
*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
1.0E-6, 
*BOUNDARY
BOT,1,1
1,2,2
7,2,2
113,2,2
*CLOAD
113,1,-275.
*CONTACT PRINT,SECONDARY=ASURF
*OUTPUT,FIELD,FREQ=10
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=ASURF
*NODE PRINT
RF,NT11,RFL
*EL PRINT
S,E,TEMP
HFL,
NFLUX, 
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,E,TEMP
NFLUX,
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
NT,RFL
*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, 2,2, 1.5
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 3 - Go to steady state - TEMP=0.68
*COUPLED TEMPERATURE,DELTMX=100.
100.,500.
*END STEP