محاسبه ضرایب تمرکز تنش در نوک ترک همراه با پدیده کینک با استفاده از روش جابجایی- ناپیوستگی مرتبه سوم

چکیده
 روش المان مرزی غیرمستقیم به علت دقت محاسبات مناسب در تحلیل ناپیوستگی های در مقیاس ریزتر مورد توجه قرار گرفته است. روش جابجایی ناپیوستگی نسبت به دیگر روش عددی تنش موهومی انعطاف پذیری بیشتری برای حل مسائل گسترش ترک دارد. برای میل به محاسبات دقیق تر المان های با مراتب بالاتر به منظور تحلیل گسترش ترک به کار گرفته شده است. پیش از این روش جابجایی ناپیوستگی مرتبه دوم برای تحلیل ترک استفاده شده است. در این تحقیق معادلات براساس مرتبه سوم به تفکیک برای المان های ساده ترک، المان کینک و المان خاص نوک ترک ارائه شده اند. المان های ذکر شده توسط نقاط وسط المان به 4 زیر المان تقسیم شده که با تعریف توابع حالت جابجایی ناپیوستگی در این نقاط محاسبه می شود و با تعریف تابع پتانسیل می توان تنش و جابجایی ها در این نقاط را بدست آورد. همچنین با به کارگیری معیار تنش مماسی بیشینه و معیار بیشترین نرخ رهایی انرژی امکان محاسبه زاویه اولیه انتشار ترک، نرخ انرژی کرنشی و ضرایب تمرکز تنش نوک تر مقدور می باشد. این پارامترها به عنوان مهمترین معیار برای تحلیل رفتار ترک در نظر گرفته می شوند. برنامه عددی نوشته شده تحت روش جابجایی ناپیوستگی می تواند با تحلیل المان به المان ترک ( TDDCKCR) مرتبه سوم برای ترک همراه با پدیده کینک احتمال رخداد پدیده کینک را بر روی سگمنت های درنظر گرفته شده در مرز ترک بررسی نماید. با درنظر گرفتن یک المان خاص نوک ترک در ابتدا و انتهای ترک می توان ضریب تمرکز تنش در نوک ترک را محاسبه نمود. برای اثبات این موضوع یک ترک خمیده تحت زاویه 45 درجه برای بررسی مورد تحلیل قرار گرفته است و ضرایب تمرکز تنش نرمال شده TDDCKCR رویداد کینک با برنامه در نوک ترک برای نسبت طول المان نوک ترک به نصف طول المان ترک برابر 0.011 با اعمال تعداد المان های مختلف حاصل شده است. باردیگر روند محاسباتی با تغییر نسبت طول المان خاص نوک ترک به طول المان و تقسیم طول ترک به 60 المان تکرار شده است. در انتها نتایج با داده های نوشته شده با روش جابجایی ناپیوستگی مرتبه دوم مقایسه شده ( TDDKCR) حاصل از برناماست که نتایج حاکی از تطابق مناسب تر با داده های تحلیلی برای روش جابجایی ناپیوستگی مرتبه سوم ( کوبیک) می باشند.

واژه های کلیدی:

روش عددی، روش جابجایی- ناپیوستگی مرتبه سوم، ترک کینک، ضریب تمرکز تن

نویسندگان :

سهند شکری  ، محمد فاتحی مرجی

پنجمین کنفرانس مکانیک سنگ ایران (1393)

ارزیابی ریسکهای ژئوتکنیکی در حفاری تونل به روش شاخص آسیبپذیری و مقایسه آن با روش ماتریس ریسک حاصل از RMR

