*HEADING
SELF CONTACT TEST, COUPLED LINEAR UNDERLYING ELEMENTS, CLOSED SURFACE CPE4RT   
PENALTY CONTACT
**
*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, REF NODE=1001
*RIGID BODY, REF NODE=1002, ANALYTICAL SURFACE =RBASE, REF NODE=1002	
**
*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,MECHANICAL CONSTRAINT=PENALTY
OUTSIDE, RBASE
*SURFACE INTERACTION, NAME=BASE
1.0, 
*FRICTION,ROUGH
*CONTACT PAIR, INTERACTION=TOP,MECHANICAL CONSTRAINT=PENALTY
OUTSIDE, RTOP
*SURFACE INTERACTION, NAME=TOP
1.0, 
*FRICTION,ROUGH
*CONTACT PAIR, INTERACTION=INSIDE,MECHANICAL CONSTRAINT=PENALTY
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