The optimization of vehicle suspension kinematic/compliance characteristics is of significant importance in the chassis development. Practical suspension system contains many uncertainties which may result from poorly known or variable parameters or from uncertain inputs. However, in most suspension optimization processes these uncertainties are not accounted for. This study explores the use of Chebyshev polynomials to model complex nonlinear suspension systems with interval uncertainties.In the suspension model, several kinematic and compliance characteristics are considered as objectives to be optimized. Suspension bushing characteristics are considered as design variables as well as uncertain parameters. A high-order response surface model using the zeros of Chebyshev polynomials as sampling points is established to approximate the suspension kinematic/compliance model. For the uncertain optimization problem, a double-loop optimization process is used, in which the inner loop is to calculate the bounds of the interval design functions and the outer loop is to achieve optimum suspension kinematic/compliance characteristics. Numerical results show the effectiveness of the proposed optimization method.