Generalized Euler equations for linked rigid bodies

We derive the equations of motion for linked rigid bodies from Lagrange mechanics and from Gauss’s principle of least constraint. The rotational motion of the subunits is described in terms of quaternion parameters and angular velocities. Different types of joints can be incorporated via axis constraints for the angular velocities. The resulting equations of motion are generalizations of the Euler equations of motion for a single rotor.