Rajesh Korabathina^{1} , Meera Saheb Koppanati^{2}
^{1, 2}Jawaharlal Nehru Technological University, Kakinada, 533003, India
^{1}Corresponding author
Mathematical Models in Engineering, Vol. 2, Issue 1, 2016, p. 2733.
Received 9 November 2015; accepted 19 January 2016; published 30 June 2016
Copyright © 2016 JVE International Ltd.
Every structure which is having some mass and elasticity is said to vibrate. Natural frequency is the one of the most important parameter associated with engineering vibration. In nature every structure has its own natural frequency. Whenever the natural frequency of the structure coincides with the frequency of external applied load, excessive deflections will occur and the structure will be failed. To avoid such condition one must be aware of the operating frequencies of the materials or structures under various conditions like simply supported, clamped and cantilever boundary conditions. There are many methods to evaluate the natural frequency of the structures. in this method the authors developed a method called “coupled displacement field method” which reduces computational efforts compared with the other methods and which is successfully applied for the HingedHinged boundary condition of a tapered (rectangular cross section) Timoshenko beam and calculated the fundamental frequency parameter values and compared the results with existing literature. The results obtained in this method are very close to the existing literature.
Keywords: linear free vibrations, coupled displacement filed method, tapered Timoshenko beams, taper ratio, slenderness ratio.
Research on vibrations of beams has been doing by so many authors for a long period of time. So far, many authors have found different methods to find the free vibration behaviour of shear flexible beams, some of them are as follows.
Abrate S. [1] analyzed the free vibration of nonuniform beams with general shape and arbitrary boundary conditions. Simple formulas are presented for predicting the fundamental natural frequency of nonuniform beams with various end support conditions.
Byoung Koo Lee et al. [2] studied on “Free vibrations of tapered Beams with general boundary condition”. This method involving the finding of the ordinary differential governing equation of beams which can be solved by numerical methods. the natural frequencies are calculated by combining the Runge Kutta method and the determinant search method.
De Rosa M. A. et al. [3] studied the dynamic behaviour of beams with linearly varying crosssection. in this method. The equation of motion is solved in terms of Bessel functions, and the boundary conditions lead to the frequency equation which is a function of four flexibility coefficients.
De Rosa M. A. et al. [4] studied the Natural vibration frequencies of Tapered beams by using EulerBernoulli beam theory in the presence of an arbitrary number of rotationally, axially and elastically flexible constraints. In this method, the dynamic analysis is performed by means of the socalled cell discretization method (CDM), according to which the beam is reduced to a set of rigid bars, linked together by elastic sections, where the bending stiffness and the distributed mass of the bars is concentrated.
FirouzAbadi R. D. et al. [5] studied the transverse free vibration of a class of variablecrosssection beams is investigated using the Wentzel, Kramers, Brillouin (WKB) approximation in which the governing equation of motion of the EulerBernoulli beam including axial force distribution is utilized to obtain a singular differential equation in terms of the natural frequency of vibration and a WKB expansion series is applied to find the solution.
Izabela Zamorska [6] used the Green’s function method for the free vibration problem of nonuniform BernoulliEuler beams, to find the Green’s function of the fourth order differential operator, occurring at the beam’s equation of motion, the power series method is proposed.
Mahmoud A. A. et al. [7] applied the differential transformation method (DTM) for free vibration analysis of beams with uniform and nonuniform cross sections. Natural frequencies and corresponding normalized mode shapes are calculated for different cases of cross section and boundary conditions. MATLAB code is designed to solve the differential equation of the beam using the differential transformation method.
Mehmet Cem Ece, et al. [8] studied the vibrations of an isotropic beam which has a variable crosssection. In this, the Governing equation is reduced to an ordinary differential equation in spatial coordinate for a family of crosssection geometries with exponentially varying width. Analytical solutions of the vibration of the beam are obtained for three different types of boundary conditions associated with simply supported, clamped and free ends. Natural frequencies and mode shapes are determined for each set of boundary conditions.
Mohamed Hussien Taha et al. [9] used a “mathematical model for vibrations of nonuniform flexural beams” they presented a simplified mathematical model for free vibrations of nonuniform viscoelastic flexural beams is presented. The mass intensity, the material damping intensity and the flexural stiffness of the beam are assumed varying as power functions along the beam. An analytical solution for the fourth order differential equation of beam vibration under appropriate boundary conditions is obtained by factorization and calculated Mode shapes and damped natural frequencies of the beam are obtained for wide range of beam characteristics.
Rossi R. E. and Laura P. A. A. [10] determined the natural frequencies and dynamic behavior vibration of linearly tapered beams subjected to different combinations of edge supports by finite element algorithmic procedures.
Stanisław Kukla, Izabela Zamojska [11] applied the Green’s function method in frequency analysis of a beam with varying cross section. The beam carries an arbitrary number of attached discrete systems. The exact solution of the problem concerns a beam with quadratically varying crosssection area.
Zhou D. and Cheung Y. K. [12] studied the vibrational characteristics of tapered beams with continuously varying rectangular crosssection of depth and breadth proportional to ${x}_{s}$ and ${x}_{t}$ respectively, where both $s$ and $t$ are arbitrary real numbers for a truncated beam and arbitrary positive numbers for a sharp ended beam and $x$ is the axial coordinate measured from the sharp end of the beam. They used the Bernoulli Euler theory of bending is used to describe the motion of the beam. They obtained the eigen frequency equation by the RayleighRitz method and The effect of the location the convergency is discussed.
In the present method the concept of Coupled Displacement Method is presented for the vibration analysis of tapered Timoshenko beams, where the transverse displacement is expressed in terms of total rotation and its derivatives depending on the transverse shear deformation theory, in which the number of un known coefficients is brought down to exactly $half\left(n\right)$ in the coupled displacement field method. To demonstrate the efficacy of the proposed method, the expressions for fundamental frequency are obtained for tapered shear flexible Timoshenko beams with hingedhinged type of boundary condition. The numerical results presented in this paper are matching very closely with the existing literature.
The concept of coupled displacement field is explained in detail. in the Couple Displacement Field Method (CDF) with the single term admissible function for $\theta $, the function for $w$ is derived using the coupling equation. The coupling equation has been derived from the Timoshenko beam theory and an admissible function for $\theta $ which satisfies all the applicable boundary conditions and the symmetric condition is assumed in the beam domain. The coupled lateral displacement $w$ distribution is evaluated from the coupling equation.
Fig. 1. Tapered Timoshenko beam with linearly varying height (depth)
Fig. 2. Tapered (rectangular cross section) Timoshenko HingedHinged beam
From The kinematics of a shear flexible beam theory (based on the Timoshenko beam theory):
where $\stackrel{\_}{u}$ is the axial displacement and $\stackrel{\_}{w}$ is the transverse displacements at an any point of the beam, $z$ is the distance of the any point from the neutral axis, $w$ is the transverse displacement and $\theta $ is the total rotation anywhere on the beam axis and $x$, $z$ are the independent spatial variables. The axial and shear strains are given by:
Now, the expressions for the strain energy ‘$U$’ and the work done $W$ by the externally applied loads are given by:
where $EI$ is the flexural rigidity, $GA$ is the shear rigidity, $k$ is the shear correction factor (taken as 5/6 in the present study), $p\left(x\right)$ is the static lateral load per unit length acting on the beam, $E$ is the Young’s modulus, $G$ is the shear modulus, $I$ is the area moment of inertia, $A$ is the area of cross section, $x$ is the axial coordinate and $L$ is the length of the beam. Applying the principle of minimization of total potential energy, as:
The following equilibrium equations can be obtained:
where $\theta $ is total rotation, $w$ is transverse displacement.
Eqs. (8) and (9) are coupled equations and can be solved for obtaining the solution for the static analysis of the shear deformable beams.
A close observation of Eq. (8) shows that it is dependent on the load term ‘$p$’ and Eq. (9) is independent of the load term ‘$p$’. Hence, Eq. (9) is used to couple the total rotation $\theta $ and the transverse displacement $w$, so that the two undetermined coefficients problem (for single term admissible function) becomes a single undetermined coefficient problem and the resulting linear free vibration problem becomes much simpler to solve.
The concept of coupled displacement field method is explained in detail. In Coupled Displacement Field method (CDF) with the single term admissible function for $\theta $, the function for $w$ is derived using the coupling Eq. (9). An admissible function for $\theta $ which satisfies all the applicable boundary conditions and the symmetric condition is assumed in the beam domain.
For hingedhinged tapered Timoshenko beam, assumed total rotation is:
where $a$ is the central lateral displacement of the beam which is also the maximum lateral displacement.
Rewriting Eq. (9):
By integrating the above equation, lateral displacement can be obtained as:
where $\gamma =EI/kGA$, $\varsigma =x/L$.
It may be noted here that because of the coupled displacement field concept $,$ the transverse displacement ($w$) distribution contains the same undetermined coefficient “$a$” as the $\theta $ distribution and satisfies all the applicable essential boundary and symmetric conditions:
For an assumed $\theta $ distribution, the coupled displacement field for the lateral displacement $w$ is evaluated, after the lateral displacement $w$ is calculated, linear free vibrations can be studied.
Using the principle of conservation of total energy at any instant of time $U+T=$ constant.
As the depth of the beam is not uniform throughout the length, the expression for $U$ and $T$ are given as:
where:
and $U$, $T$, $E$ are the strain energy, kinetic energy and young’s modulus respectively. $\rho $ is the mass density of the material of the beam, $\omega $ is the radian frequency, ${h}_{0}$ and ${h}_{L}$ are the height of the beam at left end $x=$0 and the right end $x=L$ respectively, ${A}_{o}$ and ${I}_{o}$ are cross sectional area and area moment of inertia at left end, $A$ is the area at any cross section, $I$ is the moment of inertia at any cross section, $\alpha $ is the taper ratio.
Substituting Eqs. (10), (14), (18) in Eq. (16) and after simplification:
Substituting Eqs. (10), (14), (18) in Eq. (17) and after simplification:
Applying the principle of minimization of total potential energy:
and after simplification, the nondimensional fundamental frequency parameter 𝜆 is obtained and as:
where $\lambda $ is the non dimensional fundamental frequency parameter, $\beta =L/r$ (slenderness ratio) and $r$ is radius of gyration
If the taper ratio $\alpha $ is zero, then the above equation becomes:
The above equation is applicable for uniform beam as the taper ratio is zero.
Table 1. Fundamental Frequency ($\lambda $^{1/2}) values for a HingedHinged (Rectangular cross section) tapered Timoshenko beam
Taper ratio ($\alpha $)

