Abstract：Based on the computational results from elasto-plasticity finite element analysis，an initial value problem related to a system of ordinary differential equations(ODEs) is formulated to define the potential slide lines(PSL) for two dimensional cases. A prediction-correction algorithm for the ODEs and a necessary and sufficient condition that assures the convergence of the algorithm are presented. The skills for searching PSL with presented method are discussed，and the procedures of this method have been validated against traditional limit analysis methods and limit equilibrium method.
