*HEADING
 SUBSTRUCTURE ANALYSIS
 WITH A DISTRIBUTING COUPLING CONSTRAINT IN THE SUBSTRUCTURE
*PREPRINT,MODEL=YES,HISTORY=YES
*NODE
*NODE,SYSTEM=C
*NGEN,LINE=C
*NGEN
*NSET,NSET=NY
*ELEMENT,TYPE=CPE8
*ELGEN
*ELSET,ELSET=TOP
*ELSET,ELSET=BOT
*ELSET,ELSET=CYLINDER
TOP,BOT
*SURFACE, NAME=TOPSURF
 1, S4
 2, S4
 3, S4
*COUPLING, CONSTRAINT NAME=COUPLINGDEF, REF NODE=1001, SURFACE=TOPSURF
*DISTRIBUTING
*MPC
QUADRATIC,82,81,83,85
QUADRATIC,84,81,83,85
QUADRATIC,122,121,123,125
QUADRATIC,124,121,123,125
QUADRATIC,126,125,127,129
QUADRATIC,128,125,127,129
QUADRATIC,130,129,131,133
QUADRATIC,132,129,131,133
QUADRATIC,149,129,169,209
QUADRATIC,189,129,169,209
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=A1
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
*NSET,NSET=CYLSURF,GENERATE
*STEP
*SUBSTRUCTURE GENERATE,TYPE=Z1001,OVERWRITE
*RETAINED NODAL DOFS
CYLSURF,2
*BOUNDARY
NY,1
*END STEP