*HEADING
 TEST USING DISTRIBUTING COUPLING ELEMENT WITHIN
 THE SUBSTRUCTURE;
 NEED hertzcontact_2d_coupd_substrgen.sup
*NODE
*NODE,SYSTEM=C
*NGEN,LINE=C
*NODE
*NODE,NSET=FIX
*ELEMENT,TYPE=R2D2
*ELGEN,ELSET=RIGID
*RIGID BODY,ELSET=RIGID,REF NODE=300
*NSET,NSET=BOTSURNODE,GENERATE
*SURFACE,TYPE=NODE,NAME=BOTSUR
BOTSURNODE,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=SUBSTRUCTURE,
FILE=hertzcontact_2d_coupd_substrgen
*SUBSTRUCTURE PROPERTY,ELSET=SUBSTRUCTURE,POSITION TOL=0.
*NORMAL
RIGID,FIX,,1.
*BOUNDARY
 300,1,2
 300,6
1001,6
*RESTART,WRITE
*STEP
*STATIC,DIRECT
 1.
*BOUNDARY
1001,2,2,-.4
*CONTACT PRINT ,SLAVE=BOTSUR,MASTER=RIGSUR
*CONTACT FILE ,SLAVE=BOTSUR,MASTER=RIGSUR
*PRINT,CONTACT=YES,SOLVE=YES
*NODE PRINT
U,
*OUTPUT,VAR=PRESELECT,FIELD
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT
 S ,
*NODE PRINT
 U,
*SUBSTRUCTURE PATH,LEAVE
*SUBSTRUCTURE PATH,ENTER ELEMENT=99
*OUTPUT,VAR=PRESELECT,FIELD
*EL PRINT,ELSET=TOP
E,
*EL FILE,ELSET=BOT
S,E
*NODE PRINT,NSET=CYLSURF
RF
*SUBSTRUCTURE PATH,LEAVE
*END STEP