The dynamics of a rotor system with axial contact/rub events between the disks and stator are investigated by numerical simulations. The formula for determining the contact/rub points, axial contact forces and dry friction forces are deduced. To account for their influence, the axial contact forces are substituted by equivalent forces acting at the disk centers, based on the equivalent moment rule. One-parametric model is used to estimate the contact-induced dry friction forces. The coupled equations of lateral and torsional motions of rotor and the lateral motion of disk are then established. Numerical simulations are carried out to reveal the lateral and torsional vibrations for both twodisk contact/rubs with different axial clearances, and one disk contact/rubs. Bifurcation diagrams, orbits, phase portraits, amplitude-frequency spectra and Poincaré maps are adopted to demonstrate the dynamical behaviors of the system. The results show that though both the lateral and torsional vibrations can reflect the influences of contact/rubs on rotor dynamics, the spectrum analyses of the torsional vibrations are more suitable to determine straight the extent of their effect.