*Heading
: Thermally Coupled Contact, Finite Sliding Contact With Interface Conductivity, 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.
*Node,Nset=Xtra
 998, 0., 0.
*Ngen
   1,  7,1
  11, 17,1
  101,105,1
 111,115,1
*Nset,Nset=Bot,Gen
 1,7,1
*Nset,Nset=Top,Gen
 111,115,1
**
*Element,Type=cpeg4rt,Elset=Dead
 1,1,3,13,11
 3,5,7,17,15
*Solid Section,Elset=Dead,Material=Dead,Refnode=Xtra
*Material,Name=Dead
*Elastic
 30.E6,.3
*Conductivity
 10.000, 
*Specific Heat
 0.001, 
*Density
 1000., 
**
*Element,Type=cpeg4rt,Elset=Active
 2,3,5,15,13
 101,101,103,113,111
 102,103,105,115,113
*Elset,Elset=Etop
 101,102
*Elset,Elset=Ebot,Generate
 1,3,1
*Solid Section,Elset=Active,Material=Active,Refnode=Xtra
*Material,Name=Active
*Elastic
 30.E6,.3
*Conductivity
 10., 
*Density
 1000., 
*Specific Heat
 .001, 
*Surface, Name=Asurf
 Etop,S1
*Surface, Name=Bsurf
 Ebot,S3
*Contact Pair,Interaction=Rough,Smooth=0.05
 Asurf,Bsurf
*Surface Interaction,Name=Rough
** Clearance Dependent Gap Conductance Not Present
*Gap Conductance,Pressure
 20.,100.
 20.,110.
  5.,190.
  5.,200.
*Gap Heat Generation
 0.0, 
**
*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,2,2
 1,1,1
 7,1,1
 113,1,1
 998,3,5
*Dload
 101,P3,100.
 102,P3,100.
*Dflux
 101,S3,100.
 102,S3,100.
*Film
 2,F1,0.,10.
**
*Contact Print,Secondary=Asurf
*Contact File,Secondary=Asurf,Frequency=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
*Output,Field,Freq=10
*Element Output
 S,LE,Temp,Nflux
*Node Output
 Nt,Rfl,U
*Contact Output
 Cstress, Cdisp
*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,2,2
*Boundary
 998,3,5
*Boundary,Amplitude=Slide
 113, 1,1, 1.5
*Output,Field,Freq=1
*Element Output
 S,LE,Temp,Nflux
*Node Output
 Nt,Rfl,U
*Contact Output
 Cstress, Cdisp
*End Step
**
*Step,Inc=50,Nlgeom,Unsymm=Yes
 Step 3
*Coupled Temperature-Displacement,Steady State
 100000.,100000.
*Boundary,Op=New
 Bot,2,2
 1,1,1
 7,1,1
 113,1,1
 998,3,5
*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,1,1,0.0
 998,3,5
*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, 1,1, 10.0
 998,3,5
*Boundary,Fixed
 111,2,2
 113,2,2
 115,2,2
*End Step