*HEADING
 HERTZ CONTACT PROBLEM 3D:                    
 ELEM: C3D20R, C3D20, R3D4    INTERFACE: SMALL SLIDING	
Requires substructure generation file hertzcontact_3d_gen1.inp	
*RESTART,WRITE,OVERLAY,FREQUENCY=10
*NODE
 201, 0.,-10.
1001,.0,.0,.5
1201,.0,-10.,.5
2001,.0,.0,1.
2201,.0,-10.,1.
*NODE,SYSTEM=C
209,10.,-78.
1209,10.,-78.,.5
2209,10.,-78.,1.
*NGEN,LINE=C
 201,209
1201,1209,1,1001
2201,2209,1,2001
*NSET,NSET=NCONT, GENERATE
201, 209
1201, 1209, 2
2201, 2209
*NSET,NSET=ZSYM, GENERATE
201,209,1
2201,2209,1
*NSET,NSET=NX
201,1201,2201
*NODE,NSET=RIGNODE
5000,-.5,-10,-.5
5001,-.5,-10,1.5
5002,2.5,-10,-.5
5003,2.5,-10,1.5
10000,1.,-11,1
*ELEMENT,TYPE=Z1001,ELSET=SUPER, FILE=hertzcontact_3d_gen1
999,1, 201, 202, 203, 204, 205, 206, 207, 208, 209,1201,1203,1205,1207,1209,
2201,2202,2203,2204,2205,2206,2207,2208,2209
*SUBSTRUCTURE PROPERTY,ELSET=SUPER
*ELEMENT,TYPE=R3D4,ELSET=ERIGID
100,5000,5001,5003,5002
*RIGID BODY,ELSET=ERIGID,REF NODE=10000
**
** Contact definition
** 
*surface,NAME=RIGID
ERIGID,SPOS
*surface,type=node,NAME=SECONDARYS
NCONT, 
*CONTACT PAIR, INTERACTION=HERTZ,SMALL SLIDING
SECONDARYS,RIGID
*SURFACE INTERACTION,NAME=HERTZ
*STEP,INC=10
*STATIC,DIRECT
.1,1.
*BOUNDARY
 10000,1,6,0
 1,2,,-0.4
 ZSYM,3,3
 NX,1,1
*PRINT,CONTACT=YES,SOLVE=YES
*CONTACT PRINT,MAIN=RIGID,FREQUENCY=2
*CONTACT FILE,MAIN=RIGID,FREQUENCY=10
*NODE FILE,FREQUENCY=10
U,
RF, 
*OUTPUT,FIELD,VAR=PRESELECT
*SUBSTRUCTURE PATH,ENTER ELEMENT=999
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=10
*SUBSTRUCTURE PATH,LEAVE
*END STEP