*include,input=crank_conn_gr_ref_nodes.inp
*include,input=crank_cam_nodes.inp
*include,input=crank_cam_elts.inp 
*include,input=crank_materials.inp
*** coupling nodes: 
***    19999 - the reference node of the rigid surface
***    9 - cam center (same axial position as 19999)
*include,input=crank_cam_coupling.inp
*********************************************
*********************************************
*step
*frequency,eigen=lanc
6,
*boundary
19999,1,6
9,1,6
16992,1,3
*output,field,freq=0
*output,hist,freq=0
*end step
*step
*substructure generate,type=z5,overwrite,mass=yes,libr=crank_cam,grav=yes
**recovery=yes
***
*** the two coupling nodes above and an arbitrary node for rigid body calculations 
*RETAINED NODAL DOFS,SORTED=NO
19999,1,6
9,1,6
16992,1,3
*RETAINED EIGENMODES
1,2,3,4,5,6
*substructure matrix output, stiffness=yes, mass=yes,file=crank_cam,grav=yes
**recovery=yes
*end step