Use this URL to cite or link to this record in EThOS:
Title: Numerical modelling of coupled behaviour in unsaturated soils
Author: Cabot, Martí Lloret
Awarding Body: University of Strathclyde
Current Institution: University of Strathclyde
Date of Award: 2011
Availability of Full Text:
Access through EThOS:
Access through Institution:
Unsaturated soils are present in a wide range of geotechnical applications such as slopes, highways or earth dams, in addition to many other geo-environmental applications such as underground disposal of radioactive waste or landfills. Despite a significant number of improvements in understanding unsaturated soil behaviour over the last five decades, there are still several aspects to be addressed and to be better understood, including in the areas of constitutive modelling, development of experimental procedures and improvement of numerical techniques. This work concentrates on two of these aspects: constitutive modelling of unsaturated soils accounting for the coupling between mechanical and water retention response; and associated numerical and computational aspects employed to solve engineering problems involving unsaturated soil mechanics. Wheeler et al. (2003) proposed a coupled elasto-plastic constitutive model for the mechanical and water retention behaviour of unsaturated soils and this model constitutes the basis of the work developed within this thesis. For isotropic stress states this coupled model employs as stress state variables mean Bishop’s stress p* and modified suction s*. In the model, plastic volumetric strains occur when the stress state reaches a LC yield curve and plastic changes of degree of saturation Sr occur when the stress state reaches a SD or SI yield curve. Coupled movements of the yield curves represent the influence of plastic changes of degree of saturation on mechanical behaviour and the influence of plastic volumetric strains on water retention behaviour. According to this constitutive model, during many types of loading or wetting under isotropic stress states the soil state will ultimately arrive at the corner between the LC and SD yield curves. Analysis of the model suggests that such states should fall on a unique planar normal compression surface in v:lnp*:lns* space and also on a second unique planar surface in Sr:lnp*:lns* space. The experimental results of Sivakumar (1993) from constant suction isotropic loading of compacted speswhite kaolin are presented in these spaces, and are shown to lie on planar surfaces, as predicted by the model. Presenting experimental normal compression results in this way gives an ideal method for determining the values of the key soil parameters k1, k2, λs and κs in the model. A simple extension of the Wheeler et al. (2003) model to the stress conditions of the triaxial test, by including the role of deviator stress q, has been proposed by Lloret et al. ii (2008b) and Raveendiraraj (2009). According to this extended model, critical states will commonly occur at the intersection of the LC yield surface and the SD yield surface, and such critical states should fall on a second pair of unique planar surfaces in the v:lnp*:lns* and Sr:lnp*:lns* spaces. The experimental critical state data of Sivakumar (1993) have been plotted in these spaces and the results obtained are presented and discussed. The new constitutive model has been extended to general 3D stress conditions and generalised stress-strain relationships required for numerical integration of the model are presented. Furthermore, 3D extended functions to identify the different elasto-plastic mechanisms of the model are proposed. A partial validation of this extended model is also presented using the experimental results of Sivakumar (1993). Two strain-driver algorithms for the integration of the generalised model are presented including first and second order error schemes. The complete formulation of these explicit algorithms is described and the computed results from both integration schemes is compared and discussed. Finally, the generalised fully coupled constitutive model has been implemented within the finite element program CODE_BRIGHT (Olivella et al., 1996) to solve boundary value problems involving unsaturated soil. The performance of the new implemented model is analysed and discussed by considering application to a boundary value problem involving wetting of a cylindrical soil specimen.
Supervisor: Not available Sponsor: Not available
Qualification Name: Thesis (Ph.D.) Qualification Level: Doctoral
EThOS ID:  DOI: Not available