*HEADING
 HERTZ CONTACT PROBLEM 3D WITH SUBSTRUCTURING         
 ELEM: C3D10M    INTERFACE: SMALL SLIDING
Requires substructure generation file 
hertzcontact_c3d10m_gen1.inp
*NODE
       1,          0.,        -10.
      24,    0.188987,    -9.99821
      40,    0.377907,    -9.99286
      64,    0.613052,    -9.98119
      79,    0.847859,    -9.96399
     103,     1.13969,    -9.93484
     118,     1.43054,    -9.89715
     143,     1.79112,    -9.83829
     157,      2.1493,    -9.76629
     181,     2.59144,    -9.65839
     451,          0.,        -10.,        -0.5
     471,    0.188987,    -9.99821,        -0.5
     481,    0.377907,    -9.99286,        -0.5
     501,    0.613053,    -9.98119,        -0.5
     511,    0.847859,    -9.96399,        -0.5
     531,     1.13969,    -9.93484,        -0.5
     541,     1.43054,    -9.89715,        -0.5
     561,     1.79112,    -9.83829,        -0.5
     571,      2.1493,    -9.76629,        -0.5
     591,     2.59144,    -9.65839,        -0.5
     802,          0.,        -10.,         -1.
     825,    0.188987,    -9.99821,         -1.
     841,    0.377907,    -9.99286,         -1.
     864,    0.613052,    -9.98119,         -1.
     880,    0.847859,    -9.96399,         -1.
     905,     1.13969,    -9.93484,         -1.
     919,     1.43054,    -9.89715,         -1.
     943,     1.79112,    -9.83829,         -1.
     958,      2.1493,    -9.76629,         -1.
     983,     2.59144,    -9.65839,         -1.
     430,         10.,  4.37114E-7
*NODE,NSET=RIGNODE
5000,-.5,-10,-1.5
5001,-.5,-10,.5
5002,2.5,-10,-1.5
5003,2.5,-10,.5
10000,1.,-11,1
*ELEMENT,TYPE=Z1003,ELSET=SUPER, FILE=hertzcontact_c3d10m_gen1
2999,1,24,40,64,79,103,118,143,157,181,430,
451,471,481,501,511,531,541,561,571,591,
802,825,841,864,880,905,919,943,958,983
*SUBSTRUCTURE PROPERTY,ELSET=SUPER,POSITION TOL=0.
**
** Contact definition
**
*ELEMENT,TYPE=R3D4,ELSET=ERIGID
10000,5000,5001,5003,5002
*RIGID BODY,ELSET=ERIGID,REF NODE=10000
*surface,NAME=RIGID
ERIGID,SPOS
*nset,nset=NCYLBOT
1,24,40,64,79,103,118,143,157,181,
451,471,481,501,511,531,541,561,571,591,
802,825,841,864,880,905,919,943,958,983
*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,,      
     118, 3,,      
     143, 3,,      
     157, 3,,      
     181, 3,,      
     451, 1,,      
     802, 1,,      
     802, 3,,      
     825, 3,,      
     841, 3,,      
     864, 3,,      
     880, 3,,          
     905, 3,,          
     919, 3,,          
     943, 3,,          
     958, 3,,          
     983, 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,MAIN=RIGID,SECONDARY=CYLBOT,FREQUENCY=5
*CONTACT PRINT,MAIN=RIGID,SECONDARY=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