NPTEL Mechanical Engineering Nonlinear Vibration

Joint initiative of IITs and IISc Funded by MHRD Page 1 of 104

Module 6

APPLICATIONS
NONLINEAR VIBRATION OF MECHANICAL SYSTEMS



6
Applications
31-42
SDOF Free and Forced Vibration: Duffing Equation,
van der pols Equation: Simple or primary resonance,
sub-super harmonic resonance.

Parametrically excited system- Mathieu-Hills
equation, Floquet Theory, Instability regions; Multi-
DOF nonlinear systems and Continuous system,
System with internal resonances
6 A
Free
Vibration of
nonlinear
Systems

1 Single degree of freedom Nonlinear conservative systems with Cubic
nonlinearities.
2 Single degree of freedom nonlinear conservative systems with quadratic
and Cubic and nonlinearities.
3 Single degree of freedom non-conservative systems: viscous damping,
quadratic and Coulomb damping
4 non-conservative systems: Negative damping, van der Pol oscillator,
simple pendulum with quadratic damping
6 B
Forced
nonlinear
Vibration

1 Single degree of freedom Nonlinear systems with Cubic nonlinearities:
Primary Resonance
2 Single degree of freedom nonlinear systems with Cubic nonlinearities:
Nonresonant Hard excitation
3 Single degree of freedom Nonlinear systems with Cubic and quadratic
nonlinearities and self sustained oscillations
4 Multi-degree of freedom nonlinear systems
6 C
nonlinear
Vibration
of
Parametrically
excited
system
1 Parametrically excited system: Floquet theory, Hills infinite determinant
2 Parametric Instability region: sandwich beam vibration
3 Base excited magneto-elastic cantilever beam with tip mass
4 System with internal resonance: Two-mode interaction: Base excited
cantilever beam with tip mass at arbitrary position
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 2 of 104


Module 6 Lect 1
Free Vibration of Nonlinear Conservative system
In this lecture we will learn about the free vibration response of nonlinear conservative systems.
Initially the qualitative analysis will be demonstrated and later by using one of the perturbation
methods the free vibration response of the single degree of freedom system will be illustrated using
different examples.
Qualitative Analysis of Nonlinear Systems
Consider the nonlinear conservative system given by the equation
( ) 0 u f u + = (6.1.1)
Multiplying u in Eq. (6.1.1) and integrating the resulting equation one can write as
( ) ( )
( ( ))
Or, ( )
Or, ( )
uu uf u dt h
d u d u
u dt f u dt h
dt dt
udu f u du h
+ =
+ =
+ =

(6.1.2)
Or.,
1
2
( ) , where, ( ) ( )
2
u F u h F u f u du + = =

(6.1.3)
This represents that the sum of the kinetic energy and potential energy of the system is constant. Hence,
for particular energy level h, the system will be under oscillation, if the potential energy ( ) F u is less
than the total energy h. From the above equation, one may plot the phase portrait or the trajectories for
different energy level and study qualitatively about the response of the system using the following
equation.
( ) 2 ( ) u h F u = (6.1.4)
It may be noted that velocity exists, or the body will move only when ( ) h F u > . One will obtain
equilibrium points corresponding to ( ) h F u = or when ( ) ( ) 0 F u f u = = . For minimum potential energy
a center will be obtained and for maximum potential energy a saddle point will be obtained. The
trajectory joining the two saddle points is known as homoclinic orbit. The response is periodic near the
center.
Example 6.1.1 Perform qualitative analysis of spring-mass system with a soft spring. Take mass of the
system as 1unit, linear stiffness 1 unit and stiffness corresponding to cubic nonlinear term as 0.1 unit.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 3 of 104


Solution: In this case the equation of motion of the system can be given by
3
0.1 0 x x x + = (6.1.5)
For this system
3 2 4
1 1
( ) ( ) ( 0.1 )
2 40
F x f x dx x x dx x x = = =

(6.1.6)
Figure 6.1.1 shows the variation of potential energy F(x) with x. It has its optimum values
corresponding to
( ) ( ) ( ) 0 F x f x = = 0 or 10 x = . While x equal to zero represents the system with
minimum potential energy, the other two points represent the equilibrium points with maximum
potential energy. Now by taking different energy level h , one may find the relation between the
velocity v and displacement x as
2 4
2( ( )) 2( (0.5 0.025 )) v x h F x h x x = = =

Fig. 6.1.1: Potential well (curve with blue colour) and phase portrait (red colour) showing saddle point
(S) and center (C) corresponding to maximum and minimum potential energy.
Now by plotting the phase portrait one may find the trajectory which clearly depicts that the
equilibrium point corresponding to maximum potential energy is a saddle point (marked by point S)
where the equilibrium point is unstable and the equilibrium point corresponding to the minimum
potential energy is stable center type (marked by point C). Clearly the orbit joining the points S and S is
homoclinic orbit. Depending on different initial conditions i.e the total energy of the system, near the
center one will obtain periodic orbits.
Example 6.1.2: Perform qualitative analysis for a simple pendulum.
Solution: The governing equation of motion of the system is given by (Eq. 2.1.7)
S
C
S
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 4 of 104


sin 0 or sin 0
g
ml mg
l
+ = + =

(6.1.7)
So, ( ) ( ) sin cos
g g
F f d d
l l
= = =

(6.1.8)
Hence, the potential function ( ) F has a minima corresponding to equal to zero and maxima
corresponding to equal to odd multiple of . So, the equilibrium point near 0 = is a center and near
180 =

is a saddle point. The phase portrait for different energy level h is given in Figure 6.1.3 which
is obtained by using the following relation.
2( ( )) 2 cos
g
v h F h
l
| |
= = = +
|
\ .

(6.1.9)

Figure 6.1.3: Phase portrait for the motion of a simple pendulum



Approximate solution method
C
C
C C C
S S S S
Fig. 6.1.2: Simple pendulum

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 5 of 104


Let us consider the Duffing equation with cubic nonlinear term for the free vibration study of a
nonlinear system.

2
2 3
0 2
0
d u
u u
dt
+ + = (6.1.10)
Using method of multiple scales, the solution of this equation can be given by

2 3
0 1 2
( ) u u u u o = + + +
(6.1.11)

Taking the time scale =
n
n
T t (6.1.12)
2
0 1 2
................,
n
n
d
D D D D
dt T


= + + + =

(6.1.13)
( )
2
2 2 2
0 0 1 1 0 2 2
2 2 ................
d
D D D D D D
dt
= + + + + (6.1.14)

Substituting Eqs. (6.1.11-6.1.14) in Eq. (6.1.10) and separating terms with different order of one
obtains the following equations.

Order of
0

2 2
0 0 0 0
0 D u u + = (6.1.15)

Order of
1


2 2 3
0 1 0 1 0 1 0 0
2 D u u D Du u + = (6.1.16)
Order of
2

2 2 2 2
0 2 0 2 0 1 1 0 2 0 1 0 0 1
2 2 3 D u u D Du D D u D u u u + = (6.1.17)

Solution of Eq. (6.1.15) can be given by

0 0 0 0
0 1 2 1 2
( , ,) ( , ,)
i T i T
u A T T e A T T e

= + (6.1.18)

Substituting Eq.(6.1.18) in Eq. (6.1.16) one obtains
0 0 0 0
2
0
3 2 2 3
1 0 1 0 1
2 3 [ ]

+ = + +
i T i T
i D A A D e A A u u e cc (6.1.19)
To eliminate the secular term (coefficient of
0 0
i T
e ) term marked in pink colour should be equated to
zero.

Hence,
2
0 1
2 3 0 i D A A A + = (6.1.20)
From Eq. (6.1.19), the solution of
1
u can be written as
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 6 of 104


0 0
3
3
1 2
0
8
i T
A
u e cc

= + (6.1.21)
Now Substituting
1
exp( )
2
A a i = (where a and are real number) in Eq. (6.1.20) and separating the
real and imaginary parts one obtains
1
0,
a
T

and
3
0
1
3
0
8
a a
T

+ =

(6.1.22)

Hence, a is not a function of
0
T and
1
T . So up on integration Eq. (6.1.22) can be written as

2
( ), a a T = and
2
1 0 2 2
0
3
( )
8
a T T

= + (6.1.23)
Substituting Eq. (6.1.18) and Eq. (6.1.21) in Eq. (6.1.17), one can write
0 0 0 0 0 0
Secular t
3 5 2 2 4 5
0 2 0 2 2
e
3 2
0 2
2
0 0 0
m
2
r
21 3
8
1
2
8 8
5
i T i T i T
D
A A
i u u e A Ae A e c D A c

| |
+ = + +
|
\ .

(6.1.24)

For the secular term to be zero one can write
3 2
0 2 2
0
15
2
8
A A
i D A

(6.1.25)

The solution of remaining part of Eq. (6.1.24) can be written as
0 0 0 0
5
5 3 4
2 4 4
0 0
21
64 64
i T i T
A
u e A Ae cc


= + + (6.1.26)
Now Substituting
1
exp( )
2
A a i = (where a and are real number) in Eq. (6.1.25) and separating the
real and imaginary parts one obtains

2
0,
a
T

and
4
0 2
2 0
15
256
a
T

(6.1.27)
So, a is a constant. Now using Eqs. (6.1.23, 6.1.27) one can write
4
0 2 3
0
15
256
a T

= + (6.1.28)
Hence, from Eq. (6.1.23)

2 4
1 2 3
0 0
3 15
8 256
a T a T

= + (6.1.29)
Hence the solution of the system can be written as
2 3 2 5
3
2 2 4
0 0 0
21
1 cos( ) cos3( ) cos5( ) ( )
32 32 1024
a a a
u a t t t o



| |
= + + + + + +
|
\ .
(6.1.30)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 7 of 104


where,
2 4
2 3
0 3
0 0
3 15
( )
8 256


= + +
a a
o (6.1.31)
From Eq. (6.1.31) it may be noted that the frequency of oscillation is a function of amplitude
oscillation. It may be recalled that in case of linear system frequency is independent of amplitude of
oscillation.

For, the simple pendulum taking
0
/ g L = , where g is the acceleration due to gravity and L is the
length of the pendulum (=1 m) , the variation of frequency with amplitude is shown in Fig. 6.1.4. The
phase portrait for the periodic response considering only the first order solution is shown in Fig. 6.1.5.






Exercise Problems:

1. Perform qualitative analysis for the following nonlinear systems. Using equation (6.1.30) plot the
response of the system and compare both the results.

(a) +4+
3
=0 (c) +9+0.5
2
+0.1
3
=0
(b) ++
2
=0 (d) +100+10
3
=0
(e) ++0.1
3
+0.05
3
=0





Fig. 6.1.4: Frequency amplitude
relation for a simple pendulum
Fig. 6.1.5: Phase portrait showing the periodic
motion for a simple pendulum
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 8 of 104






Module 6 Lecture 2
Free vibration of nonlinear single degree of freedom conservative systems with quadratic and
cubic nonlinearities.

In this lecture the free vibration response of a nonlinear single degree of freedom system with
quadratic and cubic nonlinearities will be discussed with numerical examples. As studied in module 3,
the equation of motion of a nonlinear single degree of freedom system with quadratic and cubic
nonlinearities can be given by

2 2 3
0 2 3
0 x x x x + + + = (6.2.1)
Here
0
is the natural frequency of the system
2
and
3
are the coefficient of the quadratic and cubic
nonlinear terms. Also is the book-keeping parameter which is less than 1. Using method of multiple
scales the solution of this equation can be written as
( ) ( ) ( ) ( ) ..... ,...... , , ,...... , , ,...... , , ;
2 1 0 3
3
2 1 0 2
2
2 1 0 1
+ + + = T T T x T T T x T T T x t x (6.2.2)
Using different time scales , ,
1 0
T T and
2
T where
n
n
T t = and Eq. (6.2.2) in Eq. (6.2.1) and separating the
terms with different order of one can write the following equations.
Order of
1

2 2
0 1 0 1
0 D x x + = (6.2.3)
Order of
2

2 2 2
0 2 0 2 0 1 1 2 1
2 D x x D D x x + = (6.2.4)
Order of
3

2 2 2 3
0 3 0 3 0 1 2 1 1 0 2 1 2 1 2 3 1
2 2 2 D x x D D x D x D D x x x x + = (6.2.5)
The solution of (6.2.3) can be written as

1 1 2 0 0 0 0
( , )exp( ) exp( ) x A T T i T A i T = + . (6.2.6)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 9 of 104


Here A is an unknown complex function and A is the complex conjugate of A. Substituting (6.2.3)
into (6.2.4) leads to

