*HEADING
 HERTZ CONTACT PROBLEM 3D WITH SUBSTRUCTURING
 ELEM: C3D10M    INTERFACE: SMALL SLIDING
Requires substructure generation file
hertzcontact_c3d10m_gen1.inp
*NODE
*NODE,NSET=RIGNODE
*ELEMENT,TYPE=Z1003,ELSET=SUPER, FILE=hertzcontact_c3d10m_gen1
*SUBSTRUCTURE PROPERTY,ELSET=SUPER,POSITION TOL=0.
**
** Contact definition
**
*ELEMENT,TYPE=R3D4,ELSET=ERIGID
*RIGID BODY,ELSET=ERIGID,REF NODE=10000
*surface,NAME=RIGID
ERIGID,SPOS
*nset,nset=NCYLBOT
*surface,type=node,NAME=CYLBOT
NCYLBOT,1.0
*CONTACT PAIR, INTERACTION=HERTZ,SMALL SLIDING
CYLBOT,RIGID
*SURFACE INTERACTION,NAME=HERTZ
*BOUNDARY
10000, 1,6
       1, 1,,
       1, 3,,
      24, 3,,
      40, 3,,
      64, 3,,
      79, 3,,
     103, 3,,
*STEP,INC=10
*STATIC,DIRECT
.1,1.
*RESTART,WRITE
*BOUNDARY
 430,2,,-0.4
*PRINT,CONTACT=YES
*EL PRINT,FREQUENCY=2
S11,E11
*NODE PRINT,FREQUENCY=10
U,
RF,
*NODE FILE,FREQUENCY=10
RF,
*CONTACT FILE,MASTER=RIGID,SLAVE=CYLBOT,FREQUENCY=5
*CONTACT PRINT,MASTER=RIGID,SLAVE=CYLBOT,FREQ=5
*OUTPUT,VAR=PRESELECT,FIELD
*SUBSTRUCTURE PATH,ENTER ELEMENT=2999
*OUTPUT,VAR=PRESELECT,FIELD
*NODE PRINT,NSET=CYLSURF
RF,
*SUBSTRUCTURE PATH,LEAVE
*END STEP