*HEADING
: THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH FRICTIONAL HEAT GENERATION, 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.
1001, 1., 0.5
1007, 7., 0.5
1101, 3., 1.5
1105, 5., 1.5
*NGEN
1,7,1
11,17,1
101,105,1
111,115,1
1001,1007,1
1101,1105,1
*NSET,NSET=BOT,GEN
1,7,1
*NSET,NSET=TOP,GEN
111,115,1
**
*ELEMENT,TYPE=CPE8T,ELSET=DEAD
1,1,3,13,11,2,1003,12,1001
3,5,7,17,15,6,1007,16,1005
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000,
*SPECIFIC HEAT
0.001,
*DENSITY
1000.,
**
*ELEMENT,TYPE=CPE8T,ELSET=ACTIVE
2,3,5,15,13,4,1005,14,1003
101,101,103,113,111,102,1103,112,1101
102,103,105,115,113,104,1105,114,1103
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.,
*DENSITY
1000.,
*SPECIFIC HEAT
.001,
**
*ELSET,ELSET=ETOP
101,102
*ELSET,ELSET=EBOT
1,2,3
*surface,NAME=ASURF
ETOP,S1
*surface, NAME=BSURF
EBOT,S3
*CONTACT PAIR,INTERACTION=ROUGH,SMOOTH=0.05
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*SURFACE BEHAVIOR, AUGMENTED
*FRICTION
0.1,
** clearance dependent gap conductance not present
*GAP CONDUCTANCE,PRESSURE
20.,100.
*GAP HEAT GENERATION
0.50,
**
*NSET,NSET=CONST
111,112,114,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
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES,Convert Sdi=YES
Step 1
*COUPLED TEMPERATURE,DELTMX=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,SUMMARY=YES,TOTALS=NO
*CONTACT FILE,SECONDARY=ASURF,F=10
*NODE PRINT
RF,NT11,RFL
*EL PRINT
S,E,TEMP
HFL,
NFLUX,
*EL FILE,F=10
S,E
*PRINT,CONTACT=YES
*ENERGY PRINT,F=1
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES,Convert Sdi=YES
Step 2 - Generate 120 units of heat, 60 units absorbed
*COUPLED TEMPERATURE,DELTMX=100.
.05,1.,,.05
*BOUNDARY,AMPLITUDE=SLIDE
113, 1,1, 1.5
*CONTACT FILE,SECONDARY=ASURF,F=10
*EL FILE,F=10
S,E,TEMP
NFLUX,
*NODE FILE,F=10
NT,RFL
*END STEP
**
*STEP,INC=50,NLGEOM,AMP=RAMP,UNSYMM=YES,Convert Sdi=YES
Step 3 - Go to steady state - TEMP=7.5
*COUPLED TEMPERATURE,DELTMX=100.
100.,500.
*CONTACT FILE,SECONDARY=ASURF,F=10
*EL FILE,F=10
S,E,TEMP
*NODE FILE,F=10
NT,RFL
*END STEP