چکیده بلند فارسی
 ریسکهای ژئوتکنیکی در تونل زنی به شرایط مخاطره آمیزی گفته میشود که میتوانند تأثیر نامطلوبی بر روی ساخت و اجرای تونل گذاشته و یا حتی سبب به مخاطره انداختن جان افراد شوند. و همچنین با صدمه رسانیدن به دستگاهها، سبب ایجاد وقفه در اجرای طرح و افزایش هزینههای پروژه گردند. این مخاطرات و ریسکهای ناشی از آنها به طور مستقیم بر روی نرخ پیشروی تأثیر گذاشته و با کاهش آن طول زمان اجرای پروژه را افزایش میدهند. تحلیل ریسک شامل فرآیندهای تشخیص، اندازهگیری و کاهش ریسک ناشی از مخاطرات مورد نظر میباشد. بررسی مخاطرات زمینشناسی و ژئوتکنیکی و ریسکهای ناشی از این مخاطرات در بحثهای مربوط به طراحی و اجرای سازههای سنگی از جمله تونلها، یکی از مهمترین مواردی است که بایستی مورد توجه و مطالعه قرار گیرد. در این مطالعه هدف بررسی ریسکهای زمینشناسی و ژئوتکنیکی موجود در حفاری تونل به روش مکانیزه و کم  یسازی این ریسکها میباشد. در این مورد مطالعاتی هدف بررسی ارتباط و میزان انطباق روش محاسبه شاخص آسیبپذیری و روش در ارزیابی ریسکهای RMR ماتریس ریسک ناپایداری حاصل از سیستم طبقهبندی تودهسنگ زمینشناسی و ژئوتکنیکی موجود در حفاری تونل به روش مکانیزه است. بر این اساس در این مطالعه ریسک حفاری در ناحیههای 21 تا 30 تونل نوسود به روش محاسبه شاخص آسیبپذیری مورد بررسی قرار گرفت و نتایج آن با روش ماتریس ریسک ناپایداری حاصل از سیستم طبقهبندی مقایسه شد. نتایج بدست آمده از این مطالعه نشان داد که بین تحلیل ریسک RMR تودهسنگ حفاری به روش شاخص آسیبپذیری و روش ماتریس ریسک ناپایداری حاصل از سیستم طبقه- تشابه بالایی وجود دارد، به طوری که در ناحیههای مورد مطالعه روند RMR بندی تودهسنگ افزایش و کاهش ریسک مشابه است. و این نشان دهنده عملکرد مطلوب این سیستم طبقهبندی درارزیابی ریسک کلی ناپایداری ناشی از حفاری مکانیزه میباشد. بر اساس نتایج دو روش محاسبه ، RMR شاخص آسیبپذیری و ماتریس ریسک ناپایداری حاصل از سیستم طبقهبندی تودهسنگ 23 ، و 26 دارای کمترین مقدار ریسک میباشند، و ناحیههای 22 و 29 دارای بیشترین ، ناحیههای 21مقدار ریسک هستند.

واژه های کلیدی:

مخاطرات زمین شناسی، تحلیل ریسک، روش شاخص آسیبپذیری، روش ماتریس ، ریسک ،RMR

نویسندگان :

مجتبی ربیعی وزیری  ، وحید حسیبی  ، امیر دریایی

پنجمین کنفرانس مکانیک سنگ ایران (1393)

تخمین جهت و مقدار تنش های برجای نیروگاه تلمبه ذخیره ای آزاد با استفاده از روش انتشار آوایی

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

واژه های کلیدی:

مغار نیروگاه، تخمین تنش برجا، روش انتشار آوایی، اثر کایزر

نویسندگان :

محمود خسروی کارشک  ، مرتضی احمدی

پنجمین کنفرانس مکانیک سنگ ایران (1393)

بررسی الگوهای درزه در نواحی گسلی

چکیده 
 نواحی گسلی و بویژه نواحی تخریب گسل، تغییرات گسترده ای از هندسه ها و الگوهای درزه را نشان می دهند. در این مقاله، طبقه بندی جامعی از ناحیه گسلی بر اساس موقعیت شان در اطراف گسل، ساختارهای داخلی و الگوهای درزه ها در داخل هر ناحیه، ارائه شده است. ناحیه گسلی به ناحیه تخریب و ناحیه هسته گسل تقسیم می شود. ناحیه هسته به لای گسلی، بِرِش گسلی و کاتاکلازیت تقسیم بندی می شود. ناحیه تخریب نیز بر اساس موقعیت در اطراف ناحیه گسلی، به نواحی تخریب راسی، اتصالی و دیوار های تقسیم شده است. متعاقباً هر یک از این نواحی نیز بر اساس الگوی درزه ها به بخ شهای کوچکتری تقسیم می شود. در مورد گسل کامل، پهنای ناحیه تخریب از چند ده سانتی متر تا چند کیلومتر تغییر می کند، درحالیکه ناحیه هسته گسلی در ناحی های 0 تا 5 متر متمرکز می شود. رابطه مستقیمی بین مقدار لغزش گسل و پهنای / باریک با پهنای تقریباً 5 نواحی هسته و تخریب وجود دارد. پهنای ناحیه تخریب در دو طرف گسل م یتوانند کاملاً نامتقارن باشند. مطالعات جامعی در خصوص تغییرات چگالی و هندسه درز هها در این مقاله صورت گرفته است. نتایج بررسی ها نشان می دهد که اولاً به عنوان یک قاعده کلی، چگالی درزه ها، مستقل از جابجایی گسل، در مجاورت گسل بیشینه است و ثانیاً چگالی درزه ها با فاصله از هسته گسل، طی یک روند نمایی، خطی یا نامنظم کاهش می یابد. همچنین درزه هایی که در طول گس لزایی تولید می شوند، زوایای بین 15 تا 42 درجه با صفحه گسل اصلی تشکیل می دهند، درحالیکه درزه هایی که با زاویه بیش از 50 درجه ( متوسط 75 درجه) با گسل قرار می گیرند، درز ههایی هستند که اساساً در واکنش به تنش تکتونیکی منطقه ای تولید می شوند و مستقیما توسط فرآیند گسل زایی تحت تاثیر قرار نمی گیرد.

