Analytical study on the vibration frequencies of tapered beams

Analytical study on the vibration frequencies of tapered beams

Mahmoud Bayat^*; Iman Pakar; Mahdi Bayat

Department of Civil Engineering, Shirvan Branch, Islamic Azad University, Shirvan, Iran

ABSTRACT

A vast amount of published work can be found in the field of beam vibrations dealing with analytical and numerical techniques. This paper deals with analysis of the nonlinear free vibrations of beams. The problem considered represents the governing equation of the nonlinear, large amplitude free vibrations of tapered beams. A new implementation of the ancient Chinese method called the Max-Min Approach (MMA) and Homotopy Perturbation Method (HPM) are presented to obtain natural frequency and corresponding displacement of tapered beams. The effect of vibration amplitude on the non-linear frequency is discussed. In the end to illustrate the effectiveness and convenience of the MMA and HPM, the obtained results are compared with the exact ones and shown in graphs and in tables. Those approaches are very effective and simple and with only one iteration leads to high accuracy of the solutions. It is predicted that those methods can be found wide application in engineering problems, as indicated in this paper.

Keywords: large amplitude free vibrations, analytical solution, tapered beam, Homotopy Perturbation Method (HPM), Max-Min Approach (MMA).

1 INTRODUCTION

Analyzing the nonlinear vibration of beams is one of the important issues in structural engineering. Applications such as high-rise buildings, long-span bridges, aerospace vehicles have necessitated the study of their dynamic behavior at large amplitudes. Many researchers have studied tapered beams, which are very important for the design of many engineering structures. The non-linear vibration of beams is governed by a non-linear partial-differential equation in space and time. Generally, finding an exact or close-form solution for nonlinear problems is very dificult. Goorman [7] is given the governing differential equation corresponding to fundamental vibration mode of a tapered beam. Evensen [6] studied on the nonlinear vibrations of beams with various boundary conditions by using the perturbation method. Pillai & Rao [24] considered the different types of solutions to the nonlinear equation of motion such as Galerkin, harmonic balance method and simple harmonic oscillations. Singh et al [32] later studied the large-amplitude vibration problem of unsymmetrically laminated beams based on classical, first-order and higher-order formulations by using the numerical integration technique introduced earlier by Singh et al [31]. Qaisi [25] used an analytical method for determining the vibration modes of geometrically nonlinear beams under various edge conditions. Rehfield [27] proposed an approximate method for nonlinear vibration problems with material nonlinear effects for various boundary conditions. Sathyamoorthy [28] developed the work on finite element method for nonlinear beams under static and dynamic loads and classical methods for the analysis of beams with material, geometric and other types of nonlinearities. Raju et al [26] studied the large amplitude vibration problem of beams and plates using Rayleigh-Ritz method by incorporating the inplane deformation as well as inertia, which were absent in the earlier studies, and also by retaining the equivalent linearization function. Klein [13] used finite element approach and Rayleigh-Ritz for analyzing the vibration of the tapered beams. A dynamic discretization technique was applied to calculate the natural frequencies of a non-rotating double tapered beam based on both the Euler-Bernoulli and Timoshenko Beam Theories by Downs [5]. Sato [29] improved the Ritz method to study a linearly tapered beam with ends restrained elastically against rotation and subjected to an axial force. Lau [14] used the exact method for studying on the free vibration of tapered beam with end mass. The Green's function method in Laplace transform domain was used to study the vibration of general elastically restrained tapered beams by Lee et al [15] for obtaining the approximate fundamental solution by using a number of stepped beams to represent the tapered beam. Junior et al [21] proposed Galerkin Method and the Askey-Wiener scheme as solutions of the stochastic beam bending problem. Oni and Awodola considered [3] the Dynamic response under a moving load of an elastically supported non-prismatic Bernoulli-Euler beam on variable elastic foundation. Boukhalfa and Hadjoui [4] analyzed the free vibration of an embarked rotating composite shaft using the hp-version of the FEM.

