*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH 
INTERFACE RADIATION, DEFORMABLE [CAX6MT]
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*NODE,NSET=ALL
  1, 1., 1.
  7, 1., 7.
***
 11, 1.5, 1.
 17, 1.5, 7.
**
 21, 2., 1.
 27, 2., 7.
**
101, 2., 3.
105, 2., 5.
**
111, 2.5, 3.
115, 2.5, 5.
**
121, 3., 3.
125, 3., 5.
*NGEN,NSET=ALL
  1,  7,1
 11, 17,1
 21, 27,1
101,105,1
111,115,1
121,125,1
*NSET,NSET=QA_TEST
ALL
*NSET,NSET=BOT,GEN
1,7,1
*NSET,NSET=TOP,GEN
121,125,1
**
*ELEMENT,TYPE=CAX6MT,ELSET=DEAD
1,1,21,23,11,22,12
2,1,23,3,12,13,2
5,5,25,27,15,26,16
6,5,27,7,16,17,6
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELEMENT,TYPE=CAX6MT,ELSET=ACTIVE
3,3,23,25,13,24,14
4,3,25,5,14,15,4
101,101,121,123,111,122,112
102,103,101,123,102,112,113
103,103,123,125,113,124,114
104,105,103,125,104,114,115
*ELSET,ELSET=ETOP
102,104
*ELSET,ELSET=EBOT, GENERATE
1,5,2
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFAN=1.E-6
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10., 
*DENSITY
1000., 
*SPECIFIC HEAT
.001, 
**
*AMPLITUDE, NAME=SLIDE, VALUE=RELATIVE, TIME=STEP,DEFINITION=SMOOTH STEP
0.0, 0.0,  20., 1.0,  40., 0.0,  60., -1.0
80., 0.0,  100., 0.0
*AMPLITUDE, NAME=RAMP, DEFINITION=SMOOTH STEP
0.0,0.0,10.,1.0
**
*INITIAL CONDITIONS,TYPE=TEMP
ALL,20.0
**
*ELSET,ELSET=EALL
ACTIVE,DEAD
**
*SURFACE, NAME=ASURF
ETOP,S1
*SURFACE, NAME=BSURF
EBOT,S2
**
*STEP
STEP 1
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.5e-3,10.
**
*CONTACT PAIR,INTERACTION=ROUGH,WEIGHT=0.0
ASURF,BSURF
*SURFACE INTERACTION,NAME=ROUGH
*GAP RADIATION
1.0, 1.0
1.0,  0.0
1.0, 10.0
*GAP HEAT GENERATION
0.0, 
**
*BOUNDARY, TYPE=VELOCITY
BOT,1,1
1,2,2
7,2,2
123,2,2
*DLOAD,AMP=RAMP
101,P2,100.
103,P2,100.
*DFLUX,AMP=RAMP
101,S2,100.
103,S2,100.
*FILM,AMP=RAMP
4,F3,20.,10.
*********************
*FILE OUTPUT,N=1
*NODE FILE,NSET=ALL
NT,
********************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
NT,
*NODE OUTPUT,NSET=ALL
U,
************
*CONTACT OUTPUT
CSTRESS,
*ELEMENT OUTPUT
S,
********************
***OUTPUT,HISTORY, FREQUENCY=500
***NODE OUTPUT,NSET=ALL
**U,NT
*******************
*OUTPUT,HISTORY, FREQUENCY=99999
*ENERGY OUTPUT
ALLAE, ALLHF, ALLMW, ALLFD, ALLIHE, ALLPW, ALLVD, ETOTAL
******************
*END STEP
**
*STEP
STEP 2
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.5e-3,100.
*BOUNDARY, TYPE=VELOCITY
TOP,1,1
*BOUNDARY,AMPLITUDE=SLIDE
123,2,2,1.5
*END STEP
**
*STEP
STEP 3 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.5e-3,5.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*BOUNDARY,OP=NEW, TYPE=VELOCITY
BOT,1,1
1,2,2
7,2,2
123,2,2
*END STEP
**
*STEP
STEP 4 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.5e-3,10.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE,NAME=RAMP2, DEFINITION=SMOOTH STEP
0.0,0.5,10.0,1.0
*DLOAD,AMP=RAMP2
101,P2,200.
103,P2,200.
*END STEP
**
*STEP
STEP 5
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT,DIRECT
0.5e-3, 40.0
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE, NAME=REMOVE, VALUE=RELATIVE, TIME=STEP, DEFINITION=SMOOTH STEP
0.0, 1.0,  30.0, 0.0,  40.0, 0.0
*DFLUX, AMP=REMOVE
101, S2, 100.0
103, S2, 100.0
*BOUNDARY,TYPE=VELOCITY
123, 2,2, 0.1
121,1,1
122,1,1
123,1,1
124,1,1
125,1,1
*END STEP