*HEADING
 SUBSTRUCTURE ANALYSIS - MPC ONLY IN SUBSTRUCTURE          1-2-3-E
Requires superelment generation file 
hertzcontact_2d_xnodes_gen1.inp
*NODE
   1, 0.,0.
  41, 0.,-3.
  81, 0.,-5.5
 101, 0.,-6.5
 121, 0.,-7.5
 201, 0.,-10.
*NGEN
 1,41,20
41,81,20
81,121,20
121,201,20
*NODE
 201, 0.,-10.
 300, 0.,-11.
*NODE,SYSTEM=C
 209,10.,-78.
*NGEN,LINE=C
 201,209
*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.
*NSET,NSET=NY
 1,21,41,61,81,101,121,141,161,181,201
*ELEMENT,TYPE=R2D2
101,301,302
*ELGEN,ELSET=RIGID
101,8,1,1
*nset,nset=NBOTSUR,GEN
201,209,1
*RIGID BODY,ELSET=RIGID,REFNODE=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=Z1001,ELSET=ZUPER, FILE=hertzcontact_2d_xnodes_gen1
 99,  1,  3,  5,  7,  9, 11, 13, 21, 41, 61, 81,101,121,141,161,
181,201,202,203,204,205,206,207,208,209
*SUBSTRUCTURE PROPERTY,ELSET=ZUPER,POSITION TO=0.0
*BOUNDARY
NY,1
 300,1,6
*EQUATION
 2, 
 3,2,-1.,1,2,1.
 2, 
 5,2,-1.,1,2,1.
 2, 
 7,2,-1.,1,2,1.
 2, 
 9,2,-1.,1,2,1.
 2, 
11,2,-1.,1,2,1.
 2, 
13,2,-1.,1,2,1.
*STEP
*STATIC,DIRECT
 1., 
*BOUNDARY
 1,2,2,-.4
*CONTACT PRINT ,SECONDARY=BOTSUR,MAIN=RIGSUR
*CONTACT FILE ,SECONDARY=BOTSUR,MAIN=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
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT,ELSET=TOP
E,
*EL FILE,ELSET=BOT
S,E
*NODE PRINT,NSET=CYLSURF
RF, 
*SUBSTRUCTURE PATH,LEAVE
*END STEP