Design of modern antenna systems has become highly dependent on computational tools, especially full-wave electromagnetic (EM) simulation models. EM analysis is capable of yielding accurate representation of antenna characteristics at the expense of considerable evaluation time. Consequently, execution of simulation-driven design procedures (optimization, statistical analysis, multi-criterial design) is severely hindered by the accumulated cost of multiple antenna evaluations. This problem is especially pronounced in the case of global search, frequently performed using nature-inspired algorithms, known for poor computational efficiency. At the same time, global optimization is often required, either due to multimodality of the design task or the lack of sufficiently good starting point. A workaround is to combine metaheuristics with surrogate modeling methods, yet a construction of reliable metamodels over broad ranges of antenna parameters is challenging. This work introduces a novel procedure for global optimization of antenna structures. Our methodology involves a simplex-based automated search performed at the level of approximated operating and performance figures of the structure at hand. The presented approach capitalizes on weakly-nonlinear dependence between the operating figures and antenna geometry parameters, as well as computationally cheap design updates, only requiring a single EM analysis per iteration. Formal convergence of the algorithm is guaranteed by implementing the automated decision-making procedure for reducing the simplex size upon detecting the lack of objective function improvement. The global optimization stage is succeeded by gradient-based parameter refinement. The proposed procedure has been validated using four microstrip antenna structures. Multiple independent runs and statistical analysis of the results have been carried out in order to corroborate global search capability. Satisfactory outcome obtained for all instances, and low average computational cost of only 120 EM antenna simulations, demonstrate superior efficacy of our algorithm, also in comparison with both local optimizers and nature-inspired procedures.