本发明公开了一种基于极限分析下限定理的三维边坡安全系数迭代求解方法，该方法以块体界面作用力的大小、方向和作用点为主要变量，以满足极限分析下限定理许可静力场要求的静力平衡和屈服准则为约束条件，将边坡稳定分析问题转化为求解安全系数最大值的非线性优化问题；同时引入过渡变量临界加速度系数η，降低由于强度折减系数F引起的高度非线性；采用二分法迭代求得当η＝0时对应的折减系数即为所求安全系数。本发明为实际复杂边坡工程的三维稳定性分析提供了高效的手段。