Modelling the behaviour of unsaturated soil using an elastoplastic constitutive model

Abstract

A numerical model of the stress–strain behaviour of unsaturated soil is presented. An extended critical state elastoplastic constitutive model based on net stress and suction is solved via the finite-element method. Numerical techniques are proposed which allow for the development of an elastoviscoplastic solution algorithm within the context of a transient finite-element analysis. The simultaneous solution of the coupled partial differential equations provides variations of net stress, suction, specific volume and degree of saturation in an unsaturated soil. Several examples are then analysed which illustrate typical features of the behaviour of unsaturated soil samples. Expected patterns of behaviour are reproduced and good correlation is obtained between numerical results and an alternative solution of the theoretical model. Comparisons of numerical predictions with results obtained from two suction-controlled experiments on compacted kaolin and sandy clay are then presented. Good correlation is again achieved, yielding confidence in the validity of the numerical algorithm. Finally, the solution of a boundary value problem is attempted, namely the prediction of seasonal ground movement. Comparison of numerical results with field observed data shows good correlation is achieved. The ability of the new approach to solve boundary problems of geotechnical significance is thus illustrated. Nous présentons une maquette numérique du comportement contrainte–déformation d'un sol non saturé. Une maquette constitutive élastoplastique d'état critique étendue basée sur une contrainte nette et une succion est résolue par la méthode d'éléments finis. Nous proposons des techniques numériques qui tiennent compte du développement d'un algorithme de solution élastoviscoplastique dans le contexte d'une analyse d'éléments finis transitoire. La solution simultanée des équations différentielles épardelles couplées donne des variations de la contrainte nette, de la succion, du volume éspécifique et du degré de saturation dans un sol non saturé. Nous analysons alors plusieurs éééexemples qui illustrent les caractéristiques type du comportement d'échantillons de sol non saturé. Nous reproduisons les modes attendus de comportement et nous obtenons une bonne corrélation entre les résultats numériques et la solution alternative du modèle théorique. Nous comparons les prédictions numériques avec les résultats de deux expériences à succion sur le kaolin compacté et Pargile sableuse. Encore une foie, une bonne corrélation est obtenue, confirmant la validité de Palgorithme numérique. Enfin, nous essayons la solution d'un probléme à valeur limite, nommément la prédiction des mouvements saisonniers du sol. La comparaison des résultats numériques avec des données observées sur le terrain montre qu'on obtient une bonne corrélation. Nous illustrons ainsi Pefficacité de la nouvelle méthode à résoudre les problèmes limites à signification géotechnique.