This paper presents a numerical model for investigating the structural dynamic response of an unbalanced rotor system supported on deep groove ball bearings. The aim of this work is to develop a numerical model for investigating the effect of the number of balls on the dynamic characteristics of the rotor ball bearing system. The fourth-order Runge-Kutta numerical integration technique has been applied. The results are presented in the form of time displacement responses and frequency spectra. The analysis demonstrates that the model can be used as a tool for predicting the nonlinear dynamic behavior of the rotor ball bearing system under different operating conditions. Moreover, the study may contribute to a further understanding of the nonlinear dynamics of rotor bearing systems.