*include,input=crank_conn_gr_ref_nodes.inp
*include,input=crank_shaft_nodes.inp
*include,input=crank_shaft_elts.inp 
*include,input=crank_materials.inp
*** coupling nodes: 
***    5 - bore center on face adjacent to CONROD
***    6 - shaft axis to be driven by connector motion
*include,input=crank_shaft_coupling.inp
*********************************************
*********************************************
*step
*frequency,eigen=lanc, SIM=NO
6,
*boundary
5,1,6
6,1,6
385,1,3
*output,field,freq=0
*output,hist,freq=0
*end step
*step
*substructure generate,type=z3,overwrite,mass=yes,libr=crank_shaft,grav=yes
**recovery=yes
***
*** the two coupling nodes above and an arbitrary node for rigid body calculations 
*RETAINED NODAL DOFS,SORTED=NO
5,1,6
6,1,6
385,1,3
*SELECT EIGENMODES
1,2,3,4,5,6
***
*substructure matrix output, stiffness=yes,mass=yes,file=crank_shaft, grav=yes
**recovery=yes
*end step