*HEADING
 TEST OF SUBSTRUCTURE MATRIX INPUT 
*NODE
 201, 0.,-10.
 300, 0.,-11.
*NODE,SYSTEM=C
 209,10.,-78.
*NGEN,LINE=C
 201,209
*NODE,NSET=FIX
 301, 0.     ,-10.
 302, 0.26177,-10.
 303, 0.52336,-10.
 304, 0.78459,-10.
 305, 1.04528,-10.
 306, 1.30526,-10.
 307, 1.56434,-10.
 308, 1.82236,-10.
 309, 2.07912,-10.
*ELEMENT,TYPE=R2D2
101,301,302
*ELGEN,ELSET=RIGID
101,8,1,1
*nset,nset=NBOTSUR,GENERATE
201,209,1
*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, 
**
**
*USER ELEMENT,TYPE=U1001,FILE=hertzcontact_3d_sub_library,OLD ELEMENT=99,STEP=1,INCREMENT=5
*ELEMENT,TYPE=U1001,ELSET=USEREL
 99,  1,201,202,203,204,205,206,207,208,209
*UEL PROPERTY,ELSET=USEREL
*BOUNDARY
 300,1,6
 201,1
*RESTART,WRITE
*STEP
*STATIC,DIRECT
 0.2 , 1.
*BOUNDARY
 1,2,2,-.4
*CONTACT PRINT ,SECONDARY=BOTSUR,MAIN=RIGSUR
*CONTACT FILE,SECONDARY=BOTSUR,MAIN=RIGSUR
*OUTPUT,FIELD
*CONTACT OUTPUT,VARIABLE=PRESELECT,SECONDARY=BOTSUR,MAIN=RIGSUR
*NODE PRINT
U,
*END STEP