*HEADING
THERMALLY COUPLED CONTACT, FINITE SLIDING CONTACT WITH 
INTERFACE CONDUCTIVITY, DEFORMABLE [C3D10MT]
**
** See ch. 1.6.7 of ABAQUS/Verification Manual
**
*Preprint, Model=yes
*NODE,NSET=ALL
  1, 1., 0.,-0.01
  2, 2., 0.,-0.01
  3, 3., 0.,-0.01
  4, 4., 0.,-0.01
  5, 5., 0.,-0.01
  6, 6., 0.,-0.01
  7, 7., 0.,-0.01
 11, 1., 1.,-0.01
 12, 2., 1.,-0.01
 13, 3., 1.,-0.01
 14, 4., 1.,-0.01
 15, 5., 1.,-0.01
 16, 6., 1.,-0.01
 17, 7., 1.,-0.01
 21, 1., 0., 1.01
 22, 2., 0., 1.01
 23, 3., 0., 1.01
 24, 4., 0., 1.01
 25, 5., 0., 1.01
 26, 6., 0., 1.01
 27, 7., 0., 1.01
 31, 1., 1., 1.01
 32, 2., 1., 1.01
 33, 3., 1., 1.01
 34, 4., 1., 1.01
 35, 5., 1., 1.01
 36, 6., 1., 1.01
 37, 7., 1., 1.01
101, 3., 1.,
102, 3.5, 1.,
103, 4., 1.,
104, 4.5, 1.,
105, 5., 1.,
 111, 3., 2.,
 112, 3.5, 2.,
 113, 4., 2.,
 114, 4.5, 2.,
 115, 5., 2.,
121, 3., 1., 1.0
122, 3.5, 1., 1.0
123, 4., 1., 1.0
124, 4.5, 1., 1.0
125, 5., 1., 1.0
 131, 3., 2., 1.0
 132, 3.5, 2., 1.0
 133, 4., 2., 1.0
 134, 4.5, 2., 1.0
 135, 5., 2., 1.0
**
501,3.,1.5,
502,3.5,1.5,
503,4.,1.5,
504,4.5,1.5,
505,5.,1.5,
**
601,3.,1.0,0.5
602,3.5,1.0,0.5
603,4.,1.0,0.5
604,4.5,1.0,0.5
605,5.,1.0,0.5
***
701,3.,1.5,1.0
702,3.5,1.5,1.0
703,4.,1.5,1.0
704,4.5,1.5,1.0
705,5.,1.5,1.0
***
801,3.,2,0.5
802,3.5,2,0.5
803,4.,2,0.5
804,4.5,2,0.5
805,5.,2,0.5
**
901,3,1.5,0.5
902,3.5,1.5,0.5
903,4.0,1.5,0.5
904,4.5,1.5,0.5
905,5,1.5,0.5
**
1001, 1., .5,-0.01
1002, 2., .5,-0.01
1003, 3., .5,-0.01
1004, 4., .5,-0.01
1005, 5., .5,-0.01
1006, 6., .5,-0.01
1007, 7., .5,-0.01
**
2001, 1.,0.0, 0.5
2002, 2.,0.0, 0.5
2003, 3.,0.0, 0.5
2004, 4.,0.0, 0.5
2005, 5.,0.0, 0.5
2006, 6.,0.0, 0.5
2007, 7.,0.0, 0.5
**
3001, 1.,0.5, 1.01
3002, 2.,0.5, 1.01
3003, 3.,0.5, 1.01
3004, 4.,0.5, 1.01
3005, 5.,0.5, 1.01
3006, 6.,0.5, 1.01
3007, 7.,0.5, 1.01
**
4001, 1.,1., 0.5
4002, 2.,1., 0.5
4003, 3.,1., 0.5
4004, 4.,1., 0.5
4005, 5.,1., 0.5
4006, 6.,1., 0.5
4007, 7.,1., 0.5
***
5001, 1.,0.5, 0.5
5002, 2.,0.5, 0.5
5003, 3.,0.5, 0.5
5004, 4.,0.5, 0.5
5005, 5.,0.5, 0.5
5006, 6.,0.5, 0.5
5007, 7.,0.5, 0.5
***********
** new node set for QA_TEST
*************************
*NSET,NSET=INTERFACE,GEN
11,17
31,37
101,105
121,125
601,603
4001,4007
************************
*NSET,NSET=QA_TEST
INTERFACE
*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
5,35,15,25,33,4005,5005,3005,34,4004,3004
6,5,25,15,3,2005,5005,1005,4,2004,1004
7,13,33,3,15,4003,5003,1003,14,4004,1004
8,23,3,33,25,2003,5003,3003,24,2004,3004
9,37,17,27,35,4007,5007,3007,36,4006,3006
10,7,27,17,5,2007,5007,1007,6,2006,1006
11,15,35,5,17,4005,5005,1005,16,4006,1006
12,25,5,35,27,2005,5005,3005,26,2006,3006
******
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
******
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
**
105,135,115,125,133,805,905,705,134,804,704
106,105,125,115,103,605,905,505,104,604,504
107,113,133,103,115,803,903,503,114,804,504
108,123,103,133,125,603,903,703,124,604,704
**
125,101,123,113,131,602,903,502,901,702,802
135,103,125,115,133,604,905,504,903,704,804
*ELSET,ELSET=DEAD5
25,35,45
*ELSET,ELSET=ACTIVE5
125,135
**
*ELSET, ELSET=DEAD1
1,5,9
*ELSET, ELSET=DEAD2
2,6,10
*ELSET, ELSET=DEAD3
3,7,11
*ELSET, ELSET=DEAD4
4,8,12
******
*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., 
**
*ELSET, ELSET=ACTIVE1
101,105
*ELSET, ELSET=ACTIVE2
102,106,
*ELSET, ELSET=ACTIVE3
103,107
*ELSET, ELSET=ACTIVE4
104,108
*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
**
*SURFACE, NAME=ASURF,TYPE=ELEMENT
ETOP,S2
*SURFACE, NAME=BSURF,TYPE=ELEMENT
EBOT,S2
**
*STEP
STEP 1
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT, DIRECT
1e-3,10.
*CONTACT PAIR,INTERACTION=ROUGH,WEIGHT=0.0
ASURF,BSURF
*SURFACE INTERACTION, NAME=ROUGH
*GAP CONDUCTANCE,PRESSURE
20.,0.
20.,100.
20.,110.
 5.,190.
 5.,200.
**
*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,VAR=PRESELECT
******************
*END STEP
**
*STEP
STEP 2
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT,DIRECT
1e-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
1e-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
1e-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
1.25e-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.
*AMPLITUDE, NAME=RAMP3, DEFINITION=SMOOTH STEP
0.0,0.0,40.,1.0
*BOUNDARY, TYPE=VELOCITY, AMP=RAMP3
113, 1,1, 0.2
133, 1,1, 0.2
TOP, 2,2
*END STEP