Consolidation
Tuesday, September 11, 2012 12:43 PM
Steven F. Bartlett, 2010
Consolidation Page 1
e vs log v curves
Tuesday, September 11, 2012 12:43 PM
Cc = compression index
Cs = recompression index
'c = preconsolidation stress
Steven F. Bartlett, 2010
Consolidation Page 2
Settlement Calculations
Tuesday, September 11, 2012 12:43 PM
Consolidation Settlement
Change in void ratio for normally
consolidated clay
Consolidation settlement for
normally consolidated clay
Change in void ratio for
overconsolidated clay below the
preconsolidation stress
Consolidation settlement for
overconsolidated clay with
increase in stress below the
preconsolidation stress
Consolidation settlement
for overconsolidated clay
with increase in stress
below the
preconsolidation stress
Steven F. Bartlett, 2010
Consolidation Page 3
Relationships - Elastic and Consolidation
Tuesday, September 11, 2012 12:43 PM
Constrained modulus
Coefficient of volume
compressibility
Also
Coefficient of volume
compressibility, where
cv is the coefficient of
vertical consolidation
(defined later).
Young's modulus
Note that cv and k are not constant but vary with void ratio or vertical
strain, hence M and E are non linear. cv is discussed later in this
lecture in time rate of consolidation section.
Steven F. Bartlett, 2010
Consolidation Page 4
Developing Elastic Model for Consolidation
Tuesday, September 11, 2012 12:43 PM
Consolidation is a non-linear process that produces a stress-strain
relationship where vertical strain is proportional to the log of the
change in stress.
If would be useful to be able to model this process with a simple
tangent modulus that varies with the vertical stress.
strain
Vary slope of this tangent
line (modulus) with stress
stress
Because the finite difference techniques uses time steps, it is a
simple matter to check the state of stress at any time step and
adjust the tangent modulus to follow this non-linear path.
Thus, we can use an algorithm that allows for the modulus to
change as a function of applied stress to mimic this non-linear
function if we can determine the relationship that expresses the
modulus in terms of the applied stress.
Steven F. Bartlett, 2010
Consolidation Page 5
Developing Elastic Model for Consolidation (cont.)
Tuesday, September 11, 2012 12:43 PM
Step 1 - Define consolidation properties in term of vertical strain,
v, instead of void ratio, e.
1
v
Cc
Log v'
Useful relationships between vertical strain and void ratio
V = H/Ho
V = eo - e / (1+eo)
Cc = Cc / (1+eo)
Step 2 - Find tangent modulus at a point (i.e., derivative)
y = log x
dy/dx = d log (x) / dx
dy/dx = 1 / (x * ln(10))
Let x be v'
Let dy = Cc
dy/dx = Cc/ (v' * ln(10))
Note that the tangent modulus as defined by dy/dx is also
called the constrained modulus, M, for a 1D compression
test.
M = Cc/ (v' * ln(10))
Steven F. Bartlett, 2010
Consolidation Page 6
Developing Elastic Model for Consolidation (cont.)
Tuesday, September 11, 2012 12:43 PM
Step 3 - Express Young's modulus E in terms of Cc and v'
M = (1-)E/[(1+)(1-2)]
ln(10)v'/Cc = (1-)E/[(1+1-2
(1+1-2v'ln(10) = Cc(1-)E
E = (1+1-2v'ln(10) / (Cc(1-))
In the equation above, we have a relationship to define Young's
modulus, E, in terms of the vertical stress, Cc, and The latter
two factors are material properties, which can be determined
from laboratory tests.
Thus we have a method to predict how Young's modulus varies
non linearly as a function of applied vertical stress for given
values of Cc, and .
Steven F. Bartlett, 2010
Consolidation Page 7
FLAC Implementation of Elastic Model
Tuesday, September 11, 2012 12:43 PM
Modeling 1D consolidation test
Axisymmetrical model with height = 2.54 cm, radius = 3.175 cm
Steven F. Bartlett, 2010
Consolidation Page 8
FLAC Implementation of Elastic Model (cont.)
Tuesday, September 11, 2012 12:43 PM
Vertical displacement vectors and displacement after applying 200
kPa, 4 tsf, stress at top of model.
Bulk modulus
Shear modulus
Young's modulus
Vertical strain
Change in moduli with vertical strain
Steven F. Bartlett, 2010
Consolidation Page 9
FLAC Implementation of Elastic Model (cont.)
Tuesday, September 11, 2012 12:43 PM
Vertical stress (Pa)
Vertical strain = 0.12 or 12%
Vertical strain
Check of FLAC results using consolidation theory
v = Cc log [(o + ) / o]
v = 0.25 log [(100 + 200)/100]
v = 0.12 or 12 percent
Steven F. Bartlett, 2010
Consolidation Page 10
FLAC Model
Tuesday, September 11, 2012 12:43 PM
config axisymmetry ats
set large; let grid deform for high strain problems
grid 10,10
gen 0,0 0,0.0254 0.03175,0.0254, 0.03175,0
define inputs ; subroutine for input values and calcs
Cce=0.25
soil_d = 2000 ; soil density
sigma_c = 100e3*(-1) ; preconsolidation stress
sigma_f = 300e3*(-1) ; final stress
P_ratio = 0.45 ; Poisson's ratio
E_ini = (1+P_ratio)*(1-2*P_ratio)*(-1)*sigma_c*2.3026/(Cce*(1-P_ratio)) ; initial Young's modulus
K_ini = E_ini/(3*(1-2*P_ratio)) ; initial bulk modulus
G_ini = E_ini/(2*(1+P_ratio)) ; inital shear modulus
;
new_E = E_ini ; initializes variables used in nonlin
new_K = K_ini ; initializes variables
new_G = G_ini ; initializes variables
end
inputs ; runs subroutine
;
model elastic
prop dens = soil_d bu = K_ini sh = G_ini
;
; boundary conditions
fix x y j 1
fix x i 11
;
; initializes preconsolidation stress in model
define preconsol
loop i (1,izones)
loop j (1,jzones)
syy(i,j) = sigma_c
endloop
endloop
end
preconsol
Steven F. Bartlett, 2010
Consolidation Page 11
FLAC Model (cont.)
Tuesday, September 11, 2012 12:43 PM
define nonlin ; subroutine to change moduli while stepping
sumstress = 0
sum_K = 0
sum_G = 0
sum_E = 0
whilestepping
loop i (1,izones)
loop j (1,jzones)
new_E = (1+P_ratio)*(1-2*P_ratio)*(-1)*syy(i,j)*2.3026/(Cce*(1-P_ratio))
new_K = new_E/(3*(1-2*P_ratio))
new_G = new_E/(2*(1+P_ratio))
bulk_mod(i,j)=new_K
shear_mod(i,j)=new_G;
sumstress = sumstress + syy(i,j)
sum_K = sum_K + new_K
sum_G = sum_G + new_G
sum_E = sum_E + new_E
avgstress = (-1)*sumstress/100 ; average vertical stress in model for 100 zones
v_strain = ydisp(6,11)*(-1)/0.0254
avg_K = sum_K/100 ; average bulk modulus in model
avg_G = sum_G/100 ; average shear modulus in model
avg_E = sum_E/100
endloop
endloop
end
;
; applies new stress at boundary (stress controlled)
set st_damping=local 2.0; required for numerical stabilitygui
apply syy sigma_f from 1,11 to 11,11 ; applies vstress at boundary
;
; applies velocity at boundary (strain controlled)
;apply yvelocity -5.0e-6 xvelocity=0 from 1,11 to 11,11 ;applies constant downward velocity to simulate
a strain-controlled test
;
;
; histories
history 1 avg_K ; creates history of bulk modulus
history 2 avg_G ; creates history of shear modulus
history 3 avg_E ; creates history of elastic modulus
history 4 avgstress ; average stess in model
history 5 unbalanced ; creates history of unbalanced forces
history 6 v_strain; vertical strain
;
solve; use this if stress is applied to top boundary (stress controlled)
;cycle 820; use this if velocity is applied to top boundary (strain controlled)
save 1d consolidation.sav 'last project state'
Steven F. Bartlett, 2010
Consolidation Page 12
Consolidation Settlement Under Strip Footing
Tuesday, September 11, 2012 12:43 PM
Stresses versus depth from elastic theory
Steven F. Bartlett, 2010
Consolidation Page 13
Consolidation Settlement Under Strip Footing (Example 1 -
Normally Consolidated Clay)
Tuesday, September 11, 2012 12:43 PM
Steven F. Bartlett, 2010
Consolidation Page 14
Consolidation Settlement Under Strip Footing (Example 1 cont.)
Tuesday, September 11, 2012 12:43 PM
Steven F. Bartlett, 2010
Consolidation Page 15
Consolidation Settlement Under Strip Footing (Example 2 cont.)
Tuesday, September 11, 2012 12:43 PM
Steven F. Bartlett, 2010
Consolidation Page 16
Consolidation Settlement Under Strip Footing (Example 2 -
Overconsolidated Clay)
Tuesday, September 11, 2012 12:43 PM
Steven F. Bartlett, 2010
Consolidation Page 17
Time Rate of Consolidation
Tuesday, September 11, 2012 12:43 PM
Governing Eq. for 1D Consolidation
We will discuss more about this equation when we cover seepage
and groundwater flow
Coefficient of Consolidation
Steven F. Bartlett, 2010
Consolidation Page 18
Time Rate of Consolidation (cont.)
Tuesday, September 11, 2012 12:43 PM
Steven F. Bartlett, 2010
Consolidation Page 19
Assignment 5
Tuesday, September 11, 2012 12:43 PM
1. Solve Figure 3.26 using the mesh shown in Figure 3.27 in the book for
an infinite strip footing. Solve this using Eq. 3.12 using elastic theory
in Excel (10 points) and solve it using FLAC (10 points). Compare the
results on a composite plot (5 points).
2. Further develop the Excel spread sheet in problem 1 to perform
consolidation settlement calculations for consolidation settlement
underneath the strip footing. Verify your spread sheet using the
examples found in the lecture notes. (10 points NC case, 10 points
OC case)
3. Develop a 2D plane strain FLAC model to perform the settlement
calculations for a strip footing placed on an 3-m clay layer overlying
bedrock. Use the nonlinear elastic method given in these lecture
notes to calculate the non-linear elastic modulus as a function of
stress level as given the lecture notes. Verify the 2D plane strain FLAC
model results with the settlement calculations for a over-
consolidated case from problem using the following (50 points):
(Hint: The total settlement should be about 3 cm)
q= applied stress at the surface from footing = 50 kPa
x = 0 (i.e., stress and center profile located below footing center)
B = 0.15 m (half width of footing)
z = 0.075 m (for spread sheet) (for FLAC use a 20 by 20 grid to
form a 3-m by 3-m box using the gen command: gen 0,-3. 0,0 3,0
3,-3 ratio 1.177 0.8
Zw = 0 (depth to water table)
sat soil = 19.62 kN/m3
p' = 19.62 kPa for the upper 2 m of the soil
The soil below 2 m is normally consolidated
Cc = 0.2
Cr = 0.02
e0 = 1
Steven F. Bartlett, 2016
Consolidation Page 20
