*HEADING
EIA3TSSR: THERMALLY COUPLED CONTACT, FINTE SLIDING CONTACT WITH INTERACE RADIATION, 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=SMOOTH,SMOOTH=0.05
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*GAP RADIATION
.74074, .74074
1.,0.
1.,1.
** converted from old data: 1.E-6, 1.E-6
*GAP HEAT GENERATION
0.0, 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, STEFAN BOLTZMANN=1.7E-6
**
*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.
1.,10.
*CONTROLS,PARAMETERS=FIELD
1.0E-6, 
*BOUNDARY
BOT,1,1
1,2,2
7,2,2
113,2,2
*DLOAD
101,P3,100.
102,P3,100.
*DFLUX
101,S3,100.
102,S3,100.
*FILM
2,F1,0.,10.
*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
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 2
*COUPLED TEMPERATURE,DELTMX=100.
100.,1000.,,100.
*BOUNDARY,FIXED
TOP,1,1
*BOUNDARY,AMPLITUDE=SLIDE
113 ,2,2, 1.5
*END STEP
**
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 3
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
100000.,100000.
*BOUNDARY,OP=NEW
BOT,1,1
1,2,2
7,2,2
113,2,2
*END STEP
**
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 4
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
100000.,100000.
*DLOAD
101,P3,200.
102,P3,200.
*BOUNDARY
113,2,2,0.0
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 5
*COUPLED TEMPERATURE-DISPLACEMENT,DELTMX=100.
10.0, 100.0
*AMPLITUDE, NAME=REMOVE, VALUE=RELATIVE, TIME=STEP
0.0, 1.0,  30.0, 0.0,  100.0, 0.0
*DFLUX, AMP=REMOVE
101, S3, 100.0
102, S3, 100.0
*BOUNDARY
113,2,2,10.0
*BOUNDARY,FIXED
111,1,1
113,1,1
115,1,1
*END STEP