*HEADING
: THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH INTERFACE CONDUCTIVITY, DEFORMABLE
*RESTART, WRITE, F=2
*NODE,NSET=ALL
*NGEN
*NSET,NSET=BOT,GEN
*NSET,NSET=TOP,GEN
**
*ELEMENT,TYPE=CAX4T,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=CAX4T,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
*surface,NAME=ASURF
ETOP,S1
*surface,NAME=BSURF
EBOT,S3
*CONTACT PAIR,INTERACTION=SMOOTH,TYPE=SURFACE
ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
** clearance dependent gap conductance not present
*GAP CONDUCTANCE,PRESSURE
*GAP HEAT GENERATION
**
*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
*AMPLITUDE, NAME=REMOVE, VALUE=RELATIVE, TIME=STEP
0.0, 1.0,  30.0, 0.0,  100.0, 0.0
*AMPLITUDE, NAME=FILM, VALUE=RELATIVE, TIME=STEP
0.0, 0.2,  10.0, 1.0
**
*INITIAL CONDITIONS,TYPE=TEMP
ALL,0.0
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 1
*COUPLED TEMPERATURE,DELTMX=100.
*CONTROLS,PARAMETERS=FIELD
1.0E-6,
*BOUNDARY
BOT,1,1
1,2,2
7,2,2
113,2,2
*DLOAD
101,P3,100.
102,P3,100.
*DFLUX
101,S3,100.
102,S3,100.
*FILM, FILM AMPLITUDE=FILM
2,F1,0.,10.
*CONTACT PRINT,SLAVE=ASURF
*PRINT,CONTACT=YES
*OUTPUT,FIELD
*ELEMENT OUTPUT,VAR=PRESELECT
*NODEOUTPUT
U,RF,NT
*CONTACT OUTPUT,SLAVE=ASURF,MASTER=BSURF
CSTRESS,CDISP
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 2
*COUPLED TEMPERATURE,DELTMX=100.
*BOUNDARY,FIXED
TOP,1,1
*BOUNDARY,AMPLITUDE=SLIDE
113,2,2,1.5
*DLOAD
101,P3,100.
102,P3,100.
*DFLUX
101,S3,100.
102,S3,100.
*FILM
2,F1,0.,10.
*NODE PRINT
RF,
NT11,RFL
*EL PRINT
S,TEMP
HFL,
NFLUX,
*PRINT,CONTACT=YES
*END STEP
**
*STEP,INC=50,NLGEOM,UNSYMM=YES
Step 3
*COUPLED TEMPERATURE-DISPLACEMENT,STEADY STATE
100000.,100000.
*BOUNDARY,OP=NEW
BOT,1,1
1,2,2
7,2,2
113,2,2
*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,2,2,0.0
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES
Step 5
*COUPLED TEMPERATURE-DISPLACEMENT,DELTMX=100.
10.0, 100.0
*DFLUX, AMP=REMOVE
101, S3, 100.0
102, S3, 100.0
*BOUNDARY
113, 2,2, 10.0
*BOUNDARY,FIXED
111,1,1
113,1,1
115,1,1
*END STEP