Three-dimensional finite element analysis should be used in stability analysis of slope because it can overcome the short advantages of two-dimensional finite element and can simulate the complex topographic and geological conditions. Based on the large-scale triaxial shear test, the modified Duncan-Chang model is established. Based on strength reduction elasto-plastic finite element, stability of high fill embankment was studied with three-dimensional finite element method considering the complex terrain conditions. Study results suggest that plastic strain and displacement mutant of slip surface node can be a sign of slope instability as a whole. At the same time calculation of three-dimensional finite element also does not converge. Therefore, it can be slope instability criterion calculate whether the finite element static analysis converges or not. On the other hand, stability safety factor of high fill embankment under three-dimensional conditions is larger than that of two-dimensional conditions, which shows that boundary conditions of high fill embankment enhance its stability.