The main objectives of this study are to use analytical methods for analyzing the free vibration of the tapered beams. Finding an exact analytical solution for nonlinear equations is extremely dificult. Therefore, many analytical and numerical approaches have been investigated. The most useful methods for solving nonlinear equations are perturbation methods. They are not valid for strongly nonlinear equations and they have many shortcomings. Many new techniques have appeared in the open literature to overcome the shortcomings of traditional analytical methods such as Variational Iteration [10, 18], Parameter-Expansion [19], Energy Balance [2, 12, 16, 22], Variational Approach [33], Iteration Perturbation [23], the Improved Amplitude-frequency Formulation [35] and differential transformation [20], etc. In this study, which is an extension of the authors' previous work [1], Max-Min Approach (MMA) and Homotopy Perturbation Method (HPM) were developed by J.H. He [8, 9] and investigated in different works [17, 30, 34], which have the following advantages, over above-mentioned methods:

1. MMA and HPM lead us to a very rapid convergence of the solution and they can be easily extended to other nonlinear oscillations.

2. In MMA and HPM, just one iteration leads us to high accuracy of solutions which are valid for a wide range of vibration amplitudes.

The Max-Min Approach (MMA) and Homotopy Perturbation Method (HPM) are used to find analytical solutions for this problem with the nonlinear governing differential equation. It is shown that the solutions are quickly convergent and their components can be simply calculated. The results of the MMA and HPM are compared with the exact one, it can be observed that MMA and HPM are accurate and require smaller computational effort. An excellent accuracy of the MMA and HPM results indicates that those methods can be used for problems in which the strong nonlinearities are taken into account.

2 TAPERED BEAM FORMULATION

In dimensionless form, the governing differential equation corresponding to fundamental vibration mode of a tapered beam is given by Goorman [7] and the schematic of a tapered beam represented by Fig. 1:

Where u is displacement and ε₁ and ε₂ are arbitrary constants. Subject to the following initial conditions:

3 OVERVIEW OF HE'S MAX-MIN APPROACH METHOD

We consider a generalized nonlinear oscillator in the form;

Where ƒ(u) is a non-negative function of u. According to the idea of the max-min method, we choose a trial-function in the form;

Where ω the unknown frequency to be further is determined.

Observe that the square of frequency, ω², is never less than that in the solution

Of the following linear oscillator

Where ƒ_min is the minimum value of the function ƒ(u).

In addition, ω²never exceeds the square of frequency of the solution

Of the following oscillator

Where ƒ_max is the maximum value of the function ƒ(u).

Hence, it follows that

According to the Chentian interpolation [8, 11], we obtain

Or

Where m and n are weighting factors, k =n/m. So the solution of Eq. (3) can be expressed as

The value of k can be approximately determined by various approximate methods [22, 23, 33]. Among others, hereby we use the residual method [11]. Substituting Eq. (12) into Eq. (3) results in the following residual:

Where ω =

If, by chance, Eq. (12) is the exact solution, then R(t; k) is vanishing completely. Since our approach is only an approximation to the exact solution, we set

Where T =2π/ω. Solving the above equation, we can easily obtain

Substituting the above equation into Eq. (12), we obtain the approximate solution of Eq. (3).

4 OVERVIEW OF HOMOTOPY PERTURBATION METHOD (HPM)

To explain the basic idea of the HPM for solving nonlinear differential equations, one may consider the following nonlinear differential equation:

That is subjected to the following boundary condition:

Where A is a general differential operator, B a boundary operator, ƒ(r) is a known analytical function, Γ is the boundary of the solution domain (Ω), and ∂u/∂t denotes differentiation along the outwards normal to Γ. Generally, the operator A may be divided into two parts: a linear part L and a nonlinear part N. Therefore, Eq. (16) may be rewritten as follows:

In cases where the nonlinear Eq. (16) includes no small parameter, one may construct the following homotopy equation

Where

In Eq. (19), p ∈ [0, 1] is an embedding parameter and u₀ is the first approximation that satisfies the boundary condition. One may assume that solution of Eq. (19) may be written as a power series in p, as the following:

