امروزه بکارگیری رویکردهای عددی به خصوص رویکرد عددی تفاضلات محدود و نرم افزار FLAC3D در ارزیابی-های مربوط به طراحی های ژئوتکنیکی یک اصل غیرقابل اغماض می باشد. توانمندی این رویکرد عددی در تحلیل شرایط پیچیده محاسباتی سبب گردیده تا بصورت گسترده در تحلیل های مهندسی پی بویژه شمع بکار برده شود. این نرم افزار یک برنامه رایانه ای سه بعدی براساس روش حل عددی برمبنای تفاضلات محدود برای مدل سازی، مش بندی، تحلیل و محاسبه تنش ها و تغییرمکان المان های سازه ای در اندرکنش با خاک است. در این پژوهش از نرم افزار FLAC3D برای ارزیابی رفتار و اندرکنش شمع-رادیه-خاک برای مدل تک شمعی منفرد، زوج شمعی و گروه شمع؛ بهره گرفته شده که عوامل مؤثر در میزان توزیع بارگذاری محوری در شمع ها و پی بصورت یک سیستم ترکیبی مورد توجه بوده است. در این راستا تغییرات هندسی، ضخامت، تعداد، نوع مقطع شمع (دایره یا مربع) و طول شمع ها به عنوان متغیرهای اصلی تحلیل مد نظر بوده است. توزیع کنتوری و پیوسته میدان تنش و متعاقب آن رخداد تغییرشکل های پلاستیک تحت مدل های رفتاری الاستیک (برای شمع) و الاستوپلاستیک (برای خاک) در این مطالعه مورد بررسی قرار گرفته است. برپایه نتایج حاصل از شبیه سازی مشخص گردیده که در مدل تک شمع، تمرکز اصلی تنش در نوک و محور اصلی شمع بوده که سبب ایجاد تغییرشکل پیوسته در خاک بیشتر خود سازه گردیده است. این مسئله برای مدل زوج شمع تعدیل داشته است و روند مشابهی دنبال می کند، اما نتایج نشان داده که در مدل زوج شمع میزان باربری با افزایش طول شمع افزوده شده و تا 85% نیز می رسد. لازم به بیان است که تمرکز تنش برای زوج شمع های 30 متری بجای نوک شمع در بدنه شمع تمرکز نشان داده که بیانگر تحمیل تغییرشکل های جانبی توسط خاک به شمع است. در مدل گروه شمع روند زوج شمع در افزایش باربری وجود داشته اما در مورد سیستم مرکزگرای گروه شمع نتیجه ای عکس بهمراه دارد. این مسئله تمرکز تنش ناشی از فشار جانبی خاک در بدنه شمع است. از سوی دیگر باربری در رادیه با شمع های با مقطع مربعی علیرغم اینکه ضخامت یکسانی با رادیه سیستم رادیه و شمع مقطع استوانه ای دارد، کمتر بوده و افزایش طول شمع ها این تفاوت را کاهش می دهد. از سوی دیگر نسبت ظرفیت باربری کل سیستم با شمع های مقطع مربعی با افزایش طول شمع ها نسبت به سیستم با شمع های استوانه ای کمتر می شود و
