This paper proposes a local search optimizer that, employed as an additional operator in multiobjective evolutionary techniques, can help to find more precise estimates of the Pareto-optimal surface with a smaller cost of function evaluation. The new operator employs quadratic approximations of the objective functions and constraints, which are built using only the function samples already produced by the usual evolutionary algorithm function evaluations. The local search phase consists of solving the auxiliary multiobjective quadratic optimization problem defined from the quadratic approximations, scalarized via a goal attainment formulation using an LMI solver. As the determination of the new approximated solutions is performed without the need of any additional function evaluation, the proposed methodology is suitable for costly black-box optimization problems.