شبیه‌سازی عددی رفتار زهکشی‌نشده خاک دانه‌ای با ذرات چندگوشه به کمک روش اجزای مجزا

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

Numerical simulation of undrained behavior of granular materials with polygonal particles by discrete element method (DEM)

In the present study, discrete element method is used to simulate undrained twodimensional behavior of soil samples containing polygonal particles. Two methods of constant volume and cylinder were used for numerical simulation of undrained behavior and the results from these two methods were compared. In the cylinder method, it is assumed that there are channels among adjacent pores and the possibility of water exchange among the centers of the pores is provided. In this method, the diameter of the cylinder represents the permeability of the soil. In the constant volume method, it is assumed that the sample volume stays constant during loading. The simulations were done by using both methods. The samples were initially subjected to confining pressure of 200 kPa and after consolidation, they were loaded under deviatoric stress. The results showed that the output obtained by the cylinder method are in good agreement with the constant volume method, and by increasing the stiffness of the water, the results of both methods are closer to each other. Then, the undrained behavior was investigated by using cylinder method in terms of shear strength and pore water pressure in a sand sample at the confining pressures of 200, 400, 800 and 1600 kPa. The trend obtained from the simulation results is in good agreement with laboratory results such that with bigger confining pressure, the shear strength and the positive pore pressure are bigger too, but by increasing confining pressure, the amount of pore pressure decreases in the specimen. The pore pressure distribution contour at the strains of 10% and 30% was presented by means of cylinder method and the effect of cylinder diameter changes on the pore pressure distribution contour was investigated. By investigating the pore pressure distribution contour at the strains of 10% and 30%, it was shown that pore pressure in the sample center has the lowest value. In addition, it was seen that by increasing the cylinder diameter, the variation of pore pressure among the voids becomes smaller, which indicates the uniformity of the pore pressure in permeable soils. By reducing the diameter of the cylinder, the variation of pore pressure becomes more through the sample.