*HEADING
 SUBSTRUCTURE ANALYSIS - SUBSTRUCTURE GENERATED EARLIER
*NODE
*NODE,SYSTEM=C
*NGEN,LINE=C
*NODE,NSET=FIX
*ELEMENT,TYPE=R2D2
*ELGEN,ELSET=RIGID
*nset,nset=NBOTSUR,GENERATE
*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=Z1001,FILE=hertzcontact_3d_sub_only,ELSET=SUPER
*SUBSTRUCTURE PROPERTY,ELSET=SUPER
*BOUNDARY
 300,1,6
 201,1
*RESTART,WRITE
*STEP
*STATIC,DIRECT
 0.2 , 1.
*BOUNDARY
 1,2,2,-.4
*CONTACT PRINT ,SLAVE=BOTSUR,MASTER=RIGSUR
*CONTACT FILE,SLAVE=BOTSUR,MASTER=RIGSUR
*NODE PRINT
U,
*OUTPUT,FIELD,VAR=PRESELECT
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT,FREQUENCY=5
 S ,
*NODE PRINT
 U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT,ELSET=TOP,FREQUENCY=5
E,
*EL FILE,ELSET=BOT,FREQUENCY=5
S,E
*NODE PRINT,NSET=CYLSURF
RF,
*SUBSTRUCTURE PATH,LEAVE
*ELEMENT MATRIX OUTPUT,ELSET=SUPER,STIFFNESS=YES,FREQUENCY=999
*END STEP