Deep-drawing operations are performed widely in industrial applications. It is very important for efficiency to achieve parts with no defects. in this work, a finite element method is developed to simulate deep-drawing operation including wrinkling. A four nodded five degree of freedom shell element is formulated. Isotropic elasto-plastic material model with Von Mises yield criterion is used. By using this shell element, the developed code can predict the bending behavior of workpiece besides membrane behavior. Simulations are carried out with four different element sizes. The thickness strain and nodal displacement values obtained are compared with results of a commercial finite element program and results of previously conducted experiments.