*HEADING
SIMPLE CONTACT TEST WITH CPS6MT ELEMENTS, 
DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
*NODE,NSET=NALL
1,1.,5.
3,2.0001,5.
7,1.,10.
9,2.0001,10.
**************
11,1.,10.0000001
15,2.,10.0000001
31,1.,15.0000001
35,2.,15.0000001
*NGEN,NSET=NALL
1,3,1
7,9,1
3,9,3
1,7,3
4,6,1
2,8,3
***
11,15,1
15,35,5
31,35,1
11,31,5
21,25,1
16,20,1
26,30,1
14,34,5
12,32,5
*NSET,NSET=BOTTOM
1,2,3
*NSET,NSET=TOP
31,32,33,34,35
*NSET,NSET=SIDE
1,4,7,11,16,21,26,31
*ELEMENT,TYPE=CPS6MT,ELSET=SOLID
1, 1,3,9,2,6,5
2, 1,9,7,5,8,4
101, 11,13,23,12,18,17
102, 13,15,25,14,20,19
103, 21,23,33,22,28,27
104, 23,25,35,24,30,29
105, 11,23,21,17,22,16
106, 13,25,23,19,24,18
107, 21,33,31,27,32,26
108, 23,35,33,29,34,28
*SOLID SECTION,ELSET=SOLID,MATERIAL=SOLMAT
*MATERIAL,NAME=SOLMAT
*ELASTIC
30.E6,.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=NELEM1
1,2,3,4,5,6,7,8,9
*NSET, NSET=NELEM2, GEN
11,35,1
*NSET,NSET=QA_TEST
NELEM2
**
** surfaces
**
*SURFACE, NAME=ASURF
2, S2
*SURFACE, NAME=BSURF
101, S1
102, 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
BOTTOM,2,2
TOP,2,2
1, 1, 1
31, 1, 1
*amplitude,name=ramp, definition=smooth step
0.,0.,100.,1.,500.,1.
*boundary,amplitude=ramp
NELEM1, 11, 11, 700.0
NELEM2, 11, 11, 500.0
************************
*FILE OUTPUT, NUM=1
*NODE FILE,NSET=NELEM2
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
11, 1, 1, 0.000025
16, 1, 1, 0.000025
21, 1, 1, 0.000025
26, 1, 1, 0.000025
31, 1, 1, 0.000025
*BOUNDARY, OP=NEW
NELEM1, 11, 11, 700
NELEM2, 11, 11, 500
*BOUNDARY, OP=NEW,type=velocity
TOP, 2, 2
BOTTOM, 2, 2
1, 1, 1
4, 1, 1
7, 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
11, 1, 1, -0.000025
16, 1, 1, -0.000025
21, 1, 1, -0.000025
26, 1, 1, -0.000025
31, 1, 1, -0.000025
*BOUNDARY, OP=NEW
NELEM1, 11, 11, 700
NELEM2, 11, 11, 500
*BOUNDARY, OP=NEW,type=velocity
TOP, 2, 2
BOTTOM, 2, 2
1, 1, 1
4, 1, 1
7, 1, 1
*END STEP