Loading...
Numerical Modeling of Unsaturated Soils Behavior (Wheeler Model) in Flac Software
Foroughi, Sina | 2015
2011
Viewed
- Type of Document: M.Sc. Thesis
- Language: Farsi
- Document No: 47015 (09)
- University: Sharif University of Technology
- Department: Civil Engineering
- Advisor(s): Ahmadi, Mohammad Mehdi
- Abstract:
- The mechanical behavior of unsaturated soils is much more complex than saturated soils and the effective stress principle cannot be initially applied to them. The soils found in most parts of Tehran, Iran, are unsaturated and this underlines the importance of having access to softwares capable of analyzing such soils; despite the fact that most available softwares are unable to do so. This thesis aims at adding one of the common unsaturated soil models to a geotechnical software and use that to model and analyze earth structures such as earth dams, foundations, embankments, etc. The Wheeler model (2003) for unsaturated soils is used in this study. This model takes advantage of two independent stress variables to model a soil’s behavior. The most prominent characteristic of this model is coupling the mechanical and hydraulic behaviors of soils in addition to considering the effect of saturation on the mechanical behavior. The geotechnical software used in this study is Flac 2D which uses the Finite Difference Method (FDM) to solve the problems. In this study, the Wheeler model is first added to the software as the constitutive model of unsaturated soils, and then the results are compared to lab data and the implemented code is thus verified. Eventually, the software is used to model an embankment and an earth foundation and the results are analyzed. The comparison of the lab results and the basic BBM model results in the literature with those of this study shows that considering the coupling of the mechanical and hydraulic behavior results in a more accurate modeling and analysis. In addition, it was observed that considering the unsaturated behavior in soils results in a lower amount of settlement. In the end, the effect of climate condition on soil settlement, the most important of which is soil collapse during the wetting process, is investigated
- Keywords:
- Unsaturated Soil ; FLAC Software ; Numerical Modeling ; Wheeler Model
- محتواي کتاب
- view
- فصل 1: مقدمه
- فصل 2: مروری بر مطالعات پیشین
- 2-1- پیشگفتار
- 2-2- مکانیک رفتاری خاکهای غیراشباع
- 2-3- مدل رفتاری خاکهای غیراشباع
- 2-4- مدلهای رفتاری موجود برای خاکهای غیراشباع
- 2-4-1- مدلهای رفتاری برای خاکهای غیر متورم شونده
- 2-4-2- مدلهای رفتاری برای خاکهای متورم شونده
- 2-4-3- مدلهای ارائهشده در سالهای اخیر:
- 2-4-3-1- مدل D. Masin & N. Khalili (2008)
- 2-4-3-2- مدل W. Fuentes & Th. Triantafyllidis (2013)
- 2-4-3-3- مدل ارائهشده توسط H.Ghasemzadeh & S.A Ghoreishian Amiri (2013)
- 2-4-3-4- مدل ارائهشده توسطG.D Nguyen & Y.Gan (2013)
- شکل (2-9) اثر نیروهای کاپیلاره بر مدل a) تماس ذرات در طول سطح فرضی b) ایدهالسازی تماس بینذرهای و نیروی کاپیلاره (G.D Nguyen & Y.Gan, 2013)
- شکل (2-10) اثر درجه اشباع بر نیروهای کاپیلاره (G.D Nguyen & Y.Gan, 2013)
- شکل (2-11) a) اثرات درجه اشباع b)مدول الاستیسیته ، بر مقاومت کششی (G.D Nguyen & Y.Gan, 2013)
- شکل (2-12) a) اثر فشار نرمال b) مدول الاستیسیته ، بر مقاومت برشی (G.D Nguyen & Y.Gan, 2013)
- 2-5- مدلسازی عددی خاکهای غیراشباع:
- 2-6- جمعبندی و نتیجهگیری:
- فصل 3: کلیات مدل ویلر
- 3-1- پیشگفتار
- 3-2- مفاهیم بهکاررفته در مدل ویلر
- شکل (3-1) هیسترزیس هیدرولیکی در منحنی مشخصه SWCC (Wheeler, 2003)
- شکل (3-2) آبهای menisci و bulk (Wheeler, 2003)
- شکل (3-3) شماتیک نیروی عمود بر سطوح (N σ)و نیروی مماس بر سطوح(N T) (Wheeler & Karbu, 1995)
- شکل (3-4) جابجایی محل تماس آب-هوا به علت مکش (Wheeler, 2003)
- شکل (3-5) هیسترزیس هیدرولیکی در منحنی مشخصه آب (Wheeler, 2003)
- شکل (3-6) سیکلهای ترشوندگی- خشک شوندگی روی رس (Alonso et al, 1995)
- شکل (3-7) سیکل ترشوندگی-خشک شوندگی بر روی خاک کائولین در شرایط ایزوتروپ (Sharma, 1998) a) حجم مخصوص b) درجهاشباع
- شکل (3-8) اثر سیکل ترشوندگی-خشک شوندگی بر رفتارهای بعدی تحت بار ایزوتروپ (Sharma, 1998) a) حجم مخصوص b) درجهاشباع
- شکل (3-9) رفتار خاک کائولین در حین بارگذاری ایزوتروپ با مکش ثابت (Sharma, 1998) a) حجم مخصوص b) درجهاشباع
- 3-3- مهمترین ویژگیهای مدل Wheeler (2003)
- 3-4- متغیرهای تنش و کرنش بهکاررفته در مدل ویلر
- 3-5- معرفی منحنیهای تسلیم مدل
- 3-5-1- منحنی تسلیم LC
- 3-5-2- منحنیهای SI و SD
- شکل (3-12) نمایش هیسترزیس هیدرولیکی بهعنوان پروسههای الاستیک-پلاستیک (Wheeler, 2003)
- شکل (3-13) حرکت مستقیم منحنیهای SI و SD به علت تسلیم روی منحنیهای SI یا SD (Wheeler, 2003)
- شکل (3-14) اثر کرنش حجمی پلاستیک بر منحنیهای ترشوندگی-خشک شوندگی اولیه (Wheeler, 2003)
- شکل (3-15) حرکت همبسته منحنیهای SI و SD به علت کرنشهای حجمی پلاستیک (Wheeler, 2003)
- شکل (3-16) منحنیهای تسلیم LC ، SD و SI در شرایط ایزوتروپ(q=0) (Wheeler, 2003)
- 3-6- معرفی پارامترهای مدل ویلر
- 3-7- فرمولبندی مدل ویلر در فضای تنشهای همسان
- 3-8- نحوه یافتن ثابتهای مدل ویلر در آزمایشگاه
- 3-9- جمعبندی و نتیجهگیری:
- فصل 4: آمادهسازی روابط ویلر برای نرمافزار Flac و صحت سنجی آن
- 4-1- پیشگفتار
- 4-2- فرمولبندی مدل برحسب مسیرهای تنش متعارف
- 4-2-1- حالت الاستیک (کرنش حجمی الاستیک و تغییر درجه اشباع برگشتپذیر)
- 4-2-2- حالت پلاستیک
- 4-2-3- حالت تسلیم روی منحنی SI یا SD (بدون رسیدن به سطح تسلیم LC)
- 4-2-4- حالت تسلیم همزمان روی دو سطح تسلیم
- (3-8) = (−) 0 ∗ 0 ∗
- (4-32) ∗ =+
- (4-33) ∗ =−
- (4-34) = ∗ + (−) ∗ (1− 1 2 )
- (4-35) =− (−) 1 ∗ (1− 1 2 )
- (4-36) = ( − ) 2 ∗ (1− 1 2 )
- (4-37) =− ∗ − ( − ) ∗ (1− 1 2 )
- (4-38) ∗ = 1+ 2 ( +1) −
- (4-39) = 1+ ∗
- (4-40) = − (1+ ) ∗
- (4-41) ∗ = +
- (4-42) ∗ =−
- (4-43) = + − 2 (1+ )
- (4-44) =1+ 2 (1+ ) −
- (4-45) ∗ = ×
- (4-46) = ( ∗ +×) (1+ )
- (4-47) = ∗ + ∗
- (4-48) =. =
- (4-49) = . , = .
- 4-3- پیشبینی مدل
- شکل (4-1) شرایط مرزی و بارگذاری یک چهارم المان (پژوهش حاضر)
- 4-3-2- بارگذاری ایزوتروپ در مکش ثابت
- 4-3-3- متورم شوندگی و فروریزش در حین تر شوندگی (مکش متغیر، بار ثابت)
- شکل (4-4) شبیهسازی مدل در شرایط ترشوندگی با تنش ایزوتروپ ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجهاشباع (Wheeler, 2003)
- شکل (4-5) مقایسه نتایج کد واردشده در نرمافزار Flac با مدل ویلر در حالت تنش ایزوتروپ ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجهاشباع (پژوهش حاضر)
- شکل (4-6) مقایسه نتایج کد واردشده در نرمافزار Flac در حالت مکش ثابت با تنشهای ایزوتروپ متفاوت a) فضای تنش اصلاحی b) حجم مخصوص c) درجهاشباع (پژوهش حاضر)
- 4-3-4- بارگذاری سیکلی (کاهش-افزایش مکش)
- 4-4- حالت چند المانی
- 4-5- مقایسه با دادههای آزمایشگاهی
- شکل (4-11) صحت سنجی مدل وارد شده در نرمافزار با شرایط آزمایشگاهی (پژوهش حاضر)
- شکل (4-12) مقایسه نتایج آزمایشگاهی و مدل BBM (Raveevdiraraj, 2009) با مطالعه حاضر، تنش خالص متغیر و مکش ثابت
- شکل (4-13) مقایسه نتایج آزمایشگاهی و مدل BBM (Raveevdiraraj, 2009) با مطالعه حاضر، تنش خالص ثابت و مکش متغیر (سیکل ترشوندگی-خشکشوندگی)
- 4-6- جمعبندی و نتیجهگیری:
- فصل 5: مدلسازی
- 5-1- پیشگفتار
- 5-2- خاکریز
- 5-2-1- مشخصات خاک و پارامترهای مدل
- 5-2-2- شبکهبندی مدل ، شرایط مرزی و فرآیند مدلسازی
- 5-2-3- نتایج بهدستآمده و تحلیل نتایج
- شکل (5-3) بردار جابجایی در خاکریز بعد از اتمام ساخت ناشی از وزن خود Zheng (2013)
- شکل (5-4) بردار جابجایی در خاکریز بعد از تر شوندگی Zheng (2013)
- شکل (5-5) بردار جابجایی در خاکریز بعد از اتمام ساخت (پژوهش حاضر)
- شکل (5-6) کنتور جابجاییها بعد از اتمام خاکریز ناشی از وزن خاکریز (مطالعه حاضر)
- شکل (5-7) بردار جابجایی در خاکریز بعد از تر شوندگی (پژوهش حاضر)
- شکل (5-8) کنتور جابجاییها بعد از تر شوندگی (مطالعه حاضر)
- شکل (5-9) مقایسه نشست محور میانی خاکریز در راستای ارتفاع خاکریز
- شکل (5-10) بردار جابجایی در خاکریز در حالت خشک (پژوهش حاضر)
- شکل (5-11) بردار جابجایی در خاکریز بعد ازتر شوندگی (مطالعه حاضر، یکلایه)
- شکل (5-12) کنتور جابجاییها بعد از تر شوندگی (مطالعه حاضر، یکلایه)
- 5-2-4- افزایش محدوده فونداسیون خاکریز در حالت قبل
- 5-3- فونداسیون
- 5-3-1- مشخصات مدل و پارامترهای خاک
- 5-3-2- شبکهبندی مدل،شرایط مرزی
- شکل (5-16) شبکهبندی فونداسیون (مطالعه حاضر)
- شکل (5-17) شرایط مرزی فونداسیون (مطالعه حاضر)
- شکل (5-18) بردار جابجایی و نشست (مطالعه حاضر)
- شکل (5-19) کنتور جابجاییها و نشست (مطالعه حاضر)
- شکل (5-20) نمودار بار-نشست (مطالعه حاضر)
- شکل (5-21) نمودار بار-نشست (Abed & Vermeer, 2008)
- شکل (5-22) مقایسه پژوهش حاضر با پژوهش Abed & Vermeer (2008)
- 5-4- جمعبندی و نتیجهگیری:
- فصل 6: نتیجهگیری و پیشنهادها