Slenderness ratio ($\beta $)


10

20

40

80

100

500


CDF Method

Ref
[10]

CDF Method

Ref
[10]

CDF Method

CDF Method

CDF Method

Ref
[10]

CDF Method

Ref
[10]


0

8.3913

8.388

9.4108

9.411

9.7471

9.8384

9.8496

9.850

9.8688

9.869

0.1

8.6917

8.683

9.8416

9.829

10.226

10.331

10.344

–

10.366

–

0.15

8.8435

–

10.059

–

10.469

10.581

10.595

–

10.619

–

0.20

8.9962

8.955

10.279

10.22

10.714

10.833

10.848

–

10.873

–

0.25

9.1496

–

10.499

–

10.960

11.086

11.102

–

11.129

–

0.30

9.3036

9.205

10.721

10.61

11.208

11.342

11.358

–

11.386

–

0.35

9.4581

–

10.944

–

11.457

11.598

11.616

–

11.646

–

0.4

9.6128

–

11.168

–

11.707

11.856

11.875

–

11.906

–

0.45

9.7676

–

11.3927

–

11.9596

12.116

12.1356

–

12.168

–

0.5

9.9226

–

11.6178

–

12.2124

12.377

12.3973

–

12.432

–

0.55

10.077

–

11.8435

–

