A dynamic model is reported herein for the study of vibrations of deep groove ball bearings having single and multiple defects on surfaces of inner and outer races. Masses of shaft, housing, races, and balls are considered in the modeling. The coupled solution of governing equations of motions is obtained using Runge–Kutta method. The model provides the vibrations of shaft, balls, and housing in time and frequency domains. Computed results from the model are validated with experimental results, which are generated using healthy and defective deep groove ball bearings. Characteristic defect frequencies and its harmonics are broadly investigated using both theoretical and experimental results. Comparison of vibration spectra for the cases having single and two defects on races reveals relatively higher velocity amplitudes with two defects. Good correlations between theoretical and experimental results are observed. Authors believe that this dynamic model can be used with confidence for the study and prediction of vibrations of healthy and defective deep groove ball bearings.