Non-linear vibration of Euler-Bernoulli beams

Non-linear vibration of Euler-Bernoulli beams

A.Barari^I,*; H.D. Kaliji^II; M. Ghadimi^III; G. Domairry^III

^IDepartment of Civil Engineering, Aalborg University, Sohngårdsholmsvej 57, 9000 Aalborg, Aalborg Denmark

^IIDepartment of Mechanical Engineering, Islamic Azad University, Semnan Branch, Semnan Iran

^IIIDepartment of Mechanical Engineering, Babol University of Technology, Babol Iran

ABSTRACT

In this paper, variational iteration (VIM) and parametrized perturbation (PPM) methods have been used to investigate non-linear vibration of Euler-Bernoulli beams subjected to the axial loads. The proposed methods do not require small parameter in the equation which is difficult to be found for nonlinear problems. Comparison of VIM and PPM with Runge-Kutta 4th leads to highly accurate solutions.

Keywords: Variational Iteration Method (VIM), Parametrized Perturbation Method (PPM), Galerkin method, non-linear vibration, Euler-Bernoulli beam.

1 INTRODUCTION

The demand for engineering structures is continuously increasing. Aerospace vehicles, bridges, and automobiles are examples of these structures. Many aspects have to be taken into consideration in the design of these structures to improve their performance and extend their life. One aspect of the design process is the dynamic response of structures. The dynamics of distributed-parameter and continuous systems, like beams, were governed by linear and nonlinear partial-differential equations in space and time. It was difficult to find the exact or closed-form solutions for nonlinear problems. Consequently, researchers were used two classes of approximate solutions of initial boundary-value problems: numerical techniques [28, 31], and approximate analytical methods [2, 26]. For strongly non-linear partial-differential, direct techniques, such as perturbation methods, were not utilized to solve directly the non-linear partial-differential equations and associated boundary conditions. Therefore first partial-differential equations are discretized into a set of non-linear ordinary-differential equations using the Galerkin approach and the governing problems are then solved analytically in time domain.

Approximate methods for studying non-linear vibrations of beams are important for investigating and designing purposes. In recent years, some promising approximate analytical solutions have been proposed, such as Frequency Amplitude Formulation [13], Variational Iteration [5, 6, 14, 17], Homotopy-Perturbation [3, 4, 7, 24], Parametrized-Perturbation [18], Max-Min [15, 19, 29], Differential Transform Method [16], Adomian Decomposition Method [22], Energy Balance [23, 30], etc.

Kopmaz et al. [20] considered different approaches to describing the relationship between the bending moment and curvature of an Euler-Bernoulli beam undergoing a large deformation. Then, in the case of a cantilevered beam subjected to a single moment at its free end, the difference between the linear and the nonlinear theories based on both the mathematical curvature and the physical curvature was shown. Biondi and Caddemi [8] studied the problem of the integration of the static governing equations of the uniform Euler-Bernoulli beams with discontinuities, considering the flexural stiffness and slope discontinuities.

The vibration problems of uniform Euler- Bernoulli beams can be solved by analytical or approximate approaches [10, 21]. Pirbodaghi et al. [25] studied non-linear vibration behaviour of geometrically non-linear Euler-Bernoulli beams subjected to axial loads using homotopy analysis method. Also, the effect of vibration amplitude on the non-linear frequency and buckling load is discussed. Burgreen [9] investigated the free vibrations of a simply supported buckled beam using a single-mode discretization. He pointed out the natural frequencies of buckled beams depend on the amplitude of vibration. Eisley [11, 12] used a single-mode discretization to investigate the forced vibrations of buckled beams and plates. He considered both simply supported and clamped-clamped boundary conditions. For a clamped-clamped buckled beam, Eisley [11, 12] used the first buckling mode in the discretization procedure. He obtained similar forms of the governing equations for simply supported and clamped-clamped buckled beams.

