*Heading
 C3D8 Elements; Static Analysis
*Node
      1,          10.,           0.,           5.
      2,     6.666667,           0.,           5.
      3,     3.333333,           0.,           5.
      4,           0.,           0.,           5.
      5,          10.,           0.,          2.5
      6,     6.666667,           0.,          2.5
      7,     3.333333,           0.,          2.5
      8,           0.,           0.,          2.5
      9,          10.,           0.,           0.
     10,     6.666667,           0.,           0.
     11,     3.333333,           0.,           0.
     12,           0.,           0.,           0.
     13,          10.,          2.5,           5.
     14,     6.666667,          2.5,           5.
     15,     3.333333,          2.5,           5.
     16,           0.,          2.5,           5.
     17,          10.,          2.5,          2.5
     18,     6.666667,          2.5,          2.5
     19,     3.333333,          2.5,          2.5
     20,           0.,          2.5,          2.5
     21,          10.,          2.5,           0.
     22,     6.666667,          2.5,           0.
     23,     3.333333,          2.5,           0.
     24,           0.,          2.5,           0.
     25,          10.,           5.,           5.
     26,     6.666667,           5.,           5.
     27,     3.333333,           5.,           5.
     28,           0.,           5.,           5.
     29,          10.,           5.,          2.5
     30,     6.666667,           5.,          2.5
     31,     3.333333,           5.,          2.5
     32,           0.,           5.,          2.5
     33,          10.,           5.,           0.
     34,     6.666667,           5.,           0.
     35,     3.333333,           5.,           0.
     36,           0.,           5.,           0.
*Element, type=C3D8
 1, 13, 14, 18, 17,  1,  2,  6,  5
 2, 14, 15, 19, 18,  2,  3,  7,  6
 3, 15, 16, 20, 19,  3,  4,  8,  7
 4, 17, 18, 22, 21,  5,  6, 10,  9
 5, 18, 19, 23, 22,  6,  7, 11, 10
 6, 19, 20, 24, 23,  7,  8, 12, 11
 7, 25, 26, 30, 29, 13, 14, 18, 17
 8, 26, 27, 31, 30, 14, 15, 19, 18
 9, 27, 28, 32, 31, 15, 16, 20, 19
10, 29, 30, 34, 33, 17, 18, 22, 21
11, 30, 31, 35, 34, 18, 19, 23, 22
12, 31, 32, 36, 35, 19, 20, 24, 23
*Elset, elset=I1, generate
  1,  12,   1
** Section: Secn
*Solid Section, elset=I1, material=Elas
1.,
** 
** MATERIALS
** 
*Material, name=Elas
*Density
25.,
*Elastic
 1e+08,0.
*step
*frequency,eigen=lanc, SIM=NO
81,
*boundary
1,1,3
4,1,3
9,1,3
12,1,3
21,1,3
25,1,3
28,1,3
33,1,3
36,1,3
*output, field,freq=0
*end step
*STEP
*substructure generate, type=z1, overwrite, 
recovery matrix=yes, mass matrix=yes,
libr=substr_contact_solid3d_gen
*RETAINED NODAL DOFS
1,1,3
4,1,3
9,1,3
12,1,3
21,1,3
25,1,3
28,1,3
33,1,3
36,1,3
*select eigenmodes,generate
1,81,1
*substructure matrix output, stiffness=yes, mass=yes, recovery=yes,
output file=user defined, file name=substr_contact_solid3d_gen
*END STEP