*HEADING
 Substructure generation file for hertzcontact_2d_unsorted.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
*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
*EQUATION
 2, 
 3,2,-1.,1,2,1.
 2, 
 5,2,-1.,1,2,1.
 2, 
 7,2,-1.,1,2,1.
 2, 
 9,2,-1.,1,2,1.
 2, 
11,2,-1.,1,2,1.
 2, 
13,2,-1.,1,2,1.
*NSET,UNSORTED,NSET=MPC1
82,84,122,124,126,128,130,132,149,189
*NSET,UNSORTED,NSET=MPC2
81,81,121,121,125,125,129,129,129,129
*NSET,UNSORTED,NSET=MPC3
83,83,123,123,127,127,131,131,169,169
*NSET,UNSORTED,NSET=MPC4
85,85,125,125,129,129,133,133,209,209
*MPC
QUADRATIC,MPC1,MPC2,MPC3,MPC4
*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,OVERWRITE
*RETAINED NODAL DOFS,SORTED=NO
CYLSURF,
1,2
*BOUNDARY
NY,1
*END STEP