12.4662

12.638

12.6602

–

12.696

–

0.6

10.232

–

12.0697

–

12.7209

12.901

12.9242

–

12.962

–

0.65

10.386

–

12.2963

–

12.9764

13.165

13.1892

–

13.229

–

0.7

10.540

–

12.523

–

13.232

13.430

13.455

–

13.497

–

0.75

10.694

–

12.750

–

13.489

13.696

13.722

–

13.765

–

0.8

10.847

–

12.977

–

13.747

13.963

13.989

–

14.035

–

0.85

11.000

–

13.204

–

14.005

14.230

14.258

–

14.305

–

0.9

11.151

–

13.431

–

14.264

14.498

14.527

–

14.577

–

0.95

11.303

–

13.658

–

14.523

14.767

14.797

–

14.849

–

1

11.453

–

13.885

–

14.783

15.036

15.068

–

15.122

–

The concept of Coupled displacement field method is used to determine the fundamental frequencies of tapered shear flexible beams, with most practically used boundary condition. The boundary conditions of the beams considered is hingedhinged, one with axially immovable ends. Suitable singleterm trigonometric functions is used to represent the total rotation $\theta $ in the present study, the corresponding coupled lateral displacement $w$ is derived using the present concept. The numerical results are obtained in terms of nondimensional fundamental frequency parameter, for various slenderness ratios ($\beta $) and taper ratios ($\alpha $).
The present results, in terms of slenderness ratios ($\beta $) and taper ratios ($\alpha $) for a tapered shear flexible hingedhinged beam, are presented in Table 1. For the sake of comparison and validation of the proposed method, the results obtained by present coupled displacement method are compared with Rossi R. E. and Laura P. A. A. and are also included in this table. The value obtained for nondimensional fundamental frequency parameter for the given uniform Timoshenko beam is exactly matching with the existing literature. This is because the admissible function assumed for the total rotation $\theta $ is exact for this classical problem.
This method is having less complexity compared to other energy methods as the computational efforts are being reduced to half compared to RayleighRitz method. This formulation can also be used for the study of the free vibration analysis of shear flexible structural members other than Timoshenko beams.