In this paper, the initial blank shape is optimized via sequential approximate optimization (SAO) with radial basis function (RBF) network. In deep drawing, the wrinkling and tearing are major defects. In addition, too small initial blank shape will results in the boundary of blank drawn into a trimming contour, which is the geometrical defect. Therefore, these three defects are taken as the design constraints. On the other hand, the flange area, which is often called the earing, is taken as the objective function. Also, variable blank holder force (VBHF) trajectory is determined. In this paper, new algorithm to evaluate the earing and the geometrical defect is also proposed for the SAO. Numerical result shows that the optimization of the initial blank shape with VBHF trajectory leads to defect-free product. In comparison with the use of square blank shape with a constant blank holder force, 60% reduction of the waste of material can be achieved.