Abstract: To deal with the problems that boundaries of control variables are difficult to determine in traditional slope stability analysis with circular sliding surface method, a new method is put forward in which arc height and X values of upper and lower points of intersection are regarded as control variables of arc. The relationships between each control variable and safety factor are analyzed. Based on these analyses, a global optimization method is presented which is validated through some engineering examples.