*include,input=crank_conn_gr_ref_nodes.inp
*include,input=crank_head_nodes.inp
*include,input=crank_head_elts.inp 
*include,input=crank_materials.inp
*** coupling nodes: 
***    1 - bore center on face adjacent to CONROD
***    2 - center of cylider top face
*include,input=crank_head_coupling.inp
*********************************************
*********************************************
*step
*frequency,eigen=lanc
6,
*boundary
1,1,6
2,1,6
5186,1,3
*output,field,freq=0
*output,hist,freq=0
*end step
*step
*substructure generate,type=z1,overwrite,mass=yes,libr=crank_head,grav=yes
**recovery=yes
***
*** the two coupling nodes above and an arbitrary node for rigid body calculations 
*RETAINED NODAL DOFS,SORTED=NO
1,1,6
2,1,6
5186,1,3
*SELECT EIGENMODES
1,2,3,4,5,6
*substructure matrix output, stiffness=yes, mass=yes,file=crank_head,grav=yes
**recovery=yes
*end step