The homotopy parameter p is also used to expand the square of the unknown angular frequency ω as follows:

Or

where ω₀ is the coeficient of u(r) in Eq. (16) and should be substituted by the right hand side of Eq. (23). Besides, ω_i (i =1, 2,...) are arbitrary parameters that have to be determined.

The best approximations for the solution and the angular frequency ω are

when Eq. (19) corresponds to Eq. (16) and Eq. (24) becomes the approximate solution of Eq. (16).

5 APPLICATIONS

5.1 Solution using MMA

We can re-write Eq. (1) in the following form

We choose a trial-function in the form

Where ω the frequency to be is determined.

By using the trial-function, the maximum and minimum values of ω² will be:

So we can write:

According to the Chengtian's inequality [8, 11], we have

Where m and n are weighting factors, k = n/m + n. Therefore the frequency can be approximated as:

Its approximate solution reads

In view of the approximate solution, Eq. (32), we re-write Eq.(??) in the form

Substituting the trial function into Eq. (34) , and using Fourier expansion series, it is obvious that:

For avoiding secular term we set b₁ =0

Solving the above equation, we can easily obtain

Substituting Eq. (37) into Eq. (32), yields

According to Eqs. (38) and (38), we can obtain the following approximate solution:

5.2 Solution using HPM

Eq. (1) can be rewritten as the following form:

To explain the analytical solution, the solution u and the square of the unknown angular frequency ω are expanded as follows:

Substituting Eqs.(41) and (42) into Eq. (40) and equating the terms with identical powers of p, the following set of linear differential equations is obtained:

Solving Eq. (43) gives: u₀ = A cos ωt. Substituting u₀ into Eq. (44), yield:

For achieving the secular term, we use Fourier expansion series as follows:

Substituting Eq. (46) into Eq. (45) yields:

Avoiding secular term, gives:

From Eq. (42) and setting p = 1, we have:

Comparing Eqs. (48) and (49), we can obtain:

Solving Eq. (50), gives:

6 RESULTS AND DISCUSSIONS

Comparisons with the analytical methods and the exact one are presented to illustrate and verify the accuracy of the Max-Min Approach (MMA) and Homotopy Perturbation Method (HPM). The exact frequency ω_e for a dynamic system governed by Eq. (1) can be derived, as shown in Eq. (52), as follows:

To demonstrate the accuracy of the MMA and HPM , the procedures explained in previous sections are applied to obtain natural frequency and corresponding displacement of tapered beams. A comparison of obtained results from the Max-Min Approach and Homotopy perturbation method and the exact one is tabulated in Table 1 for different parameters A, ε₁ and ε₂.

From Table 1, the relative error of the analytical approaches is 2.90861% for the firstorder analytical approximations, for A = 10, ε₁ = 0.1 and ε₂ = 0.5. To further illustrate and verify the accuracy of these approximate analytical approaches, a comparison of the time history oscillatory displacement response for tapered beams with exact solutions is presented in Figs. 2-4. From Figs. 2 and 3, the motion of the system is a periodic motion and the amplitude of vibration is a function of the initial conditions. Fig. 4 presents the high accuracy of both approaches with the exact one for different values of amplitude, and shows the phase-space curves ((t) versus u(t) curve) of the Eqs. (51, 38) for amplitudes, u(0)= 0.5, 1, 2 and 3. It can be observed that the phase-space curve generated from MMA and HPM are close to that of the exact curve. The phase plot shows the behavior of the oscillator when ε₁ =0.5, ε₂ =0.1. It is periodic with a center at (0, 0). This situation also occurs in the unforced, undamped cubic Dufing oscillators.

The effect of small parameters ε₁ and ε₂ on the frequency corresponding to various parameters of amplitude (A) has been studied in Figs. 5 and 6 for ε₁ and ε₂. Also, the phase plane for this problem obtained from MMA and HPM has been shown in Fig. 7. It is evident that MMA and HPM show excellent agreement with the numerical solution using the exact solution and quickly convergent and valid for a wide range of vibration amplitudes and initial conditions. The accuracy of the results shows that the MMA and HPM can be potentiality used for the analysis of strongly nonlinear oscillation problems accurately.

