شبيه سازي حالت بحراني خاك دانه اي با دانه هاي چندگوشه به كمك روش اجزاي مجزا

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

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

In this study, the numerical Discrete Element Method (DEM) is used for simulating drained and undrained behavior of 2D polygonal granular materials in order to find critical state line. For undrained behavior simulation, two methods including constant volume and cylinder methods have been used. In the constant volume method, it is assumed that the volume of the soil remains constant during the loading due to the incommensurability of water. In the cylinder method, however, a pipe is considered among adjacent pores that provides the water transformation between them. Undrained simulation was performed for sandy samples at the confining pressure of 200 kPa by both methods. Simulations show that the results obtained by cylinder method have good conformity with those of the constant volume method. As the water becomes less compressible, i.e., stiffer, the results 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 critical state, the simulations were performed until large strains where the deviatoric stress and void ratio become constant. Then the critical state line locus and also its parameters are determined. The results show that the critical state line locus does not depend on the stress path and the simulation method for undrained condition has very little impact on the critical state line locus.