The presence of wet pollution on the upper surface of a string insulator increases the electric field on the insulator surface, especially at the triple junction (pin-cement, cement-porcelain) as well as at the surrounding air of the insulator. The rise of the electric field leads to the ionization of the air surrounding the insulator. This phenomenon, called corona discharge, is accompanied by several consequences that are harmful to the electricity transmission network, such as electromagnetic interferences, energy losses, visible light, audible noise, and the destruction of materials by erosion. If favorable conditions are gathered, it may even cause the flashover of the insulator. Designing an optimal insulator shape that reduces this electric field value at the triple junction will be an important achievement in enhancing the performance of electrical grids. The objective of this paper is to evolve a hybrid algorithm based on the GWO-FEM for optimizing the shape and the electrical performance of a string insulator. To achieve this purpose, this work is structured in four parts. First, modeling of the insulator string geometry is conducted in Comsol-multiphysics, then FEM computation of the electric field on the polluted surface of the string insulator is completed, and the maximum electric field value at the triple junction is saved as the fitness function that will be sent to the GWO algorithm to be optimized (minimized). The third part of the work is devoted to the coding of the constrained (electrical and geometrical constraints) GWO algorithm in a Matlab interface, and finally, coupling the GWO code with the FEM code. This study is achieved in wet polluted conditions. The results are given in both 2D and 3D representations. From the obtained results we can confirm that the developed GWO-FEM hybrid algorithm for optimizing insulator strings is a very promising tool for designing and enhancing the shape and the electrical performance of insulators.