GT NOTES BCV502 notes (1)
GT NOTES BCV502 notes (1)
Geotechnical Engineering-I
Course Content
Module-I
Introduction: Origin of soils, formation of soils, clay mineralogy and soil structure, basic
terminology and their relations, index properties of soils. Soil classification: Particle size
distribution, use of particle size distribution curve, Particle size classification, textural
classification, HRB classification, Unified classification system, Indian standard soil
classification system, Field identification of soils. capillary tension, capillary siphoning. Stress
conditions in soil: Total stress, pore pressure and effective stress
Module-II
Seepage analysis: Seepage pressure, quick condition, Laplace equation for two –dimensional
flow, flow net, properties and methods of construction of flow net, application of flow net,
seepage through anisotropic soil and non-homogenous soil, seepage through earth dam.Inverted
filter and design of inverted filter.
Module-III
Module-IV
Shear strength of soils: Mohr’s stress circle, theory of failure for soils, determination of shear
strength (direct shear test, tri-axial compression test, unconfined compression test, van shear
test), shear characteristics of cohesionless soils and cohesive soils.
Module-V
1
Lecture Note on GTE-I
Reference Books:
Reference Books:
1. Geotechnical Engineering, T.N. Ramamurthy & T.G. Sitharam, S. Chand & Co.
2. Soil Mechanics, T.W. Lambe& Whiteman, Wiley Eastern Ltd, Nw Delhi.
Course Outcomes:
2 Solve any practical problems related to soil stresses estimation, permeability and seepage
including flow net diagram.
3 Formulate practical problems related to consolidation settlement and time rate of settlement.
5 Use stabilization techniques for soft and expansive soil by using various methods.
PO1 PO2 PO3 PO4 PO5 PO6 PO7 PO8 PO9 PO10 PO11 PO12
CO1 2 1 1 1 1
CO2 3 2 1 1 1
CO3 2 1 3
CO4 3 2 1 2 1 1
CO5 3 2 1 2 1 1 1 1 2 2 2 1
2
Lecture Note on GTE-I
3
Lecture Note on GTE-I
1.1 Physical weathering reduces the size of the parent rock material, without any change in
the original composition of the parent rock. Physical or mechanical processes taking place on the
earth's surface include the actions of water, frost, temperature changes, wind and ice. They cause
disintegration and the products are mainly coarse soils.
The main processes involved are exfoliation, unloading, erosion, freezing, and thawing. The
principal cause is climatic change. In exfoliation, the outer shell separates from the main rock.
Heavy rain and wind cause erosion of the rock surface. Adverse temperature changes produce
fragments due to different thermal coefficients of rock minerals. The effect is more for freeze-
thaw cycles.
1.2 Chemical weathering
Chemical weather not only breaks up the material into smaller particles but alters the nature of
the original parent rock itself. The main processes responsible are hydration, oxidation, and
carbonation. New compounds are formed due to the chemical alterations.
Rain water that comes in contact with the rock surface reacts to form hydrated oxides, carbonates
and sulphates. If there is a volume increase, the disintegration continues. Due to leaching, water-
soluble materials are washed away and rocks lose their cementing properties.
Chemical weathering occurs in wet and warm conditions and consists of degradation by
decomposition and/or alteration. The results of chemical weathering are generally fine soils with
altered mineral grains.
The effects of weathering and transportation mainly determine the basic nature of the soil (size,
shape, composition and distribution of the particles).
The environment into which deposition takes place, and the subsequent geological events that
take place there, determine the state of the soil (density, moisture content) and the structure or
fabric of the soil (bedding, stratification, occurrence of joints or fissures)
Transportation agencies can be combinations of gravity, flowing water or air, and moving ice. In
water or air, the grains become sub-rounded or rounded, and the grain sizes get sorted so as to
form poorly-graded deposits. In moving ice, grinding and crushing occur, size distribution
becomes wider forming well-graded deposits.
In running water, soil can be transported in the form of suspended particles, or by rolling and
sliding along the bottom. Coarser particles settle when a decrease in velocity occurs, whereas
finer particles are deposited further downstream. In still water, horizontal layers of successive
sediments are formed, which may change with time, even seasonally or daily.
Wind can erode, transport and deposit fine-grained soils. Wind-blown soil is generally
uniformly-graded. A glacier moves slowly but scours the bedrock surface over which it passes.
Gravity transports materials along slopes without causing much alteration.
1.3 BASICS OF SOIL CLASSIFICATION
Soils as they are found in different regions can be classified into two broad categories:
(1) Residual soils
(2) Transported soils
4
Lecture Note on GTE-I
(d) Melting of a glacier causes the deposition of all the materials scoured by it leading to
formation of glacial deposits.
(e) Soil particles carried by wind and subsequently deposited are known as aeolian deposits.
Soil is not a coherent solid material like steel and concrete, but is a particulate material. Soils, as they
exist in nature, consist of solid particles (mineral grains, rock fragments) with water and air in
the voids between the particles. The water and air contents are readily changed by changes in
ambient conditions and location.
As the relative proportions of the three phases vary in any soil deposit, it is useful to consider a soil
model which will represent these phases distinctly and properly quantify the amount of each
phase. A schematic diagram of the three-phase system is shown in terms of weight and volume
symbols respectively for soil solids, water, and air. The weight of air can be neglected.
5
Lecture Note on GTE-I
𝑒 = ---------------------------
𝑉𝑣
and is expressed as a decimal.
Eq.1
𝑉𝑠
6
Lecture Note on GTE-I
2. Porosity (n) is the ratio of the volume of voids to the total volume of soil (V ), and
𝑛 = -------------------------------
𝑉𝑣
is expressed as a percentage.
𝑉
Eq.2
𝑛
Void ratio and porosity are inter-related to each other as follows:
𝑒 =
1−𝑛
𝑛=
𝑒
1+𝑒
And
3. The volume of water (Vw) in a soil can vary between zero (i.e. a dry soil) and the volume
𝑆
=
𝑉𝑤
of voids. This can be expressed as the degree of saturation (S) in percentage.
𝑉𝑣
- Eq.3
For a dry soil, S = 0%, and for a fully saturated soil, S = 100%.
4. Air content (ac) is the ratio of the volume of air (Va) to the volume of voids.
𝑎
𝑎
------------------------------- Eq.4
𝑐 𝑉𝑣
n = × 100 = n ×
Va
5. Percentage air voids (na) is the ratio of the volume of air to the total volume.
a c
a
-Eq.5
V
Density is a measure of the quantity of mass in a unit volume of material. Unit weight is a
measure of the weight of a unit volume of material. Both can be used interchangeably. The units
of density are ton/m³, kg/m³ or g/cm³. The following are the basic weight relations:
1. ` The ratio of the mass of water present to the mass of solid particles is called the water
𝜔=
𝑊𝜔
content (w), or sometimes the moisture content.
----------------------- Eq.6
𝑊𝑠
Its value is 0% for dry soil and its magnitude can exceed 100%.
2. The mass of solid particles is usually expressed in terms of their particle unit
𝛾 = = 𝐺 𝛾 ----------------
𝑊𝑠
weight (𝛾𝑠) or specific gravity (Gs) of the soil grain solids.
𝑠 𝜔
Eq.7
𝑠
𝑉𝑠
where, 𝛾𝜔= unit weight of water
For most inorganic soils, the value of Gs lies between 2.60 and 2.80. The presence of organic
material reduces the value of Gs.
3. Dry unit weight (γd) is a measure of the amount of solid particles per unit volume.
𝛾--------------------
𝑠
𝑉
Eq.8
𝑠
4. Bulk unit weight (𝜸𝒐𝒓 𝜸𝒕) is a measure of the amount of solid particles plus water
per unit volume.
𝛾
= 𝑠𝛾 =𝑤
𝑤
𝑊 +𝑊
𝑉𝑠+𝑉
𝑡
7
Lecture Note on GTE-I
- ---------- Eq.9
5. Saturated unit weight (𝛾𝑠𝑎𝑡) is equal to the bulk density when the total voids is filled up
with water.
8
Lecture Note on GTE-I
6. Buoyant unit weight or submerged unit weight (γ′) is the effective mass per unit volume
𝛾; = 𝛾𝑠𝑎𝑡 − 𝛾𝑤-----------------------
when the soil is submerged below standing water or below the ground water table.
Eq.10
It is important to quantify the state of a soil immediately after receiving in the laboratory and prior to
commencing other tests. The water content and unit weight are particularly important, since they
may change during transportation and storage.
Some physical state properties are calculated following the practical measurement of others. For
example, dry unit weight can be determined from bulk unit weight and water content. The
following are some inter-relations:
1.5 Relationships Involving Porosity, Void Ratio, Degree of Saturation, Water
Content, Percent Air Voids and Air Content
𝑛 = , as a fraction
𝑉𝑣
Fig.3: Three phase diagram used for establishment of inter relationship
𝑉
= =1−
𝑉−𝑉𝑠
𝑉
--------------------------
𝑠 𝑊𝑠
=1-
𝑉 𝑉 𝐺𝛾𝜔𝑉
𝑛=1−
Eq.11
𝑊𝑑 -------------------------------------- Eq.12
𝐺𝛾𝜔𝑉
𝑒= = = −1= − 1------------------
𝑉𝑣 𝑉−𝑉𝑠 𝑉 𝑉𝐺𝛾𝜔
This may provide a practical approach to determine n
Eq.13
𝑉𝑠 𝑉𝑠 � 𝑊𝑠
𝑠
�
𝑒 = − 1---------------------
𝑉𝐺𝛾𝜔
Eq.14
𝑊𝑑
𝑛
= and 𝑒 =
𝑉𝑣 𝑉𝑣
The expression provides an approach to determine value of void ratio e
𝑉 𝑉
𝑉𝑠+𝑉𝑣 𝑠
= = = + = +1=
1 𝑉 𝑉𝑠 𝑉𝑣 1 1+𝑒
---------------------- Eq.15
𝑛 𝑉 𝑉𝑣 𝑉𝑣 𝑒 𝑒
𝑉𝑣
𝑒𝑣
So 𝑛 =
1+𝑒
------------------------------ Eq.16
a =
V
a
Vv
c
9
Lecture Note on GTE-I
and 𝑛 =
𝑉𝑣
∴ = = ------------------------------
𝑉𝑎
𝑛
Eq.17
𝑛𝑎𝑐 𝑉 �
�
By definition
10
Lecture Note on GTE-I
𝜔=
𝑊𝑤
--------------------- Eq.18
𝑊𝑠
𝑆=
𝑉𝑤
--------------------- Eq.19
𝑉𝑣
And 𝑒 =
𝑉𝑣
------------------- Eq.20
𝑉𝑠
Now , 𝑆𝑒 =
𝑉𝑤
------------------ Eq.21
𝑉𝑠
Again, 𝜔 = = = = = 𝑆𝑒/𝐺-----------
𝑊 𝑤 𝑉𝑤𝛾𝑤 𝑉𝑤𝛾𝑤 𝑉𝑤
Eq.22
𝑊𝑠 𝑉𝑠 𝑉 𝑠𝐺 𝑉𝑠𝐺
𝛾𝑠 𝛾𝑤
Hence 𝜔𝐺 = 𝑆𝑒---------------------------- Eq.23
1.6 Relationships Involving Unit Weights, Grain Specific Gravity, Void Ratio, and
Degree of Saturation:
Ww⁄
Ws(1+
Ws)
γ= = =
W Ws+Ww Eq.24
Vv
V Vs(1 ⁄
--------------------
+ )
Vs+Vv Vs
𝑊𝑠 𝑉𝑠
Again 𝛾 = = 𝐺𝛾 ----------------------------------
𝑊𝑠
Eq.25
𝑠 𝑤
Hence , 𝛾 = 𝐺𝛾
𝑠 𝑉 1+𝜔
𝑤 1+𝑒
---------------------------------- Eq.26
Or 𝛾 = (𝐺+𝐺𝜔)--------------------------------------
𝛾𝑤
(1+𝑒)
Eq.27
This is a general expression from which unit weight corresponding to the saturated and dry states
∴
of soil mass may be found out by substituting S=1 and S=0 respectively.
𝛾
=
(𝐺+𝑒) --------------------
𝛾
Eq.1.30
𝑠𝑎𝑡 (1+𝑒) 𝑤
𝛾 ----------------------------
= (1+𝑒
𝐺
𝛾𝑑
And Eq.1.31
) 𝑤
𝛾′ (𝐺+𝑒)
= 𝛾𝑠𝑎𝑡 − 𝛾𝑤-------------------------
The submerged unit weight may be written as
= 𝛾 − 𝛾 ----------------------
Eq.1.32
𝑤 𝑤
Eq.1.33
(1+𝑒)
= [(𝐺+𝑒) − 1]---------------------------
𝛾𝑤 (1+𝑒)
Eq.1.34
So 𝛾′ = 𝛾 ---------------------------
(𝐺−1)
(1+𝑒) 𝑤
Eq.1.35
11
Lecture Note on GTE-I
Example 1: A soil has void ratio = 0.72, moisture content = 12% and Gs= 2.72. Determine its
(a) Dry unit weight
12
Lecture Note on GTE-I
𝛾 = =
𝐺𝑠𝛾𝜔 2.72×9.81
Solution:
𝑑
a) =15.51 kN/m3
1+0.72
1+𝑒
𝛾 = 𝛾𝑑(1 + 𝜔)
b)
1+0.72
= × 9.81=19.62kN/m3
2.72+0.72
= ×
�
𝐺𝑠+𝑒𝑆
c)
𝛾
�
𝑠𝑎𝑡 𝜔 1+0.72
1+
𝑒
𝛾𝑠𝑎𝑡 − 𝛾
Water to be added per m3 to make the soil saturated
Example 2: The dry density of sand with porosity of 0.387 is 1600 kg/m3. Find the void ratio of
= 1000
𝑘𝑔
𝛾𝜔
the soil and the specific gravity of the soil solids. [Take ]
𝑚3
n = 0.387
= 1600 kg/m3
Solution:
𝑛 0.387
𝑒 = = = 0.631
1−𝑛 1 − 0.387
𝐺𝑠𝛾𝜔
𝛾𝑑 =
1+𝑒
𝐺 = ×𝛾 = × 1600 = 2.61
1+𝑒 1+0.631
𝑠 𝛾𝜔 1000
𝑑
Example 3
A sample of saturated soil has a water content of 35%. The specific gravity of soil solids is 2.65.
Determine its voids ratio, porosity, saturated unit weight and dry unit weight.
13
Lecture Note on GTE-I
Example 4
A sample of clay taken from a natural stratum was found to be partially saturated and when
tested in the laboratory gave the following results. Specific gravity of soil particle =2.6, wet
weight of sample =2.50 N, dry weight of sample =2.10 N and volume of sample =150 cm 3.
Determine the degree of saturation.
Sol:
Data given
14
Lecture Note on GTE-I
15
Lecture Note on GTE-I
Problem No 1 A soil sample has a porosity of 40 %. The specific gravity of soil solids is 2.70.
Calculate a) Voids ratio, b) dry density, c) unit weight if the soil is 56% saturated and d) unit
weight if the soil is completely saturated
16
Lecture Note on GTE-I
2 Engineering Properties
a. Permeability
b. Shear Strength
c. Compressibility
2.1.1 SPECIFIC GRAVITY
Definition:
The ratio of weight of soil solids to the weight of equal volume of water content at 40C standard
temperature
Laboratory Determination:
a. Pycnometer Method
17
Lecture Note on GTE-I
(𝑊4−𝑊1)
IS Sieves
18
Lecture Note on GTE-I
SIEVE ANALYSIS
Objective of sieve analysis is to classify the soil and to determine gradation.
A set of IS Sieves are taken
A suitable amount of pulverized dry soil (500g) is taken and sieved through different sieves with
the help of sieve shaker for 10 – 15 minutes.
Weight of soil retained on each sieve is determined.
Percentage weight of soil retained on each sieve is
determined. Percentage finer is determined.
Plot of Percentage finer v/s Size of the sieve on semi log sheet will give the type of gradation,
Effective size, Coefficient of uniformity and Coefficient of curvature,
TABULATION FOR SIEVE ANALYSIS
IS Sieve Weight of Soil Percentage Cumulative Percentage Finer
Size Retained(gm) Retained Percentage (100% - Cumulative Retained)
4.75mm
2mm
1mm
600μ
425μ
212μ
150μ
75μ
Pan
19
Lecture Note on GTE-I
1.
Coefficient of Uniformity (Cu)
𝐶
60
D10
----------------------- Eq.2.3
𝑢
2.
Coefficient of Curvature (Cc)
𝐷2
𝐶𝐶 30
=
𝐷60×𝐷
------------------------- Eq.2.4
10
Well graded soils have high Cu values and poorly graded soils have low Cu, values. If all the
particles of soil mass are of the same size Cu is unity.
Cc lies between 1 to 3 for well graded soil.
Uniform soil: Cu = 1
Poorly graded soil: 1< Cu<4
Well graded soil: Cu>4
2.1.5 HYDROMETER ANALYSIS
1. The analysis is carried out in two stages
• Calibration of hydrometer
• Sedimentation analysis
CALIBRATION OF HYDROMETER
A = Area of stem or graduate
20
Lecture Note on GTE-I
Rh = Hydrometer Reading
He = Effective Height
H = Depth of stem up to hydrometer reading
𝐻 = (𝐻 + + ℎ) − ℎ----------------------
ℎ 𝑉 𝑉
Vh = Volume of hydrometer in cc
𝑒
Eq.2.5
2𝐴 𝐴
2
Hydrometer is immersed in the soil suspension, hydrometer readings are noted at different time
intervals
Effective height values are calculated for corresponding ‘H’ values.
A plot of He v/s Rh is used as the calibration chart.
Note that calibration of hydrometer is done only once. The same chart can be used for any ‘n’
number of trials.
21
Lecture Note on GTE-I
TEST PROCEDURE
×
He
𝐷=
18𝜂 √ t
------------------------- Eq.2.6
(ϒ𝑠−ϒ𝑤
)
22
Lecture Note on GTE-I
𝐺−1
) × 100-----------------------
𝑅𝑐
Eq.2.7
𝑊𝑠
23
Lecture Note on GTE-I
Where,
N - % Finer
G - Specific Gravity of Solids
Rc - Corrected reading of hydrometer
Ws -Weight of solids
24
Lecture Note on GTE-I
𝛾
𝛾𝑑
Dry Density
1+𝑤
= --------------------- Eq.2.9
𝛾𝑠𝑎𝑛𝑑
SPC with sand up to W1 is taken to the field where bulk unit weight is required. The surface is
cleaned and leveled.
The tray with central hole is placed in position. A hole of approximately same size as calibrating
container is dug
The removed soil is collected on the sides of tray and its weight is noted (Wsoil).
The SPC is placed on the hole and sand is allowed to run in to hole.
Weight of SPC is noted (W4).
Weight of sand in the holeW5 =W1 –W4 +Wcone
Volume of the hole
𝑉𝑠 =
𝑊𝑆
𝛾𝑠𝑜𝑖𝑙
---------------------- Eq.2.11
𝑜𝑖𝑙
25
Lecture Note on GTE-I
26
Lecture Note on GTE-I
It is the arbitrary water content between plastic and liquid phase of soil
It is the minimum water content in the soil where it behaves as a plastic material
It is the amount of water in the soil, where the soil contains least amount of shear strength.
2.3.2 Plastic Limit
It is the minimum water content in the soil at which a thread formed by rolling on a glass plate
crumbles at 3 mm diameter.
In other words, it is the arbitrary water content between plastic and semi solid state.
2.3.3 Shrinkage Limit
It is the maximum water content at which soil behaves as a solid material.
It is the maximum water content below which, decrease in water content results in no change of
volume of soil.
In other words, it is the arbitrary water content between solid and semi solid phase of soil.
2.4 DETERMINATION OF LIQUID LIMIT
Casagrande’s Liquid Limit Apparatus
The following procedures are used
27
Lecture Note on GTE-I
1 Take 120 g of air dried soil after passing through 425μ sieve.
2 Arbitrarily mix certain amount of water and make a paste.
3 Put the soil pat in the brass cup and make a V shaped groove using grooving tool.
4 Give blows at a speed of 2 rps., till the grove closes.
5 Count the number of blows at the end of experiment.
6 Take certain amount of soil from the cup to determine water content.
7 Take at least 5 number of trials with varying number of blows from 15 to
35 The corresponding water content at 25th blow is taken as liquid limit
Procedures to be taken
1 Take about 15 g of air dried soil after passing through 425μ sieve.
2 Arbitrarily mix certain amount of water and make soil ball.
3 Roll the soil ball on glass plate with palm.
4 Make thread till it attains a thickness of 3mm.
28
Lecture Note on GTE-I
5 If the thread crumbles at 3 mm thick, then take sample to determine water content.
6 Corresponding water content is taken as PLASTIC LIMIT.
2.6 DETERMINATION OF SHRINKAGE LIMIT
Procedure to determine Shrinkage Limit
1 Take about 30 g of air dried soil after passing through 425μ sieve in an evaporating dish.
2 Soil is mixed with sufficient quantity of water so that it may flow.
3 Fill the soil mix in dish and take the weight of dish.
4 The soil pat is allowed to dry in air till the colour of pat changes from dark to light.
5 The dish is then placed in oven at 110oc till it attains constant weight.
29
Lecture Note on GTE-I
=
(𝑊1−𝑊2)−(𝑉1−𝑉2)𝛾𝑤
𝑤𝑠
Hence, Shrinkage Limit
------------------------ Eq.2.13
𝑊𝑆
It is the range of water content in soil in during which soil behaves as a plastic material
Simply it is the difference in the water content at liquid limit and plastic limit
Highly plastic soils exhibit greater plasticity values
It is the range of water content in soil in during which soil behaves as a semi solid material
Simply it is the difference in the water content at plastic limit and shrinkage limit
It is the slope of the flow curve obtained between the number of blows and water content in the
Casagrande’s test for the determination of liquid limit.
Flow index indicates the loss of shear strength of soil upon the increase of water.
Greater the slope angle less is the shear strength of soil.
Flatter slope curve indicates more shear strength in the soil.
Flow index = (𝑊1−𝑊2)------------------------
𝐼𝐹 𝑙𝑜𝑔10𝑛
Eq.2.16
2
𝑛1
2.7.4 Consistency Index (Ic)
It is the ratio of the difference between the liquid limit and natural water content of soil to its
plasticity index. This is very important to attribute the soil property
=
(𝑤𝑙−𝑤𝑛)
𝐼𝐶
Consistency Index ------------------- Eq.2.17
𝐼𝑃
30
Lecture Note on GTE-I
It is the ratio of the difference between the natural water content of soil and its plastic limit to its
plasticity index.
Liquidity index varies from 0 to 1
(𝑤𝑛−𝑤𝑝)
Liquidity Index, IL = ------------------ Eq.2.18
𝐼𝑃
2.7.6 Toughness Index (IT)
It is the ratio of plasticity index to the flow index.
This gives an idea of shear strength of soil at plastic limit
𝐼𝑃
Toughness Index, IT = -------------------------- Eq.2.19
𝐼𝐹
2.7.7 Density Index (ID)
It is the ratio of difference between the void ratio of a Cohesionless soil in the loosest state
and void ratio in its natural state to the difference between its void ratios in the loosest and
(𝑒𝑚𝑎𝑥−𝑒𝑛𝑎𝑡)
densest states.
(𝑒𝑚𝑎𝑥−𝑒𝑚𝑖𝑛)
Hence Density Index ID= ------------------------- Eq.2.20
The nature of soil can be predicted from Density Index value as given in table
Density Index (ID) Soil Description
<15 Very loose 15 - 35 Loose
35 - 65 Medium 65- 85 Dense
>85 Very dense 5/
<15 Very loose 15 - 35 Loose
35 - 65 Medium 65- 85 Dense
2.7.8 Shrinkage Ratio (SR)
It is the ratio of a given volume change in a soil, expressed as a percentage of the dry volume, to
the corresponding change in water content above the shrinkage limit.
31
Lecture Note on GTE-I
× 100-----------------------
𝑉1−𝑉2
SR=( 𝑉𝑑 )
𝑤1−𝑤2
Eq.2.21
Where
V1 = Volume of soil at water content w1 (%)
V2 = Volume of soil at water content w2 (%)
Vd = Volume of dry soil mass
At Shrinkage,
V2 = Vd
w2= ws
𝑉1−𝑉𝑑
(
Shrinkage Ratio (SR) = 𝑉𝑑
) × 100-------------------- Eq.2.22
𝑤1−𝑤𝑠
2.8 MISCELLANEOUS PROPERTIES
2.8.1 Sensitivity (St)
It is defined as the ratio of unconfined compressive strength of an undisturbed soil sample to the
unconfined compressive strength of the same soil after remolding at unaltered water content.
Sensitivity Soil Description
1-4 Normal
4-8 Sensitive
8 - 15 Extra sensitive
>15 Quick
2.8.2 Activity
It is the ratio of plasticity index of soil to the % finer fraction available in it.
Clayey soils exhibit more activity values.
Activity Classification Activity Classification
<0.75 Inactive
0.75 – 1.25 Normal
>1.25 Active
32
Lecture Note on GTE-I
Sol:
Data given
33
Lecture Note on GTE-I
34
Lecture Note on GTE-I
A soil particle may be a mineral or a rock fragment. A mineral is a chemical compound formed
in nature during a geological process, whereas a rock fragment has a combination of one or more
minerals. Based on the nature of atoms, minerals are classified as silicates, aluminates, oxides,
carbonates and phosphates.
Out of these, silicate minerals are the most important as they influence the properties of clay
soils. Different arrangements of atoms in the silicate minerals give rise to different silicate
structures.
Soil minerals are formed from two basic structural units: tetrahedral and octahedral. Considering
the valencies of the atoms forming the units, it is clear that the units are not electrically neutral
and as such do not exist as single units.
The basic units combine to form sheets in which the oxygen or hydroxyl ions are shared among
adjacent units. Three types of sheets are thus formed, namely silica sheet, gibbsite
sheet and brucite sheet.
Isomorphous substitution is the replacement of the central atom of the tetrahedral or octahedral
unit by another atom during the formation of the sheets.
The sheets then combine to form various two-layer or three-layer sheet minerals. As the basic
units of clay minerals are sheet-like structures, the particle formed from stacking of the basic
units is also plate-like. As a result, the surface area per unit mass becomes very large.
A tetrahedral unit consists of a central silicon atom that is surrounded by four oxygen atoms
located at the corners of a tetrahedron. A combination of tetrahedrons forms a silica sheet.
35
Lecture Note on GTE-I
An octahedral unit consists of a central ion, either aluminium or magnesium, that is surrounded
by six hydroxyl ions located at the corners of an octahedron. A combination of aluminium-
hydroxyl octahedrons forms a gibbsite sheet, whereas a combination of magnesium-hydroxyl
octahedrons forms a brucite sheet.
36
Lecture Note on GTE-I
Montmorillonite and illite clay minerals are the most common. A basic three-layer sheet unit
is formed by keeping one silica sheet each on the top and at the bottom of a gibbsite sheet.
These units are stacked to form a lattice as shown.
3.1.3 Montmorillonite Mineral
The bonding between the three-layer units is by vander Waals forces. This bonding is very weak
and water can enter easily. Thus, this mineral can imbibe a large quantity of water causing
swelling. During dry weather, there will be shrinkage.
Illite consists of the basic Montmorillonite units but are bonded by secondary valence
forces and potassium ions, as shown. There is about 20% replacement of aluminium with
silicon in the gibbsite sheet due to isomorphous substitution. This mineral is very stable and
does not swell or shrink.
37
Lecture Note on GTE-I
Specific surface area (SSA) is defined as the surface area of soil particles per unit mass (or
volume) of dry soil. Its unit is in m2/g or m2/m3. Clay minerals are characterized by high specific
surface area (SSA) as listed in Table 2.2. High specific surface area is associated with high soil-
water-contaminant interaction, which indicates high reactivity. The reactivity increases in the
order Kaolinite< Illite< Montmorillonite. For the purpose of comparison, SSA of silt and sand
has also been added in the table. There is a broad range of SSA values of soils, the maximum
being for montmorillonite and minimum for sand. As particle size increases SSA decreases.
For smectite type minerals such as montmorillonite, the primary external surface area amounts to
50 to 120 m2/g. SSA inclusive of both primary and secondary surface area (interlayer surface
area exposed due to expanding lattice) and termed as total surface area would be closed to 800
m2/g.For kaolinite type minerals there is possibility of external surface area where in the
interlayer surface area does not contribute much.
3.2.2 Plasticity and cohesion
Clay attracts dipolar water towards its surface by adsorption. This induces plasticity in clay.
Therefore, plasticity increases with SSA. Water in clays exhibits negative pressure due to which
two particles are held close to each other. Due to this, apparent cohesion is developed in clays.
Surface charge and adsorption
Clay surface is charged due to following reasons:
Iso-morphous substitution (Mitchell and Soga 2005): During the formation of mineral, the
normally found cation is replaced by another due to its abundant availability. For example: when
Al+3 replace Si+4 there is a shortage of one positive charge, which appears as negative charge on
38
Lecture Note on GTE-I
clay surface. Such substitution is therefore the major reason for net negative charge on clay
surface.
O-2 and OH- functional groups at edges and basal surface also induce negative charge.
Dissociation of hydroxyl ions or broken bonds at the edges is also responsible for unsatisfied
negative or positive charge. Positive charge can occur on the edges of kaolinite plates due to
acceptance of H+ in the acid pH range (Berkowitz et al. 2008). It can be negatively charged under
high Ph environment.
Absence of cations from the crystal lattice also contributes to charge formation.
In general, clay particle surface are negatively charged and its edges are positively charged.
Due to the surface charge, it would adsorb or attract cations (+ve charged) and dipolar molecules
like water towards it. As a result, a layer of adsorbed water exists adjacent to clay surface.
3.2.3 Exchangeable cations and cation exchange capacity
Due to negative charge, clay surface attracts cations towards it to make the charge neutral. These
cations can be replaced by easily available ions present in the pore solution, and are termed as
exchangeable ions. The total quantity of exchangeable cations is termed as cation exchange
capacity, expressed in milli equivalents per 100 g of dry clay. Cation exchange capacity (CEC) is
defined as the unbalanced negative charge existing on the clay surface. Kaolinite exhibits very
low cation exchange capacity (CEC) as compared to montmorillonite. Determination of CEC is
done after removing all excess soluble salts from the soil. The adsorbed cations are then replaced
by a known cation species and the quantity of known cation required to saturate the exchange
sites is determined analytically.
When two clay particles come closer to each other it experiences (a) interparticle attraction due
to weak van-der-Waal’s force (b) repulsion due to –ve charge. When particles are sufficiently
close, attraction becomes dominant active force and hence there is an edge to face configuration
for clay particles as shown in Fig. 3.5 a. Such a configuration is termed as flocculant structure.
When the separation between clay particles increase, repulsion becomes predominant and hence
the clay particles follows face to face configuration called dispersed structure. A lot of micro and
macro level behaviour of clays are associated with these arrangement of clay particles (Mitchell
and Soga 2005).
39
Lecture Note on GTE-I
Other than silica and clay, subsurface contains a variety of minerals such as oxides and
carbonates that governs the reactivity of soil and its interaction with the environment. Some of
the abundant metal oxide minerals present are iron oxides (hematite, magnetite, goethite etc.) and
aluminium oxides (gibbsite, boehmite). Other oxide minerals (such as manganese oxide, titanium
oxide) are far less than Fe and Al oxides, but because of small size and large surface area, they
would affect very significantly the geochemical properties of subsurface. These oxides are
mostly present in residual soils of tropical regions. Other major components include soluble
calcium carbonate and calcium sulphate, which has relatively high surface area. In most soils,
quartz is the most abundant mineral, with small amount of feldspar and mica present. Carbonate
minerals such as calcite and dolomite are found in some soils in the form of bulky particles,
precipitates etc. Sulphate minerals mainly gypsum are found in semiarid and arid regions.
One of the very well established methods for mineralogy characterization of fine-grained soils is
by using X-ray diffraction (XRD) analysis. Majority of the soil minerals are crystalline in nature
and their structure is defined by a unique geometry. XRD identifies minerals based on this
unique crystal structure. In XRD, characteristic X-rays of particular wave length are passed
through a crystallographic specimen. When X-ray interacts with crystalline specimen it gives a
particular diffraction pattern, which is unique for a mineral with a particular crystal structure.
The diffraction pattern of the soil specimen (according to its crystal structure), which is based on
powder diffraction or polycrystalline diffraction, is then analyzed for the qualitative and
quantitative (not always) assessment of minerals. Sample preparation method for XRD should be
done with great care as the XRD reaches only a small layer (nearly 50 µm) from the surface of
the sample. Hence, homogeneity is very important. Soil sample is initially dried and sieved
through 2 mm sieve. Sieved sample is homogenized in a tumbler mixer for 30 min. A control
mix of 30 g was taken and ground in lots of 15 g in a gyratory pulverizer. 15 percent by weight
40
Lecture Note on GTE-I
of KIO4 (internal standard) was added to 5 g of specimen and again homogenized in a mixer..
X-ray wave of monochromatic radiation (K α) is commonly obtained from copper radiation,
which is commonly known as Cu- Kα. A typical XRD output is represented by Fig. 2.5. It can be
noted from the figure that ordinate represent relative intensity of X-ray diffraction and abscissa
represents twice of angle at which a striking X-ray beam of wave length λ makes with parallel
atomic planes. Based on this diffraction pattern, the minerals can be identified by matching the
peak with the data provided by International Centre Diffraction Data (ICDD) formerly known as
Joint Committee on Powder Diffraction Standards (JCPDS).
It is understood that the area under the peak of diffraction pattern gives the quantity of each
phase present in the specimen. However, quantitative determination of mineral composition in
soils based on simple comparison of diffraction peak heights or area under peak is complex and
uncertain because of different factors such as mineral crystallinity, hydration, surface texture of
the specimen, sample preparation, non-homogeneity of soil samples, particle orientation etc. The
method of quantification will be more precise for those soils with less number of minerals. Al-
Rawas et al. (2001) have discussed about constant mineral standard method and constant clay
method for quantification of clay minerals
Fig.3.6: A typical XRD pattern with mineral identification for two different soils (modified from
Sreedeep 2006)
In the first method, increasing quantity of clay are added to the fixed mass of known standard
and the difference in X-ray diffraction intensity when the specimen changes from 100 percent
standard to 100 percent clay is noted. The peak area ratio for each component is then plotted
against percentage of clay, based on which regression equation is determined. This regression
equation is further used for mineral quantification. In the second method, known weight of pure
standard mineral is added to clay containing the same components, and the change in the
41
Lecture Note on GTE-I
reflection peak-area intensity of each component is measured to estimate the weight proportion
of that component.
42
Lecture Note on GTE-I
Since the difference between wL and wp is indicative of the compressibility of the soil we can
argue that a Na Kaolinite would be more compressible than a Ca Kaolinite. When the double
layer and mineral particle size are considered together it can be argued that the smaller the
particle (and hence the higher the specific surface) the larger the amount of water that (with
respect to the mass of dry soil) will be associated with the double layer. This feature is of great
importance when the swell of a soil, when wetted from a dry condition, is considered. Of
particular significance is the clay mineral montmorillonite: the bonding between successive
layers is weak and when an ion associated with this interlayer position hydrates, the clay mineral
stack will expand. Expansion associated with moisture can cause structural failure in some
lightly loaded structures.
3.8 STRUCTURE OF SOILS
The ‘structure’ of a soil may be defined as the manner of arrangement and state of aggregation of
soil grains. In a broader sense, consideration of mineralogical composition, electrical properties,
orientation and shape of soil grains, nature and properties of soil water and the interaction of soil
water and soil grains, also may be included in the study of soil structure, which is typical for
transported or sediments soils. Structural composition of sedimented soils influences, many of
their important engineering properties such as permeability, compressibility and shear strength.
Hence, a study of the structure of soils is important.
The following types of structure are commonly studied:
(a) Single-grained structure
(b) Honey-comb structure
(c) Flocculent structure
3.8.1 Single-grained Structure
Single-grained structure is characteristic of coarse-grained soils, with a particle size greater than
0.02mm. Gravitational forces predominate the surface forces and hence grain to grain contact
results. The deposition may occur in a loose state, with large voids or in a sense state, with less
of voids
43
Lecture Note on GTE-I
41
Lecture Note on GTE-I
42
Lecture Note on GTE-I
43
Lecture Note on GTE-I
44
Lecture Note on GTE-I
The standard recommends that when a soil possesses characteristics of two groups either in
particle size distribution or in plasticity, it is designed by combination of group symbols.
Table 3: Unified Soil Classification System (based on particle passing through 75 mm sieve)
45
Lecture Note on GTE-I
46
Lecture Note on GTE-I
𝐶𝑢 > 4
Symbol Names
𝐶𝑐between 1 &3
Coarse Gravel Clean GW Well graded When
grained soil More than gravel gravel fines are
(More than half of (Fine less GP Poor graded Not meeting all between
half of material is than 5%) gravel requirement for GW
5% to
material is larger than Gravel with GM Silty gravel Atterberg Atterberg
12%
larger than the 4.75 mm fines (Fines limit limits
the 75 – IS sieve more than below A plotting boarder
𝐼𝑝 > 7 line
line or boarder
symbol
cases GW-
requiring GP.SW-
use of SC etc
dual
symbol
𝐶𝑢 > 6
GM-GC
𝐶𝑐between 1 &3
Sands Clean sand SW Well graded
More than (Fine less sands
half of than 5%) SP Poor graded Not meeting all
material is sands requirement for SW
smaller Sands with SM Silty Sands Atterberg Atterberg
than the fines (Fines limit limits
4.75 mm IS more than below A plotting
sieve
𝐼𝑝 < 4 line with
12%) line or above A
𝐼𝑝 > 7 line
line or boarder
cases
requiring
use of
dual
symbol
SM-SC
47
Lecture Note on GTE-I
𝐼𝑝 < 7
soil (L)Liquid limit plasticity or plotting cotton soils
(More less than 35% in India
above A line
than CL Inorganic clays Atterberg limit lies partly
half of with Ip above and
with low above A line
𝐼𝑝 > 7
material plasticity or between 4 & partly
is 7 (hatched below the
smaller zone) ML- A-Line
than the CL
75 µ IS
sieve ) OL Organic silts Atterberg limit below A line Organic
with low and
plasticity inorganic
MI Inorganic silts Atterberg limit below A line soils
Intermediate with medium plotted in
compressibility plasticity same zone
(I) CI Inorganic clays Atterberg limit above A line of plasticity
Liquid limit with medium chart are
more than 35% plasticity distinguish
but less than OI Organic silts Atterberg limit below A line ed by
50% with medium colour,
plasticity odor or
MH Inorganic silts Atterberg limit below A line liquid limit
High with High test after
compressibility plasticity oven
(H) CH Inorganic clays Atterberg limit above A line drying.
Liquid limit with High Reduction
more than 50% plasticity in value
after and
OH Organic silts Atterberg limit below A line
before
with medium to
drying by
High plasticity
three fourth
indicates
clear
organic
soils
49
Lecture Note on GTE-I
SOIL WATER: Water presence in the voids of soil mass is called soil water. It can be
classified in several ways:
Structural water: It is the water chemically combined in the crystal structure of the soil mineral.
Structural water cannot be separated or removed and also not removed by oven drying at 105-
50
Lecture Note on GTE-I
110°c.It can be destroyed at higher temperature which will destroy the crystal structure.
Infiltrated water Infiltrated water is the portion of surface precipitation which soaks into ground,
moving downwards through air containing zones.
Pore Water: It is cable of moving under hydrodynamic forces unless restricted in its free
movement such as when entrapped between air bubbles or retention by capillary forces.
Gravitational and capillary water are the two types of pore water.
Solvate Water: The water which forms a hydration shell around soil grains is solvate water. it
is subjected to polar electrostatic and binding forces.
Ground water: Subsurface water that fills the voids continuously and is subjected to no force
other than gravity is known as gravitational water.
Capillary water: The minute pores of soil serve as capillary tubes through which the moisture
rise above the ground water table. Capillary water is the soil moisture located within the
interstices and voids of capillary size of the soil. Capillary water is held in the interstices of soil
due to capillary forces. Capillary action or capillarity is the phenomenon of movement of water
in the interstices of a soil due to capillary forces. The capillary forces depend upon various
factors such as surface tension of water, pressure in water in relation to atmospheric pressure,
and the size and conformation of soil pores.
The pores of soil mass may be looked upon as a series of capillary tubes, extending vertically
above water table. The rise of water in the capillary tubes, or the fine pores of the soil, is due to
the existence of surface tension which pulls the water up against the gravitational force. The
height of capillary rise, above the ground water (or free water) surface depends upon the
diameter of the capillary tube (or fineness of the pores) and the value of the surface tension.
When a capillary tube is inserted in water, the rise of water will take place up to reach the
equilibrium. At this stage the rise of water in the tube is stopped.
which water is retained in a soil mass. It decreases linearly from a maximum value of ℎ𝑐𝛾𝑤 at
pressure reduction or negative pressure in the pore water (or the pressure below atmospheric) by
the level of the meniscus to zero value at the free water surface. The pressure deficiency in the
held water is also termed as soil suction or suction pressure. Soil suction is measured by the
height hc in centimeters to which a water column could be drawn by suction in a soil mass free
from external stress. The common logarithm of this height (cm) or pressure (g/cm2) is known as
the pF value (Schofield, 1935): pF = log10 (hc) Thus, a pF value of 2 represents a soil suction of
100 cm of water or suction pressure and capillarity of 100 g /cm2.
51
Lecture Note on GTE-I
The magnitude of the pressure is the same at all height above the free water surface. The
capillarity pressure transferred from grain to grain called as inter angular or effective pressure.
Capillary action (or) capillarity: It is the phenomenon of movement of water in the interstices
of a soil due to capillary forces. The capillary forces depend upon various factors depend upon
various factors such as surface tension of water, pressure in water in relation to atmospheric
pressure and three size and conformation of soil pores.
52
Lecture Note on GTE-I
In contrast to capillary tubes the continuous voids in soils have a variable width. They
communicate with each other in all directions and constitute an intricate network of voids. When
water rises into the network from below, the lower part of the network becomes completely
saturated. In the upper part, however, the water occupies only the narrowest voids and the wider
areas remain filled with air. Fig. shows a glass tube filled with fine sand. Sand would remain
fully saturated only up to a height h' which is considerably smaller than h c. A few large voids
hc only in the smaller voids. The zone between the depths ( ℎ� 𝑐 − ℎ′ ) will remain
may effectively stop capillary rise in certain parts. The water would rise, therefore, to a height of
�
partially
saturated.
The capillary phenomenon which enables raising the water in capillary tube against the force of gravity
can also be taken place in voids of dry soil. In case of the earth dam, this may create serious
problems. Capillary forces are able to raise water against the force of gravity not only into
capillary tubes or the voids in columns of dry soil, but also into narrow open channels or V-
shaped grooves. If the highest point of the groove is located below the level to which the surface
tension can lift the water, the capillary forces will pull the water into the descending part of the
groove and will slowly empty the vessel. This process is known as capillary siphoning. The same
process may also occur in the voids of soil. For example, water may flow over the crest of an
impermeable core in a dam in spite of the fact that the elevation of the free water surface is
below the crest of the core as shown in
In an earth dam, the central impervious core provides to check the seepage of water through the body of
dam. When the water level in the reservoir reaches HFL, the portion of the upstream of the dam
will be saturated. The water level in u/s of the pervious shell will be practically the same as HFL.
Now due to capillarity, the water will raise through the height ‘h’. If the top of the impervious
core is situated at height y<h, above the HFL, the capillary force will pull the water to the
downstream of the dam. This effect will slowly empty the reservoir. This process is known as
Capillary Siphoning action in earth dam.
This will not only result in the water loss but also damage the downstream and create instability to the
earth dam. So the height of the impervious core is a major factor for design consideration.
Contact moisture. Water can also be held by surface tension round the point of contact of two particles
(spheres) capillary water in this form is known as contact moisture (or) contact capillary water.
6.0 EFFECTIVE STRESS CONCEPTS IN SOIL
At any plane in a soil mass, the total stress or unit pressure σ is the total load per unit area. This pressure
may be due to i) self weight of soil ii) over burden on the soil. The total pressure consists of two
distinct components: inter granular pressure or effective pressure and the neutral pressure or pore
pressure. Effective pressure σ' is the pressure transmitted from particle through their point of
contact through the soil mass above the plane. Such a pressure, also termed as inter
54
Lecture Note on GTE-I
granular pressure, is effective in decreasing the voids ratio of the soil mass and in mobilizing its
shear strength. The neutral pressure or the pore water pressure or pore pressure is the pressure
transmitted through the pore fluid. Therefore, this pressure is also called neutral pressure (u).
Since the total vertical pressure at any plane is equal to the sum of the effective pressure and pore
𝜎 = 𝜎′ + 𝑢-------------------
water pressure we have
Eq.6.1
55
Lecture Note on GTE-I
𝑢 = 𝛾𝑤ℎ--------------------------------
is given by
Eq.6.2
Where, h = depth below water table or overlying water surface
It is convenient to think of pore water pressure as the pressure exerted by a column of water in
an imaginary standpipe inserted at the given point.
The natural level of ground water is called the water table or the phreatic surface. Under
conditions of no seepage flow, the water table is horizontal. The magnitude of the pore water
pressure at the water table is zero. Below the water table, pore water pressures are positive.
The principle of effective stress was enunciated by Karl Terzaghi in the year 1936. This
1. At any point in a soil mass, the effective stress (represented by or 𝜎′) is related to
principle is valid only for saturated soils, and consists of two parts:
𝜎′ = 𝜎 − 𝑢-----------------------------
total stress (σ) and pore water pressure (u) as
Eq.6.3
Both the total stress and pore water pressure can be measured at any point.
2. All measurable effects of a change of stress, such as compression and a change of shearing
Compression = 𝑓1(𝜎)
resistance, are exclusively due to changes in effective stress.
56
Lecture Note on GTE-I
the capillary rise, and this depends on the grain size and the size of pores. In coarse soils, the
capillary rise is very small.
6.2 STRESSES WHEN NO FLOW TAKES PLACE THROUGH THE
SATURATED SOIL MASS
+
𝑧2)𝛾𝑤
The pressure uz is termed as the pore water pressure acting on the grains at depth z from the
weight of the submerged soil above N. If 𝛾𝑏 is the submerged unit weight of the soil, the total
surface of the sample. However, the total pressure at point N is due to the water head plus the
𝜎𝑧 = 𝑧𝛾𝑏 + (𝑧 + 𝑧2)𝛾𝑤-----------------
pressure at N is
Eq.6.4
The inter-granular or effective pressure at the point N is the difference between the total and the
𝜎′ = 𝜎 − 𝑢 =𝑧𝛾 ---------------
pore water pressures. Therefore, the effective pressure
𝑧 𝑧 𝑧 𝑏
Eq.6.5
The above equation clearly demonstrates that the effective pressure is independent of the depth
of water z2 above the submerged soil surface.
6.3 STRESSES WHEN FLOW TAKES PLACE THROUGH THE SOIL FROM TOP
TO BOTTOM
57
Lecture Note on GTE-I
𝑢𝑐 = (𝑧1 + 𝑧2 − ℎ)𝛾𝑤-------------------------
Therefore, the pore water pressure uc at the bottom level is
Eq.6.6
As per Fig.6.5(a), the pore water pressure at the bottom of container A when no flow takes place through
𝑢𝑐 = (𝑧1 + 𝑧2)𝛾𝑤-------------------
the soil sample is
Eq.6.7
𝛾𝑤when water flows through the soil sample from top to bottom. It may be understood that this
It is clear from the above expressions, that there is a decrease in pore water pressure to the extent of
decrease in pore2 water pressure is not due to velocity of the flowing water. The value of the
velocity head 𝑣 is a negligible quantity even if we take the highest velocity of flow that is
2𝑔
encountered in natural soil deposits. As shown in Fig. 6.5. (a), the total pressure 𝜎𝑐 at the bottom of the
𝜎′ = 𝜎 − 𝑢 =𝑧 + ℎ𝛾 ----------------
The effective pressure at the bottom of the container is
𝑐 𝑐 𝑐 1 𝑤
Eq.6.9
Equation indicates that in this case there is an increase in the effective pressure by ℎ𝛾𝑤 at the bottom of
the container A as compared to the earlier case. The effective pressure at the top surface of the
sample is zero as before. Therefore, the effective pressure at any depth z can be written as
58
Lecture Note on GTE-I
𝜎′ = 𝜎′ ( ) = (𝑧 𝛾 + ℎ𝛾 )( )----------------------
𝑧 𝑧
𝑧 1 𝑤
Eq.6.10
𝑐 𝑧 𝑏
=𝑧𝛾 +
𝑧1 ℎ𝑧𝛾𝑤
1
------------------ Eq.6.11
𝑧1
𝑏
ℎ𝑧𝛾𝑤
The above equation indicates that is the increase in the effective pressure as the water flows
𝑧1
from the surface to a depth z. This increase in effective pressure due to the flow of water through
the pores of the soil is known as seepage pressure. It may be noted that h is the total loss of head
as the water flows from the top surface of the sample to a depth z1.
The corresponding loss of head at depth z is (z/z1)h. Since (z/z1) = i the hydraulic gradient, the
𝜎 = 𝑧𝛾 + 𝑖𝑧𝛾 ---------------------
𝑧 𝑏 𝑤
Eq.6.12
The distribution of pore water and effective pressures are shown in Fig. 6.5.(b). In normal soil
deposits when flow takes place in the direction of gravity there will be an increase in the
effective pressure.
𝑢𝑐 = (𝑧1 + 𝑧2 + ℎ)𝛾𝑤------------------
As before the total pressure head 𝜎𝑐 at the bottom of the sample is
Eq.6.13
𝜎′ = 𝜎 − 𝑢 = 𝑧 𝛾 + (𝑧 + 𝑧 + ℎ)𝛾 − (𝑧 + 𝑧 + ℎ)𝛾
The effective pressure at the bottom of sample is, therefore,
𝑐 𝑐 𝑐 1 𝑏 1 𝑤 1 2 𝑤
2
= 𝑧1𝛾𝑏 − 𝑧𝛾𝑤------------- Eq.6.15
This equation indicates that there is a decrease in the effective pressure due to upward flow of
water. At any depth z, zγb is the pressure of the submerged soil acting downward and izγb is the
seepage pressure acting upward. The effective pressure reduces to zero when these two pressures
𝜎′ = 𝜎 − 𝑢 = 0 0r
balance. This happens when
𝑐 𝑐 𝑐
𝑧 𝛾 − 𝑧𝛾 =0 and hence 𝑖 = 𝑖 =𝛾𝑏-----------------
1 𝑏 𝑤
Eq.6.16
𝑐 𝛾𝑤
This indicates that the effective pressure reduces to zero when the hydraulic gradient attains a
maximum value which is equal to the ratio of the submerged unit weight of soil and the unit
weight of water. This gradient is known as the critical hydraulic gradient ic. In such cases,
cohesionless soils lose all of their shear strength and bearing capacity and a visible agitation of
soil grains is observed. This phenomenon is known as boiling or a quick sand condition.
Hence critical hydraulic gradient can be expressed as
𝑖𝑐 =
𝐺−1
1+𝑒
--------------- Eq.6.17
Quick conditions are common in excavations below the ground water table. This can be
prevented by lowering the ground water elevation by pumping before excavation. Quick
conditions occur most often in fine sands or silts and cannot occur in coarse soils. The larger the
particle size, the greater is the porosity. To maintain a critical gradient of unity, the velocity at
which water must be supplied at the point of inflow varies as the permeability. Therefore a quick
condition cannot occur in a coarse soil unless a large quantity of water can be supplied.
6.5 EFFECTIVE PRESSURE DUE TO CAPILLARY WATER RISE IN SOIL
Since the capillary pressure inside a soil mass is below atmospheric pressure, it draws the grains
of soils closer to each other at all points where the menisci touch the soil grains. Inter-granular
pressure of this type is called capillary pressure. The effective or inter-granular pressure at any
𝜎� ′ = 𝜎 − 𝑢𝑤-------------
point in a soil mass can be expressed by
60
Lecture Note on GTE-I
𝜎 = 𝜎 + 𝑢𝑤---------------
Fig.6.7: Effect of capillary pressure uc on soil vertical stress diagram
′
Eq.6.19
The increase in the inter-granular pressure due to capillary pressure acting on the grains leads to
greater strength of the soil mass.
At any point within the soil mass, the magnitudes of both total stress and pore water pressure are
dependent on the ground water position. With a shift in the water table due to seasonal
fluctuations, there is a resulting change in the distribution in pore water pressure with depth.
Changes in water level below ground result in changes in effective stresses below the water
table. A rise increases the pore water pressure at all elevations thus causing a decrease in
effective stress. In contrast, a fall in the water table produces an increase in the effective stress.
Changes in water level above ground do not cause changes in effective stresses in the ground
below. A rise above ground surface increases both the total stress and the pore water pressure by
the same amount, and consequently effective stress is not altered.
If both total stress and pore water pressure change by the same amount, the effective stress
remains constant.
Total and effective stresses must be distinguishable in all calculations. Ground movements and
instabilities can be caused by changes in total stress, such as caused by loading by foundations
and unloading due to excavations. They can also be caused by changes in pore water pressures,
such as failure of slopes after rainfall.
61
Lecture Note on GTE-I
Example 1: For the soil deposit shown below, draw the total stress, pore water pressure and
effective stress diagrams. The water table is at ground level.
Solution:
Total stress
At - 4m, = 1.92 x 4 = 7.68
At -11m, = 7.68 + 2.1 x 7 = 22.38
Pore water pressure
At - 4 m, u = 1 x 4 = 4
At -11 m, u = 1 x 11 = 11
Effective stress
At - 4 m , = 7.68 - 4 = 3.68
At -11m , = 22.38 - 11 = 11.38
62
Lecture Note on GTE-I
Example 2: An excavation was made in a clay stratum having = 2 T/m3. When the depth was
7.5 m, the bottom of the excavation cracked and the pit was filled by a mixture of sand and
water. The thickness of the clay layer was
10.5 m, and below it was a layer of pervious water-bearing sand. How much was the artesian
pressure in the
sand layer?
Solution:
When the depth of excavation was 7.5 m, at the interface of the CLAY and SAND layers, the
effective stress was equal to zero.
Downward pressure due to weight of clay = Upward pressure due to artesian pressure
(10.5 - 7.5) = , where h = artesian pressure head
3x2=1xh
h = 6 m = 0.6 kg/cm2 or 6 T/m2 artesian pressure
Example 3 A clay layer 3.66 m thick rests beneath a deposit of submerged sand 7.92 m thick.
The top of the sand is located 3.05 m below the surface of a lake. The saturated unit weight of
the sand is 19.62 kN/m3 and of the clay is 18.36 kN/m3 Compute (a) the total vertical pressure,
(b) the pore water pressure, and (c) the effective vertical pressure at mid height of the clay layer
(Refer to Fig. Ex.3.)
Sol:
(a) Total pressure
The total pressure cr, over the midpoint of the clay is due to the saturated weights of clay and
63
Lecture Note on GTE-I
sand layers plus the weight of water over the bed of sand, that is
σ=3.66/2 x 18.36 +7.92 x 19.62 + 3.05 x 9.81 = 33.6 + 155.4 + 29.9 = 218.9 kN/m2
(b) Pore water pressure is due to the total water column above the midpoint.
That is
u =3.66/2 x 9.81 + 7.92 x 9.81 + 3.05 x 9.81 = 125.6 kN/m2
(c) Effective vertical pressure
σ-u = 218.9-125.6 = 93.3 kN/m2
Example 4
The water table is located at a depth of 3.0 m below the ground surface in a deposit of sand 11.0
m thick (Fig. Ex. 4). The sand is saturated above the water table. The total unit weight of the
sand is 20 kN/m3. Calculate the (a) the total pressure, (b) the pore water pressure and (c) the
effective pressure at depths 0, 3.0, 7.0, and 11.0 m from the ground surface, and draw the
pressure distribution diagram.
The calculation of total and effective stresses at different depth are shown in the table
Problem No.1
A clay stratum 8.0 m thick is located at a depth of 6 m from the ground surface. The natural
moisture content of the clay is 56% and G = 2.75. The soil stratum between the ground surface
and the clay consists of fine sand. The water table is located at a depth of 2 m below the ground
64
Lecture Note on GTE-I
surface. The submerged unit weight of fine sand is 10.5 kN/m3, and its moist unit weight above
the water table is 18.68 kN/m3. Calculate the effective stress at the center of the clay layer.
The study of seepage of water through soil is important for the following engineering
problems:
65
Lecture Note on GTE-I
v = q/A = ki … Eq.7.2
Where q = discharge per unit time
A = total cross-sectional area of soil mass, perpendicular to the direction of flow
i = hydraulic gradient
k = Darcy’s coefficient of permeability
v = velocity of flow or average discharge velocity
It is defined as the average velocity of flow that will occur through the total cross-sectional are of
soil under unit hydraulic gradient. The coefficient of permeability is denoted as K. It is usually
expressed as cm/sec (or) m/day (or) feet/day.
When hydraulic gradient is unity, k is equal to V. Thus, the coefficient of permeability, or simply
permeability is defined as the average velocity of flow that will occur through the total cross-
sectional area of soil under unit hydraulic gradient. Dimensions are same as of velocity, cm/sec.
Figure 7.1 shows a soil sample of length L and cross-sectional area A. The sample is placed in a
cylindrical horizontal tube between screens. The tube is connected to two reservoirs R1 and R2 in
which the water levels are maintained constant. The difference in head between R{ and R2 is h.
This difference in head is responsible for the flow of water. Since Darcy's law assumes no
change in the volume of voids and the soil is saturated, the quantity of flow past sections AA, BB
and CC should remain the same for steady flow conditions. We may express the equation of
continuity as follows
66
Lecture Note on GTE-I
If soil mass is divided into two parts soil solid and void space, tjen the area are available for
passage of water is only through the area of voids A v. If vs is the velocity of flow in voids and v
is the average velocity across the section then we have
Or 𝑣� =
𝐴
𝑣---------------------
� 𝐴
Eq.7.5
= =
𝑣 𝐴 1
1+𝑒 𝐴𝑣 𝑛
But since, Eq.7.6
So 𝑣� = =
𝑣 1+𝑒
𝑒
-----------------
𝑛
𝑣--------------------
Since�1+𝑒𝑒is always greater than unity, 𝑣 is always greater than v.
Eq.7.7
𝑒 𝑠
7.1.3 Validity of Darcy’s Law:
In accordance with the Darcy’s Law, the velocity of flow through soil mass is directly proportion
to the hydraulic gradient for laminar flow condition only. It is expected that the flow to be
always laminar in case of fine-grained soil deposits because of low permeability and hence low
velocity of flow.
So v=ki
Where i is the hydraulic gradient and k is a constant known as coefficient of permeability.
However, in case of sands and gravels flow will be laminar up to a certain value of velocity for
each deposit and investigations have been carried out to find a limit for application of Darcy’s
law.
According to researchers, flow through sands will be laminar and Darcy’s law is valid so long as
≤ 1-------------
𝑣𝐷𝑎𝛾𝑤
Reynolds number expressed in the form is less than or equal to unity as shown below
𝜂𝑔
Eq.7.8
Where v = velocity of flow in cm/sec
Da = size of particles (average) in cm.
It is found that the limiting value of Reynolds number taken as 1 is very approximate as its actual
value can have wide variation depending partly on the characteristic size of particles used in the
equation.
7.2 FACTORS AFFECTING PERMEABILITY ARE
𝑘𝛼𝐷2
and particle size (D) can be shown from equation (1) as follows.
1
lower will be the permeability.
𝐾𝛼
𝑆𝑝𝑒𝑐𝑖𝑓𝑖𝑐 𝑆𝑢𝑟𝑓𝑎𝑐𝑒 𝐴𝑟𝑒𝑎
3. Shape of Soil Particle
Rounded Particles will have more permeability than angular shaped. It is due to specific surface
area of angular particles is more compared to rounded particles.
4. Void Ratio
In general, Permeability increases with void ratio. But it is not applicable to all types of soils. For
example, Clay has high void ratio than any other types of soil but permeability for clays is very
low. This is due to, the flow path through voids in case of clays is extremely small such that
water cannot permit through this path easily.
Clay 𝑘𝛼
𝐶𝑒3
The relation between coefficient of permeability and void ratio can be expressed as
For
1+𝑒
1+𝑒
5. Soil Structure
Structure of any two similar soil masses at same void ratio need not be same. It
varies according to the level of compaction applied. If a soil contains flocculated structure, the
particles are in random orientation and permeability is more in this case.
If the soil contains dispersed structure, the particles are in face to face orientation hence,
permeability is very low. The permeability of stratified soil deposits also varies according to the
flow direction. If the flow is parallel, permeability is more. If it is perpendicular, permeability is
less.
6. Degree of Saturation
Partially saturated soil contains air voids which are formed due to entrapped air or gas released
from the percolating fluid or water. This air will block the flow path thereby reduces the
permeability. Fully saturated soil is more permeable than partially saturated soil.
7. Water Properties
68
Lecture Note on GTE-I
Various properties of water or fluid such as unit weight and viscosity also effects the
permeability. However, unit weight of water will not affect much since it does not change much
with temperature.
But when temperature is increased viscosity decreases rapidly. From equation (1), permeability
𝛾
increase when viscosity decreases.
𝑘𝛼
𝑤
8. Temperature 𝜇
Temperature also affects the permeability in soils. The permeability is inversely proportional to
the viscosity of the fluid. It is known that viscosity varies inversely to the temperature. Hence,
Permeability is directly related to temperature.
Greater the temperature, higher will be the permeability. That is the reason; seepage is more in
1
summer seasons than in winter.
𝑘𝛼 𝛼 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑒
𝜇
9. Adsorbed Water
Adsorbed water is the water layer formed around the soil particle especially in the case of fine-
grained soils. This reduces the size of the void space by about 10%. Hence, permeability reduces.
10. Organic matter
Presence of organic matter decreases the permeability. This is due to blockage of voids by the
organic matter.
7.3 PERMEABILITY OF STRATIFIED SOIL DEPOSITS:
Natural soil deposits may exhibit stratification. Each layer may have its own coefficient of
permeability, assuming it to be homogeneous. The ‘average permeability’ of the entire deposit
will depend upon the direction of flow in relation to the orientation the bedding planes.
Two cases will be considered—the first one with flow perpendicular to the bedding planes and
the next with flow parallel to the bedding planes.
7.3.1 Flow Perpendicular to the Bedding Planes
69
Lecture Note on GTE-I
Let h1, h3, h3 ...hn be the thicknesses of each of the n layers which constitute the deposit, of total
thickness h. Let k1, k2, k3 ... kn be the Darcy coefficients of permeability of these layers
respectively.
In this case, the velocity of flow v, and hence the discharge q, is the same through all the layers,
for the continuity of flow.
∆ℎ1
The hydraulic gradients are
𝑖1 = ⁄ where i is the hydraulic gradient for deposit
ℎ1
∆ℎ2
𝑖2 = ⁄ -----------
ℎ2
∆ℎ𝑛
𝑖𝑛 = ⁄ --------------------
ℎ𝑛
Eq.7.10
Since q is the same in all the layers, and area of cross-section of flow is the same, the velocity is
the same in all layers.
This is the equation for average permeability for flow perpendicular to the bedding planes.
7.3.2 Flow Parallel to the Bedding Planes
Let the flow be parallel to the bedding planes as shown in Figure 7.3
With the same notation as in the first case, the hydraulic gradient i will be the same for all the
layers as for the entire deposit. Since v = ki, and k is different for different layers, v will be
different for the layers, say, v1, v2, ... vn
Also v1=k1i, v2=k2i and so on
Considering unit dimension perpendicular to the plane of the paper, the areas of flow for each
layer will be plane of the paper, the areas of flow for each layer will be h1, h2, ... hn respectively,
and it is h for the entire deposit.
The discharge through the entire deposit is equal to the sum of the discharge through the
individual layers. Assuming kx to be the average permeability of the entire deposit parallel to the
𝑞 = 𝑞1 + 𝑞2 + − − +𝑞𝑛------------------------
bedding planes, and applying the equation:
∴ = [
𝑘1ℎ1+𝑘2ℎ2+ −−+𝑘𝑛ℎ𝑛
]-----------------------------------
𝑘𝑥 ℎ
Eq.7.17
1 Permeameter mould, internal diameter = 100mm, effective height =127.3 mm, capacity =
1000ml.
2 Detachable collar, 100mm diameter, 60mm height
3 Dummy plate, 108 mm diameter, 12mm thick,
4 Drainage base, having porous disc
71
Lecture Note on GTE-I
5 Drainage cap having porous disc with a spring attached to the top.
6 Compaction equipment such as Proctor’s rammer or a static compaction equipment, as
specified in IS: 2720 (Part VII)-1965.
7 Constant head water supply reservoir
8 Vacuum pump
9 Constant head collecting chamber
10 Stop watch
11 Large funnel
12 Thermometer
13 Weighing balance accuracy 0.1g
14 Filter paper.
Procedure
Specimen Preparation
1. Remove the collar of the mould. Measure the internal dimensions of the mould. Weigh
the mould with dummy plate to the nearest gram.
2. Apply a little grease on the inside to the mould. Clamp the mould between the base plate
and the extension collar and place the assembly on a solid base.
3. Take about 2.5kg of the soil sample, from a thoroughly mixed wet soil, in the mould.
Compact the soil at the required dry density using a suitable compacting device.
4. Take a small specimen of the soil in a container for the water content determination.
5. Remove the collar and base plate. Trim the excess soil level with the top of the mould.
6. Clean the outside of the mould and the dummy plate. Find the mass of the soil in the
mould.
7. The mould with the sample is now placed over the permeameter. This will have drainage
and cap discs properly saturated
Test Procedure
1. Through the top inlet of the constant head reservoir, the specimen is connected.
2. The bottom outlet is opened and a steady flow is established
3. For a particular time interval, the quantity of flow can be collected.
4. Measure the difference of head (h) in levels between the constant head reservoir and the
outlet in the base.
5. For the same interval, this is repeated three times.
72
Lecture Note on GTE-I
𝑘 = 𝑄 × 𝐿/(𝐴 × 𝑡 × ℎ)-----------
Calculate the coefficient of permeability (k) using the following equation:
Eq.7.20
where: k = coefficient of permeability,
Q= quantity of water discharged,
L = distance between manometers,
A = cross-sectional area of specimen, t =
total time of discharge
h = head required for flow
7.5 FALLING HEAD PERMEABILITY TEST
Falling head permeability test is one of several techniques by which the permeability of soil is
determined. It is used to evaluate the permeability of fairly less previous soil particularly for fine
73
Lecture Note on GTE-I
grained soil sample. Permeability is the measure of the ability of soil to allow water to flow its
pores or voids.
A falling head permeameter is shown in Fig.7.5. The soil sample is kept in a vertical cylinder of
cross-sectional area A. A transparent stand pipe of cross sectional area, a, is attached to the test
cylinder. The test cylinder is kept in a container filled with water, the level of which is kept
constant by overflows. Before the commencement of the test the soil sample is saturated by
allowing the water to flow continuously through the sample from the stand pipe. After saturation
is complete, the stand pipe is filled with water up to a height of h0 and a stop watch is started. Let
the initial time be t0. The time t1 when the water level drops from h0 to h1 is noted. The hydraulic
conductivity k can be determined on the basis of the drop in head (h 0 – h1) and the elapsed time
(t1 – t0) required for the drop as explained below.
Let h be the head of water at any time t. Let the head drop by an amount dh in time dt.
𝑑𝑄 = 𝑘𝑖𝐴𝑑𝑡 = 𝑘 𝐴𝑑𝑡-------------
ℎ
The quantity of water flowing through the sample in time dt from Darcy's law is
𝐿
Eq.7.21
where, i = h/L the hydraulic gradient. The quantity of discharge dQ can be expressed as
dQ = -adh ---------------- Eq.7.22
Since the head decreases as time increases, dh is a negative quantity in Eq.7.22. & Eq. 7.21 can
= 𝑘𝐴𝑑𝑡/𝐿--------------
be equated to Eq. 7.22
−𝑎𝑑ℎ
ℎ
Eq.7.23
The discharge Q in time (t1 – t0) can be obtained by integrating Eq. 2 or 3. Therefore, Eq.4 can
𝒉𝟎
OR
𝑲 =
𝟐.𝟑𝒂𝑳
𝒍𝒐 𝟏 𝒉
𝒈
------------------- Eq.7.27
𝑨(𝒕𝟏−𝒕 𝟎 𝟏
𝟎)
Apparatus
1. Mould Assembly
The mould assembly including drainage base and drainage cap which need to be conform to IS:
11209-1985
2. Compaction Hammer
3. Set of Stand Pipes
Glass stand pipes varying in diameter from 5 to 20 mm, suitably mounted on stand or otherwise
fixed on wall
74
Lecture Note on GTE-I
75
Lecture Note on GTE-I
Suitable water reservoir capable or supplying water to the permeameter under constant head for constant
head test arrangement is required.
5. Miscellaneous Apparatus
For instance, IS sieves, mixing pan, graduated cylinder, meter scale, stop watch, 75micron wire gauge,
thermometer, and a source of de-aired water.
Take 2.5-kg soil from a thoroughly mixed air-dried or oven-dried material and evaluate
its moisture content.
Remove the collar of the mould. Measure the internal dimensions of the mould. Weigh
the mould with dummy plate to the nearest gram.
Apply a little grease on the inside to the mould.
Clamp the mould between the base plate and the extension collar and place the assembly
on a solid base.
Place soil specimen in the mould, and compact it at the required dry density using a
suitable compacting device.
Take a small specimen of the soil in a container for the water content determination.
Remove the collar and base plate. Trim the excess soil level with the top of the mould.
Clean the outside of the mould and the dummy plate.
Find the mass of the soil in the mould.
The mould with the sample is now placed over the permeameter. This will have drainage
and cap discs properly saturated
Trim the specimen in a form of cylinder not larger than 85cm in diameter, and having a
height equal to that of the mould.
Place the specimen over porous disc of the drainage base fixed to the mould.
Use impervious material like cement slurry to fill the space between mould and the
specimen.
Fix the drainage cap over the top of the mould.
Procedure
1. Connect the specimen to the selected stand-pipe through the top inlet.
2. Open the bottom outlet and record the time interval required for the water level to fall
from a known initial head to a known final head as measured above the center of the
outlet.
76
Lecture Note on GTE-I
3. Refill the stand-pipe with water and repeat the test till three successive observations give
nearly same time interval; the time intervals being recorded for the drop in head from the
same initial to final values, as in the first determination.
4. Alternatively, after selecting the suitable initial and final heads, h 1, and h2, respectively,
observe the time intervals for the head to fall from h1 to √ℎ1ℎ2 and similarly
from√ℎ1ℎ2 to h2.
5. The time intervals should be the same; otherwise the observation shall be repeated after
refilling the stand-pipe.
Data Sheet
The following values shall be recorded in the data sheet of variable head permeability test: Length of
specimen (L)
Diameter of specimen (D)
Volume of specimen (V)
Water content (ω)
Diameter of stand pipe (d)
Area of stand pipe (a)
Specific gravity of solids (Gs)
Solution
where , Q = 120 ml, t = 6 min, A = 35 cm2 , L = 20 cm, and h = 60 cm. Substituting, we have
78
Lecture Note on GTE-I
water levels in the test and observation wells remain stationary. When this condition is achieved
the water pumped out of the well is equal to the inflow into the well from the surrounding strata.
The water levels in the observation wells and the rate of water pumped out of the well would
provide the necessary additional data for the determination of k. As the water from the test well
is pumped out, a steady state will be attained when the water pumped out will be equal to the
inflow into the well. At this stage the depth of water in the well will remain constant. The
drawdown resulting due to pumping is called the cone of depression. The maximum drawdown
DQ is in the test well. It decreases with the increase in the distance from the test well. The
depression dies out gradually and forms theoretically, a circle around the test well called the
circle of influence. The radius of this circle is called the radius of influence of the depression
cone.
7.6.1 Equation for k for an Unconfined Aquifer
Figure below gives the arrangement of test and observation wells for an unconfined aquifer.
Only two observation wells at radial distances of r1 and r2 from the test well are shown. When
the inflow of water into the test well is steady, the depths of water in these observation wells are
h1 and h2 respectively.
Let h be the depth of water at radial distance r. The area of the vertical cylindrical surface of
𝐴
The= 2𝜋𝑟ℎ-----------------Eq.7.29
radius r and depth h through which water flows is given as
hydraulic gradient is 𝑖 =
𝑑ℎ
𝑑𝑟
As per Darcy's law the rate of inflow into the well when the water levels in the wells remain
stationary is
q = kiA-------- Eq.7.30
𝑞 =𝑘 2𝜋𝑟ℎ----------------
𝑑ℎ
Substituting for A and i the rate of inflow across the cylindrical surface is
𝑑𝑟
Eq.7.31
=
2𝜋𝑘ℎ𝑑ℎ
Rearranging the terms, we have
𝑑𝑟
𝑟 𝑞
------------------ Eq.7.32
79
Lecture Note on GTE-I
𝑘 = 2 2 𝑙𝑜𝑔 ------------
The equation for k2.3𝑞 𝑟2
after integration and rearranging is
𝜋(ℎ −ℎ )
Eq.7.34
2 1 10 𝑟1
Proceeding in the same way as before another equation for k in terms of r0, h0 and Ri can be
𝑘 = 𝑙𝑜𝑔 𝑖----------
2.3𝑞 to Fig.7.6)
established as (referring 𝑅
Eq.7.35
𝜋(𝐻2−ℎ2
0) 10 𝑟
0
If we write h0 = (H- D0) in the above Eq.7.35 , where D0 is the depth of maximum drawdown in
𝑙𝑜𝑔 𝑖--------------
𝑅
the test well, we have
𝑘 =
Eq.7.36
2.3𝑞
𝜋𝐷0(2𝐻− 10 𝑟
𝐷0) 0
Now the maximum yield from the well may be written from the above expression as
1
𝑞=
𝜋𝐷0(2𝐻−𝐷0)𝑘 𝑙𝑜
𝑔
𝑅 -------------
( 𝑖)
Eq.7.37
2.3
10 𝑟0
7.6.2 Equation for k in a Confined Aquifer
Figure 7.7 shows a confined aquifer with the test and observation wells. The water in the
observation wells rises above the top of the aquifer due to artesian pressure. When pumping from
such an artesian well two cases might arise. They are:
Case 1. The water level in the test well might remain above the roof level of the aquifer at steady
flow condition.
80
Lecture Note on GTE-I
81
Lecture Note on GTE-I
𝐴 = 2𝜋𝑟𝐻0-------------
area is
𝑑ℎ 𝑑𝑟
𝑞 = 𝑘𝑖𝐴 = 𝑘 2𝜋𝑟𝐻 --------------
𝑑𝑟 0
Eq.7.39
ℎ2
integration of𝑟 the equation after rearranging the terms yields
𝑑ℎ =
𝑞 2 𝑑𝑟
∫
The
∫
ℎ1
-------------------- Eq.7.40
2𝜋𝑘𝐻0 𝑟1 𝑟
− ) = 2𝜋𝑘 𝑙𝑛 -------------------
𝑞 𝑟2
(ℎ2 ℎ1
Or Eq.7.41
𝐻0 𝑟1
𝑙𝑜𝑔 -------------
𝑟2
The equation for k is
𝑘 =
Eq.7.42
2.3𝑞
2𝜋𝐻0(ℎ2−ℎ 10 𝑟
1)
1
Alternate Equations
𝑙𝑜𝑔 ----------------
𝑟1
As before we can write the following equation for determining k
𝑘 =
Eq.7.43
2.3𝑞
2𝜋𝐻0(ℎ1−ℎ 10 𝑟
0)
0
𝑙𝑜𝑔
𝑅𝑖
𝑘 =
-------------------- Eq.7.44
2.3𝑞
2𝜋𝐻0(𝐻−ℎ 10 𝑟
0)
0
𝑘 = 𝑙𝑜𝑔
2.3𝑞 𝑅𝑖
----------------------- Eq.7.45
2𝜋𝐻0 10 𝑟
0
𝐷0
Case 2. When h0 < H0
Under the condition when h0 is less than H0, the flow pattern close to the well is similar to that of
an unconfined aquifer whereas at distances farther from the well the flow is artesian. Muskat
𝑙𝑜𝑔 𝑖----------------
𝑅
(1946) developed an equation to determine the hydraulic conductivity. The equation is
𝑘 =2𝜋(𝐻𝐻
2.3𝑞
0−𝐻 −
2
Eq.7.46
10 𝑟
2
ℎ ) 0
0 0
Example No.1
A pumping test was carried out for determining the hydraulic conductivity of soil in place. A
well of diameter 40 cm was drilled down to an impermeable stratum. The depth of water above
the bearing stratum was 8 m. The yield from the well was 4 m3/ min at a steady drawdown of 4.5
82
Lecture Note on GTE-I
m. Determine the hydraulic conductivity of the soil in m/day if the observed radius of influence
was 150m.
that 𝑘 = 𝑙𝑜 𝑅𝑖
Solution
𝑔
We know
2.3𝑞
𝜋𝐷0(2𝐻−𝐷0) 10 𝑟
0
Here q = 4 m /min = 4 x 60 x 24 m3/day
3
83
Lecture Note on GTE-I
Example No.2
A field pumping test was conducted from an aquifer of sandy soil of 4 m thickness confined
between two impervious strata. When equilibrium was established, 90 liters of water was
pumped out per hour.
The water elevation in an observation well 3.0 m away from the test well was 2.1 m and another
6.0 m away was 2.7 m from the roof level of the impervious stratum of the aquifer. Find the
value of k of the soil in m/sec.
Solution:
We know that k for confined aquifer is given as
2.3𝑞 𝑟2
𝑘 =
2𝜋𝐻 (ℎ 𝑙𝑜𝑔10 1
−ℎ
0 2
) 𝑟
1
q = 90 x 103 cm3/hr = 25 x 10-6 m3/sec
H0=4 m, r1=3 m and r2=6m
Substituting the above values we get k= 1.148 x 10-6 m/s
Problem No.1
Calculate the yield per hour from a well driven into a confined aquifer. The following data are
Available, height of original piezometric level from the bed of the aquifer, H = 30m thickness of
aquifer, H0 4m, the depth of water in the well at steady state, h0 = 6 m hydraulic conductivity of
soil = 0.079 m/min, radius of well, r0 = 0.25 m, radius of influence, Ri = 150 m
this stage flow should maintain gravity flow under constant head. Water may also be allowed to
entre the pipe under some pressure head.
𝑘=
𝑞
The permeability is calculated from the expression as
5.5𝑟 --------------
ℎ
Eq.7.47
Where q= constant rate of flow
h= Difference head of flow (gravity head and pressure head)
r = radius if the casing
Packer Test
Packers are primarily used in bore holes for testing the permeability of rocks under applied
pressures. The apparatus used for the pressure test is comprised of a water pump, a manually
adjusted automatic pressure relief valve, pressure gage, a water meter and a packer assembly.
The packer assembly consists of a system of piping to which two expandable cylindrical rubber
sleeves, called packers, are attached. The packers which are provided by means of sealing a
limited section of bore hole for testing should have a length five times the diameter of the hole.
They may be of the pneumatically or mechanically expandable type. The former are preferred
since they adapt to an oversized hole whereas the latter may not. However, when pneumatic
packers are used, the test apparatus must also include an air or water supply connected, through a
pressure gage, to the packers by means of a higher pressure hose. The piping of a packer
assembly is designed to permit testing of either the portion of the hole between the packers or the
portion below the lower packer. The packers are usually set 50, 150 or 300 cm apart. The wider
spacing is used for rock which is more uniform. The short spacing is used to test individual joints
which may be the cause of high water loss in otherwise tight strata.
Two types of packer methods are used for testing of permeability. They are:
1. Single packer method.
2. Double packer method.
The single packer method is useful where the full length of the borehole cannot stand
uncased/un-grouted in soft rocks, such as soft sand stone, clay shale or due to the highly
fractured and sheared nature of the rocks, or where it is considered necessary to have
permeability values side by side with drilling. Where the rocks are sound and the full length of
the hole can stand without casing/grouting, the double packer method may be adopted. The
disadvantage of the double packer method is that leakage through the lower packer can go
unnoticed and lead to overestimation of water loss.
Single Packer Method
The method used for performing water percolation tests in a section of a drilled hole using a
single packer is shown in Fig. 3 a. In this method the hole should be drilled to a particular depth
desirable for the test. The core barrel should then be removed and the hole cleaned with water.
The packer should be fixed at the desired level above the bottom of the hole and the test
performed. Water should be pumped into the section under pressure. Each pressure should be
maintained until the readings of water intake at intervals of 5 min show a nearly constant reading
of water intake for one particular pressure. The constant rate of water intake should be noted.
85
Lecture Note on GTE-I
After performing the test the entire assembly should be removed. The drilling should then
proceed for the next test section.
𝑘=
𝑞
Reclamation, 1968)
2𝜋𝐿 𝑙𝑛
𝐿 for 𝐿 ≥
𝐻
10𝑟0
and------------------ Eq.7.48
𝑟
0
𝑘=
𝑞
2𝜋𝐿 𝑠𝑖𝑛ℎ > 𝐿 ≥ 𝑟0-------------------
−1
𝐻 𝐿
10𝑟0
for Eq.7.49
where 2𝑟0
k = hydraulic conductivity
q = constant rate of flow into the hole
L = length of the test section
H = differential head on the test section
Example: A sample in a variable head permeameter is 8 cm in diameter and 10 cm high. The
permeability of the sample is estimated to be 10x10 -4 cm/sec. If it is desired that the head in the
stand pipe should fall from 24 cm to 12 cm in 3min, determine the size the stand pipe which
should be used.
Sol:
86
Lecture Note on GTE-I
87
Lecture Note on GTE-I
𝐴𝑡 10 ℎ
2
𝜋 × 16 × 180
Substituting for various values, we get
𝑎 =
2
2.303 × 104(𝑙𝑜𝑔
1
0
2)
Hence a=1.305 cm
Diameter of stand pipe d= 1.29 cm
Solution:
2.3𝑞
𝑘 = 𝑙𝑜𝑔10𝑟2
𝜋(ℎ2 − ℎ2) 𝑟
2 70 1 1
= 𝑙𝑜𝑔 = 1.18 × 10−1mm/sec
2.3×2×100
𝜋(31.31−292)×60 10 15
88
Lecture Note on GTE-I
Consider a soil element of infinitesimally small size of dx and dz in X- and Z-directions, respectively,
through which the flow is taking place, shown in Fig. 8.1. Consider unit length of the soil
element in the Y-direction.
∑ 𝑞𝑖 = ∑ 𝑞𝑜-----------------
element is equal to the quantity of water leaving the soil element –
𝑥 𝑥 𝑥 𝑧 𝑧
𝛛𝑥 𝑧 𝛛𝑧
Eq.8.6
89
Lecture Note on GTE-I
𝑥 𝑥 𝑥 𝑧 𝑧
𝛛𝑥 𝑧 𝛛𝑧
(𝑣 ) + (𝑣 ) = 0------------------
𝛛
On simplifying, and dividing throughout with dz. dx, we have
𝛛
𝑥 𝑧
Eq.8.9
𝛛𝑥 𝛛𝑧
Let h be the total head at any point. Then the component of hydraulic gradient in the X-direction
will be
ix = – ∂h/∂x-------------- Eq.8.10
The negative sign is to indicate that the head decreases in the direction of flow. Similarly, the
component of hydraulic gradient in the Z-direction will be –
iz = -∂h/∂z------------------ Eq.8.11
By Darcy’s law, we know that –
vx = kx .ix = – kx . ∂h/∂x------------- Eq.8.12
Similarly
vz = kz .iz = – kz . ∂h/∂x-------------- Eq.8.13
𝑥 𝑧 𝛛𝑧
Eq.8.14
𝛛𝑥 𝛛
𝛛𝑥 𝑧
Assuming that the soil mass is homogeneous; permeability will be same throughout in a given
−𝑘 ( )−𝑘 ( ) = 0------------
𝛛 𝛛ℎ 𝛛 𝛛ℎ
direction. Hence
𝑥 𝛛𝑥 𝑧 𝛛𝑧
Eq.8.15
𝛛𝑧
𝛛𝑥
≡ = 0-----------------
𝛛2ℎ
+
𝛛2ℎ
𝑥
𝑘
𝑧
2 𝑘
𝛛𝑧
Eq.8.16
𝛛𝑥
+ = 0----------------
If2 the soil
𝛛 ℎ 𝛛2ℎ
is isotropic, then kx = kz = k. Then, we have
𝛛𝑥2
Eq.8.17
𝛛𝑧2
This is the Laplace equation for two-dimensional flow. It says that the change of gradient in the
x--direction plus the change of gradient in the z-direction is zero. The solution of this equation
gives a family of curves meeting at right angles to each other. One family of these curves
represents flow lines and the other equipotential lines. The graphical representation of solution to
the Laplace equation which gives two sets of mutually perpendicular curves is known as Flow
net.
8.1 Velocity Potential
Velocity potential is a scalar function of space and time such that its derivative in any direction
gives the component of velocity in that direction. Thus –
ɸ = kh ------------- Eq.8.18
∂ɸ = ∂x = vx ----------------- Eq.8.19
∂ɸ = ∂z = vz-------------------- Eq.8.20
90
Lecture Note on GTE-I
( )+ ( ) = 0----------------------------
𝛛∅ 𝛛 𝛛∅
Substituting these in Laplace Equation, we have
𝛛
𝛛 𝛛 𝛛 𝛛𝑧
Eq.8.21
𝑥 𝑥 𝑧
91
Lecture Note on GTE-I
+ = 0------------------
or2
𝛛 ∅ 𝛛2∅
𝛛𝑥2
Eq.8.22
𝛛𝑧2
This is the Laplace equation in terms of velocity potential.
8.2 Stream Function
Stream function is a scalar function of space and time such that its derivative in any direction
gives the component of velocity in the perpendicular direction clockwise. Thus a positive stream
function when derived would give negative velocity and vice versa, as shown in Fig. 8.2.
Referring to Fig. 2, we have
𝛛 ℎ + = 0----------------
𝛛2ℎ
where anisotropy exists. This equation may be written in the form
2
𝑘𝑧
𝛛𝑥
Eq.8.25
𝛛𝑧2
2
𝑘𝑥
If we consider a new coordinate variable xc measured in the same direction as x multiplied by a
constant, expressed by
𝑥
= √ -------------------
𝑘𝑧
𝑥𝑥
Eq.8.26
𝑐 𝑘𝑥
+ = 0----------------
𝛛2ℎ 𝛛2ℎ
Eq.8.27
𝛛𝑥�
2
𝛛𝑧2
�
Now this equation is a Laplace equation in the coordinates xc and z. This equation indicates that
a cross-section through an anisotropic soil can be transformed to an imaginary section which
possesses the same permeability in all directions. The transformation of the section can be
𝑘𝑧
effected by multiplying the x-coordinates by √ and keeping the z-coordinates at the natural
𝑘𝑥
scale. The flow net can be sketched on this transformed section. The permeability to be used
with the transformed section is
𝑘𝑒 = √𝑘𝑧𝑘𝑥----------------- Eq.8.28
From the transformed section, the rate of seepage can be determined using Eq. 7.1 with
exception that ke is to be substituted for k (see Fig. 8.3):
93
Lecture Note on GTE-I
94
Lecture Note on GTE-I
equipotential line k is ∆ℎ which the difference between the heads is shown by the piezometers.
This loss of head ∆ℎ is a fraction of the total head lost.
8.5 FLOW NET CONSTRUCTION
Flow nets are constructed in such a way as to keep the ratio of the sides of each block bounded
by two flow lines and two equipotential lines a constant. If all the sides of one such block are
equal, then the flow net must consist of squares. The square block referred to here does not
constitute a square according to the strict meaning of the word; it only means that the average
width of the square blocks is equal. For example, in Fig.8.4, the width a1 of block 1 is equal to its
length b1. The area bounded by any two neighboring flow lines is called a flow channel. If the
flow net is constructed in such a way that the ratio a/b remains the same for all blocks, then it can
be shown that there is the same quantity of seepage in each flow channel. In order to show this
consider two blocks 1 and 2 in one flow channel and another block 3 in another flow channel as
shown in Figure. Block 3 is chosen in such a way that it lies within the same equipotential lines
that bound the block 2. Darcy's law for the discharge through any block such as 1 per unit length
𝑏 𝑏
Eq.8.29
where ∆ℎ represents the head loss in crossing the block. The expressions in this form for each of
= =1
𝑎1 𝑎2
the blocks are all squares then
𝑏1 𝑏2
Since blocks 1 and 2 are in the same flow channel, we have∆𝑞1 = ∆𝑞2. Since blocks 2 and 3
are within the same equipotential lines we have ∆ℎ1 = ∆ℎ2 If these equations are inserted we
This proves that the same quantity flows through each block and there is the same head drop in
crossing each block if all the blocks are squares or possess the same ratio a/b. Flow nets are
constructed by keeping the ratio a/b the same in all figures. Square flow nets are generally used
in practice as this is easier to construct.
There are many methods that are in use for the construction of flow nets. Some of the important
methods are
1. Analytical method,
2. Electrical analog method,
3. Scaled model method,
4. Graphical method.
The analytical method, based on the Laplace equation although rigorously precise, is not
universally applicable in all cases because of the complexity of the problem involved. The
95
Lecture Note on GTE-I
mathematics involved even in some elementary cases is beyond the comprehension of many
design engineers. Although this approach is sometimes useful in the checking of other methods,
it is largely of academic interest.
The electrical analogy method has been extensively made use of in many important for design
problems. However, in most of the cases in the field of soil mechanics where the estimation of
seepage flows and pressures are generally required, a more simple method such as the graphical
method is preferred.
Scaled models are very useful to solve seepage flow problems. Soil models can be constructed to
depict flow of water below concrete dams or through earth dams. These models are very useful
to demonstrate the fundamentals of fluid flow, but their use in other respects is limited because
of the large amount of time and effort required to construct such models.
The graphical method developed by Forchheimer (1930) has been found to be very useful in
solving complicated flow problems. A. Casagrande (1937) improved this method by
incorporating many suggestions. The main drawback of this method is that a good deal of
practice and aptitude are essential to produce a satisfactory flow net. In spite of these drawbacks,
the graphical method is quite popular among engineers.
8.5.1 Graphical Method
The usual procedure for obtaining flow nets is a graphical, trial sketching method, sometimes
called the Forchheimer Solution. This method of obtaining flow nets is the quickest and the most
practical of all the available methods. A. Casagrande (1937) has offered many suggestions to the
beginner who is interested in flow net construction. Some of his suggestions are summarized
below:
1. Study carefully the flow net pattern of well-constructed flow nets.
2. Try to reproduce the same flow nets without seeing them.
3. As a first trial, use not more than four to five flow channels. Too many flow channels would
confuse the issue.
4. Follow the principle of 'whole to part', i.e., one has to watch the appearance of the entire flow
net and when once the whole net is found approximately correct, finishing touches can be given
to the details.
5. All flow and equipotential lines should be smooth and there should not be any sharp
transitions between straight and curved lines.
The above suggestions, though quite useful for drawing flow nets, are not sufficient for a
beginner. In order to overcome this problem, Taylor (1948) proposed a procedure known as the
procedure by explicit trials. Some of the salient features of this procedure are given below:
1. As a first step in the explicit trial method, one trial flow line or one trial equipotential
line is sketched adjacent to a boundary flow line or boundary equipotential.
2. After choosing the first trial line (say it is a flow line), the flow path between the line and the
boundary flow line is divided into a number of squares by drawing equipotential lines. These
96
Lecture Note on GTE-I
equipotential lines are extended to meet the bottom flow line at right angles keeping in view that the
lines drawn should be smooth without any abrupt transitions.
3. The remaining flow lines are next drawn, adhering rigorously to square figures.
4. If the first trial is chosen property, the net drawn satisfies all the necessary conditions.
Otherwise, the last drawn flow line will cross the bottom boundary flow line, indicating that the
trial line chosen is incorrect and needs modification.
5. In such a case, a second trial line should be chosen and the procedure repeated.
A typical example of a flow net under a sheet pile wall is given in Fig. 1. It should be
understood that the number of flow channels will be an integer only by chance. That means, the
bottom flow line sketched might not produce full squares with the bottom boundary flow line. In
such a case the bottom flow channel will be a fraction of a full flow channel. It should also be
noted that the figure formed by the first sketched flow line with the last equipotential line in the
region is of irregular form. This figure is called a singular square. The basic requirement for such
squares, as for all the other squares, is that continuous sub-division of the figures give an
approach to true squares. Such singular squares are formed at the tips of sheet pile walls also.
Squares must be thought of as valid only where the Laplace equation applies. The Laplace
equation applies to soils which are homogeneous and isotropic. When the soil is anisotropic, the
flow net should be sketched as before on the transformed section. The transformed section can be
obtained from the natural section explained earlier.
8.6 DETERMINATION OF QUANTITY OF SEEPAGE
Flow nets are useful for determining the quantity of seepage through a section. The quantity of seepage
∆𝑞 = 𝑘∆ℎ--------------
q is calculated per unit length of the section. The flow through any square can be written as
Eq.8.31
Let the number of flow channel and equipotential drops in a section be Nf and Nd, respectively. Since all
∆ℎ = -------------------
ℎ
drops are equal, we can write
Eq.8.32
𝑁𝑑
𝑞 = 𝑁𝑓∆𝑞----------------------
Since the discharge in each flow channel is the same we can write,
𝑁𝑓
𝑞 = 𝑘ℎ ----------------- Eq.8.34
𝑁𝑑
This Eq. can also be used to compute the seepage through anisotropic sections by writing ke in place of
k, where ke is equal to √𝑘𝑥𝑘𝑧 , where kx and kz are the hydraulic conductivities in the x and z
directions, respectively. The validity of this relationship can be proved as follows.
Consider a figure bounded by flow and equipotential lines in which the flow is parallel to the x
direction. In Figure 8.5 the figure in question is drawn to a transformed scale in (b) and the same to the
natural scale in (a). In Fig. 8.5 (b) the permeability has the effective value ke in both the x and z
directions and the flow through the square according to
∆𝑞 = ∆ℎ𝑎
= ∆ℎ---
𝑘𝑒 𝑘
Eq.8.35
� 𝑎
�
97
Lecture Note on GTE-I
𝑘𝑥 𝑘
----------------- Eq.8.36
𝑘𝑥
𝑎×√
𝑘𝑧
8.7 DETERMINATION OF SEEPAGE PRESSURE
Flow nets are useful in the determination of the seepage pressure at any point along the flow
path.
Consider the cubical element 1 in Fig.8.4 (a) with all the sides equal to a. Let h1 be the
piezometric head acting on the face kt and h2 on face jo.
𝑃1 − 𝑃2 = 𝑎2𝛾𝑤(ℎ1 − ℎ2)------------
The differential force acting on the element is
∆ℎ , we can write
Eq.8.39
PS = iγw-----------------
3
Where, a is the volume of the element. The force per unit volume of the element is, therefore
Eq.8.41
This force exerts a drag on the element known as the seepage pressure. It has the dimension of
unit weight, and at any point its line of action is tangent to the flow line. The seepage pressure is
a very important factor in the stability analysis of earth slopes. If the line of action of the seepage
force acts in the vertical direction upward as on an element adjacent to point x, in Fig.8.6 (a), the
force that is acting downward to keep the element stable is the buoyant unit weight of the
element. When these two forces balance, the soil will just be at the point of being lifted up, and
there will be effectively no grain-to-grain pressures.
98
Lecture Note on GTE-I
99
Lecture Note on GTE-I
100
Lecture Note on GTE-I
Solution:
Fig. gives the flow net and other details. The differential head between the bottom of trench and
the water level in the river is 14.8 ft.
Number of channels = 6
Number of equipotential drops =10
𝑁𝑓
Discharge
𝑞 = 𝑘ℎ
�𝑑
-4
From which, k= 1x10 ft/sec �
The distance between the last two equipotentials given is 2.95 ft. The calculated hydraulic
i = =0.50
∆h
gradient is
∆s
8.9 SEEPAGE FLOW THROUGH HOMOGENEOUS EARTH DAMS
In almost all problems concerning seepage beneath a sheet pile wall or through the foundation of
a concrete dam all boundary conditions are known. However, in the case of seepage through an
earth dam the upper boundary or the uppermost flow line is not known. This upper boundary is a
free water surface and will be referred to as the line of seepage or phreatic line. The seepage line
may therefore be defined as the line above which there is no hydrostatic pressure and below
which there is hydrostatic pressure. In the design of all earth dams, the following factors are very
important.
1. The seepage line should not cut the downstream slope.
2. The seepage loss through the dam should be the minimum possible.
The two important problems that are required to be studied in the design of earth dams are:
1. The prediction of the position of the line of seepage in the cross-section.
2. The computation of the seepage loss.
If the line of seepage is allowed to intersect the downstream face much above the toe, more or
less serious sloughing may take place and ultimate failure may result. This mishap can be
prevented by providing suitable drainage arrangements on the downstream side of the dam. The
section of an earth dam may be homogeneous or non-homogeneous. A homogeneous dam
contains the same material over the whole section and only one coefficient of permeability may
be assumed to hold for the entire section. In the non homogeneous or the composite section, two
or more permeability coefficients may have to be used according to the materials used in the
section. When a number of soils of different permeabilities occur in a cross-section, the
prediction of the position of the line of seepage and the computation of the seepage loss become
quite complicated.
101
Lecture Note on GTE-I
Fig.8.8: Basic parabola and the phreatic line for a homogeneous earth dam
It has been noticed from experiments on homogeneous earth dam models that the line of seepage
assumes more or less the shape of a parabola as illustrated in Fig. 8.8. In some sections a little
divergence from a regular parabola is required at the surfaces of entry and discharge of the line
of seepage. In some ideal sections where conditions are favorable the entire seepage line may be
considered as a parabola. When the entire seepage line is a parabola, all the other flow lines will
be confocal parabolas. The equipotential lines for this ideal case will be conjugate confocal
parabolas as shown in Fig.1. As a first step it is necessary to study the ideal case where the entire
flow net consists of conjugate confocal parabolas.
8.9.1 Flow net consisting of conjugate confocal parabolas
As a prelude to the study of an ideal flow net comprising of parabolas as flow and equipotential
line, it is necessary to understand the properties of a single parabola. The parabola ACV
illustrated in Fig.8.9, is defined as the curve whose every point is equidistant from a point F
called the focus and a line DG called the directrix. If we consider any point, say, A, on the curve,
we can write FA = AG, where the line AG is normal to the directrix. If F is the origin of
coordinates, and the coordinates of point A are (x, y), we can write
𝐴𝐹 = 𝑦√𝑥2 2+ 𝑦2 = 𝐴𝐺 = 𝑥 + 𝑦0------ Eq.8.44
2
−𝑦
Or 𝑥 = 0
----------------------------- Eq.8.45
2𝑦0
where, y0 = FD
Eq.8.45 is the equation of the basic parabola. If the parabola intersects the y-axis at C, we can
write
FC=CE = y0
Similarly for the vertex point V, the focal distance a0 is
𝐹𝑉 = 𝑉𝐷 𝑦0
= 𝑎0 =2
Figure 8.9 illustrates the ideal flow net consisting of conjugate confocal parabolas. All the
parabolas have a common focus F.
The boundary lines of such an ideal flow net are:
1 . The upstream face AB, an equipotential line, is a parabola.
102
Lecture Note on GTE-I
𝑎0
Or
2
Eq. 8.48
=
𝑑𝑦 ------------------------------ Eq.8.50
𝑦0
𝑑
𝑥 √2𝑥𝑦0+𝑦2
0
= = 1-----------------------------
𝑦0
For the point C which has coordinates (0, y0), the hydraulic gradient from Eq. (8.50) is
𝑑𝑦
𝑑
Eq.8.51
𝑥 √𝑦02
𝑞=𝑘 𝑦 = 𝑘𝑦 ------------------------
𝑑𝑦
Therefore, the seepage quantity across section FC is
Eq.8.52
𝑑𝑥 0 0
8.9.3 Seepage Through Homogeneous and Isotropic Earth Dams
Types of Entry and Exit of Seepage lines
The flow net consisting of conjugate confocal parabolas is an ideal case which is not generally
met in practice. Though the top flow line resembles a parabola for most of its length, the
departure from the basic parabola takes place at the faces of entry and discharge of the flow line.
The departure from the basic parabola depends upon the conditions prevailing at the points of
entrance and discharge of the flow line as illustrated in Fig. 8.10 from (a) to (e).
upstream face AB with the horizontal is less than 90°. It can be assumed in this case the coarse
material used to support the face AB is highly permeable and does not offer any resistance for
105
Lecture Note on GTE-I
flow. In such cases AB taken as the upstream equipotential line. The top flow line cannot therefore be
normal to the equipotential line. However, this line possesses zero gradient and velocity at the
point of entry. This zero condition relieves the apparent inconsistency of deviation from a normal
intersection. The conditions prevailing at the downstream toe of the dam affect the type of exit of
the flow line at the discharge face. In Fig. 8.10 (c) the material at the toe is the same as in the
other parts of the dam whereas in (d) and (e) rock toe drains are provided. This variation in the
soil condition at the toe affects the exit pattern of the flow line. The flow line will meet the
discharge face FE tangentially in 8.10 (c). This has to be so because the particles of water as they
emerge from the pores at the discharge face have to conform as nearly as possible to the direction
of gravity. But in cases where rock toe drains are provided, the top flow line becomes tangential
to the vertical line drawn at the point of exit on the discharge face as shown in (d) and (e) of Fig.
8.10.
8.10 METHOD OF LOCATING SEEPAGE LINE
The general method of locating the seepage line in any homogeneous dam resting on an impervious
foundation may be explained with reference to Fig. 8.11 (a).
As explained earlier, the focus F of the basic parabola is taken as the intersection point of the bottom
flow line BF and the discharge face EF. In this case the focus coincides with the toe of the dam.
One more point is required to construct the basic parabola. Analysis of the location of seepage
lines by A. Casagrande has revealed that the basic parabola with focus F intersects the upstream
water surface at A such that AA'= 0.3m, where m is the projected length of the upstream
equipotential line A'B on the water surface. Point A is called the corrected entrance point.
106
Lecture Note on GTE-I
The parabola APSV may now be constructed as per Eq. (2). The divergence of the seepage line
from the basic parabola is shown as AT1 and SD in Fig. 8.11 (a). For dams with flat slopes, the
divergences may be sketched by eye keeping in view the boundary requirements. The error
involved in sketching by eye, the divergence on the downstream side, might be considerable if
the slopes are steeper. Procedures have therefore been developed to sketch the downstream
divergence as explained below. As shown in Fig. 8.11 (a), E is the point at which the basic
DF as a. The values of ∆𝑎 and a + ∆a, vary with the angle, β, made by the discharge face with
parabola intersects the discharge face. Let the distance ED be designated as Aa and the distance
the horizontal measured clockwise. The angle may vary from 30° to 180°. The discharge face is
horizontal as shown in Fig. 8.10(e). Casagrande (1937) determined the ratios of ∆𝑎 for a
(𝑎+∆𝑎)
number of discharge slopes varying from 30° to 180° and the relationship is shown in a graphical
with F as the focus. With the known (a +∆𝑎) and the discharge face angle β, ∆𝑎 can be
form Fig. 8.11 (b). The distance (a +∆a) can be determined by constructing the basic parabola
determined from Fig. 8.11 (b). The point D may therefore be marked out at a distance of ∆𝑎
from
E. With the point D known, the divergence DS may be sketched by eye.
It should be noted that the discharge length a, is neither an equipotential nor a flow line, since it
is at atmospheric pressure. It is a boundary along which the head at any point is equal to the
elevation.
8.10.1 Analytical Solutions for Determining a and q
𝑎
Casagrande (1937) proposed the following equation for determining ‘a’ for β < 30°
=
𝑑
− √
𝑑2
𝑐𝑜𝑠2 − 𝑠𝑖𝑛2
𝑐𝑜𝑠𝛽 ℎ2 ---------------------- Eq. 8.53
𝛽 𝛽
L. Casagrande (1932) gave the following equation for ‘a’ when β lies between 30° and 90°.
𝑎 = √ℎ2 + 𝑑2 − √𝑑2 − ℎ2𝑐𝑜𝑡2𝛽-------------- Eq. 8.54
The discharge q per unit length through any cross-section of the dam may be expressed as
on a portion of ox near the wall become equal to the downward forces exerted by the submerged
soil, the surface of the soil rises as shown in Fig. 8.12 (a). This heave occurs simultaneously with
108
Lecture Note on GTE-I
an expansion of the volume of the soil, which causes its permeability to increase. Additional
seepage causes the sand to boil, which accelerates the flow of water and leads to complete
failure. Terzaghi determined from model tests that heave occurs within a distance of about D/2
(where D is the depth of penetration of the pile) from the sheet pile and the critical section ox
passes through the lower edge of the sheet pile.
109
Lecture Note on GTE-I
Fig.8.14: Grain size distribution curves for graded filter and protected materials
110
Lecture Note on GTE-I
9.0 COMPACTION
In construction of highway embankments, earth dams and many other engineering structures,
loose soils must be compacted to improve their strength by increasing their unit weight;
Compaction - Densification of soil by removing air voids using mechanical equipment; the
degree of compaction is measured in terms of its dry unit weight.
Compaction, in general, is the densification of soil by removal of air, which requires mechanical
energy. The degree of compaction of a soil is measured in terms of its dry unit weight. When
water is added to the soil during compaction, it acts as a softening agent on the soil particles. The
soil particles slip over each other and move into a densely packed position. The dry unit weight
after compaction first increases as the moisture content increases.
When the moisture content is gradually increased and the same compactive effort is used for
compaction, the weight of the soil solids in a unit volume gradually increases.
Beyond a certain moisture content, any increase in the moisture content tends to reduce the dry
unit weight. This phenomenon occurs because the water takes up the spaces that would have
been occupied by the solid particles. The moisture content at which the maximum dry unit
weight is attained is generally referred to as the optimum moisture content.
111
Lecture Note on GTE-I
maximum at this stage. Further addition of more water replaces the soil grains. Thus, addition of
more water after optimum moisture content reduces the unit weight as unit weight of water is
less than the unit weight of soil grains. Because of the due to the addition of water dry unit
weight of soil increases up to a certain water content value (OMC) and beyond that addition of
more water decreases the dry unit weight of soil (as shown in Fig.9.1). This can be also
explained by concept of soil structure and electrical double layer theory.
Compactive effort: As the compactive effort increases maximum dry unit weight increases in a
decreasing rate. However, the OMC value decreases as the compactive effort increases. Thus, if
the compactive effort increases the compaction curve is shifted to the top and to the left side as
shown in Figure.2. However, as the water content increases the effect of compactive effort on
dry unit weight decreases. If the peaks of the compaction curves for different compactive efforts
are joined the obtained line is called line of optimum (as shown in Figure 2). The line of
optimum is nearly parallel to line of zero air voids. Thus, the efficiency of the compaction does
not increase as the compactive effort increases.
112
Lecture Note on GTE-I
From the preceding observation and Figure 9.3, we can see that
1. As the compaction effort is increased, the maximum dry unit weight of compaction is
also increased.
2. As the compaction effort is increased, the optimum moisture content is decreased to
some extent.
The preceding statements are true for all soils. Note, however, that the degree of compaction is
not directly proportional to the compaction effort.
Type of soil: Soil type has significant effect on compaction. It is observed that cohesive soil has
higher OMC value and poorly graded sands have lower dry unit weight as compared to the
coarse-grained, well graded soil with some percentage of fines. However, excessive amount of
fines reduces the maximum dry unit weight value of the soil. Maximum dry unit weight of sand
is obtained either completely dry or saturated condition. Compaction curve has very little
importance in case of sandy soils.
113
Lecture Note on GTE-I
compaction is increased, as shown by point B, the diffuse double layers around the particles
expand, which increases the repulsion between the clay particles and gives a lower degree of
flocculation and a higher dry unit weight. A continued increase in moisture content expands the
double layers more. This expansion results in a continued increase of repulsion between the
particles and thus a still greater degree of particle orientation and a more or less dispersed
structure. However, the dry unit weight decreases because the added water dilutes the
concentration of soil solids per unit volume. At given moisture content, higher compactive effort
yields a more parallel orientation to the clay particles, which gives a more dispersed structure.
The particles are closer and the soil has a higher unit weight of compaction.
Compaction induces variations in the structure of cohesive soils. Results of these structural
variations include changes in hydraulic conductivity, compressibility, and strength. The
hydraulic conductivity, which is a measure of how easily water flows through soil, decreases
with the increase of moisture content. It reaches a minimum value at approximately the optimum
moisture content. Beyond the optimum moisture content, the hydraulic conductivity increases
slightly. The high value of the hydraulic conductivity on the dry side of the optimum moisture
content is due to the random orientation of clay particles that results in larger pore spaces.
Some expansive clays in the field do not stay compacted, but expand upon entry of water and
shrink with loss of moisture. This shrinkage and swelling of soil can cause serious distress to the
foundations of structures. Laboratory observations such as this will help soils engineers to adopt
a moisture content for compaction to minimize swelling and shrinkage.
When a soil is compacted, it changes its engineering properties and thereby behaves differently.
Some of the engineering properties which changes on application of compactive effort is briefly
described below.
1. PERMEABILITY
The effect of compaction is to decrease the permeability. In the case of fine grained soils it has
been found that for the same dry density soil compacted wet of optimum will be less permeable
than that of compacted dry of optimum.
2. COMPRESSIBILITY
In case of soil samples initially saturated and having same void ratio, it has been found that in
low pressure range a wet side compacted soil is more compressible than a dry side compacted
soil, and vice versa in high pressure range.
3. PORE PRESSURE
In undrained shear test conducted on saturated samples of clay it has been found that lower pore
pressures develop at low strains when the sample is compacted dry of optimum, compared to the
case when the sample is compacted wet of optimum. But at high strains in both types of samples
the development of pore pressure is same for same density and water content.
4. STRESS-STRAIN RELATION
Samples compacted dry of optimum produce much steeper stress-starin curves with peaks at low
strains, whereas samples compacted wet of optimum, having the same density, produce much
flatter stress-strain curves with increase in stress even at high strains.
114
Lecture Note on GTE-I
APPARATUSREQUIRED:
1. Proctor mould having a capacity of 944 cc with an internal diameter of 10.2 cm and a height
of 11.6 cm. The mould shall have a detachable collar assembly and a detachable base plate.
2. Rammer: A mechanical operated metal rammer having a 5.08 cm diameter face and a weight
of 2.5 kg. The rammer shall be equipped with a suitable arrangement to control the height of
drop to a free fall of 30 cm.
3. Sample extruder, mixing tools such as mixing pan, spoon, towel, and spatula.
PROCEDURE:
2. Weigh the proctor mould without base plate and collar. Fix the collar and base plate. Place the
soil in the Proctor mould and compact it in 3 layers giving 25 blows per layer with the 2.5 kg
rammer falling through. The blows shall be distributed uniformly over the surface of each layer.
3. Remove the collar; trim the compacted soil even with the top of mould using a straight edge
and weigh to find the bulk density,
4. Divide the weight of the compacted specimen by 944 cc and record the result as the bulk
density.
5. Remove the sample from mould and slice vertically through and obtain a small sample for
water content.
115
Lecture Note on GTE-I
6. Thoroughly break up the remainder of the material until it will pass through 4.25mm sieve as
judged by the eye. Add water in sufficient amounts to increase the moisture content of the soil
sample by one or two percentage points and repeat the above procedure for each increment of
water added. Continue this series of determination until there is either a decrease or no change in
the wet unit weight of the compacted soil.
Calculation
First, the compaction water content (w) of the soil sample is calculated using the average of the three
measurements obtained (top, middle and bottom part of the soil mass). Subsequently, the dry
unit weight (γd) is calculated as follows:
𝛾
= 1+
𝛾 --------------- Eq.9.1
𝑑 𝑤
116
Lecture Note on GTE-I
117
Lecture Note on GTE-I
2. Add suitable amount of water with the soil and mix it thoroughly. For sandy and gravelly
soil add 3% to 5% of water. For cohesive soil the amount of water to be added should be
12% to 16% below the plastic limit.
3. Weigh the mould with base plate attached, to the nearest 1g and record the weight as W1.
Attach the extension collar with the mould. Compact the moist soil into the mould in five
layers of approximately equal mass, each layer being given 25 blows, with the help of 4.9
kg rammer, dropped from a height of 450mm above the soil. The blows must be
distributed uniformly over the surface of each layer. The operator shall ensure that the
tube of the rammer is kept clear of soil so that the rammer always falls freely.
4. After completion of the compaction operation, remove the extension collar and level
carefully the top of the mould by means of straightedge. Weigh the mould with the
compacted soil to the nearest 1 g and record this weight as W2.
5. Remove the compacted soil from the mould and place it on the mixing tray. Determine
the water content of a representative sample of the specimen. Record the moisture content
as ‘M’.
6. The remainder of the soil shall be broken up and repeat Steps (iii) to (v) above, by adding
suitable increment of water to the soil. For sandy and gravelly soils the increment is
118
Lecture Note on GTE-I
generally 1% to 2% and for cohesive soils the increment is generally 2% to 4%. The total
number of determinations made shall be at least five, and the moisture contents should be
such that the optimum moisture content, at which the maximum dry density occurs, is
within that range.
7. For compacting soil containing coarse material up to 37.5 mm size, the 2250 cm3 mould
should be used. A sample weighing about 30 kg and passing the 37.5 mm IS sieve is used
for the test. Soil is compacted in five layers, each layer being given 55 blows of the 4.9
kg rammer.
NOTE 1: It is important that the water is mixed thoroughly and adequately with the soil, since
inadequate mixing gives rise to variable test results. This is particularly important with
cohesive souls when adding a substantial quantity of water to the air-dried soil. With clays of
high plasticity, or where hand mixing is employed, it may be difficult to distribute the water
uniformly through the air-dried soil by mixing alone, and it may be necessary to store the
mixed sample in a sealed container for a minimum period of about 16 hours before continuing
with the test.
NOTE 2: It is necessary to control the total volume of soil compacted, since it has been found that if the
amount of soil struck off after removing the extension is too great, the test results will be
inaccurate.
NOTE 3: The water added for each stage of the test should be such that a range of moisture contents is
obtained which includes the optimum moisture. In general, increments of 1 to 2 percent are
suitable for sandy and gravelly soils and of 2 to 4 percent for cohesive soils. To increase the
accuracy of the test it is often advisable to reduce the increments of water in the region of the
optimum moisture content.
Calculation
119
Lecture Note on GTE-I
3. Where,
γm = Bulk density of soil in g/cm3.
w = Moisture content of soil
Graph
The dry densities γd obtained in a series of determinations is plotted against the corresponding
moisture content ‘w’. A smooth curve is then drawn through the resulting points and the position
of the maximum on this curve is determined, which is called maximum dry density (M.D.D).
And the corresponding moisture content is called optimum moisture content (O.M.C.).
Fig. 9.2: Comparison of Compaction curve for Standard and Modified Proctor Test
A line showing the relation between water content and dry density at a constant degree of
saturation S may be established from the equation:
𝛾 =
𝐺𝛾𝑤
----------------- Eq.9.4
𝜔𝐺
𝑑 1+ 𝑆
Substituting S = 95%, 90%, and so on, one can arrive at γd-values for different values of water
content in %. The lines thus obtained on a plot of γd versus w are called 95% saturation line, 90%
saturation line and so on as shown in Fig.9.2.
If one substitutes S = 100% and plots the corresponding line, one obtains the theoretical
saturation line, relating dry density with water content for a soil containing no air voids. It is said
to be ‘theoretical’ because it can never be reached in practice as it is impossible to expel the pore
air completely by compaction.
We then use
𝛾 = 𝐺𝛾𝑤
---------------- Eq.9.5
𝜔𝐺
𝑑 1+100
120
Lecture Note on GTE-I
For any type of construction job which requires soil to be used as a foundation material or as a
construction material, compaction in-situ or in the field is necessary.
The construction of a structural fill usually consists of two distinct operations—placing and
spreading in layers and then compaction. The first part assumes greater significance in major
jobs such as embankments and earth dams where the soil to be used as a construction material
has to be excavated from a suitable borrow area and transported to the work site. In this phase
large earth moving equipment such as self-propelled scrapers, bulldozers, graders and trucks are
widely employed.
The thickness of layers that can be properly compacted is known to be related to the type of soil
and method or equipment of compaction. Generally speaking, granular soils can be adequately
compacted in thicker layers than fine-grained soils and clays; also, for a given soil type, heavy
compaction equipment is capable of compacting thicker layers than light equipment. Although
the principle of compaction in the field is relatively simple, it may turn out to be a complex
process if the soil in the borrow area is not at the desired optimum moisture content for
compaction. The existing moisture content is to be determined and water added, if necessary.
Addition of water to the soil is normally done either during excavation or transport and rarely on
the construction spot; however, water must be added before excavation in the case of clayey
soils. In case the soil has more moisture content than is required for proper compaction, it has to
be air-dried after excavation and compacted as soon as the desired moisture content is attained.
Soil compaction or densification can be achieved by different means such as tamping action,
kneading action, vibration, and impact. Compactors operating on the tamping, kneading and
impact principle are effective in the case of cohesive soils, while those operating on the
kneading, tamping and vibratory principle are effective in the case of cohesionless soils.
121
Lecture Note on GTE-I
122
Lecture Note on GTE-I
the travelling wheels; this facilitates the exertion of a steady pressure on uneven ground, which is
very useful in the initial stages of a fill. The weight of the roller as well as the contact pressure is
an important parameter for the performance; the latter may be varied from 0.20 to 1 N/mm2 (2 to
10 kg/cm2) through the adjustment of air pressure in the tyres. Although this type has originated
as a towed unit, self propelled units are also available. The number of passes required is similar
to that with smooth wheeled-rollers. This type is suitable for compacting most types of soil and
has particular advantages with wet cohesive materials.
(c) Sheep foot rollers: This type of roller consists of a hollow steel drum provided with
projecting studs or feet; the compaction is achieved by a combination of tamping and kneading.
The drum can be filled with water or sand to provide and control the dead weight. As rolling is
done, most of the roller weight is imposed through the projecting feet. It is generally used as a
towed assembly with the drums mounted either singly or in pairs; self-propelled units are also
available. The feet are usually club-shaped (100 × 75 mm) or tapered (57 × 57 mm), the number
on a 50 kN (5 t) roller ranging from 64 to 88. The contact pressures of the feet may range from
700 kN/m2 (7 kg/cm2) to 4200 kN/m2 (42 kg/cm2) and weight per drum from 25 kN (2.5 t) to 130
kN (13 t). Initially, the projections sink into the loose soil and compact the soil near the lowest
portion of the layer. In subsequent passes with the roller, the zone of compaction continues to
rise until the surface is reached, when the roller is said to “Walk-out”. The length of the studs,
the contact area and the weight of roller are related to the roller performance. This type of roller
is found suitable for cohesive soils. It is unsuitable for granular soils as the studs tend to loosen
these continuously. The tendency of void formation is more in soils compacted with sheeps foot
rollers.
(d) Grid rollers: This type consists of rolls made from 38 mm steel bars at 130 mm centres, with
spaces of 90 mm square. The weight of the roller ranges from 55 kN (5.5 t) to 110 kN (11 t).
This is usually a towed unit which is suitable for many types of soil including wet clays and silts.
Rammers
This type includes the dropping type and pneumatic and internal commission type, which are
also called ‘frog rammers’. They weigh up to about 1.5 kN (150 kg) and even as much as 10 kN
(1 t) occasionally. This type may be used for cohesionless soils, especially in small restricted and
confined areas such as beds of drainage trenches and back fills of bridge abutments.
123
Lecture Note on GTE-I
Vibrators
These are vibrating units of the out-of-balance weight type or the pulsating hydraulic type. Such
a type is highly effective for cohesionless soils. Behind retaining walls where the soil is
confined, the backfill, much deeper in thickness, may be effectively compacted by vibration type
of compactors.
A few of this type are dealt with below:
(a) Vibrating drum: A separate motor drives an arrangement of eccentric weights so as to cause a
high-frequency, low-amplitude, vertical oscillation to the drum. Smooth drums as well as sheep
foot type of drums may be used. Layers of the order of 1 meter deep could be compacted to high
densities.
(b) Vibrating pneumatic tyre: A separate vibrating unit is attached to the wheel axle. The ballast
box is suspended separately from the axle so that it does not vibrate. A 300 mm thick layer of
granular soil will be satisfactorily compacted after a few passes.
(c) Vibrating plate: This typically consists of a number of small plates, each of which is operated
by a separate vibrating unit. These have a limited depth of effectiveness and hence are used in
compacting granular base courses for highway and airfield pavements.
(d) Vibroflot: A method suited for compacting thick deposits of loose sandy soil is called the
‘vibro flotation’ process. The improvement of density is restricted to the surface zone in the case
of conventional compaction equipment. The vibroflotation method first compacts deep zone in
the soil and then works its way towards the surface. A cylindrical vibrator weighing about 20 kN
(2 t) and approximately 400 mm in diameter and 2 m long, called the ‘Vibroflot’, is suspended
from a crane and is jetted to the depth where compaction is to start. The jetting consists of a
water jet under pressure directed into the earth from the tip of the vibroflot; as the sand gets
displaced, the vibroflot sinks into the soil. Depths up to 12 m can be reached. After the vibroflot
is sunk to the desired depth, the vibrator is activated. The compaction of the soil occurs in the
horizontal direction up to as much as 1.5 m outward from the vibroflot. Vibration continues as
the vibroflot is slowly raised toward the surface. As this process goes on, additional sand is
continually dropped into the space around the vibroflot to fill the void created. To densify the
soil in a given site, locations at approximately 3-m spacings are chosen and treated with vibro
flotation.
9.8 CONTROL OF COMPACTION IN THE FIELD
Control of compaction in the field consists of checking the water content in relation to the
laboratory optimum moisture content and the dry unit weight achieved in-situ in relation to the
laboratory maximum dry unit weight from a standard compaction test. Typically, each layer is
tested at several random locations after it has been compacted.
Several methods are available for the determination of in-situ unit weight and moisture content
The common approaches for the determination of unit weight are the core-cutter method and
sand-replacement method. A faster method is what is known as the Proctor needle method,
which may be used for the determination of in-situ unit weight as well as in-situ moisture
content. The required density can be specified either by ‘relative compaction’ (also called
124
Lecture Note on GTE-I
‘degree of compaction’) or by the final air-void content. Relative compaction means the ratio of
the insitu dry unit weight achieved by compaction to the maximum dry unit weight obtained
from an appropriate standard compaction test in the laboratory. Usually, the relative compaction
of 90 to 100% (depending upon the maximum laboratory value), corresponding to about 5 to
10% air content, is specified and sought to be achieved. Typical values of dry unit weights
achieved may be as high as 22.5 kN/m 3 (2250 kg/m3) for well-graded gravel and may be as low
as 14.4 kN/m3 (1440 kg/m3) for clays. Approximate ranges of optimum moisture content may be
6 to 10% for sands, 8 to 12% for sand-silt mixtures, 11 to 15% for silts and 13 to 21% for clays
(as got from modified AASHO tests).
A variation of 5 to 10% is allowed in the field specification of dry unit weight at random
locations, provided the average is about the specified value.
9.8.1 Proctor Needle
The Proctor needle approach given here is an efficient and fast one for the simultaneous
determination of in-situ unit weight and in-situ moisture content, it is also called ‘penetration
needle’. The apparatus basically consists of a needle attached to a spring-loaded plunger through
a shank. An array of interchangeable needle tips is available, ranging from 6.45 to 645 mm 2, to
facilitate the measurement of a wide range of penetration resistance values. A calibration of
penetration against dry unit weight and water content is obtained by pushing the needle into
specially prepared samples for which these values are known and noting the penetration. The
penetration of the needle and the penetration resistance (load applied) may be shown on a
graduated scale on the shank and the stem of handle respectively.
The procedure for the use of the Proctor ‘plasticity’ needle, as it is called, is obvious. The spring-
loaded plunger is pressed into the compacted layer in the field with an appropriate plasticity
needle. The penetration resistance is recorded for a standard depth of penetration at a standard
time-rate of penetration. Against this penetration resistance, the corresponding values of water
content and dry unit weight are obtained from the calibration curve.
The size of the needle to be chosen depends upon the type of soil such that the resistance to be
read is neither too large nor too small. The Tennessee Valley Authority (TVA) engineers had
devised a similar device, which is called the TVA ‘Penetrometer’.
9.8.2 Special Compaction Techniques
Several special types of compaction techniques have been developed for deep compaction of in-
place soils, and these techniques are used in the field for large-scale compaction works. Among
these, the popular methods are vibroflotation, dynamic compaction, and blasting. Details of these
methods are provided in the following sections.
9.8.3 Vibro-flotation
Vibro-flotation is a technique for in situ densification of thick layers of loose granular soil
deposits. It was developed in Germany in the 1930s. The first vibroflotation device was used in
the United States about 10 years later. The process involves the use of a Vibroflot Unit.
125
Lecture Note on GTE-I
This vibrating unit has an eccentric weight inside it and can develop a centrifugal force, which
enables the vibrating unit to vibrate horizontally. There are openings at the bottom and top of the
vibrating unit for water jets.
The vibrating unit is attached to a follow-up pipe. Figure 9.4 shows the entire assembly of
equipment necessary for conducting the field compaction.
126
Lecture Note on GTE-I
Stage 1: The jet at the bottom of the Vibroflot is turned on and lowered into the ground.
Stage 2: The water jet creates a quick condition in the soil and it allows the vibrating unit to sink
into the ground.
Stage 3: Granular material is poured from the top of the hole. The water from the lower jet is
transferred to the jet at the top of the vibrating unit.
This water carries the granular material down the hole.
Stage 4: The vibrating unit is gradually raised in about 0.3 m lifts and held vibrating for about 30
seconds at each lift. This process compacts the soil to the desired unit weight.
9.9 DYNAMIC COMPACTION
Dynamic compaction is a technique that has gained popularity in the United States for the
densification of granular soil deposits. This process consists primarily of dropping a heavy
weight repeatedly on the ground at regular intervals. The weight of the hammer used varies over
a range of 80 to 360 kN, and the height of the hammer drop varies between 7.5 and 30.5 m. The
stress wave generated by the hammer helps in the densification.
The degree of compaction achieved at a given site depends on the following three factors:
1. Weight of hammer
2. Height of hammer drop
3. Spacing of locations at which the hammer is dropped
Figure 9.6, shows a dynamic compaction in progress. A site immediately after the completion is
shown in Figure 9.6 b.
127
Lecture Note on GTE-I
Leonards, Cutter, and Holtz (1980) suggested that the significant depth of influence for
compaction can be approximated by using the equation
𝐷=
1
ℎ------------------
√𝑊
Eq.9.6
2 𝐻
where D =significant depth of densification (m)
WH = dropping weight (metric ton)
h = height of drop (m)
In English units, the preceding equation takes the form
𝐷 = 0.61√𝑊𝐻ℎ------------------------- Eq.9.7
where the units of D and h are ft, and the unit of WH is kip.
In 1992, Poran and Rodriguez suggested a rational method for conducting dynamic compaction
for granular soils in the field. According to their method, for a hammer of width D having a
weight WH and a drop h, the approximate shape of the densified area will be of the type shown in
Figure 9.7, a (i.e., a semiprolate spheroid).
128
Lecture Note on GTE-I
Note that in this figure b = DI (where DI is the significant depth of densification). Figure 9.8
gives the design chart for a/D and b/D versus NWHh/Ab (D is the width of the hammer if not
circular in cross section; A = area of cross section of the hammer; and N is the number of
required hammer drops). This method uses the following steps.
Step 2: Determine the hammer weight (W H), height of drop (h), dimensions of the cross section,
and thus, the area A and the width D.
Step 4: Use Figure 6.39 and determine the magnitude of NWHh/Ab for the value of b/D obtained
in step 3.
Step 5: Since the magnitudes of W H, h, A, and b are known (or assumed) from step 2, the
number of hammer drops can be estimated from the value of NWH h/Ab obtained from step 4.
Step 6: With known values of NWH h/Ab, determine a/D and thus a from Figure 9.8.
Step 7: The grid spacing, Sg, for dynamic compaction may now be assumed to be equal to or
somewhat less than a. (See Fig. 9.9)
9.10 BLASTING
Blasting is a technique that has been used successfully in many projects (Mitchell, 1970) for the
densification of granular soils. The general soil grain sizes suitable for compaction by blasting
are the same as those for compaction by vibroflotation.
The process involves the detonation of explosive charges, such as 60% dynamite at a certain
depth below the ground surface in saturated soil. The lateral spacing of the charges varies from
about 3 to 9 m. Three to five successful detonations are usually necessary to achieve the desired
compaction. Compaction (up to a relative density of about 80%) up to a depth of about 18 m
129
Lecture Note on GTE-I
over a large area can easily be achieved by using this process. Usually, the explosive charges are
placed at a depth of about two-thirds of the thickness of the soil layer desired to be compacted.
The sphere of influence of compaction by a 60% dynamite charge can be given as follows
(Mitchell, 1970):
r = √
Wex
C
--------------------- Eq.9.8
Example 9.1: An earth embankment is compacted at a water content of 18% to a bulk density of
19.2 kN/m3. If the specific gravity of the sand is 2.7, find the void ratio and the degree of
saturation of the compacted embankment.
Sol:
Problem 1: The soil in a borrow pit has a void ratio of 0.90. A fill-in-place volume of 20,000 m3
is to be constructed with an in-place dry density of 18.84 kN/m 3. If the owner of borrow area is
to be compensated at Rs.1.50 per cubic metre of excavation, determine the cost of compensation.
130
Lecture Note on GTE-I
In the arithmetic plot as shown, as the soil compresses, for the same increase of effective
stress σ', the void ratio reduces by a smaller magnitude, from ∆e1 to ∆e2. This is on account of an
increasingly denser packing of the soil particles as the pore water is forced out. In fine soils, a
much longer time is required for the pore water to escape, as compared to coarse soils.
It can be said that the compressibility of a soil decreases as the effective stress increases. This
can be represented by the slope of the void ratio – effective stress relation, which is called
the coefficient of compressibility, av.
131
Lecture Note on GTE-I
𝑎𝑣 ∆𝜎′
For a small range of effective stress, Eq.10.1
=−
𝑑𝑒________________________________
𝑎𝑣
𝑑𝜎′
or, Eq.10.2
If the current effective stress, 𝝈′ , is equal (note that it cannot be greater than) to the pre-
consolidation stress, then the deposit is said to be normally consolidated (NC). If the current
effective stress is less than the pre-consolidation stress, then the soil is said to be over-
consolidated (OC).
The figure shows the relation of void ratio and effective stress of a clay soil as a semi-log plot.
132
Lecture Note on GTE-I
133
Lecture Note on GTE-I
The pre consolidation stress, 𝝈𝒑′ , is defined to be the maximum effective stress experienced by
along the path that it would have followed if loaded from O to R continuously.
𝒄
the soil. This stress is identified in comparison with the effective stress in its present state. For
soil at state Q or F, this would correspond to the effective stress at point P.
It may be seen that for the same increase in effective stress, the change in void ratio is much less
for an over-consolidated soil (from e0 to ef), than it would have been for a normally consolidated
soil as in path OP. In unloading, the soil swells but the increase in volume is much less than the
initial decrease in volume for the same stress difference.
The distance from the normal consolidation line has an important influence on soil behaviour.
This is described numerically by the over-consolidation ratio (OCR), which is defined as the
ratio of the pre-consolidation stress to the current effective stress.
𝜎′
𝑂𝐶𝑅
𝑝𝑐
=
𝜎′
--------------------- Eq.10.4
=
�
𝜎
---------------------- Eq.10.5
𝑐 ′2
∆𝑒
�
𝜎1
𝑙𝑜𝑔1 ′
0
134
Lecture Note on GTE-I
135
Lecture Note on GTE-I
It is possible to determine the pre-consolidation stress that the soil had experienced. The soil
sample is to be loaded in the laboratory so as to obtain the void ratio - effective stress
relationship. Empirical procedures are used to estimate the pre-consolidation stress, the most
widely used being Casagrande's construction which is illustrated as
The steps in the construction are:
The total stress increases when additional vertical load is first applied. Instantaneously, the pore
water pressure increases by exactly the same amount. Subsequently there will be flow from
regions of higher excess pore pressure to regions of lower excess pore pressure causing
dissipation. The effective stress will change and the soil will consolidate with time. This is
shown schematically.
136
Lecture Note on GTE-I
137
Lecture Note on GTE-I
the beginning of transient flow, simulating the phenomenon of consolidation. The openness of
the valve is analogous to the permeability of soil.
Referring to Fig.10.5 (d), flow has occurred to the extent of dissipating 50% of the excess pore
pressure. The pore water pressure at this instant is half the initial value, i.e., 1/2(δP/A). This
causes a corresponding increase in the stress in the spring of ½ (δP/A), the total stress remaining
constant at [(P/A) + (δP/A)]. This stage refers to that of “50% consolidation”.
Referring to Fig.10.5 (e), the final equilibrium condition is reached when the transient flow
situation ceases to exist, consequent to the complete dissipation of the pore water pressure. The
spring compresses to a final height Hf < Hi, carrying the total stress of (P + δP)/A, all by itself,
since the excess pore water pressure has been reduced to zero, the pressure in it having equalled
the atmospheric. The system has reached the equilibrium condition under the load (P + δP). This
represents “100% consolidation” under the applied load or stress increment. We may say that the
“soil” has been consolidated to an effective stress of (P + δP)/A.
In this mechanistic model, the compressible soil skeleton is characterised by the spring and the
pore water by the water in the cylinder. The more compressible the soil, the longer the time
required for consolidation; the more permeable the soil, the shorter the time required. There is
one important aspect in which this analogy fails to simulate consolidation of a soil. It is that the
pressure conditions are the same throughout the height of the cylinder, whereas the consolidation
of a soil begins near the drainage surfaces and gradually progresses inward. In may be noted that
soil consolidates only when effective stress increases; that is to say, the volume change
behaviour of a soil is a function of the effective stress and not the total stress.
Similar arguments may be applied to the expansion characteristics under the decrease of load.
An alternative mechanical analogy to the consolidation process is shown in Fig. 10.6. A cylinder
is fitted with a number of pistons connected by springs to one another. Each of the compartments
thus formed is connected to the atmosphere with the aid of standpipes. The cylinder is full of
water and is considered to be airtight. The pistons are provided with perforations through which
water can move from one compartment to another. The topmost piston is fitted with valves
which may open or close to the atmosphere. It is assumed that any pressure applied to the top
piston gets transmitted undiminished to the water and springs.
Initially, the cylinder is full of water and weights of the pistons are balanced by the springs; the
water is at atmospheric pressure and the valves may be open. The water level stands at the
elevation PP in the standpipes as shown. The valves are now closed, the water level continuing
to remain at PP. An increment of pressure∆𝜎 is applied on the top piston. It will be observed
that the water level rises instantaneously in all the stand pipes to an elevation
QQ, above PP by a height h = Δσ/γw. Let all the valves be opened simultaneously with the
application of the pressure increment, the time being reckoned from that instant. The height of
the springs remains unchanged at that instant and the applied increment of pressure is fully taken
up by water as the hydrostatic excess pressure over and above the atmospheric. An equal rise of
water in all the standpipes indicates that the hydrostatic excess pressure is the same in all
compartments immediately after application of pressure. As time elapses, the water level in the
138
Lecture Note on GTE-I
pipes starts falling, the pistons move downwards gradually and water comes out through the
open valves.
139
Lecture Note on GTE-I
3. The soil grains and water are virtually incompressible (γw is constant and volume change of
soil is only due to change in void ratio).
4. The behaviour of infinitesimal masses in regard to expulsion of pore water and consequent
consolidation is no different from that of larger representative masses (Principles of calculus may
be applied).
5. The compression is one-dimensional (u varies with z only).
6. The flow of water in the soil voids is one-dimensional, Darcy’s law being valid.
𝛛𝑣𝑦
= = 0 and 𝑣 = 𝑘
𝛛𝑣𝑥 𝛛ℎ
--------------- Eq.10.6
𝛛 𝛛 𝑧 𝑧 𝛛𝑧
𝑥 𝑦
7. Certain soil properties such as permeability and modulus of volume change are constant; these
actually vary somewhat with pressure. (k and mv are independent of pressure).
8. The pressure versus void ratio relationship is taken to be the idealised one, av is constant).
9. Hydrodynamic lag alone is considered and plastic lag is ignored, although it is known to exist.
(The effect of k alone is considered on the rate of expulsion of pore water).
The first three assumptions represent conditions that do not vary significantly from actual
conditions.
The fourth assumption is purely of academic interest and is stated because the differential
equations used in the derivation treat only infinitesimal distances. It has no significance for the
laboratory soil sample or for the field soil deposit.
The fifth assumption is certainly valid for deeper strata in the field owing to lateral confinement
and is also reasonably valid for an oedometer sample.
The sixth assumption regarding flow of pore water being one-dimensional may be taken to be
valid for the laboratory sample, while its applicability to a field situation should be checked.
However, the validity of Darcy’s law for flow of pore water is unquestionable.
The seventh assumption may introduce certain errors in view of the fact that certain soil
properties which enter into the theory vary somewhat with pressure but the errors are considered
to be of minor importance.
The eighth and ninth assumptions lead to the limited validity of the theory. The only justification
for the use of the eighth assumption is that, otherwise, the analysis becomes unduly complex.
The ninth assumption is necessitated because it is not possible to take the plastic lag into account
in this theory. These two assumptions also may be considered to introduce some errors.
Now let us see the derivation of Terzaghi’s theory with respect to the laboratory oedometer
sample with double drainage as shown in Fig. 10.7.
140
Lecture Note on GTE-I
𝑘 +𝑘 [𝑒 + 𝑆 ]------------------
2
Equation 2,𝛛for flow of𝛛2water
ℎ ℎ 1 𝛛𝑠 𝛛𝑒
through soil, holds here also,
𝑥 𝛛𝑥2 𝑧
= Eq.10.7
1+ 𝛛 𝛛𝑡
𝛛𝑧2 𝑒 𝑡
𝑘 [𝑒 + 𝑆 ]------------------
2
For 𝛛one-dimensional
ℎ 1 𝛛𝑠
flow𝛛𝑒situation, this reduces to:
𝑧
= Eq.10.8
1+ 𝛛 𝛛𝑡
𝛛𝑧2 𝑒 𝑡
During the process of consolidation, the degree of saturation is taken to remain constant at 100%,
while void ratio changes causing reduction in volume and dissipation of excess hydrostatic
𝛛𝑡1
Hence 𝑘 [ ] = [ ]------------------
𝛛2 ℎ 𝛛𝑒 𝛛 𝑒
𝑧 𝛛𝑧2
=− Eq.10.9
1+ 𝛛 𝛛𝑡 1+𝑒
𝑒 𝑡
Negative sign denoting decrease of e for increase of h.
grains are virtually incompressible, 𝛛 [ ] represents time-rate of volume change per unit
𝑒
Since volume decrease can be due to a decrease in the void ratio only as the pore water and soil
𝛛 1+𝑒
volume 𝑡
The flow is only due to the hydrostatic excess pressure, ℎ =
𝑢
𝛾𝑤
=−
2
𝑘 𝛛 𝑢 𝛛𝑉
𝛾𝑤
So -------------- Eq.10.10
𝛛𝑡
𝛛𝑧2
Here k is the permeability of soil in the direction of flow, and ∂V represents the change in
volume per unit volume. The change in hydrostatic excess pressure, ∂u, changes the intergranular
or effective stress by the same magnitude, the total stress remaining constant.
The change in volume per unit volume, ∂V, may be written, as per the definition of the modulus
141
Lecture Note on GTE-I
𝜕𝑉 = 𝑚𝑣𝜕𝜎 = −𝑚𝑣𝜕𝑢----------------
of volume change, mv
Eq.10.11
142
Lecture Note on GTE-I
Negative sign is used, since increase in stress reduces pore water pressure.
=
Differentiating both sides with respect to time,
𝛛𝑉
𝛛𝑢
−𝑚 � 𝛛𝑡
-------------------- Eq.10.12
𝛛𝑡 �
=
𝑘 𝛛2 𝑢
From Eq. 10.11 and 10.12, we get
𝛛𝑢
---------------- Eq.10.13
𝛛𝑡 𝛾𝑤𝑚𝑣 𝛛𝑧2
= 𝑐
2𝛛 𝑢
This is written as:
𝛛𝑢
𝑣 𝛛𝑧2
---------------- Eq.10.14
𝛛𝑡
= 𝛾𝑤𝑚𝑣
𝑘
𝑐𝑣
Where
= 𝛾𝑤𝑚
𝑘 ------------- Eq.10.15
𝑣 𝑎𝑣𝛾𝑤
𝑣
Equation 10.14 is the basic differential equation of consolidation according to Terzaghi’s theory
of one-dimensional consolidation. The coefficient of consolidation combines the effect of
permeability and compressibility characteristics on volume change during consolidation. Its units
can be shown to be mm2/s or L2 T-1.
The initial hydrostatic excess pressure, ui, is equal to the increment of pressure Δσ, and is the
same throughout the depth of the sample, immediately on application of the pressure, and is
shown by the heavy line in Fig. 10.7 (b). The horizontal portion of the heavy line indicates the
fact that, at the drainage face, the hydrostatic excess pressure instantly reduces to zero,
theoretically speaking. Further, the hydrostatic excess pressure would get fully dissipated
throughout the depth of the sample only after the lapse of infinite time*, as indicated by the
heavy vertical line on the left of the figure. At any other instant of time, the hydrostatic excess
pressure will be maximum at the farthest point in the depth from the drainage faces, that is, at the
middle and it is zero at the top and bottom. The distribution of the hydrostatic excess pressure
with depth is sinusoidal at other instants of time, as shown by dotted lines. These curves are
called “Isochrones”.
10.4.1 Alternative Method:
=
1 𝛛𝑢
With reference to Fig.10.7, the hydraulic gradient i1 at depth z
𝛛ℎ
--------------- Eq.10.16
𝛛𝑧 𝛾𝑤 𝛛𝑧
The hydraulic gradient i at depth z+𝜕𝑧 = ( +2𝛛 𝑑𝑧)-----------
1 𝛛𝑢 𝑢
Eq.10.17
𝛾𝑤 𝛛𝑧 𝛛𝑧2
2
Rate of inflow per unit area = Velocity at depth z = k.i1, by Darcy’s law.
143
Lecture Note on GTE-I
144
Lecture Note on GTE-I
∴ Change of volume = – mv∂(σ – u).dz, since the elementary layer of thickness dz and unit cross-
(The negative sign denotes decrease in volume with increase in pressure).
∴ Time-rate of change of 𝛛𝑢
� 𝛛𝑡
𝑑𝑧----------------- Eq.10.20
volume=+𝑚 �
𝑑𝑧---------------
𝑘 𝛛2𝑢 𝛛𝑢
Equating this to water lost per unit time
𝑣 𝛛𝑡
dz=-𝑚 Eq.10.21
𝛾𝑤𝑚𝑣 𝛛𝑧2
= 𝑐
𝛛𝑢 𝛛 𝑢
2
𝑣 𝛛𝑧2
---------------------- Eq.10.22
𝛛𝑡
= 𝛾𝑤𝑚
𝑘
𝑐𝑣
Where ---------------------------------------- Eq.10.23
𝑣
10.5 SOLUTION OF TERZAGHI’S EQUATION FOR ONE-DIMENSIONAL
CONSOLIDATION
Terzaghi solved the differential equation for a set of boundary conditions which have utility in
solving numerous engineering problems and presented the results in graphical form using
dimensional parameters.
The following are the boundary conditions:
1. There is drainage at the top of the sample: At z = 0, u = 0, for all t.
2. There is drainage at the bottom of the sample: At z = 2H, u = 0, for all t.
3. The initial hydrostatic excess pressure ui is equal to the pressure increment, Δσ u = ui = Δσ, at
t = 0.
Terzaghi chose to consider this situation where u = ui initially throughout the depth, although
solutions are possible when ui varies with depth in any specified manner. The thickness of the
sample is designated by 2H, the distance H thus being the length of the longest drainage path,
i.e., maximum distance water has to travel to reach a drainage face because of the existence of
two drainage faces. (In the case of only one drainage face, this will be equal to the total thickness
of the clay layer).
The general solution for the above set of boundary conditions has been obtained on the basis of
𝑢 = 𝑓(𝑧, 𝑛𝜋𝑧 𝑡) = ∑∞
1 2𝐻
[ 𝑢 𝑠𝑖𝑛
𝑛𝜋𝑧
𝑑𝑧] (𝑠𝑖𝑛 )𝑒
separation of variables and Fourier series expansion and is as follows:
−𝑛 𝜋 𝐶𝑣𝑡/4𝐻
2 2 2
∫
-- Eq.10.24
𝑛=1
𝐻 0 2 2𝐻
𝐻
𝑖
This solution enables the hydrostatic excess u to be computed for a soil mass under any initial
system of stress ui, at any depth z, and at any time t.
𝑛=1
Eq.10.25
145
Lecture Note on GTE-I
𝑛𝜋 2𝐻
When n is even, (1 – cosnπ) vanishes; when n is odd, this factor becomes 2. Therefore it is
146
Lecture Note on GTE-I
𝑈 = 𝑖 = (1 − )-------------------
relation𝑢to−𝑢the initial value:
𝑢
Eq.10.27
𝑢𝑖 𝑢 𝑖
𝑧
The subscript z is significant, since the extent of dissipation of excess pore water pressure is
different for different locations, except at the beginning and the end of the consolidation process.
The third dimensionless parameter, relating to time, and called ‘Time-factor’, T, is defined as
𝑇 =
𝑐𝑣𝑡
follows:
-------------------- Eq.10.28
𝐻2
Where cv is the coefficient of consolidation,
H is the length drainage path, and t is the elapsed time from the start of consolidation process.
In the context of consolidation process at a particular site, cv and H are constants, and the time
𝑀 𝐻
Eq.10.29
The one-dimensional consolidation testing procedure was first suggested by Terzaghi. This test
is performed in a consolidometer (sometimes referred to as an oedometer). The schematic
diagram of a consolidometer is shown in Figure 10.8(a). Figure 10.8(b) shows a photograph of a
consolidometer. The soil specimen is placed inside a metal ring with two porous stones, one at
the top of the specimen and another at the bottom. The specimens are usually 64 mm in diameter
and 25 mm thick. The load on the specimen is applied through a lever arm, and compression is
measured by a micrometer dial gauge. The specimen is kept under water during the test. Each
load usually is kept for 24 hours. For consolidation testing, it is generally desirable that the
applied pressure at any loading stage be double than that at the preceding stage. The test may,
therefore, be continued using a loading sequence which would successively apply stress of 0.1,
0.2, 0.4, 0.8, 1.6, 3.2, etc, kgf/cm 2 on the soil specimen. For each loading increment, after
application of load, readings of the dial gauge shall be taken using a time sequence such as 0,
0.25, 1, 2.25, 4,6·25,9, 12·25, 16, 20.25, 25, 36, 49, 64, 81, 100, 121, 144, 169, 196, 225 min etc,
up to 24. These time sequences facilitate plotting of thickness or change of thickness of specimen
147
Lecture Note on GTE-I
against square root of time or against logarithm of time. The loading increment shall be left at
least until the slope of the characteristic linear secondary compression portion of thickness
versus log time plot is apparent or until the end of pre-consolidation is Indicated on a square roof
of time plot. A period of 24 h will usually be sufficient, but longer times may be required. In
every case, the same load increment duration shall be used for all load increments during a
consolidation test.
On completion of the final loading stage, the specimen shall be unloaded by pressure decrements
which decrease the load to one-fourth of the last load. Dial gauge readings may be taken as
necessary during each stage of unloading. If desired, the time intervals used during the
consolidation increments may be adopted; usually it is possible to proceed much more rapidly.
After the completion of test, dry weight of the soil sample is taken.
Figure 1(c) shows a consolidation test in progress. The general shape of the plot of deformation
of the specimen against time for a given load increment is shown in Figure 10.9. From the plot,
we can observe three distinct stages, which may be described as follows:
Stage I: Initial compression, which is caused mostly by preloading
Stage II: Primary consolidation, during which excess pore water pressure gradually is transferred
into effective stress because of the expulsion of pore water
Stage III: Secondary consolidation, which occurs after complete dissipation of the excess pore
water pressure, when some deformation of the specimen takes place because of the plastic
readjustment of soil fabric
148
Lecture Note on GTE-I
Fig. 10.8 (b): Photograph of a consolidometer Fig. 10.8 (c): Consolidation test in progress
149
Lecture Note on GTE-I
After the time–deformation plots for various loadings are obtained in the laboratory, it is
necessary to study the change in the void ratio of the specimen with pressure. Following is a
step-by-step procedure for doing so:
𝜌𝑤 =density of water
𝐻𝑉 = 𝐻 − 𝐻𝑆----------------------
Step 2: Calculate the initial height of voids as
Eq.10.34
Where, HV = initial height of the specimen.
𝑒 = = =
Step 3:𝑉Calculate
𝐻𝑉 𝐴 𝐻𝑉
the initial void ratio, of the specimen, using the equation
𝑉
---------------------- Eq.10.35
𝐻𝑆 𝐻𝑆
0
𝑉𝑆 𝐴
∆𝑒 =
∆𝐻1
---------------- Eq.10.36
𝐻𝑆
1
∆H1 is obtained from the initial and the final dial readings for the loading).
It is important to note that, at the end of consolidation, total stress 𝜎1 is equal to effective stress
𝑒1 = 𝑒0 − ∆𝑒1------------------
Step 5: Calculate the new void ratio after consolidation caused by the pressure increment as
For the next loading, 𝜎2 (note: 𝜎2 equals the cumulative load per unit area of specimen), which
Eq.10.37
causes additional deformation ∆H2, the void ratio at the end of consolidation can be calculated
as
𝑒 = −
∆𝐻2
2 𝑒1
----------------- Eq.10.38
𝐻2
150
Lecture Note on GTE-I
At this time, 𝜎2′ effective stress, proceeding in a similar manner, one can obtain the void ratios at
the end of the consolidation for all load increments.
151
Lecture Note on GTE-I
The effective stress 𝜎′ and the corresponding void ratios (e) at the end of consolidation are
plotted on semi logarithmic graph paper. The typical shape of such a plot is shown in Figure
10.11
152
Lecture Note on GTE-I
𝑡90
---------------------------------- Eq.10.39
𝑣
where t90 is read off from Fig. 10.7 (a)
T90 is 0.848 from Terzaghi’s theory
H is the drainage path, which may be taken as half the thickness of the sample for double
drainage conditions,
153
Lecture Note on GTE-I
Fig. 10.7: Square root of time fitting method (After Taylor, 1948)
Fig. 10.8(a): Sample thickness/Dial gauge reading versus logarithm of time (Laboratory curve)
154
Lecture Note on GTE-I
Fig. 10.8 (b) Logarithm of time fitting method (After A. Casagrande, 1939)
Since the early portion of the curve is known to approximate a parabola, the corrected zero point
may be located as follows: The difference in ordinates between two points with times in the ratio
of 4 to 1 is marked off; then a distance equal to this difference may be stepped off above the
upper points to obtain the corrected zero point. This point may be checked by more trials, with
different pairs of points on the curve. After the zero and 100% primary compression points are
located, the point corresponding to 50% consolidation and its time may easily be obtained and
the coefficient of consolidation computed from:
𝑐 =
𝑇50𝐻2
𝑡50
----------------- Eq.10.40
𝑣
Where t50 is read off from Fig. 10.8 (a)
T50 = 0.197 from Terzaghi’s theory, and H is the drainage path as stated in the previous
subsection.
155
Lecture Note on GTE-I
156
Lecture Note on GTE-I
𝑒 = 𝑙𝑜𝑔1 ( )----------------------
𝑡2
−𝐶𝛼
Eq.10.42
𝑡1
0
In other words, 𝐶𝛼 may be taken to be the slope of the straight line representing the secondary
𝐶
= 1+
𝛼 --------------------- Eq.10.43
𝛼 𝑒
Generally α and 𝐶𝛼 increase with increasing stress
Some common values of 𝐶𝛼 are given below
Sl. No Nature of Soil 𝐶𝛼 – Value
1 Over consolidated days 0.0005 to 0.0015
2 Normally consolidated days 0.005 to 0.030
3 Organic soils, peats 0.04 to 0.10
average effective overburden pressure, Because of an increase of effective pressure, 𝜎′, let the
Let us consider a saturated clay layer of thickness H and cross-sectional area A under an existing
∆𝑉 = 𝑉0 − 𝑉1 = 𝐻𝐴 − (𝐻 − 𝑆𝐶)𝐴--------
primary settlement be Sc. Thus, the change in volume (Figure 10.11) can be given by
Hence ∆𝑉 = 𝑆𝐶𝐴-----------------------------
Eq.10.44
Eq.10.43
∆𝑉 = ∆𝑒𝑉𝑆------------------------------------
of void ratio, it follows that
Where, ∆e
Eq. 10.45
𝑉 = =
𝑉0
is the𝐴𝐻change of void ratio. But
𝑆
----------------------------- Eq. 10.46
1+𝑒0 1+𝑒0
Where, 𝑒0 is initial void ratio at volume V0. Thus, from Eqs. (10.43) through (10.46),
∆𝑉 = 𝑆 𝐴 = ∆𝑒𝑉 = 𝐴𝐻 ∆𝑒----------------
𝐶 𝑆
Eq.10.47
1+𝑒0
𝑆 =𝐻
∆𝑒
----------------------------------- Eq. 10.48
𝐶 1+𝑒0
158
Lecture Note on GTE-I
𝑆𝐶
𝐻𝑙𝑜𝑔10+∆𝜎
( 0 )-------------------------
′ ′
𝑐𝑐 𝜎
= 1+𝑒 𝜎
′0
Eq.10.50
0
In over-consolidated clays (see Figure 10.13), for field e-log 𝜎′ variation will be along the line
Fig.10.13: Consolidation characteristics of over consolidated clays
hj, the slope of which will be approximately equal to that for the laboratory rebound curve. The
𝑆𝐶
𝐻𝑙𝑜𝑔10+∆𝜎
( 0 )-------------------
′ ′
𝑐𝑠 𝜎
= 1+𝑒 𝜎
′0
Eq.10.52
0 𝑐 ′
If
𝑐𝑠 𝜎 𝜎′ +∆𝜎′
𝑆𝐶 𝐻𝑙𝑜𝑔10 ( ) + 𝐻𝑙𝑜𝑔10( )-
𝑐 𝑐 𝑐 0
=
Eq.10.53
1+𝑒0 𝜎0
′
1+𝑒0 � �′
� �
10.9 CORRELATIONS FOR COMPRESSION INDEX (CC)
The compression index for the calculation of field settlement caused by consolidation can be
determined by graphic construction (as shown in Figure 10.13) after one obtains the laboratory
test results for void ratio and pressure.
Skempton (1944) suggested the following empirical expression for the compression index for
𝑐𝑐 = 0.009(𝐿𝐿 − 10%)---------------
undisturbed clays:
Eq.10.54
Where LL is the liquid limit
Several other correlations for the compression index are also available. They have been
developed by tests on various clays. Some of these correlations are given in Table 1.
On the basis of observations on several natural clays, Rendon-Herrero (1983) gave the
𝑐 = 0.14𝐺 ( ) -------------
1+𝑒𝑜 2.38
relationship for the compression index in the form
Eq.10.55
𝑠 𝐺𝑠
𝑐
𝑐 = 0.2343 ( ) 𝐺 --------------
𝐿𝐿%
Nagaraj and Murty (1985) expressed the compression index as
𝑐 𝑠
Eq.10.56
100
More recently, Park and Koumoto (2004) expressed the compression index by the following
relationship:
𝑐
= 371.747−4.275 -----------------
𝜂0 Eq.10.57
𝑐
𝜂0
where 𝜂0 = in situ porosity of the soil.
𝑐 = 𝑡𝑜 𝑐 ------------------
can be 1determined
1
from laboratory tests. In most cases,
𝑠 𝑐
Eq.10.58
5 10
161
Lecture Note on GTE-I
𝑐 = 0.0463 ( ) 𝐺 ---------------
𝐿𝐿%
The swell index was expressed by Nagaraj and Murty (1985) as
𝑠 𝑠
Eq.10.59
100
Based on the modified Cam clay model, Kulhawy and Mayne (1990) have shown that
𝑐 =
𝑃𝐼
370
-------------------- Eq.10.60
𝑠
Where PI= Plasticity Index
Sol:
Sol:
148
Lecture Note on GTE-I
Sol:
Sol:
149
Lecture Note on GTE-I
Sol:
150
Lecture Note on GTE-I
Sol:
151
Lecture Note on GTE-I
Example 7:
Sol:
152
Lecture Note on GTE-I
Example 8:
Sol:
153
Lecture Note on GTE-I
Example 9:
Sol:
154
Lecture Note on GTE-I
Example 10
Sol:
155
Lecture Note on GTE-I
Example 11:
Sol:
156
Lecture Note on GTE-I
∆𝒆 𝟏 + 𝒆𝟏 + ∆𝒆
Now
𝟏+
=
∆𝑯 𝒆𝟎 𝑯𝟎
=
𝑯𝟎
∆𝑒
= 19.0
1.541+∆𝑒
3.52
Hence
0
OR ∆𝑒 = 0.350 and 𝑒0 = 0.541 + 0.350 = 0.891
In general, the relationship between ∆𝑒 and ∆𝐻 is given by
∆e 1.891
∆H 19.0
=
157
Lecture Note on GTE-I
158
Lecture Note on GTE-I
INTRODUCTION
Shearing Strength’ of a soil is perhaps the most important of its engineering properties. This is
because all stability analyses in the field of geotechnical engineering, whether they relate to
foundation, slopes of cuts or earth dams, involve a basic knowledge of this engineering property
of the soil. ‘Shearing strength’ or merely ‘Shear strength’ may be defined as the resistance to
shearing stresses and a consequent tendency for shear deformation.
Shearing strength of a soil is the most difficult to comprehend in view of the multitude of factors
known to affect it. A lot of maturity and skill may be required on the part of the engineer in
interpreting the results of the laboratory tests for application to the conditions in the field.
Basically speaking, a soil derives its shearing strength from the following:
(1) Resistance due to the interlocking of particles.
(2) Frictional resistance between the individual soil grains, which may be sliding friction, rolling
friction, or both.
(3) Adhesion between soil particles or ‘cohesion’.
Granular soils of sands may derive their shear strength from the first two sources, while cohesive
soils or clays may derive their shear strength from the second and third sources.
Highly plastic clays, however, may exhibit the third source alone for their shearing strength.
Most natural soil deposits are partly cohesive and partly granular and as such, may fall into the
second of the three categories just mentioned, from the point of view of shearing strength. The
shear strength of a soil cannot be tabulated in codes of practice since a soil can significantly
exhibit different shear strengths under different field and engineering conditions.
In granular or cohesionless soil masses, the resistance to sliding on any plane through the point
within the mass is similar to that discussed in the previous sub-section; the friction angle in this
case is called the ‘angle of internal friction’. However, the frictional resistance in granular soil
masses is rather more complex than that between solid bodies, since the nature of the resistance
is partly sliding friction and partly rolling friction. Further, a phenomenon known as
‘interlocking’ is also supposed to contribute to the shearing resistance of such soil masses, as part
of the frictional resistance.
The angle of internal friction, which is a limiting angle of obliquity and hence the primary
criterion for slip or failure to occur on a certain plane, varies appreciably for a given sand with
the density index, since the degree of interlocking is known to be directly dependent upon the
density. This angle also varies somewhat with the normal stress. However, the angle of internal
friction is mostly considered constant, since it is almost so for a given sand at a given density.
159
Lecture Note on GTE-I
Since failure or slip within a soil mass cannot be restricted to any specific plane, it is necessary
to understand the relationships that exist between the stresses on different planes passing through
a point, as a prerequisite for further consideration of shearing strength of soils.
11.2 PRINCIPAL PLANES AND PRINCIPAL STRESSES—MOHR’S CIRCLE
At a point in a stressed material, every plane will be subjected, in general, to a normal or direct
stress and a shearing stress. In the field of geotechnical engineering, compressive direct stresses
are usually considered positive, while tensile stresses are considered negative. A ‘Principal
plane’ is defined as a plane on which the stress is wholly normal, or one which does not carry
shearing stress. From mechanics, it is known that there exists three principal planes at any point
in a stressed material. The normal stresses acting on these principal planes are known as the
‘principal stresses’. The three principal planes are to be mutually perpendicular. In the order of
decreasing magnitude the principal stresses are designated the ‘major principal stress’, the
‘intermediate principal stress’ and the ‘minor principal stress’, the corresponding principal planes
being designated exactly in the same manner. It can be engineering by two-dimensional analysis,
the intermediate principal stress being commonly ignored.
Let us consider an element of soil whose sides are chosen as the principal planes, the major and
the minor, as shown in Fig. 11.1(a):
160
Lecture Note on GTE-I
=𝜎1+𝜎3 + ( ) 𝑐𝑜𝑠2𝜃-----------
𝜎1−𝜎3
2 2
Eq.11.1
𝜏 = 𝑠𝑖𝑛2𝜃--------------
𝜎1−𝜎3
2
Eq.11.2
𝜃
Fig. 11.2: Mohr’s circle for the stress conditions illustrated in Fig. 11.1
intermediate principal plane may be expressed in terms of 𝜎1 𝜎3 , and θ. Otto Mohr (1882)
Thus it may be noted that the normal and shearing stresses on any plane which is normal to the
represented these results graphical in a circle diagram, which is called Mohr’s circle. Normal
stresses are represented as abscissae and shear stresses as ordinates. If the coordinates, σ ɵ and τɵ
in Fig.11.2. This circle has its centre on the axis and cuts it at values 𝜎3 and 𝜎1. This circle is
represented by Eqs.1 and 2 are plotted for all possible values of θ, the locus is a circle as shown
𝜎𝜃 = 𝑀𝐺 = 𝑀𝐹 + 𝐹𝐺
the figure, the co-ordiantes of the point C, are established as follows:
𝜏𝜃 𝜎1 + 𝜎3
+(
2
= 𝐶𝐺 = 𝑠𝑖𝑛2𝜃
𝜎1−𝜎3
2
And
161
Lecture Note on GTE-I
𝜎1 − 𝜎3
) 𝑐𝑜𝑠2𝜃
2
162
Lecture Note on GTE-I
These are the same as in Eqs. 11.1 and 11.2, which prove our statements. In the special case
where the major and minor principal planes are vertical and horizontal respectively, or vice-
versa, the origin of planes will be D or E, as the case may be. In other words, it will lie on the σ-
axis. A few important basic facts and relationships may be directly obtained from the Mohr’s
circle:
1. The only planes free from shear are the given sides of the element which are the principal
planes. The stresses on these are the greatest and smallest normal stresses.
2. The maximum or principal shearing stress is equal to the radius of the Mohr’s circle, and
𝜏𝑚𝑎 =
it occurs on planes inclined at 45° to the principal planes.
𝜎1−𝜎3
2
------------------- Eq.11.3
𝑥
3. The normal stresses on planes of maximum shear are equal to each other and is equal to
half the sum of the principal stresses.
𝜎=
(𝜎1+𝜎3)
2
---------------- Eq.11.4
𝑐
4. Shearing stresses on planes at right angles to each other are numerically equal and are of
an opposite sign. These are called conjugate shearing stresses.
on which it is 𝜎𝜃 �as 𝜎′
𝜎𝜃 + 𝜎 �
′
= 𝜎1 + 𝜎3--------------
�
Of the two stresses 𝜎𝜃 and �𝜎′ the one which makes the smaller angle with σ1 is the greater of the
Eq.11.5
�
�
𝜎𝑟 = √𝜎2 + 𝜏2 -----------------
𝜃 𝜃
Eq.11.6
𝛽 = 𝑡𝑎𝑛−1 ( 𝜃)----------------
𝜏
and has an obliquity, β, which is equal to
Eq.11.7
𝜎𝜃
7. Stresses on conjugate planes, that is, planes which are equally inclined in different
directions with respect to a principal plane are equal. (This is indicated by the co-ordinates of C
and C1 in Fig. 11.2).
8. When the principal stresses are equal to each other, the radius of the Mohr’s circle
becomes zero, which means that shear stresses vanish on all planes. Such a point is called an
isotropic point.
9. The maximum angle of obliquity, βm, occurs on a plane inclined at
𝜃𝑐 = 450+𝛽𝑚
2
--------------- Eq.11.8
𝑟
This may be obtained by drawing a line which passes through the origin and is tangential to the
Mohr’s circle. The co-ordinates of the point of tangency are the stresses on the plane of
maximum obliquity; the shear stress on this plane is obviously less than the principal or
maximum shear stress. On the plane of principal shear the obliquity is slightly smaller than βm.
163
Lecture Note on GTE-I
maximum shear, since the criterion of slip is limiting obliquity. When 𝛽𝑚 approaches and equals
It is the plane of maximum obliquity which is most liable to failure and not the plane of
the angle of internal friction, φ, of the soil, failure will become incipient. Mohr’s circle affords
an easy means of obtaining all important relationships. The following are a few such
relationships
𝑠𝑖𝑛𝛽 = ( )-----------------
𝜎1−𝜎3
Eq.11.9
𝑚 𝜎1+𝜎3
= ( )----------------
𝜎1 1+𝑠𝑖𝑛𝛽𝑚
Eq.11.10
𝜎3 1−𝑠𝑖𝑛𝛽𝑚
In case the normal and shearing stresses on two mutually perpendicular planes are known, the
principal planes and principal stresses may be determined with the aid of the Mohr’s circle
diagram, as shown in Fig. 11.3(a). The shearing stresses on two mutually perpendicular planes
are equal in magnitude by the principle of complementary shear
stresses 𝜎𝑥 and 𝜎𝑦 on mutually perpendicular planes and shear stresses 𝜏𝑥𝑦on these planes,
Figure 11.3 (a) shows an element subjected to a general two-dimensional stress system, normal
as
164
Lecture Note on GTE-I
indicated. Fig. 11.3 (b) shows the corresponding Mohr’s circle, the construction of which is
obvious.
From a consideration of the equilibrium of a portion of the element, the normal and shearing
stress components, 𝜎𝜃 and 𝜏𝜃, respectively, on a plane inclined at an angle θ, measured counter-
𝜎𝑥+𝜎𝑦 𝜎𝑥−𝜎𝑦
𝜎 = +( ) 𝑐𝑜𝑠2𝜃 + 𝜏 𝑠𝑖𝑛2𝜃---------
clockwise with respect to the plane on which σx acts, may be obtained as follows:
𝜃 𝑥𝑦
Eq.11.11
2 2
𝜎𝑥−𝜎𝑦
𝜏 = 𝑠𝑖𝑛2𝜃 𝑥 𝑐𝑜𝑠2𝜃---------------
𝜃 − 𝜏
Eq.11.12
𝑦
2
Squaring and adding these Eqs, we obtain
2
[𝜎𝜃 𝜎𝑥+𝜎𝑦 𝜎𝑥−𝜎
−
2 + 𝜏2 = + 𝜏2 ----------
)
𝑦 2
] (
Eq.11.13
2 � 𝑥
, 0 and radius √(
� 𝜎
𝑦𝑥+𝜎𝑦 𝜎𝑥−𝜎𝑦
2
2
) + 𝜏𝑥𝑦
This represents a circle with centre
2 2
Once the Mohr’s circle is constructed, the principal stresses σ1 and σ3, and the orientation of the
principal planes may be obtained from the diagram. The shearing stress is to be plotted upward
or downward according as it is positive or negative. It is common to take a shear stress which
tends to rotate the element counter-clockwise, positive.
It may be noted that the same Mohr’s circle and hence the same principal stresses are obtained,
irrespective of how the shear stresses are plotted. (The centre of the Mohr’s circle, C, is the mid-
𝜏 = √(
1 𝜎𝑥−𝜎𝑦
2
2
) + 4𝜏𝑥𝑦
--------------- Eq.11.17
𝑚𝑎𝑥 2
2
Invariably, the vertical stress will be the major principal stress and the horizontal one the minor
principal stress in geotechnical engineering situations.
11.2.1 Mohr’s Strength Theory
We have seen that the shearing stress may be expressed as τ = σ tan β on any plane, where β is
the angle of obliquity. If the obliquity angle is the maximum or has limiting value φ, the shearing
165
Lecture Note on GTE-I
stress is also at its limiting value and it is called the shearing strength, s. For a cohesionless soil
the shearing strength may be expressed as:
s = σ tan φ ... Eq.11.18
If the angle of internal friction φ is assumed to be a constant, the shearing strength may be
represented by a pair of straight lines at inclinations of + φ and – φ with the σ-axis and passing
166
Lecture Note on GTE-I
through the origin of the Mohr’s circle diagram. A line of this type is called a Mohr envelope.
The Mohr envelopes for a cohesionless soil, as shown in Fig. 11.4, are the straight lines OA and
OA′.
𝑠 = 𝜎𝑓tan∅---------------------------------------------
we add the subscript f to σ:
Eq.11.19
It is possible to express the strength in terms of normal stress on any plane, with the aid of the
𝑠 = 𝜎 𝑡𝑎𝑛∅ = 𝑐𝑜𝑠∅---------------------------
𝜎1−𝜎3
Eq.11.20
Eq.11.21
𝑓 2
167
Lecture Note on GTE-I
The primary assumptions in the Mohr’s strength theory are that the intermediate principal stress
has no influence on the strength and that the strength is dependent only upon the normal stress on
the plane of maximum obliquity. However, the shearing strength, in fact, does depend to a small
extent upon the intermediate principal stress, density speed of application of shear, and so on.
But the Mohr theory explains satisfactorily the strength concept in soils and hence is in vogue. It
may also be noted that the Mohr envelope will not be a straight line but is actually slightly
curved since the angle of internal friction is known to decrease slightly with increase in stress.
11.2.3 Mohr-Coulomb Theory
The Mohr-Coulomb theory of shearing strength of a soil, first propounded by Coulomb (1976)
and later generalised by Mohr, is the most commonly used concept. The functional relationship
between the normal stress on any plane and the shearing strength available on that plane was
assumed to be linear by Coulomb; thus the following is usually known as Coulomb’s law:
s = c + σ tan φ ------------- Eq. 11.22
where c and φ are empirical parameters, known as the ‘apparent cohesion’ and ‘angle of shearing
resistance’ (or angle of internal friction), respectively. These are better visualised as ‘parameters’
and not as absolute properties of a soil since they are known to vary with water content,
conditions of testing such as speed of shear and drainage conditions, and a number of other
factors besides the type of soil.
Coulomb’s law is merely a mathematical equation of the failure envelope shown in Figure 11.5
(a); Mohr’s generalisation of the failure envelope as a curve which becomes flatter with
increasing normal stress is shown in Fig. 11.5 (b).
168
Lecture Note on GTE-I
In a different way, it can be said that the Mohr’s circle of stress relating to a given stress
condition would represent, incipient failure condition if it just touches or is tangent to the
strength or failure envelope (circle I); otherwise, it would wholly lie below the envelopes as
shown in circle II, Fig.11.5 (b).
The Coulomb envelope in special cases may take the shapes given in Fig. 11.6 (a) and (b); for a
purely cohesionless or granular soil or a pure sand, it would be as shown in Fig. 11.6 (a) and for
a purely cohesive soil or a pure clay, it would be as shown in Fig. 11.6 (b).
Fig. 11.6: Coulomb envelopes for pure sand and for pure clay
s = c′ + 𝜎𝑓 tan φ′ --------------
Eq. 11.23 may now be modified to read:
Eq. 11.24
Where, c′ and φ′ are called the effective cohesion and effective angle of internal friction,
respectively, since they are based on the effective normal stress on the failure plane.
169
Lecture Note on GTE-I
Collectively, they are called ‘effective stress parameters’, while c and φ of Eq. 24 are called
‘‘total stress parameters’’.
Determination of shearing strength of a soil involves the plotting of failure envelopes and
evaluation of the shear strength parameters for the necessary conditions. The following tests are
available for this purpose:
Laboratory Tests
1. Direct Shear Test
2. Triaxial Compression Test
3. Unconfined Compression Test
4. Laboratory Vane Shear Test
5. Torsion Test
6. Ring Shear Tests
Field Tests
1. Vane Shear Test
2. Penetration Test
The first three tests among the laboratory tests are very commonly used, while the fourth is
gaining popularity owing to its simplicity. The fifth and sixth are mostly used for research
purposes and hence are not dealt with here.
The principle of the field vane test is the same as that of the laboratory vane shear test, except
that the apparatus is bigger in size for convenience of field use. The penetration test involves the
measurement of resistance of a soil to penetration of a cone or a cylinder, as an indication of the
shearing strength. This procedure is indirect and rather empirical in nature although correlations
are possible. The field tests are also not considered here. The details of the test procedures are
available in the relevant I.S. codes or any book on laboratory testing, such as Lambe (1951).
11.4.1 Direct Shear Test
The direct shear device, also called the ‘shear box apparatus’, essentially consists of a brass box,
split horizontally at mid-height of the soil specimen, as shown schematically in Figure 11.7.
The soil is gripped in perforated metal grilles, behind which porous discs can be placed if
required to allow the specimen to drain. For undrained tests, metal plates and solid metal grilles
may be used. The usual plan size of the specimen is 60 mm square; but a larger size such as 300
mm square or even more, is employed for testing larger size granular material such as gravel.
The minimum thickness or height of the specimen is 20 mm. After the sample to be tested is
placed in the apparatus or shear box, a normal load which is vertical is applied to the top of the
sample by means of a loading yoke and weights. Since the shear plane is predetermined as the
horizontal plane, this becomes the normal stress on the failure plane, which is kept constant
throughout the test. A shearing force is applied to the upper-half of the box, which is zero
initially and is increased until the specimen fails.
170
Lecture Note on GTE-I
171
Lecture Note on GTE-I
carefully machined, balanced and calibrated. The deflection of the annular ring is measured with
the aid of a dial gauge set inside the ring, the causative force being got for any displacement by
means of the calibration chart supplied by the manufacturer.
The shear displacement is measured again with the aid of another dial gauge attached to the side
of the box
In both cases, a dial gauge attached to the plunger, through which the normal load is applied, will
enable one to determine the changes in the thickness of the soil sample which will help in the
computation of volume changes of the sample, if any. The strain-controlled type is very widely
used. The strain is taken as the ratio of the shear displacement to the thickness of the sample. The
proving ring readings may be taken at fixed displacements or even at fixed intervals of time as
the rate of strain is made constant by an electric motor. A sudden drop in the proving ring
reading or a leveling-off in successive readings indicates shear failure of the soil specimen.
The shear strain may be plotted against the shear stress; it may be plotted versus the ratio of the
shearing stress on normal stress; and it may also be plotted versus volume change.
Each plot may yield information useful in one way or the other. The stresses may be obtained
from the forces by dividing them by the area of cross-section of the sample. The stress-
conditions on the failure plane and the corresponding Mohr’s circle for direct shear test are
shown in Fig. 11.8 (a) and (b) respectively.
Fig.11.8: Mohr’s circle representation of stress condition for direct shear test results
The failure plane is predetermined as the horizontal plane here. Several specimens are tested
under different normal loads and the results plotted to obtain failure envelopes.
The direct shear test is a relatively simple test. Quick drainage, i.e., quick dissipation of pore
pressures is possible since the thickness of the specimen is small. However, the test suffers from
the following inherent disadvantages, which limit its application.
1. The stress conditions are complex primarily because of the non-uniform distribution of normal
and shear stresses on the plane.
2. There is virtually no control of the drainage of the soil specimen as the water content of a
saturated soil changes rapidly with stress.
3. The area of the sliding surface at failure will be less than the original area of the soil specimen
and strictly speaking, this should be accounted for.
4. The ridges of the metal gratings embedded on the top and bottom of the specimen, causes
distortion of the specimen to some degree.
172
Lecture Note on GTE-I
5. The effect of lateral restraint by the side walls of the shear box is likely to affect the results.
6. The failure plane is predetermined and this may not be the weakest plane. In fact, this is the
most important limitation of the direct shear test.
11.4.2 Triaxial Shear Test
The triaxial shear test is one of the most reliable methods available for determining shear strength
parameters. It is used widely for research and conventional testing. A diagram of the triaxial test
layout is shown in Figure 11.9.
In this test, a soil specimen about 36 mm in diameter and 76 mm (3 in.) long generally is used. The
specimen is encased by a thin rubber membrane and placed inside a plastic cylindrical chamber
that usually is filled with water or glycerine. The specimen is subjected to a confining pressure
by compression of the fluid in the chamber. (Note: Air is sometimes used as a compression
medium.) To cause shear failure in the specimen, one must apply axial stress (sometimes called
deviator stress) through a vertical loading ram.
This stress can be applied in one of two ways:
1. Application of dead weights or hydraulic pressure in equal increments until the specimen fails.
(Axial deformation of the specimen resulting from the load applied through the ram is measured
by a dial gauge.)
2. Application of axial deformation at a constant rate by means of a geared or hydraulic loading
press. This is a strain-controlled test.
The axial load applied by the loading ram corresponding to a given axial deformation is measured by a
proving ring or load cell attached to the ram.
173
Lecture Note on GTE-I
Test Procedure
The essential steps in the conduct of the test are as follows:
(i) A saturated porous stone is placed on the pedestal and the cylindrical soil specimen is placed
on it.
(ii) The specimen is enveloped by a rubber membrane to isolate it from the water with which the
cell is to be filled later; it is sealed with the pedestal and top cap by rubber ‘‘O’’ rings.
(iii) The cell is filled with water and pressure is applied to the water, which in turn is transmitted
to the soil specimen all-round and at top. This pressure is called ‘cell pressure’, ‘chamber
pressure’ or ‘confining pressure’.
(iv) Additional axial stress is applied while keeping the cell pressure constant. This introduces
shearing stresses on all planes except the horizontal and vertical planes, on which the major,
minor and intermediate principal stresses act, the last two being equal to the cell pressure on
account of axial symmetry.
(v) The additional axial stress is continuously increased until failure of the specimen occurs.
(What constitutes failure is often a question of definition and may be different for different kinds
of soils. This aspect would be discussed later on).
A number of observations may be made during a Triaxial compression test regarding the
physical changes occurring in the soil specimen:
(a) As the cell pressure is applied, pore water pressure develops in the specimen, which can be
measured with the help of a pore pressure measuring apparatus, such as Bishop’s pore pressure
device (Bishop, 1960), connected to the pore pressure line, after closing the valve of the drainage
line.
(b) If the pore pressure is to be dissipated, the pore water line is closed, the drainage line opened
and connected to a burette. The volume decrease of the specimen due to consolidation is
indicated by the water drained into the burette.
(c) The axial strain associated with the application of additional axial stress can be measured by
means of a dial gauge, set to record the downward movement of the loading piston.
(d) Upon application of the additional axial stress, some pore pressure develops. It may be
measured with the pore pressure device, after the drainage line is closed. On the other hand, if it
is desired that any pore pressure developed be allowed to be dissipated, the pore water line is
closed and the drainage line opened as stated previously.
(e) The cell pressure is measured and kept constant during the course of the test.
(f) The additional axial stress applied is also measured with the aid of a proving ring and dial
gauge.
Thus the entire triaxial test may be visualised in two important stages:
(i) The specimen is placed in the triaxial cell and cell pressure is applied during the first stage.
(ii) The additional axial stress is applied and is continuously increased to cause a shear failure,
the potential failure plane being that with maximum obliquity during the second stage.
174
Lecture Note on GTE-I
11.4.3 Area Correction for the Determination of Additional Axial Stress or Deviatoric Stress
The additional axial load applied at any stage of the test can be determined from the proving ring
reading. During the application of the load, the specimen undergoes axial compression and
horizontal expansion to some extent. Little error is expected to creep in if the volume is supposed
to remain constant, although the area of cross-section varies as axial strain increases.
The assumption is perfectly valid if the test is conducted under undrained conditions, but, for drained
conditions, the exact relationship is somewhat different.
If A0, h0 and V0 are the initial area of cross-section, height and volume of the soil specimen respectively,
and if A, h, and V are the corresponding values at any stage of the test, the corresponding
Hence 𝐴 =
𝑉0+∆𝑉
Eq.11.25
------------------------------- Eq.11.26
ℎ0−∆ℎ
But, for axial 𝑉 ∆𝑉
compression, Δh∆𝑉
is known to be negative.
0(1+ ) 𝐴0(1+
𝐴0+∆𝑉
=
𝑉0 𝑉0
=
𝑉 ---------------------------- Eq.11.27
1−𝜖𝑎
=∆ℎ
ℎ0(1− )
ℎ0−∆ℎ ℎ0
an undrained test, 𝐴 =
since the axial strain, εa = Δh/h0------------------------- Eq.11.28
For
𝐴0 ------------------------- Eq.11.29
1−𝜖𝑎
since ΔV = 0.
1
This is called the ‘Area correction’ and
1−𝜖𝑎
is the correction factor.
𝐴 =𝐴0
A more accurate expression for the corrected area is given by
=
𝑉0+∆𝑉
1−𝜖
------------------------ Eq.11.30
𝑎
ℎ0−∆ℎ
Once the corrected area is determined, the additional axial stress or the deviator stress, Δσ, is
The cell pressure or the confining pressure, σc, itself being the minor principal stress, σ3, this is
𝜎1 = 𝜎3 + ∆𝜎------------------------
constant for one test; however, the major principal stress, σ1, goes on increasing until failure.
Eq.11.31
11.4.4 Mohr’s Circle for Triaxial Test
The stress conditions in a triaxial test may be represented by a Mohr’s circle, at any stage of the
test, as well as at failure, as shown in Figure 11.10:
175
Lecture Note on GTE-I
176
Lecture Note on GTE-I
The cell pressure, 𝜎𝑐 which is also the minor principal stress is constant and 𝜎11, 𝜎12, 𝜎13, 𝜎1𝑓
are the major principal stresses at different stages of loading and at failure. The Mohr’s circle at
failure will be tangential to the Mohr-Coulomb strength envelope, while those at intermediate
stages will be lying wholly below it. The Mohr’s circle at failure for one particular value of cell
pressure will be as shown in Fig. 11.11.
∴ α = 45° + φ/2------------------
From ΔDCG, 2α = 90° + φ
Eq.11.32
Again from ΔDCG
𝜎1−𝜎3
𝑠𝑖𝑛∅ = = =
𝐷𝐶 𝐷𝐶
( )
2 𝜎 + -------------------------------
𝐺𝐶 𝑐𝑐𝑜𝑡∅+ 1
Eq.11.33
𝜎3
𝐺𝑀+𝑀𝐶 2
𝜎1 − 𝜎3 = 2𝑐𝑐𝑜𝑡∅ + (𝜎1 + 𝜎3)𝑠𝑖𝑛∅------------------------ Eq.11.34
= σ3 tan2 (45° + φ/2) + 2c tan(450 + )------------------------
∅
𝜎1
Or, Eq.11.35
2
177
Lecture Note on GTE-I
178
Lecture Note on GTE-I
The above equations define the relationship between the principal stresses at failure. This state of
stress is defined as ‘Plastic equilibrium condition’, when failure is imminent.
sets are necessary to evaluate the parameters c and ∅ conventionally, three or more such sets are
From one test, a set of σ1 and σ3 is known; however, it can be seen from that at least two such
(8) Special tests such as extension tests are also possible to be conducted with the triaxial testing
apparatus.
(9) It provides an ingenious and a symmetrical three-dimensional stress system better suited to
simulate field conditions.
In the CD test, the saturated specimen first is subjected to an all around confining pressure, 𝜎3, by
compression of the chamber fluid. As confining pressure is applied, the pore water pressure of
the specimen increases by uc (if drainage is prevented). This increase in the pore water pressure
𝐵 = ------------------
𝑢𝑐
can be expressed as a non dimensional parameter in the form
Eq.11.39
𝜎3
Where, B= Skempton’s pore pressure parameter (Skempton, 1954).
180
Lecture Note on GTE-I
For saturated soft soils, B is approximately equal to 1; however, for saturated stiff soils, the
magnitude of B can be less than 1. Black and Lee (1973) gave the theoretical values of B for
various soils at complete saturation.
Now, if the connection to drainage is opened, dissipation of the excess pore water pressure, and
thus consolidation, will occur. With time, uc will become equal to 0. In saturated soil, the change
in the volume of the specimen (∆Vc) that takes place during consolidation can be obtained from
Next, the deviator stress, ∆𝜎𝑑 on the specimen is increased very slowly (Figure 11.13b). The
the volume of pore water drained (Figure 11.13a).
drainage connection is kept open, and the slow rate of deviator stress application allows
complete dissipation of any pore water pressure that developed as a result (∆𝑢𝑑=0)
A typical plot of the variation of deviator stress against strain in loose sand and normally
over consolidated clay. The volume change, ∆ Vd, of specimens that occurs because of the
consolidated clay is shown in Figure 11.13b. Figure 1c shows a similar plot for dense sand and
application of deviator stress in various soils is also shown in Figures 11.13d and e.
Fig. 11.13: Consolidated-drained triaxial test: (a) volume change of specimen caused by
chamber-confining pressure; (b) plot of deviator stress against strain in the vertical direction for
loose sand and normally consolidated clay; (c) plot of deviator stress against strain in the vertical
direction for dense sand and over consolidated clay; (d) volume change in loose sand and
normally consolidated clay during deviator stress application; (e) volume change in dense sand
and over consolidated clay during deviator stress application
181
Lecture Note on GTE-I
Fig.11.14: Effective stress failure envelope from drained tests on sand and normally consolidated
clay
The coordinates of the point of tangency of the failure envelope with a Mohr’s circle (that is,
point A) give the stresses (normal and shear) on the failure plane of that test specimen.
𝑠𝑖𝑛∅ = ----------------------
′
For normally consolidated clay, referring to Figure 2
′ 𝑂𝐴
𝑂𝑂′
Eq.11.40
OR 𝑠𝑖𝑛∅′=𝜎
′ ′
𝜎1 3
−𝜎
′
+𝜎1
--------------------------- Eq.11.41
1 3 𝜎′−𝜎′
Hence ∅′ = 𝑠𝑖𝑛−1𝜎
( ′ 1+𝜎
3
)----------------------
1 Eq.11.42
1 3
Also, the failure plane will be inclined at an angle of 𝜃 = 45 + to the major principal plane, as
∅
2
shown in Figure 11.14.
𝜏𝑓 = 𝑐′ + 𝜎′𝑡𝑎𝑛∅′--------------------
cohesion intercept, and the shear strength equation for this branch can be written as
Eq.11.43
182
Lecture Note on GTE-I
Fig. 11.15: Effective stress failure envelope for over consolidated clay
A consolidated-drained triaxial test on a clayey soil may take several days to complete. This
amount of time is required because deviator stress must be applied very slowly to ensure full
drainage from the soil specimen. For this reason, the CD type of triaxial test is uncommon.
saturated soil specimen is first consolidated by an all-around chamber fluid pressure, 𝜎3, that
The consolidated-undrained test is the most common type of triaxial test. In this test, the
183
Lecture Note on GTE-I
deviator stress against axial strain for loose sand and normally consolidated clay; (e) deviator
stress against axial strain for dense sand and over consolidated clay; (f) variation of pore water
pressure with axial strain for loose sand and normally consolidated clay; (g) variation of pore
water pressure with axial strain for dense sand and over consolidated clay.
deviator stress, ∆𝜎𝑑, on the specimen is increased to cause shear failure (Figure 11.16c). During
After the pore water pressure generated by the application of confining pressure is dissipated, the
permitted, the pore water pressure, ∆𝑢𝑑 , will increase. During the test, simultaneous
this phase of the test, the drainage line from the specimen is kept closed. Because drainage is not
measurements of ∆𝜎𝑑, and ∆𝑢𝑑 are made. The increase in the pore water pressure, ∆𝑢𝑑can be
𝐴 = 𝑑----------------------
∆𝑢
expressed in a nondimensional form as
Eq.11.44
∆𝜎𝑑
where 𝐴̅ = Skempton’s pore pressure parameter (Skempton, 1954).
The general patterns of variation of ∆𝜎𝑑 and ∆𝑢𝑑with axial strain for sand and clay soils
are shown in Figures 4 d through 4g. In loose sand and normally consolidated clay, the pore
water pressure increases with strain. In dense sand and over consolidated clay, the pore water
pressure increases with strain to a certain limit, beyond which it decreases and becomes
negative (with respect to the atmospheric pressure). This decrease is because of a tendency of
the soil to dilate. Unlike the consolidated-drained test, the total and effective principal stresses
are not the same in the consolidated-undrained test. Because the pore water pressure at failure
is measured in this test, the principal stresses may be analyzed as follows:
• Major principal stress at failure (total):
• Major principal stress at failure (effective):
• Minor principal stress at failure (total):
• Minor principal stress at failure (effective):
𝜎1 − 𝜎3 = 𝜎′ − 𝜎 ′ --------------------------
The preceding derivations show that
1 3
Eq.11.45
Tests on several similar specimens with varying confining pressures may be conducted to
determine the shear strength parameters. Figure 11.17 shows the total and effective stress
Fig. 11.17: Total and effective stress failure envelopes for consolidated undrained triaxial tests
184
Lecture Note on GTE-I
Mohr’s circles at failure obtained from consolidated-undrained triaxial tests in sand and normally
consolidated clay. Note that A and B are two total stress Mohr’s circles obtained from two tests.
C and D are the effective stress Mohr’s circles corresponding to total stress circles A and B,
respectively. The diameters of circles A and C are the same; similarly, the diameters of circles B
and D are the same. In Figure 11.17, the total stress failure envelope can be obtained by drawing
a line that touches all the total stress Mohr’s circles. For sand and normally consolidated clays,
this will be approximately a straight line passing through the origin and may be expressed by the
𝜏𝑓 = 𝜎𝑡𝑎𝑛∅--------------------
equation
∅ =the angle that the total stress failure envelope makes with the normal stress axis, also known
as the consolidated-undrained angle of shearing resistance
∅ = 𝑠𝑖𝑛−1 ( )--------------------
𝜎1−𝜎3
For sand and normally consolidated clay, we can write
Eq.11.47
𝜎1+𝜎3
∅′
= 𝑠𝑖𝑛−1 (
𝜎1−𝜎3
And
𝜎1+𝜎3−2(∆𝑢
)-------------------- Eq.11.48
𝑑)
Fig. 11.18: Total stress failure envelope obtained from consolidated-undrained tests in over
consolidated clay
Again referring to Fig.11.17, we see that the failure envelope that is tangent to all the effective
𝜏𝑓 = 𝜎′𝑡𝑎𝑛∅′-----------------------
stress Mohr’s circles can be represented by the equation
Eq.11.48
which is the same as that obtained from consolidated-drained tests.
In over consolidated clays, the total stress failure envelope obtained from consolidated undrained
𝜏𝑓 = 𝜎𝑡𝑎𝑛∅1------------------------
tests will take the shape shown in Fig.11.18. The straight line is represented by the equation
Eq.11.49
and the straight line follows the relationship given by the above equation. The effective stress
failure envelope drawn from the effective stress Mohr’s circles will be similar to that shown in
Figure 11.15.
Consolidated-drained tests on clay soils take considerable time. For this reason, consolidated-
undrained tests can be conducted on such soils with pore pressure measurements to obtain the
drained shear strength parameters. Because drainage is not allowed in these tests during the
185
Lecture Note on GTE-I
application of deviator stress, they can be performed quickly. Skempton’s pore water pressure
∆𝑢defined
̅
𝐴̅ = 𝐴̅̅ =
𝑑𝑓
parameter was as follows. At failure, the parameter can be written as
𝑓
----------------- Eq.11.50
∆𝜎
�𝑓
�
11.5.3 Unconsolidated-Undrained Triaxial Test
application of chamber pressure 𝜎3. The test specimen is sheared to failure by the application of
In unconsolidated-undrained tests, drainage from the soil specimen is not permitted during the
deviator stress, ∆𝜎𝑑, and drainage is prevented. Because drainage is not allowed at any stage,
the test can be performed quickly. Because of the application of chamber confining pressure, 𝜎3,
the pore water pressure in the soil specimen will increase by uc. A further increase in the pore
water pressure (∆𝑢𝑑) will occur because of the deviator stress application. Hence, the total pore
𝑢 = 𝑢𝑐 + ∆𝑢𝑑--------------------------
water pressure u in the specimen at any stage of deviator stress application can be given as
𝑢 = 𝐵𝜎3 +𝐴̅∆𝜎𝑑----------------------
Now substituting above relation, we obtain
Eq.11.52
concept for cohesive soils if the soil is fully saturated. The added axial stress at failure ∆𝜎𝑑𝑓 is
This test usually is conducted on clay specimens and depends on a very important strength
practically the same regardless of the chamber confining pressure. This property is shown in
𝜏𝑓 = 𝑐 = 𝑐𝑢---------------------------- Eq.11.53
Note that the ∅ = 0 concept is applicable to only saturated clays and silts.
Where, cu is the undrained shear strength and is equal to the radius of the Mohr’s circles.
Fig. 11.19: Mohr’c circle for total and effective stress and failure envelope ∅ = 0
The reason for obtaining the same added axial stress ∆𝜎𝑑𝑓 regardless of the confining
pressure can be explained as follows. If clay specimen (No. I) is consolidated at a chamber
186
Lecture Note on GTE-I
pressure 𝜎3and
187
Lecture Note on GTE-I
then sheared to failure without drainage, the total stress conditions at failure can be represented
equal to ∆𝑢𝑑𝑓. Thus, the major and minor principal effective stresses at failure are, respectively,
by the Mohr’s circle P in Fig.11.19. The pore pressure developed in the specimen at failure is
Any value of 𝜎3 could have been chosen for testing the specimen. In any case, the deviator stress
Note that the diameters of circles P and Q are the same.
∆𝜎𝑑𝑓 to cause failure would have been the same as long as the soil was fully saturated and fully
undrained during both stages of the test.
188
Lecture Note on GTE-I
2
Eq.11.56
1
The two unknowns–c and ϕ–cannot be solved since any number of unconfined compression tests
would give only one value for σ 1. Therefore, the unconfined compression test is mostly found
useful in the determination of the shearing strength of saturated clays for which ϕ is negligible or
𝜎1 = 𝑄𝑢 = 2𝑐-----------------------
zero, under undrained conditions. In such a case, the above equation reduces to
Eq.11.57
189
Lecture Note on GTE-I
190
Lecture Note on GTE-I
𝑇 = 𝜋𝐷2𝑠 [ + ]----------------
𝐻 𝐷
If only the bottom end partake in the shearing the above equation takes as
2 12
Eq.11.62
Regarding the shearing stress distribution on the soil cylinder, it is assumed uniform on the
cylindrical surface but it is triangular over the shear end faces, varying from zero at the axis of
the vane device, to maximum at the edge, as shown in Fig. 11.23.
Fig. 11.23: Shearing Distribution on the Sides and Faces of Soil Cylinder in the Vane Shear Test
The vane shear test is particularly suited for soft clays and sensitive clays for which suitable
cylindrical specimens cannot be easily prepared.
191
Lecture Note on GTE-I
interlocking. The Mohr envelopes merely show large ordinates and steeper slopes for dense soils
than for loose ones.
The angle of internal friction is a measure of the resistance of the soil to sliding along a plane.
This varies with the density of packing, characterised by density index, particle shape and
roughness and particle size distribution. Its value increases with density index, with the
angularity and roughness of particles and also with better gradation. This is influenced to some
extent by the normal pressure on the plane of shear and also the rate of application of shear.
The ‘angle of repose’ is the angle to the horizontal at which a heap of dry sand, poured freely
from a small height, will stand without support. It is approximately the same as the angle of
friction in the loose state.
Some clean sands exhibit slight cohesion under certain conditions of moisture content, owing to
capillary tension in the water contained in the voids. Since this is small and may disappear with
change in water content, it should not be relied upon for shear strength. On the other hand, even
small percentages of silt and clay in sand give it cohesive properties which may be sufficiently
large so as to merit consideration.
Unless drainage is deliberately prevented, a shear test on a sand will be a drained one as the high
value of permeability makes consolidation and drainage virtually instantaneous. A sand can be
tested either in the dry or in the saturated condition. If it is dry, there will be no pore water
pressures and if it is saturated, the pore water pressure will be zero due to quick drainage. In
either case, the inter granular pressure will be equal to the applied stress. However, there may be
certain situations in which significant pore pressures are developed, at least temporarily, in
sands. For example, during earth-quakes, heavy blasting and operation of vibratory equipment
instantaneous pore pressures are likely to develop due to large shocks or dynamic loads. These
may lead to the phenomenon of ‘liquefaction’ or sudden and total loss of shearing strength,
which is a grave situation of lack of stability.
192
Lecture Note on GTE-I
193
Lecture Note on GTE-I
through the origin and inclined at the angle of internal friction to the normal stress axis. It is
shown in Fig. 6 (b).
194
Lecture Note on GTE-I
axial normal stress is gradually increased keeping the cell pressure constant, until failure occurs.
The value of φ is obtained by plotting the Mohr Circles and the corresponding Mohr’s envelope.
The failure envelope obtained from a series of drained triaxial compression tests on saturated
sand specimens initially at the same density index is approximately a straight line passing
through the origin, as shown in Fig. 11.27.
∅= 21° 48’
tan∅= 4/10 = 0.4
𝜃 = 45° +21048’/2
= 55°54’
195
Lecture Note on GTE-I
Fig.1
196
Lecture Note on GTE-I
Sol:
197
Lecture Note on GTE-I
Sol:
198
Lecture Note on GTE-I
Graphical Sol:
Sol:
199