Three kinds of viscoelastic damper model, which has a non-linear spring as an element is studied analytically and numerically The behavior of the damper model shows non-linear hysteresis curves which is qualitatively similar to those of real viscoelastic materials. The motion is governed by a non-linear constitutive equation and an additional equation of motion. Harmonic balance method is applied to get analytical solutions of the system. The frequency-response curves sallow that multiple solutions co-exist and that the jump phenomena can occur. In addition, it is shown that separate solution branch exists and that it can merge with the primary response curve. Saddle-node bifurcation sets explain the occurrences of such non-linear Phenomena.