*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH 
INTERFACE RADIATION, DEFORMABLE [C3D10MT]
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*NODE,NSET=ALL
  1, 1., 0.,-0.01
  7, 7., 0.,-0.01
1001, 1., .5,-0.01
1007, 7., .5,-0.01
 11, 1., 1.,-0.01
 17, 7., 1.,-0.01
2001, 1.,0.0, 0.5
2007, 7.,0.0, 0.5
 21, 1., 0., 1.01
 27, 7., 0., 1.01
3001, 1.,0.5, 1.01
3007, 7.,0.5, 1.01
 31, 1., 1., 1.01
 37, 7., 1., 1.01
4001, 1.,1., 0.5
4007, 7.,1., 0.5
***
5001, 1.,0.5, 0.5
5007, 7.,0.5, 0.5
***********
101, 3., 1.
105, 5., 1.
***
501,3.,1.5,
505,5.,1.5,
***
111, 3., 2.
115, 5., 2.
***
601,3.,1.0,0.5
605,5.,1.0,0.5
***
121, 3., 1., 1.0
125, 5., 1., 1.0
***
701,3.,1.5,1.0
705,5.,1.5,1.0
***
131, 3., 2., 1.0
135, 5., 2., 1.0
***
801,3.,2,0.5
805,5.,2,0.5
**
901,3,1.5,0.5
905,5,1.5,0.5
*NGEN, NSET=ALL
  1,  7
 11, 17
 21, 27
 31, 37
101,105
111,115
121,125
131,135
501,505
601,605
701,705
801,805
901,905
1001,1007
2001,2007
3001,3007
4001,4007
5001,5007
*NSET,NSET=QA_TEST
ALL
*NSET,NSET=BOT,GEN
1,7
2001,2007
21,27
*NSET,NSET=TOP,GEN
111,115
131,135
801,805
**
*ELEMENT, TYPE=C3D10MT
1,33,13,23,31,4003,5003,3003,32,4002,3002
2,3,23,13,1,2003,5003,1003,2,2002,1002
3,11,31,1,13,4001,5001,1001,12,4002,1002
4,21,1,31,23,2001,5001,3001,22,2002,3002
******
******
101,133,113,123,131,803,903,703,132,802,702
102,103,123,113,101,603,903,503,102,602,502
103,111,131,101,113,801,901,501,112,802,502
104,121,101,131,123,601,901,701,122,602,702
**
*ELEMENT,TYPE=C3D10MT,ELSET=DEAD5
25,5,27,17,35,2006,5007,1006,5005,3006,4006
35,3,25,15,33,2004,5005,1004,5003,3004,4004
45,1,23,13,31,2002,5003,1002,5001,3002,4002
*ELEMENT,TYPE=C3D10MT,ELSET=ACTIVE5
125,101,123,113,131,602,903,502,901,702,802
135,103,125,115,133,604,905,504,903,704,804
**
*ELGEN, ELSET=DEAD1
1,3,2,4
*ELGEN, ELSET=DEAD2
2,3,2,4
*ELGEN, ELSET=DEAD3
3,3,2,4
*ELGEN, ELSET=DEAD4
4,3,2,4
******
*ELSET, ELSET=DEAD
DEAD1,DEAD2,DEAD3,DEAD4,DEAD5
*SOLID SECTION,ELSET=DEAD,MATERIAL=DEAD
*MATERIAL,NAME=DEAD
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELGEN, ELSET=ACTIVE1
101,2,2,4
*ELGEN, ELSET=ACTIVE2
102,2,2,4
*ELGEN, ELSET=ACTIVE3
103,2,2,4
*ELGEN, ELSET=ACTIVE4
104,2,2,4
*ELSET, ELSET=ACTIVE
ACTIVE1,ACTIVE2,ACTIVE3,ACTIVE4,ACTIVE5
*SOLID SECTION,ELSET=ACTIVE,MATERIAL=ACTIVE
*MATERIAL,NAME=ACTIVE
*ELASTIC
30.E6,.3
*CONDUCTIVITY
10.000, 
*SPECIFIC HEAT
0.001, 
*DENSITY
1000., 
**
*ELSET,ELSET=ETOP, GENERATE
102,108,2
*ELSET,ELSET=EBOT, GENERATE
1,11,2
*ELSET, ELSET=ELOAD, GENERATE
101,107,2
*ELSET, ELSET=EFILM
6,8
**
*AMPLITUDE, NAME=SLIDE, VALUE=RELATIVE, TIME=STEP, DEFINITION=SMOOTH STEP
0.0, 0.0,  20.0, 1.0,  40., 0.0,  60., -1.0
80., 0.0,  100.0, 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
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFAN=1.E-6
*SURFACE, NAME=ASURF,TYPE=ELEMENT
ETOP,S2
*SURFACE, NAME=BSURF,TYPE=ELEMENT
EBOT,S2
**
*STEP
STEP 1
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.45e-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,2,3
TOP,3,3
1,1,1
7,1,1
21,1,1
27,1,1
113,1,1
133,1,1
*DLOAD, AMP=RAMP
ELOAD,P2,100.
*DFLUX, AMP=RAMP
ELOAD,S2, 100.
*FILM, AMP=RAMP
EFILM,F2,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, TIME INTERVAL=1.0
*ENERGY OUTPUT
ALLAE,ALLHF,ALLMW,ALLIE,ALLIHE,
ALLPW,ALLKE,ALLSE,ETOTAL,ALLVD
******************
*END STEP
**
*STEP
STEP 2
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.35e-3,100.
*BOUNDARY,TYPE=VELOCITY
TOP,2,2
*BOUNDARY,AMPLITUDE=SLIDE
113,1,1, 1.5
133,1,1, 1.5
*END STEP
**
*STEP
STEP 3 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT,DIRECT
0.45e-3,5.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*BOUNDARY,OP=NEW,TYPE=VELOCITY
BOT,2,3
TOP,3,3
1,1,1
7,1,1
21,1,1
27,1,1
113,1,1
133,1,1
*END STEP
**
*STEP
STEP 4 RUN TO STEADY STATE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.45e-3,10.
*FIXED MASS SCALING, ELSET=EALL, FACTOR=1.E2
*AMPLITUDE,NAME=RAMP2
0.0,0.5,10.0,1.0
*DLOAD, AMP=RAMP2
ELOAD, P2,200.
*END STEP
**
*STEP
STEP 5
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
0.175e-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
ELOAD, S2, 100.
*BOUNDARY, TYPE=VELOCITY
113, 1,1, 0.1
133, 1,1, 0.1
*BOUNDARY,FIXED
TOP,2,2
*END STEP