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