Coupled modelling of hydro‐mechanical behaviour of unsaturated compacted expansive soils

The Barcelona basic model cannot predict the mechanical behaviour of unsaturated expansive soils, whereas the Barcelona expansive model (BExM) can only predict the stress–strain behaviour of unsaturated expansive soils without the water‐retention behaviour being incorporated. Moreover, the micro‐parameters and the coupling function between micro‐structural and macro‐structural strains in the BExM are difficult to determine. Experimental data show that the compression curves for non‐expansive soils under constant suctions are shifted towards higher void ratios with increasing suction, whereas the opposite is true for expansive soils. According to the observed water‐retention behaviour of unsaturated expansive soils, the air‐entry value increases with density, and the relationship between the degree of saturation and void ratio is linear at constant suction. According to the above observation, an elastoplastic constitutive model is developed for predicting the hydraulic and mechanical behaviour of unsaturated expansive soils, based on the existing hydro‐mechanical model for non‐expansive unsaturated soil. The model takes into consideration the effect of degree of saturation on the mechanical behaviour and that of void ratio on the water‐retention behaviour. The concept of equivalent void ratio curve is introduced to distinguish the plastic potential curve from the yield curve. The model predictions are compared with the test results of an unsaturated expansive soil, including swelling tests under constant net stress, isotropic compression tests and triaxial shear tests under constant suction. The comparison indicates that the model offers great potential for quantitatively predicting the hydraulic and mechanical behaviour of unsaturated expansive soils. Copyright © 2011 John Wiley & Sons, Ltd.     

饱和度对非饱和土力学性质的影响

现在被广泛公认的由Fredlund提出的非饱和土力学的双参数理论,即净应力和吸力为非饱和土的应力状态量,不能直接考虑饱和度或含水率对非饱和土的应力-应变关系和强度的影响。在非饱和土三轴试验结果表明,即使在净应力和吸力路径相同的条件下,具有不同饱和度试样的应力-应变关系和强度也是不同的。其他条件相同时,试样饱和度越高,其应力比-应变关系曲线越高,强度越大。最新的水力-力学特性耦合的弹塑性本构模型可以定量地表示上述非饱和土的性质     

A hierarchical thermo‐hydro‐plastic constitutive model for unsaturated soils and its numerical implementation

Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities.     

A density‐dependent elastoplastic hydro‐mechanical model for unsaturated compacted soils

This paper presents a three‐dimensional elastoplastic constitutive model for predicting the hydraulic and mechanical behaviour of unsaturated soils. It is based on experimental results obtained from a series of controlled‐suction triaxial tests on unsaturated compacted clay with different initial densities. Hydraulic hysteresis in the water‐retention behaviour is modelled as an elastoplastic process, with the elastic part modelled by a series of scanning curves and the elastoplastic part modelled by the main drying and wetting curves. The effect of void ratio on the water‐retention behaviour is studied using data obtained from controlled‐suction wetting–drying cyclic tests on unsaturated compacted clay with different initial densities. The effect of the degree of saturation on the stress–strain‐strength behaviour and the effect of void ratio on the water‐retention behaviour are considered in the model, as is the effect of suction on the hydraulic and mechanical behaviour. The initial density dependency of the compacted soil behaviour is modelled by experimental relationships between the initial density and the corresponding yield stress and, thereby, between the initial density and the normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure and yield surfaces in the deviatoric stress plane are given by the Matsuoka–Nakai criterion. Model predictions of the stress–strain and water‐retention behaviour are compared with those obtained from triaxial tests with different initial densities under isotropic compression, triaxial compression and triaxial extension, with or without variation in suction. The comparisons indicate that the model accurately predicts the hydraulic and mechanical behaviour of unsaturated compacted soils with different initial densities using the same material constant. Copyright © 2006 John Wiley & Sons, Ltd.     

Temperature effects on the hydraulic behaviour of an unsaturated clay

The influence of temperature on the hydraulic properties of unsaturated clays is of major concern in the design of engineered barriers in underground repositories for high-level radioactive waste disposal. This paper presents an experimental study centred on the investigation of the influence of temperature on soil hydraulic properties related to water retention and permeability. Laboratory tests were conducted on artificially prepared unsaturated fabrics obtained from a natural kaolinitic-illitic clay. Special attention is given to the testing procedures involving controlled suction and temperature oedometer cells and the application of the vapour equilibrium technique at high temperatures. Retention curves at different temperatures show that total suction tends to reduce with increasing temperatures at constant water content. Temperature influence on water permeability is more relevant at low matric suctions corresponding to bulk water preponderance (inter-aggregate zone). Below a degree of saturation of 75% no clear effect is detected. Experimental data show that temperature dependence on permeability at constant degree of saturation and constant void ratio is smaller than what could be expected from the thermal change in water viscosity. This behaviour suggests that phenomena such as porosity redistribution and thermo-chemical interactions, which alter clay fabric and pore fluid, can be relevant.     

