Popis: |
Zemin emmesi ile su içeriği arasındaki ilişkiyi gösteren zemin-su karakteristik eğrisi (ZSKE) doygun olmayan zeminlerin drenaj özelliklerini yansıtır. Yakın geçmişte ZSKE ölçümü için geliştirilmiş olan çeşitli deneysel yöntemler çoğunlukla zaman alıcı, pahalı, ve kullanıcının becerisinden fazlaca etkilenen prosedürlerdir. Bu problemin bir çözümü ampirik veya yarı-ampirik hesaplar yoluyla temel zemin özelliklerinden dönüştürmek veya mevcut veritabanlarındaki eğrilerle benzeştirmek yoluyla ZSKE'ne bağlamaktır. Diğer bir yaygın çözüm kısmi ZSKE veri noktasına regresyon ile bir eğri oturtmaktır. İsabet ve maliyet açısından daha akılcı bir yaklaşım zemin ortamını gözenek ölçeğinde bilgisayarla modellemek olabilir. Bu çalışma toprak ortamının drenajının MATLAB platformunda hazırlanmıs bir bilgisayar yazılımıyla modellenmesi yoluyla ZSKE'nin elde edilmesini hedeflemiştir. Önerilen yaklaşımda öncelikle dane boyu dağılımı ve boşluk oranı esas alınarak zemin yapısı oluşturulur. Sonra boşluk suyu hacmi boşluk gövdesi ve sıvı köprüleri olarak ikiye ayrılır. Her üç birbirine komşu parçacığın arasında kalan boşluğun en küçük kesiti boşluk boğazı olarak varsayılır, ve her boşluk boğazının hava giriş basıncı hesaplanabilir. Bu amaç için Young-Laplace denklemi sonlu farklar yakınsaması ve Newton-Raphson (Jacobi) sayısal metodlarının birlikte kullanılması yoluyla çözülür. Bu simülasyonda suya uygulanan emme kademeli olarak arttırılarak drenaj adım adım izlenebilir. Sonunda boşluk gövdelerinde ve sıvı köprülerinde kalan su içeriği her emme kademesi için hesaplanarak ZSKE çizilir. Önerilen yöntemin kullanımı plastik olmayan zeminlerle sınırlıdır, kil danelerini modellemek için uygun değildir. Geliştirilen yöntem deneysel sonuçlarile neredeyse mükemmel tutarlılığa sahip ZSKE'lerini nisbeten hızlıca çıkarabilmektedir. Ayrıca geliştirilen yöntemin, çeşitli gerçek toprak numuneleri için literatürdeki ampirik hesap yöntemlerinin bulabildiğinden çok daha isabetli şekilde ZSKE'leri hesaplayabildiği gösterilmiştir. The soil water characteristic curve (SWCC) reflects fundamental drainage properties of partially saturated soils by revealing the relation between soil suction and water content. Although during the past decades, various experimental methods have been proposed for obtaining the SWCC, these approaches are generally time consuming, expensive and highly dependent on operator skills. One solution to this has been sought through various empirical or physico-empirical formulations that link the SWCC to common soil properties, through transformation calculations or within available SWCC databases. Another solution is to fit a curve onto partial SWCC data via regression analysis. A more reasonable approach both in terms of accuracy and cost could be to simulate the pore-scale drainage of soil medium with application of computer modeling techniques. This study is aims to model the drainage of soil medium through developing a computer program in MATLAB for simulating the drying SWCC. In the proposed approach, initially, a solid phase of soil structure is generated based on the particle size distribution as well as void ratio (e) of a given soil sample. Next, the volume of pore water is divided into two main portions: (i) bulk water and, (ii) liquid bridges. Assuming the inter-particle region of three adjacent particles with minimum hydraulic radius as a pore throat, the air entry pressure of each pore throat can be computed. To this end, the Young-Laplace equation is solved through a combination of the finite difference approximation and Newton-Raphson (Jacobian) numerical methods. In this simulation the amount of applied suction is gradually increased and resulting drainage scheme of bulk pores is visualized step by step. Eventually the final residual water content of bulk pores and liquid bridges after each suction increment are determined and plotted as SWCC. The usefulness of the developed computer program is verified by comparison of the obtained SWCC to those attained based on experimental and empirical techniques. The application of proposed method is limited to non-plastic soils, and is not suitable for modeling clay particles. The developed method is capable of providing SWCC, with an almost perfect level of conformity to the experimental results, in a timely manner. Moreover, it is shown that the SWCCs obtained by the proposed algorithm are much more accurate than results of the estimation methods found in the literature for a variety of real soils. 167 |