*HEADING
: THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH FRICTIONAL HEAT GENERATION, DEFORMABLE
*RESTART,WRITE,F=2
*NODE,NSET=ALL
*NGEN
*NSET,NSET=BOT,GEN
*NSET,NSET=TOP,GEN
**
*ELEMENT,TYPE=CPE4RHT,ELSET=DEAD
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000,
*SPECIFIC HEAT
0.001,
*DENSITY
1000.,
**
*ELEMENT,TYPE=CPE4RHT,ELSET=ACTIVE
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.,
*DENSITY
1000.,
*SPECIFIC HEAT
.001,
**
*ELSET,ELSET=ETOP
*ELSET,ELSET=EBOT,GENERATE
*surface, NAME=ASURF
ETOP,S1
*surface, NAME=BSURF
EBOT,S3
*CONTACT PAIR, INTERACTION=ROUGH,SMOOTH=0.3
ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*FRICTION,LAGRANGE
0.1,
** clearance dependent gap conductance not present
*GAP CONDUCTANCE,PRESSURE
*GAP HEAT GENERATION
**
*NSET,NSET=CONST
*EQUATION
CONST,1,1.0,113,1,-1.0
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
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 1
*COUPLED TEMPERATURE,DELTMX=100.
*SOLUTION TECHNIQUE,TYPE=SEPARATED
*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, SLAVE=ASURF
*CONTACT FILE, SLAVE=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
*ENERGY PRINT,F=1
*ENERGY FILE,F=1
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP
Step 2 - Generate 120 units of heat, 60 units absorbed
*COUPLED TEMPERATURE,DELTMX=100.
*BOUNDARY,AMPLITUDE=SLIDE
113,1,1,1.5
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP
Step 3 - Go to steady state - should go to Temp=7.5
*COUPLED TEMPERATURE,DELTMX=100.
*END STEP