*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH INTERFACE CONDUCTIVITY, DEFORMABLE
*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.
*NGEN,NSET=ALL
  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=CPE4RPT,ELSET=DEAD
1,1,3,13,11
3,5,7,17,15
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
**########################################
** Material Definitions for elset DEAD
**########################################
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
*PERMEABILITY,SPECIFIC=1.0
 1.E-4,
**
*EXPANSION, PORE FLUID
 0., 
*CONDUCTIVITY, PORE FLUID
10.0, 
*DENSITY, PORE FLUID
1000., 
*SPECIFIC HEAT, PORE FLUID
0.001, 
**#########################################
*ELEMENT,TYPE=CPE4RPT,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
**########################################
** Material Definitions for elset ACTIVE
**########################################
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10., 
*DENSITY
1000., 
*SPECIFIC HEAT
.001, 
*PERMEABILITY,SPECIFIC=1.0
 1.E-4,
*EXPANSION, PORE FLUID
 0., 
*CONDUCTIVITY, PORE FLUID
10.0, 
*DENSITY, PORE FLUID
1000., 
*SPECIFIC HEAT, PORE FLUID
0.001, 
**##################################
** Surface and contact Definitions
**##################################
*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, DEFINITION=TABULAR, 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
*INITIAL CONDITIONS,TYPE=RATIO
ALL, 1.0
**#######################################
*NSET, NSET=BOT_ELE_NODS, ELSET=EBOT
** Node set to specify pore-pressure BC 
** (includes all nodes except contact secondary nodes)
*NSET, NSET=NOD_POR_BC
BOT_ELE_NODS, TOP
**############
**   STEP 1
**############
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 1
*SOILS, CONSOLIDATION, DELTMX=100., UTOL=100.
1.,10.,,1
*CONTROLS,PARAMETERS=FIELD
1.0E-6, 
*BOUNDARY
BOT,2,2
1,1,1
7,1,1
113,1,1
NOD_POR_BC,8,8
*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
*EL FILE,F=10
S,E,TEMP
*NODE FILE,F=10
NT,RFL
*PRINT,CONTACT=YES
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,LE,TEMP
*NODE OUTPUT
NT,RFL,U,POR
*CONTACT OUTPUT
CSTRESS, CDISP
*END STEP
**############
**   STEP 2
**############
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 2
*SOILS, CONSOLIDATION, DELTMX=100., UTOL=100.
100.,1000.,,100.
*BOUNDARY,FIXED
TOP,2,2
*BOUNDARY, AMPLITUDE=SLIDE
113, 1,1, 1.5
*END STEP
**############
**   STEP 3
**############
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 3
*SOILS
100000.,100000.
*BOUNDARY,OP=NEW
BOT,2,2
1,1,1
7,1,1
113,1,1
NOD_POR_BC,8,8
*END STEP
**############
**   STEP 4
**############
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 4
*SOILS
100000.,100000.
*DLOAD
101,P3,200.
102,P3,200.
*BOUNDARY
113,1,1,0.0
*END STEP
**############
**   STEP 5
**############
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 5
*SOILS, CONSOLIDATION, DELTMX=100., UTOL=100.
10.0, 100.0,,10.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
*BOUNDARY,FIXED
111,2,2
113,2,2
115,2,2
*END STEP