*HEADING
SIMPLE CONTACT TEST WITH C3D4T ELEMENTS, 
DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
TESTs THE HEAT TRANSFER OUTPUT (SURFACE, NODE and ELEMENT)
*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=C3D4T,ELSET=SOLID
1, 109,9,103,107
2, 3,103,9,1
3, 7,107,1,9
4, 101,1,107,103
5, 1,103,9,107
**
101, 118,18,112,117
102, 12,112,18,11
103, 17,117,11,18
104, 111,11,117,112
105, 11, 112,18,117
**
106, 118,119,112,18
107, 113,112,119,13
108, 19,18,13,119
109, 12,13,18,112
110, 13, 112,119,18
**
*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,13,1
17,19,1
111,113,1
117,119,1
*NSET, NSET=QA_TEST
NELEM2
**
** surfaces
**
*SURFACE, NAME=ASURF
1, S2
3, S2
*SURFACE, NAME=BSURF
102, S2
104, S2
107, S2
109, S2
**
*** 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
17, 1, 1, 0.000025
111, 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
17, 1, 1, -0.000025
111, 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