Blade-rotor systems frequently encounter the problem of blade-to-case rubbing, which affects their safety and stability. Numerical simulation can be used to predict the steady-state response of these systems. However, such simulation is frequently computationally expensive because of the high dimensions of the dynamic model of a blade-rotor system. To overcome this problem, a new method that combines the receptance-based dimension-reduction approach with the incremental harmonic balance (IHB) method is presented in this study. First, a dynamic model of a blade-rotor system is developed using the finite element method, and the number of dimensions of the model is reduced by the receptance method. Subsequently, the steady-state response is obtained by the improved IHB method to conveniently manage the large number of super-harmonic components of the local rubbing system. Finally, the precision and efficiency of the proposed method is verified by comparing its results with those obtained by the Newmark-b method. The proposed method is found to be efficient in analyzing local rubbing blade-rotor systems with high dimensions, local nonlinearities, and rich super-harmonics.