<p>The variation in rolling bearing clearance exhibits significant sensitivity to the dynamics behavior of the system and it is often given as a constant value in the current research. However, the reconstruction phenomenon of the clearance constrained by interference assembly could significantly affect the bearing contact characteristics, the resultant nonlinear dynamics behaviors of the system have not been further investigated, especially when the High-speed centrifugal effect is taken into account. To this end. In this paper, we develop a two degree of freedom ball bearing dynamics model that accounts for clearance variations. The complete response curve of system is obtained through harmonic balance and alternating frequency/time domain technique (HB-AFT), complemented by arc-length continuation (ALC) for solution tracking. The stability of the periodic solutions is discussed based on Floquet theory. Numerical analysis illustrates that the varying compliance (VC) vibration presents complex nonlinear characteristics when the varying clearance is introduced through the interference assembly model. Comparison with the results of the Rung-Kutta method indicates that one VC period (<i>T</i>-1) solution alone cannot comprehensively describe the system's dynamic response. The system exhibits rich and complex nonlinear behaviors, including <i>T</i>-2 solution, instances where both <i>T</i>-1 and <i>T</i>-2 solution coexist and chaos.</p>