Design optimization of layered plates bonding process is conducted by considering uncertainties in a manufacturing process, in order to reduce the crack failure arising due to the residual stress at the surface of the adherent which is caused by different thermal expansion coefficients. Robust optimization is peformed to minimize the mean as well as its variance of the residual stress, while constraining the distortion as well as the instantaneous maximum stress under the allowable reliability limits. In this optimization, the dimension reduction (DR) method is employed to quantify the reliability such as mean and variance of the layered plate bonding. It is expected that the DR method benefits the optimization from the perspectives of efficiency, accuracy, and simplicity. The obtained robust optimal solution is verified by the Monte Carlo simulation.