*HEADING
 Substructure generation file for hertzcontact_2d_coupk_substr.inp 
*NODE
   1, 0.,0.
   5, 2.,0.
   9, 5.5,0.
  13, 10.,0.
  41, 0.,-3.
  45, 2.,-3.
  49, 5.2,-3.
  81, 0.,-5.5
  85, 2.,-5.5
  89, 4.,-5.5
 121, 0.,-7.5
 129, 2.,-7.5
 133, 3.5,-7.5
 101, 0.,-6.5
 103, 1.,-6.5
 105, 2.,-6.5
 107, 3.75,-6.5
 201, 0.,-10.
*NODE,SYSTEM=C
 33,10.,-15.
 53,10.,-30.
 73,10.,-37.5
 93,10.,-45.
113,10.,-55.
213,10.,-65.
209,10.,-78.
*NGEN,LINE=C
 209,213
 201,209
*NGEN
 1,5,2
 5,9,2
 9,13,2
 41,45,2
 45,49,2
 49,53,2
 81,85
 85,89,2
 89,93,2
 121,129
 129,133
 133,213,40
 131,211,40
 1,41,20
 5,45,20
 9,49,20
 41,81,20
 45,85,20
 49,89,20
 81,121,20
 121,201,20
 123,203,20
 125,205,20
 127,207,20
 129,209,20
 161,163
 163,165
 165,167
 167,169
*NSET,NSET=NY
 1,21,41,61,81,101,121,141,161,181,201
*ELEMENT,TYPE=CPE8
 1, 1,41,45,5,21,43,25,3
 7, 81,121,125,83,101,123,103,82
 8, 83,125,129,85,103,127,105,84
 9, 85,129,133,89,105,131,107,87
 10, 89,133,213,93,107,173,113,91
 11, 121,161,163,123,141,162,143,122
 19, 129,209,211,131,169,210,171,130
 20, 131,211,213,133,171,212,173,132
*ELGEN
 1,3,4,1
 1,2,40,3
 4,3,4,1
 11,4,2,1
 11,2,40,4
 15,4,2,1
*ELSET,ELSET=TOP
 1,2,3,4,5,6,7,8,9,10,11,12,13,14,19,20
*ELSET,ELSET=BOT
 15,16,17,18
*ELSET,ELSET=CYLINDER
TOP,BOT
*SURFACE, NAME=TOPSUR, TYPE=NODE
 3,1.0
 5,1.0
 7,1.0
 9,1.0
 11,1.0
 13,1.0
*COUPLING, CONSTRAINT NAME=COUPLINGDEF, REF NODE=1, SURFACE=TOPSUR
*KINEMATIC
2
*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
201,209,1
**
**
*STEP
*SUBSTRUCTURE GENERATE,TYPE=Z1001
*RETAINED NODAL DOFS
CYLSURF,2
1,2,6
*BOUNDARY
NY,1
*END STEP