Capillary collapse of loose pyroclastic unsaturated sands characterized at grain scale

The reduction in volume for unsaturated soils wetted at constant total stress is indicated as capillary collapse. Several studies conducted on standard laboratory specimens (macro-scale) outlined the role of initial void ratio, confining pressure and matric suction on collapse onset. Conversely, few observations were made at grain scale, although an important influence of soil structure has been supposed since years. This paper investigated the collapse of coarse and fine sands derived from a pyroclastic soil of Southern Italy. The X-ray computed tomography was used to identify the mechanisms acting at grain scale and to measure the local variations of soil structure. The experimental procedure consisted in preparing remoulded unsaturated specimens and reducing the matric suction until the collapse occurred under self-weight. At different stages of the process, the sample was imaged by X-ray tomography. The experimental results provided original insight into: (1) transformation of soil structure during the wetting tests; (2) variation of porosity, water content and degree of saturation for the whole specimen; and (3) local variations of those variables in several representative sub-volumes. It is worth noting that collapse of coarse sand specimen occurred before saturation. This was also emphasized by the presence of macro-voids at collapse.     

Investigation of the hydro-mechanical behaviour of GMZ bentonite pellet mixtures

Bentonite pellet mixtures are considered as one of the candidate sealing materials for deep geological disposals of radioactive waste. One of the particularities of this material is the initial heterogeneous distribution of pellets and porosity within the mixture, leading to complex hydro-mechanical behaviour. In this paper, the hydro-mechanical properties of GMZ bentonite pellet mixtures were investigated in the laboratory by carrying out water retention tests on pellet mixtures under constant-volume condition and single pellets under free swelling condition, as well as a infiltration test on a column specimen of pellet mixture. In the infiltration test, the relative humidity and radial swelling pressure were monitored at five heights, the axial swelling pressure was also recorded. The instantaneous profile method was applied to determine the unsaturated hydraulic conductivities. Results show that, in high suction range (>?10 MPa) the water retention curve of pellet mixture under constant-volume condition was comparable to that of a single pellet under free swelling condition, while in low suction range (<?10 MPa) the latter exhibits a much higher water retention capacity. Due to clogging of large pores, the unsaturated hydraulic conductivity decreases as suction decreases to around 25 MPa. However, with further suction decrease, the hydraulic conductivity increases continuously until the value at saturated state, as in the case of most unsaturated soils. The radial swelling pressure at different heights develops with local sudden increase and decrease, which was attributed to local rearrangement of pellets upon wetting. By contrast, as the axial swelling pressure was measured on the global surface of the specimen, it develops in a more regular fashion.

     

A water retention model accounting for void ratio changes in double porosity clays

This paper presents a constitutive model that predicts the water retention behaviour of compacted clays with evolving bimodal pore size distributions. In line with previous research, the model differentiates between the water present inside the saturated pores of the clay aggregates (the microstructure) and the water present inside the pores between clay aggregates (the macrostructure). A new formulation is then introduced to account for the effect of the macrostructural porosity changes on the retention behaviour of the soil, which results in a consistent evolution of the air-entry value of suction with volumetric deformations. Data from wetting tests on three different active clays (i.e. MX-80 bentonite, FEBEX bentonite, and Boom clay), subjected to distinct mechanical restraints, were used to formulate, calibrate, and validate the proposed model. Results from free swelling tests were also modelled by using both the proposed double porosity model and a published single porosity model, which confirmed the improvement in the predictions of degree of saturation by the present approach. The proposed retention model might be applied, for example, to the simulation of the hydromechanical behaviour of engineered bentonite barriers in underground nuclear waste repositories, where compacted active clays are subjected to changes of both suction and porosity structure under restrained volume conditions.

     

Retention behaviour of natural clayey materials at different temperatures

The water retention capacity of geomaterials, and especially clayey soils, is sensitive to temperature changes as the physical mechanisms of retention, such as capillarity or adsorption, are affected by it. It is therefore a major issue to be able to define temperature-dependent behaviour of materials, especially for geo-energy and geo-environmental applications involving non-isothermal conditions. This paper presents results of experiments conducted on two representative materials: a hard clay (Opalinus clay) and a plastic clay (Boom clay), both of which have been considered as buffer materials for underground radioactive waste disposal, in Switzerland and Belgium, respectively. Two new devices were developed for this purpose to permit the analysis of water retention behaviour at different temperatures. The behaviour of these two materials at ambient (20 °C) and high temperature (80 °C) was observed and described through the evolution of the degree of saturation, the water content and the void ratio with respect to suction. It appears that the retention capability of the clays reduces significantly with an increase in temperature; on the other hand, the change in temperature had less of an effect on the total volume variation.     

