Up to now, many constitutive models for clay have been proposed and studied based on the elasto-plastic or elasto-viscoplastic theory and it has been recognized that the effect of time on the loading process is a salient feature. In the present study, cyclic behavior characteristics of clay was studied with a viscoelastic-viscoplastic constitutive model for clay based on the non-linear kinematic hardening rule. In order to examine the behavior of clay several cyclic untrained triaxial tests and also their numerical simulations were performed. As results of that, it was found that the proposed model can well describe cyclic behaviors of clay such as frequency dependent characteristics, and have the high feasibility of numerical simulation for dynamic analysis.