By specializing the Hamiltonian of Nose's extended system in rigid multibody systems, we develop rigorous motion equations for rigid multibody molecular dynamics in canonical ensemble. The equations are proved to generate correct canonical equilibrium distributions. The corresponding reversible and symplectic integrator is also given in this paper.