Predicting the dependency of a degree of saturation on void ratio and suction using effective stress principle for unsaturated soils

The paper presents an approach to predicting variation of a degree of saturation in unsaturated soils with void ratio and suction. The approach is based on the effective stress principle for unsaturated soils and several underlying assumptions. It focuses on the main drying and wetting processes and does not incorporate the effects of hydraulic hysteresis. It leads to the dependency of water retention curve (WRC) on void ratio, which does not require any material parameters apart from the parameters specifying WRC for the reference void ratio. Its validity is demonstrated by comparing predictions with the experimental data on four different soils taken over from the literature. Good correlation between the measured and predicted behaviour indirectly supports applicability of the effective stress principle for unsaturated soils. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of degree of saturation on mechanical behaviour of unsaturated soils and its elastoplastic simulation

The two stress-state variable approach has been widely used in interpreting unsaturated soil behaviour. However this approach cannot take into account the effect of degree of saturation or water contents on the stress–strain behaviour and strength of unsaturated soils. The triaxial test results presented in this paper show that even if the same path of net stress and suction is followed, the stress–strain relation and strength are different due to different degrees of saturation. When other conditions are the same, the higher the degree of saturation for the soil sample is, the higher the stress ratio corresponding to a given axial strain will be. This effect can be modeled by using an elasto-plastic constitutive model coupling hydraulic and mechanical behaviour of unsaturated soils. Comparisons between the predicted and measured results are presented, which demonstrate that the model can quantitatively simulate the influence of the degree of saturation on stress–strain behaviour and strength of unsaturated soils.     

Influence of Infiltration of Sodium Chloride Solutions on SWCC of Compacted Bentonite–Sand Specimens

The degree of saturation of compacted bentonite buffer in deep geological repositories is subject to alterations from infiltration of groundwater and heat emanated from the waste canisters. The matric suction (ψ)–degree of saturation (S _r ) relations of unsaturated clays is represented by soil–water characteristics curves (SWCC) that are influenced by soil structure, initial compaction condition and stress history. Infiltration of groundwater besides increasing the degree of saturation can also alter the pore water chemistry; the associated changes in cation hydration and diffuse double layer thickness could impact the micro-structure and matric suction values. This study examines the influence of infiltrating sodium chloride solutions (1,000–5,000 mg/L) on the transient ψ–S _r relations of compacted bentonite–sand specimens. Analysis of the ψ–S _r plots, and X-ray diffraction measurements indicated that infiltration of sodium chloride solutions has progressively less influence on the micro-structure and SWCC relations of bentonite–sand specimens compacted to increasingly higher dry densities. The micro-structure and SWCC relations of specimens compacted to 1.5 Mg/m³ were most affected, specimens compacted to 1.75 Mg/m³ were less affected, while specimens compacted to 2 Mg/m³ remained unaffected upon infiltration with sodium chloride solutions.     

Experimental Set-up for Determining Soil Water Retention Curves for Granular Soils During Drying

Abstract: Soil water retention curves (SWRCs) provide an important means of describing the response of unsaturated soils during drying / wetting processes in terms of variations of degree of saturation, water content or void ratio with suction. A key consideration in generating these curves is how to measure the suction. Frequently the filter paper technique is adopted, especially when high suctions are developed, e.g. with plastic clays. As each measurement takes at least a week with this technique, it can take months or years to generate a full SWRC in drying and wetting. Developments in laboratory tensiometers now allow matrix suctions up to about 1.5 MPa to be measured. With such a device it is possible to develop SWRCs for granular soils such as silts and clays in hours or days by continuous measurement. This paper describes an experimental set-up that was developed to measure changes in volume, water content and matrix suction during drying of three granular soils. Limitations of the apparatus and usefulness of the curves are discussed.     

Shear strength of a compacted scaly clay in variable saturation conditions

Scaly clays are stiff and highly fissured clays often used as construction materials. This paper presents the results of triaxial compression tests carried out on saturated and unsaturated samples of a compacted scaly clay. Complementary investigation on the microstructural features and their evolution with the amount of water stored into the material are also presented in order to shed light on the evolution of the micro- and macroporosity with suction. The water retention behaviour of the compacted scaly clay is also addressed. The results from the controlled suction triaxial tests are used to discuss the applicability of a single-shear strength criterion to compacted double-structured clays when the effective stress concept for unsaturated soils is used. The choice of the degree of saturation to be included in the effective stress definition for obtaining a satisfactory representation of the shear strength is addressed. It is shown that the best results are obtained when the macropore degree of saturation is considered along with its evolution during the applied stress path.     

Experimental investigation on the effect of the mode of deposition on the pore-access size distribution of sand

