The modal characteristics of constrained multibody systems undergoing constant accelerated motions are investigated in this paper. Relative coordinates are employed to derive the equations of motion, which are generally nonlinear in terms of the coordinates. The dynamic equilibrium position of a constrained multibody system needs to be obtained from the nonlinear equations of motion, which are then linearized at the dynamic equilibrium position. The mass and the stiffness matrices for the modal analysis can be obtained from the linearized equations of motion. To verify the effectiveness and the accuracy of the proposed method, two numerical examples are solved and the results obtained by using the proposed method are compared with those obtained by analytical and other numerical methods. The proposed method is found to be accurate as well as effective in predicting the modal characteristics of constrained multibody systems undergoing constant accelerated motions.