*HEADING
SELF CONTACT TEST, COUPLED LINEAR UNDERLYING ELEMENTS, CLOSED SURFACE CPE4RT   
**
*NODE,NSET=NALL
       1,          3.,          0.
      73,          3.,          0.
     101,        3.25,          0.
     173,        3.25,          0.
     201,         3.5,          0.
     273,         3.5,          0.
     301,        3.75,          0.
     373,        3.75,          0.
     401,          4.,          0.
     473,          4.,          0.
*NGEN,LINE=C,NSET=NALL
1,73,1,,0.,0.,0.,0.,0.,1.
101,173,1,,0.,0.,0.,0.,0.,1.
201,273,1,,0.,0.,0.,0.,0.,1.
301,373,1,,0.,0.,0.,0.,0.,1.
401,473,1,,0.,0.,0.,0.,0.,1.
**
*NSET,NSET=QA_TEST
NALL,
*NODE, NSET=ODB_NSET
    1001,     0.00000,     5.00000
    1002,     0.00000,    -5.00000
**
**
*ELEMENT, TYPE=CPE4RT, ELSET=WHEEL
1,1,101,102,2
72,72,172,101,1
*ELGEN,ELSET=WHEEL
1,71,1,1,4,100,100
72,4,100,100
**
*SOLID SECTION, ELSET=WHEEL, MATERIAL=SIMPLE,CONTROLS=HGLASS
1.,
*SECTION CONTROLS,NAME=HGLASS,HOURGLASS=STIFFNESS
1.0,
*MATERIAL,NAME=SIMPLE
*HYPERELASTIC,N=1
1E3,,1.E-3 
**
*CONDUCTIVITY
0.0005,
*SPECIFIC HEAT
0.1,
*DENSITY
1., 
*NSET,NSET=COLD,GEN
 38, 72,1
138,172,1
238,272,1
338,372,1
438,472,1
*NSET,NSET=HOT,GEN
  1, 37,1
101,137,1
201,237,1
301,337,1
401,437,1
*INITIAL CONDITIONS,TYPE=TEMPERATURE
HOT,100.
*INITIAL CONDITIONS,TYPE=TEMPERATURE
COLD,0.   
**
*ELSET, ELSET=INSIDE, GEN
1,72,1
*ELSET, ELSET=OUTSIDE, GEN
301,372,1
**
**
*RIGID BODY, REF NODE=1001, ANALYTICAL SURFACE =RTOP
*RIGID BODY, REF NODE=1002, ANALYTICAL SURFACE =RBASE	
**
*ELSET,ELSET=ODB_ESET
WHEEL,
*NSET, NSET=ODB_NSET
NALL, NALL, ODB_NSET, COLD
HOT, 
*SURFACE,TYPE=ELEMENT, NAME=INSIDE
INSIDE,S4
*SURFACE,TYPE=ELEMENT, NAME=OUTSIDE
OUTSIDE,S2
*SURFACE, TYPE=SEGMENTS, NAME=RTOP
START, 0.0, 6.0
CIRCL, 1.0, 5.0, 0.0, 5.0
CIRCL, 0.0, 4.0, 0.0, 5.0
CIRCL,-1.0, 5.0, 0.0, 5.0
CIRCL, 0.0, 6.0, 0.0, 5.0
*SURFACE, TYPE=SEGMENTS, NAME=RBASE
START,-15.0, -4.0
LINE,   4.0, -4.0
**
**
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
 ,100.0
*FIXED MASS SCALING,ELSET=WHEEL,FACTOR=40.0
*CONTACT PAIR, INTERACTION=BASE
OUTSIDE, RBASE
*SURFACE INTERACTION, NAME=BASE
1.0, 
*FRICTION,ROUGH
*CONTACT PAIR, INTERACTION=TOP
OUTSIDE, RTOP
*SURFACE INTERACTION, NAME=TOP
1.0, 
*FRICTION,ROUGH
*CONTACT PAIR, INTERACTION=INSIDE
INSIDE,
*SURFACE INTERACTION, NAME=INSIDE
1.0, 
*GAP CONDUCTANCE
0.0005,0.
0.,0.001
0.,10.
*BOUNDARY,TYPE=VELOCITY
219,1,1
255,1,1
1002, 1,2
1002, 6,6
1001, 6,6
1001, 1,1
1001, 2,,-6.2E-02
*FILE OUTPUT, NUM=1
*NODE FILE
U,NT
***********************
*OUTPUT,FIELD,NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,NT
************************
*OUTPUT,FIELD,VARIABLE=PRESELECT,NUMBER INTERVAL=1
*OUTPUT, HISTORY, TIME INTERVAL=100
*ENERGY OUTPUT,VARIABLE=PRESELECT
*ENDSTEP
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT,EXPLICIT
 ,100.0
*FIXED MASS SCALING,ELSET=WHEEL,FACTOR=40.0
*BOUNDARY, TYPE=VELOCITY, OP=NEW
1001,1,1,-10.0E-02
1001,2,2,  0.0
1001,6,6, -8.0E-02
1002,1,2
1002,6,6
*ENDSTEP