Three constitutive models were investigated based on the concept of a non-associated flow rule. The key concept is the assumption that sand can be regarded as an isotropic work hardening elasto-plastic material. It was proved that the proposed elasto-plastic models, especially Model I could effectively represent the yield and stress-deformation characteristics of sand.