It is well known that the specimen preparation method and the resulting sand fabric significantly affect sand behaviour and sand liquefaction resistance. In this study, a simple experimental method namely tensiometric method is presented for evaluation of the effect of the samples reconstituting methods on microstructure and the pore-access size distribution of sands. Specimens were prepared using dry funnel pluviation (DFP) and wet deposition (WD) at two densities Dr = 40% (loose state) and Dr = 80% (dense state). The water-retention curve of a sand specimen is evaluated using a technique similar to the hanging column method by application of a step-by-step variation of suction and measure of the corresponding changes in water content. The pore-access size distribution is then evaluated from the water-retention curve. In this experimental study the effects of the mode of deposition and the uniformity coefficient on the pore-access size distribution are investigated. It was found that the specimen generated by the WD method contained a higher percentage of pore size, contrary to the specimen prepared by the DFP which contained a limited percentage of pore size, causing a great influence on the mechanical behaviour of the specimens.     

Role of degree of saturation in denitrification of unsaturated sand specimens

Nitrate contamination of groundwater arises from anthropogenic activities, such as, fertilizer and animal manure applications and infiltration of wastewater/leachates. During migration of wastewater and leachates, the vadose zone (zone residing above the groundwater table), is considered to facilitate microbial denitrification. Particle voids in vadose zone are deficient in dissolved oxygen as the voids are partially filled by water and the remainder by air. Discontinuities in liquid phase would also restrict oxygen diffusion and therefore facilitate denitrification in the vadose/unsaturated soil zone. The degree of saturation of soil specimen (S _r) quantifies the relative volume of voids filled with air and water. Unsaturated specimens have S _r values ranging between 0 and 100 %. Earlier studies from naturally occurring nitrate losses in groundwater aquifers in Mulbagal town, Kolar District, Karnataka, showed that the sub-surface soils composed of residually derived sandy soil; hence, natural sand was chosen in the laboratory denitrification experiments. With a view to understand the role of vadose zone in denitrification process, experiments are performed with unsaturated sand specimens (S _r = 73–90 %) whose pore water was spiked with nitrate and ethanol solutions. Experimental results revealed 73 % S _r specimen facilitates nitrate reduction to 45 mg/L in relatively short durations of 5.5–7.5 h using the available natural organic matter (0.41 % on mass basis of sand); consequently, ethanol addition did not impact rate of denitrification. However, at higher S _r values of 81 and 90 %, extraneous ethanol addition (C/N = 0.5–3) was needed to accelerate the denitrification rates.     

Microstructural deformation mechanisms of unsaturated granular soils

A discrete model for unsaturated granular soils has been developed. Three discrete entities have been defined: particles, water menisci and pores. Local interaction forces and water transfer mechanisms have been integrated into a model through the appropriate equilibrium and balance equations. The results of several numerical tests using this model have been described and discussed. Simulations include wetting and drying under load tests, the application of suction cycles and the effect of a deviatoric stress ratio on wetting‐induced collapse. The model reacts just as true granular soil samples behave in laboratory tests. The model provides a new insight into the internal mechanisms leading to large‐scale features of behaviour such as wetting‐induced collapse or the increase in soil strength provided by suction. The paper also stresses that matric suction changes acting on a granular structure are capable of explaining most of the macroscopic features of stress–strain behaviour. Copyright © 2002 John Wiley & Sons, Ltd.     

非饱和土的水力和力学特性及其弹塑性描述

简单回顾了非饱和土本构模型研究的发展历程,总结了近几年非饱和土弹塑性本构模型最新研究成果,重点介绍了能统一模拟非饱和土水力性状和力学性状耦合的弹塑性本构模型。通过对建立模型过程中的几个核心问题讨论,较详细地说明该类模型的结构、性能以及相关问题。非饱和土水力性状的滞回性用假定存在饱和度弹性区间的弹塑性过程来模拟;该类耦合模型不仅考虑了吸力对非饱和土水力性状和力学性状的影响,还考虑了饱和度对应力-应变关系和强度的影响以及土体变形对土-水特征曲线的影响。用同一套模型参数,耦合模型可统一预测在吸力控制或含水率控制下沿各种应力路径下非饱和土的水力-力学特性,并简单介绍了膨胀性非饱和土的弹塑性本构模型以及耦合模型在有限元数值计算中的应用。     

Coupled consolidation in unsaturated soils: An alternative approach to deriving the Governing Equations

The equations governing coupled consolidation in unsaturated soils are known to contain additional parameters when compared to the equations for saturated soils. Nonetheless, the variation of these parameters with suction or degree of saturation is not generally agreed upon. The paper introduces a novel approach to deriving general equations for each of these parameters and their variation, and explains that, for consistency with the constitutive and soil-water retention curve models adopted, these general equations need to be transformed into case-specific expressions. Finally, a conceptual model is presented highlighting how the behaviour of unsaturated soil reflects aspects of its water content.