واژه های کلیدی:

الگوی درز هها، گسل، ناحیه تخریب، چگالی درزه

نویسندگان :

مهدی نوروزی ، رضا کاکایی ، سیدمحمداسماعیل جلالی

پنجمین کنفرانس مکانیک سنگ ایران (1393)

شکست مصالح دان های در شرایط محصور

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

نویسندگان :

وحید محمدنژاد  ، مهرداد امامی تبریزی ، حسن افشین

پنجمین کنفرانس مکانیک سنگ ایران (1393)

توسعه مدل تصادفی شبکه درزه ها با در نظر گرفتن ویژگی آماری اندازه درزه

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

واژه های کلیدی:

مدل تصادفی، اندازه درزه، تابع چگالی احتمال، برنامه کامپیوتری.

نویسندگان :

مهدی نوروزی ، سیدمحمداسماعیل جلالی ، رضا کاکایی

پنجمین کنفرانس مکانیک سنگ ایران (1393)

ارزیابی بهبود های صورت گرفته در کدگذاری ماتریس اندرکنش در رویکرد RES

چکیده
 از جمله رو شهایی که امروزه در حل مسائل پیچیده مورد استفاده قرار می گیرد، روش های مبتنی بر مبنای تفکر سیستمی بنا ( RES) بر تفکر سیستمی می باشند. رویکرد سیست مهای مهندسی سنگ RES نهاده شده و امروزه کاربرد گسترده ای در حل مسائل ژئوتکنیکی پیدا کرده است. در رویکرد با تشکیل ماتریسی تحت عنوان ماتریس اندرکنش به بررسی و تحلیل مسئله مورد نظر پرداخته می شود. پارامترهای اصلی مرتبط با مسئله مورد نظر در امتداد قطر اصلی ماتریس لیست شده و اندرکنش های هر جفت از پارامترها در سایر درایه ها شکل می گیرند. یکی از مهم ترین بخ شها در استفاده از این روش، نحوه کدگذاری ماتریس اندرکنش است. در این مقاله به بررسی مزایا و معایب انواع روش های کدگذاری ماتریس اندرکنش که تا به امروز ارائه شده است، پرداخته م یشود. این شبکه های عصبی، ،( FCM) استفاده از نظریه گراف، طرح شناخت فازی ، CQC ، ESQ روش ها شامل با وجود ضع فهای آشکار، تا حد ESQ هوش محاسباتی و استفاده از تئوری احتمالات است. روش زیادی موفق بوده و تاکنون از بیشترین کاربرد برخوردار بوده است. از مزایای روش کدگذاری می توان به محدوده وسیع تر عددی آن اشاره کرد که از 0 تا 100 را در نظر ESQ نسبت به CQC می گیرد. هدف اصلی از بهر هگیری از نظریه گراف در کدگذاری ماتریس اندرکنش، در نظر گرفتن به منظور FCM رابطه همزمان بیش از دو پارامتر در ماتریس اندرکنش بوده است. در روش + کدگذاری اندرکنش بین پارامترها و توصیف پویایی و پیچیدگی آن ها، از اعدادی در بازه 1 - تا 1 استفاده شده است. نتایج این بررسی نشان م یدهد که استفاده از شبکه های عصبی و یا احتمالات، کمک شایانی در جهت کاهش خطاهای ناشی از قضاوت مهندسی کارشناسان می نماید.

واژه های کلیدی:

تفکر سیستمی، رویکرد سیست مهای مهندسی سنگ (RES) ،ماتریس اندرکنش

نویسندگان :

عماد خراسانی ، مسعود زارع ، سید محمد اسماعیل جلالی ، شکرالله زارع

پنجمین کنفرانس مکانیک سنگ ایران (1393)

بهسازي لرزه اي قابهاي خمشی فولادي با مهاربندهاي مختلف

خلاصه
نظر به آنکه بخش قابل توجهی از سازه هاي موجود که بر اساس استاندار 2800 ویرایش اول طرح شده اند قادر به پاسخگویی عملکرد مورد انتظار تحت اثر زلزله هاي شدید نیستند، لذا در این تحقیق اقدام به بررسی یکی از روشهاي بهسازي لرزه اي سازه فولادي قاب خمشی متوسط (در دو جهت) با مهاربندها شده است.در این تحقیق 12وفق آیین نامه 2800 ویرایش اول طرح شده و براساس آیین نامه 2800 ویرایش سوم ، سه مدل سازه فولادي قاب خمشی متوسط (در دو جهت) در تعداد طبقات 8، 3 مورد بازنگري قرار گرفته است و ملاحظه گردید که این سازه ها عملکرد لرزه اي لازم را ندارند. سپس بهسازي لرزه اي آنها به دو روش : 1- با استفاده از مهاربند هاي صورت گرفته و با انجام تحلیل غیر خطی ضریب رفتار سازه هاي بهسازي شده مورد ارزیابی قرار (BRB) -2 با استفاده از مهاربند هاي کمانش تاب .(TBB) رایج گرفت.

