In this work the robustness of residual stresses in finite element simulations with respect to deviations in mechanical parameters in castings is evaluated. Young's modulus, the thermal expansion coefficient and the hardening are the studied parameters. A 2D finite element model of a stress lattice is used. The robustness is evaluated by comparing purely finite element based Monte Carlo simulations and Monte Carlo simulations based on linear and quadratic response surfaces. Young's modulus, the thermal expansion coefficient and the hardening are assumed to be normal distributed with a standard deviation that is 10% of their nominal value at different temperatures. In this work an improved process window is also suggested to show the robustness graphically. By using this window it is concluded that least robustness is obtained for high hardening values in combination to deviations in Young's modulus and the thermal expansion coefficient. It is also concluded that quadratic response surface based Monte Carlo simulations substitute finite element based Monte Carlo simulations satisfactory. Furthermore, the standard deviation of the responses are evaluated analytically by using the Gauss formula, and are compared to results from Monte Carlo simulations. The analytical solutions are accurate as long as the Gauss formula is not utilized close to a stationary point.