In this Study, the dynamic equation for vibration analysis of mechanical systems with kinematic constraints is derived. This equations are derived in terms of small displacements of Cartesian coordinates, and are applied to compute the dynamic response and the natural modes of the suspension system of a vehicle. The results are validated through the comparison with the results from conventional nonlinear dynamic analysis and modal test.