The main purpose of this study is to obtain the analytical expression for geometrically non-linear vibration of clamped-clamped Euler-Bernoulli beams fixed at one end. Geometric non-linearity arises from non-linear strain-displacement relationships. This type of nonlinearity is most commonly treated in the literature. Sources of this type of nonlinearity include mid-plane stretching, large curvatures of structural elements, and large rotation of elements. First, the governing non-linear partial differential equation using Galerkin method was reduced to a single non-linear ordinary differential equation. It was then assumed that only fundamental mode was excited. The later equation was solved analytically in time domain using VIM and PPM. Ultimately, VIM and PPM methods are compared with Runge-Kutta 4th method.

2 DESCRIPTION OF THE PROBLEM

Consider a straight beam on an elastic foundation with length L, a cross-section A, a mass per unit length µ, moment of inertia I, and modulus of elasticity E that subjected to an axial force of magnitude as shown in Fig. 1. It is assumed that the cross-sectional area of the beam is uniform and its material is homogenous. The beam is also modeled according to the Euler-Bernoulli beam theory. Planes of the cross sections remain planes after deformation, straight lines normal to the mid-plane of the beam remain normal, and straight lines in the transverse direction of the cross section do not change length. The first assumption ignores the in plane deformation. The second assumption ignores the transverse shear strains and consequently the rotation of the cross section is due to bending only. The last assumption, which is called the incompressibility condition, assumes no transverse normal strains. The last two assumptions are the basis of the Euler-Bernoulli beam theory [27].

The equation of motion including the effects of mid-plane stretching is given by:

Where C is the viscous damping coefficient, is a foundation modulus and U is a distributed load in the transverse direction.

Assume the non-conservative forces were equal to zero. Therefore Eq. (1) can be written as follows:

For convenience, the following non-dimensional variables are used:

Where R = (I/A)^0.5 is the radius of gyration of the cross-section. As a result, Eq. (2) can be written as follows:

Assuming W(X, t) = Φ(X)ψ(t)whereΦ(X)is the first eigenmode of the beam [32] and applying the Galerkin method, the equation of motion is obtained as follows:

Where α = α₁ +α₂F + K and α₁ , α₂ and β are as follows:

The Eq. (5) is the governing non-linear vibration of Euler-Bernoulli beams. The center of the beam subjected to the following initial conditions:

where A denotes the non-dimensional maximum amplitude of oscillation.

3 BASIC IDEA OF VARIATIONAL ITERATION METHOD

To illustrate the basic concepts of the VIM, we consider the following differential equation:

Where L is a linear operator, N a nonlinear operator and g(t) an inhomogeneous term. According to VIM, we can write down a correction functional as follows:

Where λ is a general Lagrange multiplier which can be identified optimally via the variational theory [17]. The subscript n indicates the nth approximation and _n is considered as a restricted variation [17], i.e. δ

_n = 0.

4 APPLICATION OF VARIATIONAL ITERATION METHOD

To solve Eq. (5) by means of VIM, we start with an arbitrary initial approximation:

From Eq. (5), we have:

Integrating twice yields:

Equating the coefficients of cos(ωt) in ψ₀ and ψ₁, we have:

And therefore,

Where δ

_n = 0 is considered as restricted variation.

Its stationary conditions can be obtained as follows:

Therefore, the multiplier, can be identified as

As a result, we obtain the following iteration formula:

By the iteration formula (20), we can directly obtain other components as:

Where ω is evaluated from Eq. (13).

In the same manner, the rest of the components of the iteration formula can be obtained.

5 APPLICATION OF PARAMETRIZED PERTURBATION METHOD

Equation of motion, which reads:

We let

By substituting Eq. (23) in Eq. (22):

We suppose that the solution of Eq. (24) and the constant a, can be expressed in the forms:

Substituting Eqs. (25) and (26) into Eq. (24) and equating coefficients of same powers of e yields the following equations:

Solving Eq. (27) we obtain:

Therefore, Eq. (28) can be re-written as:

Avoiding the presence of a secular terms needs:

Substituting Eq. (31) into Eq. (26)

Solving Eq. (30), we obtain:

Its first-order approximation is sufficient, and then we have:

Where the angular frequency can be written by Eq. (32).

