Considering the influences of harmonic moment and roll eccentricity, a dynamic model of horizontal vibration for work roll was established, and the amplitude-frequency characteristic equations of horizontal principal resonance, super-harmonic resonance and sub-harmonic resonance of work roll under the influences of harmonic moment and roll eccentricity were obtained by multi-scale approximate analytic method. Then, the amplitude-frequency characteristics of the system under the changes of damping coefficient, nonlinear coefficient of rolling force and harmonic moment were analyzed. The results show that the amplitude and resonance region of the amplitude-frequency characteristic curve decrease with the increasing of the damping coefficient of the system, the curvature and amplitude of the amplitude-frequency characteristic curve increase with the increasing of the nonlinear coefficient of rolling force and the harmonic moment, and the stability of the system is reduced. Finally, the model was numerically solved, and the motion rules of the system under the influences of harmonic moment and roll system eccentricity were analyzed by the methods such as phase trajectory and bifurcation chaos diagrams. It is found that the chaotic motion range of work roll in the horizontal direction increases when the harmonic moment and the roll eccentricity force work together, making the system more prone to instability. By comparing the numerical solution of the model with the approximate solution, it is further verified that the harmonic moment and the roll eccentricity force reduce the stability of work roll. The results provide a theoretical basis for improving the operation stability of work roll.