Limit equilibrium methods of slope stability analysis have been widely adopted mainly due to their simplicity and applicability. However, the conventional methods may not give reliable and convincing results for various geological conditions such as nonhomogeneous and anisotropic soils. Also, they do not take into account soil slope history nor the initial state of stress, for example excavation or fill placement. In contrast to the limit equilibrium analysis, the analysis of deformation and stress distribution by finite element method can deal with the complex loading sequence and the growth of inelastic zone with time. This paper proposes a technique to determine the critical slip surface as well as to calculate the factor of safety for shallow failure on partially saturated soil slope. Based on the effective stress field in finite element analysis, all stresses are estimated at each Gaussian point of elements. The search strategy for a noncircular critical slip surface along weak points is appropriate for rainfall-induced shallow slope failure. The change of unit weight by seepage force has an effect on the horizontal and vertical displacements on the soil slope. The Drucker-Prager failure criterion was adopted for stress-strain relation to calculate coupling hydraulic and mechanical behavior of the partially saturated soil slope.