Abstract This paper presents an efficient damage model based on the strong coupling of both anisotropic elasto-plasticity and isotropic ductile damage. The elasto-plasticity is modeled with mixed non-linear work hardening considering both isotropic and kinematic effects. The developed model is then coupled with the finite element method using ABAQUS finite element code. First, the formulation of the proposed damage model is presented. A comparison between numerical and experimental results is presented, in the context of the square cup deep-drawing benchmark test of Numisheet 1993. The results demonstrate the capability of the model to predict where and when the damaged zones will appear in the workpiece during the forming operation.