An Elasto-Viscoplastic Model and Multiphase Coupled FE Analysis for Unsaturated Soil

Rate sensitivity is an important characteristic of geomaterials for both saturated and unsaturated soils. However, many constitutive models for unsaturated soil have been constructed within the framework of the rate independent theory. The present study addresses an elasto-viscoplastic constitutive model which considers the effect of suction for unsaturated clayey soil and a soil-water-air three-phase coupled analysis using the elasto-viscoplastic model. The proposed constitutive model adopts the average skeleton stress for the effective stress from the viewpoint of the mixture theory. Hence, it has become possible to construct a model for unsaturated soil starting with a model for saturated soil by substituting the average skeleton stress for the effective stress and introducing the suction effect into the constitutive model. Furthermore, the collapse behavior, which is brought about by a decrease in suction, is described by the shrinkage of the overconsolidation boundary surface, the static yield surface, and the viscoplastic potential surface. A numerical analysis for multiphase materials is conducted within the framework of a continuum mechanics approach through the use of the theory of porous media. The theory is a generalization of Biot's two-phase mixture theory for saturated soil. A soil-water-air three-phase coupled finite element method is developed in the present study using the governing equations for multiphase soil based on the non-linear finite deformation theory. The average skeleton stress is defined as the difference between the total stress and the average pressure of the two fluids and is used in the proposed elasto-viscoplastic constitutive model. A van Genuchten (1980) type of equation is employed as the constitutive equation between the liquid saturation and the suction pressure. Numerical simulations of unexhausted-undrained compression with different strain rates are conducted under plane strain conditions, and the applicability of the proposed method is evaluated with respect to strain localization and the effect of suction.