ZSOIL.PC: A UNIFIED APPROACH TO STABILITY, BEARING CAPACITY, CONSOLIDATION, CREEP AND FLOW FOR TWO AND THREE-DIMENSIONAL SIMULATIONS IN GEOTECHNICAL PRACTICE.1/
Th. Zimmermann 1,2, A. Truty 2, A. Urbanski 3,1, S. Commend 2 & K. Podles 3,1
1 - Zace Services Ltd, Lausanne Switzerland
2 - Swiss Federal Institute of Technology Lausanne, Switzerland
3 - Cracow University of Technology, Poland

4. Two-phase medium and flow

The coupled problem of the equilibrium of a soil in presence of transient underground flow conditions is solved using a mixed formulation in which displacement and total water pressure are the main variables. The definition of effective stresses (Eq.4.1) is modified to account for non saturated media:

(4.1)

where the saturation ratio S is a function of the water pressure.

The two-phase coupled equilibrium equations are:

The first equation expresses equilibrium of the two-phase medium in total stresses and bi = gi, with = nS w +(1-n) s and n, the porosity; the second equation expresses transient continuity in the pore fluid and accounts for consolidation and flow. kk is the volumetric strain time derivative in the skeleton, q the averaged (Darcy) fluid velocity and c the specific storage coefficient.

Fig. 4.1 - Boundary and initial conditions

An appropriate set of boundary conditions must be prescribed. These include imposed surface loads in total stress t and imposed displacements (as on 1). In addition, constant pressure values are prescribed where pressure is constant (as on 2), flow across the boundary can be prescribed where necessary ( = 0 on 1) and seepage boundary conditions are imposed on 3 and 4 to account for the fact that the free surface and H2 are functions of time. This means that seepage elements are introduced along these boundaries and external, time-dependent pressures are prescribed. This type of boundary condition will switch automatically to the appropriate one depending on the saturation in the medium.

Appropriate initial conditions must also be prescribed for the fluid in case of transient flow. In the situation described above, these correspond to the hydrostatic pressure in the medium, i.e to a water table following the upper surface of the medium.