Uncertainties in engineering design may lead to low reliable solutions that also exhibit high sensitivity to uncontrollable variations. In addition, there often exist several conflicting objectives and constraints in various design environments. In order to obtain solutions that are not only multi-objectively optimal, but also reliable and robust, a probabilistic optimization method was presented by integrating six sigma philosophy and multi-objective genetic algorithm. With this method, multi-objective genetic algorithm was adopted to obtain the global Pareto solutions, and six sigma method was used to improve the reliability and robustness of those optimal solutions. Two engineering design problems were provided as examples to illustrate the proposed method.