Abstract: According to Bishop’s method and non linear programming theory, this paper renders a rather complete numerical model of optimization for method of circular sliding surface applied widely in engineering, and uses obstacle function method to solve it. For avoiding the risk to get in "part optimization solution", a global optimization method is established through applying a system of randomly obtained points in the paper. The analysis of some engineering examples shows satisfactory results.