*HEADING
: THERMALLY COUPLED CONTACT, FINTE SLIDING CONTACT WITH INTERACE RADIATION, DEFORMABLE
*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=CGAX6MT,ELSET=DEAD
1,3,1,11,2,1001,1002
3,7,5,15,6,1005,1006
4,11,13,3,12,1003,1002
6,15,17,7,16,1007,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=CGAX6MT,ELSET=ACTIVE
2,5,3,13,4,1003,1004
5,13,15,5,14,1005,1004
101,103,101,111,102,1101,1102
103,111,113,103,112,1103,1102
102,105,103,113,104,1103,1104
104,113,115,105,114,1105,1104
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10., 
*DENSITY
1000., 
*SPECIFIC HEAT
.001, 
**
*ELSET,ELSET=ETOP
102,101
*ELSET,ELSET=EBOT
4,5,6
*surface,NAME=ASURF
ETOP,S1
*surface,NAME=BSURF
EBOT,S1
*CONTACT PAIR,INTERACTION=SMOOTH,SMOOTH=0.05
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*GAP RADIATION
.74074, .74074
1.,0.
1.,1.
*GAP HEAT GENERATION
0.0, 
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16, STEFAN BOLTZMAN=1.7E-6 
**
*AMPLITUDE, NAME=SLIDE, VALUE=RELATIVE, TIME=STEP
0.0, 0.0, 200., 1.0, 400., 0.0, 600., -1.0
800., 0.0,  1000., 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
BOT,5,5
1,2,2
7,2,2
113,2,2
113,5,5
*DLOAD
103,P1,100.
104,P1,100.
*DFLUX
103,S1,100.
104,S1,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,U
*PRINT,CONTACT=YES
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 2
*COUPLED TEMPERATURE,DELTMX=100.
10.,1000.,,100.
*BOUNDARY,FIXED
TOP,1,1
*BOUNDARY,AMPLITUDE=SLIDE
113, 2,2, 1.0607
113, 5,5, 0.5304
*END STEP
**
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 3
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
100000.,100000.
*BOUNDARY,OP=NEW
BOT,1,1
BOT,5,5
1,2,2
7,2,2
113,2,2
113,5,5
*END STEP
**
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 4
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
100000.,100000.
*DLOAD
103,P1,200.
104,P1,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
103, S1, 100.0
104, S1, 100.0
*BOUNDARY
113,2,2,10.0
*BOUNDARY,FIXED
111,1,1
113,1,1
115,1,1
*END STEP