*HEADING
SIMPLE CONTACT TEST WITH C3D10MT ELEMENTS, 
DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
*RESTART,WRITE,NUM=1
*NODE
       1,          1.,          5.
       2,         1.5,          5.
       3,          2.,          5.
       4,          1.,         7.5
       5,          2.,         7.5
       6,          1.,         10.
       7,         1.5,         10.
       8,          2.,         10.
       9,          1.,          5.,         0.5
      10,          2.,          5.,         0.5
      11,          1.,         10.,         0.5
      12,          2.,         10.,         0.5
      13,          1.,          5.,          1.
      14,         1.5,          5.,          1.
      15,          2.,          5.,          1.
      16,          1.,         7.5,          1.
      17,          2.,         7.5,          1.
      18,          1.,         10.,          1.
      19,         1.5,         10.,          1.
      20,          2.,         10.,          1.
      21,         1.5,          5.,         0.5
      22,         1.5,         10.,         0.5
      23,         1.5,         7.5,         
      24,         1.5,         7.5,          1.
      25,          2.,         7.5,         0.5
      26,          1.,         7.5,         0.5       
     121,          1.,        12.5
     122,          1.,        11.25
     123,          1.,        10.,
     124,        1.25,        12.5,
     125,        1.25,        10.,
     126,         1.5,        12.5
     127,         1.5,        11.25
     128,         1.5,        10.
     129,          1.,        12.5,         0.5
     130,          1.,        10.,          0.5
     131,         1.5,        12.5,         0.5
     132,         1.5,        10.,          0.5
     133,          1.,        12.5,          1.
     134,          1.,        11.25,         1.
     135,          1.,        10.,           1.
     136,        1.25,        12.5,          1.
     137,        1.25,        10.,           1.
     138,         1.5,        12.5,          1.
     139,         1.5,        11.25,         1.
     140,         1.5,        10.,           1.,
     144,        1.75,        12.5,
     145,        1.75,        10.
     146,          2.,        12.5
     147,          2.,        11.25
     148,          2.,        10.
     151,          2.,        12.5,          0.5
     152,          2.,        10.,           0.5
     156,        1.75,        12.5,          1.
     157,        1.75,        10.,           1.
     158,          2.,        12.5,          1.
     159,          2.,        11.25,         1.
     160,          2.,        10.,           1.
     161,        1.25,        10.,           0.5
     162,        1.25,        12.5,          0.5
     163,        1.25,        11.25,
     164,        1.25,        11.25,         1.
     165,        1.5,         11.25,         0.5
     166,          1.,        11.25,         0.5
     167,        1.75,        10.,           0.5
     168,        1.75,        12.5,          0.5
     169,        1.75,        11.25, 
     170,        1.75,        11.25,         1.
     171,          2.,        11.25,         0.5
*NSET,NSET=BOTTOM
 1,2,3,9,10,13,14,15,21
*NSET,NSET=TOP
121,124,126,144,146,129,
162,131,168,151,133,136,
138,156,158
*NSET,NSET=BACK
 1,2,3,
121,124,126,144,146
*ELEMENT, TYPE=C3D10MT, ELSET=SOLID
  1,    15,    3,    1,    8,   10,    2,   21,
        25,    5,   23   
  2,     1,   13,   15,   18,    9,   14,   21,
        26,   16,   24
  3,     8,   20,   18,   15,   12,   19,   22,
        25,   17,   24
  4,    18,    6,    8,    1,   11,    7,   22,
        26,    4,   23
  5,    15,    8,    1,   18,   25,   23,   21, 
        24,   22,   26
  6,   140,  128,  123,  126,  132,  125,  161,
       165,  127,  163
  7,   123,  135,  140,  133,  130,  137,  161,
       166,  134,  164
  8,   126,  138,  133,  140,  131,  136,  162,  
       165,  139,  164
  9,   133,  121,  126,  123,  129,  124  ,162,
       166,  122,  163
 10,   140,  126,  123,  133,  165,  163,  161,
       164,  162,  166
 11,   160,  148,  128,  146,  152,  145,  167,
       171,  147,  169
 12,   128,  140,  160,  138,  132,  157,  167,
       165,  139,  170
 13,   146,  158,  138,  160,  151,  156,  168,
       171,  159  ,170
 14,   138,  126,  146,  128,  131,  144,  168,
       165,  127,  169
 15,   160,  146,  128,  138,  171,  169,  167,
       170,  168,  165