Secular term
2 2 2
0 2 0 2 2 0 0 0 1 0 0
exp(2 ( ) ) 2 exp D x x A i T AA cc i D A i T
(
+ = + +

(6.2.7)
Here cc denotes the complex conjugate of the preceding terms. To have a bounded solution one should
eliminate the secular term and hence
1
1
0
dA
D A
dT
= = (6.2.8)
Therefore A must be independent of
1
T . With
1
0 D A = the particular solution of (6.2.7) can be written
as

2
2 2
2 0 0 2 2
0 0
exp(2 )
3
A
x i T AA cc

= + (6.2.9)
Substituting the expression for
1
x and
2
x from equation (6.2.6) and (6.2.9) into (6.2.5) and recalling
that
1
0 D A = we obtain
( )
2 2
2 2 3 0
0 2 0 0
2
0
2 2
0 3 0 3
2 2
3 3 0 2
0 0
2
0
Secular Te
3 2
exp
10 9
2 exp
3

rm
(3 )
3
i D A A A i D x x
A i
T
T cc

+ =
+
+
(

(6.2.10)
To eliminate the secular terms from
3
x , we must put

2 2
2 3 0 2
0 2 2
0
9 10
2 0
3
i D A A A

+ = (6.2.11)
Using
1
exp( )
2
A a i = where a and are real function of
2
T in Eq. (6.2.11) and separating the result
into real and imaginary parts, one obtains
0 a = and
2 2
3 2 3 0
0 2
0
10 9
0
24
a a

+ = (6.2.12)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 10 of 104


where the prime denotes the derivative with respect to
2
T . As 0 a = , a is a constant and
2 2
3 2 3 0
2
0 0
10 9
24
a
a

= or
2 2
2 3 0 2
2 0 3
0
9 10
24
a T

= + (6.2.13)



Here
0
is a constant. Now using
2
2
T t = one may write

2 2
2 2 3 0 2
0 3
0
9 10 1
exp
2 24
A a i a t i

(
= +
(

(6.2.14)
Substituting Eq. (6.2.14) in the expressions for
1
x and
2
x in Eqs. (6.2.6), (6.2.9) and (6.2.2), one
obtains
2 2
3 2
0 0 2
0
1
cos( ) 1 cos(2 2 ) ( )
2 3
a
x a t t O

(
= + + +
(

(6.2.15)
Here
2 2
2 2 3 3 0 2
0 4
0
9 10
1 ( )
24
a O

(
= + +
(

(6.2.16)
This solution is in good agreement with the solution obtained using the Lindstedt-poincare procedure.
[ Nayfeh and Mook, 1979].
Example 6.2.1: Taking

=k=100,
2
1 = and
3
1.5 = in eq
n
6.2.1find the response of the
system.
Using Eq. (6.2.16) the variation of frequency with amplitude is shown in figure (6.2.1). Taking two
values of a (viz., a=0.009 and a=.029) the time response has been plotted in figure (6.2.2). It may be
noted
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 11 of 104



Figure 6.2.1: Variation of amplitude with frequency for

Figure 6.2.2.(a): Time response (b) Phase portrait corresponding to initial amplitude a=0.009 and
a=.029
By changing the quadratic nonlinear terms
2
from .5 to 2.5 and keeping all other parameter same
figure (6.2.3) shows the variation of the frequency with amplitude of oscillation.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 12 of 104



Figure 6.2.3: Variation of frequency with amplitude for different values of
2


Figure 6.2.4: Variation of frequency with amplitude for different values of
3

2
.5 =
2
1.5 =
2
1 =
2
2.5 =
2
2 =
3
.5 =
3
1 =
3
1.5 =
3
2 =
3
2.5 =
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 13 of 104


In Figure 6.2.3, by varying the amplitude form -1 to +1 , the frequency of the system changes and it
decreases from 1.005 to .98 by increasing
2
from .5 to 2.5. One may observe the reverse phenomenon
by increasing
3
. Also the orientation of the frequency response curve changes from left to right when
3
is greater than 1 or
2
1 . It may be noted that for a linear system the response frequency of the
system does not depend on the amplitude of the response. But in case of the nonlinear system it
depends on the response amplitude. A Matlab code is given below which may be used to obtain the
frequency response, time response, and phase portrait for different system parameters.
%NPTEL WEB MODULE6 L2--Free Vibration: Duffing Oscillator
% Response plot using the method of multiple scales
clear all
clc
omega=1;
alpha1=omega^2;
alpha2=1;
alpha3=1.5;
ep=0.1;
p1=(9*alpha3*alpha1-10*alpha2^2)/(24*alpha1^2)
i=1
for a=0:0.001:1
om=sqrt(alpha1)*(1+p1*ep^2*a.^2)
s(i,1)=a;
s(i,2)=om;
i=i+1;
end
figure(1)
plot(s(:,2),s(:,1),s(:,2),-s(:,1),'linewidth',2)
grid on
set(gca,'FontSize',15) % For changing fontsize of tick no
xlabel('\bf frequency','Fontsize',15)
ylabel('\bf amplitude','Fontsize',15)
n1=length(s)
bt0=0;
for ii=10:20:40
om1=s(ii,2)
a1=s(ii,1)
T=2*pi/om1;
jj=1;
for t=0:T/1000:10*T
p2=ep^2*a1^2*alpha2/(2*alpha1);
p3=1-(1/3)*cos(2*om1*t+2*bt0);
x(jj)=ep*a1*cos(om1*t+bt0)-p2*p3;
tt(jj)=t;
jj=jj+1;
end
n2=jj-1
for k=2:n2
xt(k-1)=(x(k)-x(k-1))/(tt(k)-tt(k-1));
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 14 of 104


end
figure(2)
plot(tt,x,'linewidth',2)
hold on
grid on
set(gca,'FontSize',15) % For changing fontsize of tick no
xlabel('\bf Time','Fontsize',15)
ylabel('\bfx','Fontsize',15)
figure(3)
plot(x(1:n2-1),xt,'linewidth',2)
hold on
end
grid on
set(gca,'FontSize',15) % For changing fontsize of tick no
xlabel('\bf displacement','Fontsize',15)
ylabel('\bf velocity','Fontsize',15)

Exercise problem 6.2.1:
Find the frequency response of a single degree of freedom system with mass=1 kg, stiffness=100 N/m,
nonlinear quadratic and cubic stiffness parameter equal to 20 N/m
2
and 10 N/m
3
respectively. Vary the
book-keeping parameter and study the variation of frequency with amplitude.




Module 6 Lecture 3
FREE VIBRATION OF NONLINEAR SINGLE DEGREE OF FREEDOM
NONCONSERVATIVE SYSTEMS

In this lecture discussion on the vibration of a linear single degree of freedom system with viscous,
Coulomb damping, quadratic dumping and will be carried out using method of averaging.

System with viscous damping

Let us consider a single degree of freedom system with viscous damping. The equation of motion of
this system with massm, stiffness k and damping c can be written as

0 mu ku cu + + = . (6.3.1)
The same system can be written using the term natural frequency
n
, damping ratio as
2
2 0
n n
u u u + + = (6.3.2)
( )
2
Or, , 2 2
n n
u u f u u u u + = = = (6.3.3)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 15 of 104


By using Krylov-Bogoliubov method of averaging for an under damped system ( 1 < ) the solution can
be written as
sin( )
n
u a t = + (6.3.4)
where
( ) ( )
2
0
sin cos , sin
2
n
n
a f a a d

(6.3.5)
( ) ( )
2
0
cos cos , sin
2
n
n
f a a d
a

(6.3.6)
Substituting expression for
( ) , f u u from Eq. (6.3.3) in Eq. (6.3.5) and Eq. (6.3.6), one obtains
2
2
0
sin
a
a d a

= =

(6.3.7)
and sin cos 0 d

= =

. (6.3.8)
Solving Eq. (6.3.7) and Eq. (6.3.8) yields
( ) ( )
0 0 0
exp exp ,
n
a a t a t = = = (6.3.9)
Here
0
a and
0
are the initial displacement and phase of the response. Substituting Eq. (6.3.9) in Eq.
(6.3.4) one obtains the following equation.
( ) ( ) ( )
0 0
exp cos
n n
u a t t O = + + (6.3.10)

This equation is same as the expression one may obtain by finding the complementary function of the
differential equation (6.3.2). Using the
0
u and
0
u as the initial displacement and velocity respectively,
one may write Eq. (6.3.10) as
( ) ( ) ( )
0 0 0
s exp cos i / n
n d n d d
u t u t u u t
(
= + +

(6.3.11)
Where the damped natural frequency
2
1
d n
=
For over damped ( 1 > ) system one may use the following expressions for the response.
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
2 2 2
0 0
2 2 2
0 0
/ 2 exp
/ 2 ex
1 1 1
1 p 1 1
n n n
n n n
u u u t
u u t


= + + + +
+

+

(6.3.12)
For critically damped ( 1 = ) system one may write the response as follows.
( ) ( ) ( )
0 0 0
exp
n n
u u u u t t = + + (6.3.13)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 16 of 104


A Matlab code is given below to plot the under damped, critically damped and over damped response
of a system as shown in figure 6.3.1.

Figure 6.3.1(a): Time response of a linear single degree of freedom with viscous damping.




Figure 6.3.1(b): Time response of the system with linear damping.
( )
0 0
3, 1, 0.5, 0.09, 3.15
n
a = = = = =

Using Eq. (6.3.10) the time response is shown in Fig. 6.3.1(b). It may be noted that the response decreases exponentially.
The corresponding Matlab code is given below

Under damped
Critically damped
Over damped
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 17 of 104


Matlab code 6.3.1:
%Fr ee Vi br at i on r esponse of a l i near si ngl e degr ee of f r eedomsyst em

x0=0. 1;
xt 0=0. 001;
wn=2;
zet a=1. 5;
t =0: 0. 001: 20;

%over damped

z1=- zet a+sqr t ( zet a^2- 1) ;
z2=- zet a- sqr t ( zet a^2- 1) ;
z3=2*wn*sqr t ( zet a^2- 1) ;
A=( xt 0- z2*wn*x0) / z3;
B=( - xt 0+z1*wn*x0) / z3;
x1=A*exp( z1*wn*t ) +B*exp( z2*wn*t ) ;

%cr i t i cal l y damped

x2=( x0+( xt 0+wn*x0) *t ) . *exp( - wn*t )

%under damped

zt =0. 2 %Dampi ng f act or
wd=wn*sqr t ( 1- zt ^2) ;
x3=exp( - zt *wn*t ) . *( ( ( xt 0+zt *wn*x0) / wd) . *si n( wd*t ) +x0*cos( wd*t ) ) ;

pl ot ( t , x1, ' r ' , t , x2, ' b' , t , x3, ' g' , ' l i newi dt h' , 2)
gr i d on

set ( gca, ' Font Si ze' , 15) %For changi ng f ont si ze of t i ck no
xl abel ( ' \ bf Ti me' , ' Font si ze' , 15)
yl abel ( ' \ bf x' , ' Font si ze' , 15)

Matlab code 6.3.2:
%pl ot t i ng of l i near dampi ng ( Eq. 6. 3. 10) .
cl c
cl ear al l
a0=3;
ep=. 5;
mu=. 09;
t =0: 0. 1: 100;
omega=1;
bet a=- 3. 15;
a=a0*exp( - ep*mu*t ) ;
u=a0*exp( - ep*mu*t ) . *cos( omega*t +bet a) ;
pl ot ( t , u, t , a, ' - - ' , t , - a, ' - - ' )
%t i t l e( ' SYSTEM WI TH LI NEAR DAMPI NG' )
set ( f i ndobj ( gca, ' Type' , ' l i ne' ) , ' Col or ' , ' b' , ' Li neWi dt h' , 2) ;
set ( gca, ' Font Si ze' , 14)
xl abel ( ' t ' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
yl abel ( ' u' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
gr i d on
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 18 of 104




Single degree of freedom system with quadratic damping.
Here, the equation of motion of the system can be written as
( )
2
, | |
n
u u f u u u u + = = (6.3.14)
Similar to viscous damping here also using KB method the solution can be written as
sin( )
n
u a t = + (6.3.15)
Here, a and can be given by Eq. (6.3.5) and (6.3.6). Now using the expression for
( ) , f u u in Eq.
(6.3.5) and (6.3.6), one may write


( ) ( )
2 2 2
2
0 0
sin cos , sin sin sin
2 2
n
n
n
a
a f a a d d




= =


=
2 2
3 3 2
0
4
sin sin
2 3
n
n
a
d d a

(
=
(


(6.3.16)
and
( ) ( )
2 2
0 0
cos cos , sin sin cos sin 0
2 2
n
n
n
a
f a a d d
a




= = =

(6.3.17)
Solving Eq. (6.3.16) and Eq. (6.3.17) one may write Eq. (6.3.15) as

( ) ( )
0
0
0
cos
4
1
3
n
n
a
u t O
a
t

= + +
+
(6.3.18)

Figure 6.3.2: Time response of the system with quadratic damping.
( )
0 0
2, 1, 0.1, 3.15
n
a = = = =
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 19 of 104



Using Eq. (6.3.18) the time response is shown in Fig. 6.3.2. It may be noted that unlike the linear system the response does
not decreases exponentially but decreases algebraically. The corresponding Matlab code is given in Matlab code 6.3.3.

Matlab code 6.3.3:
%pl ot t i ng of quadr at i c dampi ng. (Eq. 6.3.18)
cl c
cl ear al l
a0=2;
ep=. 1;
t =0: 0. 1: 80;
omega=1;
bet a=- 3. 15;
a=a0. / ( 1+( 4*ep*omega*a0*t ) / ( 3*pi ) ) ;
u=a. *cos( omega*t +bet a) ;
pl ot ( t , u, t , a, ' - - ' , t , - a, ' - - ' )
%t i t l e( ' SYSTEM WI TH QUADRATI C DAMPI NG' )
set ( f i ndobj ( gca, ' Type' , ' l i ne' ) , ' Col or ' , ' b' , ' Li neWi dt h' , 2) ;
set ( gca, ' Font Si ze' , 14)
xl abel ( ' t ' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
yl abel ( ' u' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
gr i d on

System with Coulomb damping

In this case the equation of motion of the system can be given by

0
for 0
sgn( )
for 0
c
c
mx kx F
N x
F N x
N x

+ + =
>
= =

<

(6.3.19)

Using method of averaging this equation can be written as
2
0
for 0
/
for 0
c
g x
x x f F m
g x

>
+ = = =

<

(6.3.20)
The solution of the above equation can be written as
0
cos( ) x a t = +
where
( ) ( )
2 2
0 0
2
sin cos , sin sin sin
2 2
n
n n n
g g
a f a a d d d

(
= = =
(


(6.3.21)
( ) ( )
2 2
0 0
cos cos , sin cos cos 0
2 2
n
n n
g
f a a d d d
a

(
= = =
(

(6.3.22)
Integrating Eq. (6.3.21) and Eq. (6.3.22) one may obtain
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 20 of 104


0 0
2
,
n
g
a a t

= = (6.3.23)
Substituting Eq. (6.3.23) in Eq. (6.3.20) the response of the system can be written as
( ) ( )
0 0
2
cos
n
n
g
x a t t O

| |
= + +
|
\ .
(6.3.24)


Figure 6.3.3: Time response of the system with Coulomb damping.
( )
0 0
2, 1.3, .1, 0.4, 3.15
n
a g = = = = =

Here it may be noted that the response of the system decreases linearly. A Matlab code is given below.


Matlab code 6.3.4:
%pl ot t i ng of t i me r esponse f or syst emwi t h Coul omb dampi ng. (Eq.6.3.24)
cl c
cl ear al l
a0=2;
ep=. 1;
mug=. 4;
t =0: 0. 1: 100;
omega=1. 3;
phi =omega;
bet a=- 3. 15;
a=a0- ( ( 2*ep*mug*t ) / ( pi *omega) ) ;
u=a. *cos( omega*t +phi ) ;
pl ot ( t , u, t , a, ' - - ' , t , - a, ' - - ' )
%t i t l e( ' SYSTEM WI TH COULOMB DAMPI NG' )
set ( f i ndobj ( gca, ' Type' , ' l i ne' ) , ' Col or ' , ' b' , ' Li neWi dt h' , 2) ;
set ( gca, ' Font Si ze' , 14)
xl abel ( ' t ' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
yl abel ( ' u' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
gr i d on

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 21 of 104


Exercise problem :

1. Find the response of a single degree of freedom system with mass 1 kg, stiffness 100 N/m and
damping factor 10 N.s/m. Plot the time response and phase portrait. Also plot the phase portrait
considering coulomb damping and quadratic damping. Develop a Matlab code for finding the time
response and phase portrait by using second order governing differential equation of motion (Use
Runge-Kutta method).
Hints-The Matlab code for the system with viscous damping is given below
Matlab code 6.3.5:

%Use Runge-Kutta method to obtain the response of a sdof vibrating system

m=input( mass of the system in kg = )
k=input( Stiffness of the system in N/m = )
c=input( damping factor of the system in N.S/m= )
u0=input(initial Displacement in m=
v0=input(initial velocity in m=
omega_n=sqrt(k/m),
zeta=c/(2*m*omega_n);

if (zeta>1)
display(over damp system)
u=
end
if(zeta==1)
display(critically damped system)

u=
end

if(zeta<1)
display(under damped system)
u=u0*sin(omega_n*t)+(v0/omega_n)*cos(omega_n*t)
end

plot(t,u)

[T,u]=ode45(@ex631f,[0,20],[0.1,0.01]

Xlabel
Ylabel
Title
function
Ex631f

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 22 of 104


2. Find the response of a single degree of freedom system with Hysteretic damping. The equation of
motion in this case is given by.


2
0
x x f + =

2 and 2
s b b c
s c d
s d a d
s d a
c b s a d s
x x x x x x
x x x x
f
x x x x x x
x x x x
x x x x x x
+

=
`



)
= = +


3. Find the response of a single degree of freedom system with material damping by considering (a)
Maxwell model (spring and dashpot) in series, (b) Kelvin-Viogot Model (spring and dashpot) in
parallel. Consider soft spring with cubic nonlinearity in both the cases.







Module 6 Lecture 4
FREE VIBRATION OF SYSTEMS WITH NEGATIVE DAMPING

In this lecture initially the system with negative damping will be discussed. Then the free vibration
response of systems similar to van der Pol type of oscillator will be discussed with the help of
numerical examples. Finally the nonlinear response of a simple pendulum with viscous damping will be
illustrated.
There are many systems which can be modeled as a system with negative damping. This type of system
particularly occurs in control system where the derivative gains if not properly adjusted will give rise to
negative damping. Also this type of damping can be found in the high voltage transmission lines. The
equation of motion for this type of system for example that of a Rayleigh oscillator can be given by
( )
2 3
0
u u f u u + = = (6.4.1)
There are many systems which can be modeled as a system similar Rayleigh oscillator.
Using KB method the response amplitude a and phase of the system can be given by
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 23 of 104


( ) ( ) ( )
2 2
2 2 2 4
0 0
sin cos , sin sin sin
2 2
n n
n
a
a f a a d a d




= =


2 2
1 3
1
2 4
n
a a
| |
=
|
\ .
(6.4.2)
( ) ( ) ( )
2 2
2 2 2
0
0 0
cos cos , sin 1 sin sin cos 0
2 2
n
n
f a a d a d
a




(
= = =
(

(6.4.3)
Solving Eq. (6.4.2) one can write
( )
2
2 0
2 2 2 2
0 0 0 0
3 3
1 exp
4 4
a
a
a a t
=
| |
+
|
\ .
(6.4.4)
The time responses obtained by using Eq. (6.4.4) are shown in figure 6.4.1 for large and small initial
disturbance condition. It is observed that while with large initial disturbance the response decreases to

Figure 6.4.1: Time response of the system with Rayleigh damping. (a) Large initial disturbance(
0
2 a = )
(b) Small initial disturbance (
0
0.2 a = );
( )
0
1.3, .1, 0.4, 3.15
n
g = = = =
attain the steady state periodic response, and in case of small initial disturbance the response grows to
attain the steady state periodic response. Matlab code 6.4.1 may be used for plotting the time response
of the system with Rayleigh damping.


Matlab code 6.4.1:

%pl ot t i ng of t i me r esponse of t he syst emwi t h Rayl ei gh dampi ng. Equat i on no- 6. 4. 4.
cl c
cl ear al l
a0=2;
%a0=. 2;
ep=. 1;
t =0: 0. 1: 80;
omega=1;
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 24 of 104


bet a=- 3. 15;
a1=0. 75*( omega*a0) ^2;
a2=a0. ^2. / ( a1+( 1- a1) *exp( - ep*t ) ) ;
a=sqr t ( a2) ;
u=a. *cos( omega*t +bet a) ;
pl ot ( t , u, t , a, ' - - ' , t , - a, ' - - ' )
%t i t l e( ' SYSTEM WI TH Rayl ei gh DAMPI NG' )
set ( f i ndobj ( gca, ' Type' , ' l i ne' ) , ' Col or ' , ' b' , ' Li neWi dt h' , 2) ;
set ( gca, ' Font Si ze' , 14)
xl abel ( ' t ' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
yl abel ( ' u' , ' f ont si ze' , 14, ' f ont wei ght ' , ' b' ) ;
gr i d on

THE VAN DER POL OSCILLATOR

There are many systems which can be modeled as a system similar to van der Pols oscillator which is
named after the Dutch physicist Balthasar van der Pol (27 J anuary 1889 6 October 1959). Mostly this
equation is used in electrical circuits but it can also be used in some mechanical system where self
oscillation takes place due to negative damping.
The van der Pols equation can be written as
2
2
2
(1 )
d u du
u u
dt dt
+ = (6.4.5)
Using method of multiple scales the solution of this equation can be given by
( ) ( ) ( )
0 0 1 1 0 1
; , , ....... u t u T T u T T = + + (6.4.6)
Where , 0,1,2,
n
n
T t n = = . Using Eq. (6.4.6) in Eq. (6.4.5) and separating the terms with different
order of , one obtains the following equations.
2
0 0 0
0 D u u + =
(6.4.7)
2 2
0 1 1 0 1 0 0 0 0
2 (1 ) D u u D Du u D u + = +
(6.4.8)

2 2 2 2
0 2 2 0 1 1 1 0 0 2 0 0 0 1 0 1 0 0 1 0 0
2 2 (1 ) (1 ) 2 D u u D Du D u D D u u D u u Du u u D u + = + + (6.4.9)
The solution of Eq.(6.4.7) can be written as
0 0
0 1 2 1 2
( , ,) ( , ,)
iT iT
u A T T e A T T e

= +
(6.4.10)
Substituting Eq. (6.4.10) in Eq. (6.4.8)
0 0
3 2 2 3
0 1 1 1
secular term
(2 ]
iT iT
D u u i D A A A A e iA e cc + = + +
(6.4.11)
Eliminating the secular term marked in Eq. (6.4.11) one can write
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 25 of 104


2
1
2D A A A A =
(6.4.12)
Now the solution of Eq.(6.4.11) can be written as
0 0
3 3
1 1 2
1
( , )
8
iT iT
u B T T e iA e cc = + +

(6.4.13)

( )
1 2 1 2
1
( , )exp ( , )
2
A a T T i T T = (6.4.14)
Substituting Eq. (6.4.14) in Eq. (6.4.12) one can write
( ) ( ) ( ) ( ) ( )
2
1 1
1 1 1 1 1
2 exp exp exp exp exp
2 2 2 2 2
a
i ia i a i a i a i
T T

| | | |
+ =
| |

\ .
\ .
(6.4.15)
Separating the real and imaginary terms one can write

2
1 1
1 1
0, 1
2 4
a
a a
T T
| |
= =
|

\ .
(6.4.16)
Hence,
1
2
2
2
4
( ), and
1 ( )
T
T a
c T e


= =
+
(6.4.17)
So the first order solution of the system can be given by
cos ( ) u a t o = + (6.4.18)

Where
( )
2
2
0
4
4
1 1 exp
a
t
a

=
| |
+
|
\ .
(6.4.19)
To obtain the second or higher order solution one may use the expression for
0 1
and u u in Eq. (6.4.9)
which yields the following equation.
0 0
2
0 2 2 1 2 1 2
( , ,) ( , ,) NST
iT iT
D u u Q T T e Q T T e + = + + (6.4.20)

Where NST contains non-secular terms.
3 2
2 2 2
1 2 1 1 1
2 (1 2 ) 2 (1 2 )
8
A A
Q iD B i AA B iA B iD A D A AA D A A D A = + + + (6.4.21)
Now again substituting Eq. (6.4.14) in the secular term of Eq. (6.4.20) and separating the real and
imaginary parts of the resulting equation one obtains the following equation.

3
2
0, ( )

= =

a
a a T
T
(6.4.22)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 26 of 104


2
1 1 2 1
2 1 7 1
2 2
16 16 4
b da d da
b a a
T a dT dT dT
| | | |
= + +
| |

\ .
\ .
(6.4.23)
1
2
1 7 1
16 32 8
b d
d dT a da
a dT a
| |
| | | |
= + +
| | |
\ . \ .
\ .
(6.4.24)
Integrating one obtains
3
1 0 2
2
1 7 1
ln ( )
16 64 8
d
b a T a a a ab T
dT
| |
= + + +
|
\ .
(6.4.25)
In order that the solution to be bounded for all
1
T , the coefficient of
1
T in the above equation for b must
vanish. Hence one obtains
2
1
0
16
d
dT
| |
+ =
|
\ .
.

So ,
2 0
1
16
T = + (6.4.26)
Here
0
is a constant. Now using the expression for
0 1
and u u the second order solution can be given by

2 2
0 0
2 2
0
3 2
0
7 1 1
ln sin 1
64 8 16
1
cos 1 ( )
16
1 1
sin3 1
32 16
a a ab t
u a t o
a t



( | | | |
+ +
| | (
\ . \ . ( | |
= + +
`
| (
\ . ( | |

+ +
| (

\ . )
(6.4.27)
Here a can be given by the Eq. (6.4.19). Equation (6.4.27) can also be written as

3 2
1
cos( ) sin3( ) ( )
32
u a t a t o = + (6.4.28)
where
2 2
0
1 1 7
ln
16 8 64
t a a = + + (6.4.29)

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 27 of 104



Figure 6.4.2: Time response of the system with van der Pols equation
( )
0
1, .1, 3.15
n
= = =
(a) Small initial disturbance(
0
0.2 a = ) (b) Large initial disturbance (
0
3.5 a = ).

Example: 6.4.2 Simple pendulum with quadratic damping
The equation of motion of a simple pendulum with quadratic damping can be given by
2 2 3
0 0
1
2 sin 2 0
6

| |
+ + = + + =
|
\ .

(6.4.30)
Using method of multiple scales with different time scales , 0,1,2,
n
n
T t n = = and writing the motion
of the pendulum as follows
( ) ( ) ( ) ( ) ... , , , , , , ;
2 1 0 3
3
2 1 0 2
2
2 1 0 1
+ + + = T T T T T T T T T t

(6.4.31)
and substituting Eq. (6.4.31) in (6.4.29) and separating the terms with different order of one obtains
the following equations.
0
1
2
0 1
2
0
= + D (6.4.32)
1 1 0 2
2
0 2
2
0
2 D D D = + (6.4.33)
3
1 1 0 1 0 1
2
1 1 2 0 2 1 0 3
2
0 3
2
0
6
1
2 2 2 + = + D D D D D D D D (6.4.34)
The solution of Eq. (6.4.32) can be written as
( ) ( ) cc T i T T A + =
0 0 2 1 1
exp , (6.4.35)
Using Eq. (6.4.35) in (6.4.33) it can be written as
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 28 of 104


( ) ( )
2 2
0 2 0 2 0 1 1 0 1 0 0
2 2 exp D D D i D A i T cc + = = + (6.4.36)
As the terms in the right hand sides are secular terms, it will be eliminated if
1
0 D A = . Hence, Ais not a
function of
1
T . After eliminating the secular term, the solution of Eq. (6.4.36) will contain only
auxiliary part and hence,
2
may be dropped. So one can write
( ) ( )
1 2 0 0
exp A T i T cc = + (6.4.37)
Now substituting (6.4.37) in (6.4.34) one obtains
( ) ( ) ( ) ( )
( ) ( ) ( ) ( )
( ) ( )
' '
0 0 0 0 0 0 0 0 0 0 0
2 2
0 0 0 0 0
2 2
0 3 0 3
3 3
0 0 0 0
0 0 0 0 0
2 exp exp 2 exp exp
1 1
.exp exp ex
non secular
p exp
1 1
exp 3 exp 3
2 2

6 6
terms
D
A i T
i A i T A i T i A i T i A i T
i A i T i A i T A A i T AA i T
A i T

( (

+ =
+ +
+
+

(6.4.38)

Now substituting ( ) i a A exp
2
1
= in Eq. (6.4.38) one can write


2 2 ' ' 3 3 2 2
0 3 0 3 0 0 0
1 1
2 sin 2 cos cos cos3 2 sin sin
8 24
D a a a a a + = + + + + (6.4.39)
where
0
t = + . It may be noted that the damping term is periodic and can be expanded in Fourier
series as
1
sin sin sin
n
n
f n

=
=

(6.4.40)
Where
8
3
n
f

= . The secular terms will be eliminated from () if

0
16
, 0
3
8
0
2
' 2 0 '
= + = +

a
a a (6.4.41)
Solving Eq. (6.4.41) one may write
( )
0
2 0 0
2
0
0
2
2 0 0
0
8 3 128
9
,
8 3
3

+
+
=
+
=
T a
a
T a
a
a (6.4.42)
Hence, the response of the system can be given by
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 29 of 104


( )
( )
3
2
0
0
0 0
2
0
0
2
0
0 0
0
128
3
8 3 128
9
cos
8 3
3

O
t
t
t
+
(

+
+
+
= (6.4.43)

EXERCISE PROBLEM:
1. Obtain the response of a simple pendulum with quadratic damping. The equation of motion for this
system can be written as follows.
(a)
( )
2
0
2 sin 0 + =

, (b)
( )
2
0
2 sin 0 + + =


Plot the phase portrait and discuss about the equilibrium solution.
2. Obtain the response of a simple pendulum with viscous damping. The equation of motion for this
system can be written as
2 2 3
0 0
1
2 sin 2 0
6

| |
+ + = + + =
|
\ .


Plot the phase portrait and discuss about the equilibrium solution. Compare the results with that
obtained in problem 6.4.1.

















Module 6 Lecture 5
Forced Vibration of single degree of freedom system with cubic nonlinearities
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 30 of 104


In this lecture, the response of a nonlinear single degree of freedom system with cubic nonlinearies will
be discussed considering a weak forcing function. The simplest form of this equation can be given by
the forced Duffing equation as follows.
2 3
0
2 cos u u u u f t + + + = (6.5.1)
As discussed in the previous lectures, in the absence of external force, the free vibration response
amplitude of such system is a function of the natural frequency
0
of the system. Similar to linear
vibration of the system here we may consider the behaviour of the system near the resonance condition,
i.e., when the external frequency is equal to the natural frequency of the system. This condition is
known as the primary resonance condition (
0
). In case of multi-degree of freedom system one
may reduce the system into a number of single degree of freedom system and follow a procedure as
outlined here.
To study the behaviour of the system near the primary resonance condition, one may use the detuning
parameter which represents the nearness of the external frequency to that of the natural frequency.
Hence one may write
0
= + (6.5.2)
Using method of multiple scales the solution of Eq. (6.5.1) can be written as
( ) ( ) ( )
0 0 1 1 0 1
; , , ....... u t u T T u T T = + + (6.5.3)
Where
n
n
T t = . Substituting Eq. (6.5.3) in Eq. (6.5.1) and separating terms with different order of ,
one obtains the following equations.
2 2
0 0 0
0 D u + =

(6.5.4)
2 2 3
0 1 0 1 0 1 0 0 0 0 0 0 1
2 2 cos( ) D u u D Du D u u f T T + = + +

(6.5.5)
The solution of Eq. (6.5.4) can be written as
( ) ( ) ( ) ( )
0 1 0 0 1 0 0
exp exp u A T i T A T i T = +

(6.5.6)
Substituting Eq. (6.5.6) in Eq. (6.5.5) one obtains the following equation.
( ) ( ) ( ) ( )
2 2 ' 2 3
0 1 0 1 0 0 0 0 0 0 0 1
Secular term
Mixed Secular term
1
2 3 exp exp 3 exp
2
D u u i A A A A i T A i T f i T T cc
(
+ = + + + + + (

(6.5.7)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 31 of 104


To eliminate the secular and near secular terms from Eq. (6.5.7), one can write
( ) ( )
2
0 1
1
2 ' 3 exp 0
2
+ + = i A A A A f i T .

(6.5.8)
Now substituting ( )
1
exp
2
A a i = in Eq.(6.5.8) and separating the real and imaginary parts following
reduced equations are obtained.

( )
1
0
1
' sin
2

= +
f
a a T

(6.5.9)
( )
3
1
0 0
3 1
' cos
8 2

=
f
a a T

(6.5.10)
To write these two equations in its autonomous form one may use
1
T = and obtained the
following equations.
0
1
' sin
2

= +
f
a a
(6.5.11)
3
0 0
3 1
' cos
8 2

= +
f
a a a

(6.5.12)
One should solve these two equations to obtain a and and can write the first order solution of the
system in the following form
( ) ( )
( ) ( )
0 0 1
0
cos ( ) cos ( )
cos ( ) cos ( )
u a t O a t T O
a t t O a t O


= + + = + +
= + + = +

(6.5.13)
Now for steady state as and a equals to 0, one can write Eq. (6.5.11-12) as

0
1
sin
2
f
a

=

(6.5.14)
3
0 0
3 1
cos
8 2
f
a a

=

(6.5.15)
Now eliminating from the above equations, one obtains
2
2
2 2 2
2
0 0
3
8 4
f
a a

(
| |
( + =
|
(
\ .


(6.5.16)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 32 of 104


Eq. (6.5.16) is a 6
th
order Polynomial in a , but quadratic Polynomial in . Hence, by solving this
quadratic equation, one can write the expression for the frequency response curve as follows.
1
2
2
2 2
2 2
0 0
3
8 4
f
a
a

| |
=
|
\ .

(6.5.17)
Hence, for a particular value of detuning parameter one can get different amplitude of the response and
the phase. Now to check the stability of the obtained steady state response, one can perturb the Eq.
(6.5.11) and (6.5.12) to obtain the following J acobian matrix.
2
0
0
0
2
0
0 0
3
8
9 1
8
a
a
J
a
a

| |

|
\ .
=
| |

|
\ .
(6.5.18)
Now to find the stability of the steady state response one can find the eigenvalues by finding the
determinant of the J I matrix. This leads to expression
2 2
2 2 0 0
0 0
3 9
2 0
8 8
a a


| || |
+ + + =
| |
\ .\ .
(6.5.19)
The system will be unstable when the real part of at least one of the eigenvalue becomes positive. This
gives rise to the following relation.
2 2
2 0 0
0 0
3 9
0
8 8
a a


| || |
= + <
| |
\ .\ .

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 33 of 104



Fig:6.5.1 Frequency response curves for =1,f=10, =10,=0.1


Fig:6.5.2 Frequency response curves with different values of





NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 34 of 104


References for further reading
1. J . A. Goitwald, L. N. Virgin and E. H. Dowell, Experimental mimicry of Duffing's equation,
J ournal of Sound and Vibration 158(3), 441461 , (1992)
2. R. E. Mickens, Mathematical and numerical study of the Duffing-harmonic oscillator, J ournal
of Sound and Vibration, 244, no.(3), 563567, 2001
3. S. B. Tiwari, B. N. Rao, N. S. Swamy, K. S. Sai, and H. R. Nataraja, Analytical study on a
Duffing-harmonic oscillator, J ournal of Sound and Vibration, 285, (4-5), 12171222, 2005.
4. H. Hu, Solution of a quadratic nonlinear oscillator by the method of harmonic balance,
J ournal of Sound and Vibration,293, (1-2), 462468, 2006
5. C. W. Lim, B. S. Wu, and W. P. Sun, Higher accuracy analytical approximations to the
Duffing-harmonic oscillator, J ournal of Sound and Vibration, 296, (4-5), 10391045, 2006.
6. T. zi and A. Yildirim, Determination of the frequency-amplitude relation for a Duffing-
harmonic oscillator by the energy balance method, Computers & Mathematics with
Applications, 54, (7-8), 11841187, 2007.
7. A. Belndez, C. Pascual, E. Fernandez, C. Neipp, and T. Belendez, Higherorder approximate
solutions to the relativistic and Duffing-harmonic oscillators by modified He's homotopy
methods, Physica Scripta,77, Article ID 065004, 14 pages, 2008.
8. J . Cai, X. Wu, and Y. P. Li, An equivalent nonlinearization method for strongly nonlinear
oscillations, Mechanics Research Communications, 32, 553560, 2005.





NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 35 of 104




Module 6 Lect 6
System with non resonant hard excitations
In the previous lecture a single degree of freedom nonlinear system is considered when the amplit
ude
of the external excitation (f) is one order less than the linear term (i.e.
2
0
). In the present lecture
the forcing term is assumed to be of same the order as that of the linear term. So the equation of motion
considered in this case is
2 3
0
2 cos u u u u f t + + + =
(6.6.1)
Following similar procedure of method of multiple scales, one may write
( ) ( ) ( )
0 0 1 1 0 1
; , , ... u t u T T u T T = + +
(6.6.2)

Now separating the terms with different order of one obtains the following equations.

2 2
0 0 0 0 0
cos D u u f T + =
(6.6.3)

2 2 3
0 1 0 1 0 1 0 0 0 0
2 2 D u u D Du D u u + =
(6.6.4)
The solution of Eq. (6.6.3) can be written as
( ) ( ) ( )
0 1 0 0 0
exp exp u A T i T i T cc = + +
(6.6.5)
Where
( )
2 2
0
2
f

=


It may be noted that unlike the previous lecture, where only the complementary part of the solution was
present, in this case both complimentary and particular integral parts are present in the solution of
0
u .
Now substituting Eq. (6.6.5) in Eq. (6.6.4) one obtains the following equation.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 36 of 104


( ) ( )
( ) ( ) ( )
( )
2 2 ' 2 2
0 1 0 1 0 0 0
3 3 2
0 0 0 0 0
2
0 0
Secular term
Mixed Secular term 1
Mixed
2 6 3 exp
exp 3 ex
Secular ter
p 3 3 exp 2
3 exp
m
2
2
D u u i A A A A A i T
A i T i T A i T
A i T

(
+ = + + +

+ + + (

+ (

( ) ( )
( )
2 2
0 0 0 0
2
0
Mixed Secular term 3
Mixed Secular t
3 exp 2 3 exp 2
2 3 6 ex
4
p
erm
A i T A i T
i AA i T cc





`
+ + + ( (



)
( + + +

(
6.6.6
)



It may be noted from Eq. (6.6.6) that when the exponent terms of the marked mixed secular terms are
equal to
0
a resonance condition will occur. Hence, resonance will be observed in the system when

0
0 0 0
0 0
0 0
(Primary resonance)
(Sub harmonic resonan

2 3
1
3
ce)
(Super harmonic resonance)
3
2 0

=
= =
= =
= =



For the non resonant case, i.e., when the external frequency is away from 0,
0
,
0
1
3
or 3
0
, from Eq.
(6.6.6) eliminating the secular terms yield the following equation.
( )
' 2 2
0
2 6 3 0 i A A A A A + + + =
(
6.6.7
)
Now using ( )
1
exp
2
A a i = in Eq.(6.6.7) and separating the real and imaginary parts following reduced
equations are obtained.
1 0 1 0
exp( ) exp( ) Da a a a T a t = = = (6.6.8)


and
2 2 2 2 2 2
0 1 1 0 0
1 1 1
3 3 3
8 8 8
aD a a a T a t
| | | | | |
= + = + + = + +
| | |
\ . \ . \ .
(6.6.9)


NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 37 of 104


( )
( )
( )
( )
( )
0
2 2
0
2 2
0 0 0
2 2
0
cos cos
1
exp( )cos 3 cos
8
f
u a t t O
f
a t a t t O

= + + +

| | | | | |
= + + + + +
| | |
\ . \ . \ .
(
6.6.10
)
The free oscillation solution decays with time and hence the steady state response consists of forced
solution only similar to the linear case.
Superharmonic Resonance (
0
1
3
)
To express the nearness of the external excitation frequency to one third of the natural frequency one
may use the detuning parameter ( ) as follows.
0
3 = +
(6.6.11)
Now to include the mixed secular (or nearly secular or small divisor) term 1in Eq. (6.6.6) in this
resonance condition one may write

( )
0 0 0 0 0 0 0 0 1
3 T T T T T T = + = + = +
(6.6.12)

Now to eliminate the secular and near secular terms from Eq. (6.6.6) one can write
( ) ( )
' 2 2 3
0 1
2 6 3 exp 0 i A A A A A i T + + + + =
(
6.6.13
)

Using ( )
1
exp
2
A a i = in Eq.(6.6.13) and separating the real and imaginary parts following reduced
equations are obtained.
( )
3
1
0
sin a a T

= (6.6.14)


( )
3
2 2
1
0 0
3 1
cos
8
a a a T



| |
= + +
|
\ .

(6.6.15)
Now to express the above equations in their autonomous form one may use the following
transformation.
1
T =
(6.6.16)
Hence, Eq. (6.6.14-15) can be written as

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 38 of 104


3
0
' sin a a

=
(6.6.17)

2 3
' 3
0 0 0
3 3
cos
8
a a a



| |
=
|
\ .
(6.6.18)
By solving the above two equations, one can obtain a and , and then can write the solution of the
system as
( )
( )
( )
0
2 2
0
cos 3 cos
f
u u a t t O

= = + +

(6.6.19)
For steady state as the time derivative terms vanish, Eq. (6.6.17) and (6.6.18) can be written as
3
0
sin a

=
(6.6.20)

2 3
3
0 0 0
3
3 cos
8
a a



| |
=
|
\ .
(6.6.21)
Now eliminating

from the above two equations, one can obtain a closed form equation which can be
used for finding the frequency response of the system.
2
2 2 6
2 2 2
2
0 0 0
3
3
8
a a



(
| |
( + =
|
(
\ .

(6.6.22)

Solving Eq. (6.6.22) one may write the relation between the detuning parameter and amplitude of the
response as follows.
1/2
2 2 6
2 2
2 2
0 0 0
3
3
8
a
a



| |
= +
|
\ .
(6.6.23)
Hence, in this resonance condition, the free oscillation term does not decay to zero inspite of the
presence of damping. Moreover, the nonlinearity adjusts the frequency of the free oscillation term to
exactly three times the frequency of the excitations so that the response is periodic. Since the frequency
of the free oscillation term is 3 times the frequency of excitation, such resonances are called super
harmonic resonances or overtones.
From Eq. (6.6.20) the peak amplitude of the free oscillation term is given by
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 39 of 104


3
0
p
a

= (6.6.24)


2 2 4
2 2
0 0
3
1
8
p

| |
= +
|
\ .
(6.6.25)

Subharmonic Resonance
0
3 =
When the external frequency is nearly 3 times the natural frequency of the system, using detuning
parameter one can write

0
3 , = +
(6.6.26)

Or, ( )
0 0 0 0 0
2 T T T = +
(6.6.27)

Using a similar procedure of method of multiple scales, to eliminate the secular terms from Eq. (6.6.6)
one can write
( ) ( )
' 2 2 2
0 1
2 6 3 exp 0 i A A A A A A i T + + + + = .
(
6.6.28
)
Using ( )
1
exp
2
A a i = in Eq.(6.6.13) and separating the real and imaginary parts following reduced
equations are obtained.
( )
2
1
0
3
' sin 3
4
a a a T

=
(6.6.29)

( )
' 2 2 2
1
0 0
3 1 3
cos 3
8 4
a a a a T



| |
= + +
|
\ .
(6.6.30)
Now to express the above equations in their autonomous form one may use the following
transformation.
1
3 T =
(6.6.31)

2
0
3
' sin
4
a a a

=
(6.6.32)

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 40 of 104


2
' 3 2
0 0 0
9 9 9
cos
8 8
a a a a



| |
=
|
\ .
(6.6.33)

( ) ( ) ( )
1
2 2
0
1
cos cos
3


(
= + +
(

u a t f t O
(6.6.34)
For steady state one can write

2
0
3
sin
4
a a

=
(6.6.35)
2
3 2
0 0 0
9 9 9
cos
8 4
a a a



| |
=
|
\ .
(6.6.36)
Now eliminating

from the above two equations, one can obtain a closed form equation which can be
used for finding the frequency response of the system.

2
2 2 2
2 2 2 4
2
0 0 0
9 9 81
9
8 16
a a a



(
| |
( + =
|
(
\ .

(6.6.37)
This shows either the system will have a trivial state response (i.e., a
=0
) and a non trivial response
which can be obtained by solving the following equation.
2
2 2 2
2 2 2
2
0 0 0
9 9 81
9
8 16
a a



(
| |
( + =
|
(
\ .

(6.6.38)
This equation is quadratic in
2
a and hence the solution can be written as
( )
1/2
2 2
a p p q =
(6.6.39)

Where
2 0
8
6
9
p

=
(6.6.40)

and
2
2 2
2 0
2
0
64 9
9
81
q



(
| |
( = +
|
(
\ .

(6.6.41)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 41 of 104


As q is always positive, so the nontrivial free oscillation amplitudes occur when
2
0 and p p q > . For
these conditions one should have
2 2
2 2 0
0 0
4 63
and 2 0
27 8



| |
<
|
\ .
(6.6.42)
So and should have same sign. From Eq. (6.6.42), for a given nontrivial solutions can exist only
if
2 2 2
0
2
0
2 63
8

(6.6.43)
Similarly for a given

, nontrivial solution can exist only if

1/2 1/2
2 2 2
2 2
0
63
63 63
4


| | | |
+
| |
\ . \ .
(6.6.44)
In the

plane the boundary of the region where nontrivial solutions can exist can be given by
1/2
2 2
2
0
63
63
4


| |
=
|
\ .
(6.6.45)
Exercise Problems
1. Study the response of the single degree of freedom system with both quadratic and cubic
nonlinearities. Consider primary, subharmonic and superharmonic resonance conditions. Use either
method of multiple scales or the method of normal forms.
References:
1. A. H. Nayfeh and D. T. Mook: Nonlinear Oscillations, Wiley, 1979
2. A. H. Nayfeh, Method of Normal Forms, Wiley,1993
Exercise problems
1. 1. Using numerical techniques plot the time response of the system given by Eq. (6.6.1) with
primary, subharmonic and superharmonic resonance. Compare the results with those obtained
from the perturbation method.
2. Study different resonance conditions for the Duffing equation with two frequency excitation
terms as given by the following equation.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 42 of 104


( ) ( )
2 3
0 1 1 1 2 2 2
2 cos cos u u u u f t f t + + + = + + +

Ans: One may observe following resonance conditions







Q No. 3 Plot the frequency response curves for the system with sub harmonic and super harmonic
resonance condition using equation 6.6.25 and 6.6.45.


Module 6 Lecture 7
Forced vibration Single-Degree of freedom system
In this lecture briefly following analysis will be carried out
Forced vibration of Single-Degree of freedom system with cubic and quadratic nonlinearies
System with self sustained oscillation

System with cubic and quadratic nonlinearies
In this case considering damping, cubic nonlinearity and the forcing parameters to be one order less than
quadratic nonlinearity which is one order less than the linear term, the equation of motion can be written as
2 2 2 2 3 2
0 2 3
2 cos u u u u u f t + + + + =
(6.7.1)
Following similar procedure as in the previous two lectures one may study the primary, super harmonic
and sub harmonic resonance conditions.
( )
( )
0 0 0 0
0 1 2 0 1 0 2
0 2 1 1 2 0 2 1 0 2 1
2 1 0 2 1 0
2 1 0 2 1
2 1 0 2 1 0
2 1 0 2 1 0
1 1
3 , ,, 1 2
3 2
1 1
3 3 ,
3 3
1
2 2 , 2 ,
2
9 3 , 3
1 5
, 5
3 3
7 2
7 , 2
3 3
7 5
7 , 5
3 3
n n m n m n
or or
or







+






NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 43 of 104


In case of primary resonance, taking
2
0
= + and
( ) ( ) ( ) ( ) , , , , , , ,
2
0 0 1 2 1 0 1 2 2 0 1 2
u t u T T T u T T T u T T T = + + (6.7.2)
and applying the usual procedure of method of multiple scales, one can obtain the following reduced equations
(Nayfeh and Mook, 1979).

0
' sin
2
f
a a

=
(6.7.3)

2 2
' 3 3 0 2
3
0 0
9 10
cos
24 2
f
a a a

= +
(6.7.4)
Comparing these two equations with those obtained from Duffing equation with only cubic
nonlinearities for primary resonance (Eq. 6.5.11 and 12), one may observe that both sets of equations
are identical if
2
2
3 2
0
10
9

= (6.7.5)
When
2
2
3 2
0
10
0
9

= = which is negative, the system will show softening effect. Also if

2
2
3 2
0
10
, 0
9

< < one obtains softening effect in which the frequency response curves bends towards the
lower frequencies irrespective of the sign of
2
.
When
2
2
3 2
0
10
,
9

= 0 = , the system behaviour will be similar to that of a linear system. Here the effect of
cubic nonlinearities will be cancelled by that of the quadratic nonlinearity. Similar to the expression in the cubic
nonlinear system here the following equation can be used to find the frequency response.
2
2
2 2 2
2
0 0
3
8 4
f
a a

(
| |
( + =
|
(
\ .

(6.7.8)
The expression for the response can be given by
( ) ( ) ( )
2 2 2 2 2
2 2
0 0
1 1
cos cos 2 2
2 6
u a t a a t O



= + + + (6.7.9)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 44 of 104


Comparing this expression with that of the system with only cubic nonlinearity, it can be observed that the
oscillating motion is not centered at 0 u = and there is a drift or shift of the steady state part by an amount of
2 2
2
0
1
2
a

.
Superharmonic resonance
This resonance condition may occur when one consider, the forcing term is same order that of the linear part.
Also considering the damping and quadratic nonlinearity of the order of and cubic nonlinearity of the order
of
2
one may write the equation of motion as (Nayfeh and Mook, 1979)
2 2 2 3
0 2 3
2 cos u u u u u f t + + + + =
(6.7.10)
Applying method of multiple scales while eliminating the secular term, one may observe that a
resonance condition will occur when
0
2 = + . In this case the steady state solution can be written as
( )
( ) ( )
( ) ( )
2
2 2
2 2 2 1/2
2 2 2 2
0
0 0
cos sin 2
4
f f
u t t O

= +

(6.7.11)
1
tan

| |
=
|
\ .

Subharmonic resonance
For the system given by equation (6.7.10), subharmonic resonance will occur when external frequency is nearly
equal to
0
2 = + . In this case one may obtain the expression for the frequency of resulting oscillation as
1/2
2 2 2
2
2
0
4

| |
=
|
\ .
(6.7.12)
Hence if
2 2
2 2
2
0
4

> , the motion is oscillatory and decays with time

If
2 2 2 2
2 2 2 2
2 2
0 0
4 4 4



| |
> >
|
\ .
, the response decays without oscillating
If
2 2
2 2 2
2
0
4 4

| |
>
|
\ .
, the system becomes unstable as the response grows

Systems with self-sustained oscillations
Let us consider the van der Pols oscillator with soft harmonic excitation which is given by the
following equation.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 45 of 104


3 2
0
1
cos
3
f t u u u u
| |
=
|
\ .
+
(6.7.13)
Let us consider the primary resonance case in which the external frequency is assumed to be near the
linear system frequency
0
. So using detuning parameter one may write
0
t = + (6.7.14)
Now to solve Eq. (6.7.14), taking
( ) ( ) ( ) ( ) , , , u t u T T u T T o

= + +
2
0 0 1 1 0 1
(6.7.15)
one will obtain the following reduced equations
2 2
0
0
1 1
' sin 1
2 2 4
f
a a a

| |
= +
|
\ .
(6.7.16)

'
0
cos
2
f
a a

= +
(6.7.17)
The steady state solution can be given by
( ) ( ) cos u a t O = + (6.7.18)

For steady state taking ' a =
'
=0 and eliminating form the equations (6.7.16-17) one can obtain

( )
2
2
2
4
1
4
f

+ =

(6.7.19)
Where
2 2
0
1
4
a = (6.7.20)
Figure 6.7.1 shows the frequency response obtained using equation 6.7.20 for different values of forcing
parameter.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 46 of 104



Fig. 6.7.1: Frequency response curves for primary resonance for van der Pols oscillator

Exercise Problems:
1. Using numerical techniques plot the time response of the system given by Eq. (6.7.1) with primary,
subharmonic and superharmonic resonance. Compare the results with those obtained from the
perturbation method.
2. Find the frequency response curves for van der Pols oscillator considering strong forcing term.
Study the subharmonic and superharmonic resonance conditions.
3. Find the expressions for frequency response curves for a single degree of freedom system for
different resonance conditions when subjected to 2 and 3 frequency excitations.


0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6

Response of the vander pol oscillator for primary resonance

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 47 of 104



Module 6 Lecture 8
Free and Forced vibration of Multi-Degree of freedom system
In this lecture briefly we will discuss about the free vibration of multi-degree of freedom nonlinear
systems. Initially the system with quadratic nonlinearity and then the system with cubic nonlinearity
will be considered.
Free vibration of the system with quadratic nonlinearities
Let us consider a two degree of freedom system where the equation of motion can be given by
2
1 1 1 1 1 1 1 2
0 2 u u u u u = + +
(6.8.1)
2 2
2 2 2 2 2 2 1
2 0 u u u u + + =
(6.8.2)
To use method of multiple scales one may assume the solution of these equations using different time
scales
( )
0 1
,
n
n
T T T t = as follows.
2
1 11 0 1 12 0, 1
( , ) ( ) ... u u T T u T T = + + (6.8.3)
2
2 21 0 1 22 0, 1
( , ) ( ) ... u u T T u T T = + + (6.8.4)
Substituting Eq. (6.8.3) and (6.8.4) in Eq. (6.8.1) and (6.8.2) and separating the terms with different
order of one obtains
Order of
2 2
0 11 1 11
0 D u u + =
(6.8.5)
2 2
0 21 2 21
0 D u u + =
(6.8.6)
Order of
2

2 2
0 12 1 12 0 1 11 1 11 1 11 21
2 ( ) D u u D Du u u u + + + + (6.8.7)
2 2 2
0 22 2 22 0 1 21 2 21 1 11
2 ( ) D u u D Du u u + + + + (6.8.8)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 48 of 104


The solution of Eq. (6.8.5) and (6.8.6) can be given by
11 1 1 1 0
( )exp( ) u A T i T cc = +
(6.8.9)
21 2 1 2 0
( )exp( ) u A T i T cc = +
(6.8.10)
Substituting Eq. (6.8.9) and (6.8.10) in Eq. (6.8.7) and (6.8.8) yields
| | { }
2 2
0 12 1 12 1 1 1 1 1 1 0 1 1 2 1 2 0
Secul
2
ar term
( )exp( ) exp ( ) D u u D A A i T A A i T cc + = + + + +

(6.8.11)
2 2 2
0 22 2 22 2 1 2 2 2 2 0 2 1 1 0 1 1
Secular
2 ( )exp( ) exp(2
em
)
T r
D u u D A A i T A i T A A cc ( + = + + + +

(6.8.12)
To eliminate the secular term from Eq. (6.8.11) and Eq. (6.8.12) one can write
1 1 1 1
0 D A A + =

and
1 2 2 2
0 + = D A A (6.8.13)
Solving Eq. (6.8.13) one can write
1 1 1 1
exp( ) A a T =

and
2 2 2 1
exp( ) A a T =

(6.8.14)
Substituting Eq. (6.8.14) in (6.8.3) and (6.8.4), the first order solution of the system can be given by
( ) ( )
2
1 1 1 1
exp( ) exp u t a i t cc O = + + (

(6.8.15)
( ) ( )
2
2 2 2 2
exp( ) exp u t a i t cc O = + + (

(6.8.16)
Hence, for steady state as time tends to infinity, both the response
1 2
0 u u = =
.
Resonant case (System with internal resonance)
Considering internal resonance of 1:2, i.e., when the second frequency is nearly equal to twice the first
frequency one can write
2 1
2 = +
(6.8.17)
So,
1 0 2 0 0 2 0 1
2 T T T T T = =
and ( )
2 1 0 1 0 0 1 0 1
= + = + T T T T T (6.8.18)
Following similar procedure of method of multiple scales to eliminate the secular terms one obtains the
following equations.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 49 of 104


( ) ( )
'
1 1 1 1 1 2 1 1
2 exp 0 i A A A A i T + + = and
( ) ( )
' 2
2 2 2 2 2 1 1
2 exp 0 i A A A i T + + = (6.8.19)
Now using
( ) ( )
1 1 1 2 2 2
1 1
exp and exp
2 2
A a i A a i = = and introducing
2 1 1
2 T = +
and separating
the real and imaginary parts the following reduced equations are obtained.
' 1
1 1 1 1 2
1
sin
4
a a a a

= + and
' 2 2
2 2 2 1
2
sin
4

= a a a (6.8.20)
' 1
1 1 1 2
1
cos
4
a a a

= and
' 2 2
2 2 1
2
cos
4
a a

= (6.8.21)
Or,
2 2
1 2 2 1
2 2
1 2
cos
2 4
a a
a a



| |
= +
/ |
\ .
(6.8.22)
For steady state using
' ' '
1 2
0 a a = = = , one may write
1
1 1 1 2
1
sin 0
4
a a a

+ = (6.8.23)
2 2
2 2 1
2
sin 0
4
a a

= (6.8.24)
2 2 1 2
2 1 2
1 2
cos 0
2 4
a a a



| |
+ =
|
\ .
(6.8.25)
Eliminating from Eq. (6.8.23) and (6.8.24) one obtains the following equations.
2 2 2 2 1
1 2
1 1 2
0 a a


+ = (6.8.26)
If
1 2
and are of different sign,
1 2
and a a can differ from zero, that is the system may have nontrivial
response. Hence, in the presence of internal resonance, even though there is no external forcing energy
will transfer from 1
st
mode to second mode and self sustained oscillation will continue.
Forced vibration Multi-Degree of freedom system with quadratic nonlinearities
Let us consider a two degree of freedom system with quadratic nonlinearies given by the following
equations.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 50 of 104


( )
1 1
2
1 1 1 1 1 1 1 2
cos 2 f t u u u u u = + + + +
(6.8.27)
( )
2 2 2
2 2 2 2 2 2 1 2 2
2 cos u u u u f t + + + = +
(6.8.28)
Now to solve these equations let us use method of multiple scales by considering
2
1 11 0 1 12 0, 1
( , ) ( ) ... u u T T u T T = + + and

2
2 21 0 1 22 0, 1
( , ) ( ) ... u u T T u T T = + +
(6.8.29)
Substituting Eq. (6.8.29) in Eqs. (6.8.27-6.8.28) and equating coefficients of like power of one
obtains
( )
2 2
0 11 1 11 1 0 1
cos D u u f T + = + (6.8.30)

2 2
0 21 2 21
0 D u u + = (6.8.31)
2 2
0 12 1 12 0 1 11 1 11 1 11 21
2 ( ) D u u D Du u u u + = + (6.8.32)

2 2 2
0 22 2 22 0 1 21 2 21 2 11 2 0 2
2 ( ) cos( ) D u u D Du u u f T + = + + + (6.8.34)

Now, solution of Eq. (6.8.30) and Eq. (6.8.31) can be given by
( )
11 1 1 1 0 0 1
( )exp( ) exp u A T i T i T cc = + + + (

(6.8.35)
21 2 1 2 0
( )exp( ) u A T i T cc = +

(6.8.36) Here,
( )
2 2
1 1
/ 2 f = Now
substituting these two equations in Eq. (6.8.33) and (6.8.34) one can write
( )
2 2 '
0 12 1 12 1 1 1 1 1 0 1 2 1 2 1 0
secular term
2 ( )exp( ) exp D u u A A i T A A i T + = + + (



( ) ( )
2 1 2 1 0 2 2 0 1
1
2
2 1
near secular in case of internal resonance
exp exp A A i T A i T i

| | + + + ( (

|
|
\ .

( ) ( )
1
1 2 2 0 1 1 0 1
1 2
near secular term if near secular term if
exp 2 exp A i T i i i T cc


+
+ + + ( (


(6.8.37)
And

( )
2 2 ' 2
0 22 2 22 2 2 2 2 2 0 2 1 1 0
Secular term
2 ( )exp( ) exp 2 D u u A A i T A i T + = +


( ) ( )
2
2 1 1 1 1 0 1
2 exp A A A i T i + + + + (


NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 51 of 104



( ) ( )
2
1 1 0 1 0 1
2
may be near secular term considering internal resonance
2 exp exp 2 A i T i i T

| | + + + ( (

|
\ .

( )
2
2 0 2
near secular term if
1
exp
2
f i T cc

+ + + (

(6.8.38)
From Eq. (6.8.37) and (6.8.38) one may observe that one may get many resonance conditions such as
(a)
1
, (b)
2
(c)
1 2
+ . Also, some resonancecondition occurs when one consider
internal resonance i.e.
2 1
2 .
Let us consider when the external frequency is nearly equal to the second mode frequency i.e.,
2 1
= + (6.8.39).
Hence, without considering internal resonance, to eliminate secular terms one can write
'
1 1 1
0 A A + =
(6.8.40)
and
( ) ( )
'
2 2 2 2 2 1 1 2
1
2 exp
2
i A A f i T + = + (

(6.8.41)
Now solution of Eq. (6.8.40) can be written as
1 1 1 1 1
1
exp( )
2
A a T i = +
(6.8.42)
( )
( )
2
2 2 2 1 2 1 1 2
2 2 1
1
exp( ) exp
2 4
if
A a T i i T
i


= + + (

+
(6.8.43)


For steady state as t tends to 0, one obtains
( )
( )
2
1 2 1 1 2
2 2 1
0, and exp
4
if
A A i T
i


= = + (

+
(6.8.44)
Hence one obtains,

( )
( ) ( )
2 1
1 1
2 2
1
cos
F
u t O

= + +

(6.8.45)

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 52 of 104


( )
( ) ( )
2 2
2 2 0
2 2
2 2 1
sin
2
F
u t O

= + +
+
(6.8.46)

where
1 1
0
2
tan

| |
=
|
\ .
(6.8.47)

Now if one considers the system with internal resonance i.e.,
2 1
2 or
2 1 2
2 = + , then to
eliminate secular term one can write
( ) ( )
' '
1 1 1 1 1 2 1' 2 1
2 exp 0 i A A A A i T + = (6.8.48)

( ) ( ) ( )
' 2
2 2 2 2 2 1 2 1 2 1 1 2
1
2 exp exp 0
2
i A A A i T f i T ( + + + =

(6.8.49)
Substituting ( )
1
exp
2
n n n
A a i = where n =1, 2, in Eqs. (6.8.48) and (6.8.49) one obtains
' 1
1 1 1 1 2 2
1
sin
4
a a a a

= (6.8.50)
' 2 1 2
2 2 2 1 2 2 2 1
2
1
sin sin
4 2
a a a f

= + + (6.8.51)
' 1
1 1 2 1 2
1
cos
4
a a a

= (6.8.52)
' 2 2
2 2 1 2 2 1
2 2
1
cos cos
4 2
a a f

= (6.8.53)
Where,
1 1 1 2 2
T = + and
2 2 1 2 1
2 T = (6.8.54)
For steady state response,
'
1
a ,
'
2
a ,
'
1

'
2
are zero and there are two possibilities. In the first case one will
obtain the linear solution similar to Eq. 6.8.44. In the second case
1/2
1/2
2 2
1 1 2 2
1
2
4
a P f P
(
| |
=
(
|
\ .
(

(6.8.55)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 53 of 104


( )
1/2
2
2
2 1 1 1 2
2 4 a
(
= +

(6.8.56)
Here,
( )
1 1 2 1 2 1 1 2
2 2 P ( = +


( )
2 1 2 1 1 2 2 1
2 2 P ( =

(6.8.57)
Substituting these values in
1
u and
2
u one can obtain the response of the system. It may be noted that the
second mode amplitude
2
a is not a direct function of the external excitation. Using these equations one
may plot the frequency response curves and observe many resonance phenomena similar to those
observed in case of single degree of freedom system.
Exercise Problems:
1. Plot the frequency response curves using Eq. (6.8.55) and (6.8.56) taking
1
=
2
=-1. Observe
different nonlinear phenomena. Study the stability of the system.
2. Study the resonant and non resonant free vibration of two degree of freedom system with cubic
nonlinearities. The equation of motion of this system can be given by
2 3 2 2 3
1 1 1 1 1 1 1 2 1 2 3 1 2 4 2
2 0 u u u u u u u u u + + + + + + =
2 3 2 2 3
1 1 1 1 1 1 2 1 2 3 1 2 4 2
2 0 u u u u u u u u u + + + + + + =











NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 54 of 104






Module 6 Lect 9
Parametrically excited system
In this lecture, the parametrically excited system will be considered and using Floquet theory the
conditions for instability regions will be determined. It may be noted that a simple parametrically
excited system can be represented by the following equation.

1 2
( ) ( ) 0 u p t u p t u + + = (6.9.1)
Where, the term
1
( ) p t and
2
( ) p t are periodic function of time. It may be noted that as these time varying
terms are coefficients of the response and its derivative, this equation is called the equation of a
parametrically excited system. One can have many variation of this equation by including different
nonlinear terms and forcing terms. One can consider a single or multi degree of freedom system also.
Eq. (6.9.1) can also be written in the following form by substituting
1
1
exp ( )
2
u x p t dt
| |
=
|
\ .

in Eq. (6.9.1). The resulting equation can be written as
( ) 0 x p t x + = (6.9.2)
where
2
2 1 1
1 1
( )
4 2
p t p p p = (6.9.3)
Equation (6.9.2) is called the Hills equation who studied this system in 1886 (Nayfeh and Mook 1979).
Now by substituting ( ) 2 cos2 p t t = + in Eq. (6.9.2) one can write
( 2 cos2 ) 0 x t x + + = (6.9.4)
This equation is known as Mathieus equation. It may be noted that the basics of parametrically excited
systems are based on Mathieu-Hill types of equation. As pointed out in module 4, one may use Floquet
theory to study the stability of these systems.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 55 of 104


Example 6.9.1: Study the stability of the Hills equation given in Eq. (6.9.2) by taking the following
initial condition.
1 1 2 2
(0) 1, (0) 0, (0) 0, (0) 1 u u u u = = = = (6.9.5)
Solution
Writing the fundamental sets of solution as
1 11 1 12 2
2 21 1 22 2
( ) ( ) ( )
( ) ( ) ( )
u t T a u t a u t
u t T a u t a u t
+ = +
+ = +
(6.9.6)
1 11 1 12 2
2 21 1 22 2
( ) ( ) ( )
( ) ( ) ( )
u t T a u t a u t
u t T a u t a u t
+ = +
+ = +


(6.9.7)
one can obtain
11 1 21 2 12 1 22 2
( ), ( ), ( ), ( ) a u T a u T a u T a u T = = = = (6.9.8)
Or
1 1
2 2
( ) ( )
( ) ( )
u T u T
A
u T u T
(
=
(

(6.9.9)
Finding the determinant of A I matrix one may write
2
2 0 + = (6.9.10)
where
1
1 2 1 2 1 2 2
[ ( ) ( )], ( ) ( ) ( ) ( ) u T u T u T u T u T u T = + = (6.9.11)
The parameter is known as the Wronskian determinant of
1
( ) u T and
2
( ) u T .
In case of Hills equation the Wronskian determinant can also be obtained as follows.
As,
1
( ) u t and
2
( ) u t are the fundamental set of solution, hence they must satisfy Eq. (6.9.2). Hence, one can
write
1 1
2 2
( ) 0
( ) 0
u p t u
u p t u
+ =
+ =

(6.9.12)
Or,
1 2 1 2
0 u u u u = (6.9.13)
Which can be integrated to obtain,
1 2 1 2
( ) ( ) ( ) ( ) ( ) t u t u t u t u t = =constant (6.9.14)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 56 of 104


At t =0, ( ) t = 1. So the roots of the Eq. (6.9.10) can be given by
2
1,2
1 = (6.9.15)
Or,
1 2
1 = (6.9.16)
From Eq.(6.9.15) it may be noted that

>1 one of the root will be greater than unity while the other root is
less than one. Hence one of the normal solutions is unbounded and the other is bounded. When

<1 both
the roots will be complex conjugate and their absolute value will be less than one. Hence they will be in the
unit circle. So the solutions will be bounded. It may be observed that the transistion from stable to unstable will
takes place when

=1. This corresponds to a periodic solution of period T when

1 2
1 = = and a periodic
solution of period 2T when
1 2
1 = = .
Example 6.9.2: Study the stability of the Mathieus equation using same initial condition given in
example 6.9.1.
Solution: In case of Mathieu equation
( 2 cos2 ) 0 x t x + + = (6.9.17)
is a function of and . The values of and for which

>1 are called unstable values while those for

which

=1 are called transition values. The locus of transition values separates the plane into regions of
stability and instability as shown in the fallowing figure.

Fig 6.9.1 Transition curve for Matheiu equation
Hills Infinite Determinat
1.9 2 2.1 2.2 2.3 2.4 2.5
0
0.5
1
1.5
2

Stable
Stable
Unstable
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 57 of 104


One may also use hills infinite determinant method to find the transition curve which is explained
below.
( ) 2 cos2 0 u t u + + = (6.9.18)
Using Floquet theory one may assume the solution of the equation (6.9.18) as
( ) ( ) exp u t t = (6.9.19)
Where
( ) ( ) t t T = + . One may expand
( ) t in a Fourier series to obtain the following equation.

( ) exp 2
n
n
u in t

=
( = +

(6.9.20)
Where
n
is constant. Substituting Eq. (6.9.20) in (6.9.18) one obtains
( ) ( )
{ }
( ) ( )
{ }
2
2 exp 2 exp 2 1 exp 2 1 0
n n
n n
in in t t i n t t i n t

= =
(
( ( ( + + + + + + + + =



(6.9.21)
Equating each of the coefficients of the exponential functions to the zero one can obtain the following infinite
set of linear, algebraic, homogenous equations for
n

( ) ( )
2
1 1
2 0
n n n
in
+
(
+ + + + =

(6.9.22)
For a non trivial solution the determinant of the coefficient matrix must vanish. Since the determinant is infinite
one may divide the m
th
row by
2
4m for convergence considerations to obtain the following Hills
determinant.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 58 of 104


( )
( )
( )
( )
( )
2
2 2
2
2 2 2
2
2
2 2 2
2
2 2 2
... . . . . . . . . . ...
4
... 0 0 0 0 0 0 0 ...
4 4
2
... 0 0 0 0 0 0 ...
2 2 2
... 0 0 0 0 0 0 ...
2
... 0 0 0 0 0 0 ...
2 2 2
4
... 0 0 0 0 0 0 ...
4 4 4
i
i
i
i
+

+


+ =

+


+



(6.9.23)
The determinant can be rewritten as (Whittaker and Watson 1962, Nayfeh and Mook 1979)
( ) ( )
|
.
|

\
|
|
.
|

\
|
=

2
1
sin
2
1
sin
0
2
2
i
(6.9.24)
Since the characteristic exponents are solution of
( ) =0, they are given by
( )
2
1
2 1
6
2
1
sin 0 sin
2
(

|
.
|

\
|
=

i
(6.9.25)
Once is known
n
can be related to
0
using equation (6.9.22).
One may also consider the central three rows and columns to approximate the characteristic equation as
follows
( )
( )
( )
2
2
2
2 0
0
0 2
i
i
+
= = +
+ +


Or ( ) ( ) ( ) ( ) ( )
2 2 2 2
2 2 2
2 2 2 2 0 i i i i
( ( ( (
+ + + + + + + =

(6.9.26)
The transition curve separating stability from instability correspond to =0 (i.e., periodic solution with period
) or i = (i.e., periodic motion with period 2).
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 59 of 104


When =0 Eq. (6.9.26) gives the transition curves
2
1
2
= and
2
1
4
2
= + (6.9.27)
When i = Eq. (6.9.26) gives the transition curves
1 = and
2
1
9
8
= + (6.9.28)

Exercise problem
Problem 6.9.1: Use Floquet theory to study the stability of the periodic motion corresponding to
primary resonance of Duffing equation.
Problem 6.9.2: Use Method of Multiple Scales to determine the equations for transition curves for
Matheiu equation. Plot these transition curves near
2
= =1, 4. Taking few points in the stable and
unstable regions and by using numerical method to solve the Mathieus equation, plot the time
responses to check whether the marked instability regions are correct.


Module 6 Lect 10
Multi-degree-of freedom parametrically excited system
Case study: Instability Region of a Sandwich beam
In this lecture a case study has been taken by considering a three-layered, soft-cored, symmetric
sandwich beam subjected to a periodic axial load. For completeness purpose the derivation of the
governing equation of motion is given here and then the parametric instability regions for simple and
combination resonances are investigated for simply supported and clamped-free end conditions by
modified Hsus method.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 60 of 104



Figure 6.10.1 Symmetric three-layered soft cored sandwich beam with periodic axial load
Figure 6.10.1 shows a simply supported, symmetric, three-layered sandwich beam of length L and
width b with a flexible soft core. The top, core and bottom layer thickness are d
t
, c and d
b
,
respectively. The upper and lower layers (face layer) of the beam are of the same elastic material and
the core is of soft viscoelastic material. The sandwich beam is subjected to an axial periodic load
( ) , cos
1 0
t P P t P + = being the frequency of the applied load, t being the time and
0
P and
1
P are
the amplitudes of static and dynamic load, respectively.
Figure 6.10.2 shows the geometry of the sandwich beam, the load and internal forces and moments in
different layers and the deflection in x and z directions before and after deformations. Here,
xx
Q is the
shear force,
xx
N is the axial force and
xx
M is the bending moment. Superscripts and t b represent the
top and bottom layer, respectively. The assumptions made for deriving the governing equations are
similar to that by Frostig [1-4] and are (i) the face sheets of the sandwich beam are modeled as Euler-
Bernoulli beams (ii) the transversely flexible core layer is considered as a two dimensional elastic
medium with small deformations where its height may change under loading, and its cross section does
not remain planar. The longitudinal (in-plane) stresses in the core are neglected and (iii) the interface
layers between the face sheets and the core are assumed to be infinitely rigid and provide perfect
continuity of the deformations at the interfaces.








NA
NA
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 61 of 104
















Figure 6.10.2 (a) Geometry, (b) load, internal forces and moments (c) displacement pattern through depth of
section. N.A is the neutral axis.
The internal potential energy (U ) in terms of direct stresses () and shear stress and strains (, ) is given by,
,
xx xx xx xx c c zz zz
core core top
bot
v v v v
U dv dv dv dv = + + +

(6.10.1)
where,
top
v ,
bot
v , and
core
v are the volume of the top, bottom and core layer, respectively. One may note that, as
the core is taken to be flexible, deformation takes place in the transverse direction ( z direction) and the last term
of equation (1) takes care of that effect. The kinetic energy T can be given by
(b)
(c)
(a)









CORE




NA
NA
Lower Skin
Upper Skin





NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 62 of 104


( ) ( ) ( )
2 2 2 2 2 2
0 0
12
L L
t t t b b b c c c c
core core
v v
T m u w dx m u w dx u dv w dv


= + + + + +
`

)

. (6.10.2)
Here
,
and
t b
m m are the mass per unit length of the top and bottom layer, respectively and
c
is the density of
the core material;
t
u and
b
u are the displacement at the neutral axis of the top and bottom layer along x
(longitudinal) direction, respectively;
t
w and
b
w are the displacement at top and bottom layer along z (vertical)
direction, respectively (Fig. 6.10.2(c)). Also,
c
u and
c
w are the displacement of the core along x and z
directions, respectively and can be given by [1]
( ) { }
, ,
,
( / 2) ( / 2)
( / 2)
b b b x t t t x
c t t t x
u d w u d w z
u u d w
c
+
= + ,
( )
b t
c t
w w z
w w
c

= + (6.10.3)
Here ()
,x
represents the differentiation with respect to x and subscripts t , b, c represent top, bottom and core
layer, respectively.
The non-conservative work done due to the applied load can be given by
( )
2 2
, ,
0 0
12
L L
nc t x b x
W Pw dx Pw dx
| |
= +
|
\ .

. (6.10.4)
The following non-dimensional parameters are used in this analysis.
2
0 0
(2 )
q q
P P L E I = ,
2
1 1
(2 )
q q
P PL E I = ,
2
=
c c c
G A L E ,
2
t t t
E AL E = ,
2
, =
b b b
E A L E
2
c c c
E A L E = ,
( )
2 2
= / ( / )( / ) ( / )( / )
c t t t b b b
g G E c d L d E c d L d + ,
0
/ , t t t = / x x L = , /
q q
u u L = , /
q q
w w L = ,
/
q q
m m m = , /
c c
m m m = . (6.10.5)
Here,
0
P and
1
P are, respectively, the non-dimensional static and dynamic load amplitudes; , and
q q q
E I A are
the Youngs modulus, moment of inertia and the area of cross-section of the qth layer (q equal to t for top layer
and b for bottom layer);
t t b b
E E I E I = + ;
c
E ,
c
A and
c
m are the Youngs modulus, area of cross-section and
mass per unit length of the core, respectively. The non-dimensional time,
4 (1/2)
0
( / ) t mL E = , where mis the total
mass per unit length. The complex shear modulus of the viscoelastic core is given by
*
(1 )
c c c
G G j = + , where
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 63 of 104


c
G is the phase shear modulus, 1 j = and
c
is the core loss factor. The non-dimensional term g is known
as the shear parameter of the system.
Using equations (6.10.1-6.10.5), the governing non-dimensional equations of motion and the boundary
conditions are derived by applying the extended Hamiltons principle. These resulting governing equations of
motion are as follows.
2 2 4
, ,
2 4
, ,
( /3) ( /12)( / ) ( / ) ( /576)( / ){1 ( / )}( / ) ( / )
( / 24)( / )( / )( / ) ( /576)( / ){1 ( / )}( / ) ( / )
( /6) ( /6)( / )( / )
t c t c t t xx c t t c c t xxxx
c t b b xx c b t c c b xxxx
c b c t
m m w m d c c L w m d c d c c L w
m d c d c c L w m d c d c c L w
m w m d c c L u


+ + +
+ + +
+ +



3
, ,
3
, ,
2 2 2
, ,
(1/ 48)( /6){1 ( / )}( / ) ( / )
( /12)( / )( / ) (1/ 48)( /6){1 ( / )}( / ) ( / )
( / ) ( / 4){1 ( / )} ( / ) ( / 4){1 ( / )}{1 ( / )}
t x t c t c c t xxx
c t b x b c t c c b xxx
c t c t t xx c b c t b b xx
m m d c c L u
m d c c L u m m d c c L u
L c w d c w L c w d c d c w



+ +
+ + + +
+ + +
+

3
, ,
3
, ,
2 2 2
, ,
( / 2)( / ){1 ( / )} ( / ){1 ( / )}( / 48)( / )
( / 2)( / ){1 ( / )} ( / ){1 ( / )}( / 48)( / )
( /12)( / ) ( / ) ( / ) 0
c t t x c c t t t xxxxx
c t b x c c t b b xxxxx
t t t xxxx t xx
L c d c u d c c L u
L c d c u d c c L u
d c c L w PL E w

+ + +
+ +
+ + =
(6.10.6)

2 4
, ,
2 2 4
, ,
( /6) ( / 24)( / )( / )( / ) ( /576)( / ){1 ( / )}( / ) ( / )
( /3) ( /12)( / ) ( / ) ( /576)( / ){1 ( / )}( / ) ( / )
( /12)( / )( / )
c t c t b t xx c t b c c t xxxx
b c b c t b xx c b b c c b xxxx
c b
m w m d c d c c L w m d c d c c L w
m m w m d c c L w m d c d c c L w
m d c c L


+ + +
+ + + +

3
, , ,
3 2
, ,
2 2
,
(1/ 48)( /6){1 ( / )}( / ) ( / ) ( / 6)( / )( / )
(1/ 48)( /6){1 ( / )}( / ) ( / ) ( / ) ( / 4){1 ( / )}{1 ( / )}
( / ) ( / 4){1 ( / )}
t x t c b c c t xxx c b b x
b c b c c b xxx c b c t b t xx
c t c b b xx
u m m d c c L u m d c c L u
m m d c c L u L c w d c d c w
L c w d c w



+ +
+ + + + + +
+ +

,
3
, ,
3 2 2 2
, , ,
( / 2)( / ){1 ( / )}
( / ){1 ( / )}( / 48)( / ) ( / 2)( / ){1 ( / )}
( / ){1 ( / )}( / 48)( / ) ( /12)( / ) ( / ) ( / ) 0
c b t x
c c b t t xxxxx c b b x
c c b b b xxxxx b b b xxxx b xx
L c d c u
d c c L u L c d c u
d c c L u d c c L w PL E w

+
+ + +
+ + + =
(6.10.7)
3
, ,
3
, ,
2
,
2
( /6)( / )( / ) ( / 288)( / )( / ) ( / )
( /12)( / )( / ) ( / 288)( / )( / ) ( / )
(1/ 24)( /6)( / ) ( / ) ( /3)
(1/ 24)( /6)( / ) ( / )
c t t x c t c c t xxx
c b b x c b c c b xxx
t c c c t xx t c t
b c c c
m d c c L w m d c c L w
m d c c L w m d c c L w
m m c L u m m u
m m c L

+ + +
+

, ,
2 2
, , ,
2 2
,
( /6) ( / 2){1 ( / )}( / )
( / 2)( / ){1 ( / )} ( / ) ( / )( / 24)( / )
( / ) ( / )( / 24)( / ) 0
b xx c b c t t x
c b b x t t xx c t c c t t xxxx
c b c c b b xxxx
u m u d c L c w
L c d c w u L c u c L u
L c u c L u

+ +
+ + +
+ + =

(6.10.8)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 64 of 104



3
, ,
3
, ,
2
,
( /12)( / )( / ) ( / 288)( / )( / ) ( / )
( /6)( / )( / ) ( / 288)( / )( / ) ( / )
( /6) (1/ 24)( /6)( / ) ( / ) ( /3)
(1/ 24)( /6)( /
c t t x c t c c t xxx
c b b x c b c c b xxx
c t t c c c t xx b c b
b c
m d c c L w m d c c L w
m d c c L w m d c c L w
m u m m c L u m m u
m m c



+
+
+ +
+ +



2
, ,
2 2
, ,
2 2
, ,
) ( / ) ( / 2){1 ( / )}( / )
( / 2)( / ){1 ( / )} ( / ) ( / )( / 24)( / )
( / ) ( / )( / 24)( / ) 0
c c b xx c t t x
c b b x c t c c t t xxxx
b b xx c b c c b b xxxx
L u d c L c w
L c d c w L c u c L u
u L c u c L u



+
+ + +
+ =

(6.10.9)
As the above equations of motion (6.10.6-6.10.9) are in space and time co-ordinates, generalized
Galerkins principle is used to reduce these equations to their temporal form. For multi-mode
discretization one may take
2 3
4
1 1 2 1
3 1
( ) ( ), ( ) ( ), ( ) ( ) and
( ) ( ).
N N N
t p p b q q t r r
N
b s s
p q N r N
s N
w f t w x w f t w x u f t u x
u f t u x
= = + = +
= +
= = =
=

(6.10.10)
Here, N is a positive integer representing the number of modes taken in the analysis, and ( )
p
f t , ( ) t f
q
, ( ) t f
r
and ( ) t f
s
are the generalized co-ordinates and ( )
p
w x , ( ) x w
q
, ( ) x u
r
and ( ) x u
s
are the shape
functions chosen to satisfy as many as the boundary conditions. The resulting equation of motion
becomes

| |{ } | |{ } | |{ } { }
1
cos M f K f P t H f + =

. (6.10.11)
Here, ( ) ( ) t d d = ,
{ } { } { } { } { }
{ }
,
T
T T
T T
p q r s
f f f f f = and | | | | | |
1 0
= K K P H ,
where
| |
| | | | | | | |
| | | | | | | |
| | | | | | | |
| | | | | | | |
11 12 13 14
21 22 23 24
31 32 33 34
41 42 43 44
M M M M
M M M M
M
M M M M
M M M M
(
(
(
=
(
(
(

,
| |
| | | | | | | |
| | | | | | | |
| | | | | | | |
| | | | | | | |
11 12 13 14
21 22 23 24
1
31 32 33 34
41 42 43 44
K K K K
K K K K
K
K K K K
K K K K
(
(
(
=
(
(
(



NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 65 of 104


| |
| | | | | | | |
| | | | | | | |
| | | | | | | |
| | | | | | | |
11
22
H
H
H




(
(
(
=
(
(
(

, { } and | | are null matrices.
The elements of the various sub matrices are given below.
( ) ( ) ( )( )
{ }
( )( )( )( ) ( )
{ }
1 1
2
11
0 0
1
4
0
3 12
576 1 / /
t c i j c t i j
ij
c t t c c i j
M m m ww dx m d L ww dx
m d c d c c L ww dx
| | | |
= + +
| |
\ . \ .
| |
+ +
|
\ .

( ) ( ) ( )( ) { }
( )( )( )( ) ( )
{ }
1 1
2
12
0 0
1
4
0
6 24
576 1
c i j c t b i j
ij
c b t c c i j
M m ww dx m d d l ww dx
m d c d c c L ww dx
| | | |
=
| |
\ . \ .
| |
+ +
|
\ .

( ) ( )( )( )( ) ( )
{ }
( )( ) { }
1 1
3
13
0 0
6 148 1 6
t c t c c i j c t i j
ij
M m m d c c L wu dx m d L wu dx
| | | |
= + + +
| |
\ . \ .


( ) ( )( )( ) ( )
{ }
( )( ) { }
1 1
3
14
0 0
(148) 6 1 /12
b c t c c i j c i j
ij
M m m d c L ww dx m L wu dx
| | | |
= + + +
| |
\ . \ .

( ) ( ) ( )( ) { }
( )( )( )( ) ( )
{ }
1 1
2
21
0 0
1
4
0
6 24
576 1
c i j c t b i j
ij
c t b c c i j
M m ww dx m d d L ww dx
m d c d c c L ww dx
| | | |
=
| |
\ . \ .
| |
+ +
|
\ .

( ) ( ) ( )( )
{ }
( )( )( )( ) ( )
{ }
1 1
2
22
0 0
1
4
0
3 12
576 1 / /
b c i j c i j
ij
c b b c c i j
M m m ww dx m c L ww dx
m d c d c c L ww dx
| | | |
= + +
| |
\ . \ .
| |
+ +
|
\ .

( ) ( )( )( ) ( )
{ }
( )( ) { }
1 1
3
23
0 0
(148) 6 1 12
t c b c c i j c i j
ij
M m m d c c L wu dx m c L wu dx
| | | |
= + +
| |
\ . \ .


NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 66 of 104


( ) ( )( )( ) ( )
{ }
( )( )
1 1
3
24
0 0
(148) 6 1 / /6 /
b c b c c i j c b i j
ij
M m m d c c L wu dx m d L wu dx
| | | |
= + +
| |
\ . \ .

( ) ( )( )( ) ( )
{ }
( )( ) { }
1 1
3
31
0 0
288 + 6 /
c t c c i j c t i j
ij
M m d c c L u w dx m d L u w dx
| | | |
=
| |
\ . \ .


( ) ( )( ) { } ( )( )( ) ( )
{ }
1 1
3
32
0 0
12 288 /
c b i j c b c c i j
ij
M m d L u w dx m d c c L u w dx
| | | |
=
| |
\ . \ .

( ) ( )( ) ( )
{ }
( )
1 1
2
33
0 0
(124) 6 3
t c c c i j t c i j
ij
M m m c L u u dx m m u u dx
| | | |
= + +
| |
\ . \ .

( ) ( ) { } ( )( ) ( )
{ }
1 1
2
34
0 0
6 (124) 6 /
c i j b c c c i j
ij
M m u u dx m m c L u u dx
| | | |
= + +
| |
\ . \ .

( ) ( )( ) { } ( )( )( ) ( )
{ }
1 1
3
41
0 0
12 288
c t i j c t c c i j
ij
M m d L u w dx m d c c L u w dx
| | | |
= +
| |
\ . \ .

( ) ( )( ) { } ( )( )( ) ( )
{ }
1 1
3
42
0 0
6 / 288
c b i j c b c c i j
ij
M m d L u w dx m d c c L u w dx
| | | |
= +
| |
\ . \ .

( ) ( ) { } ( )( ) ( )
{ }
1 1
2
43
0 0
6 (124) 6 /
c i j t c c c i j
ij
M m u u dx m m c L u u dx
| | | |
= + +
| |
\ . \ .

( ) ( ) { } ( )( ) ( )
{ }
1 1
2
44
0 0
3 (124) 6 /
b c i j b c c c i j
ij
M m m u u dx m m c L u u dx
| | | |
= + +
| |
\ . \ .

( ) ( )( )
{ }
( )( )
{ }
( )
{ }
1 1 1
2 2 2
11
0 0 0
14 1 112 /
t c i j t c i j c i j
ij
K d c ww dx d L ww dx L c ww dx
| | | | | |
= + + +
| | |
\ . \ . \ .

( ) ( )
{ }
( )( )( ) { }
( )( ) ( )( )
{ }
1 1
2
12
0 0
1
3
0
14 1 1
148 1
c i j t b c i j
ij
t c c b i j
K L c ww dx d c d c ww dx
c L d c wu dx


| | | |
= + + +
| |
\ . \ .
| |
+
|
\ .

( ) ( )( )( ) { } ( )( ) ( )( )
{ }
1 1
3
13
0 0
12 / 1 148 1
t c i j t c c t i j
ij
K L c d c wu dx c L d c wu dx
| | | |
= + +
| |
\ . \ .

( ) ( )( )( ) { } ( )( ) ( )( )
{ }
1 1
3
14
0 0
12 / 1 148 1
t c i j t c c b i j
ij
K L c d c wu dx c L d c wu dx
| | | |
= + + +
| |
\ . \ .

( ) ( )
{ }
( )( )( ) { }
1 1
2
21
0 0
14 1 1
c i j t b c i j
ij
K L c ww dx d c d c ww dx
| | | |
= + + +
| |
\ . \ .


NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 67 of 104


( ) ( )( )
{ }
( )( )
{ }
( )
{ }
1 1 1
2 2 2
22
0 0 0
14 1 112
b c i j b c i j c i j
ij
K d c ww dx d L ww dx L c ww dx
| | | | | |
= + + +
| | |
\ . \ . \ .

( ) ( )( )( ) { } ( )( ) ( )( )
{ }
1 1
3
23
0 0
12 1 148 / 1
b c i j b c c t i j
ij
K L c d c wu dx c L d c wu dx
| | | |
= + +
| |
\ . \ .

( ) ( )( )( ) { } ( )( ) ( )( )
{ }
1 1
3
24
0 0
12 / 1 148 1
b c i j b c c b i j
ij
K L c d c wu dx c L d c wu dx
| | | |
= + + +
| |
\ . \ .

( ) ( )( )( ) { }
1
31
0
12 1
t c i j
ij
K L c d c u w dx
| |
= +
|
\ .

( ) ( )( )( ) { }
1
32
0
12 1
b c i j
ij
K L c d c u w dx
| |
= +
|
\ .

( ) ( ) ( )
{ }
( )( ) ( )
{ }
1 1 1
2 2
33
0 0 0
/ 124
t i j c i j c c t i j
ij
K u u dx L c u u dx c L u u dx
| | | | | |
=
| | |
\ . \ . \ .


( ) ( )
{ }
( )( ) ( )
{ }
1 1
2 2
34
0 0
124
c i j c c b i j
ij
K L c u u dx c L u u dx
| | | |
= +
| |
\ . \ .


( ) ( )( )( ) { }
1
41
0
12 / 1
t c i j
ij
K L c d c u w dx
| |
= +
|
\ .

( ) ( )( )( ) { }
1
42
0
12 1
b c i j
ij
K L c d c u w dx
| |
= +
|
\ .

( ) ( )
{ }
( )( ) ( )
{ }
1 1
2 2
43
0 0
/ 124
c i j c c t i j
ij
K L c u u dx c L u u dx
| | | |
= +
| |
\ . \ .


( ) ( ) ( )
{ }
( )( ) ( )
{ }
1 1 1
2 2
44
0 0 0
/ 124
b i j c i j c c b i j
ij
K u u dx L c u u dx c L u u dx
| | | | | |
=
| | |
\ . \ . \ .


( )
1
11
0
i j
F ww dx =

( )
1
22
0
i j
F ww dx =

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 68 of 104


In the above equations () () x = .
Equation (6.10.11) is a set of coupled Mathieu-Hill equations with complex coefficients. For numerical
calculations following shape functions (Ray and Kar [7]) are considered.
For simply supported beam
( ) ( ) sin
p
w x p x = , ( ) ( ) x q x w
q
sin = , ( ) ( ) x r x u
r
cos = , and ( ) ( ) x s x u
s
cos = . (6.10.12)
These shape functions satisfy all the boundary conditions. Here 1,2 , p N = ,
( ), q q N = ( ) 2 r r N = and ( ) 3 s s N = . (6. 10.13a)
For clamped-free beam, the shape functions are as follows (Ray and Kar [7])
{ } { } { }
{ } { }
( 1)
2 2 1 2
2 ( 3)
2 1 2
( 1)
( ) ( 3)( 2) ( 2)( 1) 2( 3)( 1) ( 2) ( 1) ( 2)( 1)
( 2)( 1) ( 1) ( 1) / ( 3)( 2)( 1) ( 3) ,
( ) ( 1)
i
i
i
k k
k
w x i i i i x i i i i i i i i
i i i i i i i i i i x
u x k x kx
+
+
+
( = + + + + + + + + + + +

+ + + + + + + + + +

= +
(6.10.13b)
Here and i k are same as the previous boundary conditions.
If
| |
L is a normalized modal matrix of | | | |,
1
K M

then the linear transformation
{ } | |{ } f L U = , (6.10.14)
transforms equation (6.10.11) to,

4
2
1
( ) 2 cos 0
N
q q q qp p
p
U U t b U

=
+ + =

, = q 1.4N; (6.10.15)
where ( )
2

q
are the distinct eigen values of | | | | K M
1
and

qp
b are the elements of
| | | | | | | || |
1 1
B L M H L

= . Also, < = 2
1
P 1 for the present analysis. The complex frequency and forcing
parameters in terms of real and imaginary parts are given by
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 69 of 104



, , q q R q I
j

= + and
, , qp qp R qp I
b b jb

= + . (6.10.16)
The boundaries of the regions of instability for simple and combination resonances are obtained by the modified
Hsus [9] method. When the system is excited at a frequency nearly equal to twice the natural frequencies
principal parametric resonance and when it is excited near a frequency, which is equal to the sum or differences
of any two modal frequencies combination resonances of sum or difference types take place. Following relations
are used to obtain the boundaries of the regions of instability for simple and combination resonances [7].
(1) Simple resonance case

,
( / 2)
R
<
1
4

, =1, 2..4N (6.10.17)
where

=
2 2 2
, , 2
,
2
,
4 (
16
R I
I
R
b b

(
+
(
(

. (6.10.18)
(2) Combination resonance of sum type

, ,
( )
R R
+ <

(6.10.19)
when damping is present,

=
1/2
2
, , , , , ,
, ,
1/2
, ,
, ,
( ) 4 ( )
16 ,
4( )
I I R R I I
I I
R R
I I
b b b b








(
+ +
(
(

(6.10.20)
and for the undamped case,

1/2
, ,
, ,

.

R R
R R
b b

(
=
(
(

, , = 1, 2..4N (6.10.21)
(3) Combination resonance of difference type
, ,
( )
R R
< > , , = 1, 2..4N (6.10.22)
when damping is present,
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 70 of 104


1/2
2
, , , , , ,
, ,
1/2
, ,
, ,
( ) 4 ( )
16
4( )
I I I I R R
I I
R R
I I
b b b b








(
+
( =
(

, (6.10.23)
and for the undamped case,
1/2
, ,
, ,

.

R R
R R
b b

(
=
(
(

(6.10.24)
Numerical Results and Discussions
Here the parametric instability regions of a three-layered symmetric sandwich beam with simply
supported, and clamped-free boundary conditions have been determined numerically using MATLAB.
For visco-elastic materials, core loss factor (
c
) is a measure energy dissipation capacity and the shear
parameter
( )
2
= / 2 ( / )( / )
c t t t
g G E c d L d is a measure of stiffness of the material and is important in
determining how much energy gets into the visco-elastic material. So these two parameters are varied
in determining the instability regions for the parametrically excited beams. Also the effects of core and
skin thickness on the instability regions are studied for all these boundary conditions. In the parametric
instability regions shown in the following figures, the regions enclosed by the curves are unstable and
the regions outside the curves are stable. Here the ordinate
1
P is the amplitude of non-dimensional
dynamic load and the abscissa is the non-dimensional forcing frequency. Following physical
parameters are taken for the numerical analysis. The span of the beam, L=300 cm, width, b=50 mm, the
top and bottom face thickness 2 = =
b t
d d mm and the core thickness, c = 30 mm. The non-
dimensional static load amplitude
0
0.1 P = for all the figures except it is specifically mentioned. The
top and bottom faces are of steel and the core is of soft plastic foam (Divinycell H60). The mechanical
properties of steel and Divinycell H60 are given in Table (1).
Table 1:Material properties of sandwich beam [10]

Material
Youngs
modulus
E, Gpa
Shear
modulus
G, Gpa
Poissons
ratio

Density
, kg/m
3


Steel

210

81

0.3

7900
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 71 of 104




Divinycell H60


0.056

0.022

0.27

60

Simply supported beam
Using the shape functions given in equation (6.10. 12) the instability regions for the simply supported
beam are determined and shown in Figures (6.10.3-6.11.6) for the first three modes. Figure 6.10.3
shows the parametric instability regions obtained using both the higher-order theory and classical
theory [10] for simple resonances. One may observe that for all the three modes, the region of
instability starts at a lower frequency for higher order theory in comparison to the classical theory,
which is due to the fact that, the core is considered to be more flexible in higher order theory than in
case of classical theory. Also, it is clearly observed from these figures that the instability region is
wider in case of higher order theory as compared to the classical theory. With change in
0
P (say
0
P
=0.1), while instability region with higher order theory remains almost unchanged, it is observed that
for lower value of
1
P , the instability regions with classical theory shifts towards left.

Figure 6.10.3: Comparison of instability regions using higher-order and classical theories,
0
0.5,
c
P = = 0.1;
g=0.05; , higher-order theory; ------, classical theory.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 72 of 104



Figure 6.10.4: Effect of shear parameter on instability regions for =
c
0.0. , g =0.05; -----, g =0.1; +++
, g =0.5.

Figure 6.10.5: Effect of shear parameter on instability regions for =
c
0.3. , g =0.05; -----, g =0.1; +++,
g =0.5.
Figures (6.10.4-5) show the influence of core loss factor (
c
) and the shear parameter (g) upon the instability
region obtained by using higher order theory. It is clearly observed that increase in core loss factor improves the
stability by shifting the instability zones upwards and reducing the area of instability, which is similar to those,
obtained by classical theory. It is also observed that with increase in shear parameter stability of the system
improves. From the above figures it is clearly understood that to get a more stable system one may go for higher
value of core loss factor (
c
) and shear parameter (g).
Clamped-Free beam.
Using the shape functions (equation (6.10.15)) for the clamped-free beam, the instability regions for the first
three modes are determined and shown in Figures (6.10.6-8).
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 73 of 104



Figure 6.10.6: Comparison of instability regions using higher order and classical theories for
0
0.5,
c
P = = 0.1,
g =0.05. , higher-order theory; ------, classical theory.
Using higher order theory and classical theory the parametric instability regions for simple resonances
are shown in Figure 6.10.6. Here also, higher order theory gives a conservative design for lower modes.

Figures (6.10.7-8) show the influence of core loss factor (
c
) and the shear parameter (g) upon the instability
region and it is observed that with increase in core loss factor and shear parameter stability of the system
improves.

Figure 6-10.7: Effect of shear parameter on instability regions for =
c
0. , g =0.05;
-----, g =0.1; +++, g =0.5.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 74 of 104



Figure 6.10.8: Effect of shear parameter on instability regions for =
c
0.3. , g =0.05; -----, g =0.1; +++, g =
0.5.
For all the boundary conditions the system is always found to be stable at combination resonances of sum and
difference type. In these cases, for simple resonances it is observed that with increase in shear parameter the
instability plot moves upward implying that there exists critical forcing amplitude below which the system is
always stable. For example, when a cantilevered sandwich beam with
c
=0.3 and g=0.1 is excited near twice the
first natural frequency ( 8.2 ), the system will not vibrate if the forcing amplitude is less than 0.485 (point P
c

on figure 6.10.8). But for the same
c
and g =0.05, with same amplitude of forcing, the system will vibrate at a
slightly less frequency (say 7.8 ). Again with increase in shear parameter, the instability region shifts
towards right and hence, for same forcing amplitude, the system becomes unstable at a higher frequency. As the
shear parameter
( )
2
= / 2 ( / )( / )
c t t t
g G E c d L d , is a function of dimension and material properties of both skin
and core material, using the above stability charts, a designer will be able to construct sandwich beams having
very less or vibration free structures.
REFERENCES
1. FROSTIG, Y., and BARUCH, M., Bending of sandwich beams with transversely flexible core.
American Institute of Aeronautics and Astronautics Journal, 27, 523-531, 1990.
2. FROSTIG, Y., and BARUCH, M., Free vibrations of sandwich beams with a transversely flexible core:
A higher order approach. Journal of Sound and Vibration, 176(2), 195-208, 1994.
3. FROSTIG, Y., Buckling of sandwich panels with flexible core-high order theory. International Journal
of Solids and Structures, 35(3-4), 183-204,1998.
P
c
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 75 of 104


4. FROSTIG, Y., and THOMSEN, O.T., High-order free vibration of sandwich panels with a flexible core.
International Journal of Solids and Structures, 41, 1697-1724, 2004.
5. SAITO, H., and OTOMI, K., Parametric response of viscoelastically supported beams. Journal of Sound
and Vibration, 63, 169-178, 1979.
6. KAR, R.C., and SUJ ATA, T., Dynamic stability of a tapered symmetric sandwich beam. Computers &
Structures, 40, 1441-1449, 1991.
7. RAY, K., and KAR, R.C., Parametric instability of a sandwich beam with various boundary conditions.
Computers & Structures, 55, 857-870, 1995.
8. RAY, K., and KAR, R.C., Parametric instability of multi-layered sandwich beams. Journal of Sound
and Vibration, 193(3), 631-644, 1996.
9. HSU, C.S., On the parametric excitation of a dynamic system having multiple degrees of freedom.
Journal of Applied Mechanics, ASME, 30, 367-372, 1963.
10. Dwivedy S. K., Sahu K. C. and Babu Sk., Parametric instability regions of three layered soft-cored
sandwich beam using higher order theory, Journal of Sound and Vibration, 304, 326-344, 2007.






Module 6 Lecture 11
Parametrically excited continuous system
Case study: Nonlinear Vibration of a Magneto-Elastic Cantilever Beam With Tip Mass
In this work the effect of the application of alternating magnetic field on the large transverse vibration
of a cantilever beam with tip mass is investigated. The governing equation of motion is derived using
the D Alemberts principle which is reduced to its non-dimensional temporal form by using the
generalize Galerkins method. The temporal equation of motion of the system contains the
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 76 of 104


nonlinearities of geometric and inertial type along with parametric excitation and non-linear damping
term. Method of multiple scales is used to determine the instability region and frequency response
curves of the system. The influences of the damping, tip mass, amplitude of magnetic field strength,
permeability and conductivity of the beam material on the frequency response curves are investigated.
These perturbation results are found to be in good agreement with those obtained by numerically
solving the temporal equation of motion and experimental results.








Fig. 6.11.1: Schematic diagram of a flexible single-link cantilever beam with tip mass.
Figure 6.11.1 shows a flexible cantilever beam with a tip mass M. The beam is subjected to a harmonic
transverse magnetic field
0
cos
m
B B t = where
m
B and are respectively, the amplitude and
frequency of the magnetic field strength. In this work, the flexible cantilever beam with tip mass is
modeled as an Euler-Bernoulli beam with a tip mass. For the purpose of completeness a brief derivation
s


v
u
M
X
Y
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 77 of 104


of the equation of motion using d Alemberts principle is given below. The bending moment ( ) M s of
the beam at a distance s from the fixed end (Fig. 1) can be expressed as [7,11]
( )
2
1
2
ss s ss
M s E I v v v
| |
+
|
\ .
. (6.11.1)
Here, v is the transverse displacement of the beam. ( )
s
is the first derivative with respect to s. One
may write the inextensibility condition of the beam in terms of longitudinal displacement ( ) t u , and
transverse displacement ( ) t v , as [7]
( )
2
2
1 1
s s
v u + + = . or,
( ) ( )
1
2
2
, 1
0
u t v d

. (6.11.2)
Here , are the integration variables. Considering the inertia forces , , and Au Av M u M v , and
using the d Alemberts principle, one may write Eq. (6.11.1) as follows
( ) ( ) ( ) 0
L
M s M s M s

= . (6.11.3)
Here ( ) M s

is the summation of the moment due to inertia force of beam and the moment due to
external magnetic force a distance from the roller support and the couple due to magnetic field.
( )
L
M s is the moment due to inertia force for the pay load at the tip of the manipulator. The
expressions for these moments are given below.

( ) ( ) sin cos
sin cos ,
d
L L
M s Au d d Av C d d
s s s s
L
p d d c d
s s s
v

= +

(6.11.4)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 78 of 104


Here, p and c are the body force and body couple of the beam due to the magnetic field
0
B which are
expressed as [3, 5, 6, 9]
2 2 2 2
0 0
0
1 1
1 , and
2 2
0
s s s s s s
r
m
p hdB v v v v v v d c hdB v
s
| | | |
= =

| |
\ . \ .




. (6.11.5)
Also,
( ) sin cos
L L
M s M u d M v d
L
s s
=

. (6.11.6)
By differentiating Eq. (6.11.3) twice with respect to s and applying the Leibnitzs rules one may obtain
the following governing differential equation of motion.
( ) ( )
2 3
2

1
3
2
s
ssss s ssss s ss sss ss s s ss
0
v v v d

EI v v v v v v v Av v v
| |
+ |
\ .
+ + +

+ +
( )
( )
( ) ( )
2 2
L
L
d s ss ss
s
s 0 0
d Av C v M v v v d d M d v A v v v v v v

+ + +

| |
+ +
|
\ .


( ) ( )
( ) ( ) ( )
2 2 2
1
1 1
1 1 1 0
2 2
2
L
ss s s s ss s s d
s
dc
v Av C v pd pv v v v v c
ds
v + + + =

| | | |
| |
\ . \ .

. (6.11.7)
To obtain the temporal equation of motion, one may discretize the governing equation of motion
(6.11.7) by using following assumed mode expression.
( ) ( ) ( )
, v s t r s q t = . (6.11.8)
Here, r is the scaling factor; q(t) is the time modulation and ( ) s is the eigen-function of the cantilever
beam with tip mass, which is given by [7]
( ) ( )
sin sinh
( ) cos cosh sin sinh
cos cosh
L L
s s s s s
L L
| | +
= +
|
+
\ .
. (6.11.9)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 79 of 104


One may determine L from the following equation.
( ) 1 cos cosh cos sinh sin cosh 0 L L m L L L L L + + = . (6.11.10)
Following non-dimensional parameters are used in this analysis.
L
s
x = ,
e
t = ,
e

,
r
r
L
= ,
M
m
AL
=

,
4
AL
EI

= . (6.11.11)
Substituting Eq. (6.11.9) into Eq. (6.11.6) and using the generalized Galerkins method, one may obtain
the resulting non-dimensional temporal equation of motion, which can be expressed as
( ) ( )
( ) ( )
3 2 2
1 2 3 1
2
1
2 cos 2
- 1 cos 2 0.
q q q q q q q q f q
k qq
+ + + + +
+ =

. (6.11.12)
The expressions for the coefficients (
1 2 3 1 1
i.e. , , , , , f k ) in this equation are given below.
The natural frequency of the lateral vibration of an elastic beam

2
15 1 1
4 2
14 0 14
2
2
m m
r
B
e
hd h E I h
mL h mL h
| |
=
|

\ .
,

( )
2
2
15 1 1
4
14 0 1 1
2
2
1 1
m m
m
r
B
L
hd L h E I h
B
mL h E I h
| | | |
= =
| |
|

\ . \ .
. (A1)
Here,
1 1
4
14
2
L
E I h
mL h
= ,
2
15
0 1 1
2
and .
2
m
r
m
m
B hd L
h
B
E I h
| |
=
|

\ .

Damping ratio due to the viscous damping to the system,
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 80 of 104


2
d L
e e
C
m
| |
= = = |
|

\ .
, (A2)
Coefficient of the nonlinear geometric term
3
q =
3 2 4
1
14 14 14
2
4 2
3
2
e
h h h EI r
h h h
mL
| |
= + +
|

\ .
, (A3)
Coefficient of the nonlinear inertia term
2
q q =
5 6 7 8 9 10
2
14 14 14 14 14 14
2
h h h h h h r
m m
h h h h h h
| |
= + +
|

\ .
, (A4)
Coefficient of the nonlinear inertia term
2
q q =
13 11 12
3
14 14 14
2
h h h r
m
h h h
| |
=
|

\ .
, (A5)
Coefficient of the parametric excitation termcos(2 )q =
0
1 2
2
2 2
r
e c
B
f
f
B
| |
= =
|

\ .
, where
2
2
2
m
r
B
B = and
2
0 1 1
14
2 r
m
c
E IL h
B
hd h
| |
=
|

\ .
. (A6)
Coefficient of the nonlinear damping terms ( )
2
1 cos(2 ) q q +
=
2
2 16 17
1
14 14
2
m
e
B hd h h
k r
m h h
| |
= +
|

\ .
. (A7)
Here,
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 81 of 104


( )
( )
4 1
1
0
d x
h x dx
dx

,
( )
( )
3
2 1
2
0
d x
h x dx
dx
| |
=
|
\ .

,
( ) ( )
( )
4 1
3
0
d x d x
h x dx
dx dx

=

( ) ( ) ( )
( )
2 3 1
4 2 3
0
d x d x d x
h x dx
dx dx dx

=

,
( ) ( )
2
1
5
0 0
( )
x
d
d x
h d x dx
dx d
| |

| =
|

\ .

,
( ) ( )
( ) ( )
2 1 1
6 2
0 x
d x d x
h d x dx
dx dx

=

,
( ) ( )
( )
2 1
2
7 2
0
( )
d x d x
h x dx
dx dx

=

,
( ) ( )
( )
2
2 1 1
8 2
0 0 x
d x d x
h d d x dx
dx d

| |
=
|

\ .

,
( ) ( )
( )
2
2 1
9 2
0 0
x
d
d x
h d x dx
dx d
| |

| =
|

\ .

,
( )
( ) ( )
2
1
2
10
0
d x
h x dx
dx
| |
=
|
\ .

,
( ) ( )
( )
2
1
11
0 0
x
d
d x
h d x dx
dx d
| |

| =
|

\ .

,
( ) ( )
( )
2
2 1 1
12 2
0 0 x
d
d x
h d d x dx
dx d

| |

| =
|

\ .

,
( ) ( )
( )
2
2 1
13 2
0 0
x
d
d x
h d x dx
dx d
| |

| =
|

\ .

,
( ) ( )
1
2
14
0
h x dx =

,
( )
( )
2 1
15 2
0
d x
h x dx
dx
(
=
(

,
1
16
0 0
( ) ( )
( )
x
d x d
h d x dx
dx d
| |

=
|

\ .

,
and
1 1 2
17 2
0 0
( ) ( )
( )
x
d x d
h d d x dx
dx d

| |

=
|
|

\ .

.
Here one may observe that the non-dimensional temporal Eq. (6.11.12) has parametric term
( )
1
cos 2 f q and nonlinear damping term
( ) ( )
2
1
1 cos 2 k qq + , along with cubic geometric (
1
3
q )
and inertial (
2 3
2 2
q q q q + ) nonlinear terms. Hence, it may be noted that the temporal equation of
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 82 of 104


motion Eq. (6.11.12) contains many nonlinear terms and it is very difficult to find the exact solution.
Hence one may go for the approximate solution by using the perturbation method. Here method of
multiple scales is used.
In method of multiples scales, the displacement q can be represented in terms of different time scales
0 1
( , ) T T and a book keeping parameter as follows.

( ) ( ) ( )
2
0 0 1 1 0 1
( ; ) , , q q T T q T T O = + + . (6.11.13)
Here,
0
T = , and
1
T = . The transformation of first and second time derivatives are given by
) (
2
1 0

O D D
d
d
+ + = , (6.11.14)
) ( 2
2
1 0
2
0
2
2

O D D D
d
d
+ + = . (6.11.15)
where,
0
0
D
T

, and
1
1
D
T

. Substituting Eqs. (6.11.13- 6.11.15) into Eq. (6.11.12) and equating the
coefficient of like powers of , yields the following equations.
Order :
0

2
0 0 0
0 D q q + = , (6.11.16 )
Order :
1

2
0 1 1 0 1 0 0 0
2 2 D q q D D q D q + =
( ) ( )
2
3 2 2
1 0 2 0 0 0 3 0 0 0
q D q q D q q
( ) ( ) ( )( )
2
1 0 0 1 0 0 0 0
cos 2 1 cos 2 f T q k T D q q + + + . (6.11.17)
General solutions of Eq. (6.11.16) can be written as
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 83 of 104



0 1 2 0 2 0
( , )exp( ) ( , )exp( ) q A T T iT A T T iT = + . (6.11.18)
Substituting Eq. (6.11.18) into Eq. (6.11.17) leads to
( )
2 2
0 1 1 0 0 1 2 3 1 0
2 exp( ) 2 exp( ) 3 3 exp( ) D q q i A iT i A iT ik A A iT + = +
( ) ( ) ( )
3 3 1
1 2 3 1 0
+ exp 3 + exp 3 + exp (2 1) exp (2 1)
2
0 0
f
A iT ik A iT A i A i T ( + + +



3 2 2 1 1 1
0 0 0
+ exp (2 3) exp (2 1) exp (2 1)
2 2 2
ik ik ik
A i T A A i T A A i T + + +

3 1
0
+ exp (3 2 )
2
ik
A i T cc + . (6.11.19)
Here, cc is the complex conjugate of the preceding terms One may observe that any solution of Eq.
(6.11.19) will contain secular or small divisor terms when non-dimensional frequency of magnetic field
strength ( ) is nearly equal to 1 which may be called as simple resonance case. In this case, one may
use detuning parameter to express the nearness of to 1, as
( ) 1 O and , 1 = + = . (6.11.20)
Substituting Eq. (6.11.20) into Eq. (6.11.19), one may obtain the following secular or small divisor
terms.
( )
2
0 0 1 2 3 1 0
2 exp( ) 2 exp( ) 3 3 exp( ) i A iT i A iT ik A A iT +

( ) ( ) ( )
3 2 1 1 1
1 1 1
exp 2 exp 2 exp 2 0.
2 2 2
f k k
A T i A T i A A T + + = (21)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 84 of 104


Putting A equal to
( )
( )
1
1
1

2
i T
a T e

and
1
2 2 T = into Eq. (6.11.21) and separating the real and
imaginary terms, one may a set of reduced equations as given below.

3
1 1
sin
8 4
k f
a a a a = + + , (6.11.22)

3 3
3 1
1 2 1
1 3 1
2 sin cos
4 3 4 2
f
a a a a k a
| | | |
= + + +
| |

\ . \ .
. (6.11.23)
One may observe from the Eqs. (6.11.21)-( 6.11.22) that the system possesses both trivial and
nontrivial responses. Hence one may determine both responses by solving Eqs. (6.11.22, 6.11.23)
simultaneously. To find the stability of the steady state responses, one may perturb the above Eqs.
(6.11.22, 6.11.23), by substituting
1 o
a a a = + and
0 1
= + where
0 0
, a are the singular points, and
then investigate the eigenvalues of the J acobian matrix (J) which is given by

2
1 1 1
0 0 0 0
3
1 2 0 0 1 0 2
1
0 1 0 0
3
sin cos
8 4 4
3 1
+ sin
1
cos sin 2 3 2
4 2

k f f
a a
J
a a k
f
a k
(
+ +
(
(
= | |
(
+
|
(

\ .
(
(

. (6.11.24)
It may be noted that the system will be stable if and only if all the real parts of the eigen-values are
negative.
For trivial state instability region, one may use the following expression for the transition curve for
simple resonance case which has been obtained by finding the eigen-values of J acobian matrix (J)
given in Eq.( 6.11.24).
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 85 of 104



( )
2
2 2 0
4 .
2 4
L
f
O
| |
= = +
|

\ .
(6.11.25)
Here, the expression for
0
, , and f are given in appendix [7].
It may be noted that this simple expression has been obtained by using first order method of multiple
scales is different from the expression given in the work of Pratiher and Dwivedy [7], which was
obtained by using the second order method of multiple scales.
Now the first order non-trivial steady state response of the cantilever beam with a tip mass can be given
by
( )
1
cos
2
q a
| |
=
|
\ .
. (6.11.26)
Here for numerical simulations, a steel beam similar to that considered in the work of Wu [5, 6] with
length L =0.5 m, width d =0.005 m, depth h =0.001 m, Youngs Modulus E =
11
10 94 1 . N/m
2
, mass
of the beam per unit length m =0.03965 kg, and the permeability of the vacuum,
0
=
6
10 26 . 1

Hm
-1

have been considered. Using these parameters, the reduced Eqs. (6.11.22, 6.10.23), have been solved
numerically to obtain the instability regions and the frequency response curves. In the instability plot,
the regions bounded by the curves are unstable and regions outside the curves are stable. In the
frequency response curves dotted and solid lines represent, respectively the unstable and stable
response of the system. The effect of the amplitude of magnetic field strength (
m
B ), damping (
d
C ), tip
mass (M), material conductivity (), and relative permeability of the material (
r
) on the frequency
response have been investigated.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 86 of 104



Fig. 6.11.2: The region of instability of a cantilever beam with tip mass subjected to magnetic field,
-.-.-.-.-. Moon and Paos theoretical result, - - - - - Moon and Paos experimental result, and ---- the
present result.
For simple resonance case, the beam is subjected to a transverse magnetic field with a frequency nearly
equal to the natural frequency of the system. Here, the instability regions are plotted in ( )
2
/
L
Vs
( )
2
/
r c
B B plane similar to the work of Moon and Pao [2], Wu [5, 6], and Pratiher and Dwivedy [7]. The
experimental and theoretical results of Moon and Pao [2] are also being plotted in Fig.6.11.2 for
comparison with the present result. It is found from Fig.6.11.2 that the result obtained in the present
work is in good agreement with the experimental results Moon and Pao [2]. The accuracy of the
instability region obtained by using the first order method of multiple scales can be verified by
numerically solving the temporal equation of motion (6.11.12) and plotting the time response (Fig.
6.11.3) for two different points A and B as marked in Fig.6.11.2. Figure 6.11.3(i) clearly shows that the
response is stable and Fig.6.11. 3(ii) shows that the response is unstable which are in good agreement
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 87 of 104


with the result shown in Fig.6.11. 2. Hence, one may use the first order closed form solution (Eq.
6.11.25) for finding the instability region instead of going for a second order solution as reported in the
work of Pratiher and Dwivedy [7]. But when more accurate result is required, one may use the
expression given in the work of Pratiher and Dwivedy [7].


Fig.6.113. (a) Time response for the point A and (b) time response for point B marked in Fig.6.11.2.
Figure 6.11.4 shows the frequency response curve for four different values of amplitude of magnetic
field strength
m
B . From Fig.6.11.4, it may be noted that with increase in
m
B , though the maximum
response amplitude remain unchanged, the trivial state becomes unstable which is similar to that shown
in Fig.6.11.2. The trivial state becomes unstable by the sub-critical pitchfork bifurcation at
1
R , which
ends with a super-critical pitchfork bifurcation at
2
R . Here, one may observe that the system has a
tendency to jump up from the unstable trivial state at
1
R to the stable non-trivial state at
/
1
R .
Figures 6.11.5(a) and (b) show the transient and steady state response for point C marked in Fig.
6.11.4(c). The solid line and dotted line respectively represent the response of the system with and
without magnetic field. In the presence of magnetic field, it clearly shows that the steady state response
q
q
Time () Time ()
(a) (b)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 88 of 104


has zero response amplitude. Also it may be noted that the free vibration response of the beam shown
as dotted line in Fig. 6.11.5, is reduced by applying the magnetic field.



Fig.6.11.4. Effect of the magnetic field strength (
m
B ) on the frequency response curves for 0.02 kg M = ,
2
0.01 N-s/m
d
C = , 3000
r
= ,
7
10 = Vm
-1
(a) 0.20
m
B = Am
-1
(b) 0.25
m
B = Am
-1
(c) 0.30
m
B = Am
-1

(d) 0.35
m
B = Am
-1
.
(a) (b)
a a

(c)
(d)
a
a


A BA


C
R
2
R
1


NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 89 of 104



Fig.6. 11. 5. (a) Transient response and (b) steady state time response for the point C with and without magnetic
field.
The effect of damping
d
C on the response curves is shown in Fig.6.11.6 and it has been observed that
with increase in
d
C , while the non-trivial response amplitude remains unchanged the trivial state
unstable region decreases, the sub-critical pitchfork bifurcation point occurs at a higher value of and
the corresponding jump length decreases.

Fig.6.11.6: Influence of damping on frequency response curve for 0.02 kg M = , 3000
r
= ,
7
10 = Vm
-1
,
0.30
m
B = Am
-1
(a)
2
0.02 N-s/m
d
C = ,(b)
2
0.03 N-s/m
d
C = .
Similarly one can study the influence of effect of relative permeability (
r
), material conductivity ()
and mass ratio on the frequency response curves of the system.
(a)
(b)
q q
Time () Time ()
(a) (b)
a
a

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 90 of 104


Reference
[1] Dwivedy S. K., and Eberhard P., 2006, Dynamic Analysis of Flexible Manipulators, Mechanism
and Machine Theory, 41, pp 749-777.
[2] Moon, F.C., and Pao, Y. H., 1969, Vibration and Dynamic Instability of a Beam-Plate in a
Transverse Magnetic Field, J ournal of Applied Mechanics, 36, pp. 92100.
[3] Wu, G.Y., Tsai, R., and Shih, Y.S., The Analysis of Dynamic Stability and Vibration Motions of a
Cantilever Beam with Axial Loads and Transverse Magnetic Fields, J ournal of the Acoustical Society
of ROC, 4, pp. 4055.
[4] Chen, C.C., and Yah, M.K., 2001, Parametric Instability of a Beam under Electromagnetic
Excitation, J ournal of Sound and Vibration, 240, pp. 747764.
[5] Wu, G.Y., 2005, The Analysis of Dynamic Instability and Vibration Motions of a Pinned Beam
with Transverse Magnetic Fields and Thermal Loads, J ournal of Sound and Vibration, 284, pp. 343
360.
[6] Wu, G.Y., 2007, The Analysis of Dynamic Instability on the Large Amplitude Vibrations of a
Beam with Transverse Magnetic Fields and Thermal Loads, J ournal of Sound and Vibration, 302,
pp.167-177.
[7] Pratiher, B., and Dwivedy, S. K., 2007, Parametric Instability of a Cantilever Beam with Magnetic
Field and Periodic Axial Load, J ournal of Sound and Vibration, 305, pp.904-917.
[8] Kojima, H., and Nagaya, K., 1985, Nonlinear Forced Vibration of a Beam with a Mass Subjected
to Alternating Electromagnetic Force, Bull. of the J apan Society of Mechanical Engineers, 28, pp.
468474.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 91 of 104


[9] Lu, Q. S., To, C.W.S., and Huang, K.L., 1995, Dynamic Stability and Bifurcation of an Alternating
Load and Magnetic Field Excited Magneto-Elastic Beam, J ournal of Sound and Vibration, 181, pp.
873891.
[10] Nayfeh, A.H., and Mook, D.T., 1995, Nonlinear Oscillations, Wiley, New York.
[16] Cuvalci, O., 2000, The Effect of Detuning Parameters on the Absorption Region for a Coupled
System: a Numerical and Experimental Study, J ournal of Sound and Vibration, 229, pp. 837-857.


Module 6 Lect 12
Parametrically excited System with internal Resonance
In this lecture a case study is taken for a parametrically excited system with internal resonance. The
system considered is a uniform cantilever beam of length carriying a mass at an arbitrary position
from the fixed end and subjected to base motion ( ) t Z t z = cos
0
as shown in Fig. 6.12.1.

Fig. 6. 12. 1. Base Excited Cantilever beam with attached mass at arbitrary position

The equation of motion of the beam is given by Kar and Dwivedy[1999] as
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 92 of 104


( )
{ }
( )
{ }
( )( ) ( )
2 3 2
0
1 1
3 1
2 2
0
ssss s ssss s ss sss ss s tt t
L
s sss tt t s s
tt
s
s
EI m s d c
m d c d j s d N


| |
( + + + + + +
`
|

) \ .
( ( + + + =

(6.12.1)
Subject to the boundary conditions
( ) ( ) ( ) ( ) 0 , , 0 , , 0 , 0 , 0 , 0 = = = = t L t L t t
sss ss s
(6.12.2)
Where
( ) ( ) ( ) ( )
( ) ( ) ( )
2 2
0
2 2
0
1 1
2 2
1
1
2
s L
s s tt
tt tt
L s s
s
tt stt s s st
L
N d d m d d d m z g
s
d d L z g J s d
L





= + +
` `

) )
| |
+ +
`
|
\ . )

(6.12.3)

Here, ( )
( )
( )
( )
s t
s t

= ,
Here, E, I and are, respectively, the Youngs modulus, the second moment of area of the cross-section
of the beam and mass per unit length of the beam;
0
is the moment of inertia of the concentrated mass
about its centroidal axis perpendicular to the X-Y plane; is the lateral displacement of the beam;
, and are, respectively, the acceleration due to gravity, the coefficient of viscous damping and the
displacement of the base; and is the Dirac delta function. Assuming a solution of Eq. (6.12.1) in the
form
( ) ( ) ( ) ,
n n
n
s t r s u t

=
=

1
(6.12.4)
Wherer is a scaling factor,
n
s is the shape function of the n th mode, and

is the time
modulation of the n
th
mode. Applying generalized Galerkins method and using the following non-
dimensional parameters,
, , , ,
1
1

n
n
t
L
d
L
s
x = = = =
1
2
0 0
, , , ,

= = = = =
Lr
J
J
Z
Z
T
L
m
L
r
r
(6.12.5)

Eq. (6.12.1) reduces to
{ } 2
1 1 1 1
0,
2 cos
n n n
klm k l m klm k l m klm k l m
n n n n n nm m
m k l m
u u u u u u u u u
u u u f u



= = = =
+ + =
+ + +


(6.12.6)
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 93 of 104


Here, 1,2,... n = , ( ) ( ) d d /
.
= . For details of the coefficients one may refer (Kar and Dwivedy 1999).
The small dimensionless parameter is the book-keeping parameter to indicate the smallness of
damping, non-linearities and excitation. So, we have number of coupled equations with cubic
geometric and inertial non-linearities, where represents the number of modes participating in the
resulting oscillation. Due to the absence of any internal and external excitation for 3 n , the amplitude
of these higher modes die out in the presence of damping and hence two mode discretization in the
Galerkins method is sufficient in this particular system.
The approximate solution of Eq. (6.12.6) can be obtained using the method of multiple scales. Let
( ) ( ) ( ) ... , , ;
1 0 1 1 0 0
+ + = T T u T T u u
n n n
(6.12.7a)
= = = ,..., 2 , 1 , ,
1 0
n T T (6.12.7b)
Substituting Eqs. (6.12.7a) and (6.12.7b) into Eq. (6.12.6) and equating the coefficients of
0
and to
zero, yields
, 0
0
2
0
2
0
= +
n n n
u u D (6.12.8)

=
+ = +
1 ,
0 0 1 0 0 0 1
2
1
2
0
cos 2 2 [
m n
m nm n n n n n n
u f u D D u D u u D
0
0
2
0 0 0 0 0 0 0 0 0 0 0
= + + +

klm
m l k
n
klm m l k
n
klm m l k
n
klm
u D u u u D u D u u u u (6.12.9)
Where
0 0
/ T D = and
1 1
/ T D = . The solution of Eq. (9) is given by
( ) ( ) cc T i T A u
n n n
+ =
0 1 0
exp (6.12.10)
Where cc indicates the complex conjugate of the preceding terms and
n
A is determined in the following
section. Considering the Principal parametric resonance ( )
1
2 , to express the nearness of to
1
2 the detuning parameter
1
is introduced. Also, to account for the internal resonance, the detuning
2
is used. Hence, we have
, = + = +
1 1 2 1 2
2 3 (6.12.11)
Substituting Eqs. (6.12.10) and (6.12.11) into Eq. (6.12.9) and eliminating the secular terms, we get for
= 1
( ) ( ) ( ) { } | |
0 1 2 2 12 0 1 1 11
'
1 1 1 1
exp exp
2
1
2 T i A f T i A f A A i + +
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 94 of 104


( )

=
= + +
1
0 2
2
1 2 12 1 1
0 exp
j
j j j e
T i A A Q A A A (6.12.12)
For n=2,
( ) ( ) { }
0 2 1 1 21
'
2 2 2 2
exp
2
1
2 T i A f A A i +
( )

=
= + +
1
0 2
3
1 2 21 1 2
0 exp
j
j j j e
T i A A Q A A A (6.12.13)
For 3 n ,
( )

=
= + +
1
'
0 2
j
n j j enj n n n n
A A A A A i (6.12.14)
Where a prime denotes the derivative with respect to
1
. Since the higher modes ( ) 3 n are neither
directly excited by external excitation nor indirectly excited by internal resonance, from Eq. (6.12.14) it
can be shown that the response amplitude of these modes die out due to the presence of damping.
Letting ( ) ( ) { }
1 1
exp
2
1
T i T a A
n n
= (where
n
a and
n
are real) in Eqs. (6.12.12) and (6.12.13) and then
separating into real and imaginary parts, one obtains
( ) ( ) { }
2 1 2 12 1 1 11
'
1 1 1 1
sin 2 sin
2
1
2 + + a f a f a a ( ) 0 3 sin 25 . 0
2 1
2
1 2 12
= + a a Q (6.12.15a)
( ) { }
( )
'
1 1 1 1 11 1 1 12 2 1 2
2
2 2
1 1 12 2 1 1 2
1
1 1
2 cos2 cos
2 2
1 1
cos 3 0
4 4
e j j
j
a f a f a
a a Q a a


=
| |
+ +
|
\ .
+ =

(6.12.15b)
( ) ( ) { } ( )
' 3
2 2 2 2 21 1 2 1 21 1 2 1
1 1
2 sin sin 3 0
2 4
a a f a Q a + + = (6.12.15c)
( ) ( ) ( )
2
' 2 3
2 2 2 2 1 21 1 2 1 2 2 21 1 2 1
1
1 1 1
2 1.5 cos cos 3 0
2 4 4
e j j
j
a f a a a Q a
=
+ + + =

(6.12.15d)
Where
( )
1 1 1 1 2 2 1 2 1
1
, and 1.5
2
T T = + = +
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 95 of 104


The above equations are known as the reduced equations. For steady state, 0
'
2
'
2
'
1
'
1
= = = = a a . So,
now we have a set of non-linear algebraic equations which is solved numerically to obtain the fixed
point response of the system. The first-order solution of the system can be given by
( ) { }
1 1 1 1 1
2 / cos + = a u (6.12.16a)
( ) { } | |
2 2 1 2 2 2
5 . 1 cos + = a u (6.12.16b)
Stability equations of steady-state response
By directly perturbing the reduced equations, one can study the stability of the non-trivial steady sate
solution. But, as the reduced Eqs. (6.12.15a-d) have the coupled terms
'
1 1
a and
'
2 2
a , the perturbed
equations will not contain the perturbations
'
1
or
'
2
for trivial solutions and hence the stability of
the trivial state cannot be studied by directly perturbing these equations. To circumvent this difficulty,
normalization method is adopted by introducing the transformation
2 , 1 , sin , cos = = = i a q a p
i i i i i i
(6.12.17)
Into equations (6.12.15) to obtain the following normalized reduced equations or the Cartesian form of
modulation equations:
( )
2 12 1 11 1 1 1 1
'
1 1
2
1
2
1
2 q f q f p p + |
.
|

\
|
+ +
( ) { } ( )

=
= + + +
2
1
2 2
1 1 1 2 1
2
1
2
1 2 12
0
4
1
2
4
1
j
j j j e
q p q q p p p q q Q (6.12.18a)
( )
2 12 1 11 1 1 1 1
'
1 1
2
1
2
1
2 p f p f q q + |
.
|

\
|
+ +
( ) { } ( )

=
= + + +
2
1
2 2
1 1 2 1 1
2
1
2
1 2 12
0
4
1
2
4
1
j
j j j e
q p p q q p q p p Q (6.12.18b)
( ) ( )
2 2 1 2 1 21 2 2
'
2 2
2 3
2
1
2 q q f p p + + +
( ) ( )

=
=
2
1
2 2
2 2
2
1
2
1 1 21
0
4
1
3
4
1
j
j j j e
q p q q p q Q (6.12.18c)
( ) ( )
2 2 1 2 1 21 2 2
'
2 2
2 3
2
1
2 p p f q q +
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 96 of 104


( ) ( )

=
= + + +
2
1
2 2
2 2
2
1
2
1 1 21
0
4
1
3
4
1
j
j j j e
q p p q p p Q (6.12.18d)
Now perturbing the above equations, one obtains
{ }
T
q p q p
'
2
, '
2
'
1
'
1
, , , | |{ }
T
c
q p q p J
2 2 1 1
, , , = (6.12.19)
Where T is the transpose and [

] is the J acobian matrix whose eigenvalues will determine the stability

and bifurcation of the system.
The stability boundary for the linear system (i.e. the trivial state) can be obtained from the eigen values
of the matrix [

] by letting 0
2 1 2 1
= = = = q q p p .
The first order solution of the system in terms of ( ) 2 , 1 , = i q p
i i
can be given by
, sin cos
1 1 1 1 1
q p u + = (6.12.20a)

1 2 1 2 2
3 sin 3 cos q p u + = , (6.12.20b)
where
1 1 1
2
1
+ = (6.12.21)
If the external frequency
m n
= where
n
is the nth natural frequency of the system one will
obtain combination resonance of sum( )
m n
= + or difference ( ) -
m n
= type for which one
may refer to the work of Dwivedy and Kar (1999). Also an exhaustive list of literature is given for the
interested reader.
Numerical Results and Discussion
Following Zavodney and Nayfeh [7] and keeping internal resonance in view, a metallic beam is
considered with the following properties:
L=125.4 mm, I=0.04851
4
, E=0.20936

2
,

= 1 , = 0.1 . /
2
,
= 0.03332

, = 3.68979, J =0.959, = 0.25

The roots of the characteristics equation are found numerically to be
1
=1.80097,
2
=3.2836 and the
corresponding non-dimensional natural frequencies are
1
=1 and
2
=3.33179. The book keeping
parameter and scaling factor are taken as 0.001 and 0.1, respectively. The coefficients of damping
( )
n
, excitation ( )
mm
f and non-linear terms ( )
n
klm
n
klm
n
klm
, , are found to be of the same order. The
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 97 of 104


values of other required parameters expressed in the appendix (Kar and Dwivedy 1999) are calculated
to be

11
=2.54149,
12
=-12.2027,
21
=-6.63699,
22
=-195.55,
12
=14.62282,
21
=7.84674,

11
f =0.0655762,

12
f =0.0122118,

21
f =0.04249,

22
f =0.1699298,

1
=0.0118963,

2
=0.0045865
Figure 6.12.2 shows the trivial state instability regions for the system with principal parametric
resonance for different damping parameters. While the region bounded by the curves is unstable the
regions outside the curves are stable. Clearly due to the presence of internal resonance, in addition to
the main unstable region near =2, additional alternate zones of stable and unstable trivial branches
exists. With increase in damping and forcing amplitude these additional zones get merged with the
main unstable region. Here it may be noted that while with increase in damping the instability region
decreases, with increase in forcing amplitude, the instability region increases.

Fig. 6.12.2: A typical principal instability region
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 98 of 104



Fig. 6.12.3: Frequency response curve 8, 0.001

A typical frequency response curve is shown in Figure 6.12.3 for both the first mode (lines without
bullet point) and the second mode (lines with bullet point). The stable branches are shown by solid lines
and the unstable branches are shown by dotted line. One may observe multi stable regions for a wide
range of frequency of the system. The nontrivial response amplitude of the first mode is observed to be
larger than the second mode. While supercritical and subcritical pitchfork bifurcations are observed in
the trivial state, both saddle node and Hopf bifurcations are observed in the nontrivial state. Due the
presence of Hopf bifurcation stable periodic response occurs in the trivial unstable region. With
decrease in forcing amplitude and damping parameter figure 6.12.4 shows the frequency response curve
for both the modes for 5, 0.01 . In addition to the other phenomenon described in the previous
figure, here one may clearly observe (in the insert) the additional alternate stable and unstable trivial
states near the main unstable region. While increasing the frequency one may observe jump up
phenomena and while decreasing the frequency one may observe the jump down phenomena in the
system.

NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 99 of 104



Fig. 6.12.4: Frequency response curve for 5, 0.01

Jump up
Jump down
Jump up
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 100 of 104


Figure 6.12.5 shows the force response curve =2.0, =0.001.


Figure 6.12.6 shows the force response curve =1.75, =0.001.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 101 of 104



Figure 6.12.7 Phase portraits: cascade of period doubling leading to chaos for =2.13, 8.0 , (a) =
8.5 (periodic), (b) =8.4(2T periodic), =8.3 (chaotic orbits).

Figure 6.12.8 Poincare showing period doubling route to chaos for =2.13
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 102 of 104


Figure 6.12.5 shows the force response curve for the perfectly tuned system at =2.0, =0.001.
Herethe nontrivial fixed point becomes unstable with saddle nodde bifurcation points at 0.9 and
1.6.And Hopf bifurcation points at 5.75 . Here the trivial branch istotally unstable except at 0
.Similarly Figure 6.12.6 shows the force response curve of the system at =1.75, =0.001. Though
the trivial response loses its stability at 7.55 through supercritical Pitchfork bifurcation, and a
Hopf bifurcation is observed at 7.95 , the system will fail through blue sky catastrophe if the
amplitude of excitation is increased beyond the
turning point at 10.75 . Figure 6.12.7(a) shows the periodic response originating from the Hopf
bifurcation for =2.13, 8.0 , =8.5. With decrease in the damping parameter to 8.42 one may
observe a response with double period (Figure 6.12.7(b)). This period doubling phenomena continues
with further decrease in damping parameter and finally a chaotic response (Figure 6.12.7(c)) is
observed. Figure 6.12.8 shows the Poincare section depicting cascade of period-doubling leading to
chaos.

Reference for further study
1. K.L. Hadoo, V. Sundarajan, Parametric instability of a cantilevered column with end mass,
J ournal of Sound and Vibration 18, 45-53, 1971

2. E.G. Tezak, D.T. Mook, A.H. Nayfeh, Nonlinear analysis of the lateral response of coloums to
periodic loads, ASME J ournal of Machine Design, 100, 651-655, 1978.
3. C.W.S. To, Vibration of a cantilever beam with a base excitation and tip mass, J ournal of Sound
and Vibration 83, 445-460, 1982.
4. A.H. Nayfeh, The response of multi degree freedom systems with quardtic non-linearities to a
harmonic parametric resonance, J ournal of Sound and Vibration, 90, 237-244, 1983.
5. A.H. Nayfeh, A.E.S. J ebrill, The resonance of two degree of freedom systems with quardtic
non-linearities to parametrtic excitation, J ournal of Sound and Vibration 88, 547-557, 1983.
6. L. D. Zavodney, A theoretical and experimental investigation of parametrically excited non-
linear mechanical system. PhD Dissertation, Dept. of. Engineering Science and Mechanics,
Virginia Polytechnic and State University, Blacksburg, December 1987.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 103 of 104


7. L. D. Zavodney and A. H. Nayfeh, The nonlinear response of a slender beam carrying a lumped
mass to a principal parametric excitation: theory and experiment, International J ournal of Non-
linear Mechanics, 24, 105-125, 1989.
8. A.H. Nayfeh, P.F. Pai, Nonlinear non-planar parametric responses of an in-extensional beam,
International J ournal of Non-linear Mechanics, 24, 139-158, 1989.
9. A.H. Nayfeh, B. Balachandram, Modal interactions in dynamical and structrual systems, ASME
Applied Mechanics Review, 42, s175-s201, 1989.
10. H.P. Lee, Stability of a cantilever beam with tip mass subjected to axial sinusoidal excitations,
J ournal of Sound and Vibration 183, 91-98, 1995.
11. B. Banarjee, A.K. Bajaj, P. Davies, Resonant dynamics of an autoparametric systems: a study
using higher order avearging, International J ournal Mechanical, 31, 21-39, 1996.
12. E.H.K. Fung, Z.X. Shi, Vibration frequencies of a constrained flexible arm carrying an end
mass, J ournal of Sound and Vibration 203, 373-387, 1997.
13. M.N. Hamdan, M.H.F. Dado, Large amplitude free vibrations of a uniform cantilever beam
carrrying an intermediate lumped mass and rotary inertia, J ournal of Sound and Vibration 206,
151-168, 1997.
14. H. Erol, M. Grgze, On The Effect of an Attached Spring-mass System On The Frequency
Spectrum Of Longitudinally Vibrating Elastic Rods Carrying A Tip Mass, J ournal of Sound and
Vibration, 227, (2) 471-477, 1999.
15. S. Inceo , M. Grgze, Bending Vibrations Of Beams Coupled By Several Double Spring-
mass Systems, J ournal of Sound and Vibration, 243 (2), 370-379, 2001.
16. H Erol, M Grgze, Longitudinal vibrations of a double-rod system coupled by springs and
dampers, J ournal of Sound and Vibration, 276 (12), 419-430, 2004.
17. M. Grgze, H. Erol , On laterally vibrating beams carrying tip masses, coupled by several
double springmass systems, J ournal of Sound and Vibration, 269 (12), 431-438, 2004.
18. M. Grgze On the representation of a cantilevered beam carrying a tip mass by an equivalent
springmass system J ournal of Sound and Vibration, 282, (12), 538-542, 2005.
19. S. K. Dwivedy and R. C. Kar, Dynamics of a slender beam with an attached mass under
combination parametric and internal resonances, Part II: Periodic and Chaotic response. J ournal
of Sound and Vibration, 222, no 2, 281-305, 1999.
NPTEL Mechanical Engineering Nonlinear Vibration


Joint initiative of IITs and IISc Funded by MHRD Page 104 of 104


20. S. K. Dwivedy and R. C. Kar, Dynamics of a slender beam with an attached mass under
combination parametric and internal resonances, Part I: steady state response. J ournal of Sound
and Vibration, 221, no 5, 823-848, 1999
21. R. C. Kar and S. K. Dwivedy, Non-linear dynamics of a slender beam carrying a lumped mass
with principal parametric and internal resonances. International J ournal of Nonlinear
Mechanics, 34, (3), 515-529, 1999.
22. S. K. Dwivedy and R. C. Kar, Simultaneous combination, principal parametric and internal
resonances in a slender beam with a lumped mass: three mode interactions. J ournal of Sound
and Vibration Vol. 242, no 1 pp 27-46, 2001.
23. S. K. Dwivedy and R. C. Kar, Non-linear dynamics of a slender beam carrying a lumped mass
under principal parametric resonances with three-mode interactions. International J ournal of
Nonlinear Mechanics, Vol. 36, no. 6, pp. 927-945, 2001.
24. S. K. Dwivedy and R. C. Kar, Nonlinear response of a parametrically excited system using
higher order multiple scales. Nonlinear Dynamics, 20, 115-130, 1999.
25. S. K. Dwivedy and R. C. Kar, Nonlinear dynamics of a cantilever beam carrying an attached
mass with 1:3:9 internal resonances. Nonlinear Dynamics, 31(1), 49-72, 2003.
26. S. K. Dwivedy and R. C. Kar, Simultaneous combination and 1:3:5 internal resonances in a
parametrically excited beam-mass system. International J ournal of Nonlinear Mechanics, 38,
(4), 585-596, 2003.
27. S. K. Dwivedy and R. C. Kar, Nonlinear response of a parametrically excited slender beam
carrying a lumped mass with 1:3:9 internal resonance, Advances in Vibration Engineering. 2,
(3), 1-9, 2003.