Stability of a Rigid Body With an Oscillating Particle: An Application of MACSYMA

This problem is a generalization of the classical problem of the stability of a spinning rigid body. We obtain the stability chart by using: (i) the computer algebra system MACSYMA in conjunction with a perturbation method, and (ii) numerical integration based on Floquet theory. We show that the form of the stability chart is different for each of the three cases in which the spin axis is the minimum, maximum, or middle principal moment of inertia axis. In particular, a rotation with arbitrarily small angular velocity about the maximum moment of inertia axis can be made unstable by appropriately choosing the model parameters. In contrast, a rotation about the minimum moment of inertia axis is always stable for a sufficiently small angular velocity. The MACSYMA program, which we used to obtain the transition curves, is included in the Appendix.