In the present study, in order to apply a cyclic viscoelastic-viscoplastic constitutive model to multi-layered ground conditions during large earthquake, the numerical simulations of the 1995 Hyogoken Nanbu Earthquake at Port Island, Kobe, Japan, were performed by the seismic response analysis. From the seismic response analysis, it was found that the acceleration calculated from the cyclic elasto-viscoplastic model and cyclic viscoelastic-viscoplastic models for clay was in close agreement with the recorded accelerations at the Port Island down-hole array, and the cyclic elastic-viscoplastic and viscoelastic-viscoplastic constitutive models showed little different behavior characteristics near clay layer. Thus, the propriety of viscoplastic model for clay was convinced. Therefore, it can be concluded that a cyclic viscoelastic-viscoplastic constitutive model can give a good description of the amplification and also it showed accurate damping characteristics of clay during large event which induces plastic deformation in large strain range.