Lecture 1623481301
Lecture 1623481301
LECTURE NOTE
DYNAMICS OF SOILS AND FOUNDATIONS
SECOND SEMESTER
M.TECH (GTE)
1
Lecture Notes on DSF 2021
DISCLAIMER
These lecture notes are being prepared and printed for the use in training
the students. No commercial use of these notes is permitted and copies
of these will not be offered for sale in any manner. This lecture notes
will provide the basic knowledge about the course. The students are
advised to refer Text and Reference books to gather more knowledge.
The readers are encouraged to provide written feedbacks for
improvement of course materials.
2
Lecture Notes on DSF 2021
Subject Name: DYNAMICS OF SOILS AND FOUNDATIONS MCEGT201
Course Content
Module-I
Fundamentals of vibrations: single, two and multiple degree of freedom systems,
vibration isolation, vibration absorbers, vibration measuring instruments.
Module-II
Wave propagation: elastic continuum medium, semi-infinite elastic continuum
medium, soil behaviour under dynamic loading.
Module-III
Liquefaction of soils: liquefaction mechanism, factors affecting liquefaction, studies by
dynamic tri-axial testing, shake table and blast tests, assessment of liquefaction potential.
Module-IV
Dynamic elastic constants of soil: determination of dynamic elastic constants, various
methods including block resonance tests, cyclic plate load tests, wave propagation tests,
oscillatory shear box test.
Module-V
Theory of Vibration of Foundation: Vertical, sliding, torsional and rocking oscillation of
footing resting on Elastic half space. Oscillation of rigid circular footing supported by an
elastic layer. Introduction of bearing capacity of dynamically loaded shallow foundation.
Reference Books:
Das, B.M., “Fundamentals of Soil Dynamics”, Elsevier, 1983.
Steven Kramer, “Geotechnical Earthquake Engineering”, Pearson, 2008.
Prakash, S., Soil Dynamics, McGraw Hill, 1981.
Kameswara Rao, N.S.V., Vibration analysis and foundation dynamics, Wheeler
Publication Ltd., 1998.
Richart, F.E. Hall J.R and Woods R.D., Vibrations of Soils and Foundations, Prentice Hall
Inc., 1970.
Prakash, S. and Puri, V.K., Foundation for machines: Analysis and Design, John Wiley &
Sons, 1998
COURSE OUTCOME
Students can interpret theory of vibration and resonance phenomenon, dynamic
amplification.
Students can investigate propagation of body waves and surface waves through soil.
Students can predict dynamic bearing capacity and assess liquefaction potential of any site.
Student exposed to different methods for estimation of dynamic soil properties required for
design purpose.
Students apply theory of vibrations to design machine foundation based on dynamic soil
properties and bearing capacity
3
Lecture Notes on DSF 2021
1.0 FUNDAMETALS OF VIBRATION
In order to understand the behaviour of a structure subjected to dynamic load lucidly, one must
study the mechanics of vibrations 'caused by the dynamic load. The pattern of variation of a
dynamic load with respect to time may be either periodic or transient. The periodical motions can
be resolved into sinusoidally varying components e.g. vibrations in the case of reciprocating
machine foundations. Transient vibrations may have very complicated non-periodic time history
e.g. vibrations due to earthquakes and quarry blasts.
A structure subjected to a dynamic load (periodic or transient) may vibrate in one of the
following four ways of deformation or a combination there-of:
(i) Extensional
(ii) Bending
(iii) Shearing
(iv) Torsional
The forms of vibration mainly depend on the mass, stiffness distribution and end conditions of
the system.
To study the response of a vibratory system, in many cases it is satisfactory to reduce it to an
idealized system of lumped parameters. In this regard, the simplest model consists of mass,
spring and dashpot. This chapter is framed to provide the basic concepts and dynamic analysis of
such systems. Actual field problems which can be idealized to mass-spring-dashpot systems,
have also been included.
1.1 Important Definition
Vibrations: If the motion of the body is oscillatory in character, it is called vibration.
Degrees of Freedom: The number of independent co-ordinates which are required to define the
position of a system during vibration, is called degrees of freedom (Fig. 1)
Periodic Motion: If motion repeats itself at regular intervals of time, it is called periodic motion.
Free Vibration: If a system vibrates without an external force, then it is said to undergo free
vibrations. Such vibrations can be caused by setting the system in motion initially and allowing it
to move.
Natural Frequency: This is the property of the system and corresponds to the number of free
oscillations made by the system in unit time.
Forced Vibrations: Vibrations that are developed by externally applied exciting forces are called
forced vibrations. These vibrations occur at the frequency of the externally applied exciting
force.
Forcing Frequency: This refers to the periodicity of the external forces which acts on the system
during forced vibrations. This is also termed as operating frequency.
Frequency Ratio: The ratio of the forcing frequency and natural frequency of the system is
referred as frequency ratio.
Amplitude of Motion: The maximum displacement of a vibrating body from the mean position is
amplitude of motion.
Time Period: Time taken to complete one cycle of vibration is known as time period.
4
Lecture Notes on DSF 2021
Resonance: A system having n degrees of freedom has n natural frequencies. If the frequency of
excitation coincides with anyone of the natural frequencies of the system, the condition of
resonance occurs. The amplitudes of motion are very excessive at resonance.
6
Lecture Notes on DSF 2021
The Eq. (1.1) is plotted as function of time in Fig.3. The various terms of this equation are as
follows:
Z = Displacement of the rotating mass at any time t
A = Displacement amplitude from the mean position, sometimes referred as single amplitude.
The distance 2A represents the peak-to-peak displacement amplitude, sometimes referred to as
double amplitude, and is the quantity most often measured from vibration records.
ω= Circular frequency in radians per unit time. Because the motion repeats itself after 2π radians,
the frequency of oscillation in terms of cycles per unit time will be 𝜔⁄2𝜋. It is denoted by f
θ= Phase angle. It is required to specify the time relationship between two quantities having the
same frequency when their peak values having like sign do not occur simultaneously. In Eq. (1)
the phase angle is a reference to the time origin.
The time period, T is given by
1 2𝜋
𝑇= = ------------------- Eq.1.2
𝑓 𝜔
The velocity and acceleration of motion are obtained from the derivatives of Eq. (1.1)
𝑑𝑍
Velocity = = 𝐴𝜔𝑐𝑜𝑠(𝜔𝑡 − 𝜃)----------------- Eq.1.3
𝑑𝑡
=𝐴𝜔sin(𝜔𝑡 − 𝜃 + 𝜋⁄2)
𝑑2 𝑍
Acceleration = = 𝜔2 𝐴𝑠𝑖𝑛(𝜔𝑡 − 𝜃)--------------- Eq.1.4
𝑑𝑡 2
2
=𝜔 𝐴𝑠𝑖𝑛(𝜔𝑡 − 𝜃 + 𝜋)
Equations (1.3) and (1.4) show that both velocity and acceleration are also harmonic and can be
represented by vectors ωA and 𝜔2 𝐴,which rotate at the same speed as A, i.e. ω rad/unit time.
These, however, lead the displacement and acceleration vectors by 𝜋⁄2and π respectively. In
Fig.4 vector representation of harmonic displacement, velocity and acceleration is presented
considering the displacement as the reference quantity (θ = 0)
7
Lecture Notes on DSF 2021
1.3 VIBRATIONS OF A SINGLE DEGREE FREEDOM SYSTEM
The simplest model to represent a single degree of freedom system consisting of a rigid mass m
supported by a spring and dashpot is shown in Fig. 1.5 a. The motion of the mass m is specified
by one co-ordinate, Z. Damping in this system is represented by the dashpot, and the resulting
damping force is proportional to the velocity. The system is subjected to an external time
dependent force F (t).
𝑉
And 𝐴𝑍 = √𝑍02 + (𝜔0 )2 ---------------------- Eq.1.16
𝑛
The displacement, velocity and acceleration of mass as expressed in above eqs can be
graphically shown as
9
Lecture Notes on DSF 2021
1 𝑔
It can be shown that 𝑓𝑛 = 2𝜋 √𝛿 ----------- Eq.1.18
𝑠𝑡
10
Lecture Notes on DSF 2021
𝑍 = 𝐴1 𝑒 𝜆1 𝑡 + 𝐴2 𝑒 𝜆2 𝑡 ------------------- Eq.1.24
𝐶
The physical significance of this solution depends upon the relative magnitudes of (2𝑚)2 and
(K/m), which determines whether the exponents are real or complex quantities.
𝐶
Case I: ( )2 > (𝐾/𝑚
2𝑚
The roots λ1 and λ2 are real and negative. The motion of the system is not oscillatory but is an
exponential as shown in Fig.1.7).
11
Lecture Notes on DSF 2021
The system in this condition is known as critically damped system and C c is known as critical
damping constant.' The ratio of the actual damping constant to the critical damping constant is
defined as damping ratio:
𝐶
Damping ratio, 𝜉 = 𝐶
𝑐
𝐶
By substituting this value of' 𝜉 = 𝐶 in Eq. (1.23), we get
𝑐
𝐶 𝐾
Case III: ( )2 < ( 𝑚 )
2𝑚
The roots λ1 and λ2 are complex and are given as
𝜆1,2 = (−𝜉 ± 𝑖√1 − (𝜉)2 )𝜔𝑛 ------------------ Eq.1.28
The complete solution to the Eq.27, gives
2)𝜔 𝑡 2)𝜔 𝑡
𝑍 = 𝐴1 𝑒 (−𝜉+𝑖√1−𝜉 𝑛 + 𝐴2 𝑒 (−𝜉−𝑖√1−𝜉 𝑛 ----------------- - Eq.1.29
(𝑖√1−𝜉 2)𝜔𝑛 𝑡 (−𝑖√1−𝜉 2 )𝜔𝑛 𝑡
Or 𝑍 = 𝑒 −𝜉𝜔𝑛 𝑡 [𝐴1 𝑒 + 𝐴2 𝑒 ]--------------- Eq.1.30
The above equation can be written as
𝑍 = 𝑒 −𝜉𝜔𝑛 𝑡 [𝐶1 sin(𝜔𝑛 √1 − 𝜉 2 𝑡) + 𝐶2 cos(𝜔𝑛 √1 − 𝜉 2 𝑡)]-------------- Eq.1.31
Or 𝑍 = 𝑒 −𝜉𝜔𝑛 𝑡 [𝐶1 sin(𝜔𝑛𝑑 𝑡 + 𝐶2 cos(𝜔𝑛𝑑 𝑡]-------------- Eq.1.32
2
Where 𝜔𝑛𝑑 = 𝜔𝑛 (1 − 𝜉 )is known as damped natural frequency
The motion of the system is oscillatory (Fig.1.9) and the amplitude of vibration goes on
decreasing in an exponential fashion.
12
Lecture Notes on DSF 2021
Thus, damping of a system can be obtained from a free vibration record by knowing the
successive amplitudes which are one cycle apart.
If the damping is very small, it may be convenient to measure the differences in peak amplitudes
for a number of cycles, say n, as
1 𝑍
𝜉 = 2𝜋𝑛 𝑙𝑛 𝑍0 ------------------ Eq.1.36
𝑛
Therefore, a system is
Over damped if ξ> 1;
Critically damped if ξ = 1 and
Under damped if ξ< 1
term only.
The system will vibrate harmonically, with the same frequency as the forcing and the peak
amplitude is given by
𝐹0⁄
𝐾
𝐴𝑍 = -------------------------- Eq.1.53
𝜔2 𝜔
√(1− 2 )2+4𝜉 2(𝜔 )2
𝜔𝑛 𝑛
The quantity 𝐹0⁄𝐾 equals to the static deflection of the mass under force F0. Dynamic
magnification factor M is the ratio of the dynamic amplitude Az to the static deflection and is
given by
1
𝑀= ----------------------- Eq.1.54
𝜔2 𝜔
√(1− 2 )2 +4𝜉 2 (𝜔 )2
𝜔𝑛 𝑛
𝜔
It would be seen that the frequency ratio near (𝜔 =η=) 1, the value of frequency is maximum.
𝑛
This is called resonance and the forcing frequency at which this occurs is called as the resonant
frequency.
Differentiating Eq. (1.53) with respect to η and equating to zero, it can be shown that resonance
will occur at a frequency ratio given by
𝜂 = √1 − 2𝜉 2 ------------------- Eq.1.55
which is approximately equal to unity for small values of ξ
Now 𝜔𝑛𝑑 = 𝜔𝑛 √1 − 2𝜉 2 ------------------------- Eq.1.56
This is known as damped resonance frequency.
Maximum value of magnification factor can be obtained as
1
𝑀𝑚𝑎𝑥 = ------------------------- Eq.1.57
2𝜉√1−𝜉 2
14
Lecture Notes on DSF 2021
Example:1
An unknown weight W is attached to the end of an unknown spring k and natural frequency of
the system was found to be 90 cpm. If 1 kg weight is added to W, the natural frequency reduced
to 75 cpm. Determine the unknown weight W and spring constant k.
Sol:
𝜔𝑛 = 90𝑐𝑝𝑚
When 1 kg added to the weight W, the natural frequency reduced to 75 cpm
𝜔𝑛 = 90 cpm
Or f= 90/60=1.5 cps
𝜔 = 2𝜋𝑓 = 2𝜋 × 1.5𝑟/𝑠
𝐾𝑔⁄
𝜔2 = 𝐾⁄𝑚 = 𝑊 = 88.92--------------1
Again, f= 75/60=1.25 cps
And
𝜔 = 2𝜋𝑓 = 2𝜋 × 1.25 = 61.88
𝐾𝑔
⁄(𝑊 + 1) = 61.88----------------2
Solving for 1 and 2, we get W=2.27kg
And Spring constant K=201 kg/cm
Example 2:
A spring and dashpot are attached to a body weighing 140 N. The spring constant is 3.0 kN/m.
The dashpot has a resistance of 0.75 N at a velocity of 0.06 m/s. Determine the following for free
vibration:
(i) whether the system is over damped, under damped or critically damped
Sol:
Given:
K=3 kN/m, Damping force = 0.75 N at a velocity of 0.06 m/s
Hence damping coefficient C= 0.75/0.06=12.5 N.s/m
We know:
K =3kN/m
for over damped vibration
𝐶 𝐾
( )2 > ( )
2𝑚 𝑚
140 N For critical damping
𝐶 𝐾
( )2 = ( )
2𝑚 𝑚
For under damped
𝐶
(2𝑚)2 < (𝐾/𝑚)
𝐶 12.5×9.81
Now checking for damping condition, we have 2𝑚 = = 0.437
2×140
𝐾 3000×9.81
Again, √ = √ = 14.5
𝑚 140
15
Lecture Notes on DSF 2021
So the system is under damped.
Example 3:
An SDF system is excited by a sinusoidal force. At resonance the amplitude of displacement was
measured to be 2 mm. At an exciting frequency of one-tenth of the natural frequency of the
system, the displacement amplitude was measured to be 0.2 mm. Estimate the damping ratio of
the system.
Sol:
Given:
Umax= 2 mm
u=0.2mm at the exciting frequency of one-tenth of the natural frequency (At small frequency)
We know that
𝐹0⁄
𝑢= 𝐾
𝜔2 2 𝜔
√(1 − ) + 4𝜉 2 ( )2
𝜔𝑛2 𝜔𝑛
𝑢
At low frequency ratio 𝐹0 =1
⁄𝑘
𝑢𝑚𝑎𝑥 1 1
And 𝐹0 = ~
⁄𝑘 2𝜉√1−𝜉 2 2𝜉
0.2
Hence 𝐹0 =1
⁄𝐾
𝐹0
So =0.2
𝐾
2 1 2 1
Now 𝐹0 = , which gives =
⁄𝐾 2𝜉 0.2 2𝜉
0.2
Hence ξ= = 0.05 or 5%
4
Example 4:
A body weighing 600 N is suspended from a spring which deflects 12 mm under the load. It is
subjected to a damping effect adjusted to a value 0.2 times that required for critical damping.
Find the natural frequency of the un-damped and damped vibrations, and in the latter case,
determine the ratio of successive amplitudes.
Sol:
𝑊 600
𝐾= = = 5 × 104 𝑁⁄𝑚
𝛿 12 × 10−3
m=60 kg
𝐶
Damping ratio 𝜉 = 𝐶 = 0.2
𝑐
𝐾 5×104
Natural Frequency=√𝑚=√ =28.86 rpm
60
16
Lecture Notes on DSF 2021
𝑍
So 2𝜋 × 0.2 = 𝑙𝑛 𝑍1
2
𝑍1 2𝜋×0.2
Or 𝑍 = 𝑒 = 3.51
2
Problem No.1
For a machine foundation, given weight = 60 kN, spring constant = 11,000 kN/m, and c = 200
kN-s/m, determine
(a) whether the system is overdamped, underdamped, or critically damped,
(b) the logarithmic decrement, and
(c) the ratio of two successive amplitudes.
Problem No.2
For Problem No.1, determine the damped natural frequency.
Problem No. 3
A machine and its foundation weight 140 kN. The spring constant and the damping ratio of the
soil supporting the soil may be taken as 12 × 104 kN/m and 0.2, respectively. Forced vibration of
the foundation is caused by a force that can be expressed as Q (kN) = Q0 sin ωt
Q0 = 46 kN,ω = 157 rad/s
Determine
(a) the undamped natural frequency of the foundation,
(b) amplitude of motion, and
(c) maximum dynamic force transmitted to the sub-grade.
1.4 TWO DEGREES OF FREEDOM SYSTEMS
1.4.1 Undamped free vibration
Figure 1.10 shows a mass-spring system with two degrees of freedom.
The solutions of Eq. (1.58) and (1.59) will be of the following form
𝑍1 = 𝐴1 sin(𝜔𝑛 𝑡)--------------- -------------------- Eq.1.60
𝑍2 = 𝐴2 sin(𝜔𝑛 𝑡)----------------------------------- Eq.1.61
Substitution of Eqs. (1.20) and (1.61), into Eqs. (1.58) and (1.59) yields:
(𝐾1 + 𝐾2 − 𝑚1 𝜔𝑛2 )𝐴1 − 𝐾2 𝐴2 = 0-------------- Eq.1.62
(𝐾2 + 𝐾3 − 𝑚2 𝜔𝑛2 )𝐴2 + 𝐾2 𝐴1 = 0-------------- Eq.1.63
For nontrivial solutions of 𝜔𝑛 in Eqs. (1.62) and (1.63),
𝐾 + 𝐾2 − 𝑚1 𝜔𝑛2 −𝐾2
| 1 | = 0-------- Eq.1.64
−𝐾2 𝐾2 + 𝐾3 − 𝑚2 𝜔𝑛2
Or
𝐾1+𝐾2 𝐾2 +𝐾3 𝐾1 𝐾2 +𝐾2 𝐾3 +𝐾3 𝐾1
𝜔𝑛4 − ( + ) 𝜔𝑛2 + = 0----------- Eq.1.65
𝑚1 𝑚2 𝑚1 𝑚2
Equation (1.65) is quadratic in 𝜔𝑛2 , and the roots of this equation are:
1 𝐾1 +𝐾2 𝐾2 +𝐾3 𝐾1+𝐾2 𝐾2 +𝐾3 2 4𝐾22
𝜔𝑛2 − 2 [ + ] ± √( − ) +𝑚 --------------- Eq.1.66
𝑚1 𝑚2 𝑚1 𝑚2 1 𝑚2
From Eq.(9),two values of natural frequencies (𝜔𝑛1 )and (𝜔𝑛2 ) can be obtained.
𝜔𝑛1 , is corresponding to the first mode and 𝜔𝑛2 is of the second mode of vibration
The general equation of motion of the two masses can now be written as
𝑍1 = 𝐴11 𝑠𝑖𝑛𝜔𝑛1 𝑡 + 𝐴21 𝑠𝑖𝑛𝜔𝑛2 𝑡------------------- Eq.1.67
1 2
𝑍2 = 𝐴2 𝑠𝑖𝑛𝜔𝑛1 𝑡 + 𝐴2 𝑠𝑖𝑛𝜔𝑛2 𝑡------------------- Eq.1.68
The superscripts in A represent the mode.
The relative values of amplitudes A1 and A2 for the two modes can be obtained using Eqs.1.62
and 1.63. Thus
𝐴11 𝐾2 2
𝐾2 +𝐾3 −𝑚2 𝜔𝑛1
=𝐾 = -------------- Eq.1.69
𝐴12 2
1 +𝐾2 −𝑚1 𝜔𝑛1 𝐾2
𝐴21 𝐾2 2
𝐾2 +𝐾3 −𝑚2 𝜔𝑛2
=𝐾 = -------------- Eq.1.70
𝐴22 2
1 +𝐾2 −𝑚1 𝜔𝑛2 𝐾2
1.4.2 Undamped forced vibrations
Consider the system shown in Figure 1.11 with excitation force
F0 sin (ω t ) acting on mass m1. In this case, equations of motion will be:
𝑚1 𝑍1̈ + 𝐾2 𝑍1 + 𝐾2 (𝑍1 − 𝑍2 ) = 𝐹0 𝑠𝑖𝑛𝜔𝑡------------ Eq.1.71
AND
𝑚2 𝑍2̈ + 𝐾3 𝑍2 + 𝐾2 (𝑍2 − 𝑍1 ) =0---------------- Eq.1.72
For steady state, the solutions will be as
𝑍1 = 𝐴1 𝑠𝑖𝑛𝜔𝑡------------ Eq.1.73
AND
𝑍2 = 𝐴2 𝑠𝑖𝑛𝜔𝑡-------- Eq.1.74
18
Lecture Notes on DSF 2021
Substituting Eqs. (1.73) and (1.74) in Eqs. (1.71) and (1.72), we get
(𝐾1 + 𝐾2 − 𝑚1 𝜔2 )𝐴1 − 𝐾2 𝐴2 = 𝐹0 --------------- Eq.1.75
AND
−𝐾2 𝐴1 + (𝐾2 + 𝐾3 − 𝑚2 𝜔2 )𝐴2 = 0---------- Eq.1.76
The above Two equations give steady state amplitude of vibration of the two masses
respectively, as a function of ω. The denominator of the two equations is same. It may be noted
that:
(i) The expression inside the bracket of the denominator of Eqs.1.77 and 1.78 is of the same type
as the expression of natural frequency given by Eq. (1.66). Therefore at 𝜔 = 𝜔𝑛1 and 𝜔 =
𝜔𝑛2 values of A1 and A2 will be infinite as the denominator will become zero.
(ii) The numerator of the expression for Al becomes zero when
𝐾1+𝐾3
𝜔=√ -------------------- Eq.1.79
𝑚2
Thus it makes the mass m1 motionless at this frequency. No such stationary condition exists for
mass m1. The fact that the mass which is being excited can have zero amplitude of vibration
under certain conditions by coupling it to another spring-mass system forms the principle of
dynamic vibration absorbers which will be discussed latter on.
19
Lecture Notes on DSF 2021
1.5 SYSTEM WITH n DEGREES OF FREEDOM
1.5.1 Undamped free vibrations
Consider a system shown in Figure 1.12 having n-degree of freedom.
If Z1, Z2, Z3 ... Zn are the displacements of the respective masses at any instant, then equations of
motion are:
𝑚1 𝑍1̈ + 𝐾1 𝑍1 + 𝐾2 (𝑍1 − 𝑍2 ) = 0--------------------- Eq.1.80
𝑚2 𝑍2̈ − 𝐾2 (𝑍1 − 𝑍2 ) + 𝐾3 (𝑍2 − 𝑍3 ) = 0------------ Eq.1.81
-----------------------------------
-------------------------------------
𝑚𝑛 𝑍𝑛̈ − 𝐾𝑛 (𝑍𝑛−1 − 𝑍𝑛 ) = 0--------------------------- Eq.1.82
The solution of Eqs. (1.80) to (1.82) will be of as follows;
𝑍1 = 𝐴1 𝑠𝑖𝑛𝜔𝑛 𝑡------------------ Eq.1.83
𝑍2 = 𝐴2 𝑠𝑖𝑛𝜔𝑛 𝑡------------------- Eq.1.84
---------------------
𝑍𝑛 = 𝐴𝑛 𝑠𝑖𝑛𝜔𝑛 𝑡------------------ Eq.1.85
Substitution of Eqs. (1.83) to (1.85) into Eqs. (1.80) to (1.82), yields:
[(𝐾1 + 𝐾2 ) − 𝑚1 𝜔𝑛 2 ]𝐴1 − 𝐾2 𝐴2 = 0-------------------------------- Eq.1.86
2
−𝐾2 𝐴1 + [(𝐾2 + 𝐾3 ) − 𝑚2 𝜔𝑛 ]𝐴2 − 𝐾3 𝐴3 = 0-------- Eq.1.87
2
−𝐾3 𝐴2 + [(𝐾2 + 𝐾4 ) − 𝑚3 𝜔𝑛 ]𝐴3 − 𝐾4 𝐴4 = 0-------------- Eq.1.88
-------------------------------
-------------------------------
2
−𝐾𝑛 𝐴𝑛−1 + [𝐾𝑛 − 𝑚𝑛 𝜔𝑛 ]𝐴𝑛 = 0------------------- Eq.1.89
Equation (1.90) is of nth degree in 𝜔𝑛2 and therefore gives n values of 𝜔𝑛 corresponding to n
natural frequencies. The mode shapes can be obtained from Eq. (1.86 to 1.89) by using, at one
time, one of the various values of 𝜔𝑛 obtained from Eq. (1.90).
When the number of degrees of freedom exceeds three, the problem of forming the frequency
equation and solving it for determination of frequencies and mode shapes becomes tedious.
Numerical techniques are more useful in such cases. ,
Holzer's numerical technique is a convenient method of solving the problem for an idealized
system
20
Lecture Notes on DSF 2021
.
Fig. 1.13: An idealized multiple degree of freedom system
Inertia force at a level below mass 𝑚𝑖−1 = ∑𝑖−1 ̈
𝑗=1 𝑚𝑗 𝑍𝑗 -------------------- Eq.1.91
Spring force at that level corresponding to the difference of adjoining masses
21
Lecture Notes on DSF 2021
𝐾𝑖−1 (𝑍𝑖 − 𝑍𝑖−1 )------------------------- - Eq.1.92
Equating the above eqs, we obtain
∑𝑖−1 ̈
𝑗=1 𝑚𝑗 𝑍𝑗 = 𝐾𝑖−1 (𝑍𝑖 − 𝑍𝑖−1 )------------------- Eq.1.93
Putting 𝑍𝑖 = 𝐴𝑖 𝑠𝑖𝑛𝜔𝑛 𝑡 in Eq.1.93, we get
∑𝑖=1 2
𝑗=1 𝑚𝑗 (−𝐴𝑖 𝜔𝑛 𝑠𝑖𝑛𝜔𝑛 𝑡) = 𝐾𝑖−1 (𝐴𝑖 𝑠𝑖𝑛𝜔𝑛 𝑡 − 𝐴𝑖−1 𝑠𝑖𝑛𝜔𝑛 𝑡)--------------- Eq.1.94
2
𝜔𝑛
Or 𝐴𝑖 = 𝐴𝑖−1 − 𝐾 ∑𝑖−1 ̈
𝑗=1 𝑚𝑗 𝑍𝑗 ------------------- Eq.1.95
𝑖−1
Equation (1.95) gives a relationship between any two successive amplitudes. Starting with any
arbitrary value of Ai amplitude of all other masses can be determined. A plot of An+1 versus 𝜔𝑛2
would have the shape as shown in Figure 1.14. Finally An+1 should worked out to zero because
of base fixity.
The intersection of the curve with (𝜔𝑛2 ) axis would give various𝜔𝑛2 . The mode shape can be
obtained by substituting the value of 𝜔𝑛2 in Eq. (1.95).
𝐹0 𝜔2 𝜔2 𝑒
= 2𝑚𝑒 𝑒 = 2𝑚𝑒 𝑒 𝑚𝜔2 = (2𝑚𝑒 𝑚)𝜉 2 ---------- Eq.1.109
𝑘 𝐾 𝑛
2𝜂𝜉
𝜃 = 𝑡𝑎𝑛 −1 (1−𝜂2 )--------------- ---------------------- Eq.1.110
Differentiating Eq. (1.111) with respect to η and equating to zero. It can be shown that resonance
will occur at a frequency ratio given by:
1
𝜂= ------------------ ---------------- Eq.1.112
√1−2𝜉 2
𝜔𝑛
Or 𝜔𝑑 = -------------------------- ------ Eq.1.113
√1−2𝜉 2
By substituting Eq. (1.113) in Eq. (1.111), we get:
24
Lecture Notes on DSF 2021
𝐴𝑍 1
(2𝑚𝑒 𝑒⁄ = ----------------- Eq.1.114
𝑚)𝑚𝑎𝑥 2𝜉√1−𝜉 2
1
= For small damping
2𝜉
2𝑚𝑒 𝑒𝜔2⁄
𝐾
= sin(𝜔𝑡 − 𝜃)
√(1−𝜂 ) +4𝜉2 (𝜂)2
2 2
−1 2𝜂𝜉
Where, 𝜃 = 𝑡𝑎𝑛 [ ]-------------- Eq.1.117
1−𝜂2
The only force which can be applied to the foundation is the spring force KZ and the damping
force, 𝐶𝑍̇; hence the total force transmitted to the foundation during steady state forced vibration
is
𝐹𝑡 = 𝐾𝑍 + 𝐶𝑍̇---------------------------- --------------------------------------------Eq.1.118
Now substituting Eq. (1.116) in Eq. (1.118), we get
2𝑚𝑒 𝑒𝜔 2⁄
2𝑚𝑒 𝑒𝜔 2 𝐾
𝐹𝑡 = sin(𝜔𝑡 − 𝜃) + Cω cos(𝜔𝑡 − 𝜃)---------- Eq.1.119
√(1−𝜂 2 )2+4𝜉 2 (𝜂) 2 2
𝜔2 𝜔 2
√(1− 2 ) +4𝜉 2 ( )
𝜔𝑛 𝜔𝑛
25
Lecture Notes on DSF 2021
Where β is the phase difference between the exciting force and the force transmitted to the
foundation and is given by,
𝐶𝜔
𝛽 = 𝜃 − 𝑡𝑎𝑛−1 [ 𝐾 ]--------------------- ---------------------------------------------Eq.1.121
Since the force 2𝑚𝑒 𝑒𝜔2 is the force which would be transmitted if springs are infinitely rigid, a
measure of the effectiveness of the isolation mounting system is given by,
𝐹𝑡 √1+(2𝜂𝜉)2
𝜇 𝑇 = 2𝑚 = --------------- ----------------------------------------Eq.1.122
𝑒 𝑒𝜔2
√(1−𝜂2 )2+4𝜉2(𝜂)2
It will be noted from the figure that for any frequency ratio greater than√2, the force transmitted
to the foundation will be less than the exciting force. However in this case, the presence of
damping reduces the effectiveness of the isolation system as the curves for damped case are
above the undamped ones for η>√2. A certain amount of damping, however, is essential to
maintain stability under transient conditions and to prevent excessive amplitudes when the
vibrations pass through resonance during the starting or stopping of the machine. Therefore, for
the vibration isolation system to be effective η should be greater than√𝟐.
26
Lecture Notes on DSF 2021
27
Lecture Notes on DSF 2021
√1+(2𝜂𝜉)2
𝑍𝑚𝑎𝑥 = 𝑌0 ------------------------------------------------------------- Eq.1.127
√(1−𝜂2 )2 +(2𝜂𝜉)2
Equation (1.128) is the same expression as Eq. (1.122) obtained earlier. Transmissibility of such
system can also be studied from the response curves shown in Fig.1.18. It is again noted that for
the vibration isolation to be effective, it must be designed in such a way that η>√𝟐.
Materials used for vibration isolation are rubber, felt, cork and metallic springs. The
effectiveness of each depends on the operating conditions.
i) Rubber: Rubber is loaded in compression or in shear; the latter mode gives higher
flexibility. With loading greater than about 0.6 N per sq mm, it undergoes much faster
deterioration. Its damping and stiffness properties vary widely with applied load,
temperature, shape factor, excitation frequency and the amplitude of vibration. The
maximum temperature up to which rubber can be used satisfactorily is about 65°c. It
must not be used in presence of oil which attacks rubber. It is found very suitable for
high frequency vibrations.
ii) Felt: Felt is used in compression only and is capable of taking extremely high loads.
It has very high damping and so is suitable in the range of low frequency ratio. It is
mainly used in conjunction with metallic springs to reduce noise transmission.
28
Lecture Notes on DSF 2021
iii) Cork: Cork is very useful for acoustic isolation and is also used in small pads placed
underneath a large concrete block. For satisfactory working it must be loaded from 10
to 25 N/sq mm. It is not affected by oil products or moderate temperature changes.
However, its properties change with the frequency of excitation.
iv) Metallic springs: Metallic springs are not affected by the operating conditions or the
environments. They are quite consistent in their behaviour and can be accurately
designed for any desired conditions. They have high sound transmissibility which can
be reduced by loading felt in conjunction with it. It has negligible damping and so is
suitable for working in the range of high frequency ratio.
The purpose of a vibration measuring instrument is to give an output signal which represents, as
closely as possible, the vibration phenomenon. This phenomenon may be displacement, velocity
or acceleration of the vibrating system and accordingly the instrument which reproduces signals
proportional to these are called vibrometers, velometers or accelerometers.
There are essentially two basic systems of vibration measurement. One method is known as the
directly connected system in which motions can be measured from a reference surface which is
fixed. More often such a reference surface is not available. The second system, known as
“Seismic System" does not require a fixed reference surface and therefore is commonly used for
vibration measurement.
𝜔
Where 𝜂 = 𝜔 = Frequency ratio
𝑛
𝜉 =damping ratio
2𝜂𝜉
𝜃 = 𝑡𝑎𝑛 −1 ( 2 )------------------- --------------------------- Eq.1.132
1−𝜂
Equation (1.131)can be rewritten as
𝑿 = 𝜼𝟐 µ𝒀𝟎 𝐬𝐢𝐧(𝝎𝒕 − 𝜽)------------ -------------------------- Eq.1.133
Where
1
µ= --------------------------------------------- Eq.1.134
√(1−𝜂2 )2 +(2𝜂𝜉)2
30
Lecture Notes on DSF 2021
It is seen when η is large, 𝜂2 µ is approximately equal to 1 and θ is approximately equal to 180°.
Therefore to design a displacement pickup, η should be large which means that the natural
frequency of the instrument itself 'should be low compared to the frequency to be measured. Or
in other words, the instrument should have a soft spring and heavy mass. The instrument is
sensitive, flimsy and can be used in a weak vibration environment. The instrument cannot be
used for measurement of strong vibrations.
1.8.2 Acceleration Pickup (Accelerometer)
Equation (1.133) can be rewritten as
1
𝑋̈ = 𝜔2 µ𝜔2 𝑌0 sin(𝜔𝑡 − 𝜃)----------------------------------------- Eq.1.135
𝑛
The output of the instrument will be proportional to the acceleration of the structure if µ is
constant. It is seen that µ is approximately equal to unity for small values of η. Therefore to
design an acceleration pick up, it should be small which means that the natural frequency of the
instrument itself should be high compared to the frequency to be measured. In other words, the
instrument should have a stiff spring and small mass. The instrument is less sensitive and
suitable for the measurement of strong motion. The instrument size is small.
1.8.3 Velocity Pickup
Equation (1.133) can be rewritten as
1
𝑋̇ = 𝜔 𝜂µ𝑌0 𝜔sin(𝜔𝑡 − 𝜃)----------------------------------------- Eq.1.136
𝑛
1
The output of the instrument will be proportional to velocity of the structure if 𝜂µ is a
𝜔𝑛
constant.
At η= 1, Eq. (1.136) can be written as
1 1 1
𝑋́ = 𝜔 𝑌 𝜔sin(𝜔𝑡 − 𝜃)-----------------------------------------Eq.1.137 as at η=1,µ = 2𝜉
𝑛 2𝜉 0
Since 𝜔𝑛 and ξ are constant, the instrument will measure the velocity at η= 1.
It may be noted that the same instrument can be used to measure displacement, acceleration and
velocity in different frequency ranges.
𝑋𝛼𝑌, 𝑖𝑓𝜂 ≫ 1, Displacement pickup (Vibrometer)
𝑋𝛼𝑌, 𝑖𝑓𝜂 ≪ 1, Acceleration pickup (Accelerometers)
𝑋𝛼𝑌, 𝑖𝑓𝜂 = 1, Velocity pickup (Velometers)
Displacement and velocity pickups have the disadvantage of having rather a large size if motions
having small frequency of vibration are to be measured. Calibration of these pickups is not
31
Lecture Notes on DSF 2021
simple. Further corrections have to be made in the observations as the response is not flat in the
starting regions. From the point of view of small size, flat frequency response, sturdiness and
ease of calibration, acceleration pickups are to be favored. They are relatively less sensitive and
this disadvantage can easily be overcome by high gain electronic instrumentation.
1.8.4 Transducer
A transducer is a device for converting the mechanical motion of vibration into an electrical
signal, commonly called pickup.
There are three kinds of transducers: Displacement, Velocity, and Acceleration
It is the most common type of transducer which is operated on the eddy current principle. It sets
up a high-frequency electric field in the gap between the end of the Proximity Probe and the
metal surface that is moving. It senses the change in the gap and measures relative displacement
not absolute displacement.
Proximity Probe
It is sensitive to shaft surface defects such as scratches, dents and vibrations in conductivity and
permeability.
Senses shaft run out, and it is very difficult to distinguish vibration from run out.
The practical maximum frequency of proximity probes is about 1500Hz.The minimum frequency
is zero. It can also measure static displacement. A useful application of proximity probes is to
measure very slow relative movement like thermal expansion. It is useful in situations where the
vibrating part cannot tolerate the mass of the pickup.
32
Lecture Notes on DSF 2021
1.8.6 Velocity Transducer
Typically accelerometer has very high natural frequency, typically 25000 Hz Its response is
linear for about 1/3 of this range. It has a useful frequency range of from about 5 to
approximately 100000 Hz depending on its size. The primary considerations in selecting an
accelerometer are sensitivity and frequency response.
If high-amplitude motions are to be measured, i.e. greater than 10g, such as in shock
measurement, then a low-sensitivity accelerometer is appropriate 10 mV/g or less.
If the level motion is to be measured, such as building or structural motions at low frequencies
then a high sensitivity accelerometer should be chosen 1000 mv/g.
For most machinery monitoring,100 mV/g sensitivity accelerometer provide the right balance of
sensitivity and frequency response. Other considerations in accelerometer selection or transducer
are Temperature exposure
Linearity - It is expressed as the percent deviation from a constant value of the sensitivity.
Transverse Sensitivity is the ability of the transducer to detect motion in directions perpendicular
to its sensitive axis.
Damping is very low in piezoelectric accelerometer but can be significant in other types, such as
piezo-resistive accelerometer. Strain sensitivity is the ability of the transducer to generate a
signal when the base is distorted, such as when it is clamped against a non flat surface.
34
Lecture Notes on DSF 2021
2.0 WAVE PROPAGATION; BASIC ELASTIC PROPERTIES AND RELATIONSHIP
2.1 Elastic Constants
An elastic material is one which obeys Hook's law of proportionally between stress and strain.
For an isotropic elastic material subjected to normal stress 𝜎𝑥 in the x-direction, the strains in x, y,
z directions are given as
𝜎𝑥
𝜀𝑥 = --------------- Eq.2.1
𝐸
𝜎𝑥
𝜀𝑦 = 𝜀𝑧 = −𝜇 ---------- Eq.2.2
𝐸
In this section, the propagation of stress waves in an infinite, homogeneous, isotropic medium
presented in Figure.2.1 shows the stresses acting on a soil element with sides dx, dy, dz. For
obtaining the differential equations of motion, the sum of the forces acting parallel to each axis is
considered.
In the x-direction the equilibrium equation is given as
36
Lecture Notes on DSF 2021
In the above expressions, ρ is the mass density of the soil; u, v and ω are displacements in the x,
y, and z directions respectively. To express the right hand sides of these Eqs., the relationship for
an elastic medium given is used. The equations for strains and rotations of elastic and isotropic
materials in terms of displacements are as follows:
37
Lecture Notes on DSF 2021
Equations (2. 27) are the equations of motion of an infinite homogeneous, isotropic, and
elastic medium. On differentiating these equations with respect to x, y and z, respectively, and
adding, we get
𝜕2 𝜕𝑢 𝜕𝑣 𝜕𝜔 𝜕 2𝜀̅ 𝜕 2 𝜀̅ 𝜕 2 𝜀̅ 𝜕𝑢 𝜕𝑣 𝜕𝜔
𝜌 𝜕𝑡 2 [𝜕𝑥 + 𝜕𝑦 + ] = (𝜆 + 𝐺 ) [ + 𝜕𝑦 2 + 𝜕𝑧 2 ] + 𝐺∇2 (𝜕𝑥 + 𝜕𝑦 + )------------- Eq.2. 28
𝜕𝑧 𝜕𝑥 2 𝜕𝑧
𝜕 2 𝜀̅
𝜌 𝜕𝑡 2 = (𝜆 + 𝐺 )(∇2 𝜀̅) + (𝐺∇2 𝜀̅)----------------------- Eq.2. 29
𝜕 2 𝜀̅
Hence, 𝜌 𝜕𝑡 2 = (𝜆 + 2𝐺 )(∇2 𝜀̅)--------------- Eq.2. 30
0r
𝜕 2 𝜀̅ (𝜆+2𝐺)
= (∇2 𝜀̅) = 𝑉𝑝2 ∇2 𝜀̅------------------ Eq.2. 31
𝜕𝑡 2 𝜌
(𝜆+2𝐺)
Where 𝑉𝑝2 = ----------------------------- Eq.2.32
𝜌
Vp is the ve1ocity of compression waves which is also referred as primary wave or, P-wave. It is
important to note the difference in the wave velocities for an infinite elastic medium with those
(𝜆+2𝐺)
obtained for an elastic rod is, 𝑉𝑐 = √𝐸⁄𝜌: but in the infinite medium, 𝑉𝑝 = √ . This means
𝜌
that Vp >Vc, that is compression wave travels faster in infinite medium. It is due to the fact that
in infinite medium, there are no lateral displacements, while in the elastic rod lateral
displacements are possible.
2.2.3 Shear-Waves
Differentiating Eq. (2.27,b) with respect to z and Eq. (2.27,c) with respect to y, we get
𝜕2 𝜕𝑣 𝜕𝜀̅ 𝜕𝑣
𝜌 𝜕𝑡 2 ( 𝜕𝑧) = (𝜆 + 𝐺 ) (𝜕𝑦)(𝜕𝑧) + 𝐺∇2 𝜕𝑧----------------- Eq. 2.33
𝜕2 𝜕𝜔 𝜕𝜀̅ 𝜕𝜔
𝜌 𝜕𝑡 2 ( 𝜕𝑦 ) = (𝜆 + 𝐺 ) (𝜕𝑦)(𝜕𝑧) + 𝐺∇2 𝜕𝑦 --------------- Eq.2.34
Subtracting Eq.2.34 from Eq.2.33, we get
𝜕2 𝜕𝜔 𝜕𝑣 𝜕𝜔 𝜕𝑣
𝜌 𝜕𝑡 2 ( 𝜕𝑦 − 𝜕𝑧) = 𝐺∇2 ( 𝜕𝑦 − 𝜕𝑧 )----------------------- Eq.2.35
FromEq.(2.22,a)
𝜕𝜔 𝜕𝑣
2𝜔
̅̅̅̅𝑥 = −
𝜕𝑦 𝜕𝑧
Therefore,
𝜕 2 ̅̅̅̅
𝜔𝑥
𝜌 = 𝐺∇2 ̅̅̅̅---------------------
𝜔𝑥 Eq.2.35
𝜕𝑡 2
Or,
𝜕 2 ̅̅̅̅
𝜔𝑥 𝐺
= 𝜌 ∇2 ̅̅̅̅
𝜔𝑥 = 𝑉𝑠2 ∇2 ̅̅̅̅---------------------
𝜔𝑥 Eq.2.36 (a)
𝜕𝑡 2
Similar expression can be obtained for 𝜔
̅𝑦 𝑎𝑛𝑑𝜔
̅̅̅̅𝑧 as
𝜕 2 ̅̅̅̅
𝜔𝑦 𝐺
= 𝜌 ∇2 ̅̅̅̅
𝜔𝑦 = 𝑉𝑠2 ∇2 ̅̅̅̅----------------------
𝜔𝑦 Eq.2. 36 (b)
𝜕𝑡 2
𝜕 2 ̅̅̅̅
𝜔𝑧 𝐺
= ∇2 ̅𝜔
̅̅̅𝑧 = 𝑉𝑠2 ∇2 ̅𝜔
̅̅̅---------------------
𝑧 Eq.2. 36 (c)
𝜕𝑡 2 𝜌
The above expressions indicate that the Rotation is propagated with velocity Vs which is equal to
38
Lecture Notes on DSF 2021
√𝐺⁄𝜌. Shear wave is also referred as distortion wave or S-wave. It may be noted that shear wave
propagates at the same velocity in both the rigid elastic medium like rod or bar and the infinite-
medium.
2.3 WAVEPROP AGATION IN ELASTIC HALF-SPACE
In an elastically homogeneous ground, stressed suddenly at a point 'S' near its surface as shown
in (Figure 2.2), three elastic waves travel outwards at different speeds. Two are body waves;
which are propagated as spherical, fronts affected only a minor extent by the free surface of the
ground, and the third is a surface wave which is confined to the region, near the free surface.
39
Lecture Notes on DSF 2021
40
Lecture Notes on DSF 2021
𝜕2
(𝜕𝑡 2 ) {𝐹(𝑧) exp[𝑖(𝜔𝑡 − 𝑓𝑥 )]}=𝑉𝑝2 ∇2 {𝐹(𝑧) exp[𝑖(𝜔𝑡 − 𝑓𝑥 )]}---------- Eq.2.52
Or −𝜔2 𝐹(𝑧) = 𝑉𝑝2 [𝐹 " (𝑧) − 𝑓 2 𝐹 (𝑧)]--------------------------------------- Eq.2.53
Where F(z) and G(z) are functions of depth
Similarly, substituting Eq. (2.51) into Eq. (2.49) results in
−𝜔2 𝐺 (𝑧) = 𝑉𝑠2 [𝐺 " (𝑧) − 𝑓 2 𝐺 (𝑧)]---------------- Eq.2.54
Where
𝜕 2 𝐹(𝑧)
𝐹 " (𝑧 ) = ---------------------------------------- Eq.2.55
𝜕𝑧 2
and
𝜕 2 𝐺(𝑧)
𝐺 " (𝑧 ) = ------------------------------------------------- Eq.2.56
𝜕𝑧 2
Now Equations (2.53) and (2.54) can be rearranged to the form
𝐹 " (𝑧) − 𝑞 2 𝐹(𝑧) = 0-------------------------------- Eq.2.57
𝐺 " (𝑧) − 𝑠 2 𝐹(𝑧) = 0--------------------------------- Eq.2.58
Where
𝜔2
𝑞 2 = 𝑓 2 − 𝑉 2 ----------------------------------------------- Eq.2.59
𝑝
𝜔2
𝑠 2 = 𝑓 2 − 𝑉 2 ------------------- Eq.2.60
𝑠
2 2
4𝐺𝑓 2 𝑠𝑞 = (𝑠2 + 𝑓 ) [(𝜆 + 2𝐺)𝑞2 − 𝜆𝑓 ]------- Eq.2.74
The term 𝛼 2 can also be expressed in terms of Poisson’s ratio. From the relations given in Eq.
(2.81),
2𝜇𝐺
---------------------- Eq.2.81
1−2𝜇
Substitution of this relation in Eq. (2.80) yields,
𝐺 (1−2𝜇)
𝛼 2 = 𝜆+2𝐺 = 2(1−𝜇)--------------- Eq.2.82
Again, substituting Eqs. (2.78), (2.79), and (2.80) into Eq. (2.76),
16(1 − 𝛼 2 𝑉 2 )(1 − 𝑉 2 ) = (2 − 𝑉 2 )(2 − 𝑉 2 )2
Or, 𝑉 6 − 8𝑉 4 − (16𝛼 2 − 24)𝑉 2 − 16(1 − 𝛼 2 ) = 0------------- Eq.2.83
42
Lecture Notes on DSF 2021
Equation (2.83) is a cubic equation in V2. For a given value of Poisson’s ratio, the proper value
of V2 can be found and, hence, so can the value of Vr in terms of Vp or Vs.
Example 1:
Given μ = 0.25, determined the value of the Rayleigh wave velocity in terms of Vs
Solution:
𝑉 6 − 8𝑉 4 − (16𝛼 2 − 24)𝑉 2 − 16(1 − 𝛼 2 ) = 0
For µ=0.25
1 − 2𝜇
𝛼2 = = 1/3
2 − 2𝜇
1 1
𝑉 6 − 8𝑉 4 − (16 × − 24) 𝑉 2 − 16 (1 − ) = 0
3 3
3𝑉 6 − 24𝑉 4 + 56𝑉 2 − 32 = 0
(𝑉 2 − 4)(3𝑉 4 − 12𝑉 2 + 8) = 0
2 2
Therefore, 𝑉 2 = 4, 2 + ,2−
√3 √3
𝑆2
If V2=4, 𝑓2 = 1 − 𝑉 2 = 1 − 4 = −3
2
So S/f is imaginary. This is also the case for V2=2 +
√3
It can be seen that when q/f and s/f are imaginary, it does not yield the primary and secondary
waves as discussed.
2 𝑣
For V2=2 − , 𝑣 = 𝑣𝑟 =0.9194
√3 𝑠
Or 𝑣𝑟 =0.9194𝑣𝑠
43
Lecture Notes on DSF 2021
3.0 LIQUEFACTION OF SOIL
Previous earthquake devastation was an illustration of catastrophic damages to structures and
resulting in loss of life which was due to liquefaction phenomenon. Liquefaction is defined as a
condition where a soil will undergo continuation of deformation at a constant low residual stress
or with no residual resistance, due to the build-up and maintenance of high pore water pressure
which reduces the effective confining pressure to a very low value. The pore pressure so build-up
leading to true liquefaction of this type may be due either to static or cyclic stress applications.
It denotes a condition where, during the course of cyclic stress applications, the residual pore water
pressure on completion of any full stress -cycle becomes equal to the applied confining stress.
MECHANISM OF LIQUEFACTION:
The strength of sand is primarily due to internal friction. In saturated state it may be expressed as
𝑠 = ̅̅̅tanφ------
𝜎𝑛 Eq.3.1
Where S= Shear strength of sand
𝜎𝑛 Effective normal stress on any plane at a depth of z
̅̅̅=
Φ= Angle of internal friction
45
Lecture Notes on DSF 2021
imparted to the deposit. Further, horizontal vibrations have more severe effect than
vertical vibrations. Multi directional shaking is more severe than one directional loading
(Seed at al.1977), as the pore water pressure build up is much faster and the stress ratio
required is about 10 percent less than that required for unidirectional shaking.
e) Location of Drainage and Dimension of Deposit: Sands are more pervious than fine
grained soil. However, if a pervious deposit has large dimensions, the drainage path
increases and the deposit may behave as un-drained, thereby, increasing the chances of
liquefaction of such a deposit. The drainage path is reduced by the introduction of drains
made out of highly pervious material.
f) Surcharge Load: If the surcharge load, i.e. the initial effective stress is large, then
transfer of stress from soil grains to pore water will require higher intensity vibrations or
vibration for a longer duration. If the initial stress condition is not isotropic as in field,
then stress condition causing liquefaction depends upon K0 (coefficient of earth pressure
at rest) and for K0> 5, the stress condition required to cause liquefaction increases by at
least 50%.
g) Method of Soil Formation: Sands unlike clays do not exhibit a characteristics structure.
But recent investigations show that liquefaction characteristics of saturated sands under
cyclic loading are significantly influenced by method of sample preparation and by soil
structure.
h) Period under Sustained Load: Age of sand deposit may influence liquefaction
characteristics. A 75% increase in liquefaction resistance has been reported on
liquefaction of undisturbed sand compared to its freshly prepared sample which may be
due to some form of cementation or welding at contact points of sand particles and
associated with secondary compression of soil.
i) Previous Strain History: Studies on liquefaction characteristics of freshly deposited
sand and of similar deposit previously subjected to some strain history reveal, that
although the prior strain history caused no significant change in the density of the sand, it
increased the stress that causes liquefaction by a factor of 1.5.
j) Trapped Air: If air is trapped in saturated soil and pore pressure develop, a part of it is
dissipated due to the compression of air, hence trapped air helps to reduce the possibility
of liquefaction.
k) Groundwater Table: The most conducive condition to liquefaction is near the surface of
ground water table. Unsaturated soil located above the groundwater table will not liquefy.
At the location where groundwater table significantly fluctuates, the liquefaction will also
fluctuate.
At a depth below the ground surface, liquefaction will occur if shear stress induced by
earthquake is more than the shear stress predicted. By comparing the induced and predicted shear
stresses at various depths, liquefaction zone can be obtained.
46
Lecture Notes on DSF 2021
In a sand deposit consider a column of soil of height h and unit area of cross section subjected to
maximum ground acceleration Qmax(Fig.3.2).
Assuming the soil column to behave as a rigid body, the inertia force F can be obtained as
𝑊
𝐹= 𝑎𝑚𝑎𝑥 --------- Eq.3.6
𝑔
𝛾𝑧 𝑎𝑚𝑎𝑥
Or = 𝑔 𝑎𝑚𝑎𝑥 = 𝜎0 -------------- Eq.3.7
𝑔
The maximum shear stress 𝜏𝑚𝑎𝑥 at a depth h is given by
𝐹 𝑎𝑚𝑎𝑥
𝜏𝑚𝑎𝑥 = 𝐴=𝜎0 ----------------- Eq.3.8
𝑔
As base area of soil column is taken as unity
Where g = Acceleration due to gravity and γ = Unit weight of soil
Since the soil column behaves as a deformable body, the actual shear stress at depth h, (𝜏𝑚𝑎𝑥 ) is
taken as
𝛾ℎ
𝜏𝑎𝑐𝑡 = 𝑟𝑑 𝜏𝑚𝑎𝑥 = 𝑟𝑑 ( 𝑔 )𝑎𝑚𝑎𝑥 ------ Eq.3.9
Where 𝑟𝑑 = Depth reduction factor
If linear variation is assumed between reduction factor and depth, than rd can be taken as
𝑟𝑑 = 1 − 0.012𝑧------- Eq.3.10
The above relation is valid for depth up to 15 m.
According to Seed and Idriss (1971), the average equivalent uniform shear stress 𝜏𝑎𝑣𝑔 is about
65 percent of the maximum shear stress 𝜏𝑚𝑎𝑥 .Therefore
𝛾𝑧
𝜏𝑎𝑣𝑔 = 0.65 𝑔 𝑎𝑚𝑎𝑥 𝑟𝑑 -------- Eq.3.11
The corresponding number of significant cycles N s for 𝜏𝑎𝑣𝑔 is given in table Table.3.1
Table 3.1: Significant cycles N s corresponding to 𝜏𝑎𝑣𝑔
Earthquake magnitude, M on Richter’s scale Ns
7 10
7.5 20
8.0 30
47
Lecture Notes on DSF 2021
In order to felicitate liquefaction analysis one non dimensional parameter known as Cyclic Stress
Ratio (CSR) or Seismic Stress Ratio (SSR), can be defined as
𝜏𝑎𝑣𝑔
𝐶𝑆𝑅𝑜𝑟𝑆𝑆𝑅 = ---------- Eq.3.12
𝜎0′
𝜎 𝑎𝑚𝑎𝑥
Thus, 𝐶𝑆𝑅𝑜𝑟𝑆𝑆𝑅 = 0.65𝑟𝑑 (𝜎0′ ) ------ Eq.3.13
0 𝑔
Seed and Idriss (1971) suggested the value of cyclic stress ratio values Cr as given in Table 3.2
It was observed that up to a relative density of 80%, the peak pulsating shear stress causing
liquefaction increases almost linearly with the increase in relative density. Keeping this fact in
view, the following general relation is suggested:
𝜏 𝜎 𝐷
(̅𝜎̅̅ℎ̅)𝑓𝑖𝑒𝑙𝑑𝐷𝑅 = (2𝜎𝑑 )𝑡𝑟𝑖𝑎𝑥,50 𝐶𝑟 50𝑅 ---------- Eq.3.14
𝑣 3
𝜏
Where, (̅𝜎̅̅ℎ̅)𝑓𝑖𝑒𝑙𝑑𝐷𝑅 = Cyclic shear stress ratio in field at relative density of D R percentage
𝑣
𝜎𝑑
(2𝜎 )𝑡𝑟𝑖𝑎𝑥,50 =Stress ratio obtained from triaxial test at relative density of 50%. It can be
3
determined from Fig.3.3
Fig. 3.3: Stress condition causing liquefaction of sands (Seed and Idriss, 1971)
48
Lecture Notes on DSF 2021
ii) Using Eq. 11 determine 𝜏𝑎𝑣𝑔 at depth h below ground surface.
𝜎𝑑
iii) Using Fig.2, determine the value of for given value of D50 of soil and number of
2𝜎3
equivalent cycles NS for the relative density of 50%.
𝜏
iv) Using Eq. 12, determine the value of (̅𝜎̅̅ℎ̅)𝑓𝑖𝑒𝑙𝑑𝐷𝑅 for the relative density of DR of the
𝑣
𝜏ℎ
soil at site. Multiplying (̅𝜎̅̅̅)𝑓𝑖𝑒𝑙𝑑𝐷𝑅 with effectives stress at depth h, we can obtain
𝑣
the value of shear stress 𝜏ℎ required for causing liquefaction.
v) At depth h, liquefaction will occur if 𝜏𝑎𝑣𝑔 > 𝜏ℎ
vi) Repeat steps (ii) to (iv) for other values of h to locate the zone of liquefaction.
𝜏𝑎𝑣𝑔 and 𝜏ℎ can be plotted in a graph to identify the zone of liquefaction.
Problem No.1
At a given site, a boring supplemented with standard penetration tests was done up to 15.0m
depth. The results of the boring are as given below:
Depth Classification of soil D50 N-Value DR Remarks
(m) (mm) (%)
1.5 SP 0.18 3 19 Position of ground water
3.0 SP 0.2 5 30 lies 1.5 m below the ground
4.5 SM 0.12 6 35 surface
6.0 SM 0.14 9 40
7.5 SM 0.13 12 45 𝛾𝑚𝑜𝑖𝑠𝑡 =19 kN/m3
9.0 SP 0.16 17 52 𝛾𝑠𝑢𝑏 =10 kN/m3
10.5 SW 0.2 20 52
12.0 SW 0.22 18 46
13.0 SW 0.22 24 60
15.0 SW 0.24 30 65
The site is located in seismically, active region, and is likely to be subjected by an earthquake of
Magnitude 7.5. Determine the zone of liquefaction using Seed and ldriss (1971) method.
Fig. 3.4: Chart for correction of N value for over burden pressure
3.7 FOLLOWING PROCEDURE FOR LIQUEFACTION ANALYSIS IS USED
i) Establish the design earthquake, and obtain the peak ground acceleration a max. Also
obtain number of significant cycles corresponding to the magnitude of earthquake
using Table 3.1.
ii) Using Eq. 3.11 determine 𝜏𝑎𝑣𝑔 at depth h below ground surface.
iii) Determine the value of standard penetration resistance value (N) at depth h below
ground surface. Obtain corrected Ncorrected value after applying overburden correction
to N using Fig.3.4
𝜏
iv) Using Fig.3.3, determine ̅𝜎̅̅ℎ̅for the given magnitude of earthquake and N corrected value
𝑣
𝜏
obtained in step (iii). Multiplying 𝜎̅̅̅ℎ̅ with effective stress at depth h below ground
𝑣
surface, obtain the value of shear stress 𝜏ℎ required for causing liquefaction.
v) At depth h, liquefaction will occur if
𝜏𝑎𝑣𝑔 > 𝜏ℎ
vi) Repeat steps (ii) to (v) for other values of h to locate the zone of liquefaction.
Example No.1
At a given site boring supplement with SPT was done up to 20 m depth. The results of the boring
are given below. Water table lies 2 m below the ground surface. Take 𝛾𝑠𝑢𝑏 =10 kN/m3. The site is
located in seismically active zone and the likely to be subjected by an earthquake of magnitude
7.5 and maximum ground acceleration is 0.15g. Find the zone of liquefaction if any.
Depth(m) 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 20.0
Nfield 4 4 5 7 9 10 12 14 16 18
Solution:
The effective stress, 𝜎 ′ = 𝛾𝑠𝑢𝑏 𝑍, Reduction factor 𝑟𝑑 = 1 − 0.012𝑧
The calculation are tabulated as below
50
Lecture Notes on DSF 2021
Depth N 𝜎′ 𝐶𝑁 𝑁𝑐𝑜𝑟𝑟𝑒𝑐𝑡 𝜏ℎ 𝜎0 𝑟𝑑 𝑎𝑚𝑎𝑥 𝜏ℎ 𝜏𝑎𝑣𝑔
0.65
(m) =NCN 𝜎′ 𝑔
100 × 𝑟𝑑
=√ ′
𝜎
2 4 20 2.24 8.96 0.09 20 0.976 0.0952 1.8 1.9
4 4 40 1.58 6.32 0.07 60 0.952 0.0928 2.8 5.56
6 5 60 1.29 6.45 0.075 100 0.928 0.0905 4.5 9.05
8 7 80 1.12 7.84 0.08 140 0.904 0.0881 6.4 12.33
10 9 100 1.0 9.0 0.1 180 0.88 0.0858 10 15.44
12 10 120 0.91 9.1 0.1 220 0.856 0.0835 12 18.37
14 12 140 0.84 10.8 0.12 260 0.832 0.0811 16.8 21.08
16 14 160 0.79 11.06 0.13 300 0.808 0.0788 20.8 23.64
18 16 180 0.74 11.84 0.14 340 0.784 0.0764 25.2 25.97
20 18 200 0.707 12.73 0.16 380 0.76 0.0741 32 28.15
From the above calculation, it is found that up to 18 m depth, 𝜏𝑎𝑣𝑔 > 𝜏ℎ . Hence liquefaction can
occur up to depth of 18 m from the ground surface.
51
Lecture Notes on DSF 2021
Fig.3.5: Chart to determine the cyclic resistance ratio for clean and silty sand for M=7.5
earthquake (After Seed et al.1975)
Example No.2
The sand deposit of fine sand (finer ≤ 5%) of finite thickness is located at a depth of 3.0 m from
the ground surface and ground water table is located at 1.5m below the ground surface. This is
located in seismic prone area where the anticipated GPA is 0.40g.The standard penetration test
was performed at depth of 3.0m. The corrected N value is 8. The unit weight of sand is taken
as18.4 kN/m3. Calculate the factor of safety against liquefaction for the saturated sand.
Sol:
Ncorrected = 8, Unit weight of sand = 18.4 kN/m3
Submerged unit weight = 8.59 kN/m3
PGA=0.4g
Effective stress 𝜎0′ = 𝜎0 − 𝑢=18.4x1.5+1.5x8.59=40.485 kN/m3
Total stress 𝜎0 = 𝜎0′ + 𝑢=40.485+1.5x9.81=55.2 kN/m3
Now using the linear relationship for the stress, depth reduction factor can be computed as
𝑟𝑑 = 1 − 0.012𝑧=1-0.012x3=0.964
52
Lecture Notes on DSF 2021
𝜎 𝑎𝑚𝑎𝑥 55.2
So CSR=0.65x𝑟𝑑 (𝜎0′ ) ( )=0.65 × 0, .964 × 40.485 × 0.4=0.342
0 𝑔
Using Fig.4 with Ncorrected=8, CRR can be obtained as 0.09
0.09
Hence 𝐹𝑆 = 0.342=0.263
So based on the calculation of factor of safety against liquefaction, the sand deposit liquefy.
Problem No 1:
A 10 m thick loose sand deposit (Dr=42%, finer≤5%) is saturated below a depth of 4 m. The
sand layer region is highly prone to liquefaction. Estimate the ground acceleration that would be
required to produce sand soils in a M=7.5 earthquake.
53
Lecture Notes on DSF 2021
4.0 DYNAMIC SOIL PROPERTIES
4.1 Laboratory Method
The soil properties which are needed in analysis and design of a structure subjected to dynamic
loading are:
(a) Dynamic moduli, such as Young's modulus E, shear modulus G, and bulk modulus K
(b) Poisson's ratio µ
(c) Dynamic elastic constants, such as coefficient of elastic uniform compression Cu, coefficient
of elastic uniform shear, 𝐶𝜏 , coefficient of elastic non-uniform Compression 𝐶∅ and coefficient of
elastic non-uniform shear 𝐶𝜑
(d) Damping ratio, ξ
(e) Liquefaction parameters, such as cyclic stress ratio, cyclic deformation and pore pressure
response.
f) Strength-deformation characteristics in terms of strain rate effects.
Since the dynamic properties of soils are strain dependent various laboratory and field techniques
have been developed to measure these properties over a wide range of strain amplitudes.
4.2 LABORATORY TECHNIQUES
The laboratory methods used for determining the dynamic properties of soils are:
i) Resonant column test,
ii) Ultrasonic pulse test,
iii) Cyclic simple shear test,
iv) Cyclic torsional simple shear test, and
v) Cyclic triaxial compression test
(iii) With the added mass removed and with the specimen cap, specimen and all apparatus,
determine the resonant frequency, 𝑓𝑛0 . The value of mass polar moment of inertia of the rigid
mass, J0 can be computed using Eq.4.2.
𝐾
𝐽0 = 4𝜋𝑓02 ---------------------- Eq.4.2
𝑛0
Now at resonance cut off the power and record the decay curve for the vibration, From the decay
curve compute the logarithmic decrement for the apparatus, as follows
1 𝐴
𝛿 = 𝑛 𝑙𝑜𝑔𝑒 𝐴1 ---------------- Eq.4.3
2
Under steady state vibrations, the apparatus damping constant, D is given by
𝛿
𝐷 = 𝜋 √𝐾0 𝐽0 -------------------- Eq.4.4
The procedure of obtaining G and ξ has been explained in the following steps:
Calculate the mass density of the specimen, ρ, from Eq. (4.5),
4𝑊
𝜌= -------------------------- Eq.4.5
𝜋𝑑2 𝑙𝑔
Where W =Total weight of specimen
55
Lecture Notes on DSF 2021
l= Length of specimen
d = Diameter of specimen
g =Acceleration due to gravity,
(ii) Calculate the inertia of the specimen about its axis J, as follows:
𝜋
𝐽 = 𝜌 32 𝑑 4 𝑙-------------------- Eq.4.6
(iii) Calculate the system factor, T as follows:
𝐽0 𝐾
𝑇= − 4𝜋2 𝑓02 𝐽-------------------- Eq.4.7
𝐽 𝑛𝑅
Where J0 = Mass polar moment of inertia of the apparatus
K0 = Rotational spring constant,
J = Inertia of the specimen
fnR = Resonant frequency of the complete system.
iv) To measure the torque current constant, Kt excite the apparatus successively at frequencies
(√2, 2), 𝑓𝑛0 √2 and 2𝑓𝑛0 , during the steady state vibration at each of these frequencies measure
the current flowing through the coils, C in amperes, and the displacement amplitude of vibration,
θ in radians. For each frequency compute the torque-current constant Kt as follows
𝐾 𝜃
0
𝐾𝑡 = 𝐶𝑀 ---------------- Eq.4.8
𝑓
Time of loading was defined as the time between the beginning of test and the point at which the
maximum compressive stress is reached (Figure 4.3). The pendulum loading apparatus (Figure 4.
4) utilizes the energy of a pendulum which, when released from a selected height, strikes a spring
connected to the piston rod of a hydraulic lower cylinder. This lower cylinder is connected
hydraulically to an upper cylinder, which is mounted on a loading frame.
57
Lecture Notes on DSF 2021
load from the beam to the specimen (Fig. 4.4 b). A small beam mounted above the
yoke counter-balances the weight of the beam.
Fig. 4.4 (a) Pendulum loading apparatus (Casagrande and Shannon. 1948)
58
Lecture Notes on DSF 2021
For measuring load, a load gage of rectangular or cylindrical shape is used, with four strain gages
mounted on the inside face. For measuring deflection, a thin flexible steel spring cantilever is
used with strain gages mounted on the cantilever, the base of which is clamped to the loading
piston.
A simultaneous plot of stress and strain versus time from an unconfined compression test with a
time of loading of 0.02 s on cambridge clay is shown in Figure 4.5. Similar plots were prepared
for other times of loading on Manchester sand. Using this data, stress-strain plots were obtained
as shown in Figs.6 a and b. In these figures, stress-strain curves for corresponding static tests are
59
Lecture Notes on DSF 2021
also shown. Typical plots of maximum compressive stress versus time of loading (or unconfined
and confined transient tests on Cambridge clay are shown in Fig. 4.7 a and b respectively. A
typical plot in terms of principal stress ratio a failure and time of loading for Manchester sand is
shown in Fig. 4.8.
Fig.4.5: Time Vs stress and strain in an unconfined transient test on Cambridge clay (Casagrande &
Shannon, 1948)
60
Lecture Notes on DSF 2021
Fig. 4.7(a): Maximum Compressive Stress (unconfined) versus Time of Loading for Cambridge clay
Fig.4.7 (b): Maximum Compressive Stress (confined) versus Time of Loading for Cambridge
clay
Fig. 4.8: Maximum principal stress ratio versus time of loading for Manchestor sand
61
Lecture Notes on DSF 2021
Modulus of deformation is defined as the slope of a line drawn from the origin through the point
on the stress-deformation curve and corresponding to stress of one-half the strength. It is found
that in case of clays, modulus of deformation in fast transient tests was about two times that
obtained in static tests. In case of sands, modulus of deformation was found independent of the
time of loading.
4.2.3 Summary of Cyclic Tests
In the preceding sections, various types of laboratory test methods were presented, from which
the fundamental soil properties such as the shear modulus, modulus of elasticity, and damping
ratio are determined. These parameters are used in the design and evaluation of the behavior of
earthen, earth-supported, and earth-retaining structures. As was discussed in the preceding
sections, the magnitudes of G and ξ are functions of the shear strain amplitude γ′. Hence, while
selecting the values of G and ξ for a certain design work, it is essential to know the following:
a. Type of test from which the parameters can be obtained
b. Magnitude of the shear strain amplitude at which these parameters needs to be measured For
example, strong ground motion and nuclear explosion can develop large strain amplitudes
whereas some sensitive equipment such as electron microscopes may be very sensitive to small
strain amplitudes.
Figure 9 provides is a useful reference table for geotechnical engineers; as it gives the amplitude
of shear strain levels, type of applicable dynamic tests, and the area of applicability of these test
results. Despite the fact that laboratory testing is not ideal, it will continue to be important
because soil conditions can be better controlled in the laboratory. Parametric studies necessary
for understanding the soil behaviour of soils under dynamic loading conditions must be
performed in the laboratory conditions. Table 4.3 provides a comparison of the relative qualities
(what property can be measured and what is the degree of quality of the measured property) of
various laboratory techniques for measuring dynamic soil properties. Similarly, Table 4.4 gives a
summary of the different engineering parameters that can be measured in different dynamic or
cyclic laboratory tests
62
Lecture Notes on DSF 2021
Field methods generally depend on the measurement of velocity of waves propagating through
the soil or on the response of soil structure systems to dynamic excitation. The following
methods are in use for determining dynamic properties of soil:
1. Seismic cross-bore hole survey
2. Seismic up-hole survey
3. Seismic down-hole survey
4. Seismic refraction survey
5. Vertical block resonance test
63
Lecture Notes on DSF 2021
6. Horizontal block resonance test
7. Cyclic plate load test
8. Standard penetration test
4.3.1 Seismic Cross-borehole Survey
This method is based on the measurement of velocity of wave propagation from one borehole to
another. Figure 4.10 shows the essentials of seismic cross-hole method outlined by Stoke and
Woods (1972).A source of seismic energy is generated at the bottom of one borehole and the
time of travel of the shear wave from this borehole to another at known distance is measured.
Shear wave velocity is then computed by dividing the distance between the boreholes by the
travel time.
As discussed above, seismic cross-borehole survey can be done using two boreholes one has the
source for causing wave generation and another having geophone for recording travel time.
However, for extensive investigations and better accuracy, three or more boreholes arranged in a
straight line should be used.
In this case the wave velocities can be calculated from the time intervals between succeeding
pairs of holes, eliminating most of the concern over triggering the timing instruments and the
effects of borehole casing and backfilling (Stokoe and Hour, 1978). Also this arrangement of
bore holes in a straight line overcomes problems of site anisotropy by examining one direction
only at a time.
64
Lecture Notes on DSF 2021
al., 1977).This method gives the average value of wave velocity for the soil between the
excitation and the receivers if one receiver is used, or between the receivers.
The major disadvantage in seismic up-hole survey is that it is more difficult to generate waves of
the desired type.
65
Lecture Notes on DSF 2021
67
Lecture Notes on DSF 2021
From the value of Cu obtained from Eq. (3) for the test block of contact area A the value of Cu1
for the actual foundation having contact area A1 may be obtained from Eq. (4.14)
𝐴
𝐶𝑢1 = 𝐶𝑢 √𝐴 ---------------- Eq.4.14
1
The Eq. (4.14) is valid for base areas of foundations up to 10 m2. For areas larger than 10 m2,
the value Cu obtained for 10 m2 is used.
The value of damping ratio ξ is determined using Eq. (4.15) as
𝑓2 −𝑓1
𝜉= ---------------- Eq.4.15
2𝑓𝑛𝑧
𝐴𝑚𝑎𝑥
Where, f1, f2= Two frequencies at which amplitudes is equal to
√2
Amax= Maximum amplitude
fnz = Resonant frequency
The coefficient of elastic uniform compression (Cu) is related to the elastic Young's modulus (E)
by Eq (6) which is in the form of Boussinesq relationship for the elastic settlement of a surface
footing.
𝐸 𝐶𝑠
𝐶𝑢 = 1−𝜇2 ------------------- Eq.4.16
√𝐵𝐿
where 𝜇 = Poisson's ratio
B = Width of base of the block
L = Length of base of the block
Cs = Coefficient depending on 𝐿⁄𝐵 ratio
Barkan (1962) recommended the values of Cs for various 𝐿⁄𝐵 ratios as listed in Table 4.5
Table 4.5: values of Cs for various 𝐿⁄𝐵 ratios
𝐿⁄ Cs
𝐵
1.0 1.06
1.5 1.07
2.0 1.09
3.0 1.13
5.0 1.22
10.0 1.41
68
Lecture Notes on DSF 2021
A vertical vibration test was conducted on a 1.5m x 0.75 m x 0.70 m high concrete block in an
open pit having depth 2.0 which is equal to the anticipated depth of actual foundation. The test
was repeated at different settings (e) of eccentric masses.
The data obtained from the tests are given below:
Sl.No θ(degree) 𝑓𝑛𝑧 Amplitude at Resonance
(Micron)
1 36 41.0 13.0
2 72 40.0 24.0
3 108 34.0 32.0
4 144 31.0 40.0
The soil is sandy in nature having angle of internal friction is 35° and saturated density is 20 kN/
m3). The water table lies at a depth of 3.0 m below the ground surface. Probable size of the actual
foundation is 4.0 x 3.0 x 3.5 m high. Determine the values of Cu, E and G to be adopted for the
design of actual foundation.
Limiting vertical amplitude of the machine is 150 microns.
Sol:
1 Area of Block=1.5x0.75=1.125 m2
Mass of Block=1.125x0.75x2400=1890 kg
Mass of oscillator and motor=100 kg (assumed)
Mass of block, oscillator and motoa=1890+100=1990 kg
4𝜋2 𝑓𝑛𝑧 2𝑚 4𝜋2 ×𝑓𝑛𝑧 2 ×1990
2 𝐶𝑢 = = = 69.84𝑓𝑛𝑧 2 𝑘𝑁/𝑚2
𝐴 1.125×100
The calculated values of 𝐶𝑢 for different observed resonance frequencies are tabulated as
shown in Table 1
𝐸 𝐶𝑠
Cu can be evaluated from Eq.4.6 as𝐶𝑢 = for L/B=2 CS=1.09
1−𝜇 2 √𝐵𝐿
√1.125(1−0.352 )
Assume µ=0.35, 𝐸 = × 𝐶𝑢 =0.854𝐶𝑢 𝑘𝑁/𝑚2
1.09
𝐸 0.854𝐶
𝑢
𝐺 = 2(1+𝜇) = 2(1+0.35)) =0.316𝐶𝑢 𝑘𝑁/𝑚2
For different values of 𝐶𝑢 , E and G values are calculated and tabulated as shown in table 1
3 Correction for confining pressure and area
The mean effective confining pressure ̅̅̅̅at𝜎01 depth of' one -half the width below the centre
of block is given by
(1 + 2𝐾0 )
𝜎01 = ̅̅̅
̅̅̅̅ 𝜎𝑣
3
Where ̅̅̅𝜎𝑣 = ̅̅̅̅
𝜎𝑣1 + ̅̅̅̅
𝜎𝑣2
𝜎𝑣1 =Effective overburden pressure at the depth under consideration
̅̅̅̅
𝜎𝑣2
̅̅̅̅=Increase in vertical pressure due to the weight of block
Assuming that the top 2.0 m soil has a moist unit weight of 18 kN/m3, and the next 1.0 m soil
i.e. up to water table is saturated then
0.70
𝜎𝑣1 = 18 × 2.0 + 20 ×
̅̅̅̅ = 43𝑘𝑁/𝑚2
2
69
Lecture Notes on DSF 2021
4𝑞 2𝑚𝑛√𝑚2 + 𝑛2 + 1 𝑚2 + 𝑛2 + 2 −1
2𝑚𝑛√𝑚2 + 𝑛2 + 1
𝜎𝑣2 =
̅̅̅̅ [ 2 × + 𝑠𝑖𝑛 ]
4𝜋 𝑚 + 𝑛2 + 1 + 𝑚2 𝑛2 𝑚2 + 𝑛2 + 1 𝑚2 + 𝑛2 + 1 + 𝑚2 𝑛2
𝐿/2 1.5/2
𝑚= = 0.7/2=2.14
𝑍
𝐵/2 0.75/2
𝑛= = =1.07
𝑍 0.7/2
Where q=24x0.7=16.8 𝑘𝑁/𝑚2
Substituting for the values of m and n, we get
𝜎𝑣2 = 13.44𝑘𝑁/𝑚2
̅̅̅̅
Now ̅̅̅
𝜎𝑣 = ̅̅̅̅ 𝜎𝑣2 = 43 + 13.44 = 56.44𝑘𝑁/𝑚2
𝜎𝑣1 + ̅̅̅̅
𝐾0 = 1 − 𝑠𝑖𝑛∅ = 0.426
(1 + 2𝐾0 ) 1 + 2 × 0.426
𝜎01 = ̅̅̅
̅̅̅̅ 𝜎𝑣 = 56.44 × = 34.84𝑘𝑁/𝑚2
3 3
𝑘𝑁
For the actual foundation,̅̅̅̅̅=18x2.0+20x1.0+
𝜎𝑣1 (20-10)x0.5= 61 𝑚2
4.0/2
𝑚=
= 1.33
3.0/2
3.0/2
𝑛= = 1.0
3.0/2
𝑘𝑁
𝑞 = 24 × 3.5 = 84 2
𝑚
Substituting the values of m, n and q we get
1+2×0.426 𝑘𝑁
𝜎02
̅̅̅̅=124.76[ ] = 77.01 𝑚2
3
Area of actual foundation =4.0x0.3=12.0 m2 (> 10𝑚2 )
𝐶𝑢2 𝐸2 𝐺2 𝜎02 0.5
̅̅̅̅̅ 𝐴 0.5 77.0 0.5 1.125 0.5
Hence, = = =( ) × ( 1) =( ) ×( ) = 0.4986
𝐶𝑢1 𝐸1 𝐺1 𝜎01
̅̅̅̅̅ 𝐴2 34.84 10
For actual foundation Cu =0.4986x Cu for block
Table 1:
Sl. θ 𝑓𝑛𝑧 Amplitude For Test Block For Actual Foundation
No in at
Cu E G Cu E G
deg. Resonance 4
10 kN/m 2
104kN/m2 4
10 kN/m 2 4
10 kN/m 2 4
10 kN/m 2
104kN/m2
(micron)
1 36 41 13 11.74 10.03 3.71 5.85 5.00 1.85
2 72 40 24 11.17 9.54 3.53 5.57 4.76 1.77
3 108 34 32 7.15 6.11 2.26 3.56 3.05 1.13
4 144 31 40 6.71 5.73 2.12 3.35 2.86 1.06
150×10−6
Strain in Actual foundation= = 0.5 × 10−4
3.0
The value of Cu, E and G corresponding to actual strain level of foundation can be obtained by
interpolation as
0.5−0.427
𝐶𝑢 = [3.56 − (3.56 − 3.35) 0.533−0.427] × 104 = 2.31x104 kN/m2
0.5−0.427
𝐸 = [3.05 − (3.05 − 2.86) 0.533−0.427 ] × 104 =1.97x104 kN/m2
0.5−0.427
𝐺 = [1.13 − (1.13 − 1.06) 0.533−0.427] × 104 = 0.73x104 kN/m2
71
Lecture Notes on DSF 2021
𝑀
Where 𝑟 = 𝑀 𝑚
𝑚0
𝑓𝑛𝑥 = Horizontal resonant frequency of block soil system
𝐴
𝐴0 = 𝑀 -
A = Contact area of block with soil
M = Mass of block, oscillator and soil
𝐼
𝐼0 = 3.46 𝑀 ------------------------------------- Eq.4.20
𝑚𝑜
Mm = Mass moment of inertia of block, oscillator, motor, etc. about the horizontal axis passing
through the centre of gravity of block and perpendicular to the direction of vibration
Mmo = Mass moment of inertia of the block, oscillator; motor etc. about the horizontal axis
passing through centre of contact area of block and soil and perpendicular to the direction of
vibration.
I = Moment of inertia of the foundation contact area about the horizontal axis passing through
the centre of gravity of area and perpendicular to the direction of vibration.
In Eq. (4.19), negative sign is taken when the system vibrates in first mode and positive sign
when the system vibrates in second mode. For the size of the block recommended in IS 5249-
1977 and for first natural frequency, the Eq. (4.19) reduces to
2
𝐶𝜏 = 92.3𝑓𝑛𝑥 ----------------- Eq.4.21 unit of 𝐶𝜏 in this equation is kN/m2
The coefficient of elastic uniform shear (Cτ1) for actual area of foundation (A1) is given by
𝐴
𝐶𝜏1 = 𝐶𝜏 √𝐴 -------------------- Eq.4.22
1
The cyclic plate load test is performed in a test pit dug up to the proposed base level of
foundation. The equipment is same as used in static plate load test. Circular or square bearing
plates of mild steel not less than 25 mm thickness and varying in size from 300 to 750 mm with
grooved bottom are used. The test pit should be at least five times the width of the plate. The
equipment is assembled according to details given in IS 1988-1982. A typical set up is shown in
Figure 4.16.
73
Lecture Notes on DSF 2021
74
Lecture Notes on DSF 2021
𝑃
𝐶𝑢 = 𝑆 (𝑘𝑁⁄𝑚3 )------------------ Eq.4.26
𝑒
Where P=Load Intensity in kN/m2
Se= Elastic rebound corresponding to P in mm
Fig. 4.17: Load Intensity versus Settlement Fig.4.19: Load Intensity versus Elastic rebound
75
Lecture Notes on DSF 2021
Example 1
The plot of q versus s (settlement) obtained from a cyclic plate load test is shown in Figure
below. The area of the plate used for the test was 0.3 m2. Calculate
a. Kplate, and
b. Shear modulus G (assume μ = 0.35).
𝑞𝐴
𝐾𝑝𝑙𝑎𝑡𝑒 = = 142.86 × 0.3
𝑆𝑒
= 42.86𝑀𝑁/𝑚
Following two types of corrections are normally applied to the observed SPT values (N) in
cohesionless soils:
i) Corrections due to dilatancy:
In very fine, or silty, saturated sand, Terzaghi and Peck (1967) recommend that the observed N-
values be
Corrected to N' if N was greater than 15 as
1
𝑁 ′ = 15 + 2 (𝑁 − 15)------------------------Eq.4.30
Bazaraa (1967) recommended the correction as
N' = 0.6N (for N > 15) -------------------- Eq.4.31
This correction is introduced with the view that in saturated dense sand (N > 15); the fast rate of
application of shear through the blows of drop hammer, is likely to induce negative pore
pressures and thus temporary increase in shear strength will occur. This will lead to a N-value
higher than the actual one. Since sufficient experimental evidence is not available to confirm this
correction, many engineers are not applying this correction. However this correction has also
been recommended in IS: 2131-1981.
ii) Correction due to overburden pressure:
On the basis of field tests, corrections to the N-value for overburden effects were proposed by
many investigators (Gibbs and Holtz 1957; Teng 1965; Bazaraa 1967; Peck, Hanson and
Thornburn 1974). The methods which are normally used are:
Bazaraa (1967)
77
Lecture Notes on DSF 2021
4𝑁
For 𝜎̅0 < 75𝑘𝑃𝑎, 𝑁 ′ = ------------------- Eq.4.32
̅̅̅0̅
1+0.04𝜎
4𝑁
For 𝜎̅0 > 75𝑘𝑃𝑎, 𝑁 ′ = ------------------ Eq.4.34
̅̅̅0̅
3.25+0.01𝜎
where 𝜎̅0 =effective over burden pressure, kPa
Peck, Hanson and Thornburm (1974) recommended
2000
𝑁 ′ = 0.77𝑁𝑙𝑜𝑔10 ---------------------- Eq.4.35
̅0
𝜎
Figure 1 gives the correction factor based on Eq.4.35. Use of this figure has been recommended
in IS: 2131-1981. In this figure,
2000
𝐶𝑁 = 𝐶𝑜𝑟𝑟𝑒𝑐𝑡𝑖𝑜𝑓𝑎𝑐𝑡𝑜𝑟 = 0.77𝑁𝑙𝑜𝑔10 𝜎̅ ------------- Eq.4.36
0
Hardin and Black (1968) have given the following factors which influence the shear modulus,
elastic modulus and elastic constants:
78
Lecture Notes on DSF 2021
(i) Type of soil including grain characteristics, grain shape, grain size, grading and mineralogy;
(ii) Void ratio
(iii) Initial average effective confining pressure;
(iv) Degree of saturation;
(v) Frequency of vibration and number of cycles of load
(vi) Ambient stress history and vibration history
(vii) Magnitude of dynamic stress; and
(viii) Time effects
Soil behavior over a wide range of strain amplitudes is nonlinear and on unloading follows a
different stress-strain path forming a hysteresis loop as shown in Figure 4.21. The area inside this
loop represents the energy absorbed by the soil during its deformation and is a measure of the
internal damping within the soil.
At very low strain amplitudes (≪ 0.0001 %) the soil acts essentially as a linear elastic material
with little or no loss of energy. The shear modulus under these conditions is maximum, but as the
strain amplitude is increased, the shear modulus decreases and the damping within the soil
increase.
79
Lecture Notes on DSF 2021
5.0 DESIGN OF MACHINE FOUNDATIONS
5.1 Categories of machine foundations
Reciprocating machines:
It produces periodic unbalanced force and operating frequency is 600rpm. For designing
unbalanced force is taken as varying sinusoidally.
Impact machines:
It produces impact loads at an operating frequency of 60-150 blows/min. Dynamic load attends
the peak within short duration and then die out quickly. Designed as over tuned.
Rotary machines:
These are high speed machines with high operating frequency. Hence the foundations are
designed as under tuned.
Block type
Caisson type
Frame type
81
Lecture Notes on DSF 2021
The second criteria says minimum all-around clearance of 150 mm. must be provided as per IS
codal provision.
A third criterion is that, the foundation block should be placed deep enough on good bearing
strata.
The combined centre of gravity of the machine plus foundation block should be as far below the
top of the foundation as possible.
5.4.2 Vibration Check
Foundation which is having natural frequency either much higher or lower than the operating
frequency of the machine, is called under tuned or over tuned respectively.
If the ratio of operating frequency to the natural frequency is less than or equal to 0.5 that can be
designed as Under tuned criterion.
If the operating frequency is much higher than the natural frequency and frequency ratio
𝜔
𝜉= > 2 for important machine and
𝜔𝑛
𝜔
𝜉 = 𝜔 > 1.5for less important machine
𝑛
For both types of machine design criteria is Over tuned.
• To design machine for which operating frequency is very high, over tuned criteria is used
because an under tuned type of foundation design result will provide negative value of
mass or no mass because of the value of k which has to be excessively high.
• To design machine for which operating frequency is low, under tuned criteria is used.
• However for the range of say 1000 rpm or even in the range of 600 rpm, it is always
better to check your design for both over tuned and under tuned.
5.4.3 Displacement Criteria
• The amplitude of permissible dynamic displacement should be less than or equals to 0.2
mm, If it exceeds, foundation is to be redesigned.
• The permissible displacement should be checked using Richart’s chart. So, that it should
not become annoying to the workers or adjacent structures.
• Y-axis of the chart shows dynamic displacement amplitude and X axis represents
operating frequency of the machine.
82
Lecture Notes on DSF 2021
5.5 LINEAR ELASTIC WEIGHTLESS SPRING MSD MODEL
Vertical
𝐾𝑧
𝜔𝑛𝑧 =√ Vibration
𝑚
Horizontal
𝐾𝑥
𝜔𝑛𝑥 =√ Vibration
𝑚
Rocking Mode
𝐾𝑛∅
𝜔𝑛∅ =√
𝑀𝑚0
Yawing Mode
𝐾𝑛𝜑
𝜔𝑛𝜑 =√
𝑀𝑚𝑧
𝑃𝑥
Maximum Amplitude of Vibration = 2 −𝜔2 ) ------------------ Eq.5.14
𝑚(𝜔𝑛𝑧
83
Lecture Notes on DSF 2021
5.4.3 LINEAR ELASTIC THEORY
(Based on Elastic Half Space Theory)
In 1904, Lamb studied the problem of vibration of single vibrating force acting at a point on the
surface of an elastic half-space. This study included cases in which the oscillating force R acts in
the vertical direction and in the horizontal direction, as shown in Figure 5.3 a and b. This is
generally referred to as the dynamic Boussinesq problem.
In 1936, Reissner analyzed the problem of vibration of a uniformly loaded flexible circular area
resting on an elastic half-space. The solution was obtained by integration of Lamb’s solution for
a point load. Based on Reissner’s work, the vertical displacement at the center of the flexible
loaded area (Figure 5.4 a) can be given by
𝑄𝑂 𝑒 𝑖𝜔𝑡
𝑍= (𝑓1 + 𝑖𝑓2 )----------------------- Eq.5.15
𝐺𝑟0
Where Q0=amplitude of the exciting force acting on the foundation
Z= periodic displacement at the centre of the loaded area
𝜔= circular frequency of the applied load
𝑟0 = radius of the loaded area
G= shear modulus of the soil
𝑓1 , 𝑓2 =Reissner’s displacement functions
84
Lecture Notes on DSF 2021
Fig.5.4: a) Vibration of uniformly loaded circular flexible area, b) Flexible circular area
subjected to force vibration
The displacement functions𝑓1 and 𝑓2 are related to the Poisson’s ratio of the medium and the
frequency of the exciting force. Now, consider a flexible circular foundation of weight W (mass
= m = W/g) resting on an elastic half-space and subjected to an exciting force of magnitude of
(𝑄0 𝑒 𝑖(𝜔𝑡+𝛼) as shown in Figure 5.4b. (Note: α is the phase difference between the exciting force
and the displacement of the foundation.)
Using the displacement relation given in Eq. (5.15) and solving the equation of equilibrium of
force, Reissner obtained the following relationships:
𝑄
𝐴𝑧 = 𝐺𝑟0 𝑍---------------------- Eq.5.16
0
Where Az = the amplitude of vibration
Z= dimensionless amplitude
𝑓2 +𝑓2
=√(1−𝑏𝑎2 𝑓 1)2 +(𝑏𝑎
2
2 𝑓 ))2--------------------- Eq.5.17
0 1 0 2
Vs= velocity of shear wave in the elastic material on which the foundation is resting
The classical work of Reissner was further extended by Quinlan (1953) and Sung (1953). As
mentioned before, Reissner’s work related only to the case of flexible circular foundations where
the soil reaction is uniform over the entire area (Figure 5.5a). Both Quinlan and Sung considered
the cases of rigid circular foundations, the contact pressure of which is shown in Figure 5.3b,
flexible foundations (Figure 5.5a), and the types of foundations for which the contact pressure
distribution is parabolic, as shown in Figure 5.5c. The distribution of contact pressure q for all
three cases may be expressed as follows.
For flexible circular foundations
𝑄0 𝑒 𝑖(𝜔𝑡+𝛼)
𝑞= for 𝑟 ≤ 𝑟0 --------------------- Eq.5.20
𝜋𝑟02
For rigid circular foundations
𝑄0 𝑒 𝑖(𝜔𝑡+𝛼)
𝑞= (𝑓𝑜𝑟𝑟 ≤ 𝑟0 )--------------------- Eq.5.21
2𝜋𝑟0 √𝑟02−𝑟 2
Quinlan derived the equations only for the rigid circular foundation; however, Sung presented
the solutions for all the three class described. For all cases, the amplitude of motion can be
expressed in a similar form to Eqs. (5.2 to 5.5). However, the displacement functions 𝑓1 and
𝑓2 will change, depending on the contact pressure distribution.
𝜌
Where 𝑎0 = 𝜔𝑟0 √𝐺 -------------------------------- Eq.5.25
𝐺𝑎02
𝜔2 = -------------------------------------- Eq.5.26
𝜌𝑟02
Substituting Eq. 5.26 into Eq. 5.24 we get
𝑚1𝑒𝑎02 𝑓12 +𝑓22 𝑚1 𝑒
𝐴𝑧 = √ 2 = 𝑍 ′ ------------- Eq.5.27
𝜌𝑎03 (1−𝑏𝑎02 𝑓1) +(𝑏𝑎02𝑓2 ))2 𝜌𝑟03
𝑓2 +𝑓2
Where 𝑍 ′ = 𝑎02 √(1−𝑏𝑎2𝑓 1)2 +(𝑏𝑎
2
2 𝑓 ))2-------------------- Eq.5.28
0 1 0 2
Figures 5.7 and 5.8 show the plots of the variation of the dimensionless amplitude with a0
(Richart, 1962) for rigid circular foundations (for μ = Poisson’s ratio = 0.25 and b = 5, 10, 20,
and 40).
87
Lecture Notes on DSF 2021
Fig. 5.7: Plot of Z versus a0 for rigid circular Fig. 5.8: Variation of Z′ with a0 for
foundation , Richart, 1962) rigid circular foundation (redrawn
after Richart, 1962)
88
Lecture Notes on DSF 2021
Fig. 5.9: Effect of contact pressure distribution variation of Z’ with a0 (redrawn after Richart and
Whitman, 1967)
Fig.5.10: Effect of Poisson’s ratio on the variation of Z’ with a0 (redrawn after Richart and
Whitman, 1967)
90
Lecture Notes on DSF 2021
4𝐺𝑟0
𝐾𝑧 = 𝑠𝑡𝑎𝑡𝑖𝑐𝑠𝑝𝑟𝑖𝑛𝑔𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑓𝑜𝑟𝑟𝑖𝑔𝑖𝑑𝑐𝑖𝑟𝑐𝑢𝑙𝑎𝑟𝑓𝑜𝑢𝑛𝑑𝑎𝑡𝑖𝑜𝑛 = --------- Eq.5.32
1−𝜇
And
3.4𝑟02
𝐶𝑧 = √𝐺𝜌---------------------------- Eq.5.33
1−𝜇
Fig.5.12: Plot of F1 and -F2 against a0 for rigid circular foundation subjected to vertical
vibration (after Lysmer and Richart, 1966)
In Eqs. (5.32) and (5.33) the relationships for Kz and Cz are frequency independent. Equations
(5.31 to 5.33) are referred to as Lysmer’s analog.
91
Lecture Notes on DSF 2021
the mathematical expressions presented in earlier section. The general procedure is outlined
next.
A. Resonant Frequency
1. Calculation of natural frequency. as
1 𝐾 1 4𝐺𝑟
𝑓𝑛 = 2𝜋 √ 𝑚𝑧 = 2𝜋 √(1−𝜇)𝑚
0
-------------------- Eq.5.34
2. Calculation of damping ratio ξ. As given
0 4𝐺𝑟 4𝐵𝑧 𝜌𝑟03
Critical damping 𝐶𝑐 = 2√𝐾𝑧 𝑚 = 2√(1−𝜇) 𝑚, substituting for m as 𝑚 =
(1−𝜇)
0 4𝐺𝑟
𝑧 0 𝐵 𝜌𝑟 3
0 8𝑟 2
= 4√((1−𝜇) ) ((1−𝜇) ) = (1−𝜇) √𝐺𝜌𝐵𝑧 ------------------------ Eq.5.35
3.4𝑟2
0
𝐶 1−𝜇
√𝐺𝜌 0.425
Now , 𝜉 = 𝐶 = 8𝑟02 = -------------------------- Eq.5.36
𝐶 √𝐺𝜌𝐵𝑧 √ 𝐵𝑧
(1−𝜇)
3. Calculation of the resonance frequency (that is, frequency at maximum displacement). For
constant force-type excitation,
𝑓𝑚 = 𝑓𝑛 √1 − 2𝜉 2
1 04𝐺𝑟 0.425 2
= 2𝜋 √(1−𝜇)𝑚 × √1 − 2( ) --------------------- Eq.5.37
√𝐵𝑧
It has also been shown by Lysmer that, for Bz ≥ 0.3, the following approximate relationship can
be established:
1 𝐺 1 𝐵𝑧 −0.36
𝑓𝑚 = 2𝜋 √(𝜌 ) (𝑟 ) × √ ----------------------- Eq.5.38
0 𝐵𝑍
92
Lecture Notes on DSF 2021
𝑄0
⁄𝐾
𝑧
𝐴𝑧 = ------------------- Eq.5.43
2
𝜔2 𝜔 2
√(1− 2 ) +4(𝜉𝑧 )2 ( )
𝜔𝑛 𝜔𝑛
𝐴 𝜔 𝜔
Figure 5.13 shows the plot of 𝑄0 𝑧 versus (𝜔 ). So, with known values of 𝜉𝑧 and (𝜔 ), one can
⁄𝐾 𝑛 𝑛
𝑧
𝐴𝑧
determine the value of 𝑄0 and, from that, Az can be obtained.
⁄𝐾
𝑧
In a similar manner, for rotating mass-type excitation, Eq. (2.95) can be used to determine the
amplitude of vibration, or
𝑚 𝑒 2
( 𝑒 ⁄𝑚 )(𝜔⁄𝜔𝑛 )
𝐴𝑧 = ----------------------------- Eq.5.44
2
𝜔2 𝜔 2
√(1− 2 ) +4(𝜉𝑧 )2 ( )
𝜔𝑛 𝜔𝑛
Fig.5.13: Plot of various non-dimensional parameters against ( 𝜔⁄𝜔𝑛 ) for constant force-type
vibrator (Note: ξ = ξz for vertical vibration, ξ = ξθ for rocking, ξ = ξx for sliding; ξ = ξα for
torsional vibration.)
93
Lecture Notes on DSF 2021
Fig.5.14: Plot of various non-dimensional parameters against ( 𝜔⁄𝜔𝑛 ) for rotating mass-type
vibrator (Note: ξ = ξz for vertical vibration, ξ = ξθ for rocking, ξ = ξx for sliding; ξ = ξα for
torsional vibration.)
The procedure here described relates to a rigid circular foundation having a radius of r0. If a
foundation is rectangular in shape with length L and width B, it is required to obtain an
equivalent radius, which can then be used in the preceding relationships as discussed in above.
This can be done by equating the area of the given foundation to the area of an equivalent circle.
Thus,
𝜋𝑟02 = 𝐿𝐵
where r0 = radius of the equivalent circle.
It is obviously impossible to eliminate vibration near a foundation.
However, an attempt can be made to reduce the vibration problem as much as possible. Richart
(1962) compiled guidelines for allowable vertical vibration amplitude for a particular frequency
of vibration, and this is given in Figure 5.15.
94
Lecture Notes on DSF 2021
Fig.5.15:
The data presented in Figure 5.15 refer to the maximum allowable amplitudes of vibration.
These can be converted to maximum allowable accelerations by
Maximum acceleration = (maximum displacement)𝜔2
Arnold, Bycroft, and Wartburton (1955) have provided theoretical solutions for sliding vibration
of rigid circular foundation (Figure 5.16) acted on by a force, 𝑄 = 𝑄0 𝑒 𝑖𝜔𝑡 . Hall (1967)
95
Lecture Notes on DSF 2021
developed the mass-spring-dashpot analog for this type of vibration. According to this analog,
the equation of motion of the foundation can be given in the form
𝑚𝑥̈ + 𝐶𝑥 𝑥̇ + 𝐾𝑥 𝑥 = 𝑄0 𝑒 𝑖𝜔𝑡 ------------------------ Eq.5.45
where m= mass of the foundation
Same figures 5.13 & 5.14 are used to calculate various non-dimensional parameters against
(𝜔⁄𝜔𝑛 ) for constant force and rotating mass-type vibrator respectively.
5.8 TORSIONAL VIBRATION OF FOUNDATIONS
Figure 5.21a shows a circular foundation of radius r0 subjected to a torque 𝑇 = 𝑇0 𝜃 𝑖𝜔𝑡 about an
axis z-z. The vibration problem of this type was solved by Reissner (1937) solved considering a
linear distribution of shear stress 𝜏𝑧𝜃 (shear stress zero at center and maximum at the periphery
of the foundation), as shown in Figure 5.17b which represents the case of a flexible foundation.
97
Lecture Notes on DSF 2021
In 1944 Reissner and Sagoli solved the same problem for the case of a rigid foundation
considering a linear variation of displacement from the center to the periphery of the foundation.
Similar to the cases of vertical, rocking, and sliding modes of vibration, the equation for the
torsional vibration of a rigid circular foundation can be written as
Where, Jzz= mass moment of inertia of the foundation about the axis z-z
Cα = dashpot coefficient for torsional vibration
16
Kα = static spring constant for torsional vibration = 𝐺𝑟03
3
α = rotation of the foundation at any time due to the application of a torque 𝑇 = 𝑇0 𝑒 𝑖𝜔𝑡
The damping ratio ξα for this mode of vibration has been determined as (Richart, Hall, and
Wood, 1970) given below
0.5
𝜉𝛼 = ----------------------- Eq.5.58
1+2𝐵𝛼
Where
𝐽
Bα = the dimensionless mass ratio for torsion at vibration= 𝜌𝑟𝑧𝑧5---------------- Eq.5.59
0
Calculation Procedure for Foundation Response Using Eq. (5.51)
Resonant Frequency
1. Calculate the natural frequency fn
1 𝐾
𝑓𝑛 = 2𝜋 √𝐽 𝛼 --------------------- Eq.5.60
𝑧𝑧
2. Calculate the damping ratio Bα using Eq. (5.59) and damping ratio 𝜉𝛼 by Eq. 5.58
3. For constant force excitation (that is, T0 = constant), calculate
𝑓𝑚 = 𝑓𝑛 √1 − 2𝜉𝛼2 --------------------------Eq.5.61
98
Lecture Notes on DSF 2021
4. For rotating mass type excitation, calculate
𝑓𝑛
𝑓𝑚 = --------------------------------------Eq.5.62
√1−2𝜉𝛼2
For constant force excitation, calculate 𝜔⁄𝜔𝑛 and then refer to Figure 5.13 to obtain α/(T0/Kα).
For rotating mass-type excitation, calculate 𝜔⁄𝜔𝑛 and then refer to Figure 5.14 to obtain α/[m1
e(x/2)/Jzz].
For a rectangular foundation with dimensions B × L, the equivalent radius may be given by
𝐵𝐿(𝐵 2+𝐿2 )
𝑟0 = √ -------------------- Eq.5.67
6𝜋
The torsional vibration of foundations is uncoupled motion and hence can be treated
independently of any vertical motion. Also, Poisson’s ratio does not influence the torsional
vibration of foundations
A theoretical solution for foundations subjected to rocking vibration was presented by Arnold,
Bycroft, and Wartburton (1955) and Bycroft (1956). Rocking mode of vibration for rigid
circular foundations is shown in Figure 5.18.
99
Lecture Notes on DSF 2021
A mass-spring-dashpot model for rigid circular foundations was developed by Hall (1967) in the
same manner as Lysmer and Richart (1966) developed for vertical vibration. According to Hall,
the equation of motion for a rocking vibration can be given as
3(1−𝜇) 𝐼
Inertia ratio, 𝐵𝜃 = ( ) 𝜌𝑟05--------------------------- Eq.5.72
8 0
100
Lecture Notes on DSF 2021
A. Resonant Frequency
1 𝐾
1. Calculate the natural frequency:𝑓𝑛 =
2𝜋
√ 𝐼 𝜃 ------------------------ Eq.5.73
0
𝑓𝑛
𝑓𝑚 = -----------(for rotating mass-type excitation)---------------- Eq.5.77
√1−2𝜉𝜃2
𝑚1 𝑒𝑍 ′ 1
𝜃𝑟𝑒𝑠𝑜𝑛𝑎𝑛𝑐𝑒 = --------------------- Eq.5.79
𝐼𝜃
2𝜉𝜃 √1−𝜉𝜃2
101
Lecture Notes on DSF 2021
𝑚 𝑒𝑧 , 2
( 1 ⁄𝐼 )(𝜔⁄𝜔𝑛 )
0
𝜃= ---------------------- Eq.5.81
2
𝜔2 𝜔 2
√(1− 2 ) +4(𝜉𝛼 )2 ( )
𝜔𝑛 𝜔𝑛
In the case of rectangular foundation, the preceding relationships can be used by determining the
equivalent radius as
4 𝐵𝐿3
𝑟0 = √ 3𝜋 -------------- Eq.5.82
Q0=27 kN (not frequency dependent). The operating frequency is 150 cpm. For the soil
supporting the foundation, unit weight = 19.5 kN/m3, shear modulus = 45000 kPa. , and Poisson’
ratio = 0.3. Determine:
a. resonant frequency,
b. the amplitude of vertical vibration at resonant frequency, and
c. the amplitude of vertical vibration at the operating frequency
Sol:
The machine imparts a vertical vibrating force Q = Q0 𝑠𝑖𝑛𝜔𝑡 where Q0=27 kN
The operating frequency = 150 cpm=2.5 Hz
Equivalent radius 𝑟0 = 1.25𝑚
Total weight of the machine and the foundation =270 kN
Mass ratio
(1−𝜇)𝑚 (1−𝜇)𝑊 (1−0.3)×270
𝐵𝑧 = = = 4×19.5×1.253 =1.24
4𝜌𝑟03 4𝛾𝑟03
Damping factor
102
Lecture Notes on DSF 2021
0.425 0.425
𝜉𝑍 = = = 0.382
√𝐵𝑧 √1.24
Spring constant
4𝐺𝑟0 4×4500×1.25
𝐾𝑧 = = =32142.86 kN/m2
1−𝜇 1−0.3
10 4𝐺𝑟 1 4×4500×1.25
Natural frequency 𝑓𝑛 = 2𝜋 √(1−𝜇)𝑚 = 2𝜋 √(1−0.3)×270⁄ =5.44 Hz
9.81
𝑄0 27⁄
⁄𝐾
𝑧 32142.86
𝐴𝑧 = = =
2 2 2
𝜔2 𝜔 2 2
√(1− 2 ) +4(𝜉𝑧 )2 ( ) √(1− 2.5 2 ) +4(0.382)2 ( 2.5 )
𝜔𝑛 𝜔𝑛 5.44 5.44
Example No.2:
A radar antenna foundation is shown below. For torsional vibration of the foundation, given
T0 = 250 kN-m (due to inertia)
T0 = 83 kN-m (due to wind)
Mass moment of inertia of the tower about the axis z-z = 13 × 106 kg·m2, and the unit weight of
concrete used in the foundation = 24 kN/m3. Calculate
a) the resonant frequency for torsional mode of vibration; and
b) angular deflection at resonance.
Damping factor
104
Lecture Notes on DSF 2021
105
Lecture Notes on DSF 2021
5.10 DYNAMIC BEARING CAPACITY OF SHALLOW FOUNDATION
During the application of single pulse dynamic loads which may be in vertical or horizontal
directions, the foundation may get excessive settlement. Horizontal dynamic loads on
foundations are due mostly to earthquakes. These types of loading may induce large
permanent deformation in foundations. Isolated column footings, strip footings, mat footings,
and even pile foundations all may fail during seismic events. Such failures are generally
attributed to liquefaction. However, a number of failures have occurred where field
conditions indicate there was only partial saturation or a dense soil and therefore liquefaction
alone is a very unlikely explanation. Rather, the reason for the seismic settlements of these
foundations seems to be that the bearing capacity was reduced (Richards, Elms and
Budhu,1993).
During the analysis of the time dependent motion of a foundation subjected to dynamic
loading or estimating the bearing capacity under dynamic conditions several factors need to
be considered.
Most important of these factors are
a) Nature of variation of the magnitude of the loading pulse,
b) Duration of the pulse, and
c) Strain-rate response of the soil during deformation
However, when load is applied rapidly to a foundation to cause failure, the ultimate bearing
capacity changes by somewhat. This fact has been shown experimentally by Vesic, Banks,
and Woodward (1965), who conducted several laboratory model tests with a 101.6 mm
diameter rigid rough model footing placed on the surface of a dense river sand (i.e., Df = 0),
both dry and saturated. The rate of loading to cause failure was varied in a range of 2.54 ×
10-4 mm/s to over 254 mm/s. Hence, the rate was in the range of static (2.54 × 10-4 mm /s) to
impact (254 mm/s) loading conditions. All but the four Based on the experimental results
available, the following general conclusions regarding the ultimate dynamic bearing capacity
of shallow foundations in sand can be drawn:
1. For a foundation resting on sand and subjected to an acceleration level of amax ≤ 13g, it is
possible for general shear type of failure to occur in soil (Heller, 1964).
106
Lecture Notes on DSF 2021
2. For a foundation on sand subjected to an acceleration level of amax > 13g, the nature of soil
failure is by punching (Heller, 1964).
3. The difference in the nature of failure in soil is due to the inertial restrain of the soil
involved in failure during the dynamic loading. The restrain has almost a similar effect as the
overburden pressure as observed during the dynamic loading which causes the punching
shear type failure in soil.
4. The minimum value of the ultimate dynamic bearing capacity of shallow foundations on
dense sands obtained between static to impact loading range can be estimated by using a
friction angle ∅𝑑𝑦 , such that (Vesic, 1973)
∅𝑑𝑦 = ∅ − 20 ------------------ Eq.5.84
The value of ∅𝑑𝑦 can be subsequently used to find various bearing capacity factors.
However, if the soil strength parameters with proper strain rate are known from laboratory
testing, they should be used instead of the approximate equation.
5. The increase of the ultimate bearing capacity at high loading rates is due to the fact that the
soil particles in the failure zone do not always follow the path of least resistance. This results
in at higher shear strength of soil, which leads to a higher bearing capacity.
6. In the case of foundations resting on loose submerged sands, transient liquefaction effects
may exist (Vesic, 1973). This may results in unreliable prediction of ultimate bearing
capacity.
7. The rapid increase of the ultimate bearing capacity in dense saturated sand at fast loading
rates is due to the development of negative pore water pressure in the soil.
The dynamic bearing capacity problem attracted attention of the investigators in 1960 when the
performance of foundations under transient loads became of concern to the engineering
profession (Wallace, 1961; Cunny and Sloan, 1961; Fisher, 1962; Johnson and Ireland, 1963;
Mckee and shenkman 1962: White, 1964; Chummar, 1965; Triandafilidis, 1965).
All analytical approaches are based on the assumption that soil rupture under transient
loads occurs along a static rupture surface. In this section the salient features of the analysis
developed by Triandafilidis (1965) and Wallace (1961) for transient vertical load; and by
Chummar (1965) for transient horizontal load have been presented.
107
Lecture Notes on DSF 2021
5.10.2 Triandafilidis’s Solution:
Fig. 5.19: Illustrations of mode of failure and dynamic equilibrium of moving soil mass
TriandafiIidis's Analysis
Triandafilidis (1965) has presented a solution for dynamic response of continuous surface
footing supporting by saturated cohesive soil (φ=0 condition) and subjected to vertical transient
load. The analysis is based on the following assumptions:
(i) The failure surface of soil is cylindrical for evaluation of bearing capacity under static
condition (Fig.5.19).
(ii) The saturated cohesive soil (φ= 0) behaves as a rigid plastic material (Fig. 5.20).
(iii) The forcing function is assumed to be an exponentially decaying pulse (Fig.5.21)
(iv) The influence of strain rate on the shear strength is neglected.
(v) The dead weight of the foundation is neglected.
108
Lecture Notes on DSF 2021
Fig. 5.20 Assumed Stress-strain relationship Fig. 5.21 Transient Vertical Load
Analysis
Let the transient stress pulse be expressed in the form
𝑞𝑑 = 𝑞0 𝑒 −𝛽𝑡 = 𝜆𝑞𝑢 𝑒 −𝛽𝑡 -------- Eq.5.85
Where 𝑞𝑑 = Stress at any time t
𝛽= Decaying function
𝑞𝑢 = Static bearing capacity of continuous footing
𝑞0 =instantaneous peak intensity of stress pulse
𝑞
𝜆= Over load factor = 𝑞0
𝑢
The rupture surface is shown in Fig. 5.19 with centre of rotation at point0 located at a height of
0.43 B above the ground surface.
The equation of motion is written by equating the moment of the disturbing and restoring forces
taken about the point O. The only disturbing and restoring force is an externally applied dynamic
pulse. The restoring forces consist of shearing resistance along the rupture surface, the inertia of
the soil mass taken in the motion and the resistance caused by the displacement of centre of
gravity of soil mass.
Disturbing moment 𝑀𝑑𝑝 due to applied dynamic pulse is given as
1
𝑀𝑑𝑝 = 2 𝑞𝑑 𝐵2 --------- Eq. 5.86
109
Lecture Notes on DSF 2021
The static bearing capacity of a continuous footing along the failure surface (Fellenius, 1948) is given as
𝑞𝑢 = 5.54𝑐𝑢
Where cu is the un-drained shear strength of soil
Now Resisting moment Mrs due to shear strength is taken as
1
𝑀𝑟𝑠 = 2 𝑞𝑢 𝐵2 ------------- Eq.5.87
Due to the application of pulse, the soil mass is subjected to an acceleration. So the resisting
moment 𝑀𝑟𝑖 due to the rigid body motion of the failed soil mass is given as
𝑀𝑟𝑖 = 𝐽0 𝜃̈ ------- Eq.5.88
𝑊𝐵 2
𝐽0 =Polar mass moment of inertia= 1.36𝑔
The displaced position of the soil mass generates a restoring moment 𝑀𝑟𝑤 which may be expressed as
By equating the moments of disturbing forces to those of the restoring forces, the following
equation of motion is obtained
Mdp = Mrs + Mri + Mrw ------ Eq.5.91
Substituting for moments and rearranging, we get
3𝑔 0.68𝑔
𝜃̈ + 𝜋𝐵 𝜃 = [ 𝑊 ] 𝑞𝑢 [𝜆𝑒 −𝛽𝑡 − 1]------ Eq.5.92
Equation (5.92) is a second order, non-homogeneous, linear differential equation with constant
coefficients. The natural frequency and the time period of the system are given by
3𝑔
𝜔𝑛 = √𝜋𝐵--------- Eq.5.93
110
Lecture Notes on DSF 2021
The above re1ation can be used to trace the history of motion of the foundation. For
determination of the maximum angular deflection𝜃, Eq. (5.95) can be differentiated with respect
to time. Thus
𝑊 2𝑇𝜋 𝛽2 𝑇 2 2𝜋𝑡 𝛽𝜆𝑇 2𝜋𝑡 𝛽𝜆𝑇
(𝜃̇ ) = [{𝜆 − 1 − } sin ( )+ cos ( )− 𝑒 −𝛽𝑡 ]------- Eq. 5.96
0.68𝑔𝑞𝑢 4𝜋2 +𝛽2 𝑇 2 4𝜋2 𝑇 2𝜋 𝑇 2𝜋
For obtaining the critical time 𝑡 = 𝑡𝑐𝑟 which corresponds to 𝜃 = 𝜃𝑚𝑎𝑥 the right-hand side of Eq.
2𝜋𝑇
(5.96) is equated to zero. Since 4𝜋2 +𝛽2 𝑇 2 cannot be zero,
By using small increments of time t in Eq. (5.97), the value of tcr can be obtained. This value of t
= tcr can then be substituted in to Eq. (5.96) with known values of λ,𝛽𝑎𝑛𝑑𝐵 to obtain
𝑊
𝜃𝑚𝑎𝑥 = 𝐾,is the dynamic load factor.
0.68𝑔𝑞𝑢
Figures 5.22 to 5. 24 give the values of K (s2) for B = 0.6, 1.5 and 3.0 m, respectively, with λ =
1-5 and 𝛽= 0-50 s-1
Fig.5.22 Relationship between overload ratio and dynamic load factor for continuous footings
0.6 m wide
111
Lecture Notes on DSF 2021
(i) The failure of the footing occurs with the application of horizontal dynamic load
acting at a certain height above the base of the footing.
(ii) The resulting motion in the footing is of a rotatory nature. The failure surface is a
logarithmic spiral with its centre on the base corner of the footing, which is also the
centre of rotation as shown in (Fig.5.25).
(iii) The rotating soil mass is considered to be a rigid body rotating about a fixed axis.
(iv) The soil exhibits rigid plastic, stress-strain characteristics.
Fig.5.25: Transient horizontal load on a continuous strip footing resting on ground surface
112
Lecture Notes on DSF 2021
Analysis: The static bearing capacity of the footing is calculated by assuming that the
footing fails when acted upon by a vertical static load, which causes rotation of the logarithmic
spiral failure. The ultimate static bearing capacity qu is given by:
1
𝑞𝑢 = 𝑐𝑁𝑐 + 2 𝛾𝑁𝛾 𝐵 Eq. 5.98
Where c =Cohesion
B= Footing width and equal to the initial radius of spiral curve.
γ=Unit weight of the soil.
𝑁𝑐 and𝑁𝛾 are bearing capacity factors for the assumed type of failure.
Now considering moment of the forces about0, the centre of rotation:
(𝑒 2𝜋𝑡𝑎𝑛∅ −1)
Moment due to cohesion c, 𝑀𝑅𝐶 = 𝛹𝑐𝐵2 , where = 𝛹----------- Eq.(5.99)
2𝑡𝑎𝑛∅
𝑡𝑎𝑛∅(𝑒 3𝜋𝑡𝑎𝑛∅ +1)
Moment due to weight W of soil wedge, 𝑀𝑅𝑊=𝜀𝛾𝐵3 , where = 𝜀-----Eq.(5.100)
9𝑡𝑎𝑛 2 ∅+1
∅is the angle of internal friction
𝐵2
Moment due to qu about point O is given as qu 2
𝐵2
Under equilibrium condition, we get qu 2 =𝑀𝑅𝐶 + 𝑀𝑅𝑊 , which gives
𝑐 2𝛾𝐵𝑡𝑎𝑛∅(𝑒 3𝜋𝑡𝑎𝑛∅ +1)
𝑞𝑢 = 𝑡𝑎𝑛∅ (𝑒 2𝜋𝑡𝑎𝑛∅ − 1) + --- Eq. (5.101)
9𝑡𝑎𝑛 2 ∅+1
Combining Eq.(5.98) and (5.101), yields
4𝑡𝑎𝑛∅(𝑒 3𝜋𝑡𝑎𝑛∅ +1)
𝑁𝛾 = --------------- Eq. (5.102)
9𝑡𝑎𝑛 2 ∅+1
𝑒 2𝜋𝑡𝑎𝑛∅ −1
𝑁𝑐 = --------------------------------- Eq. (5.103)
9𝑡𝑎𝑛 2 ∅+1
By considering a suitable factor of safety F, the static vertical force on the foundation per unit
length can be given as
𝐵 1
𝑄 = 𝐹 (𝑐𝑁𝑐 + 2 𝛾𝐵𝑁𝛾 )------------------------- Eq. (5.104)
The variation of dynamic force in the above analysis is considered as
𝑄𝑑(max) = 𝜆𝑄--------------------------------- Eq. (5.105)
113
Lecture Notes on DSF 2021
Where, 𝑄𝑑(max) is the maximum value of horizontal transit load per unit length acting at height H
above base of the footing and 𝜆is over load factor.
For considering of the dynamic equilibrium of the foundation with the horizontal transient load,
the moment of each of the forces (per unit length) about the centre of the log-spiral needs to be
considered:
1
1. Moment due to the vertical force Q,𝑀1 = 2 𝑄𝐵-------------------- Eq. (5.106)
2. Moment due to the horizontal force Qd at any time t,
𝑄𝑑(𝑚𝑎𝑥)𝐻𝑡 𝑀𝑑(𝑚𝑎𝑥)𝑡
𝑀2 = 𝑄𝑑 𝐻 = = -------------------------------------- Eq. (5.107)
𝑡𝑑 𝑡𝑑
3. Moment due to the cohesive force acting along the failure surface is given by Eq.(5.88)
4. Moment due to weight of soil mass in the failure wedge is given by Eq. (5.89)
5. Moment of the force due to displacement of the centre of gravity of the failure wedge (as
shown in Fig.5.26) from its initial position:
M3 =W X ---------------------------- Eq. (5.108)
Where W is the weight of the failure wedge, and given by
(𝑒 2𝜋𝑡𝑎𝑛∅ −1)𝛾𝐵 2
𝑊= ------------------------- Eq. (5.109)
4𝑡𝑎𝑛∅
∆𝑋̅ = 𝑅𝑐𝑜𝑠(𝜂 − 𝛼 ) − 𝑅𝑐𝑜𝑠𝜂------------------ Eq. (5.110)
When α is small, Eq.(5.110) can be written as ∆𝑋̅ = (𝑅𝑠𝑖𝑛𝜂)𝛼-------------- Eq. (5.111)
But 𝑅 = √𝑥̅ 2 + 𝑧̅ 2 and
−4𝐵𝑡𝑎𝑛 2 ∅(𝑒 3𝜇𝑡𝑎𝑛∅ +1)
𝑥̅ = (9𝑡𝑎𝑛2 ∅+1)(𝑒 2𝜋𝑡𝑎𝑛∅ −1)-------------------------------------- Eq. (5.112)
Example No.2.1
A 2.5m wide continuous surface footing is subjected to a horizontal transient load of duration
0.4s applied at a height of 4.0 m from the base of footing. The properties of the soil are
γ=17kN/m3, c=30kN/m3 and φ =32°.Determine the value of the maximum horizontal load that
can be applied on the footing. Also compute the angular rotation at time equal to 0.6 s.
Sol: Given that
γ=17kN/m3, c =30kN/m3and φ =32°.H= 4.0 m, td=0.4 s
(𝑒 2𝜋𝑡𝑎𝑛∅ −1) 4𝑡𝑎𝑛∅(𝑒 3𝜋𝑡𝑎𝑛∅ +1)
Nc= =79.4, Nγ= =200
𝑡𝑎𝑛∅ 9𝑡𝑎𝑛 2 ∅+1
𝐵 1
Q= 2 (𝑐𝑁𝑐 + 2 𝛾𝐵𝑁𝛾 )=8290 Taking a suitable value of factor of safety as 2.0
i) Determination of various parameters
(𝑒 2𝜋𝑡𝑎𝑛∅ −1) 𝑡𝑎𝑛∅(𝑒 3𝜋𝑡𝑎𝑛∅ +1)
ψ= =39.7, ϵ= = 50
2𝑡𝑎𝑛∅ 9𝑡𝑎𝑛 2 ∅+1
(𝑒 3𝜇𝑡𝑎𝑛∅ +1) 𝑒 4𝜋𝑡𝑎𝑛∅ −1
=56.6, 𝜇𝑐 = = 256
3(√9𝑡𝑎𝑛 2 ∅+1) 16𝑡𝑎𝑛∅
−4𝐵𝑡𝑎𝑛 2 ∅(𝑒 3𝜇𝑡𝑎𝑛∅ +1) 4𝐵𝑡𝑎𝑛 2 ∅(𝑒 3𝜇𝑡𝑎𝑛∅ +1)
𝑥̅ = (9𝑡𝑎𝑛2 ∅+1)(𝑒 2𝜋𝑡𝑎𝑛∅ −1) =-2.52 B, 𝑧̅ = 3(√9𝑡𝑎𝑛2 ∅+1)(𝑒 2𝜋𝑡𝑎𝑛∅ −1)=2.85 B
𝑧̅
And, 𝑠𝑖𝑛𝜂=√𝑥̅ 2 =0.75
+𝑧̅ 2
ii) Determination of K, A and E
𝑔𝛽𝑠𝑖𝑛𝜂 𝑔
𝐾=√ = 807, 𝐴 = 𝛾𝐵4 𝜇 = 0.0000577, and 𝐸 = 𝜓𝑐𝐵2 + 𝜖𝛾𝐵3 = 20700 kN
𝜇𝐶 𝐵 𝐶
iii) Determination of Md(max) in terms of λ
Md(max)=HQd(max)=HQλ
= 4x8290λ=33160λ
iv) Determination of λcr which corresponds to α=0
115
Lecture Notes on DSF 2021
𝐴 1 𝐴 𝑀𝑑(𝑚𝑎𝑥) 𝐴 𝑀𝑑(𝑚𝑎𝑥) 1
𝛼= (𝐸 − 𝑄𝐵) cos(𝐾𝑡) − 3 sin(𝐾𝑡) + 2 ( + 𝑄𝐵 − 𝐸)
𝐾2 2 𝐾 𝑡𝑑 𝐾 𝑡𝑑 2
For t=td equals to 0.4s
0.0000577 1 0.0000577 33160𝜆
𝛼= 2
(20700 − 8290 × 2.5) cos(0.8070 × 0.4) − sin(0.8070 × 0.4)
0.807 2 0.8073 0.4
0.0000577 33160𝜆 1
+ ( + 8290 × 2.5 − 20700)
0.8072 0.4 2
=0.9159 cos(0.3228) 9.10𝜆 sin(0.3228) + 2.94𝜆0.9181 − 1.834
=0.05λ−0.0474
For α=0, λ=0.948=𝜆𝑐𝑟
v) Determination of Md(max) for λ=𝜆𝑐𝑟
Md(max)=33160λcr= 33160 0.948=31436 kNm
vi) Determination of G1 and G2
𝐴 1 𝐴 𝑀𝑑(𝑚𝑎𝑥) 𝐴𝑀𝑑(𝑚𝑎𝑥)
𝐺1 = 𝐾2 (𝐸 − 2 𝑄𝐵) cos(𝐾𝑡𝑑 ) − 𝐾3 sin(𝐾𝑡𝑑 ) + 𝐾2 =0.9159
𝑡𝑑
𝐴 1 𝐴 𝑀𝑑(𝑚𝑎𝑥) 𝐴𝑀𝑑(𝑚𝑎𝑥)
𝐺2 = 𝐾 (𝐸 − 2 𝑄𝐵) sin(𝐾𝑡𝑑 ) − cos(𝐾𝑡𝑑 ) + 𝐾2𝑡 = −4.05
𝐾2 𝑡𝑑 𝑑
iv) The dynamic load applied to the footing is initially peak triangular force pulse (Fig.5.28).
v)The footing is assumed to be weightless and to impart uniform load to the soil surface
Analysis:
The applied load is assumed to be an initial-peak triangular force which decays to zero at time td
as shown in Fig.5.28. The peak load q is expressed in pressure units. Since the function is
discontinuous at time td, two equations are necessary
𝑡
For 0 ≤ 𝑡 ≤ 𝑡𝑑 Loading function =𝑞𝐵 (1 − 𝑡 )--------- Eq. (5.130)
𝑑
The bearing capacity factors depend on∅ and K, K being 2 (Distance OA)/H, Fig. 5.27. The
value of K locates the centre of the spiral which is the centre of rotation. Obviously the correct
value of K is that which yields the minimum value of the bearing capacity. It is obtained by trial
and error for each set of problem parameters. The values of 𝑁𝑐 , 𝑁𝑞 , 𝑁𝛾 for various values of ∅ and
K are given in Column 3, 4 and 5 of Table 5.2.
Any acceleration of the soil mass ACDBA due to the downward movement of the footing will
cause inertial forces which will resist such movement. The inertial forces are directly
proportional to the acceleration of each individual soil mass and thereby dependent on
displacements. The effective total inertial force is obtained by combining the inertial forces on
each separate mass using energy considerations.
𝑑2 ∆
The inertial force is given by 𝐼𝑓 = 𝑁𝐼 𝛾𝐵 𝑑𝑡 2 ------------ Eq. (5.135)
Where, ∆= Displacement at any time t
𝑁𝐼 =Coefficient of dynamic inertial shear resistance.
The coefficient 𝑁𝐼 depends on ∅ and K, and its values are listed in column no. 6 of Table 5.2
Displacement of the soil mass within the failure surface due to downward movement of the
footing will increase the restoring moment about the point 0, and the increase in moment will be
proportional to the displacement provided the rotation is not excessive. It is expressed as
𝑅𝐹 = 𝑁𝑅 𝐵𝛾∆--------- Eq. (5.136)
The coefficient NR also depends on ∅and K. Its values are listed in column no. 7 of Table 5.2.
The differential equations are established by equating the four vertical forces to zero. There must
be separate equations for before and after time t d
For 0 ≤ 𝑡 ≤ 𝑡𝑑
𝑑2 ∆ 𝑡
𝑁𝐼 𝛾𝐵 𝑑𝑡 2 + 𝑁𝑅 𝛾𝐵∆ + 𝑞𝑢 𝐵 − 𝑞𝐵 (1 − 𝑡 ) = 0------- Eq. (5.137)
𝑑
Or
For 𝑡 ≥ 𝑡𝑑
𝑑2 ∆
𝑁𝐼 𝛾𝐵 𝑑𝑡 2 + 𝑁𝑅 𝛾𝐵∆ + 𝑞𝑢 𝐵 = 0----------- Eq. (5.138)
𝑑2 ∆ 𝑁 𝑞
Or, 𝑑𝑡 2 + 𝑁 𝑅𝐵 ∆= − 𝑁 𝛾𝐵
𝑢
------------------ Eq. (5.139)
𝐼 𝐼
The solution of the differential equations will yield equations of footing displacement versus
time. The forms of the particular solutions of Eq. 5.140 are found to be
𝑞−𝑞𝑢 𝑞
∆= 𝐶1 cos(𝐾 , 𝑡) + 𝐶2 sin(𝐾 , 𝑡) + ( 𝑁 ) − (𝑁 ) 𝑡---------- Eq. (5.140)
𝑅 𝛾 𝑅 𝛾𝑡𝑑
𝑞
And ∆= 𝐶3 cos(𝐾 , 𝑡) + 𝐶4 sin(𝐾 , 𝑡) − (𝑁 𝑢𝛾)--------- Eq. (5.141)
𝑅
2𝑁𝑅
In which 𝐾 ′ = √ and 𝐶1 , 𝐶2 , 𝐶3 𝑎𝑛𝑑𝐶4 are coefficient of integration. The C1, C2 are evaluated
𝑁𝐼 𝑅
by the initial conditions. The coefficients C3 and C4 are obtained by conditions of displacement
118
Lecture Notes on DSF 2021
and velocity at t d as defined by Eq. (5.140). After finding the solution and substitution of the
coefficients yield non-dimensional Eqs. 5.140 & 141 as
For 0 ≤ 𝑡 ≤ 𝑡𝑑
𝑞
𝑁𝑅 𝛾 𝑞 ′ 𝑞𝑢
( 𝑞 ) ∆= (𝑞 − 1) [1 − cos(𝐾 𝑡)] + [sin(𝐾 ′ 𝑡) − 𝐾 ′ 𝑡]---------- Eq. (5.142)
𝑢 𝑢 𝑡𝑑 𝐾′
For 𝑡 ≥ 𝑡𝑑
𝑞 𝑞
𝑁𝑅 𝛾 𝑞 𝑞𝑢 ′ ′ 𝑞𝑢
( ) ∆= [(1 − )+ 𝑠𝑖𝑛(𝐾 𝑡𝑑 )] cos(𝐾 𝑡𝑑 ) + [ (1 − cos(𝐾 ′ 𝑡𝑑 ))] sin(𝐾 ′ 𝑡𝑑 ) − − − − −Eq. (5.143)
𝑞𝑢 𝑞𝑢 𝑡𝑑 𝐾 ′ 𝑡𝑑 𝐾 ′
119
Lecture Notes on DSF 2021
Fig. 5.29: Non-dimensional maximum displacement
Table 5.2: Bearing Capacity Factors
Φ(dgree) K 𝑁𝛾 𝑁𝑐 𝑁𝑞 𝑁𝐼 𝑁𝑅 𝑁𝑅
√
𝑁𝐼
0 -0.05 0.000 5.7277 1.00 0.0633 2.0125 5.6366
0.00 0.000 5.7124 1.00 0.0631 1.9723 5.5887
+0.05 0.000 5.7258 1.00 0.0633 1.9433 5.5394
-0.65 0.0454 79.6255 7,9664 0.03755 8.9076 4.8709
-0.60 0.1445 29.8163 3.6986 0.2280 6.4362 5.3126
-0.55 0.1481 18.9958 2.6619 0.1579 5.0332 5.6460
-0.50 0.1553 14.3469 2.2552 0.1213 4.1699 5.8636
-0.45 0.1655 11.8179 2.0339 0.1011 3.6088 5.9750
-0.40 0.1786 10.2699 1.8985 0.0897 3.2299 6.0020
-0.35 0.1945 9.2580 1.8100 0.0833 2.9674 5.9698
5 -0.30 0.2131 8.5723 1.7500 0.0799 2.7828 5.9005
-0.25 0.2344 8.1007 1.7087 0.0786 2.6523 5.8108
-0.20 0.2585 7.7778 1.6805 0.0785 2.5604 5.7116
-0.15 0.2855 7.5629 1.6617 0.0793 2.4969 5.6099
-0.10 0.3154 7.4291 1.6500 0.0809 2.4547 5.5096
-0.05 0.3483 7.3580 1.6437 0.0829 0.4288 5.4128
0.00 0.3843 7.3366 1.6419 0.0853 2.4155 5.3205
0.05 0.4233 7.3553 1.6435 0.0881 2.4122 5.2330
-0.60 0.5700 53.9491 10.5127 0.1120 5.7922 7.1922
-0.55 0.5588 28.9945 6.1125 0.0935 4.8411 7.1948
-0.50 0.5645 20.5266 4.6194 0.0833 4.2238 7.1228
-0.45 0.5832 16.3539 3.8837 0.0779 3.8095 6.9932
-0.40 0.6127 13.9337 3.4569 0.0757 3.5264 6.8293
-0.35 0.6521 12.4031 3.1870 0.0755 3.3323 6.6445
-0.30 0.7008 11.3881 3.0080 0.0767 3.2008 6.4587
10 -0.25 0.7586 10.7004 2.8868 0.0790 3.1147 6.2781
-0.20 0.8253 10.2345 2.8046 0.0821 3.0625 6.1071
-0.15 0.9012 9.9267 2.7503 0.0858 3.0360 5.9474
-0.10 0.9863 9.7361 2.7167 0.0901 3.0294 5.7994
-0.05 1.0807 9.6352 2.6990 0.0948 3.0386 5.6676
0.00 1.1848 9.6049 2.6936 0.0999 3.0604 5.5360
+0.05 1.2986 9.6313 2.6983 0.1053 3.0923 3.4187
-0.55 1.5462 46.5473 13.4724 0.0707 5.2677 8.6324
-0.50 1.5198 30.2759 9.1124 0.0696 4.7177 8.2310
-0.45 1.5342 23.2038 7.2175 0.0707 4.3564 7.8481
15 -0.40 1.5806 19.3483 6.1844 0.0734 4.1189 7.4903
-0.35 1.6540 16.9964 5.5542 0.0773 3.9669 7.1622
-0.30 1.7520 15.4722 5.1458 0.0823 3.8766 6.8645
-0.25 1.8730 14.4550 4.8732 0.0881 3.8322 6.5961
-0.20 2.0166 13.7730 4.6905 0.0947 3.8232 6.3542
120
Lecture Notes on DSF 2021
-0.15 2.1825 13.3257 4.5706 0.1020 3.8418 6.1361
Φ(dgree) K 𝑁𝛾 𝑁𝑐 𝑁𝑞 𝑁𝐼 𝑁𝑅 𝑁𝑅
√
𝑁𝐼
15 -0.10 2.3710 13.0501 4.4968 0.1101 3.8825 5.9388
-0.05 2.5823 12.9048 4.4579 0.1183 3.9413 5.7596
0.00 2.8168 12.8613 4.4462 0.1282 4.0149 5.5961
+0.05 3.0750 12.8991 4.4563 0.1383 4.1008 5.4463
-0.50 3.6745 46.2884 17.8477 0.0673 5.6658 9.1768
-0.45 3.6419 33.8986 13.3381 0.0728 5.3067 8.5380
-0.40 3.6943 27.6099 11.0492 0.0796 5.0886 7.9941
-0.35 3.8151 23.9213 9.7067 0.0877 4.9684 7.5267
-0.30 3.9952 21.5875 8.8572 0.0970 4.9199 7.1214
20
-0.25 4.2298 20.0542 8.2992 0.1076 4.9258 6.7672
-0.20 4.5161 19.0369 7.9289 0.1194 4.9764 6.4552
-0.15 5.8533 18.3742 7.6877 0.1325 5.0582 6.1783
-0.10 5.2413 17.9678 7.5398 0.1470 5.1704 5.9309
-0.05 5.0864 17.7542 7.4620 0.1629 5.3068 5.7084
0.00 6.1717 17.6903 7.4368 0.1802 5.4638 5.5072
+0.05 6.7161 17.4757 7.4589 0.1989 5.6486 5.3243
-0.50 8.5665 73.8778 35.4499 0.0732 7.2346 9.9384
-0.45 8.3599 51.2706 24.9079 0.0835 6.8363 9.0503
-0.40 8.3728 40.7056 19.9814 0.0954 6.6214 8.3291
-0.35 8.5541 34.7663 17.2119 0.1094 6.5339 7.7297
-0.30 8.8760 31.1015 15.5029 0.1254 6.5404 7.2223
-0.25 9.3230 28.7315 14.3977 0.1437 6.6199 6.7864
25
-0.20 9.8871 27.1750 13.6720 0.1646 6.7584 6.4075
-0.15 10.5646 26.1681 13.2024 0.1882 6.9462 6.0748
-0.10 11.3542 25.5533 12.9157 0.2148 7.1761 5.7803
-0.05 12.2569 25.2309 12.7654 0.2445 7.4429 5.5178
0.00 13.2745 25.1345 12.7205 0.2775 7.7423 5.2825
+0.05 14.4095 25.2180 12.7594 0.3139 8.0710 5.0704
-0.45 19.3095 80.8644 47.6872 0.1064 9.3123 9.3540
-0.40 19.1315 62.4470 37.0539 0.1267 9.0899 8.4705
-0.35 19.3718 52.5548 31.3426 0.1506 9.0494 7.7518
-0.30 19.9400 46.6067 27.9084 0.1787 9.1446 7.1533
-0.25 20.1887 42.8208 25.7226 0.2116 9.3473 6.6458
30
-0.20 21.9566 40.3597 24.3017 0.2500 9.6392 6.2095
-0.15 23.3512 38.7778 23.3884 0.2944 10.0081 5.8303
-0.10 24.9984 37.8159 22.8330 0.3456 10.4452 5.4979
-0.05 26.8993 37.3127 22.5425 0.4041 10.9441 5.2044
0.00 29.0580 37.1624 22.4558 0.4706 11.4998 4.9436
+0.05 31.4810 37.2926 22.5309 0.5457 12.1084 4.7107
121
Lecture Notes on DSF 2021
35 -0.45 46.2942 134.3023 95.0397 0.1527 13.4981 9.4021
Φ(dgree) K 𝑁𝛾 𝑁𝑐 𝑁𝑞 𝑁𝐼 𝑁𝑅 𝑁𝑅
√
𝑁𝐼
-0.40 45.4427 100.66099 71.4837 0.1887 13.2639 8.3844
-0.35 45.6687 83.4477 59.4308 0.2323 13.3114 7.5703
-0.30 46.7356 73.3676 52.3727 0.2849 13.5708 6.9017
-0.25 48.5145 67.0529 47.9511 0.3481 14.0015 6.3419
-0.20 50.9356 62.9887 45.1052 0.4237 14.5786 5.8661
35
-0.15 53.9640 60.3926 43.2874 0.5133 15.2895 5.4569
-0.10 57.8568 58.8199 42.1864 0.6191 16.1127 5.1018
-0.05 61.8051 57.9989 41.6113 0.7428 17.0515 4.7911
0.00 66.6196 57.7539 41.4398 0.8868 18.0970 4.5175
+0.05 72.0773 57.9662 41.5884 1.0526 19.2451 4.2753
-0.40 115.7097 172.8231 146.0161 0.3229 20.8738 8.0404
-0.35 115.5504 141.1002 119.3973 0.4107 21.1138 7.1701
-0.30 117.6386 123.0124 104.2199 0.5195 21.7125 6.4650
-0.25 121.5875 111.8576 94.8599 0.6536 22.6077 5.8817
-0.20 127.1879 104.7472 88.8935 0.8175 23.7619 5.3914
40
-0.15 134.3346 100.2323 85.1051 1.0168 25.1570 4.9741
-0.10 142.9868 97.5069 82.8181 1.2572 26.7775 4.6152
-0.05 153.1451 96.0866 81.6263 1.5450 28.6173 4.3038
0.00 164.839 95.6630 81.2709 1.8870 30.6724 4.0317
+0.05 178.1176 96.0303 81.5791 2.2904 32.9409 3.7924
-0.40 327.6781 322.2748 323.2752 0.6576 36.2961 7.4295
-0.35 325.4943 259.1345 260.1349 0.8611 37.0113 6.5559
-0.30 329.9752 224.0769 225.0772 1.1194 38.3965 5.8568
-0.25 339.8627 202.7837 203.7840 1.4447 40.3468 5.2846
45 -0.20 354.4804 189.3358 190.3361 1.8515 42.8070 4.8083
-0.15 373.4971 180.8450 181.8452 2.3565 45.7496 4.4062
-0.10 393.7473 175.7358 176.7361 2.9784 49.1634 4.0628
-0.05 424.2605 173.0775 174.0778 3.7386 53.0475 3.7669
0.00 456.1177 172.2851 173.2853 4.6607 57.4067 3.5096
+0.05 492.4763 172.9729 173.9732 5.7709 62.2499 3.2843
122
Lecture Notes on DSF 2021
The shallow foundations may fail during seismic events. Published studies relating to the
bearing capacity of shallow foundations in such instances are rare. In 1993, however,
Richards et al. developed a seismic bearing capacity theory to find seismic bearing capacity
of granular soil. Figure 5.30 shows a failure surface in soil assumed for the subsequent
analysis, under static conditions. Similarly, Figure 5.31 shows the assumed failure under
earthquake conditions. Note that, in the two figures,
𝛼𝐴 ,𝛼𝐴𝐸 = inclination angles for active pressure conditions
And
𝛼𝑃 ,𝛼𝑃𝐸 = inclination angles for passive pressure conditions
Fig.5.30: Assumed failure surface in soil for static bearing capacity analysis
Fig.5.31: Assumed failure surface in soil for seismic bearing capacity analysis
According to this theory, the ultimate bearing capacities for continuous foundations in
granular soil are
1
𝑞𝑢 = 𝑞𝑁𝑞 + 2 𝛾𝐵𝑁𝛾 , For static condition--------------------- Eq. (5.144)
1
𝑞𝑢𝐸 = 𝑞𝑁𝑞𝐸 + 2 𝛾𝐵𝑁𝛾𝐸 , For Earthquake conditions------------ Eq. (5.145)
Where, 𝑁𝑞 , 𝑁𝛾 , 𝑁𝑞𝐸 , 𝑁𝛾𝐸 are bearing capacity factors
123
Lecture Notes on DSF 2021
Again, 𝑁𝑞 , 𝑁𝛾 = 𝑓(∅′ )
And 𝑁𝑞𝐸 , 𝑁𝛾𝐸 = 𝑓(∅′ , 𝑡𝑎𝑛𝜃)
𝑘ℎ
Where 𝑡𝑎𝑛𝜃 = ---------------- Eq. (5.146)
1−𝑘𝑣
𝑘ℎ = Horizontal coefficient of earthquake acceleration
𝑘𝑣 = Vertical coefficient of earthquake acceleration
Using the failure surface shown in Figure 5.31, Richards, Elms and Budhu (1993) provided
the values of bearing capacity factors, 𝑁𝑞 , and 𝑁𝛾 ,. They are given in Table 5.3
Table 5.3: Bearing capacity factors
The variations of 𝑁𝑞 and 𝑁𝛾 with ∅′ are shown in Figure 5.32. Figure 5.33 shows the variations
𝑁 𝑁
of 𝛾𝐸⁄𝑁 and 𝑞𝐸⁄𝑁 with 𝑡𝑎𝑛𝜃 and the soil angle ∅′ based on this analysis.
𝛾 𝑞
Fig. 5.32: Variation of 𝑁𝑞 and 𝑁𝛾 based on failure surface assumed in Figure 5.30
124
Lecture Notes on DSF 2021
𝑁 𝑁
Fig. 5.33: Variation of 𝛾𝐸⁄𝑁 and 𝑞𝐸⁄𝑁 with 𝑡𝑎𝑛𝜃
𝛾 𝑞
Under static conditions, bearing capacity failure can lead to a substantial sudden downward
movement of the foundation. However, bearing capacity related settlement in an earthquake is
𝑘
important and it takes place when the ratio 𝑡𝑎𝑛𝜃 = ℎ⁄(1 − 𝑘 )reaches the critical value. The
𝑣
𝑘
critical value can be expressed as [(1−𝑘ℎ )] 𝑐𝑟𝑖𝑡𝑖𝑐𝑎𝑙 becomes equal to 𝑘ℎ∗ when 𝑘𝑣 =0
𝑣
Table 5.4: Variation of 𝛼𝐴𝐸 with 𝑘ℎ∗ and soil friction angle ∅′
Figure 5.34 shows the variation of 𝑘ℎ∗ (for 𝑘𝑣 = 0 ) with the factor of safety (FS) applied to the
𝐷
ultimate static bearing capacity [Eq. 5.144], with ∅′ , and with 𝑓⁄𝐵 (for ∅′ = 300 and 400 ).
125
Lecture Notes on DSF 2021
The settlement of a strip foundation due to an earthquake using a sliding block approach can be
estimated (Richards, Elms and Budhu, 1993) as
−4
𝑉2 𝑘∗
𝑆𝐸𝑞 (𝑚) = 0.174 𝐴𝑔 [ 𝐴ℎ ] 𝑡𝑎𝑛𝛼𝐴𝐸 ---------------------- Eq. (5.147)
Where, V = peak velocity for the design earthquake (m/sec)
A = acceleration coefficient for the design earthquake
g = acceleration due to gravity (9.81 m/s2)
The values of 𝑘ℎ∗ and 𝛼𝐴𝐸 can be obtained from Figure 5.34 and Table 5.4, respectively. This
approach can be used to design a footing based on limiting seismic settlements.
Problem No.1
A rectangular foundation has a length L of 2.5 m. It is supported by medium dense sand with a
unit weight of 17 kN/m3. The sand has an angle of friction of 36º. The foundation may be
subjected to a dynamic load of 735 kN increasing at a moderated rate. Using a factor of safety
equal to 2, determine the width of the foundation. Use Df = 0.8 m.
126
Lecture Notes on DSF 2021
DESIGN OF RECIPROCATING MACHINE FOUNDATION
Reciprocating Machine should be placed over suitable Vibration Absorber to reduce the
magnitude of displacement due to dynamic unbalanced force produced due the Machine.
Q.1 Determine the stiffness of the absorber to be kept between a reciprocating machine and
foundation to bring the vibration amplitude to less than 0.02 mm. The weight of the machine is
18 kN. It produces an unbalanced force of 4 kN, when operated at speed of 600rpm. Shear
modulus of foundation soil G=20MN/m2 and µ= 0.35.
Sol: Mass of the block foundation of dimension 4x3m2 with 1.5 m height
24×4×3×1.5
𝑚1 = = 44kg
9.81
127
Lecture Notes on DSF 2021
18
Mass of machine 𝑚2 = 9.81 = 1.83kg,
4×3
Equivalent radius of circular footing 𝑟 = √ =1.95 m
𝜋
4𝐺𝑟
Stiffness in vertical direction, 𝐾𝑧 = 1−µ=240000 kN/m
240000
Natural frequency of vibration, 𝜔𝑛𝑙1 = √44+1.83=72.36 r/sec
2𝜋×600
Operating frequency of machine, 𝜔 = = 62.83 r/sec
60
𝑍1 0.02
Efficiency of absorber= 𝑍 =0.068=0.333
𝑚𝑎𝑥
𝜔𝑛𝑙1 𝑚
Frequency ratio, 𝑎1 = =1.15, Mass ratio, 𝜂 = 𝑚1 =0.0415
𝑤 2
2 𝑎 2 (1+0.0415)(1.152 −1)
Hence, 0.333=[1−(1+0.0415)(1.152 +𝑎 2−1.152 𝑎 2 )]
2 2
So 𝑎2 =1.31
𝜔𝑛𝑙2
But 𝑎2 = , hence 𝜔𝑛𝑙2 =81.97 r/s
𝜔
[(1+𝑛)𝑎12+𝑛𝑎22 −1]𝐹0 [(1+0.0415)1.152 +0.0415×1.312 −1]×4
Z2=𝑚 = = −58.47mm
2 𝜔2 [1−(1+𝑛)(𝑎 2+𝑎 2 −𝑎 2𝑎 2 )] 1.83×62.832 [1−(1+0.0415)(1.152 +1.312 −1.152 ×1.312 )]
1 2 1 2
Piles divided into two groups: End bearing and Friction Piles
Longitudinal Vibration of Short Elastic Bar:
128
Lecture Notes on DSF 2021
i) End Bearing Pile: Mass of Pile is negligible
𝜕 2𝑢 𝜕2𝑢
Wave equation: =𝑣𝑐2 𝜕𝑥 2 -------- Eq.(1)
𝜕𝑡 2
Where u is the displacement function and vc is shear wave velocity
Solution to above second order differential equation is expressed as
u(x,t)=U(x)(A1sinωnt+A2cosωnt)--- Eq.(2)
Where A1 and A2 are two constants, ωn is the natural frequency and U(x) is amplitude of
displacement along the length of the bar and is independent of time “t”
Eq. 1 can be rewritten by putting value from Eq.2 as
𝜕 2 𝑈(𝑥) 𝜌
+ 𝐸 𝜔𝑛2 𝑈 (𝑥 ) = 0----- Eq.(3)
𝜕𝑥 2
Now general solution to the above Eq.3 can be given as
𝜔 𝑥 𝜔
U(x) =𝐵1 sin( 𝑣𝑛 ) +𝐵2 cos( 𝑣𝑛𝑥 )----- Eq. (4)
𝑐 𝑐
Using end conditions a) Fixed
At x=0, U(x) =0 which gives B2=0
𝑑𝑈(𝑥) 𝐵 𝜔𝑛𝐿
At x=L, = 0, gives 𝑣1 𝜔𝑛 𝑐𝑜𝑠 =0,
𝑑𝑥 𝑐 𝑣𝑆
𝜔𝑛 𝐿 𝜔𝑛 𝐿 (2𝑛−1)𝜋
Hence, 𝑐𝑜𝑠 = 0, or =
𝑣𝑐 𝑣𝑐 2
1 𝑣𝑐 𝜋 𝑣𝑐
Which gives, 𝜔𝑛 = 2 (2𝑛 − 1)𝜋 𝐿 , for n=1 𝜔𝑛 = 2 𝐿
𝜔𝑛 1 𝐸𝑝
So natural frequency 𝑓𝑛 = = √𝜌 --------------- Eq.(5)
2𝜋 4𝐿 𝑝
129
Lecture Notes on DSF 2021
𝜔𝑛 𝐿
At x=L ,𝐴𝐸 = 𝑚𝜔𝑛 𝑣𝑐 tan( )---------- Eq.(13)
𝑣𝑐
𝜔𝑛 𝐿
Or, 𝐴𝜌𝑣𝑐2 = 𝑚𝜔𝑛 𝑣𝑐 tan( )-------------- Eq.(14)
𝑣𝑐
𝐸
as𝑣𝑐 = √𝜌
Now consider a non-dimensional parameter,
𝐴𝐿𝛾 𝜔𝑛 𝐿 𝜔𝑛 𝐿
= tan( )--------------- Eq. (15)
𝑊 𝑣𝑐 𝑣𝑐
which may be expressed as 𝛽 = 𝛼𝑡𝑎𝑛𝛼
𝐴𝐿𝛾 𝜔𝑛 𝐿
Where 𝛽 = and 𝛼 =
𝑊 𝑣𝑐
The above relation can be placed in tabular form as Table1
T1: Coefficients for natural frequency of piles
𝛽 0.1 0.3 0.5 0.7 1.0 2.0 4.0 10.0
𝛼 0.32 0.53 0.66 0.75 0.86 1.08 1.27 1.43
Assumptions:
1 Pile is vertical and circular in cross section, if not equivalent circular radius is taken
2 The pile is floating the soil foundation, (Not restricted)
3 The pile is perfectly connected to soil.
4 The soil above the pile tip behaves as an infinitesimal thin independent linearly
elastic layer.
5 The dynamic stiffness and damping of pile material can be described interms of a
complex stiffness matrix as proposed by Novak and EI-Sharnouby 1983.
Where fz1 and fz2 are two factors can be taken from Fig. 1 (a),(b), proposed by Novak et al.
130
Lecture Notes on DSF 2021
Fig. 1:Stiffness and Damping Parameters of Vertical Response for: a) End Bearing Piles and b)
Floating Piles (reprinted from Novak and El Sharnouby 1983, © ASCE)
Where 𝛼𝑟 = Interaction factor and is obtained from the table T:3.6 page 49 (Hand book of
Machine Foundation by Srinivasulu &Vaidyanatham)
Where r0 is the equivalent radius of pile cap and Df is the depth of pile cap from ground surface
S1 and S2 are constants and may be taken as 2.7 and 6.7 respectively.
131
Lecture Notes on DSF 2021
𝐾𝑧 (𝑇)
𝜔𝑛 = √ ------------------------------------------------ Eq.(29)
𝑚
Q.1 A group of four piles having dimension 0.3x0.3 m2 is supported a machine foundation as
shown in Fig. Determine total stiffness and damping coefficient, given E p=2.1x107 kN/m2, unit
weight of soil18.9 kN/m2, Poisson’s ratio 0.5 and shear modulus Gs=28120 kN/m2.
GL
1.5m 2.0 m
12m
B
A
2.1 m
D C
1.5 m
(1.5𝑋0.6)2
Solution: Equivalent radius of pile cap r0=√ = 1.18 m
𝜋
132
Lecture Notes on DSF 2021
(0.3𝑋0.3)2
Equivalent radius of pile R= √ = 0.17 m
𝜋
𝐴𝐸 2.1𝑥107 ×0.32
Vertical stiffness of pile, 𝐾𝑧 = ( 𝑅 )𝑓𝑧1 = × 𝑓𝑧1 , 𝑓𝑧1 is found from Fig 1(b) as 0.035
0.17
3 2
= 389.1x10 kN/m
E 2.1×107
Now = = 746.8= say 750
G 28120
𝐿 12
= 0.17=70.58, From Fig. 1(b)𝑓𝑧2 =0.06
𝑅
𝐸𝐴 2.1×107 ×0.32
Hence Damping Coefficient is calculated as 𝐶𝑧 = ( 𝐺)fz2 = 28120×9.81 x0.06=938.7 kN.s/m
√
𝜌
√ 18.9
Calculation of Interaction factor:
𝑆 1.5
= =4.41 for pile A to B, A to D
𝐷 0.34
A A B C D
1 0.54 0.48 0.54
B 0.54 1.0 0.54 0.48
C 0.48 0.54 1.0 0.54
D 0.54 0.48 0.54 10
∑ αr 2.56 2.56 2.56 2.56
If the sum of interaction factors is not same, than take average value of ∑ 𝛼𝑟
∑𝐾 4×389100
Kzr(g)=∑ 𝛼𝑧 = =607968.75kN/m
𝑟 2.56
∑𝐶
𝐶𝑍 (𝑔) = ∑ 𝑧 =1446.78 kNs/m
𝛼𝑟
Problem 1 A machine is supported by four pre-stressed concrete piles driven into a bed rock.
The length of each pile is 8o ft long and is 12x12 in2 in cross section. The weight of the machine
and foundation is 300x103 lbs, unit weight 150 lb/ft 3. The Young’s modulus is 3.5x106 lb/in2.
Determine the natural frequency of pile foundation system.
134
Lecture Notes on DSF 2021
135
Lecture Notes on DSF 2021
𝑊 𝐿3 2𝐾+1
𝑚
𝛿1 = 96𝐸𝐼 ,
𝑏 𝐾+2
𝑞𝐿4 5𝐾+2
𝛿2 =
384𝐸𝐼𝑏 𝐾+2
3 𝑞𝐿
𝛿3 = (𝑊𝑚 + )
5𝐸𝐴𝑏 2
ℎ 𝑊𝑚 𝐿
𝛿4 = (𝑊𝐿 + )
𝐸𝐴𝑐 2
Where
𝐴𝑏 =Cross sectional area of beam, 𝐴𝑐 = Cross sectional area of column
𝐼𝑐 =Moment of inertia of beam about the axis of bending
𝐼 ℎ
E= Young’s modulus of concrete, K=Relative stiffness factor= 𝑏 ×
𝐼𝑐 𝐿
𝑏
Fig. 2: values of a and b for a frame with haunches. Fig.3: α versus 𝑙
0
136
Lecture Notes on DSF 2021
The natural frequency of a transverse frame in vertical vibrations is given by
𝐾𝑧 𝑔
𝜔𝑛𝑧 = √ ------- Eq.3
𝑊
137
Lecture Notes on DSF 2021
∑(𝐾𝑥 )𝑔
𝜔𝑛𝑥𝑎 = √ --------------- Eq.8
𝑊𝑇
138
Lecture Notes on DSF 2021
𝑊4 =Weight of the tranverse beam
The stiffness K1 of both the column of a transverse frame is given by
2𝐸𝐴𝑐
𝐾1 = --------------- Eq.13
ℎ
𝐾
𝜔𝑛𝑙2 = √𝑚2 ------------ Eq.22
2
𝑚
𝜂 = 𝑚1 --------- Eq.23
2
Now the two natural frequencies of the system for forced vibration condition can be obtained by
considering the equations motion as
𝑚1 𝑍1̈ + 𝐾1 𝑍1 + 𝐾2 (𝑍1 − 𝑍2 ) = 0-------- Eq.24
139
Lecture Notes on DSF 2021
𝑚2 𝑍2̈ + 𝐾2 (𝑍2 − 𝑍1 ) = 𝐹𝑧 𝑠𝑖𝑛𝜔𝑡------- Eq.25
By solving these equations the amplitude of vertical vibration can be obtained as
2
𝐹𝑧 𝜔𝑛𝑙2
𝐴𝑍1 = 𝑚 4 2 2 2 2 2 ------- Eq.26
1 ⌊𝜔 −(1+𝜂)(𝜔𝑛𝑙1 +𝜔𝑛𝑙2 )𝜔 +(1+𝜂)𝜔𝑛𝑙1 𝜔𝑛𝑙2 ⌋
2 2
𝐹𝑧 ⌊(1+𝜂)𝜔𝑛𝑙 +𝜂𝜔𝑛𝑙2 −𝜔2 ⌋
𝐴𝑍2 = 𝑚 4 −(1+𝜂)(𝜔2 +𝜔2 )𝜔2 +(1+𝜂)𝜔2 𝜔2 ⌋ ------- Eq.27
2 ⌊𝜔 𝑛𝑙1 𝑛𝑙2 𝑛𝑙1 𝑛𝑙2
Problem No. 1
Plan of deck slab with loading position is shown. A reinforced concrete frame with vertical loads
at bvarious points are also shown. The details of these loads are
1 and 2 = 5t each, 3,4,5 and 6 = 2t each
Ec = 3x106 t/m2 and unit weight of concrete = 2.24t/m3.
Calculate the natural frequency of horizontal vibration in the longitudinal direction by treating
the frame vas single degree freedom system.
6 2 4
2.0m
4.8 2.8m
m
5 1 3
2.0m
1.2m 1.2m
1.0 m 3.1m
140
Lecture Notes on DSF 2021
0.6m
1.0 m 0.4 m
0.4 m
0.9m 4.8 m
5.4 m
m
1.0 m
1.0 m 0.4 m
0.4 m
0.6 m 3.6 m
1.0 m
6.8 m
Q.2 A spring and dashpot are attached to a body weighing 140 N. The spring constant is 3.0
kN/m. The dashpot has a resistance of 0.75 N at a velocity of 0.06 m/s. Determine the following
for free vibration:whether the system is over damped, under damped or critically damped.
Q.3 A counter rotating eccentric mass exciter is used to produce forced oscillation of a spring
supported mass. By varying the speed of rotation, resonant amplitude of 5 mm was recorded.
When the speed of rotation was increased considerably beyond the resonant frequency, the
amplitude appeared to approach a constant value of 0.6 mm. Determine the damping factor of the
system.
Q.5 A body weighing 600 N is suspended from a spring which deflects 12 mm under the load.
It is subjected to a damping effect adjusted to a value 0.2 times that required for critical damping.
Find the natural frequency of the un-damped and damped vibrations, and in the latter case,
determine the ratio of successive amplitudes.
Q.6 In a cyclic plate load test on a plate of 0.60 m × 0.60 m size settles 0.65 mm under a
pressure of 20 kN/m2. On unloading observed plate settlement was 0.60 mm. Determine the
value of coefficient of elastic uniform compression of the soil.
Q.7 A mass attached to a spring of 5 N/mm has a viscous damping device. When the mass
was displaced and released, the period of vibration was found to be 2 s and ratio of the
consecutive amplitudes was 10/3. Determine the damping factor and natural frequency of the
system. Determine also the amplitude of motion when a force of 3 sin 4t N acts on the system.
Q.8 A machine of weight 17.5 kN and operating frequency 400 rpm has to be installed on
ground which has properties G= 40 MN/m2 , 𝛾𝑠 =20 kN/m3 and μ=0.3. The machine contains a
unbalanced rotating parts which produce an eccentric moment of 18 Nm in vertical direction.
The permissible amplitude of vibration for the system is 0.2 mm and the equivalent diameter of
the foundation required to install the machine is 1.6 m. Design the foundation.
Q.9 Determine the stiffness of the absorber to be kept between a reciprocating machine and
foundation to bring the vibration amplitude to less than 0.02 mm. The weight of the machine is
142
Lecture Notes on DSF 2021
25 kN. It produces an unbalanced force of 4 Kn in the vertical direction, when operated at speed
of 750 rpm. Shear modulus of foundation soil G=2.5x104kN/m2 and Poisson’s ratio 0.3.
Q.10 A foundation is subjected to a constant force type vibration. Given that the total weight of
machine and foundation block is 1500 N. Unit of foundation soil γ= 15 kN/m3, Shear modulus
and poison’s ratio are 15 MN/m2 and 0.4 respectively. The amplitude of vibrating force F0 =
1500N.Operating frequency of machine is 80 cpm. The size of foundation block is 10m long and
3 m wide. Determine
Q.11 A machine is supported by four pre-stressed concrete piles driven into a bed rock. The
length of each pile is 80 ft and they are 12x12in2 in cross section. The weight of the machine and
foundation is 300x103 lbs, Unit of concrete is 150 lb/ft 3. Young’s modulus 3.5x106 lb/ft3.
Determine the natural frequency of pile foundation system.
Q.12 Following are the field standard penetration test number (N) in a deposit of sand. Ground
water table is encountered at a depth of 3m below the ground surface. Soil properties of sand are,
dry unit weight 18.5kN/m3 and saturated unit weight 20.6 kN/m3. Determine for an earth quake
magnitude of 7.5, whether liquefaction will occur? Assume ground acceleration as a max= 0.15g.
Q.13 A horizontal piston type compressor is placed on a block type foundation as shown In
Fig.1
The operating frequency is 600 cpm. The amplitude of the horizontal unbalanced force of
compression is 30kN and it produces a rocking motion of the foundation about point O. The
mass moment of inertia of the compressor assembly about the axis BOB’ is 16x10 5kg.m2.
Determine
i) The resonance frequency
ii) The amplitude of rocking vibration at resonance.
Q.14 A concrete bock foundation of a machine has the following dimension, L=4m, B=3m and
height H=1.5 m. The foundation is subjected to a sinusoidal horizontal force from the machine
having amplitude of 10 kN at a height of 2.0 m from the base of the foundation as shown in
Fig.2.
The soil supporting the foundation is sandy clay with G=30,000 kN/m2, μ=0.25and unit weight
γ=17 kN/m3. Determine
i) The resonance frequency for sliding and rocking mode of vibration of the
foundation (Independent mode analysis)
143
Lecture Notes on DSF 2021
ii) Total horizontal displacement at the top of the foundation block.
Fig.1
10 kN
2.0 m
1.5 m
3.0 m
+
1.5 m 4.0m
3.0 m
Fig. 2 Dimension of Block Foundation
Q.15 A machine is supported by four pre-stressed concrete piles driven into a bed rock. The
length of each pile is 8o ft long and is 12x12 in2 in cross section. The weight of the machine and
foundation is 300x103 lbs, unit weight 150 lb/ft 3. The Young’s modulus is 3.5x106 lb/in2.
Determine the natural frequency of pile foundation system.
144
Lecture Notes on DSF 2021
MULTIPLE CHOICE QUESTIONS
Q.1 A spring and dashpot are attached to a body weighing 200 N. The spring constant is 2.0 kN/m.
The dashpot has a resistance of 0.65 N at a velocity of 0.05 m/s. If the system is set to a motion, the
system is:
(a) Over damped (b) Under damped (c) Critically damped (d) None of above
Q.2 A counter rotating eccentric mass exciter is used to produce forced oscillation of a spring
supported mass. By varying the speed of rotation, resonant amplitude of 6 mm was recorded. When the
speed of rotation was increased considerably beyond the resonant frequency, the amplitude appeared to
approach a constant value of 0.6 mm. The damping factor (in %) of the system is:
Q.3 An SDF system is excited by a sinusoidal force. At resonance the amplitude of displacement was
measured to be 4 mm. At an exciting frequency of one-tenth of the natural frequency of the system, the
displacement amplitude was measured to be 0.3 mm. The damping factor of the system (in %) is:
Q.4 A body weighing 600 N is suspended from a spring which deflects 10 mm under the load. It is
subjected to a damping effect adjusted to a value 0.4 times that required for critical damping. The un-
damped natural frequency of the vibrations (in radian per second) is:
Q.5 The damped natural frequency of the vibrations (in radian per second) of the system in Q 4 is:
Q.6 The ratio of successive peak amplitudes of the vibrations for the system of Q 4 is:
Q.7 A mass attached to a spring of 9 N/mm has a viscous damping device. When the mass was
displaced and released, ratio of the consecutive amplitudes was 11/6. The damping factor of the system
(in %) is:
Q.8 As per Indian standard code IS 2974, IV, the permissible amplitude for Rotating machine speed
>1500 rpm is
Q.9 If ω and ωn are operating and natural frequency of system respectively, which one is correct for
the design condition of dynamically loaded foundation
145
Lecture Notes on DSF 2021
Q.10 An unknown weight W is attached to the end of an unknown spring k and natural frequency of
the system was found to be 92cpm. If 1 kg weight is added to W, the natural frequency reduced to 80
cpm. Find the unknown weight
Q.11 The value of equivalent radius of circular footing for rocking mode of vibration of a rectangular
foundation of size 4mx6m is
Q.14 Modified mass ratio as per Lysmer Analysis for a block type foundation of size 4x6 m2 subjected
to vertical mode of vibration due to total weight of 150 kg resting on foundation soil having unit weight
and poison’s ratio as 17.2 kN/m3 and 0.25 respectively is
Q.16 According to the Richart, the maximum operating frequency up to which no noticeable amplitude
of vibration can be identified by person
(a) 500 rpm (b) 750 rpm (c) 1000 rpm (d) 2000 rpm
Q.17 According the IS 5249, the relation between Coefficient of uniform elastic compression for
various mode of excitation can be given as
(a) CU=1.73 Cτ, Cφ=2CU (b)CU=2Cτ, Cφ=1.75CU (c)CU=2.25Cτ, Cφ=2.2CU (d)None of above
Q.18 The value of coefficient of elastic uniform compression C u of soil obtained from block vibration
test of contact area of 10 m2 is 12 kN/m2. The of Cu for a base area of foundation of 12 m2 is
(a)15 kN/m2 (b) 10.65 kN/m2 (c) 12 kN/m2 (d) None of above
Q.19 Coefficients for natural frequency of end bearing pile of length 12m resting on soil having shear
velocity of 350 m/s is 0.43. The natural frequency of vibration is
(a) 14.56 rpm (b) 18.98 rpm (c) 12.54 rpm (d) None of above
Q.20 The efficiency of absorber η used in foundation for reciprocating engine when maximum
displacement is 0.034 mm and the permissible displacement of 0.02 mm is