In this paper, the nonlinear dynamics of a novel model based on optically pumped spin-polarized vertical-cavity surface-emitting lasers (spin-VCSELs) with optical feedback is investigated numerically. Due to optical feedback being the external disturbance component, the complex nonlinear dynamical behaviors can be enhanced and the regions of different nonlinear dynamics in size can be extended with appropriate parameters of spin-VCSELs. According to the equations of the modified spin-flip model (SFM), the comparison of bifurcation diagrams is first presented for the clear presentation of different routes to chaos. Meanwhile, numerous bifurcation diagrams in color are illustrated to demonstrate the rich dynamical regimes intuitively, and the crucial effects of optical feedback strength, feedback delay, linewidth enhancement factor, and spin-flip relaxation rate on the region evolvement of complex dynamics of the proposed model are revealed to investigate the dependence of dynamical behaviors on external and internal parameters when the optical feedback scheme is introduced. These parameters play a remarkable role in enhancing the mechanism of complex dynamic oscillations. Furthermore, utilizing combination with time series, power spectra, and phase portraits, the various dynamical behaviors observed in the bifurcation diagram are simulated numerically. Correspondingly, the powerful measure 0–1 test is employed to distinguish between chaos and non-chaos.