Mesh generation is the key to the numerical simulation in CFD with complex aerodynamic shapes. The algorithm of mesh generation is required to be fast and robust. First, this framework established efficient initial tetrahedron arithmetic based on viewed-able face which is fit to the arbitrary polyhedron. Then, Delaunay generation method was improved to refine meshes and the effect of different mesh quality criterion was discussed and the divergence with Delaunay face swapping was developed to optimize the mesh generation process, Finally, a vertex spring model for mesh optimization was established. The results indicate that the method established in the research is efficient and stable.