In the demonstration phase of a missile scheme, to obtain the optimum proposal, designers need to modify the parameters of the overall structure frequently and significantly, and perform the structural analysis repeatedly. In order to reduce the manual workload and improve the efficiency of research and development, a general integrated method of missile structure modeling, analysis and optimization was proposed. First, CST (Class and Shape transformation functions) parametric method was used to describe the general structure of the missile. The corresponding software geometric modeling and FEM (Finite Element Method) analyzing of the missile were developed in C/C++ language on the basis of the CST parametric method and UG (Unigraphics) secondary development technology. Subsequently, a novel surrogate model-based optimation strategy was proposed to obtain a relatively light mass missile structure under existing shape size. Eventually, different missile models were used to verify the validity of the method. After executing the structure modeling, analysis and optimization modules, satisfactory results can be obtained that demonstrated the stability and adaptability of the proposed method. The method presented saves plenty of time comparing to the traditional manual modeling and analysis method, which provides a valuable technique to improve the efficiency of research and development.3D modeling program based on the UG secondary development technology, which has improved the modeling efficiency. Wang et al. [4] used the technology of UG/KF secondary development for the automatic modeling of the wind turbine blades, and the aerodynamic characteristics of the impeller were analyzed by the Fluent software. Moreover, He et al. [5] implemented the parametric modeling of the blended-wing-body underwater glider structure by using UG secondary development technology and studied its structural performance with ANSYS. Therefore, UG is considered as an excellent tool that can realize the integration of modeling and simulation.Meanwhile, airfoils have been parameterized in a wide variety of ways. Among them, CST is considered to be a more effective method, which can represent the different types of geometries in a generic way with fewer parameters. Kulfan and Bussoletti [6] proposed a parameterization method that based on class and shape transformation functions, which has the advantages of high precision and fewer variables. The CST method is a powerful and versatile means to describe complex aeronautical shapes ranging from 2D airfoils to 3D geometries of aircraft using analytical functions. Since it has been introduced, this parameterization method has been the subject of several studies [7]. Marco Ceze of the University of Michigan studied the characteristics of CST parameterization and found that the ill-conditioning of parameterization when fitting airfoil with high-order Bernstein polynomials [8]. Straathof et al. proposed a parameterization method that used a combination of Bernstein polynomials and B-splines to allow...