In a finite element analysis of the metal forming processes having large plastic deformation, largely distorted elements are unstable and hence they influence upon the result toward negative way so that adaptive remeshing is required to avoid a failure in the numerical computation. Therefore automatic mesh generation and regeneration is very important to avoid a numerical failure in a finite element analysis. In case of generating quadrilateral mesh, the automation is more difficult than that of triangular mesh because of its geometric complexity. However its demand is very high due to the precision of analysis. Thus, in this study, an automatic quadrilateral mesh generation and regeneration method using grid-based approach is developed. The developed method contains decision of grid size to generate initial mesh inside a two dimensional domain, classification of boundary angles and inner boundary nodes to improve element qualities in case of concave domains, and boundary projection to construct the final mesh.