In this paper, we introduce a general framework to efficiently optimize phased array antennas integrated with frequency selective surfaces (FSS) and feeding elements. The framework consists of two main parts: the first being parametric models that enable fast evaluation of the full-wave problem for a range of frequency, scan angle and geometry parameters. The second part consists of a number of optimization methods to generate an optimum set of geometry parameters by minimization of a penalty function. Results are shown for parametric models and the optimization of various examples. (7 pages)