6 RESULTS AND DISCUSSIONS

The behavior of ψ(A, t) obtained by VIM and PPM at α = π and β = 0.15 is shown in Figs. 2 and 3. Influence of coefficients β and α on frequency and amplitude has been investigated and plotted in Figs. 4 and 5, respectively. The comparison of the dimensionless deflection versus time for results obtained from VIM, PPM and Runge-Kutta 4th order has been depicted in Fig. 6 for α = π and β = 0.15, with maximum deflection at the center of the beam equal to five (A=5). The solutions are also compared for t=0.5 in Table 1. It can be observed that there is an excellent agreement between the results obtained from VIM and PPM with those of Runge-Kutta 4th order method [1].

7 CONCLUSIONS

In this paper, nonlinear responses of a clamped-clamped buckled beam are investigated. Mathematically, the beam is modeled by a partial differential equation possessing cubic non-linearity because of mid-plane stretching. Governing non-linear partial differential equation of Euler-Bernoulli's beam is reduced to a single non-linear ordinary differential equation using Galerkin method. Variational Iteration Method (VIM) and Paremetrized Perturbation Method (PPM) have been successfully used to study the non-linear vibration of beams. The frequency of both methods is exactly the same and transverse vibration of the beam center is illustrated versus amplitude and time. Also, the results and error of these methods are compared with Runge-Kutta 4th order. It is obvious that VIM and PPM are very powerful and efficient technique for finding analytical solutions. These methods do not require small parameters needed by perturbation method and are applicable for whole range of parameters. However, further research is needed to better understanding of the effect of these methods on engineering problems especially mechanical affairs.

Received 25 Oct 2010

In revised form 22 Feb 2011

