*HEADING
**  Test thermal contact with analytical rigid surface
**
*NODE
      1,           0.,          0.2,           0.
      2,    -0.533333,          0.2,           0.
      3,     -1.06667,          0.2,           0.
      4,         -1.6,          0.2,           0.
      5,           0.,          0.2,          0.5
      6,    -0.533333,          0.2,          0.5
      7,     -1.06667,          0.2,          0.5
      8,         -1.6,          0.2,          0.5
      9,           0.,          0.2,           1.
     10,    -0.533333,          0.2,           1.
     11,     -1.06667,          0.2,           1.
     12,         -1.6,          0.2,           1.
     13,           0.,           0.,           0.
     14,    -0.533333,           0.,           0.
     15,     -1.06667,           0.,           0.
     16,         -1.6,           0.,           0.
     17,           0.,           0.,          0.5
     18,    -0.533333,           0.,          0.5
     19,     -1.06667,           0.,          0.5
     20,         -1.6,           0.,          0.5
     21,           0.,           0.,           1.
     22,    -0.533333,           0.,           1.
     23,     -1.06667,           0.,           1.
     24,         -1.6,           0.,           1.
     25,         -2.3,          0.4,          0.5
*ELEMENT, TYPE=C3D8RT, ELSET=SOLID
1, 13, 14, 18, 17,  1,  2,  6,  5
2, 14, 15, 19, 18,  2,  3,  7,  6
3, 15, 16, 20, 19,  3,  4,  8,  7
4, 17, 18, 22, 21,  5,  6, 10,  9
5, 18, 19, 23, 22,  6,  7, 11, 10
6, 19, 20, 24, 23,  7,  8, 12, 11
*SURFACE, TYPE=CYLINDER, NAME=RIGIDSURFACE
-1.6,0.2,1., 1.,0.2,1.
-1.6,0.2,0.
START,         -0.7,          0.2
 LINE,         -0.9,          0.2
*SURFACE, TYPE=ELEMENT, NAME=SOLID
SOLID, S2
*NSET,NSET=REF
25,
*RIGID BODY, REF NODE=REF, ANALYTICAL SURFACE=RIGIDSURFACE, ISOTHERMAL=YES
*MATERIAL, NAME=MATERIAL-1
*CONDUCTIVITY
1.,
*DENSITY
2500.,
*ELASTIC
 3E+06, 0.3
*SPECIFIC HEAT
100.,
*EXPANSION
0.,
*SOLID SECTION, ELSET=SOLID, MATERIAL=MATERIAL-1
*BOUNDARY
24, 1, 3
16, 1, 2
13, 2, 2
21, 2, 2
REF, 2, 6
*SURFACE INTERACTION, NAME=INTPROP-1
*GAP CONDUCTANCE
1000.,0.
 0.,0.5
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
, 1.
*BOUNDARY, TYPE=VELOCITY
REF, 1, 1, 2.
*BOUNDARY
REF, 11, 11, 100.
*CONTACT
*CONTACT INCLUSION, ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
SOLID,RIGIDSURFACE,INTPROP-1
*NSET,NSET=QA_TESTN,GEN
1,25,1
*OUTPUT, FIELD, NUM=1
*NODE OUTPUT, NSET=QA_TESTN
U, NT
*OUTPUT, FIELD, NUM=10
*NODE OUTPUT
U, NT
*CONTACT OUTPUT
HFL,HTL,CSTRESS
*END STEP