A new method for the determination of the critical slip surfaces of slopes is proposed in this paper. In this paper, the original critical slip field method is extended in terms of the total residual moment, values of residual work as well as the unbalanced thrust force at the exit point for a given non-circular slip surface. The most critical slip surface with the maximum representative value for a prescribed factor of safety will be optimized and located using the harmony search algorithm. The prescribed factor of safety is modified with certain tiny interval in order to find the critical slip surface where the maximum representative value is zero. The aforementioned approach to the location of the critical slip surface is greatly different from the traditional limit equilibrium procedure. Three typical soil slopes are evaluated by use of the proposed method, and the comparisons with the classical approaches have illustrated the applicability of the proposed method.
