For optimization of real-world arch dams, it is unavoidable to consider two or more conflicting objectives. This paper employs two multi-objective differential evolution algorithms (MoDE) in combination of a parallel working MATLAB-APDL code to obtain a set of Pareto solutions for optimal shape of arch dams. Full dam-reservoir interaction subjected to seismic loading is considered. A benchmark arch dam is then examined as the numerical example. The numerical results are compared to show the performance of the MoDE methods.