In this paper, we present an improved and efficient analytical solution for evaluating the stability of 3D slopes based on the minimum potential energy method. The equilibrium equation of the landslide is employed to determine the shear stress on the failure surface, and a new formulation for shear potential energy is introduced. The vector-combined direction of the force acting on the landslide is considered equivalent to the sliding direction (SD). Two classic examples are reanalyzed to validate the accuracy and performance of the proposed method. The analysis demonstrates that the obtained results align well with reference solutions and tend to be conservative in engineering applications. Parameter investigations are conducted to explore the influences of SD, shear strength parameters, slope angle, and soil unit weight on the safety factor (SF) of 3D slopes. The findings indicate that a larger SD corresponds to a less stable slope. Finally, the effects of optimization parameters, landslide volume, and surcharge on the location of the failure surface, SF, and SD are also examined. The developed analytical approach offers a novel and straightforward solution with sufficient accuracy and eliminates the need for iteration in determining the SF, making it a valuable complementary method for assessing the stability of 3D slopes.