• [1] D. A. Anderson and J.C. Tannehill. Computational Fluid Mechanics and Heat Transfer Hemisphere Publishing Corp, 1984.
• [2] L. Azrar, R. Benamar, and R.G. White. A semi-analytical approach to the non-linear dynamic response problem of S-S and C-C beams at large vibration amplitudes part I: general theory and application to the single mode approach to free and forced vibration analysis. J. Sound Vib, 224:183207, 1999.
• [3] A. Barari, B. Ganjavi, M. Ghanbari Jeloudar, and G. Domairry. Assessment of two analytical methods in solving the linear and nonlinear elastic beam deformation problems. Journal of Engineering, Design and Technology, 8(2):127 145, 2010.
• [4] A. Barari, A. Kimiaeifar, G. Domairry, and M. Moghimi. Analytical evaluation of beam deformation problem using approximate methods. Songklanakarin Journal of Science and Technology, 32(3):207326, 2010.
• [5] A. Barari, M. Omidvar, D.D. Ganji, and A. Tahmasebi Poor. An approximate solution for boundary value problems in structural engineering and fluid mechanics. Journal of Mathematical Problems in Engineering, pages 113, 2008. Article ID 394103.
• [6] A. Barari, M. Omidvar, Abdoul R. Ghotbi, and D.D. Ganji. Application of homotopy perturbation method and variational iteration method to nonlinear oscillator differential equations. Acta Applicandae Mathematicae, 104:161 171, 2008.
• [7] M. Bayat, M. Shahidi, A. Barari, and G. Domairry. On the approximate analysis of nonlinear behavior of structure under harmonic loading. International Journal of Physical Sciences, 5(7):10741080, 2010.
• [8] B. Biondi and S. Caddemi. Closed form solutions of euler-bernoulli beams with singularities. International Journal of Solids and Structures, 42:30273044, 2005.
• [9] D. Burgreen. Free vibration of a pin-ended column with constant distance between pin ends. Journal of Applied Mechanics, 18:135139, 1984.
• [10] A. Dimarogonas. Vibration for engineers Prentice-Hall, Inc., 2nd edition, 1996.
• [11] J. G. Eisley. Large amplitude vibration of buckled beams and rectangular plates. AIAA Journal, 2:22072209, 1964.
• [12] J. G. Eisley. Nonlinear vibration of beams and rectangular plates. ZAMP, 15:167175, 1964.
• [13] A. Fereidoon, M. Ghadimi, A. Barari, H.D. Kaliji, and G. Domairry. Nonlinear vibration of oscillation systems using frequency-amplitude formulation. Shock and Vibration, 2011. DOI: 10.3233/SAV-2011-0633.
• [14] F. Fouladi, E. Hosseinzadeh, A. Barari, and G. Domairry. Highly nonlinear temperature dependent fin analysis by Variational Iteration Method. Journal of Heat Transfer Research, 41(2):155165, 2010.
• [15] S.S. Ganji, A. Barari, and D.D. Ganji. Approximate analyses of two mass-spring systems and buckling of a column. Computers & Mathematics with Applications, 61(4):10881095, 2011.
• [16] S.S. Ganji, A. Barari, L.B. Ibsen, and G. Domairry. Differential transform method for mathematical modeling of jamming transition problem in traffic congestion flow. Central European Journal of Operations Research, 2011. DOI: 10.1007/s10100-010-0154-7.
• [17] J.H. He. Variational iteration method a kind of non-linear analytical technique: Some examples. International Journal of Nonlinear Mechanics, 34:699708, 1999.
• [18] J.H. He. Some asymptotic methods for strongly nonlinear equations. International Journal of Modern Physics B, 20(10):11411199, 2006.
• [19] L.B. Ibsen, A. Barari, and A. Kimiaeifar. Analysis of highly nonlinear oscillation systems using He's max-min method and comparison with homotopy analysis and energy balance methods. Sadhana, 35:116, 2010.
• [20] O. Kopmaz and Ö. Gündogdu. On the curvature of an euler-bernoulli beam. International Journal of Mechanical Engineering Education, 31:132142, 2003.
• [21] L. Meirovitch. Fundamentals of vibrations McGraw-Hill, 2001. International Edition.
• [22] H. Mirgolbabaei, A. Barari, L.B. Ibsen, and M.G. Sfahani. Numerical solution of boundary layer flow and convection heat transfer over a flat plate. Archives of Civil and Mechanical Engineering, 10:4151, 2010.
• [23] M. Momeni, N. Jamshidi, A. Barari, and D.D. Ganji. Application of He's energy balance method to Duffing harmonic oscillators. International Journal of Computer Mathematics, 88(1):135144, 2011.
• [24] 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:127135, 2010.
• [25] T. Pirbodaghi, M.T. Ahmadian, and M. Fesanghary. On the homotopy analysis method for non-linear vibration of beams. Mechanics Research Communications, 36:143148, 2009.
• [26] M.I. Qaisi. Application of the harmonic balance principle to the nonlinear free vibration of beams. Appl. Acoust, 40:141151, 1993.
• [27] S. S. Rao. Vibraton of Continuous Systems John Wiley & Sons, Inc., Hoboken, New Jersey, 2007.
• [28] B.S. Sarma and T.K. Varadan. Lagrange-type formulation for finite element analysis of nonlinear beam vibrations. J. Sound Vib, 86:6170, 1983.
• [29] M. G. Sfahani, S.S. Ganji, A. Barari, H. Mirgolbabaei, and G. Domairry. Analytical solutions to nonlinear conservative oscillator with fifth-order non-linearity. Earthquake Engineering and Engineering Vibration, 9(3):367374, 2010.
• [30] M.G. Sfahani, A. Barari, M. Omidvar, S.S. Ganji, and G. Domairry. Dynamic response of inextensible beams by improved Energy Balance Method. In Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, volume 225(1), pages 6673, 2011.
• [31] Y. Shi and C.A. Mei. Finite element time domain model formulation for large amplitude free vibrations of beams and plates. J. Sound Vib, 193:453465, 1996.
• [32] F. S. Tse, I. E. Morse, and R. T. Hinkle. Mechanical Vibrations. Theory and Applications Allyn and Bacon Inc, second edition, 1978.

Publication Dates

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

History

• Received
  25 Oct 2010
• Reviewed
  22 Feb 2011