*include,input=crank_conn_gr_ref_nodes.inp
*include,input=crank_conrod_nodes.inp
*include,input=crank_conrod_elts.inp 
*include,input=crank_materials.inp
*** coupling nodes: 
***    3 - bore center on face adjacent to shafts
***    4 - bore center on face adjacent to head
*include,input=crank_conrod_coupling.inp
*********************************************
*********************************************
*step
*frequency,eigen=lanc
6,
*boundary
3,1,6
4,1,6
16145,1,3
*output,field,freq=0
*output,hist,freq=0
*end step
*step
*substructure generate,type=z2,overwrite,mass=yes,libr=crank_conrod,grav=yes
**recovery=yes
***
*** the two coupling nodes above and an arbitrary node for rigid body calculations 
*RETAINED NODAL DOFS,SORTED=NO
3,1,6
4,1,6
16145,1,3
*RETAINED EIGENMODES
1,2,3,4,5,6
*substructure matrix output, stiffness=yes, mass=yes,file=crank_conrod,grav=yes
**recovery=yes
*end step