The uncertainties of design variables, noise parameters, and metamodel are important factors in simulation-based robust design optimization. Most conventional metamodel construction methods only consider one or two uncertainties. In this paper, a new surrogate modeling method simultaneously measuring all the uncertainties is proposed for simulation-based robust design optimization of complex product. The effect of metamodel uncertainty on product performance uncertainty is quantified through uncertainty propagation analysis among design variables uncertainty, noise parameters uncertainty, metamodel uncertainty, and performance uncertainty. Then, the sampling points are selected and the metamodel is constructed based on the predictive interval of product performance and mean square error of the Kriging metamodel. The constructed metamodel is applied to robust design optimization considering multiple uncertainties. Results of two mathematical examples show that the proposed metamodel considering multiple uncertainties increases the result accuracy of robust design optimization. Finally, the proposed algorithm is applied to robust design optimization of a heat exchanger, and the total heat transfer rate is enhanced under uncertainties of fin structural parameters, operation conditions parameters and simulation metamodel.