The objective of this study is to evaluate the effect of constitutive equations on the prediction accuracy for springback in cold stamping with various deformation modes. In this study, two types of yield functions—Hill’48 and Yld2000-2d—were considered to describe yield behavior. Isotropic and kinematic hardening models based on the Yoshida–Uemori model were also adopted to describe hardening behavior. Various material tests (such as uniaxial tension, tension- compression, loading-unloading, and hydraulic bulging tests) were carried out to determine the material parameters of the models. The obtained parameters were implemented in the finite element (FE) simulation to predict springback, and the results were compared with experimental data. U-bending and T-shape drawing were employed to evaluate the springback prediction accuracy. Obviously, the springback prediction accuracy was greatly influenced by constitutive equations. Therefore, it is important to choose appropriate constitutive equations for accurate description of material behaviors in FE simulation.