• The multi-disk rub-impact analysis of a two-spool aero-engine dual-rotor model is performed using the approximate time variational method. • The significance of multi-disk rub-impact studies compared to the single-disk rub-impact problems are discussed based on the response and stability analysis. • Limit point and Neimark-Sacker bifurcations are observed in the response indicating sudden jump and origin of quasi-periodic motion respectively. • The influences of rubbing parameters on the dynamic characteristics of the dual-rotor model are critically evaluated. • The onset of quasi-periodic motion is happening early as the coefficient of friction and contact stiffness are increased. [Display omitted] The main aim of this paper is to propose a numerical procedure for capturing the nonlinear dynamic characteristics of a two-spool aero-engine rotor system undergoing multi-disk rub-impact. In aero-engines, the possibility for the multi-disk rub-impact is high during the fan blade-out (FBO) event and subsequent windmilling action. It intensifies the nonlinear effects on the rotor vibrations and leads to certain undesired circumstances in the engine. The dual-rotor model consists of multi-stage compressors and single-stage turbines that undergo rubbing whenever their deflection exceeds the clearance. The dynamic model of the dual-rotor system is constructed using the tapered Timoshenko beam elements, rigid disks and rolling contact bearings. A proper model reduction technique based on component mode synthesis coupled with the Craig-Bampton substructuring is utilized to reduce the size of the finite element model. A semi-analytic technique called the approximate time variational method is employed to investigate the steady-state response of the system under multi-disk rub impact. Based on the proposed method, the response characteristics of the model are obtained and are verified with the numerical integration results. Compared to the single-disk rub-impact, the nonlinearities are intensified and significant variations are observed in the response characteristics and stability of the system. Period-5, quasi-periodic, and dry friction backward whirl motions are observed in the response for different values of the system parameters. During quasi-periodic motion, some unknown fractional components such as 0.716 ω 1 , 0.766 ω 1 , 0.916 ω 1 and 0.964 ω 1 are appeared in the response. Moreover, the dry friction backward whirl happened with a very large amplitude and it contains a superharmonic frequency component in the response. [ABSTRACT FROM AUTHOR]