A simplified numerical procedure to predict drawbead restraining forces(DBRF) has been developed based on the hybrid membrane/bending method which superposes bending effects onto membrane solutions. As a semi-analytical method, the new approach is especially useful to analyze the effects of various constitutive parameters. The present model can accommodate general anisotropic yield functions along with non-linear isotropic-kinematic hardening under the plane strain condition. For the preliminary results, several sensitivity analyses for the process and material effects such as friction, drawbead depth, hardening behavior including the Bauschinger effect and yield surface shapes on the DBRF are carried out.