شبیه سازی حالت بحرانی خاک دانه‌ای با دانه های چندگوشه به کمک روش اجزای مجزا

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

Simulation of Critical State Behavior of Granular Soils with Polygonal Particles Using Discrete Element Method (DEM)

In this study, the numerical Discrete Element Method (DEM) was applied for simulating both drained and undrained behavior of granular materials with twodimensional polygonal particles in order to find a critical state line. For undrained behavior simulation, two methods including constant volume and cylinder methods were utilized. In the constant volume method, it was assumed that the volume of the soil remains constant during the loading due to the incompressibility of water. In the cylinder method, however, a pipe was considered among adjacent pores that provide the water transformation between them. In other words, the transmission of water among the voids can be taken into account. An undrained simulation was performed for sandy samples at the confining pressure of 200 kPa by both methods. Simulations showed that the results obtained by the cylinder method have good conformity with those of the constant volume method. A parametric study on the water compressibility was done. As the water becomes less compressible, i.e., stiffer, the stressstrain paths of both methods become closer. Also, the effect of confining pressure on the drained and undrained behavior by constant volume and cylinder methods was investigated. The results of the simulations showed that by increasing the confining pressure, the deviatoric stress and the contraction tendency increase in drained and undrained simulations. To achieve a critical state in the soil samples, the simulations were performed with a large strain level where both deviatoric stress and void ratio become constant. Then the critical state line locus, as well as its parameters, are determined. The results show that the critical state line locus does not depend on the stress path. Furthermore, the simulation method for the undrained condition has very little impact on the critical state line locus.