In this paper a phenomenological material model for simulation of sheet metal forming processes was introduced. This model is able to describe the anisotropic behavior of sheet metals in both yield stresses and plastic strain ratios (r-values) by using the non-associated flow rule and quadratic yield and potential functions. Additionally, to reproduce an accurate prediction of cyclic plastic deformation phenomena, a two-surface mixed isotropic-nonlinear kinematic hardening model was combined with the quadratic non-associated anisotropic formulation. This mixed kinematic hardening model is amongst the most sophisticated models with acceptable degree of complexity and minimum requirement of experimental tests for material coefficients. The main advantage of the model over the more complex nonquadratic yield and potential functions along with associated or non-associated flow rules is its simplicity and computational efficiency. The plasticity foundation of the model was introduced and then a general return mapping algorithm for numerical stress integration of the constitutive model was developed in order to implement it into a finite element code. Finally, the model was used to simulate both forming stage and subsequent springback of a deep drawing problem. The results showed that the model can accurately predict springback as well as earing phenomenon of the stamped part.