Please enable javascript in your browser.
Page
of
0
مدل سازی عددی رفتار خاک های غیر اشباع (مدل ویلر)در نرم افزار Flac
فروغی، سینا Foroughi, Sina
Cataloging brief
مدل سازی عددی رفتار خاک های غیر اشباع (مدل ویلر)در نرم افزار Flac
پدیدآور اصلی :
فروغی، سینا Foroughi, Sina
ناشر :
صنعتی شریف
سال انتشار :
1393
موضوع ها :
خاک غیراشباع Unsaturated Soil نرم افزار FLAC FLAC Software مدل سازی عددی Numerical...
شماره راهنما :
09-47015
Find in content
sort by
page number
page score
Bookmark
فصل 1: مقدمه
(12)
1-1- پیشگفتار
(13)
1-2- اهداف پژوهش
(14)
1-3- ساختار پژوهش
(15)
فصل 2: مروری بر مطالعات پیشین
(16)
2-1- پیشگفتار
(17)
2-2- مکانیک رفتاری خاک‌های غیراشباع
(17)
2-2-1- مکش در خاک و روش‌های اندازه‌گیری آن
(17)
2-2-2- مفهوم مکش و مشخصات فیزیکی مرتبط با آن
(18)
(2-1) =− .ln()
(18)
(2-2) =+
(18)
(2-3) − = 2
(19)
2-2-3- مکش ساختاری
(19)
شکل (2-1) مویینگی و اثر آن بر مکش ساختاری(Mitchell,1993)
(20)
شکل (2-2) جذب آب توسط ذرات رسی(Jotisankasa, 2005)
(20)
2-2-4- روش‌های اندازه‌گیری مکش ساختاری
(20)
2-2-5- منحنی مشخصه آب-خاک (SWCC)
(21)
شکل (2-3) منحنی مشخصه آب- خاک (Brooks and Corey, 1966)
(22)
2-3- مدل رفتاری خاک‌های غیراشباع
(23)
2-3-1- متغیرهای تنش
(23)
(2-4) ′ = − +( − )
(24)
(2-5) − و −
(24)
(2-6) − و ( − )
(24)
(2-7) − و (− )
(24)
(2-8) − − =
(24)
(2-9) − + − − =
(26)
2-3-2- تغییرات حجم
(26)
2-3-3- مقاومت برشی
(26)
(2-10) = − tan ′ + − tan + ′ = − tan ′ +
(27)
2-4- مدل‌های رفتاری موجود برای خاک‌های غیراشباع
(27)
2-4-1- مدل‌های رفتاری برای خاک‌های غیر متورم شونده
(27)
2-4-1-1- مدل BBM
(27)
(2-11) تنش میانگین خالص − = 1 3 1 +2 3 −
(27)
(2-12) تنش انحرافی = 1 − − 3 − = 1 − 3
(27)
(2-13) مکش = −
(27)
2-4-1-2- مدل Wheeler & Sivakumar (1995)
(28)
2-4-1-3- مدل ارائه‌شده توسط Bolzon (1996)
(29)
(2-14) ′ =− +
(29)
(2-15) =1− tanh(.)
(29)
2-4-2- مدل‌های رفتاری برای خاک‌های متورم شونده
(29)
2-4-2-1- مدل Wheeler (2003)
(29)
(2-16) ∗ = 0 ∗
(30)
(2-17) ∗ = ∗
(30)
(2-18) ∗ = ∗
(30)
(2-19) = ∗ ∗
(30)
(2-20) = (−) 0 ∗ 0 ∗
(30)
(2-21) =0
(30)
(2-22) =− ∗ ∗
(30)
(2-23) =− − ∗ ∗ =− − ∗ ∗
(30)
(2-24) =0
(30)
(2-25) = (−) (1− 1 2 ) 0 ∗ 0 ∗ − 1 ∗ ∗
(30)
(2-26) =− ( − ) (1− 1 2 ) ∗ ∗ − 2 0 ∗ 0 ∗
(31)
2-4-2-2- مدل Santagiuliana (2006)
(31)
(2-27) تجزیه نمو کرنش به دو بخش الاستیک و پلاستیک = +
(31)
(2-28) رابطه بین نمو تنش و کرنش الاستیک ′ =
(31)
(2-29) قانون جریان =
(31)
(2-30) = 1 ′ + 1
(31)
(2-31) = + +
(31)
(2-32) = 1 1+ 2 { ; }
(31)
(2-33) =(1+)( −)
(31)
(2-34) = 1 1+ 2 { ; }
(31)
(2-35) =(1+)( −)
(31)
(2-36) = 0 ′ ( + )
(32)
(2-37) 0 = 1+ 0 0 −
(32)
(2-38) =1+
(32)
(2-39) = 0 1 exp(− 0 )
(32)
2-4-2-3- مدل De'an Sun (2007)
(32)
(2-40) ′ = ′ ( 0 ′ ) 0 − −
(32)
(2-41) 0 = ′ ( ′ ′ ) − 0 −
(32)
(2-42) = 0 + +
(32)
(2-43) ′ = ′ 0 0 + ′
(32)
(2-44) ′ 0 = 0 − − ( 0 ′ ) 0 − −
(32)
(2-45) ′ = ′ ( 0 −) ( − ) 2 ( + ) 2 ′ 0
(32)
(2-46) = ′ (1+) ′
(32)
(2-47) = ( 0 −) 0 (1+) 0
(33)
(2-48) = ( 0 −) (1+) 0 ′ − ′ / ′ 0
(33)
(2-49) =− −
(33)
(2-50) =− −
(33)
2-4-3- مدل‌های ارائه‌شده در سال‌های اخیر:
(33)
2-4-3-1- مدل D. Masin & N. Khalili (2008)
(33)
(2-51) = (ℒ:+ )
(33)
(2-52) = +.
(34)
(2-53) χ= 1 ≥ ( ) <
(34)
(2-54) = + 1− .
(34)
(2-55) ln 1+ = − ∗
(35)
(2-56) =+ ln( )
(35)
(2-57) ∗ = ∗ + ln( )
(35)
(2-58) = ℒ:+ +
(35)
(2-59) =
(35)
(2-60) = () − ()
(35)
(2-61) = ℒ:+ +
(36)
(2-62) = () − () < >
(36)
2-4-3-2- مدل W. Fuentes & Th. Triantafyllidis (2013)
(36)
شکل (2-4) شبیه‌سازی مدل برای شاخه‌های اصلی منحنی SWCC (Brooks & Corey, 1966) (a،(Van (b Genuchten, 1980)
(37)
شکل (2-5) منحنی توسعه‌یافته (Van Genuchten, 1980) برای مواد شکل‌پذیر
(38)
شکل (2-6) شبیه‌سازی تنش ایزوتروپ در شرایط مکش ثابت a) درصدتخلخل b) درجه‌اشباع W.Fuentes & Th.Triantafyllidis, 2013
(38)
(2-63) مدل هیدرولیکی = −(1− ) − (1+)
(38)
(2-64) مدل مکانیکی =( − + 3 <− >), ℎ = +
(39)
2-4-3-3- مدل ارائه‌شده توسط H.Ghasemzadeh & S.A Ghoreishian Amiri (2013)
(39)
شکل (2-7) منحنی SWCC مدل Wheeler (2003)
(39)
شکل (2-8) سطوح تسلیم برای شرایط تنش سه محوری (H.Ghasemzadeh & S.A Ghoreishian Amiri, 2013)
(40)
(2-65) ∗ = 0 ∗ , = 0 , = 0
(40)
(2-66) 0 ∗ = 0 − 1
(40)
(2-67) = − (1− ) 2
(40)
2-4-3-4- مدل ارائه‌شده توسطG.D Nguyen & Y.Gan (2013)
(40)
شکل (2-9) اثر نیروهای کاپیلاره بر مدل a) تماس ذرات در طول سطح فرضی b) ایده‌ال‌سازی تماس بین‌ذره‌ای و نیروی کاپیلاره (G.D Nguyen & Y.Gan, 2013)
(41)
شکل (2-10) اثر درجه اشباع بر نیروهای کاپیلاره (G.D Nguyen & Y.Gan, 2013)
(42)
شکل (2-11) a) اثرات درجه اشباع b)مدول الاستیسیته ، بر مقاومت کششی (G.D Nguyen & Y.Gan, 2013)
(42)
شکل (2-12) a) اثر فشار نرمال b) مدول الاستیسیته ، بر مقاومت برشی (G.D Nguyen & Y.Gan, 2013)
(42)
2-5- مدلسازی عددی خاک‌های غیراشباع:
(43)
2-6- جمعبندی و نتیجهگیری:
(43)
فصل 3: کلیات مدل ویلر
(47)
3-1- پیشگفتار
(48)
3-2- مفاهیم به‌کاررفته در مدل ویلر
(48)
شکل (3-1) هیسترزیس هیدرولیکی در منحنی مشخصه SWCC (Wheeler, 2003)
(49)
شکل (3-2) آب‌های menisci و bulk (Wheeler, 2003)
(49)
شکل (3-3) شماتیک نیروی عمود بر سطوح (N σ)و نیروی مماس بر سطوح(N T) (Wheeler & Karbu, 1995)
(50)
شکل (3-4) جابجایی محل تماس آب-هوا به علت مکش (Wheeler, 2003)
(51)
شکل (3-5) هیسترزیس هیدرولیکی در منحنی مشخصه آب (Wheeler, 2003)
(51)
شکل (3-6) سیکلهای ترشوندگی- خشک شوندگی روی رس (Alonso et al, 1995)
(52)
شکل (3-7) سیکل ترشوندگی-خشک شوندگی بر روی خاک کائولین در شرایط ایزوتروپ (Sharma, 1998) a) حجم مخصوص b) درجه‌اشباع
(53)
شکل (3-8) اثر سیکل ترشوندگی-خشک شوندگی بر رفتارهای بعدی تحت بار ایزوتروپ (Sharma, 1998) a) حجم مخصوص b) درجه‌اشباع
(54)
شکل (3-9) رفتار خاک کائولین در حین بارگذاری ایزوتروپ با مکش ثابت (Sharma, 1998) a) حجم مخصوص b) درجه‌اشباع
(55)
3-3- مهم‌ترین ویژگی‌های مدل Wheeler (2003)
(55)
3-4- متغیرهای تنش و کرنش به‌کاررفته در مدل ویلر
(56)
(3-1) ∗ = −[ + 1− ]
(56)
(3-2) ∗ =− −(1− )
(57)
(3-3) ∗ =.
(57)
3-5- معرفی منحنی‌های تسلیم مدل
(57)
3-5-1- منحنی تسلیم LC
(57)
شکل (3-10) رابطه بین مکش و اضافه نیروی نرمال ناشی از آب menisci (Wheeler, 2003)
(57)
شکل (3-11) منحنی تسلیم LC a) حرکت مستقیم منحنی به علت تسلیم شدن روی منحنی LC b) حرکت همبسته سطح تسلیم به علت ایجاد درجه اشباع برگشت ناپذیر (Wheeler, 2003)
(58)
3-5-2- منحنی‌های SI و SD
(59)
شکل (3-12) نمایش هیسترزیس هیدرولیکی به‌عنوان پروسه‌های الاستیک-پلاستیک (Wheeler, 2003)
(59)
شکل (3-13) حرکت مستقیم منحنی‌های SI و SD به علت تسلیم روی منحنی‌های SI یا SD (Wheeler, 2003)
(60)
شکل (3-14) اثر کرنش حجمی پلاستیک بر منحنی‌های ترشوندگی-خشک شوندگی اولیه (Wheeler, 2003)
(60)
شکل (3-15) حرکت همبسته منحنی‌های SI و SD به علت کرنش‌های حجمی پلاستیک (Wheeler, 2003)
(61)
شکل (3-16) منحنی‌های تسلیم LC ، SD و SI در شرایط ایزوتروپ(q=0) (Wheeler, 2003)
(61)
3-6- معرفی پارامترهای مدل ویلر
(62)
شکل (3-17) مدلی برای منحنی مشخصه آب (Wheeler, 2003)
(62)
3-7- فرمول‌بندی مدل ویلر در فضای تنش‌های همسان
(62)
3-7-1- سطوح تسلیم
(62)
(3-4) ∗ = 0 ∗
(62)
(3-5) ∗ = ∗
(62)
(3-6) ∗ = ∗
(62)
3-7-2- کرنش حجمی
(62)
(3-7) = ∗ ∗
(63)
(3-8) = (−) 0 ∗ 0 ∗
(63)
(3-9) =0
(63)
3-7-3- درجه اشباع
(63)
(3-10) =− ∗ ∗
(63)
(3-11) =− − ∗ ∗ =−( − ) ∗ ∗
(63)
(3-12) =0
(63)
3-7-4- وابستگی منحنی‌های تسلیم
(63)
(3-13) 0 ∗ 0 ∗ = 1 ∗ ∗ = 1 ∗ ∗
(63)
(3-14) ∗ ∗ = ∗ ∗ = 2 0 ∗ 0 ∗
(63)
(3-15) 0 ∗ 0 ∗ = (−) − 1 ( − )
(63)
(3-16) ∗ ∗ = ∗ ∗ =− − + 2 (−)
(64)
(3-17) = (−) (1− 1 2 ) 0 ∗ 0 ∗ − 1 ∗ ∗
(64)
(3-18) =− ( − ) (1− 1 2 ) ∗ ∗ − 2 0 ∗ 0 ∗
(64)
3-8- نحوه یافتن ثابت‌های مدل ویلر در آزمایشگاه
(64)
3-8-1- پارامترهای رفتار مکانیکی برای شرایط تنش همسان( و )
(64)
3-8-2- پارامترهای مرتبط با رفتار هیدرولیکی ( و )
(64)
3-8-3- پارامترهای همبستگی (k2 و k1 )
(65)
3-9- جمعبندی و نتیجهگیری:
(65)
فصل 4: آماده‌سازی روابط ویلر برای نرم‌افزار Flac و صحت سنجی آن
(67)
4-1- پیشگفتار
(68)
4-2- فرمول‌بندی مدل برحسب مسیرهای تنش متعارف
(68)
4-2-1- حالت الاستیک (کرنش حجمی الاستیک و تغییر درجه اشباع برگشت‌پذیر)
(68)
(4-1) ∗ = − + +
(69)
(4-2) ∗ =−
(69)
(4-3) ≠0 و =0
(69)
(4-4) d ∗ = d( p – u ) + sd
(69)
(4-5) ∗ =−
(69)
(4-6) ∗ =+
(69)
(4-7) ∗ =−
(69)
(4-8) = (+ ) (+ )
(69)
(4-9) =− (− )
(70)
(4-10) =( × ∗ ×× )/(1− × 2 ∗ ×× )×
(70)
(4-11) = ∗
(70)
(4-12) =0 و ≠0
(70)
(4-13) d ∗ = ds + sd
(70)
(4-14) ∗ =−
(70)
(4-15) = ( + ) (+ )
(70)
(4-16) =− (− )
(70)
(4-17) = − (− ∗ / 1− × 2 ∗ ×
(71)
4-2-2- حالت پلاستیک
(71)
4-2-2-1- حالت تسلیم روی منحنی LC (بدون رسیدن به سطوح تسلیم SI یا SD)
(71)
(4-18) ∗ =+
(72)
(4-19) ∗ =−
(72)
(4-20) 0 ∗ = ∗ و 0 ∗ = ∗
(72)
(4-21) = − 2 0 ∗ / ∗ − (−) 2 2 0 ∗
(72)
(4-22) ∗ = +
(73)
(4-23) ∗ =−
(73)
(4-24) = ( (−) ∗ −) ∗ − 2 (−) ∗
(73)
4-2-3- حالت تسلیم روی منحنی SI یا SD (بدون رسیدن به سطح تسلیم LC)
(73)
(4-25) ∗ =+
(74)
(4-26) ∗ =−
(74)
(4-27) ∗ = ∗ و ∗ = ∗
(74)
(4-28) = 2 ∗ ∗ − − 2 2 ∗
(75)
(4-29) ∗ = +
(75)
(4-30) ∗ =−
(75)
(4-31) = ( 2 ∗ −) ∗ ( − ) − 2 2 ∗
(76)
4-2-4- حالت تسلیم همزمان روی دو سطح تسلیم
(76)
(3-8) = (−) 0 ∗ 0 ∗
(76)
(4-32) ∗ =+
(76)
(4-33) ∗ =−
(76)
(4-34) = ∗ + (−) ∗ (1− 1 2 )
(77)
(4-35) =− (−) 1 ∗ (1− 1 2 )
(77)
(4-36) = ( − ) 2 ∗ (1− 1 2 )
(77)
(4-37) =− ∗ − ( − ) ∗ (1− 1 2 )
(77)
(4-38) ∗ = 1+ 2 ( +1) −
(77)
(4-39) = 1+ ∗
(77)
(4-40) = − (1+ ) ∗
(77)
(4-41) ∗ = +
(78)
(4-42) ∗ =−
(78)
(4-43) = + − 2 (1+ )
(79)
(4-44) =1+ 2 (1+ ) −
(79)
(4-45) ∗ = ×
(79)
(4-46) = ( ∗ +×) (1+ )
(79)
(4-47) = ∗ + ∗
(79)
(4-48) =. =
(79)
(4-49) = . , = .
(79)
4-3- پیش‌بینی مدل
(79)
شکل (4-1) شرایط مرزی و بارگذاری یک چهارم المان (پژوهش حاضر)
(80)
4-3-2- بارگذاری ایزوتروپ در مکش ثابت
(80)
شکل (4-2) شبیه‌سازی مدل برای تنش ایزوتروپ در مکش ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (Wheeler, 2003)
(81)
شکل (4-3) مقایسه نتایج کد واردشده در نرم‌افزار Flac با مدل ویلر a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (پژوهش حاضر)
(83)
4-3-3- متورم شوندگی و فروریزش در حین تر شوندگی (مکش متغیر، بار ثابت)
(83)
شکل (4-4) شبیه‌سازی مدل در شرایط ترشوندگی با تنش ایزوتروپ ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (Wheeler, 2003)
(84)
شکل (4-5) مقایسه نتایج کد واردشده در نرم‌افزار Flac با مدل ویلر در حالت تنش ایزوتروپ ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (پژوهش حاضر)
(85)
شکل (4-6) مقایسه نتایج کد واردشده در نرم‌افزار Flac در حالت مکش ثابت با تنش‌های ایزوتروپ متفاوت a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (پژوهش حاضر)
(87)
4-3-4- بارگذاری سیکلی (کاهش-افزایش مکش)
(87)
شکل (4-7) شبیه‌سازی مدل برای یک سیکل ترشوندگی و خشک‌شوندگی a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (Wheeler, 2003)
(87)
شکل (4-8) مقایسه نتایج کد وارد شده در نرم‌افزار Flac برای یک سیکل ترشوندگی-خشک شوندگی با مدل ویلر(2003) a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (پژوهش حاضر)
(89)
4-4- حالت چند المانی
(89)
شکل (4-9) مقایسه حالت چند المانی با تک المانی و مدل ویلر در حالت مکش ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (پژوهش حاضر)
(90)
شکل (4-10) مقایسه حالت چند المانی با تک المانی و مدل ویلر در حالت تنش ایزوتروپ ثابت a) فضای تنش اصلاحی b) حجم مخصوص c) درجه‌اشباع (پژوهش حاضر)
(91)
4-5- مقایسه با داده‌های آزمایشگاهی
(92)
شکل (4-11) صحت سنجی مدل وارد شده در نرم‌افزار با شرایط آزمایشگاهی (پژوهش حاضر)
(93)
جدول (4-1) پارامترهای آزمایشگاهی و مورداستفاده برای مدل BBM (Raveevdiraraj, 2009)
(94)
جدول (4-2) ثابت‌های خاک (مطالعه حاضر)
(94)
جدول (4-3) موقعیت اولیه سطوح تسلیم (مطالعه حاضر)
(95)
شکل (4-12) مقایسه نتایج آزمایشگاهی و مدل BBM (Raveevdiraraj, 2009) با مطالعه حاضر، تنش خالص متغیر و مکش ثابت
(95)
شکل (4-13) مقایسه نتایج آزمایشگاهی و مدل BBM (Raveevdiraraj, 2009) با مطالعه حاضر، تنش خالص ثابت و مکش متغیر (سیکل ترشوندگی-خشک‌شوندگی)
(96)
4-6- جمعبندی و نتیجهگیری:
(96)
فصل 5: مدلسازی
(98)
5-1- پیشگفتار
(99)
5-2- خاکریز
(99)
5-2-1- مشخصات خاک و پارامترهای مدل
(99)
جدول (5-1) پارامترهای مدل BBM مورد استفاده در پژوهش Zheng (2013)
(100)
جدول (5-2) پارامترهای مدل ویلر (پژوهش حاضر)
(100)
جدول (5-3) شرایط اولیه سطح تسلیم (پژوهش حاضر)
(100)
5-2-2- شبکه‌بندی مدل ، شرایط مرزی و فرآیند مدلسازی
(101)
شکل (5-1) هندسه مدلZheng, 2013
(101)
شکل (5-2) شبکه‌بندی و شرایط مرزی مدل (Zheng, 2013)
(102)
5-2-3- نتایج به‌دست‌آمده و تحلیل نتایج
(103)
شکل (5-3) بردار جابجایی در خاکریز بعد از اتمام ساخت ناشی از وزن خود Zheng (2013)
(103)
شکل (5-4) بردار جابجایی در خاکریز بعد از تر شوندگی Zheng (2013)
(103)
شکل (5-5) بردار جابجایی در خاکریز بعد از اتمام ساخت (پژوهش حاضر)
(104)
شکل (5-6) کنتور جابجایی‌ها بعد از اتمام خاکریز ناشی از وزن خاکریز (مطالعه حاضر)
(104)
شکل (5-7) بردار جابجایی در خاکریز بعد از تر شوندگی (پژوهش حاضر)
(105)
شکل (5-8) کنتور جابجایی‌ها بعد از تر شوندگی (مطالعه حاضر)
(105)
شکل (5-9) مقایسه نشست محور میانی خاکریز در راستای ارتفاع خاکریز
(107)
شکل (5-10) بردار جابجایی در خاکریز در حالت خشک (پژوهش حاضر)
(107)
شکل (5-11) بردار جابجایی در خاکریز بعد ازتر شوندگی (مطالعه حاضر، یک‌لایه)
(109)
شکل (5-12) کنتور جابجایی‌ها بعد از تر شوندگی (مطالعه حاضر، یک‌لایه)
(109)
5-2-4- افزایش محدوده فونداسیون خاکریز در حالت قبل
(110)
شکل (5-13) بردار جابجایی در خاکریز بعد از اتمام ساخت (مطالعه حاضر، افزایش محدوده فونداسیون)
(111)
شکل (5-14) بردار جابجایی در خاکریز بعد از تر شوندگی (مطالعه حاضر، افزایش محدوده فونداسیون)
(111)
5-3- فونداسیون
(112)
5-3-1- مشخصات مدل و پارامترهای خاک
(112)
شکل (5-15) هندسه مدل و مش بندی Abed & Vermeer (2008)
(113)
جدول (5-4) پارامترهای مدل BBM بر طبق پژوهش Abed & Vermeer (2008)
(113)
جدول (5-5) پارامترهای مدلسازی (مطالعه حاضر)
(113)
جدول (5-6) شرایط اولیه سطوح تسلیم (مطالعه حاضر)
(114)
5-3-2- شبکه‌بندی مدل،شرایط مرزی
(114)
شکل (5-16) شبکه‌بندی فونداسیون (مطالعه حاضر)
(115)
شکل (5-17) شرایط مرزی فونداسیون (مطالعه حاضر)
(115)
شکل (5-18) بردار جابجایی و نشست (مطالعه حاضر)
(116)
شکل (5-19) کنتور جابجایی‌ها و نشست (مطالعه حاضر)
(116)
شکل (5-20) نمودار بار-نشست (مطالعه حاضر)
(117)
شکل (5-21) نمودار بار-نشست (Abed & Vermeer, 2008)
(117)
شکل (5-22) مقایسه پژوهش حاضر با پژوهش Abed & Vermeer (2008)
(118)
5-4- جمعبندی و نتیجهگیری:
(119)
فصل 6: نتیجهگیری و پیشنهادها
(120)
6-1- نتیجهگیری
(121)
6-2- پیشنهادها
(122)
