Unit 5 Part 4-Modified
Unit 5 Part 4-Modified
Effective Stress
A 1-D flow condition is the one where the velocity vectors are all
parallel and of equal in magnitude .
In other words, the water always moves parallel to some
axis and through a constant c/s area . 1-D flow
Steady downward flow occurs when
water is pumped from an underground
aquifer.
Steady upward flow occurs as a result of
artesian pressure when a less permeable layer
is underlain by a permeable layer which is
connected through the ground to a water
source providing pressures higher than local
Flow through hydrostatic pressures.
a confined aquifer
Different types of fluid flow in soils 2-D flow
It occurs when all the velocity vectors are confined to a single plane but vary in direction
and magnitude within that plane.
Long excavation
E.g.
Flow towards a water well
Multi-dimensional fluid flow in soils
As of now, we have considered one -dimensional flow in soils, where all fluid is flowing
in the same direction.
In most cases, however, fluid in different regions will be flowing in different directions. →
Multi-dimensional flow
The concept of flow net is based on Laplace’s equation of continuity, which describes the steady
flow condition for a given point in the soil mass.
Seepage
Equation of continuity
Details of small
z q+ z qy
element at A → dqz
qx + dqx
qx
x
qz + dqz
H1 h H2 y
qz
Flow at point A
Datum
A
Seepage
Equation of continuity
Flows entering the soil prism in x , y, and z directions can be given from Darcy’ s law
as:
h
q x = k x i x A x = k x (dydz)
x
h h = hydraulic
q =k i A =k (dxdz)
y head at point A
y
y y y y
h
q z = k z i z A z = k z (dxdy)
z
Equation of continuity
The respective flows leaving the soil prism in x, y, and z
directions can be given from Darcy’s law as:
h 2 h
q x + dq x = k x (ix + dix ) Ax = k x + 2 dx (dydz)
x x
h + 2 h dy (dxdz)
q y + dq y = k y (i y + di y ) Ay = k y
y y
2
h + 2 h dz (dydx)
q z + dq z = k z (iz + diz ) Az = k z
z z
2
Seepage
Equation of continuity
qin = qout
Seepage
Equation of continuity
h h h
2 2 2
+ + =0 kx = ky = kz = k
dx2 dy2 dz2 k 0
Simplified Laplace ’s equation of continuity
2h
2
= 0 A = area of
dz C/ S of the L h z
specimen
Change in i z
gradient/unit distance
in z-direction
i.e. vx = vy = 0
Ix =0; iz = constant.
Possible methods for solving the Laplace Equation
The flow net solution is a graphical method of solving the 2-D Laplace’s equation:
h h
2 2
h2
+ =0 =0
x2 z2 x2
2-D flow 1-D flow
h
= v x = −k (A)
x x
h
= v z = −k (B) z
z z
Seepage - Flow net solution
By differentiating and substituting in Laplace equation of continuity:
We get
2 2
+ =0
x 2
z 2
h ( x , z ) = − C − ( x , z )
1 -- (C)
k
If h (x, z) represents a constant h 1, equation (C) represents a
curve in (x, z) plane.
Seepage - Flow net solution
= vx = vz
It gives: x z
d = dx + dz
x z dx dz vx
d = vxdx + vzdz vz
h
= v x = −k 1
z x 2
h
− = v z = −k z
x z Inverse of potential function
Combining equations and substituting in 2 2
Laplace equation: → →
So, (x, z) also satisfies the + =0
Laplace equation x 2
z 2
Seepage - Flow net solution
From the definition of partial differentiation and combining equations:
It gives:
− = vz = vx
x z
Total differential of the (x, z) is: vz
dx
d = dx + dz dz
x z
vx
d = − vzdx + vxdz
For a given flow line, if is a constant,
d z = v z
d = 0 dx v
x
Seepage - Flow net solution
It is apparent that the flow lines and equipotential lines intersect each
other at right angles (or orthogonal).
dz vx
= −
d x vz
d z = v z
dx v
x
Pore pressure in steady state seepage conditions
u = i w ( s )
Seepage - Flow net solution
Number of
Potential drops
Nd = 6
Seepage pressure
Seepage- Flow net solution : Computation of discharge q
Aspect ratio: a / b
hL / N
i= d
Hydraulic gradient i = h/ b b
Equipotential drops between two flow lines h =hL /Nd
From Darcy’s law, flow in each channel is:
hL / N a q = Total discharge
q = k d
per unit width
b
a
N a
q = khL f
N d b b
q ? and p w at B
Nf = 4
Nd = 12
h =
2/12
Drawing Flow Nets
• To calculate the flow and pore pressures in the ground a flow net must be drawn.
• The flow net must be comprised of a family of orthogonal lines (preferably defining
a square mesh) that also satisfy the boundary conditions.
Common boundary conditions
a. Submerged soil boundary - Equipotential
Water
H-z
H uw
h = +z
w
z now
uw = w ( H − z) (7)
Datum so
( H − z) w
h = +z= H
w
Common boundary conditions
b. Impermeable soil boundary - Flow Line
Permeable Soil
vn=0
Flow Line
vt
Impermeable Material
Common boundary conditions
c. Line of constant pore pressure - eg. phreatic surface
uw
Head is given by h = +z
w
u w
and thus h = + z
w
Step 1 - Draw the system to scale (no-flow boundaries are indicated by gray zones),
Step 2 - Draw equipotential lines to coincide with head boundaries (black lines),
Step 3 - Draw flow lines to coincide with no-flow boundaries (blue arrow following the no-flow boundary
of the dam
Drawing a Flow Net for Flow Beneath an Impermeable Dam
Step 4: Draw flow lines along paths where you envision groundwater flowing (blue arrows), ensuring they are
perpendicular to equipotential lines on the boundaries. Do not be concerned about getting them right initially.
Drawing flow nets is a trial-and-error process. As experience with flow nets grows, intuition improves and it
becomes easier to place flow lines in nearly the right position on the first attempt
Drawing a Flow Net for Flow Beneath an Impermeable
Dam
Step 5: Draw equipotential lines between lines drawn at constant head boundaries (black lines), ensuring they are
perpendicular to no-flow boundaries, perpendicular to flow lines and attempting to form curvilinear squares.
Drawing a flow net by hand is a trial-and-error process because the equipotential lines and flow lines are adjusted
until curvilinear squares are formed. It is useful to sketch round shapes within and touching the boundaries of the
space formed by the equipotential lines and flow lines. If the shapes are not circular, as in this first attempt to draw
the flow net, then the lines should be adjusted.
Drawing a Flow Net for Flow Beneath an Impermeable Dam
Step 6: Creating shapes with a constant aspect ratio is a requirement when drawing a flow net. The best way to
achieve that is by drawing curvilinear squares. Sketching a circle within the shapes can help discern whether the
shapes are curvilinear squares. A slight misfit is not important. The misfits need to be large enough such that it is
necessary to add or delete flow or equipotential lines in order to obtain the near-curvilinear squares because
achieving the proper ratio of the number of flow lines and equipotential lines is key to drawing a valid flow net.
The number of flow lines is same in this figure and previous Figure , but the number of equipotential lines differ
indicating the redrawing was necessary to obtain a flow net that can be used to calculate flow through the system.
Value of head on equipotentials
H
h = (9)
Number of potential drops
Phreatic line
15 m
Datum
h = 15m h=0
h = 3m
h = 12m h = 9m h = 6m
Calculation of flow
Phreatic line
15 m
h = 15m h=0
h = 3m
h =12m h = 9m h = 6m
h = 15m h=0
5m
P
h = 3m
h = 12m h = 9m h = 6m
5m uw
Pore pressure from Ph = +z (11a)
w
At P, using dam base
as datum u w = [12 − ( −5)] w (11b)
Example problem
1
Flow net around the horizontal drain
2-D problem Datum
9
9
3m 1
7 3 7
6m
5 5
k = 5 x 10-7 m/s 3
Sand
3 2
6
Long horizontal drain
Flow net around the horizontal drain
Example problem
a) Find the discharge through the drain in m3 /day per metre
length of drain
➢ Consider datum at free water surface
➢ This makes total head at top of the ground is 3 m
5 1 0 − 7 ( 3 ) 2 ( 2 4 3 6 0 0 )
6
9
0.1426m3 / day Per m run
b) Find the PWP at X, 1.5 m into the soil, directly above the drain
Potential drop between two EQP lines = 3/9
Total head at X = (5.6/9) x 3 = 1.867 m
Elevation head at X = - 1.5 m
If a flow occurs parallel to the bedding planes the flow lines are also parallel to
these planes and the hydraulic gradient has at every point of every bed the same
value i, which is independent of the coefficient of permeability of the beds.
khorizontal =
Flow nets in Anisotropic Material
On the other hand, if flow occurs at right angles to the bedding planes every
water particle passes in succession through every one of the layers. Since the
flow is continuous the discharge velocity v must be the same in every layer,
whereas the coefficient of permeability of the layers is different. Hence the
hydraulic gradient must also be different in the different layers.
kvertical =
Flow nets in Anisotropic Material
In developing procedure for plotting flow nets, we assumed that
the permeable layer is isotropic, i.e.
khorizontal = kvertical = k
Consider the case of constructing flow nets for seepage through soils that
show anisotropy with respect to permeability
2h 2h Is not a Laplacian
kx 2 + kz 2 = 0
x z equation, since this equation
applies for a soil mass in
which stratification occurs.
◆Also results in distorted flow net (not necessary that
flow lines and Eqp lines are orthogonal)
Flow nets in Anisotropic Material
k x 2h 2h
By rewriting: + = 0 -- A
k z x z
2 2
c= kz
kx
2h 2h
Equation (B) is of the same form as + = 0
x z
2 2
X= kz
kx
x
Flow nets in Anisotropic Material
Q Direction of flow Q
kz kt
D
kx
L kz
Flow through natural L kx
anisotropic section
Assumption: Anisotropy merely Flow through transformed
causes linear distortion of the (isotropic) section
flow net
Flow nets in Anisotropic Material
Starting with the natural section, the flow quantity per unit length is:
Q = (L x 1) kz i ---- (C)
Within the transformed section, the same quantity of flow is calculated
as follows:
Q= (L kz
kx
x 1 ) kt i ---- (D)
=
kz kvertical
2. Determine
kx khorizontal
Nf
q = k x k z (h)
Nd
Example problem:
1
c= kz
kx
=
4
kt = k x k z = 4x10−8 m / s
Piping Failure In Hydraulic Structures
When the water flows in the soil in horizontal direction,
it exerts seepage force on the particles in the horizontal
direction and when it flows vertically, the force is
exerted in vertical direction.
For the safety of the structure the exit gradient must be very low from the critical
hydraulic gradient. This piping phenomenon leads to the whole structure failure and it is
so serious that for safety against piping a factor of safety of at least 6 is recommended.
If the dam that holds the water is an earth dam, that is it is made up of earth material such as rock and soil,
this backward erosion piping may also occur in the body of structure.
2. Heave piping failure
Heave meaning is to lift or raise something with great force.
This kind of failure may also occur on the downstream side of a
hydraulic structure when water is coming out vertically upwards.
It occurs when the upward seepage force acting on the particles over
an area exceeds the downward force because of submerged weight
of the soil above that area.
This condition destabilizes the soil and if this uplift force acts
opposite the force of gravity due to the weight of the dam and if it is
high enough, it may destabilize the whole hydraulic structure.
Safety of hydraulic structures against piping
In Terms of forces
To prevent backward erosion piping or heave piping failure the upward seepage force exerted by
water on soil particles should be less than the downward force due to weight of soil or structure
or any other arrangement.
So factor of safety against piping in the soil can be written as weight of everything above which
is downward force divided by upward seepage force applied by water.
If value of downward force that is weight of overlying material is twice the value of upward
force we receive the factor of safety as 2. This means the structure is twice safe. Downward force
is greater than uplift force so structure is safe and is twice strong.
Safety of hydraulic structures against piping
In Terms of hydraulic gradient
Exit gradient is the hydraulic gradient at the downstream side of the flow line
where percolating water leaves the soil mass and emerges as free water.
If Δh is the head loss across the last flow field in the flow net occurring over distance of L,
then
ie= Δh /L
To avoid piping ie<icr
FOS= icr/ie
FOS against piping: 6-7 for fine sands, 4-5 for coarse sand.
Methods for preventing the piping failure
The loaded filter increases the downward force to counter the upward
seepage force at the point where water emerges out from the soil. Hence the
loaded filter increases the factor of safety against piping and the factor of
safety is given by weight of soil plus weight of filter divided by upward
seepage force.
Filters are so designed that they are more permeable than the soil they protect
and yet their particles are not very big or not so coarse that the fine particles
of the soil they protect can move through the voids in the filter.
Filter Design
When seepage water flows from a soil with relatively
fine grains into a coarser material, there is danger that the
fine soil particles may wash away into the coarse
material.
Over a period of time, this process may clog the void
spaces in the coarser material. Hence, the grain-size
distribution of the coarse material should be properly
manipulated to avoid this situation.
A properly designed coarser material is called a filter.
Figure 8.15 shows the steady-state seepage condition in
an earth dam which has a toe filter. For proper selection
of the filter material, two conditions should be kept in
mind:
Condition 1: The size of the voids in the filter material should be small enough to hold the larger particles of
the protected material in place.
Condition 2: The filter material should have a high hydraulic conductivity to prevent buildup of large seepage
forces and hydrostatic pressures in the filters.
Filter Design
https://www.soilmanagementindia.com/soil/seepage-analysis/how-
to-construct-a-flow-net-for-seepage-analysis-laplaces-equation-
and-methods/16514
Filter design
ic
FOS =
According to Harza (1935)
3
6
5
a) From flow net; Nf = 3; Nd = 8; hL = 3 m;
Nf
q = k(hL ) = 0.836 m3 / day per metre
Nd
EH = -9 m
PH = 10.5 m ; → PWP = 105 kN/ m 2
c) Head loss per EQP drop h = 3/8 = 0.375 m
Max. exit hydraulic gradient = 0.375/2.6 =
0.144