Most structures in the real life are subject to multiple load cases. This paper aims at extending the evolutionary structural optimization (ESO) algorithm to optimal contact shape design for elastic bodies under the multiple load cases. To evaluate the reference stresses of each contact node in a finite element framework, an extreme stress criterion (the worst case design) and a weighted average criterion (Pareto design) are presented. In the extreme stress method, the highest nodal contact stress under all load cases is adopted as the reference level. In the weighted average method, the weighted sum of nodal contact stresses over all the load cases is regarded as the reference. It is found that these two criteria can produce different resu...