In this paper, a new approach is proposed and addressed for designing vehicle suspension systems based on the scheme of multi-objective programming. For complicated vehicle random vibration, a linear model is used to describe the dynamic behavior of vehicles running on randomly profiled roads. The road irregularity is regarded as a Gaussian random process. Pesudo excitation method has been used to solve the dynamic responses. And a Kriging model is introduced to build the approximate mapping relationship between the design variables and the responses. Optimal solutions are derived by means of the method of centers for structural optimization with multiple objectives. Numerical examples are given, and compared with other optimization methods.