*HEADING
SIMPLE CONTACT TEST WITH CPE4RT ELEMENTS, 
DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
*NODE
1,1.,5.
3,2.0001,5.
9,2.0001,10.
7,1.,10.
11,1.,10.0000001
13,2.,10.0000001
19,2.,15.0000001
17,1.,15.0000001
*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
*NSET,NSET=BOTTOM
1,2,3
*NSET,NSET=TOP
17,18,19
*NSET,NSET=SIDE
1,4,7,11,14,17
*ELEMENT,TYPE=CPE4RT,ELSET=SOLID
1, 1,3,9,7
2, 11,12,15,14
3, 12,13,16,15
4, 14,15,18,17
5, 15,16,19,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
*NSET, NSET=NELEM2, GEN
11,19,1
*NSET,NSET=QA_TEST
NELEM2
**
** surfaces
**
*SURFACE, NAME=ASURF
1, S3
*SURFACE, NAME=BSURF
2, S1
3, 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
17, 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
*************
*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
*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
*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
*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
*END STEP