*HEADING
HERTZ CONTACT PROBLEM WITH SUBSTRUCTURING AND CPE6M ELEMENTS 
SUBSTRUCTURE ANALYSIS
Requires substructure generation file 
hertzcontact_cpe6m_gen1.inp
*NODE
       1,          0.,        -10.
      21,    0.130871,    -9.99914
      29,    0.261719,    -9.99657
      38,    0.433182,    -9.99061
      49,    0.604518,    -9.98171
      61,    0.828872,    -9.96559
      78,     1.05281,    -9.94443
      96,     1.34557,    -9.90906
     115,     1.63716,    -9.86508
     133,     2.01724,    -9.79443
     147,      2.3943,    -9.70914
     272,         10.,          0.
300,0.,-11.
*NODE,NSET=FIX
301,0.,     -10.
302,0.26177,-10.
303,0.52336,-10.
304,0.78459,-10.
305,1.04528,-10.
306,1.30526,-10.
307,1.56434,-10.
308,1.82236,-10.
309,2.07912,-10.
*ELEMENT,TYPE=R2D2
201,301,302
*ELGEN,ELSET=RIGID
201,8,1,1
*nset,nset=nbotsur
1,21,29,38,49,61,78,96,115,133,147
*RIGID BODY,ELSET=RIGID,REF NODE=300
*surface,type=node,NAME=BOTSUR
nbotsur,1.0
*surface,NAME=RIGSUR
RIGID,SPOS
*CONTACT PAIR,INTERACTION=NO_FRI,SMALL SLIDING
BOTSUR,RIGSUR
*SURFACE INTERACTION,NAME=NO_FRI
1.0, 
**
**
*ELEMENT,TYPE=Z1002,ELSET=SUPER, FILE=hertzcontact_cpe6m_gen1
99,1,21,29,38,49,61,78,96,115,133,
147,272
*SUBSTRUCTURE PROPERTY,ELSET=SUPER,POSITION TOL=0.
*NORMAL
RIGID,FIX,,1.
*BOUNDARY
 300,1,6
 1,1
*RESTART,WRITE
*STEP
*STATIC,DIRECT
 1., 
*BOUNDARY
 272,2,2,-.4
*CONTACT PRINT ,SLAVE=BOTSUR,MASTER=RIGSUR
*CONTACT FILE ,SLAVE=BOTSUR,MASTER=RIGSUR
*PRINT,CONTACT=YES,SOLVE=YES
*NODE PRINT
U,
*OUTPUT,VAR=PRESELECT,FIELD
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT
S,
*NODE PRINT
U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*NODE PRINT,NSET=CYLSURF
RF, 
*SUBSTRUCTURE PATH,LEAVE
*END STEP