Constitutive and numerical modelling of unsaturated soils (original) (raw)
Related papers
Géotechnique, 54(4): 293-295, 2004
The authors should be commended for presenting such an innovative, elaborate elasto-plastic framework for predicting constitutive unsaturated soil response under isotropic stress states. One of the key aspects of the work is the authors' deep knowledge of the physical interactions that lie behind the mechanical behaviour of unsaturated soils. The proposed constitutive model and their associated conceptual framework provide an important improved understanding of soil behaviour. The present review is intended to offer some comments that further highlight the relevance of the work and its potential for further refinement and tuning.
International Journal for Numerical and Analytical Methods in Geomechanics, 2003
This paper presents a complete finite-element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress-strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time-stepping scheme which automatically adjusts the time-step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non-linearity caused by suction-dependent plastic yielding, suction-dependent degree of saturation , and saturation-dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress-strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction.
Géotechnique, 53(1): 123-135, 2003
The paper presents an elasto-plastic model for unsaturated soils that takes explicitly into account the mechanisms with which suction affects mechanical behaviour as well as their dependence on degree of saturation. The proposed model is formulated in terms of two constitutive variables directly related to these suction mechanisms: the average skeleton stress, which includes the average fluid pressure acting on the soil pores, and an additional scalar constitutive variable, 'csi', related to the magnitude of the bonding effect exerted by meniscus water at the inter-particle contacts. The formulation of the model in terms of variables closely related to specific behaviour mechanisms leads to a remarkable unification of experimental results of tests carried out with different suctions. The analysis of experimental isotropic compression data strongly suggests that the quotient between the void ratio, e, of an unsaturated soil and the void ratio es, corresponding to the saturated state at the same average soil skeleton stress, is a unique function of the bonding effect due to water menisci at the inter-particle contacts. The same result is obtained when examining critical states at different suctions. Based on these observations, an elastoplastic constitutive model is developed using a single yield surface the size of which is controlled by volumetric hardening. In spite of this simplicity, it is shown that the model reproduces correctly many important features of unsaturated soil behaviour. It is especially remarkable that, although only one yield surface is used in the formulation of the model, the irreversible behaviour in wetting–drying cycles is well captured. Because of the behaviour normalisation achieved by the model, the resulting constitutive law is economical in terms of the number of tests required for parameter determination.
Finite element formulation and algorithms for unsaturated soils. Part I: Theory
2003
This paper presents a complete finite-element treatment for unsaturated soil problems. A new formulation of general constitutive equations for unsaturated soils is first presented. In the incremental stress-strain equations, the suction or the pore water pressure is treated as a strain variable instead of a stress variable. The global governing equations are derived in terms of displacement and pore water pressure. The discretized governing equations are then solved using an adaptive time-stepping scheme which automatically adjusts the time-step size so that the integration error in the displacements and pore pressures lies close to a specified tolerance. The non-linearity caused by suction-dependent plastic yielding, suction-dependent degree of saturation , and saturation-dependent permeability is treated in a similar way to the elastoplasticity. An explicit stress integration scheme is used to solve the constitutive stress-strain equations at the Gauss point level. The elastoplastic stiffness matrix in the Euler solution is evaluated using the suction as well as the stresses and hardening parameters at the start of the subincrement, while the elastoplastic matrix in the modified Euler solution is evaluated using the suction at the end of the subincrement. In addition, when applying subincrementation, the same rate is applied to all strain components including the suction.
Theoretical Evaluation of the Mechanical Behavior of Unsaturated Soils
Geotechnical and Geological Engineering, 2011
An evaluation method for the mechanical behavior of unsaturated soils is studied in this paper. Although the mechanical behavior of unsaturated soils is complicated, a simple modeling is preferable in practice. This is because the soil properties are not homogeneous and ground data is limited when structures are being designed. In addition, in order to evaluate the reliability of the design, the physical meanings of the parameters applied in the prediction model should be clear. Firstly, the authors study the relationship between compaction curves and compression indexes in the unsaturated state that is used in the proposed constitutive model. Based on the constitutive model, the stress paths for constant volume shear tests are formulated under a constant void ratio condition and the stress paths for undrained shear tests are calculated under a constant water content condition. In the case of unsaturated specimens, the volume of these specimens changes with the shear deformation and the stress paths depend on the initial degree of saturation. The results of the calculation qualitatively describe the test results by considering the changes in effective confining pressure in the undrained condition and the water retention curves.
Normally, suction-controlled triaxial tests are used to characterize soil behavior in constitutive modeling of unsaturated soils. However, this type of tests requires sophisticated equipment and is time-consuming. This has been one of the major obstacles to the implementation and dissemination of unsaturated soil mechanics beyond the research context. In contrast to suction-controlled triaxial tests, the suction-controlled oedometer test requires simpler equipment and a shorter testing period. Oedometer tests represent the at-rest earth pressure (K0) condition, which is an important stress state in any simulation. The major disadvantage of the oedometer test is that its lateral stress is controlled by the condition of zero lateral strain and remains unknown during the testing process. At present, no well-established, simple, and objective methods are available that take advantage of oedometer test results for constitutive modeling purposes. This paper derives an explicit formulation of the at-rest coefficient for unsaturated soils and develops an optimization approach for simple and objective identification of material parameters in elasto-plastic models for unsaturated soils using the results from suction-controlled oedometer tests. This is achieved by combining a modified state surface approach (MSSA), recently proposed to model the elasto-plastic behavior of unsaturated soils, with the quasi-Newton method to simultaneously calibrate all parameters governing virgin behavior in elasto-plastic models. The Barcelona Basic Model (BBM) is used to demonstrate the application of the proposed explicit formulation and calibration method. Results predicted using obtained parameters are compared with laboratory test results for the same stress paths in order to evaluate the simplicity and objectivity of the proposed method.
On the choice of stress-strain variables for unsaturated soils and its effect on plastic flow
Geomechanics for Energy and the Environment, 2018
The net stress plus suction and the average skeleton stress plus modified suction are two alternative sets of energetically consistent stress variables for modelling the hydro-mechanical behaviour of unsaturated soils. When used in conjunction with their work-conjugate strains, both sets of stress variables correctly calculate the first-order term of the hydro-mechanical work input into a soil element subjected to infinitesimal changes of deformation and water content. They therefore also correctly calculate the increment of internal energy along a given stress-strain path, that is the integral of the first-order term of the infinitesimal work input. This paper shows, however, that the above two sets of stress variables lead to different values of the second-order term of the hydro-mechanical work input. They are therefore no longer equivalent with respect to other aspects of material behaviour governed by the second-order work such as the flow rule of elasto-plastic models. The flow rule assumes the normality between plastic strains and equipotential surfaces defined in the conjugate stress-strain space. This normality is however lost when an elasto-plastic model originally formulated in terms of net stress plus suction is recast in terms of average skeleton stress plus modified suction (or vice versa) by means of standard mapping relationships between stress variables. To restore normality in both stress spaces, it is necessary to impose specific forms of elastic and plastic behaviour.
A new insight into modelling the behaviour of unsaturated soils
International Journal for Numerical and Analytical Methods in Geomechanics, 2012
Understanding the response of partially saturated soils under different loads is important for the design and construction of economical and safe geotechnical engineering structures. This paper presents a coupled elastoplastic constitutive model for predicting the hydraulic and stress-strain-strength behaviour of unsaturated soils. The model proposed is built according to the following principle. A constitutive relation is given for each phase (solid, liquid and gas) and coupling relations between each phases are also derived. In the present case, we assume that each phase is not miscible and that pressure in voids not filled by water remains more or less constant, which is reasonable for most geotechnical problems. Therefore, the model is written in a classical manner with a non associated elastoplastic model for the granular skeleton behaviour; an incompressible liquid phase; a water retention description; and an assumption of the existence of an effective stress concept defined by Bishop.
A constitutive model for partially saturated soils
Géotechnique, 1990
The Paper presents a constitutive model for describing the stress-strain behaviour of partially saturated soils. The model is formulated within the framework of hardening plasticity using two iodependent sets of stress variables: the excess of total stress over air pressure and the suction. The mode1 is able to represent, in a consistent and unified manner, many of the fundamental features of the behaviour of partially saturated soils which had been treated separately by previously proposed models. On reaching saturation, the mode1 becomes a conventional critical state model. Because experimental evidence is still limited, the model has been kept as simple as possible in order to provide a basic framework from which extensions are possible. Tbe mode1 is intended for partially saturated soils which are slightly or moderately expansive. After formulating the model for isotropic and biaxial stress states, typical predictions are described and compared, in a qualitative way, with characteristic trends of the behaviour of partially saturated soils. Afterwards, the results of a number of suction-controlled laboratory tests on compacted kaolin and a sandy clay are used to evaluate the ability of the model to reproduce, quantitatively, observed behaviour. The agreement between observed and computed results is considered satisfactory and confirms the possibilities of reproducing the most important features of partially saturated soil behaviour using a simple general framework. KEYWORDS: compaction; coustitutive relatioas; partial saturation; plasticity. L'article presente un modele constitutif pour d&crire le comportement contrainte/deformation des sols partiellement sat&s.