فهرست مطالب
مجله مهندسی عمران شریف
سال بیست و نهم شماره 2 (تابستان 1392)
- تاریخ انتشار: 1392/08/20
- تعداد عناوین: 13
-
-
صفحات 3-12برای سری های فوریه و چبیشف، تحلیل خطا در تقریب ضرایب انتگرالی کوادراتیک گاوسین از راه ساده ی تحلیل همانندسازی ممکن است. برای ایجاد کدهای شبه طیفی مناسب و پایدار با حذف این خطا در مسائل مکانیک سیالات استفاده از فرمولاسیون «قانون 3/2» پیشنهاد می شود. در این نوشتار به بررسی تاثیر الگوهای همانندسازی غیرمحذوف، حذف همانندسازی در راستای قائم، حذف همانندسازی در راستای طولی و عرضی، و حذف همانند سازی در 3 راستا (کامل) در حذف خطای همانندسازی در جریان یک کانال پرداخته شده است. تحلیل جریان با حل مستقیم معادلات ناویر استوکس راکم ناپذیر D N S به روش شبه طیفی و الگوریتم گام زمانی متغیر و ساختار چرخشی جمله جابجایی انجام شده است. نتایج حاصله نشان می دهد علی رغم افزایش تعداد دفعات عملیات در الگوریتم های اصلاحی حذف همانندسازی در راستای طولی و عرضی و حذف همانندسازی در 3 راستا (کامل)، استفاده از این روش ها منجر به افزایش دقت و کاهش هزینه های محاسباتی در حدود 33٪ می شود. تفاوت مقادیر شدت آشفتگی میان این 2 روش نیز بیشینه ی 5٫5٪ برآورد می شود، در حالی که اختلاف مقادیر شدت آشفتگی در روش همانندسازی غیرمحذوف و الگوهای فوق بین 7 تا 9 برابر است. درنهایت با توجه به مقادیر آشفتگی محاسبه شده، نتیجه گیری شد که روش های بهبودیافته با حذف همانندسازی از منظر اقتصادی عملکرد مطلوب تری دارند.
کلیدواژگان: شبیه سازی عددی مستقیم D N S، روش شبه طیفی $t a u$، جریان کانال، خطای همانندسازی، گام زمانی متغیر -
صفحات 13-22
در این نوشتار، مقایسه یی بین روش های تحلیلی عددی مشخصه و احتمالاتی انجام، و اثر ناهمگونی تصادفی ویژگی های خاک بر روی ضریب اطمینان پایداری شیروانی های رسی بررسی شده است. ضریب اطمینان پایداری شیب رسی در شرایط زهکشی نشده با لحاظ تغییرپذیری فضایی ویژگی های مقاومتی محاسبه و قابلیت اعتماد به محاسبات ضریب اطمینان با تحلیل های اعتباری بررسی شده است. با استفاده از تئوری میدان تصادفی در ترکیب با روش تفاضل محدود و در چارچوب شبیه سازی های مونت کارلو، توزیع آماری ضریب اطمینان پایداری شیروانی محاسبه و مشاهده شد که میانگین ضریب اطمینان تحلیل های تصادفی در مقایسه با تحلیل های مشخصه، کاهش قابل ملاحظه یی دارد. نتایج به دست آمده از این مطالعه، اهمیت تحلیل های اعتباری را مشخص کرده است و لزوم محاسبات ضریب اطمینان پایداری شیروانی ها را در یک چارچوب احتمالاتی نشان می دهد.
کلیدواژگان: ناهمگونی، میدان تصادفی، تحلیل تصادفی، ضریب اطمینان، ضریب اعتباری، شبیه سازی مونت کارلو -
صفحات 23-27در این نوشتار، الگوریتم تولید خودکار هندسه ی مدل و پارامترهای موردنیاز جهت شبیه سازی بتن با استفاده از روش های گسسته با درنظرگرفتن توزیع تصادفی سنگ دانه ها ارائه شده است. برای تولید یک بافت نامنظم، از دیاگرام V o r o n o i برای شبکه بندی نمونه استفاده شده است. سنگ دانه ها با یک دانه بندی استاندارد و در شکل های متفاوت، به طور تصادفی در بافت بتن توزیع می شوند. باتوجه به اندازه ی المان ها و سنگ دانه ها، هر سنگ دانه می تواند از یک یا چند المان تشکیل شود، که بدین ترتیب امکان بررسی شکستگی سنگ دانه در بتن های با مقاومت بالا فراهم شده است. در این الگوریتم با تعیین ابعاد و تعداد المان، نمونه ی بتنی موردنظر با تمام پارامترهای موردنیاز جهت شبیه سازی بتن ساخته می شود. در پایان، چگونگی حل نمونه ی حاصل با استفاده از روش المان گسسته ارائه شده است.
کلیدواژگان: نمونه ی بتنی، دیاگرام V o r o n o i، سنگ دانه ها، ماتریس اتصال، روش المان گسسته -
صفحات 29-35خاک ها در معرض آلودگی های فراوان در اثر نشت از مخازن نفتی، چاه های نفت، و سدهای باطله هستند، و این آلودگی در خصوصیات فیزیکی و ژئوتکنیکی خاک موثر است. در این تحقیق، با انجام آزمایش های تحکیم ادئومتر، خصوصیات تراکم پذیری خاک های رسی آلوده به سیال آلی بررسی شده است. خاک رس کائولونیت و بنتونیت به منزله ی مصالح پایه و آلاینده های آلی شامل نفت سفید و نفت گاز در درصدهای وزنی صفر، 3، 6 و 9 به آن اضافه شده است. آزمایش ها در دو تراکم نسبی 50٪ و 70٪ انجام شده است. نتایج آزمایش ها نشانگر ازدیاد تورم اولیه ی خاک آلوده به سیال آلی در هنگام اشباع شدن است. به علاوه روغن کاری ذرات خاک به دلیل حضور سیال نفتی نیز موجب افزایش تراکم پذیری خاک خواهد شد، که دلیل آن ویژگی های شیمیایی سیال منفذی، ازدیاد فعالیت کانی های رسی و افزایش ضخامت لایه ی جذب سطحی است.
کلیدواژگان: خاک آلوده، سیال نفتی، تحکیم یک بعدی، تراکم نسبی، ویژگی های شیمیایی -
صفحات 37-48در این تحقیق روش اصلاح شده ی تفاضل ماتریسی با استفاده از توابع زمان بسامد مربعی برای شناسایی آسیب در پایه ی پل ها و محل آن پیشنهاد شده است. برای تایید روش پیشنهادی از مدل پل قطور به طول 448 متر و برای پردازش سیگنال های پاسخ لرزه یی پل ها از توابع زمان بسامد شامل طیف تبدیل فوریه ی زمان کوتاه، طیف تبدیل موجک، توزیع های ویگنر ویل، چوی ویلیامز، شبه ویگنر ویل هموارشده، و تداخلی کاهش یافته استفاده شده است و تابع زمان بسامد تداخلی کاهش یافته به منزله ی تابع مطلوب برای پردازش سیگنال های پاسخ لرزه یی پل ها شناسایی شده است. پلان های زمان بسامد محاسبه و مشخصه های دینامیکی سیستم استخراج شده است. روش پیشنهادی براساس پردازش سیگنال های لرزه یی خروجی پایه ریزی شده است، و لذا از مزیت های این روش بی نیازی به تعریف مدل تحلیلی و اندازه گیری نیروی ورودی و همچنین عدم نیاز به تحریک اجباری هارمونیک پل پس از وقوع زلزله برای تشخیص آسیب لرزه یی است.
کلیدواژگان: شناسایی مشخصه، تشخیص آسیب، توابع زمان، بسامد، پل، تفاضل ماتریسی اصلاح شده، همبستگی، فاصله ی مربعی کمینه -
صفحات 49-59در این نوشتار، مشخصات جریان به دست آمده از نتایج مدل آزمایشگاهی شکست سد با درنظرگرفتن پایه ی پل استوانه یی شکل واقع شده در میان کانال مستطیلی شکل بررسی شده است. تغییرات سطح آب توسط حسگرهای آلتراسونیک و دوربین با سرعت بالا و مولفه های سرعت با دستگاه سرعت سنج (A D V) برداشت شده است. مشاهدات جریان، نشان دهنده ی تغییر مسیر جریان برای عبور از پایه ی پل و ایجاد جریان های عرضی در پایین دست است؛ به طوری که پشت پایه ی پل تا فاصله یی به اندازه ی 13 برابر قطر آن، ناحیه ی برخاستگی جریان مشاهده می شود. همچنین مقایسه ی تغییرات سطح آب به دست آمده با مقادیر حل تحلیلی ریتر، بیانگر وجود درصد اختلاف های کمتر (4 الی 12 درصد) در محدوده ی بین مخزن سد تا محل پایه ی پل است؛ در حالی که در فاصله ی پشت پایه ی پل تا 13 برابر قطر آن به علت وجود گردابه های برخاستگی، مقادیر آزمایشگاهی از حل تحلیلی ریتر فاصله می گیرد.
کلیدواژگان: شکست سد، پایه ی پل، حسگرهای التراسونیک، دوربین C C D، دستگاه A D V -
صفحات 61-70در این تحقیق، ریزساختار بتن به کمک یک مدل اجزای محدود سه بعدی شبیه سازی شده است. این مدل دارای دو فاز ملات و سنگدانه های درشت است. شکل سنگدانه ها به صورت کروی و رفتار آن ها به صورت ارتجاعی خطی در نظر گرفته شده است. از منحنی دانه بندی فولر جهت بیان توزیع سنگدانه ها استفاده می شود. رفتار ملات در ساختار خمیری خسارت بیان شده است. در ابتدا، نحوه ی پیاده سازی مدل خمیری خسارت ارائه و در گام بعد، صحت پیاده سازی مدل خمیری خسارت برای یک المان بررسی می شود. در بخش پایانی اثر درصد حجمی سنگدانه، اندازه ی بزرگ ترین سنگدانه و مدول ارتجاعی سنگدانه بر روی مقاومت فشاری مشخصه ی بتن مطالعه و نتایج با دستاوردهای سایر پژوهشگران مقایسه شده است.
کلیدواژگان: ریزساختار، مقاومت فشاری بتن، روش اجزای محدود، خمیری، خسارت -
صفحات 71-76رفتار پی رادیه شمع در خاک رس نرم به دلیل پاره یی نگرانی ها نظیرظرفیت باربری ناکافی خاک و احتمال نشست زیاد پی، مسئله یی است که در سالیان اخیر مورد توجه قرار گرفته است. سهم شمع ها نسبت به رادیه در تحمل بار وارد بر پی مسئله یی مهم است. در نوشتار حاضر با مدل سازی عددی به بررسی تاثیر مشخصات هندسی پی رادیه شمع بر رفتار آن در خاک رس نرم پرداخته و تحلیل ها در قالب مطالعات پارامتریک انجام شده است. تحلیل های انجام شده مشخص کرد که نحوه ی آرایش شمع ها در زیر رادیه تاثیر بسزایی بر عملکرد پی دارد. افزایش طول و قطر شمع ها موجب کاهش نشست پی و بالارفتن سهم شمع ها در تحمل بار وارده می شود. افزایش ضخامت رادیه منجر به کاهش نشست متوسط و تفاضلی پی می شود، اما تغییر ضخامت رادیه اثری بر سهم باربری رادیه و یا شمع ها ندارد.
کلیدواژگان: رادیه، شمع، رس نرم، مدل عددی، مطالعات پارامتریک -
صفحات 77-85کاهش پاسخ سازه با اعمال نیروی کنترلی را می توان مطابق الزامات آیین نامه های طراحی براساس عملکرد با استفاده از روش های تحلیل غیرخطی بهینه ساخت. در این نوشتار با استفاده از سیستم کنترل فعال، رفتار سازه در جهت رسیدن به عملکرد مورد نظر بهبود بخشیده می شود. نیروی کنترل با الگوریتم بهینه ی لحظه یی به قاب های 4، 8 و 12 طبقه با فرض رفتار غیرخطی سازه اعمال می شود، تا جابه جایی که سازه را در محدوده ی جابه جایی های حاصل از تحلیل استاتیکی غیرخطی در سطح عملکردی هدف نگه دارد. درواقع، نقطه ی عملکرد سازه معیار مقدار نیروی کنترلی اعمالی است. نتایج نشان می دهند الگوریتم بهینه ی لحظه یی قادر به کاهش زیاد جابه جایی سازه در نسبت 9E+1c/r= است؛ که این نسبت به نیروی کنترل قابل ملاحظه یی نیاز دارد. برای سازه های مورد مطالعه نسبت c/r بین 9E1 تا 9E5، سازه را در محدوده ی عملکردی I O حفظ می کند. این نسبت برای سازه های مشابه نیز قابل استفاده است.
کلیدواژگان: کنترل فعال، طراحی بر اساس عملکرد، جابه جایی هدف، الگوریتم کنترل بهینه ی لحظه یی، سطح عملکردی -
صفحات 87-98بسیاری از سنگ های موجود در طبیعت، سخت و شکننده هستند. لذا استفاده از مدل ساختاری مناسبی که این رفتار شکننده را به خصوص در تونل های عمیق شبیه سازی کند، بسیار ضروری است. در این تحقیق، رفتار گسیختگی این سنگ ها با استفاده از روش های حل بسته و عددی (اجزاء محدود و تفاضل محدود) با پیش بینی مقادیر جابجایی دیواره و ناحیه ی خمیری در دو تونل عمیق با مقطع دایره یی و نعل اسبی بررسی می شود. نتایج در حالت رفتار خمیری کامل بین روش های عددی و حل بسته همخوانی بسیار خوبی نشان می دهد؛ ولی با فرض رفتار شکننده، بین روش های اجزاء محدود، تفاضل محدود و حل بسته، بعضا اختلافات زیادی دیده می شود. تحلیل های انجام شده حاکی از آن است که با فرض رفتار شکننده، مقادیر تغییرمکان و ناحیه ی خمیری شدیدا به پارامترهای مقاومتی پس از اوج حساس اند و با کاهش مقاومت پسماند به شدت افزایش می یابند.
کلیدواژگان: تونل عمیق، رفتار شکننده، سنگ سخت، P h a s e^2، F l a c -
صفحات 99-106یکی از مشکلات عمده در سیستم های کنترل فعال، وجود تاخیر زمانی است. تاخیر زمانی می تواند تاثیر قابل ملاحظه یی در عملکرد سازه در هنگام زلزله داشته باشد، به طوری که به ناپایداری کلی سازه نیز بیانجامد. در این نوشتار، ابتدا سازه ی موردنظر نسبت به تاخیر زمانی تحلیل حساسیت شده و سپس مقدار بحرانی آن تعیین می شود.
بررسی ها نشان دادند تاخیر زمانی بحرانی سازه به نوع زلزله ی ورودی حساس نیست و فقط به خصوصیات سازه و سیستم کنترلی ارتباط دارد. لذا برای مهار این پدیده و افزایش بازدهی سیستم کنترل، از ایده ی افزایش ظرفیت میراکنندگی سازه استفاده شد. افزایش میرایی به گونه یی انجام شد که فقط ظرفیت جذب انرژی سازه تغییر کند. الگوریتم مورد استفاده، الگوریتم کنترل گسسته ی بهینه ی آنی است. ارزیابی ها نشان داده اند افزایش میرایی می تواند تا اندازه ی زیادی حساسیت سازه و سیستم کنترل را نسبت به تاخیر زمانی محدود کند. به این شیوه می توان سازه های کنترلی را که با گذشت زمان کارآیی خود را از دست داده اند، مجددا احیا کرد.
کلیدواژگان: کنترل فعال، نیروی کنترلی، تاخیر زمانی، میرایی، الگوریتم کنترل -
صفحات 107-116در این مطالعه، الگوریتم بهینه سازی P S O با مدل جامع شبیه ساز M O D S I M به منظور حل مسئله ی تخصیص منابع آب حوضه ی رودخانه ی اترک تلفیق شده است. تابع هدف مدل منتج P S O-M O D S I M، با مجموعه یی از متغیرهای تصمیم بهره برداری و طراحی مخازن سدها، براساس بهبود شرایط تخصیص آب به نیازها نسبت به شرایط پیش از توسعه ی سیستم کالیبره شد. در ادامه، به منظور بهبود سرعت اجرا و تنظیم موثر پارامترهای مدل، نسبت به ساخت یک فرامدل جایگزین از نوع ماشین های بردار پشتیبان به جای مدل M O D S I M که در آن ارزیابی تابع هدف با استفاده از مدل تقریبی جایگزین انجام می شود، اقدام شد. با توسعه ی مدل جایگزین (P S O-M O D S I M$s i m$S V M) حل این مسئله ی حوضه یی با نتایج حاصل از مدل اصلی P S O-M O D S I M و یک مدل مبتنی بر فرامدل شبکه های عصبی مصنوعی (A N N) مقایسه شد. نتایج حاصل از فرامدل سازی با روش S V M از منظر کیفیت جواب ها و بار محاسباتی، بیانگر عملکرد مطلوب این روش و برتری آن بر فرامدل از نوع A N N است.
کلیدواژگان: تخصیص آب، حوضه ی رودخانه، بهینه سازی جایگزین، ماشین های بردار پشتیبان -
صفحات 117-126پل های خاکی -فولادی با ترکیب صفحات فولادی موج دار و خاک، سیستم جدیدی را برای پل سازی ارائه کرده اند. در تحلیل این نوع پل ها عموما تنش ها و تغییرشکل های ایجادشده در لایه های خاک و دیواره ی سازه تحت اثر بارهای اجرایی و بهره برداری مورد توجه قرار می گیرد. در این نوشتار، پس از طراحی یک نمونه پل خاکی فولادی مطابق آئین نامه ی کانادا با استفاده از نرم افزار اجزاء محدود P L A X I S، بارهای اجرائی بر پل موردنظر اعمال شده و نتایج باربری با پروفیل پیشنهادی آئین نامه یی مقایسه شده اند. این مقایسه نشان می دهد که کنترل تنش و تغییرشکل مذکور با پروفیل پیشنهادی آئین نامه یی جوابگوی تمامی ضوابط طراحی است، ولیکن دقت در بارگذاری و تحلیل گام به گام، استفاده از مدل های رفتاری غیرخطی، و درنظرگرفتن لغزش بین پوشش خاکی و جداره می تواند منجر به انتخاب پروفیل اقتصادی تری در مقایسه با آئین نامه های طراحی مبتنی بر روش های کشسان شود.
کلیدواژگان: پل های خاکی فولادی، تحلیل استاتیکی، روش اجزاء محدود، بارگذاری کامیون
-
Pages 3-12T h e o b j e c t i v e o f t h i s w o r k i s t o p e r f o r m a l i a s e d a n d d e-a l i a s e d d i r e c t n u m e r i c a l s i m u l a t i o n o f a t u r b u l e n t c h a n n e l f l o w w i t h a v a r i a b l e t i m e-s t e p p i n g a l g o r i t h m. T h e t i m e a d v a n c e m e n t i s c a r r i e d o u t b y a n i m p l i c i t 3r d o r d e r b a c k w a r d d i f f e r e n t s c h e m e f o r l i n e a r t e r m s، a n d e x p l i c i t f o r w a r d E u l e r f o r n o n l i n e a r t e r m s (S B D F3)، i n T=300s (a b o u t 12000 t i m e s t e p). A $t a u$ p s e u d o s p e c t r a l m e t h o d i s e m p l o y e d o n 128$t i m e s$65$t i m e s$128 c o m p u t a t i o n a l g r i d s i n $O m e g a$= [0. 4$p i$] $t i m e s$ [-1. 1] $t i m e s$ [0. 2$p i$]. T h e F o u r i e r s e r i e s i n s t r e a m-w i s e a n d s p a n-w i s e d i r e c t i o n s a n d C h e b y c h e v p o l y n o m i a l e x p a n s i o n i n a n o r m a l d i r e c t i o n a r e e m p l o y e d. T h e f r i c t i o n R e y n o l d s n u m b e r i s $R e_ {t a u} $=175 a n d i n i t i a l v i s c o s i t y i s s e t t o $n u$=1/4000. T h e e f f e c t o f a l i a s i n g e r r o r s i s e x a m i n e d i n t h e N o D e a l i a s i n g s c h e m e. A l s o، t h e p r o d u c t''s a l i a s i n g e r r o r s a r e e l i m i n a t e d b y t r u n c a t i n g t h e F o u r i e r c o e f f i c i e n t s o f t h e v e l o c i t y a n d t h e p r o d u c t، u s i n g t h e s o-c a l l e d ``2/3 r u l e''، i n t h r e e s c h e m e s o f: D e a l i a s Y، D e a l i a s X Z a n d D e a l i a s X Y Z. A c o m p a r i s o n o f t u r b u l e n c e i n t e n s i t i e s i s m a d e o v e r t h e i n t e r v a l f r o m $T_ {0} =100$ t o $T_ {1} =300$. T h e r e s u l t s o f t u r b u l e n c e s t a t i s t i c s، w a l l-s h e a r v e l o c i t y، s t a n d a r d d e v i a t i o n، o v e r a l l C P U t i m e a n d t o t a l n o r m a l i z e d e n e r g y o f i n s t a n t a n e o u s v e l o c i t i e s h a v e a l s o b e e n c o l l e c t e d. A l s o، t h e r e s u l t s o f t h e m e a n v e l o c i t y p r o f i l e a r e c o m p a r e d t o t h e s c a l i n g l a w s f o r v i s c o u s a n d i n e r t i a l l a y e r s i n a l l s c h e m e s. T h e p r e s e n t r e s u l t s o f t u r b u l e n c e s t r e n g t h s h o w s i g n i f i c a n t d i s c r e p a n c i e s b e t w e e n t h e r e s u l t s o f N o D e a l i a s i n g a n d D e a l i a s Y s c h e m e s a n d t h e D e a l i a s X Z a n d D e a l i a s X Y Z s c h e m e s a t t h e o t h e r s i d e. A l s o m t h e r e a r e n o c o n s i d e r a b l e d i f f e r e n c e s b e t w e e n D e a l i a s X Z a n d D e a l i a s X Y Z r e s u l t s. T h e d i f f e r e n c e s b e t w e e n c o m p u t e d v a l u e s o f t u r b u l e n t q u a n t i t i e s v a r y f r o m 96% i n t o t a l n o r m a l i z e d e n e r g y t o 3% i n d i s s i p a t i o n. C o m p a r e d t o t h e N o D e a l i a s i n g s c h e m e، m o d i f i e d D e a l i a s X Y Z r e d u c e s t h e c o m p u t a t i o n a l t i m e a r o u n d 33%. T h e r e f o r e، f r o m a p r a c t i c a l s t a n d p o i n t، t h e c h o i c e b e t w e e n t h e c o s t l y، o r d i n a r y، a n d d e-a l i a s e d r o t a t i o n a l f o r m r e d u c e s t o t h e l a s t o n e، i f t h e f i n a l d e c i s i o n r e s t s o n e c o n o m i c s.Keywords: Direct numerical simulation (DNS), pseudo spectral method, channel flow, aliasing errors, variable time, stepping algorithm
-
Pages 13-22
Slope stability analysis is a geotechnical engineering problem characterized by many sources of uncertainty. Some of the uncertainties are related to the variability of soil properties involved in the analysis. The influence of spatial variation on the mean of the safety factor is discussed in both deterministic and stochastic numerical analysis contexts. A numerical procedure for a probabilistic slope stability analysis, based on Monte Carlo simulation, which considers the spatial variability of soil properties, is presented to assess the influence of randomly distributed undrained shear strength. In the proposed method, the commercially available finite difference numerical code, FLAC 5.0, is merged with random field theory generated using the Cholesky decomposition technique. For each realization, the random field of $c_{u}$ is generated by Cholesky decomposition at the center of the element level, and the in-situ stresses in the slope are generated by applying gravity loading. This method considers the spatial variability of soil properties, and applies a strength reduction method to estimate the safety factor of a slope. Considering the case of a 5.0 m high cohesive soil slope of $54.5^0$ (a range of coefficients of variation (COV) from 10% to 50% in cohesion values), parametric studies are performed to study the effect of stochastic soil cohesion on the statistics of the safety factor, in comparison to the deterministic solution available for the uniformly constant cohesion property. Effects of the statistical range of variability of $c_{u}$ on the mean of the safety factor are also presented and give an important insight into their effect on reliability analyses. The results obtained in this study are useful to understand the role of cohesion variations in slope stability analysis under different slope conditions and material properties. The coefficient of variation of undrained shear strength was proven to have a significant effect on the reliability of safety factor calculations. The observations made from this study help to explain the requirement for slope stability calculations in a probabilistic framework. It is further demonstrated that the variability of soil properties translates into a substantial reduction in the safety factor (in an average sense), compared to the corresponding deterministic (homogene deterministic (homogeneous soil) case.
-
Pages 23-27Concrete is perhaps the most available manufactured material. The low cost, wide availability, ease of use and high durability of concrete has led to its continually increasing usage. It can be a hastily prepared, low-grade mixture for simple applications, or can be a firmly controlled, engineering material for high-performance structures. Complex physical and chemical interactions exist in cement, which plays a main role in the properties of concrete. This complicated structure leads to a complexity in fracture processes. Considering this reason, most mechanical properties of concrete and studies on its behavior are based on experimental results. Since experiments require time and money, providing mathematical models to simulate the behavior of concrete is necessary. In general, modeling fracture and damage within concrete, and other quasi-brittle materials, has been classified as either continuum or discrete approaches. Continuum models provide an average description of material behavior for a representative volume element. Because the width of the fracture process zone (FPZ) in concrete can be sizeable (roughly several times the maximum aggregate size), simulation of concrete fracture at meso-scale, with continuum approaches, is not suitable. Use of discrete micromechanical models is motivated by the need for fundamental knowledge, to improve material behavior. If the material structure (e.g at micro/meso scale of observation) is explicitly represented, the models provide a direct way for studying crack patterns; mechanisms of softening in post-peak branches and size effect/scaling phenomena. In this paper, two-dimensional geometrical models for concrete are generated, taking the random distribution of aggregates at mesoscale into consideration. The generation procedure is based upon the Voronoi diagram method. The aggregate particles are constructed by several polygons and then placed into the concrete model, in such a way that there is no intersection between them. In this method, simulation of the fracture of aggregate in high strength concrete is feasible. The generated model can be used for modeling concrete with the discrete method. Finally, the model analyzed using the discrete elemen method.Keywords: Concrete specimen, Voronoi diagram, aggregates, random distribution, connectivity matrix, discrete element method
-
Pages 29-35Soil is generally subjected to different contaminants that are being used and excreted into the environment. Leakage of these materials can influence its physical and mechanical characteristics. A review of technical literature in this field shows some contradictions between the results of different research, which reveals the impact of soil and contaminant type on the behavior of contaminated clayey soil. In order to investigate the effects of organic contamination on the physical and mechanical characteristics of clayey soil, in the present study, a series of index tests, besides one-dimensional consolidation (Oedometer tests), has been performed. Two types of clay, Kaolinite and Bentonite, with low and high plastic indexes are used as the base soils. Also kerosene and gasoil are applied as the organic contaminant in four different percentages of 0, 3, 6 and 9 to the base soil. All samples were kept for seven days in polyethylene bags to perform all chemical reactions before the tests. After that, water is used to saturate all samples normally within 48 hours. Results of physical tests reveal that Atterberg limits increase, by an increase in degree of contamination, for both considered soils. Also, maximum dry density increases with enhancement of the degree of contamination. However, according to the Proctor compaction tests, the optimum moisture content reduces with degree of contamination. This is illustrated due to the lubrication of grains during contamination with organic material. Consolidation tests are conducted on samples in two distinct relative compactions of 50% and 70%, and the results are investigated in e-Log (p) coordinates. Based on the results of one dimensional consolidation tests, compressibility coefficient increases, with an increase in contaminant percentage for both soils. Also, the effects of contamination on the consolidation behavior of clay increase with an enhancement in relative compaction. As a result, it is negligible in soil with low relative compaction. Using the results of tests, an empirical relation is suggested between the compressibility coefficient of contaminated soil and its relative compaction and liquid limit. Comparison of the suggested equation with ordinary relations in soil mechanics texts shows a great difference, which proves the necessity for new suggested forms in contaminated soil. Finally, it was concluded that chemical and organic contaminants, like kerosene and gas oil, which are not soluble in water, are more effective in changing the physical characteristics of soil, like Atterberg limits and relative compaction. However, compared to the first category, soluble organic contaminants are more effective on the mechanical and consolidation characteristics of clayey soil.Keywords: Bridge system identification, damage detection, time, frequency representation, modified matrix subtraction
-
Pages 37-48In this research, a new modified matrix subtraction method is proposed using square time-frequency representations to detect bridge pier damage and its location. The new proposed method is confirmed, using the seismic response of the 448 meter long Ghotour Bridge model, and compared with correlation and least square distance methods. Time frequency plans of damaged and non-damaged bridge matrix elements are calculated using reduced interference distribution. The difference matrix is calculated by subtraction of corresponding matrix elements. The possibility of an existing damage index is calculated by summation of all elements of difference matrices and then normalizing them with a maximum value. Linear time history analyses and earthquake acceleration records of San Fernando, Loma Prieta and Northridge earthquakes are used for seismic response analysis. In previous response analyses of the Ghotour Bridge, as stated in reference [23], the third pier from the left is considered more vulnerable to seismic damage. Therefore, for the purpose of this study, a reduction of thirty percent in the stiffness of that pier is considered, to simulate the seismic damage in the damaged analytical model. It is shown that the proposed method could satisfactorily identify the damage location. For the seismic response of the top of the bridge piers, the maximum error in locating the damage is 6 percent, while, for the seismic response of the middle of the bridge piers, it is 14 percent.The time-frequency representations used include: the short time Fourier transform spectrum, wavelet transform spectrum, Wigner-Ville Distribution, Choi-Williams distribution, smoothed pseudo Wigner-Ville distribution and reduced interference distribution, which are finally identified as optimal performance time-frequency epresentation for bridge seismic response signal processing. Reduced interference distribution and Wigner-Ville distribution are both in the Cohen class, but reduced interference distribution methods are more appropriate for processing seismic bridge transient nonstationary response signals. Time-frequency planes have been calculated and dynamic specifications of the system have been estimated. The proposed algorithm is a seismic output-only method. Therefore, it has the advantage of not needing to define the bridge analytical model and measuring seismic input loading, and, also, not needing to use harmonic forced vibration analysis after earthquake occurrence for bridge seismic damage detection, as is general in some other methods.Keywords: Bridge system identification, damage detection, time, frequency representation, modified matrix subtraction
-
Pages 49-59In this research, a new modified matrix subtraction method is proposed using square time-frequency representations to detect bridge pier damage and its location. The new proposed method is confirmed, using the seismic response of the 448 meter long Ghotour Bridge model, and compared with correlation and least square distance methods. Time frequency plans of damaged and non-damaged bridge matrix elements are calculated using reduced interference distribution. The difference matrix is calculated by subtraction of corresponding matrix elements. The possibility of an existing damage index is calculated by summation of all elements of difference matrices and then normalizing them with a maximum value. Linear time history analyses and earthquake acceleration records of San Fernando, Loma Prieta and Northridge earthquakes are used for seismic response analysis. In previous response analyses of the Ghotour Bridge, as stated in reference [23], the third pier from the left is considered more vulnerable to seismic damage. Therefore, for the purpose of this study, a reduction of thirty percent in the stiffness of that pier is considered, to simulate the seismic damage in the damaged analytical model. It is shown that the proposed method could satisfactorily identify the damage location. For the seismic response of the top of the bridge piers, the maximum error in locating the damage is 6 percent, while, for the seismic response of the middle of the bridge piers, it is 14 percent. The time-frequency representations used include: the short time Fourier transform spectrum, wavelet transform spectrum, Wigner-Ville Distribution, Choi-Williams distribution, smoothed pseudo Wigner-Ville distribution and reduced interference distribution, which are finally identified as optimal performance time-frequency representation for bridge seismic response signal processing. Reduced interference distribution and Wigner-Ville distribution are both in the Cohen class, but reduced interference distribution methods are more appropriate for processing seismic bridge transient nonstationary response signals. Time-frequency planes have been calculated and dynamic specifications of the system have been estimated. The proposed algorithm is a seismic output-only method. Therefore, it has the advantage of not needing to define the bridge analytical model and measuring seismic input loading, and, also, not needing to use harmonic forced vibration analysis after earthquake occurrence for bridge seismic damage detection, as is general in some other methods.Keywords: Bridge system identification, damage detection, time, frequency representation, modified matri
-
Pages 61-70In the present paper, the detailed mesostructure of concrete is geometrically generated and its compressive strength is numerically estimated using the 3D finite element method. The models contain two phases of mortar and coarse aggregates. The FE models of concrete are cubic in shape, with a side length of 80 mm. Aggregates are assumed to be spherical and behave in a linear elastic manner. The famous Fuller formula is utilized for the aggregate grading curve, and the simple sequential inhibition (SSI) technique is employed to fill the concrete cubes with the particles. Only aggregates bigger than 4.75 mm in diameter (gravel) are modeled, i.e., the particles smaller than 4.75 mm in diameter (sand) are not considered individually and assumed to be part of the homogenized nonlinear cement paste. A modified version of the plastic-damage model, proposed by Lee and Fenves [J. Lee, G.L. Fenves, International Journal for Numerical Methods in Engineering 50 (2001) 487-506], has been adopted to simulate the inelastic response of the mortar. This constitutive model incorporates two independent hardening variables, namely; equivalent tensile and compressive plastic strains, and, thus, is capable of tracing damage evolution due to both tensile cracking and compressive crushing. In the first stage, the numerical implementation of the plastic-damage model is presented and then its validity is examined in a 3D FE element. Next, the effects of aggregate volume fraction, aggregate maximum diameter, and aggregate elastic modulus on concrete compressive strength are studied. It is shown that: (1) compressive strength remains constant for specimens with aggregate volume fractions of up to 50%, and then increases significantly with grain content, (2) for the range of aggregate volume fractions studied in this paper, the maximum aggregate size has little influence on compressive strength, and (3) any increase in the elastic modulus of aggregates accentuates stress concentration near the aggregates, and, thus, reduces the compressive strength of concrete samples. Finally, the results are satisfactorily compared with those presented by other researchers.Keywords: Mesostructure, compressive strength, concrete, Finite Element method, plastic, damage
-
Pages 71-76Combined piled-raft foundations (CPRF), is an innovative idea proposed since the 1970s. Adding piles to raft foundations contributes to settlement reduction and an increase in the bearing capacity of the raft. Generally speaking, CPRF helps in considerably reducing the number of piles and pile lengths, improving serviceability in both total and differential settlements, and minimizing tilt and instability probabilities. Reduction in the number and length of piles also contributes to considerable saving in construction cost. The applicability of the construction of piled-raft foundations on soft clay has been under greater attention in recent years. Lack of sufficient stiffness for soil and, thus, excessive settlement, allowing mobilization of pile loads, is a major concern in this regard. This paper presents a numerical investigation into the effects of geometrical variations of piled-raft components on soft clay, as a parametric study. A 3D FEM numerical model is developed using ABAQUS. The model was calibrated by comparing physical and numerical modeling of other researchers. An elaborate parametric study was then conducted to evaluate the effect of various parameters on CPRF applicability on soft clay. In this study, attention was focused on the drained (long-term) response of a piled-raft resting on a soft clay layer, so the clay was idealized using drained shear strength parameters, c' and phi Thus, consolidation effects were neglected. In all analyses vertical uniform loading was applied to the top of the raft surface. The results indicate that pile arrangements under the raft play an important role in piled-raft performance. An increase in pile length and diameter results in settlement reduction and an increase in the pile contribution ratio. An increase in raft thickness reduces differential settlements, but has no considerable effect on the pile and/or raft bearing ratio. Overall, the construction of pile-rafts on soft clay is promising, as long as the total settlement of the structure is not imposing restrictions, such as the common 20 mm allowable settlement. But, instead, if the allowable settlements are in the range of 100 to 150 mm, for example, in the case of rigid steel tanks, the method will be applicable, with considerable cost and co tanks, the method will be applicable, with considerable cost and construction time savings.Keywords: Piled, raft, soft clay, numerical modeling, parametric study
-
Pages 77-85During recent earthquakes, the behavior of structures has demonstrated that the main reason for building collapse is the big deformation during strong seismic events. Therefore, the structural designers would like to reduce the lateral displacement of structures under a specified level. Using the shear walls and bracing and increasing the members size are the common ways which can protect the buildings against the earthquake or wind load. In two recent decades, a new method has been invented which is called Active Structural Control. Applying the external force, when it is needed; not only can reduce the response displacement but also is able to make a high reduction in response velocity and acceleration. Predict the building behavior during an earthquake event has an important role in order to have a safe structure. The structure engineers usually use two methods to compute the structure response: i) Dynamic Analysis or ii) Nonlinear Static Analysis. Since the dynamic analysis depends on several items such as frequency content of earthquake record, the most structure specialists prefer to employ the nonlinear static analysis (pushover analysis) to compute the target displacement in IO, LS and CP levels. The essential goal of this research is to create a combination between active structural control concept and nonlinear static analysis. Three reinforced concrete buildings with 4, 8 and 12 floors are selected. All buildings have three bays and are analyzed using the pushover method and for each case a performance level is found. The main purpose of using active structural control equipment in the selected buildings is to keep the building performance in that limit which is found at the analytical part. An important aspect of this investigation is the incorporation of the possibility of inelastic deformation of the members into the control algorithm. In analytic investigation, in order to solve the control formula, the Instantaneous Optimal Control algorithm is employed which is conducted to introduce c/r'' ration. This ratio shows the importance of the amount of control force against the importance of reduction in structural response. The results show that the best c/r ratio for the selected structures is between 1$times$$10^9$ and 5$times$$10^9$ which can keep the structures in IO performance level. Generally, the optimal amount of c/r ratio to get the target performance level depends on different items such as the primary strength of structure (the initial performance level) and probable seismic hazard, therefore; it seems that the hazard risk analyses is necessary as an additional analysis.Keywords: Active control, performance based design, target displacement, instantaneous optimal control algorithm, performance level
-
Pages 87-98Estimation of stress and deformation in rock mass around deep tunnels is of great concern in geotechnical engineering, rock mechanics and mining engineering. Due to large stresses in these tunnels, the surrounding rock mass may experience failure in some points. Hence, in these tunnels, the failure process is of great importance. There are different rheological models to simulate the failure behavior of rocks, including perfect plastic, strain softening and brittle models. Since many rocks in nature are hard, with brittle behavior, it is required to select a suitable mechanical model for simulating such behavior. In this research, the post-peak behavior of hard rocks is studied using the numerical methods of FEM and FDM, with $Phase^2$ and FLAC codes. The effect of brittle behavior is evaluated on deformation and the extent of the plastic zone in deep tunnels excavated in hard rock. Then, these results are compared with those obtained based on perfect plastic assumption. Due to the fact that the rock strength after failure decreases to a residual level, different residual values are considered, in order to show the importance of selecting the correct residual parameters. Firstly, a circular tunnel is analyzed by both numerical and closed form methods and the obtained stresses and displacements are compared. In addition to tunnel wall deformation, the radius of the plastic zone is evaluated. In this circular tunnel, deformation of the tunnel wall with brittle assumption is obtained twice greater than the results based on perfect plastic assumption. This ratio is decreased at higher radial distances. Moreover, radial stresses in brittle condition are lower than values obtained in perfect plastic condition. Also, tangential stresses in the yielded zone are lower than those stresses in the elastic region, because of stress redistribution. In this paper, a parametric study is accomplished on a horse-shoe tunnel, with both Mohr-Coulomb and Hoek-Brown failure envelopes. The rock is assumed homogenous and isotropic with a non-associated flow rule. The Cai et. al. approach was adopted to estimate the peak and residual strength values of the rock mass. The importance of selecting correct residual strength is again presented. This study illustrates that numerical and closed form methods submit the same results in perfect plastic assumption, but different results in brittle assumption. Specifically, the obtained deformations and plastic zone are highly dependent on post-peak residual parameters. The strain softening behavior, as an intermediate condition between perfect plastic and brittle conditions, was considered as well. In this research, it was observed that the available numerical codes usually exhibit non-symmetric results in very low residual parameters, even for symmetric loading. This can be attributed to the weak points of these numerical codes.Keywords: keywords{Deep tunnel, brittle behaviour, hard rock, $PHASE^2$, FLAc
-
Pages 99-106Time delay is a major problem in active structural control, which may lead to an inactive control performance or even result in instability of the whole structure. Sensitivity of structures to time delay may refer to characteristics of the building, the control system, and/or external loadings. If a control system is designed based on the ideal parameters of the building and control mechanism, it may be inefficient, because it does not take note of the practical problems encountered in reality. But, if a control system is perfectly designed, it may suffer from some changes in structure and/or control parameters. Hence, time delay may always affect the performance of a controlled structure. In this paper, in order to evaluate the influence of this effect, sensitivity analysis of an eight-story model building, with respect to time delay, is carried out. The control system is an active mass damper/driver (AMD) installed on the top floor. The control algorithm is the discrete instantaneous optimal control (DIOC) method, previously presented by the second author. This algorithm is powered using stable matrices to overcome the time step problems that the instantaneous optimal control methods confront. First, a critical time delay is recognized, in which the controlled structure is unstable. By changing input excitation (strong ground motion), the effect of input seismic load on critical time delay is examined. Using different earthquake loads: El Centro, Parkfield, and San Fernando, results in the same critical time delay for the mentioned buildings. Then, the stiffness of the building stories are increased without changing the mass of the floors. It is recognized that there is a relation between critical time delay and the first period of the controlled building. The proposed coefficient is about $frac{1}{7}$ to $frac{1}{6}$ of the fundamental period of the controlled building. Second, by introducing passive dampers in such a manner that the capacity of the energy absorption of the building is increased without increasing its stiffness, critical time delay is evaluated, again. There is no change in critical time delay, but extensive results show that increasing damping to a large extent may alleviate the sensitivity of the controlled building to time delay. In other words, increasing the damping of the first mode of the controlled building may extensively increase the efficiency of the control system. Hence, this may be a good solution for older active control buildings, which suffer from time delay effects.Keywords: Active control, time delay, passive damping, control algorithm, stability, discrete instantaneous optimal control, DIOC
-
Pages 107-116Basin-scale river basin simulation models are vital tools in assessing and quantifying large-scale water resource systems performances under different scenarios. However, simulation models are not capable of doing multiperiod optimization for optimal planning and management purposes. Linking river basin simulation models to multiperiod optimization algorithms could be a remedy to such shortcoming. Nevertheless, large-scale river basin models, when used as the simulation model of an optimization algorithm, may face difficulties, in terms of the needed time for combined optimization-simulation models to run. The more hydrologic and socio-economic processes to be addressed in the simulation model, the more the computational burden of the combined model would be. Meta-modeling is one of the useful approaches to dealing with this difficulty, where a fast-running approximate model, called meta-model, replaces the exact simulation model. This study presents application of the particle swarm optimization (PSO) algorithm, linked to the MODSIM decision support system (DSS), as a river basin simulation model, resulting in the integrated PSO-MODSIM optimization-simulation model, for solving basin-scale water allocation problems. The developed model was applied to the problem of optimal water resource development plans, as well as optimal water allocation of the Atrak River Basin, as a real case study in north-east Iran. Where there are serious competitions between upstream and downstream provinces of the basin in utilization of available water resources. The objective function of the model, which consists of design variables; i.e., reservoir capacities, and operational; i.e., priority numbers of reservoir target levels, was calibrated based on improving water allocation conditions compared to a benchmark scenario, in which development projects were not implemented. Since the integrated model was time consuming, support vector machines (SVMs) were used as a meta-model to develop the PSO-MODSIM$sim$ SVM model, in which an SVM-based surrogate model replaces MODSIM. The performance of the PSO-MODSIM, PSO-MODSIM$sim$ SVM and another model, with a meta-model of artificial neural network (ANNs) type, were analyzed and compared through their application to optimal water allocation and planning water resource development projects in the basin. Before using SVMs as a meta-model, we first showed how well they perform in approximating benchmark ultidimensional mathematical functions. Then, basin-scale optimal water allocation problems were solved using surrogate optimization techniques. Both SVM and ANN models were able to approximate and represent the MODSIM DSS reasonably. However, the results show that SVM-based surrogate optimization has performed satisfactorily in terms of both the quality of the solutions and the saving in computational burden compared to the model employing ANN as the meta-model.Keywords: Water allocation, river basin, surrogate optimization, support vector machines
-
Pages 117-126Soil-steel structures are composite structures made by steel rings and a soil envelope, which are constructed in spans of, maximum, thirty meters. Soil-steel interaction causes flexible steel plates to interact with the urrounding backfill for suitable distribution of external loads in a radial direction. Internal forces and deflections are of the most concern in analytical study of soil-steel structures. But, only precise modeling makes results reliable. The current methods of analyzing these types of structures are principally based on the elasticity theory, in which nonlinear interaction of soil-structure and stage construction of local effects are not taken into account. These are the main concerns of the present research. In this paper, introducing soil-steel structures for bridge construction, the behavior of such structures is investigated under stage construction loads. As a case study, a single span soil-steel bridge has been analyzed and designed, based on the Canadian highway bridge design code (CHBDC) and results are compared with the PLAXIS finite element code. Stage construction consists of 1) the initial phase where a natural trench has been idealized; 2) filling and compaction of soil under foundations; 3) construction of foundation and installation of plates; 4) backfilling and compaction of both sides of the steel structure up to crown level; 5) backfilling the soil cover up to a minimum height of cover for construction loads; 6) backfilling and compaction of soil cover up to the final level for performance loads. In addition to investigate the deformations, axial and shear forces, bending moments and stresses in each of the above construction stages, the stability of the natural trench before construction, the interface between corrugated steel plates and the soil, the assumption of the elastic behavior of backfilling, the effect of unbalanced backfilling, and the effect of the asymmetry of construction loads, have also been studied. The comparison indicates that horizontal and vertical displacements and the compressive stresses obtained from FE stage construction are 64, 43 and 56 percent, respectively, greater under standard limitations. Therefore, a more optimum and economic design, in comparison to elasticity based analysis methods, can be achieved by accounting for stage construction and the nonlinear behavior of soil-steel bridges.Keywords: Soil, steel bridges, static analysis, Finite Element method, truck load