حوادثی که سالانه به دلیل بروز پدیده¬ی گسیختگی شیب های خاکی در سراسر جهان به وقوع می پیوندد، موجب شده تا مدل های تعیین پایداری این شیب ها مورد توجه محققین قرار گیرد. در مدل های موجود با قرار دادن مقادیر متوسط پارامترهای اثر گذار، تنها یک مقدار به عنوان ضریب اطمینان پایداری شیب تعیین می گردد. این در حالی است که در روش های احتمالاتی با در نظر گرفتن دامنه ای از مقادیر محتمل برای هر یک از پارامترهای ورودی، به تمامی مقادیر محتمل بر آن پارامتر شانس وقوع داده می شود. در این پژوهش تلاش شده با استفاده از روش احتمالاتی تحلیلی ترکیب منحنی تابع توزیع متغیر های تصادفی و بر اساس روش های تعادل حدی قطعه ای ساده شده ی بیشاپ، ساده شده ی جانبو، اسپنسر و مورگشترن-پرایس، قابلیت اطمینان پایداری شیب های خاکی محدود مورد ارزیابی قرار گیرد. برای این منظور پارامتر های مقاومتی خاک از قبیل چسبندگی (c)، زاویه اصطکاک داخلی (?) و وزن مخصوص (?)، بصورت غیر قطعی و با استفاده از توزیع های احتمالاتی نرمال و لوگ نرمال و با فرض ضریب همبستگی توصیه شده در نظر گرفته و تابع چگالی احتمال ضریب اطمینان پایداری شیب روش های مذکور با یکدیگر قیاس گردیده است. همچنین قابلیت اطمینان پایداری لرزه ای شیب های خاکی با استفاده از روش شبه-استاتیکی بیشاپ نیز مورد بررسی قرار گرفت. در هر مرحله سطح گسیختگی بحرانی شیب های مورد مطالعه نیز با استفاده از الگوریتم بهینه یابی تجمع ذرات تعیین شده اند. از جمله دستاورد های این تحقیق می توان به تعیین مدل های احتمالاتی برای حالات تنش کل و لرزه ای اشاره کرد. با استفاده از این مدل ها می توان، احتمال گسیختگی شیب های خاکی را مستقیماً بر مبنای ضریب اطمینان قطعی و بدون نیاز به طی مراحل زمانبر و پیچیده ی متداول تعیین کرد.