As a new type of suspension bearing, Magnetic-Liquid Double Suspension Bearing (MLDSB) is mainly supported by electromagnetic suspension and supplemented by hydrostatic supporting. Its bearing capacity and stiffness can be greatly improved. Because of the small liquid film thickness (it is smaller 10 times than air gap), the eccentricity, crack, bending of the rotor, and the assembly error, it is easy to cause a clearance-rubbing fault between the rotor and stator. The coating can be worn and peeled, the operating stability can be reduced, and then it is one of the key problems of restricting the development and application of MLDSB. Therefore, the clearance-rubbing dynamic equation of 2-DOF system of MLDSB is established and converted into Taylor Series form and the nonlinear components are retained. Dimensionless treatment is carried out by dimensional normalization method. Finally, the rotor displacement response under different rotor eccentricity ratio and rotating speeds is numerically simulated. The studies show that the trajectory of the rotor is periodic elliptic without clearance-rubbing phenomenon when the eccentricity ratio is less than 0.2, while the rotor is greatly affected by the rotation speed and a variety of motions, such as single-period, quasi-period, double-period and chaos, are presented when greater than 0.3. Within the largest range of rotating speed and eccentricity ratio, the rotor presents the single-period trajectory, and then the number of Poincare mapping point is 1, without a clearance-rubbing fault. When the rotational speed is in the scope of (9, 13) krpm and the eccentricity ratio is in the scope of (0.27, 0.4), the number of Poincare mapping point is more than one, the maximum dimensionless rubbing force is −5.7, and then clearance-rubbing fault occurs. The research can provide a theoretical basis for the safe and stable operation of MLDSB.