کلمات کلیدي:

قاب خمشی متوسط،مهاربند رایج (TBB)، مهاربندکمانش تاب (BRB)،آنالیز استاتیکی-غیر خطی،ضریب رفتار

نویسندگان :

غلامرضا عبداله زاده ، هادي فقیه ملکی

دومین کنفرانس ملی سازه - زلزله - ژئوتکنیک (1391)

بررسی عملکرد قابهاي مهاربندي برون محور فولادي تحت زلزله هاي نزدیک گسل

خلاصه
قاب مهاربندي شده واگرا، نمایانگر یک سیستم قاب بندي شده ي مرکب است که هم سخت و هم شکل پذیر می باشد. در این قابها رفتار جانبی لرزه اي سازه ترکیبی از عملکرد خمشی _ برشی تیرها و ستون هاي دهانه مهاربندي شده و عملکرد کششی  فشاري مهاربندها می باشد. در این سیستم ایجاد شده است. به فاصله ایجاد شده در e در انتهاي دو عضو مهاري بین عضو مهاري با ستون و یا بین عضو مهاري با تیر طبقه، یک فاصله به اندازه این سیستم، پیوند اطلاق می گردد. طبق فلسفه طراحی این نوع قاب، پیوند ضعیفترین بخش قاب می باشد تا در برابر نیروهاي جانبی ناشی از زلزله شدید وارد ناحیه غیر خطی شده تا مقدمه استهلاك انرژي ورودي زلزله را فراهم آورد. در این مقاله عملکرد قابهاي مهاربندي برون محور فولادي در ساختمان هاي 6 و 12 طبقه تحت زلزله هاي نزدیک گسل بررسی شده اند. تحلیل استاتیکی غیرخطی و دینامیکی غیرخطی روي ساختمان هاي مدل شده ، انجام شده است. نتایج بیانگر آنست که اگر تیر پیوند طوري طراحی شود که رفتار آن برشی باشد، عملکرد آنها با استفاده از نتایج تحلیل استاتیکی غیر خطی، مطلوب بوده و در حد ایمنی جانی جوابگوي نیاز لرزه اي ما می باشند. طبق نتایج به دست آمده از تحلیل غیر خطی سازه ها، در برابر زلزله هاي نزدیک گسل و تعداد خرابی ها در سازه ها به این نتیجه می رسیم که با افزایش تعداد طبقات، سازه ها آسیب پذیر تر بوده و در طراحی آنها باید دقت کرد.

کلمات کلیدي:

عملکرد، نزدیک گسل، مهاربند واگرا، تحلیل غیر خطی، نیاز لرزه اي

نویسندگان :

یاسر مظفري جویباري ، یاسر محمودي ، شهرام مقصودیان

دومین کنفرانس ملی سازه - زلزله - ژئوتکنیک (1391)

نیروي استاتیکی معادل و مقایسه ي آیین نامه ها در طراحی لرزه اي سازه ها

خلاصه
هرچند که طی سالیان اخیر آیین نامه 2800 ایران تغییرات عمده اي نسبت به ویرایش هاي گذشته داشته است ولی مقایسه آن با دیگر آیین نامه هاي معتبر دنیا می تواند نقاط ضعف و قوت آن را مشخص نموده و در ویرایش هاي بعدي موثر باشد. در این تحقیق به مقایسه چهار آیین نامه اروپا) تحت یک شتاب مبنا و نوع خاك یکسان پرداخته شده است. سپس طیف ) Eurocode آمریکا) و 8 ) UBC- هند)، 97 ) IS ،( 2800 (ایران پاسخ هر آیین نامه، برش پایه، تغییر مکان نسبی طبقات و زمان تناوب تجربی و تحلیلی براي هر قاب به صورت جداگانه مورد بررسی قرار گرفته است و در نهایت با توجه به خروجی هاي ذکر شده نتایج مطلوبی از مقایسه ي این آیین نامه ها استخراج گردید.

کلمات کلیدي:

طیف طراحی، قاب خمشی بتنی متوسط، آیین نامه لرزه اي، برش پایه، تغییر مکان نسبی

نویسندگان :

محسن گرامی ، فاطمه عبدوس ، نوید سیاه پلو  ،محمدهادي تقدیسی

دومین کنفرانس ملی سازه - زلزله - ژئوتکنیک (1391)