An approach for robustness analysis of non-dominated solutions to a multi-objective optimization model of an energy management system aggregator (EMSA) in face of uncertainty is presented. The EMSA is an intermediary entity between households and the System Operator (SO), capable of contributing to balance load and supply, and therefore coping with the intermittency of renewable energy sources (RES) and facilitating a load follows supply strategy in a Smart Grid environment. Household clusters provide load flexibility to satisfy system services requested by the SO, involving decreasing or increasing load in specific time slots. The EMSA multi-objective optimization model considers the maximization of profits and the minimization of the imbalance between the amounts of load flexibility provided by the end-user clusters to satisfy SO requests, taking into account revenues from the SO and payments to the clusters. A hybrid evolutionary approach combining Genetic Algorithms (GA) with Differential Evolution (DE) has been designed to deal with this model, and its behaviour subject to different scenarios of uncertainty is evaluated. The robustness analysis of non-dominated solutions produced by the hybrid evolutionary approach is based on the degree of robustness concept, taking into account the changes in the performance of the objective functions when small perturbations of the model nominal coefficients occur.