*HEADING
 Substructure generation file for hertzcontact_2d_substr_sload.inp
*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
*EQUATION
*MPC
QUADRATIC,900000082,900000081,900000083,900000085
QUADRATIC,900000084,900000081,900000083,900000085
QUADRATIC,900000122,900000121,900000123,900000125
QUADRATIC,900000124,900000121,900000123,900000125
QUADRATIC,900000126,900000125,900000127,900000129
QUADRATIC,900000128,900000125,900000127,900000129
QUADRATIC,900000130,900000129,900000131,900000133
QUADRATIC,900000132,900000129,900000131,900000133
QUADRATIC,900000149,900000129,900000169,900000209
QUADRATIC,900000189,900000129,900000169,900000209
*SOLID SECTION,ELSET=CYLINDER,MATERIAL=A1
*MATERIAL,NAME=A1
*ELASTIC
 30.E6,0.3
*NSET,NSET=CYLSURF,GENERATE
**
**
*BOUNDARY
NY,1
*STEP
*SUBSTRUCTURE GENERATE,TYPE=Z1001,OVERWRITE
*RETAINED NODAL DOFS
CYLSURF,
*END STEP