7 CONCLUSIONS

In this paper, the MMA and HPM were employed to solve the governing equations of nonlinear oscillations of tapered beams. The analytical solutions yield a thoughtful and insightful understanding of the effect of system parameters and initial conditions. Also, Analytical solutions give a reference frame for the verification and validation of other numerical approaches. MMA and HPM are suitable not only for weak nonlinear problems, but also for strong nonlinear problems. The most significant feature of those methods is their excellent accuracy for the whole range of oscillation amplitude values. Also, it can be used to solve other conservative truly nonlinear oscillators with complex nonlinearities. The MMA and HPM solutions are quickly convergent and its components can be simply calculated. Also, compared to other analytical methods, it can be observed that the results of MMA and HPM require smaller computational effort and only the one iteration leads to accurate solutions. The successful implementation of the MMA and HPM for the large amplitude nonlinear oscillation problem considered in this paper further confirms the capability of those methods in solving nonlinear oscillation problems.

Received 23 Nov 2010

In revised form 21 Dec 2010

• [1] M. Bayat, M. Shahidi, A. Barari, and G. Domairry. The approximate analysis of nonlinear behavior of structure under harmonic loading. International Journal of the Physical Sciences, 5(7):10741080, 2010.
• [2] M. Bayat, M. Shahidi, A. Barari, and G. Domairry. Analytical evaluation of the nonlinear vibration of coupled oscillator systems. Zeitschrift für Naturforschung A - A Journal of Physical Sciences, 66(1-2):6774, 2011.
• [3] Abdelkrim Boukhalfa and Abdelhamid Hadjoui. Free vibration analysis of an embarked rotating composite shaft using the hp-version of the FEM. Latin American Journal of Solids and Structures, 7:105141, 2010.
• [4] Claudio Roberto Avila da Silva Junior, Andre Teolo Beck, and Edison da Rosa. Solution of the stochastic beam bending problem by Galerkin method and the Askey-Wiener scheme. Latin American Journal of Solids and Structures, 6:5172, 2009.
• [5] B. Downs. Transverse vibrations of cantilever beam having unequal breadth and depth tapers. ASME J. Appl. Mech., 44:737742, 1977.
• [6] D.A. Evensen. Nonlinear vibrations of beams with various boundary conditions. AIAA J., 6:370372, 1968.
• [7] D.J. Goorman. Free vibrations of beams and shafts. Appl. Mech., ASME, 18:135139, 1975.
• [8] J. H. He. Max-Min approach to nonlinear oscillators. International Journal of Nonlinear Sciences and Numerical Simulation, 9(2):207210, 2008.
• [9] J.H. He. Homotopy perturbation technique. Computer Methods in Applied Mechanics and Engineering, 178:257262, 1999.
• [10] J.H. He. Variational iteration method-a kind of non-linear analytical technique: Some examples. International Journal of Non-Linear Mechanics, 34(4):699708, 1999.
• [11] J.H. He and H. Tang. Rebuild of king fang 40 bc musical scales by He's inequality. Appl Mech Comput, 168(2):909914, 2005.
• [12] N. Jamshidi and D.D. Ganji. Application of energy balance method and variational iteration method to an oscillation of a mass attached to a stretched elastic wire. Current Applied Physics, 10(A202):484486, 2010.
• [13] L. Klein. Transverse vibrations of non-uniform beam. J. Sound Vibr., 37:491505, 1974.
• [14] J.H. Lau. Vibration frequencies of tapered bars with end mass. ASME J. Appl. Mech., 51:179181, 1984.
• [15] S.Y. Lee and Y.H. Kuo. Exact solution for the analysis of general elastically restrained non-uniform beams. ASME J. Appl. Mech., 59:205212, 1992.
• [16] I. Mehdipour, D.D. Ganji, and M. Mozaffari. Application of the energy balance method to nonlinear vibrating equations. Current Applied Physics, 10(1):104112, 2010.
• [17] Mo Miansari, Me Miansari, A. Barari, and G. Domairry. Analysis of blasius equation for flat-plate flow with infinite boundary value. International Journal for Computational Methods in Engineering Science and Mechanics, 11(2):7984, 2010.
• [18] S. T. Mohyud-Din. Modified variational iteration method for integro-differential equations and coupled systems. Z. Naturforsch, 65a:277285, 2010.
• [19] S. T. Mohyud-Din, M. A. Noor, and K. I. Noor. Parameter-expansion techniques for strongly nonlinear oscillators. Int. J. Nonlin. Sci. Num. Sim., 10(5):581583, 2009.
• [20] M. Omidvar, A. Barari, M. Momeni, and D.D. Ganji. New class of solutions for water infiltration problems in unsaturated soils. Geomechanics and Geoengineering: An International Journal, 5(2):127135, 2010.
• [21] S. T. Oni and T. O. Awodola. Dynamic response under a moving load of an elastically supported non-prismatic Bernoulli-Euler beam on variable elastic foundation. Latin American Journal of Solids and Structures, 7:320, 2010.
• [22] T. Ozis and A. Yildirim. Determination of the frequency-amplitude relation for a dufing-harmonic oscillator by the energy balance method. Computers and Mathematics with Applications, 54:11841187, 2007.
• [23] Turgut Özis and Ahmet Yildirim. Generating the periodic solutions for forcing van der Pol oscillators by the iteration perturbation method. Nonlinear Analysis Real world Applications, 10:19841989, 2009.
• [24] S.R.R. Pillai and B.N. Rao. On nonlinear free vibrations of simply supported uniform beams. J.Sound Vib., 159:527531, 1992.
• [25] M.I. Qaisi. Application of the harmonic balance principle to the nonlinear free vibration of beams. Appl. Acoust., 40:141151, 1993.
• [26] I.S. Raju, G.V. Rao, and K.K. Raju. Effect of longitudinal or inplane deformation and inertia on the large amplitude flexural vibrations of slender beams and thin plates. J. Sound Vib, 49:415422, 1976.
• [27] L.W.A. Rehfield. Simple, approximate method for analysing nonlinear free vibrations of elastic structures. J. Appl. Mech., ASME, 42:509511, 1975.
• [28] M. Sathyamoorthy. Nonlinear analysis of beams, Part II: Finite-element methods. Shock Vib. Dig., 14:718, 1982.
• [29] K. Sato. Transverse vibrations of linearly tapered beams with ends restrained elastically against rotation subjected to axial force. Int. J. Mech. Sci., 22:109115, 1980.
• [30] M. Shaban, D.D. Ganji, and M.M. Alipour. Nonlinear fluctuation, frequency and stability analyses in free vibration of circular sector oscillation systems. Current Applied Physics, 10(5):12671285, 2010.
• [31] G. Singh, G.V. Rao, and N.G.R. Iyengar. Reinvestigation of large amplitude free vibrations of beams using finite elements. J. Sound Vib., 143:351355, 1990.
• [32] G. Singh, G.V. Rao, and N.G.R. Iyengar. Analysis of the nonlinear vibrations of unsymmetrically laminated composite beams. AIAA J., 29:17271735, 1991.
• [33] L. Xu. Variational approach to solitons of nonlinear dispersive K(m,n) equations. Chaos, Solitons & Fractals, 37:137143, 2008.
• [34] Elçin Yusufoğlu. An improvement to homotopy perturbation method for solving system of linear equations. Computers & Mathematics with Applications, 58(11-12):22312235, 2009.
• [35] Hu-Li Zhang. Application of he's frequency-amplitude formulation to an force nonlinear oscillator. Int. J. Nonlin. Sci. Numer. Simul., 9(3):297300, 2008.

Publication Dates

• Publication in this collection
  15 June 2011
• Date of issue
  June 2011

History

• Reviewed
  21 Dec 2010
• Received
  23 Nov 2010