*HEADING
SIMPLE CONTACT TEST WITH C3D8T ELEMENTS, 
DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
**
** See ch. 1.6.2 of ABAQUS/Verification Manual
**
*NODE
1,1.,5.
3,2.,5.
9,2.,10.
7,1.,10.
11,1.,10.0000001
13,2.,10.0000001
19,2.,15.0000001
17,1.,15.0000001
101,1.,5.,1.
103,2.,5.,1.
109,2.,10.,1.
107,1.,10.,1.
111,1.,10.0000001,1.
113,2.,10.0000001,1.
119,2.,15.0000001,1.
117,1.,15.0000001,1.
*NGEN,NSET=NALL
1,3
3,9,3
7,9
1,7,3
2,8,3
11,13
13,19,3
17,19
11,17,3
12,18,3
101,103
103,109,3
107,109
101,107,3
102,108,3
111,113
113,119,3
117,119
111,117,3
112,118,3
*NSET,NSET=BOTTOM
1,2,3
101,102,103
*NSET,NSET=TOP
17,18,19
117,118,119
*NSET,NSET=SIDE
1,4,7,11,14,17
101,104,107,111,114,117
*ELEMENT,TYPE=C3D8T,ELSET=SOLID
1, 1,3,9,7,101,103,109,107
2, 11,12,15,14,111,112,115,114
3, 12,13,16,15,112,113,116,115
4, 14,15,18,17,114,115,118,117
5, 15,16,19,18,115,116,119,118
*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, 3, 7, 9
101, 103, 107, 109
*NSET, NSET=NELEM2, GEN
11,19,1
111,119,1
*NSET,NSET=QA_TEST
NELEM2
**
** surfaces
**
*SURFACE, NAME=ASURF
1, S5
*SURFACE, NAME=BSURF
2, S3
3, S3
**
*** 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
17, 1, 1
101,1, 1
117,1, 1
101, 3, 3
117, 3, 3
*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
*************
*CONTACT OUTPUT
HFLA, HFL, HTL,HTLA,SFDR,SFDRA,SFDRT,SFDRTA,
CSTRESS,
*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
14, 1, 1, 0.000025
17, 1, 1, 0.000025
111, 1, 1, 0.000025
114, 1, 1, 0.000025
117, 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
7, 1, 1
101, 1, 1
107, 1, 1
101, 3, 3
117, 3, 3
*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
14, 1, 1, -0.000025
17, 1, 1, -0.000025
111, 1, 1, -0.000025
114, 1, 1, -0.000025
117, 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
7, 1, 1
101, 1, 1
107, 1, 1
101, 3, 3
117, 3, 3
*END STEP