1.14.1 The Terzaghi consolidation problem

Product: ABAQUS/Standard  

This one-dimensional problem has a well-known linear solution (see Terzaghi and Peck, 1948) and, thus, provides a simple verification of the consolidation capability in ABAQUS. The analysis of saturated soils requires solution of coupled stress-diffusion equations, and the formulation used in ABAQUS is described in detail in Analysis of porous media, Section 2.8 of the ABAQUS Theory Manual, and Plasticity for non-metals, Section 4.4 of the ABAQUS Theory Manual. The coupling is approximated by the effective stress principle, which treats the saturated soil as a continuum, assuming that the total stress at each point is the sum of an “effective stress” carried by the soil skeleton and a pore pressure in the fluid permeating the soil. This fluid pore pressure can change with time (if external conditions change, such as the addition of a load to the soil), and the gradient of the pressure through the soil that is not balanced by the weight of fluid between the points in question will cause the fluid to flow: the flow velocity is proportional to the pressure gradient in the fluid according to Darcy's law. A typical case is a consolidation problem. Here the addition of a load (usually an overburden) to a body of soil causes pore pressure to rise initially; then, as the soil skeleton takes up the extra stress, the pore pressures decay as the soil consolidates. The Terzaghi problem is the simplest example of such a process. For illustration purposes, the problem is treated with and without finite-strain effects. The small-strain version is the classical case discussed by Terzaghi and Peck (1948), and the finite-strain version has been analyzed numerically by a number of authors, including Carter et al. (1979).

Problem description

Geometry and models

ABAQUS contains no one-dimensional elements for effective stress calculations. Therefore, we use a two-dimensional plane strain mesh, with one element only in the x-direction. Element type CPE4P is chosen to perform the finite-strain analysis, and element type CPE8P is chosen for the small-strain analysis. We recommend the use of linear elements for applications involving finite strain, impact, or complex contact conditions and second-order elements for problems where stress concentrations must be captured accurately or where geometric features such as curved surfaces must be modeled. In this particular example the linear and second-order elements yield almost identical results.

The soil is assumed to be linear elastic, with a Young's modulus of 689.5 GPa (108 lb/in2) and Poisson's ratio of 0.3. The specific weight of the pore fluid is assumed to be 276.8 × 103 N/m3 (1 lb/in3). The permeability is assumed to vary linearly with the void ratio, with a value of 8.47 × 10–8 m/sec (2.0 × 10–4 in/min) at a void ratio of 1.5 and a value of 8.47 × 10–9 m/sec (2.0 × 10–5 in/min) at a void ratio of 1.0. The void ratio is assumed to be 1.5 initially throughout the sample. ABAQUS uses effective permeability, which is permeability divided by the specific weight of the pore fluid. Therefore, the fluid in this problem is assigned the value 276.8 × 103 N/m3 (1lb/in3) for the specific weight (water, for example, has a specific weight of 9965 N/m3, 0.036 lb/in3) and the permeability is scaled accordingly.

The boundary conditions are as follows. On the bottom and two vertical sides, the normal component of displacement is fixed (0 on the bottom and 0 on the sides), and no flow of pore fluid through the walls is permitted. This latter is the natural boundary condition in the fluid mass conservation equation, so no explicit specifications need to be made (as with zero tractions in the equilibrium equation). On the top surface a uniform downward load (an overburden) is applied suddenly. The magnitude of this load is taken to be 689.5 GPa (108 lb/in2). This large load will cause considerable deformation, thus illustrating the difference between the small- and large-strain solutions. This surface allows perfect drainage so that the excess pore pressure is always zero on this surface.

Time stepping

The problem is run in two steps. The first step is a single increment of a *SOILS, CONSOLIDATION analysis with an arbitrary time step, with no drainage allowed across the top surface (the natural boundary condition in the mass conservation equation governing the pore fluid flow). This establishes the initial solution: uniform pore pressure equal to the load throughout the body, with no stress carried by the soil skeleton (zero effective stress). The actual consolidation is then done with a second *SOILS, CONSOLIDATION step, using automatic time stepping.

The accuracy of the time integration for the second *SOILS, CONSOLIDATION procedure, during which drainage is occurring, is controlled by the UTOL parameter. This parameter specifies the allowable pore pressure change per time step. Even in a linear problem UTOL controls the accuracy of the solution, because the time integration operator is not exact (the backward difference rule is used). In this case UTOL is chosen as 344.8 GPa (5.0 × 107 lb/in2), which is a relatively large value and, so, should only give moderate accuracy: this is considered to be adequate for the purposes of the example.

An important issue in such consolidation problems is the choice of initial time step. As the governing equations are parabolic, the initial solution (immediately after the sudden change in load) is a local, “skin effect,” solution. In this one-dimensional case the form of the initial solution is sketched in Figure 1.14.1–2 for illustration purposes. With a finite element mesh of reasonable size for modeling the solution at a later time (when the changes in pore pressure have diffused into the bulk of the body soil), this initial solution will be modeled poorly. With smaller initial time steps the difficulty becomes more pronounced, as sketched in Figure 1.14.1–2. As in any transient problem, the spatial element size and the time step are related to the extent that time steps smaller than a certain size give no useful information. This coupling of the spatial and temporal approximations is always most obvious at the start of diffusion problems, immediately after prescribed changes in the boundary values. For this particular case the issue has been discussed in detail by Vermeer and Verruijt (1981), who suggest the simple criterion

where is a characteristic element size near the disturbance (that is, near the draining surface in our case), E is the elastic modulus of the soil skeleton, k is the soil permeability, and is the specific weight of the permeating fluid. For our model we choose  254mm (10 in); and we have  689.5  GPa (108 lb/in2), 8.47 × 10–8 m/s (2.0 × 10–4 in/min), 2.768 × 105 N/m3 (1.0lb/in3), which gives .05 s (0.833 × 10–3 min). Based on this calculation, an initial time step of .06 sec (0.001 min) is used. This gives an initial solution with no “overshoot” at all, as expected.

In this case we wish to continue the analysis to steady-state conditions. This is defined by asking ABAQUS to stop when all pore pressure change rates fall below 11.5 KN/m2/s (100 lb/in2/min).

Results and discussion

Input files

References

Figures

Figure 1.14.1–1 Terzaghi consolidation problem definition.

Figure 1.14.1–2 Solutions at very early times.

Figure 1.14.1–3 Pore pressure at various times.

Figure 1.14.1–4 Pore pressure and effective stress at elevation 1.905 m (75 in).

Figure 1.14.1–5 Degree of consolidation versus time factor.

Figure 1.14.1–6 Comparisons of finite- and small-strain solutions to the Terzaghi consolidation problems.