*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH FRICTIONAL HEAT GENERATION, 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, TYPE=SURFACE TO SURFACE
ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*FRICTION
0.1, 
** clearance dependent gap conductance not present
*GAP CONDUCTANCE,PRESSURE
20.,100.
*GAP HEAT GENERATION
.50, 
**
*NSET,NSET=CONST
111,115
*EQUATION
2, 
CONST,1,1.0,113,1,-1.0
2, 
CONST,2,1.0,113,2,-1.0
**
*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
*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
**
**############
*BOUNDARY, FIXED
NOD_POR_BC, 8,8,
**
**############
**   STEP 1
**############
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 1
*SOILS, CONSOLIDATION, DELTMX=100., UTOL=100.
10.,10.
*CONTROLS,PARAMETERS=FIELD,FIELD=TEMPERATURE
5.0E-3, 
*BOUNDARY
BOT,2,2
1,1,1
7,1,1
113,1,1
*CLOAD
113,2,-200.
**
*CONTACT PRINT, SECONDARY=ASURF
*CONTACT FILE, SECONDARY=ASURF,F=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
*ENERGY PRINT,F=1
*ENERGY FILE,F=1
*OUTPUT, FIELD, FREQUENCY=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
Step 2 - Generate 120 units of heat, 60 units absorbed
*SOILS, CONSOLIDATION, DELTMX=100., UTOL=100.
.05,1.,,.05
*BOUNDARY,AMPLITUDE=SLIDE
113,1,1,1.5
*END STEP
**############
**   STEP 3
**############
*STEP,INC=50,NLGEOM,AMP=RAMP
Step 3 - Go to steady state - should go to Temp=7.5
*SOILS, CONSOLIDATION, DELTMX=100., UTOL=100.
100.,500.,,100.
*END STEP