International audienceIn the robust shape optimization context, the evaluation cost of numerical models is reduced by the use of a response surface. Multi-objective methodologies for robust optimization that consist in simultaneously minimizing the function and a robustness criterion (the second moment) have already been developed. However, efficient estimation of the robustness criterion in the framework of time-consuming simulation has not been greatly explored. A robust optimization procedure based 15 on the prediction of the function and its derivatives by kriging is proposed. The second moment is replaced by an approximated version using Taylor expansion. A Pareto front is generated by a genetic algorithm named NSGA-II with a reasonabl...