*SOLID SECTION, ELSET=SOLID, MATERIAL=SOLMAT
*MATERIAL, NAME=SOLMAT
*ELASTIC
 30.E6,0.3
*EXPANSION
 10.E-6,
*CONDUCTIVITY
 43.0, 
*DENSITY
 7700., 
*SPECIFIC HEAT
 600., 
**
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFANB=5.6697E-8
**
*** nset definitions for prescribing temperatures.
**
*NSET,NSET=NELEM15,GEN
 1,26
*NSET, NSET=QA_TEST
NELEM15,
*NSET, NSET=NELEM620
 122,123,125,127,128,130
 132,134,135,137,139,140
 141,142,143,145,147,148,149,150
 152,153,154,155,157,159,160
 161,163,164,165,166,167,169,170
** surfaces
**
*surface, NAME=ASURF
 3, S1
 4, S1
*surface, NAME=BSURF
 6, S1
 7, S1
 11,S1
 12,S1
**
*** Apply different temperatures to the nodes above and
*** below the contact surface
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
 , 500.
*fixed mass scaling, elset=solid,factor=100.0
*CONTACT PAIR, SMALL SLIDING, INTERACTION=GAP, WEIGHT=0.0
BSURF, ASURF
*SURFACE INTERACTION, NAME=GAP
 1.0,
*FRICTION
0.01,  
*GAP CONDUCTANCE
 1.,0.
 1.,1.
*GAP RADIATION
 3.4664e-2, 3.4664e-2
 1.,0.
 1.,1.
*GAP HEAT GENERATION
 0.5,0.5
*BOUNDARY,type=velocity
 BACK,3,3
 BOTTOM,2,2
 TOP,2,2
 1, 1, 1
 9, 1, 1
 13,1, 1
*AMPLITUDE,NAME=RAMP, DEFINITION=SMOOTH STEP
0.,0.,100.,1.,500.,1.
*BOUNDARY,AMPLITUDE=RAMP
NELEM15, 11, 11, 700.0
NELEM620, 11, 11, 500.0
************************
*FILE OUTPUT, NUM=1
*NODE FILE,NSET=NELEM620
NT,U, 
************************
*OUTPUT,FIELD,NUM=1
*NODE OUTPUT,NSET=QA_TEST
NT,U,
*************
*NODE OUTPUT
U,NT
*ELEMENT OUTPUT
S
**********************************
*OUTPUT,HISTORY,TIME INTERVAL=500.0
*ENERGY OUTPUT, VAR=ALL
**********************************
*END STEP
**
**** STEP 2. Push top block to the right.
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
 ,500.0
*fixed mass scaling, elset=solid,factor=100.0
**
*amplitude,name=ramp2,definition=smooth step
0.,0., 100.,1.0, 400.,1.0, 500.,0.
*BOUNDARY, OP=NEW, amplitude=ramp2,type=velocity
 123, 1, 1, 0.000025
 130, 1, 1, 0.000025
 135, 1, 1, 0.000025
 122, 1, 1, 0.000025
 134, 1, 1, 0.000025
 166, 1, 1, 0.000025
*BOUNDARY, OP=NEW
 NELEM15, 11, 11, 700
 NELEM620, 11, 11, 500
*BOUNDARY, OP=NEW,type=velocity
 BACK,3,3
 TOP, 2, 2
 BOTTOM, 2, 2
 1, 1, 1
 9, 1, 1
 13,1, 1
 4,1,1
 16,1,1
 6,1,1
 11,1,1
 18,1,1
 26,1,1
*END STEP
**
*** STEP 3. Push top block back into position.
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
 ,500.0
*fixed mass scaling, elset=solid,factor=100.0
*BOUNDARY, OP=NEW, amplitude=ramp2,type=velocity
 123, 1, 1, -0.000025
 130, 1, 1, -0.000025
 135, 1, 1, -0.000025
 122, 1, 1, -0.000025
 134, 1, 1, -0.000025
 166, 1, 1, -0.000025
*BOUNDARY, OP=NEW
 NELEM15, 11, 11, 700
 NELEM620, 11, 11, 500
*BOUNDARY, OP=NEW,type=velocity
 BACK,3,3
 TOP, 2, 2
 BOTTOM, 2, 2
 1, 1, 1
 9, 1, 1
 13,1, 1
 4,1,1
 16,1,1
 6,1,1
 11,1,1
 18,1,1
 26,1,1
*END STEP