This paper develops theoretical foundations for extending Gauss-Hermite quadrature to robust design with computer experiments. When the proposed method is applied with m noise variables, the method requires 4m + 1 function evaluations. For situations in which the polynomial response is separable, this paper proves that the method gives exact transmitted variance if the response is a fourth-order separable polynomial response. It is also proven that the relative error mean and variance of the method decrease with the dimensionality m if the response is separable. To further assess the proposed method, a probability model based on the effect hierarchy principle is used to generate sets of polynomial response functions. For typical populations of problems, it is shown that the proposed method has less than 5% error in 90% of cases. Simulations of five engineering systems were developed and, given parametric alternatives within each case study, a total of 12 case studies were conducted. A comparison is made between the cumulative density function for the hierarchical probability models and a corresponding distribution function for case studies. The data from the case-based evaluations are generally consistent with the results from the model-based evaluation.