In order to maintain the diversity of non-dominated solutions in multi-objective optimization algorithms efficiently the authors have proposed a multi-objective artificial physics optimization algorithm based on virtual force sorting (VFMOAPO). Adopting quick-sort idea, the individuals in non-dominated solutions set were sorted by the total virtual force exerting on the other individuals. So the non-dominated solution set was pruned and the individual with the maximal sum of virtual force exerting on the other individuals was selected as the global best solution. Some benchmark functions were tested for comparing the performance of VFMOAPO with MOPSO, NSGA and SPEA. The simulation results show the algorithm is feasible and competitive.