Stabilizing pile is a widely used method to reduce the development of large-scale landslides. Optimizing the pile geometry is a great challenge in the design of stabilizing piles with the purpose of cost-effectiveness, especially for soil strength parameters with large uncertainty. The objective of this study is to propose a robust and efficient method of designing piles for landslide stabilization with the consideration of the safety of slope, uncertainty of soil parameters, and cost of stabilizing piles. A new response surface, which incorporates soil parameters and stabilizing force into a quadratic polynomial function, is first proposed. Unknown coefficients of the quadratic polynomial function are solved with a numerical method at typical sampling points. Based on the solved quadratic polynomial function, the mean and standard deviation of factor of safety (FOS) of the pile-stabilized slope as well as the signal-to-noise factor are then calculated in order to evaluate the design robustness. A framework based on the concept of robust geotechnical design is presented, and its feasibility is illustrated by two cases of soil slopes. The results indicate that the proposed robust geotechnical design method could be used to optimize the design of landslide-stabilizing piles.