Unit - 7 Vibration of Continuous System
Unit - 7 Vibration of Continuous System
Continuous systems are tore which have continuously distributed mass and Elasticity. These continuous
systems are assumed to be homogeneous and isotropic obeying hook law within the elastic limit. Since to
specify the position of every point in the continuous systems an infinite number of coordinates is required
therefore a continuous systems is considered to have infinite number of degrees of freedom. Thus there will
be infinite natural frequencies
Free vibration of continuous system is sum of the principal or normal modes. For the normal mode
vibration, every particle of the body performs simple harmonic motion at the frequency corresponding to
the particular root of the frequency equation, each particles passing simultaneously through its respective
equilibrium position. If the elastic curve resulting due to vibration motion starts coinciding exactly with
one of the normal modes, only that normal mode will be produced. If the elastic curve of the system under
which the vibration is started, is identical to any one of the principal mode shapes, then the system will
vibrate only in that principal mode.
Vibration of String.
A string stretched between two support is shown in figure-1 is an infinite degree of freedom system.
T
Y θ + (δθ
(δθ/δx)
/δ dx
y
θ dx
x dx X
T
Figure-1
Let the tension T in the string be large so that for small displacement or amplitude the tension in the
string remains constant throughout the string. Consider an element of length dx at a distance ‘x’ from the
left end. At any instant of time the element of the string be displaced through a distance ‘y’ from the
equilibrium position. Then the tension at both the ends of this element is ‘T’. If θ the angle at the left end
of the element makes with the horizontal x- axis then the angle at the right end of the element is
θ + (δθ/δx)dx
The components of these two tensions at the ends of the element along x-axis balance each others.
The components along Y – axis
sin Since θ is small
1
If ‘ρ' is mass per unit length, then the mass of the element is ‘ρdx’ and the differential equation of motion
along Y-axis according to Newton’s second law of motion is,
Therefore
But we know slope then we have
or or
⁄
--- (2)
Where !⁄ ------ (3) is the velocity of wave propagation along the string. Solution of equation (2)
can be obtained by assuming Y through separation of variables to be a product of two functions in the
form. " , $ % & $ 4
, - - 0 -
& ./ ---- (6)
,
- ,
-
Left hand side of equation (6) is a function of x alone and the right hand side is a function of t alone. These
two can only be equal if each one of these expression is a constant. This constant may be positive, zero or
negative. If it is a positive constant or zero then there is no vibratory motion. Hence this constant has to be
2
negative constant and equal to -
0 -
Then ) 7
,
-
)
()" 1 ()&
2) 0 * 2) 0
% ( ) & ($ )
()" 4 ()&
) % 0 * 2) 0 8
( ) ($ )
4 4
" , $ A sin 0
: ;.< 0
D sin$ @;.<$ 11
Since 2 can have infinite values, being an infinite degree of freedom system the general solution can be
written as
" , $ ∑∞
EF BAisin 9 : ;.< 9 DDisin$ @;.<$ 11
8C 8C
The arbitrary constant A, B, D and E have to be determined from the boundary condition and the initial
conditions.
For the string of length ‘L’ stretched between two fixed points. The boundary conditions are
y(0,t) = y(L, t) = 0
4
" , $ A sin 0
D sin$ @;.<$ 12
The condition y (L, t) = 0 them fields the equation 0 BA sin H DD sin$ @;.<$
8
9
BQR/P
9
DH L or 2H/O -------- (16)
Each ‘n’ represents normal mode vibration with natural frequency determined from the equation
S
2HM/; therefore M ;/2H = )T √/V for n = 1,2,3…….. ------ (17)al with the
distribution % B sin L
H
D ----- (18)
∞
S\
" , $ WI sin t @ cos t sin T
19
1
S\
With 20
T
With the initial condition of Y(x,0) and Y (x, 0) the constants Dn and En can be evaluated.
The bar/rod considered for the analysis is assumed to be thin and a uniform cross sectional area throughout
its length. Due to axial force the displacement ‘u’ along the bar will be a function of both position ‘x’ and
time ‘t’. The bar/rod has infinite number of natural frequencies and modes of vibration
Consider an element of length dx at a distance ‘x’ from the left end of the bar as shown in the
figure -2.
At any instant of time during vibration let ‘P’ be the axial force at the left end of the element. Then the
axial force at the right end of the element is P+(
P+(∂p/∂x)dx
x)dx and If ‘u’ is the displacement at x, then the
displacement at x + dx will be u+(∂u/
u/∂x)dx.
x)dx. Then the element dx in the new position has changed in length
x)dx and hence the unit strain is ∂u/∂x.
by an amount of (∂u/∂x)dx
But from Hook’s law the ratio of unit stress to the unit strain is equal to the Modulus of Elasticity ‘E’,
‘E
Hence we have ∂u/∂xx = P/AE ------- (1)
In which ‘A’ is the cross sectional area of the bar/rod. Differentiating the equation (1) with
w respect to x we
have AE ∂2u/∂x2 = ∂P/∂x ------ (2) But (∂P/∂x) dx is the unbalanced force
Considering the Dynamic Equilibrium of the element from Newton’s second law of motions me have.
ass per unit volume, substituting equation (2) into (3) for ∂P/∂x we have
Where ρ is the density of rod/bar mass
Where c = √(E/r) ------ (6) is the velocity of propagation of the displacement or stress wave in the rod/bar.
Solution of equation (5) can be obtained by assuming u through separation of variables to be a product of
two functions in the form u(x, t) = U(x) G(t) ------ (7)
Left hand side of equation (9) is a function of x alone and the right hand side is a function of t alone. These
two can only be equal if each one of these expression is a constant.
This constant may be positive, zero or negative. If it is a positive constant or zero then there is no vibratory
motion. Hence this constant have to be negative constant and equal to - w2 . Then we have
With the general solution U = Asin(ω/c)x + Bcos(ω/c)x ----- (12) G = Dsinωt + Ecosωt ------ (13)
Since ω can have infinite values, as the system being an infinite degree of freedom the general solution can
be written as
The arbitrary constant A, B, D and E have to be determined from the boundary condition and the initial
conditions.
Example:
Determine the natural frequencies and mode shape of a bar/rod when both the ends are free.
free
Solution: When the bar/rod with both the ends being free then the stress and strain at the ends are zero
(∂u/∂x)x=0=(Aω/c)[Dsinwt + Ecosωt]
t] = 0 --- (4)
Since the equation (4) and (5) must be true for any time ‘t’
From equation (4) A must be equal to zero. Since B must be finite value in order to have vibration then
sin(ωnL/c) = 0 or ωnL/c =sin-1(0) = nπ
equation (5) is satisfied only when sin(
Each ‘n’ represents normal mode vibration with natural frequency ‘fn’ determined from the
equation (7).
Torsional
sional vibration of circular shaft.
Consider an element of length dx at a distance x from one end of the shaft as shown in the figure -3
7
Dept. Of ME, ACE 7
Mechanical Vibrations[10ME72]
At any given instant of time during vibration let T be the torque at left end of the element. Then the torque
at the right end of the element is given by T+(∂T/∂x)dx.
If θ is the angular twist of the shaft at the distance x, then θ+(∂θ/∂x)dx is the angular twist of the
shaft at the distance x+dx. Therefore the angular twist in the element of length dx is (∂θ/∂x)dx.
If ‘J’ is the polar moment of inertia of the shaft and ‘G’ the modulus of rigidity then the angular twist in the
element of length dx is given by torsion formula
In which GJ is the torsional stiffness of the shaft. Since the torque on the two face of the element being T
and T+(∂T/∂x)dx the net torque on the element will be (∂T/∂x)dx --- (4)
Substituting equation (3) for ∂T/∂x in to equation (4) we have (∂T/∂x)dx = GJ(∂2θ/∂x2)dx ------- (5)
Considering the dynamic equilibrium of the element one can obtain the equation of motion by equating the
product of mass moment of inertia ‘Jρdx’ and the acceleration ‘∂2θ/∂t2’ to the net torque acting on the
element we have.
∂2θ/∂t2 = (G/ρ)(∂2θ/∂x2) or
Where c = √(G/ρ) ------ (7) is the velocity of wave propagation in which r is the density of the shaft in
mass per unit volume.
Solution of equation (6) can be obtained by assuming q through separation of variables to be a product of
two functions in the form. θ(x, t) = Θ(x) G(t) ------ (7)
Left hand side of equation (9) is a function of ‘Θ’ alone and the right hand side is a function of ‘t’ alone.
These two can only be equal if each one of these expression is a constant. This constant may be positive,
zero or negative. If it is a positive constant or zero then there is no vibratory motion. Hence this constant
has to be negative constant and equal to - ω2
With the general solution Θ = Asin(ω/c)x + Bcos(ω/c)x ----- (12) G = Dsinωt + Ecosωt ------ (13)
Since ω can have infinite values, as the system being an infinite degree of freedom the general solution can
be written as
The arbitrary constant A, B, D and E have to be determined from the boundary condition and the initial
conditions.
To derive the differential equation of motion for lateral vibration of beams one has to consider the forces
and moments acting on the beam. Let V and M are the shear forces and bending moments respectively
acting on the beam with p(x) represents the intensity of lateral loading per unit length of the beam.
Consider a beam of length ‘L’ and moment of inerter ‘I’ subjected to distributed lateral load of p(x) N per
unit length. At any section x-x at a distance x from one end consider an element of length dx as shown in
the figure-4. The force and the moments acting on the element of length dx are as shown in the figure-5
For static equilibrium the summation of the forces and moment should be equal to zero. Thus summing the
vertical force in the Y – direction and equating to zero we have
Which state that the rate of change of shear force along the length of the beam is equal to the load per unit
length.
Summation of the moment about any point should be equal to zero. Thus summing the moments about the
right bottom corner of the element considering the clockwise moment as positive we have
Which state that the rate of change of bending moment along the length of the beam is equal to the shear
force.
From equation (1) and (2) we have ∂2M/∂x2 = ∂V/∂x = p(x) ----- (3)
But we known that the bending moment is related to curvature by flexure formula and for the coordinates
indicated we have EI(d2y/dx2) = M ----- (4)
Substituting equation (4) into equation (3) we have ∂2[EI(d2y/dx2) ]/∂x2 = p(x)
For the dynamic equilibrium when the beam is having transverse vibrations about its static equilibrium
position under its own weight and due to the distributed load per unit length must be equal to the inertia
load due its mass and acceleration. By assuming harmonic motion the inertia force ρωy2 is in the same
direction as that of the distributed load p(x).
In which ρ is the mass density per unit length of the beam. Thus substituting equation (6) into equation (5)
for p(x) we have EI(d4y/dx4) = ρω2y or EI(d4y/dx4) - ρω2y = 0 ------ (7)
If the flexural rigidity EI is constant then equation (7) can be written as d4y/dx4 - (ρω2/EI)y = 0 ---- (8)
The above equation (10) is a fourth order differential equation for the lateral vibration of a uniform cross
section beam. The solution of the above equation (10) can be obtained by assuming the displacement ‘y’ of
the form y = eax ---- (11)
Which will satisfy the fourth order differential equation (10) when a = ± β and a = ± iβ ---- (12)
Since e±βx = coshβx ± sinhβx ---- (13) and e±ibx = cosβx ± sinβx ---- (14)
Then the solution of the above equation (10) will be in the form
The natural frequencies of lateral vibration of beams are found from equation (9) given by
ωn = (βnL)2√(EI/ρL4) ------ (16) Where the number βn depends on the boundary conditions of the beam.
The following table lists the numerical values of (βnL)2 for typical end boundary conditions of the beam.
Dept. Of ME , ACE 11
Mechanic al Vibrations [10ME72]
Example-1: Determine the first three natural frequencies of a Rectangular Cantilever Beam for the
following data Length L = 1 m Breath B = 0.03 m Depth D = 0.04 m Young's Modulus
E = 2x1011 N/m2 Density r = 7850.0 Kg/m3
Solution: For a cantilever beam the boundary conditions are at x= 0 displacement and slope are zero
i.e. y = 0 and dy/dx = 0 ------ (1) and at x= L the shear force ‘V’ and the bending moment ‘M’ is zero
i.e. d2y/dx2 = 0 and d3y/dx3 = 0 -------- (2)
Since β2 ≠ 0 and A = -C and B = -D and we have A(coshβL + cosβL) + B(sinhβL + sinβL) = 0 --- (6)
Again Since β2 ≠ 0 and A= -C and B= -D and we have A(sinhβL - sinβL) + B(coshβL + cosβL) = 0 ---- (8)
From equation (7) and (9) we have A(coshβL+cosβL) + B(sinhβL+sinβL) = 0 --- (10)
Equation (11) is satisfied by a number of values of βL, corresponding to each normal mode of vibration,
which for the first three modes these values are tabulated in the table as
Circular Frequency
Natural Frequencies
Example-2: Determine the first three natural frequencies of a Rectangular Cantilever Beam for the
following data Length L = 1 m Breath B = 0.03 m Depth D = 0.04 m Young's Modulus E = 2.00E+11 N/m2
Density r = 7850.0 Kg/m3
Solution: For a cantilever beam the boundary conditions are at x= 0 displacement and slope are zero
and at x= L the shear force ‘V’ and the bending moment ‘M’ is zero
Equation (11) is satisfied by a number of values of βL, corresponding to each normal mode of vibration,
which for the first three modes these values are tabulated in the table as
Circular Frequency
Natural Frequencies