در این پایان نامه سعی بر آن بوده است، که بهینه یابی پارامترهای توزیع منحنی های شکنندگی که به صورت ابتدایی وابسته به تک متغیر تصادفی تحریک هستند علاوه بر روش کلاسیک، با استفاده از الگوریتم بهینه یابی جستجوی هماهنگی محاسبه گردد. مزیت اضافه نمودن روش بهینه یابی جستجوی هماهنگی در بهینه یابی تابع بیشینه درست نمایی امکان استفاده از یک الگوریتم فرا ابتکاری را برای دخیل نمودن متغیرهای تصادفی دیگری فراهم می آورد. برای بدست آوردن توابع شکنندگی به داده های پاسخ سازه نیاز است. با استفاده از روش تحلیل دینامیکی فزاینده، می توان به مجموعه ای از داده های پاسخ سازه تحت شدت مشخص زلزله دست یافت. این داده ها بیانگر همه سطوح شدت نخواهند بود. بنابراین با برازش توزیع های احتمالاتی موجود مانند توزیع لوگ نرمال امکان پیش بینی پاسخ در همه سطوح فراهم می گردد.