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