*HEADING
 SUBSTRUCTURE ANALYSIS - USING THE SLOAD OPTION       1-2-3-F
 Requires substructure generation file
 hertzcontact_2d_sload_gen1.inp
*NODE
*NODE
*NODE,SYSTEM=C
*NGEN,LINE=C
*NODE,NSET=FIX
*ELEMENT,TYPE=R2D2
*ELGEN,ELSET=RIGID
*nset,nset=NBOTSUR,GEN
*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=SUPER, FILE=hertzcontact_2d_sload_gen1
*SUBSTRUCTURE PROPERTY,ELSET=SUPER
*BOUNDARY
 300,1,6
 201,1
*STEP
*STATIC,DIRECT
 1.,
*BOUNDARY
900000001,2,2,-.4
*CONTACT PRINT ,SLAVE=BOTSUR,MASTER=RIGSUR
cstress,cdisp
carea,
CFN,CFS
CMN,CMS
CFT,CMT
XN,XS,XT
*OUTPUT, FIELD
*CONTACT OUTPUT, SLAVE=BOTSUR, MASTER=RIGSUR
cstress,cdisp
** *CONTACT FILE ,SLAVE=BOTSUR,MASTER=RIGSUR
** cstress,cdisp
** carea,
** CFN,CFS
** CMN,CMS
** CFT,CMT
** XN,XS,XT
*PRINT,CONTACT=YES,SOLVE=YES
*NODE PRINT
U,
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*EL PRINT
 S ,
*NODE PRINT
 U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*EL PRINT,ELSET=TOP
E,
*OUTPUT, FIELD
*ELEMENT OUTPUT, ELSET=BOT
 S, E
** *EL FILE,ELSET=BOT
** S,E
*NODE PRINT,NSET=CYLSURF
RF,
*SUBSTRUCTURE PATH,LEAVE
*END STEP