*Heading
: Thermally Coupled Contact, Finite Sliding With Interace Radiation, 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 
*Node,Nset=Refnd
 9998,0.,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=Cpeg6mt,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,Refnode=Refnd
*Material,Name=Dead
*Elastic
 30.E6,.3
*Conductivity
 10.000, 
*Specific Heat
 0.001, 
*Density
 1000., 
**
*Element,Type=Cpeg6mt,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,Refnode=Refnd
*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
*Gap Radiation
 .74074,.74074
 1.,0.
 1.,1.
*Gap Heat Generation
 0.0, 
*Physical Constants, Absolute Zero=-273.16,Stefanb=1.7e-6
**
*Amplitude, Name=Slide, Value=Relative, Time=Step
 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
*Boundary
 Refnd,3,5
**
*Step,Inc=50,Nlgeom,Amp=Ramp,Unsymm=Yes
Step 1
*Coupled Temperature,Deltmx=100.
 1.,10.
*Solution Technique,Type=Separated
*Controls,Parameter=Field,Field=Temperature
 5.0e-3, 
*Boundary
Bot,2,2
 1,1,1
 7,1,1
 113,1,1
*Dload
 103,P1,100.
 104,P1,100.
*Dflux
 103,S1,100.
 104,S1,100.
*Film
 2,F1,0.,10.
**
*Contact Print,Secondary=Asurf,Summary=No,Totals=Yes
*Contact File,Secondary=Asurf,F=10
*Node Print
 Rf,
 Nt11,
 Rfl,
*El Print
 S,
 E,
 Temp,
 Hfl,
 Nflux, 
*El File,F=10
 S,E,Temp
 Nflux, 
*Node File,F=10
 Nt,Rfl
*Print,Contact=Yes
*Output, Field,Freq=1
*Node Output
 Rf,U,Nt,Rfl
*Element Output
 S,E
*Contact Output
 Cdisp, Cstress
*End Step
**
*Step,Inc=50,Nlgeom,Amp=Ramp,Unsymm=Yes
 Step 2
*Coupled Temperature,Deltmx=100.
 10.,1000.,,100.
*Boundary,Fixed
 Top,2,2
*Boundary,Amplitude=Slide
 113, 1,1, 1.5
*End Step
**
*Step,Inc=50,Nlgeom,Unsymm=Yes
 Step 3
*Coupled Temperature-Displacement,Steady State
 10000.,10000.
*Boundary,Op=New
 Bot,2,2
 1,1,1
 7,1,1
 113,1,1
 Refnd,3,5
*End Step
**
*Step,Inc=50,Nlgeom,Unsymm=Yes
 Step 4
*Coupled Temperature-Displacement,Steady State
 10000.,10000.
*Dload
 103,P1,200.
 104,P1,200.
*Boundary
 113,1,1,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,1,1,10.0
*Boundary,Fixed
 111,2,2
 112,2,2
 113,2,2
 114,2,2
 115,2,2
*End Step