SciELO - Scientific Electronic Library Online

vol.14 número9Dynamic Compressive Strength and Failure of Natural Lake Ice Under Moderate Strain Rates at Near Melting Point TemperatureNumerical Investigation of Structural Response of Corrugated Blast Wall Depending on Blast Load Pulse Shapes índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados




Links relacionados


Latin American Journal of Solids and Structures

versão impressa ISSN 1679-7817versão On-line ISSN 1679-7825

Lat. Am. j. solids struct. vol.14 no.9 Rio de Janeiro set. 2017 


Nonlinear Vibration Analysis of Euler-Bernoulli Beams by Using Continuous Galerkin-Petrov Time-Discretization Method

M. Sabeel Khana  * 

H. Kaneeza 

a Department of Applied Mathematics, Institute of Space Technology, 44000 Islamabad, Pakistan


In this paper, we present a new numerical method for nonlinear vibrational analysis of Euler-Bernoulli beams. Our approach is based on the continuous Galerkin-Petrov time discretization method. The Euler-Bernoulli beam equation which governs its vibrations is transformed into set of ordinary differential equations and the presented method is employed in order to investigate the vibrational response. A comparison is made between present method and different other methods available in literature. It is observed that the obtained results are in strong agreement with other results in literature. We conclude that the present method has a great potential to deal with nonlinear vibration analysis problems of beams and related structures like rods and shafts.

Keywords: Nonlinear vibration; Euler-Bernoulli beams; Time discretization; Numerical method


In the design and fabrication of many engineering structures and machines, vibrational analysis and dynamical response of structures like beams (or rods or shafts) are important factor to investigate in order to increase the performance of these structures. For the analysis of complex engineering structures like bridges, tall buildings, vehicle guide-ways, huge cranes, turbines and compressor blades, beams can be used as simple model. The dynamical response of beams is governed by linear as well as nonlinear differential equations both in space and time. To study their behavior it is therefore important to design methods for the numerical solutions of these partial differential equations. In this respect, recently some approximate analytical techniques (see for instance, Barari et al., 2010; Bayat et al., 2010; Baghari et al., 2014; Jafari et al., 2014; Bayat et al., 2011;, Ganji, 2012; He, 1999; He, 2006; Khan et al., 2012; Liu and Gurram, 2009; Mirgolbabaei et al., 2010; Sfahani et al., 2011) as well as numerical methods by Lai et al. (2008) (see also, Boukhalfa et al., 2010; Da Silva et al., 2009; Ganji et al., 2011) and references therein, have been proposed to investigate the nonlinear vibrations of beams for designing and fabrication purpose.

The theory of Euler-Bernoulli beams which is based on the assumption that beam curvature is related to its bending moment provides a good explanation of long isotropic bar’s bending behavior. Several methods have been proposed for obtaining solution to the nonlinear equations of motions governing the Euler-Bernoulli beam’s response. Nonlinear vibrations of the Euler-Bernoulli beam has been investigated by Barari et al. (2011) using the Variational iteration and parameterized perturbation methods. Nikkar et al. (2014) studied the nonlinear vibration response of Euler-Bernoulli beam by approximate analytical techniques where they utilized Variational He’s approach and Laplace iteration technique in order to solve the respective nonlinear governing equations. Pakar and Bayat (2013) applied a max-min approach, also by He, to investigate the nonlinear vibrational response of Euler-Bernoulli beams subjected to an axial load. Nonlinear behavior of Euler-Bernoulli beam with different end conditions has also been studied by Rafieipour et al. (2014) using Laplace iteration method. Lindstedt-Poincare techniques for non linear vibrational analysis have been employed to double-clamped and simply-supported beams subjected to an axial load by Ahmadian et al. (2009). Javanmard et al. (2013) used He’s Energy method in order to get solution of nonlinear vibration problem of Euler-Bernoulli beam subjected to axial loads. A similar study has been performed earlier by Pirbodaghi et al. (2009) but with homotopy analysis method. The nonlinear equation of motion for large amplitude free vibrational analysis of Euler-Bernoulli beam resting on variable elastic foundation was solved by Mirzabeigy and Madoliat (2016) with the application of second order homotopy perturbation method. Moreover, the nonlinear vibrational response of Euler-Bernoulli beams with geometric nonlinarity and subjected to axial loads was computed by Johnson et al. S. (2014) where they used a differential transform and two auxiliary parameter based homotopy analysis techniques.

Here, in this article, we present a new approach to solve equations of motions governing the dynamics of Euler-Bernoulli clamped-clamped-beam which is fixed from one end. This article is organized as follows. In Section II, a mathematical framework governing the nonlinear vibrations of Euler-Bernolulli beam is presented. In Section III, continuous Galerkin-Petrov time-discretization method is developed. In Section IV, the developed method afterwards applied to the governing equations of Euler-Bernoulli beam presented in Section II. In Section V, simulations data and numerical results are discussed. In last Section, conclusions are drawn based on the obtained simulation results in Section V.


Let us consider a straight beam composed of homogeneous material having length L, placed on an elastic foundation with modulus of elasticity as E. The beam is experiencing an axial force F as shown in Figure 1, below. Let ‘A’ be the cross-sectional area of beam which is assumed to be uniform throughout its length. Further, suppose that there is no in plane deformation i.e. the planes of the cross section remains plane after deformation. Now by ignoring the transverse normal and shear strains one can write the equation of motion of the Euler-Bernoulli beam including the mid-plane stretching effect as follows

EI4Yx4+μ2Yt2+F2Yx2+dYt+KYEA2L2Yx20L(Yx)2dx=q(x,t), (1)

where, d is the coefficient of viscous damping, K is the stiffness of foundation and q(x, t) is the transverse directional load. In the absence of non-conservative forces above equation (1) reduced to

EI4Yx4+μ2Yt2+F2Yx2+KYEA2L2Yx20L(Yx)2dx=0. (2)

Figure 1 A schematic representation of an Euler-Bernoulli beam fixed at one end and subjected to axial load. 

In order to undimensionalize the above equation we choose the following variables

X˜=XL,Y˜=YrG,t˜=tEIμL4,F˜=FL2EIandK˜=KL4EI, (3)

where, rG is the gyration radius of the cross-section of the beam. Using equation (3) and omitting the tilde sign from the variables equations (2) thus can be written as follows

4Yx4+2Yt2+F2Yx2+KY122Yx201(Yx)2dx=0. (4)

Now, by assuming a product solution Y (x, t) = φ(x) η(t) where φ(x) is the beam’s eigen mode, and by applying the Galerkin method, an equation governing the beam dynamics is obtained as

d2η(t)dt2+εη(t)+λη3(t)=0 (5)

Where ε and λ can be determined by the relations

ε=K+01ϕ(iv)ϕ dx01ϕ2dx+F01ϕϕ dx01ϕ2dx,andλ=01(ϕ01ϕ2dx)ϕdx2.01ϕ2dx.

Moreover, in order to complete the above problem description, the beam is subjected to following initial conditions

η(0)=A,dηdt(0)=0, (6)

where A is the amplitude of oscillation. Equation (5) along with initial conditions in Equation (6) governs the nonlinear vibrational response of Euler-Bernoulli beams.


In this section, we present a time discretization method to solve nonlinear equations governing the dynamics of Euler-Bernoulli beams. It is easy to see that Equation (5) can be transformed into a system of nonlinear differential equations of the following form

ddt{U(t)}=f(t,U(t)). (7)

Where, U(t) = (η, ξ)T and f denotes the vector of unknown-state variables and vector comprising non-linear functions in both time and space, respectively. Here, ξ is the first order time derivative of η. In order to develop a numerical method to solve such dynamical systems, here, we use the concept of a Galerkin-type formulation and time-discretization. Let S be a space of all possible solutions to set of nonlinear differential equations in Equation (7). We search a vector of unknown states U: [0, τ] → S such that

ddt{U(t)}=f(t,U(t)),t(0,τ),whereU(0)=U0. (8)

By choosing ℑ as a test space the problem in Equation (8) can be formulated as follows: Find US such that

0τVT(t)ddtU(t)dt=0τVT(t)f(t,U(t)) dt, V .WithU(0)=U0, (9)

Equation (9) is the weak form of problem in Equation (8). Now, to solve it over the time interval I = (0, τ) let us discretized the interval I into further sub-intervals with each sub-interval In = ( tn -1, tn ) where n belongs to the set {1, ...,N} Now, each interval is having the following property

n=1N(tn1, tn)= (0,τ).

Moreover, it also satisfy

(tn1, tn) (tm1, tm)={0},for alln m n,m=1,...,N.

In order to compute the vector of state variables U(t) at time interval In , we choose Ψ n,j (t) ∈ P k ( In ) as a set of basis functions. Here, P k ( In ) represents the k-th degree polynomials space defined on interval In by

Ρ(Ιn,SH):={U (t): ΙnSH;U(t) = j=0kujtj,for allt Ιn, uj SH}. (10)

The vector of state variables U(t) in Equation (5) is thus approximated by

U(t) Uσ(t)=j=0kunj ψn,j(t), for all tIn (11)

Where, unj are vector elements in Hilbert space SH Global vector of state variables U σ(t): In SH In is in the discrete-solution space Sσk S where,

Sσk:= {U(t) C(Ιn,SH):U(t)|Ιn Ρk (Ι¯n ,SH) ;for alln =1,...,N },

with Īn = In ⋃ { tn -1, tn }. The symbol σ above represents a discretization parameter. The test function V(t) is taken from discrete test space σk where,

σk:={V(t) L2 (Ιn,SH):V(t)|Ιn Ρk-1(Ιn,SH); n =1,...,N}.

By denoting V σ(t) as a discrete approximation of the test function V(t), then the weak problem in Equation (9) takes the following Variational form which states:

Find Uσ(t) Uo +Sσ,0k such that

0τVσT(t) ddtUσ(t) dt = 0τVσT(t) f(t,Uσ(t)) dt;Vσσk. (12)

Above sσ,0k:= SσkSo denotes the subspace of Sσk with zero as initial condition. Define a polynomial φn,i P k-1( Īn ) on the time interval I zero everywhere in the set I, except at the time interval Īn = [ tn -1, tn ], and v SH as an arbitrary scalar function then the variational problem on time interval In in Equation (12) reads

Ιn(vT(t) ddtUσ(t)) ϕn,i(t) dt = Ιn(vT(t) f (t,Uσ(t))) ϕn.i(t) dt;v SH. (13)

Or, equivalently, we can write by using Equation (11)

In(vT(t)j=0kunj ddtψn,j(t))ϕn,i(t) dt=In(vT(t) f (t,j=0kunj ψn,j (t)))ϕn.i(t) dt;v SH. (14)

To define Basis functions in Equation (14), a mapping Mn:IΙ¯n is used which maps the time interval I=[1,1] from natural time interval I to physical time interval Īn . Now, for every tI there exists t Īn such that

t = Mn(t^):= tn-1+ tn2 + σn2t^ Ι¯n;n {1,...,N}.

The basis functions defined on physical time interval Īn can now be calculated in natural time interval as

ψn,j(t):= ψj (Mn1(t));j =0,...,k,φn.i(t):= φ^i(Mn1(t));i =0,...k-1.

Where, ψj Ρk(Ι^), and φ^i Ρk-1(Ι^) are satisfying the property

ψj(1) = δ0,j,ψ^j(1) = δk,j, (15)

δk,j Denote the Kronecker’s delta. To get a second-order approximation of discrete vector of state variables U σ(t) we require to find the coefficients unj ; j ∈{0, 1, 2}. By using property of basis function from Equation (15) and with the application of initial condition the coefficients un0 can be determined from the following

un0 = { U0, n=1un-1k,n>1. (16)

Since, it is easy to deal with natural-time interval Î, thus Equation (14) is reconstructed into

j=0kαi,j(vTunj) = σn2Ι^vTf(Mn(t^),j=0kunj Ψ^j(t^)) ϕ^i (t^) dt^;v SH (17)

where, i = 0, ..., k - 1 and

αi,j:= Ι^ ddt^ψ^j (t^) ϕ^i (t^) dt^. (18)

Now, with the application of (k + 1)-point Gauss-Lobatto formula in Equation (18) we write

j=0kαi,j (vTunj) l=0kw^l {vTf(Mn (t^l),j=0kunj ψ^j(t^l))ϕ^i(t^l);v SH (19)

Where, Ŵ l represents the weights t^0 =-1 ,t^k =1 , and t^l Ι^;l=1,...,k1 are roots of the Legendre polynomial P k-1(t). Using the abbreviation

tn,l:=gn(t^l),βi,l:= w^l ϕ^ (t^l)andηj.l:= ψ^j (t^l). (20)

Hence, we finally arrived at the following system of coupled equations

j=0kαi,j (vTunj) = σn2 l=0kβi,l (vTf (tn,l ,j=0kunj ηj,i));v SH,i =0,...,k-1 andj=0,...,k, (21)

with initial conditions as described in Equation (16).

3.1 Computation of Parameters α i,j , β i,j and η j,l

To compute unknown parameters in Equation (21) let us select the following set of test functions

ψ^i (t^j) = δi,j,φ^0 (t^j) = 34 t^j,φ1 (t^j) =1for alli,j {0,1,2},

Where, each ψ^j(t^)Ρk(Ι^),φ^i (t^)Ρk-1(Ι^), and t^0 =-1 ,t^1 =0 , and t^2=1 . By this choice of test functions αi,j , βi,j and ηj,l are calculated and are given as under

{α0,0=12,α0,1=1=α1,2,α0,2=α0,0,α1,0=α0,1,α1,1=0,β0,0=14,β0,1=0,β0,2=β0,0,β1,1=3β0,0,β1,0=β1,2=14β1,1ηi,j = δi,j ,for alli,j {0,1,2}.,


The method presented in Section 3 is applied to motion equation in Equation (5) and the obtained discrete set of equations are solved using numerical algorithm outlined as under

Numerical Algorithm:

Step 1: Initialize the known parameters ε and λ

Step 2: Choose a time step size σ n

Step 3: Set the initial conditions at time tn as

ηnm1=A, ξnm1=0

Step 4: First solve the following equations for unknowns ηnm,ξnm

ηnm=12ηnm1+12ηnm+1+σn8 (ξnm1ξnm+1)ξnm=12ξnm1+12ξnm+1+σn8 [ε(ηnm1ηnm+1)λ{ (ηnm1)3(ηnm+1)3}]

Step 5: Using the computed values for ηnm and ξnm from Step 4, Solve the following nonlinear equations for ηnm+1 and ξnm+1.

ηnm+1=ηnm1+σn6 (ξnm1+4ξnm+ξnm+1)ξnm+1=ξnm1σn6{ε(ηnm1+4ηnm+ηnm+1)+λ((ηnm1)3+4(ηnm)3+(ηnm+1)3)}

Step 6: Use the computed values of ηnm+1 and ξnm+1 from Step 5 to get the numerical values of ηnm and ξnm from Step 4 by back substitution.

Step 7: Compute the solution at time step tn by equation (11)

Step 8: Update the solution for the next time iteration by


Step 9: Update the time step tn tn + σ n

Step 10: Go to Step 4.


In this Section, nonlinear vibrational response of a clamped-clamped Euler-Bernoulli beam fixed at one end is presented by a new numerical method. The results computed by the presented method are compared with the other existing methods in literature (for instance, with the results of Barai et al. 2011, Nikkar et al. 2014 and Johnson et al. 2014). For the numerical computations a time step size of 0.05 is chosen for the discretization of time domain by the present method. The non-dimensional deflection η(t) is calculated by the presented method at the center of the beam and is shown in Figure 2, with the variation of non-dimensional maximum amplitude of oscillation over a wide range of time.

Figure 2 Non-dimensional deflection with varying values of amplitude A and with parametric values ( = λ = 1.0. 

In Table 1 and Table 2, a comparison is made between the presented method and four different methods, with semi-analytical approach by Barari et al. (2011), with a semi-analytical method by Nikkar et al. (2014), with semi-analytical approach by Johnson et al. (2014) and with numerical method RK-4. It is seen that the computed values of deflection at the center of the beam for different values of maximum amplitude of oscillation are in strong agreement with the results by these methods. In Table 3, relative errors are provided which are computed by the formula

(Solution by Present Method - Solution by Existing Method)Solution by Existing Method×100%,

Where in, Error-1 is calculated by taking RK-4 solution as an existing solution. In Error-II the existing solution is taken from the semi-analytical approach by Barari et al. (2011). Error-III is obtained by taking the existing solution form the method by Nikkar et al. (2014). In all of these calculations the errors are shown for three different amplitude of oscillation of beam over a wide range of time. From Table 3, it can be clearly seen that the results by present method converges closely to other results in literature and are in strong agreement.

Table 1 Comparison of time marching solution of equation (5) with parametric values ( = λ = 1.0. 

A RK-4 Barari et al. (2011) Nikkar et al. (2014) Present
0.01 0.008775827082332 0.008775710513259 0.008775710509477 0.008653113482869
0.1 0.087759719142174 0.087643206097855 0.087642827100616 0.086406655501930
0.2 0.175528204796203 0.174597455961528 0.174585253116176 0.172060842262657
0.3 0.263314167649526 0.260180487515907 0.260086892907199 0.263974531478741
0.4 0.351126207698248 0.343723299721559 0.343323578255470 0.349072483074928
0.5 0.438972761263577 0.424576467617775 0.423336627126820 0.417180214867134
0.6 0.526862050755839 0.502116086386647 0.498973411668122 0.492625565690737
0.7 0.614802037962082 0.575749233260282 0.568819355464365 0.563875657455006
0.8 0.702800381522887 0.644918846456226 0.631123805468134 0.630358772839251
0.9 0.790864399137227 0.709107932868351 0.683724780583061 0.691556439573670
1.0 0.879001034896931 0.767843030699686 0.723980197798699 0.766656852847480

Table 2 Differences between the time marching solution of equation (5) by presented method and different other schemes in literature. The parameters are chosen to be ( = λ = 1.0. 

A RK-4 present | Barari et al. (2011 ) present | Nikkar et al. (2014) present | Johnson et al. (2014) present |
0.01 0.000122713599463 0.0001225970303900 0.00012259702660799 0.0003821333327880
0.1 0.001353063640244 0.0012365505959250 0.00123617159868600 0.0019824623866750
0.2 0.003467362533546 0.0025366136988710 0.00252441085351901 0.0008865790820170
0.3 0.000660363829215 0.0037940439628340 0.00388763857154201 0.0052723348189849
0.4 0.002053724623320 0.0053491833533690 0.00574890481945800 0.0086853097383650
0.5 0.021792546396443 0.0073962527506410 0.00615641225968599 0.0087110167343409
0.6 0.034236485065102 0.0094905206959099 0.00634784597738497 0.0002395324508129
0.7 0.050926380507076 0.0118735758052760 0.00494369800935901 0.0063966672702809
0.8 0.072441608683636 0.0145600736169750 0.00076503262888294 0.0161331435918892
0.9 0.099307959563557 0.0175514932946811 0.00783165899060900 0.0078649950295989
1.0 0.112344182049451 0.0011861778522060 0.04267665504878090 0.0153845773856570

Table 3 Relative errors between the time marching solution of equation (5) by presented method with different other schemes in literature. The parameters are chosen to be ( = λ = 1.0. 

A time Error-I Error-II Error-III
0.01 5 8.39406344436937 1.43006947021705 1.43006930029566
10 1.64295821637604 1.65335398475738 0.46126609193088
25 -0.40714155296337 -0.30013048066676 -0.30013048391456
50 0.34564381067195 0.70609180050204 0.70609178786218
100 0.89942811653737 1.32438720734438 1.32438716831965
0.50 5 -0.01611628030638 -2.71860640853757 -3.34003971902284
10 -2.19636228576523 0.94834463220716 -0.15087284130191
25 1.33999271743426 0.18473004818300 -0.63051410694749
50 4.46407178480245 -0.21897390156171 -1.08460097418586
100 -5.91449682727209 -0.55016306148146 -0.57899624097327
1.00 5 -5.825758083900410 0.33425814206559 -1.25501908713661
10 -0.411754581578781 -1.42949464499397 -0.67810815791013
25 -0.614810964667922 3.07084443826624 0.52366254911929
50 0.200159632081500 -0.86954946868848 0.65475889510081
100 1.993160046364360 -0.27554338389560 0.96750105622486

In Figure 3, time marching deflection response over a wide range of time is shown where the result obtained from present method is graphically compared with the result computed from semi-analytical approach by Nikkar et al. (2014). In Figure 4, the time marching response of deflection by present method is compared with the solution of deflection response of center of the beam by Barari et al. (2011). It can be observed from these two Figures that the center of the beam shows same deflection behavior for a large time period as depicted by the semi-analytical approaches in literature. A parametric study is carried out to observe time marching response of the central deflection of beam. In Figure 5 and Figure 6, the central deflection of beam is observed with varying the parameters ( and λ The initial maximum amplitude of vibration was chosen to be 1.0 and a time step size of 0.005 for the calculation by the presented method. By varying the parameters ( and λ an increase in the maximum amplitude of vibration is seen as the time marches. In Table 4, a comparison is made between the accumulated errors by classical numerical scheme RK-4 and present method. Also the accuracy of solution by both of the numerical methods is shown for two different time step sizes. In all the computations of accumulated errors in Table 4 the reference solution by Johnson et al. (2014) is taken into consideration. It is observed that our present numerical method is more accurate in predicting the solution over the classical RK-4 method. This is highly due to the fact that present method does not involve iterative time integration and therefore have less accumulated errors as shown in Table 4. Moreover, it is observed that that the accumulated error decreases with an increase in the time step size for both the classical RK-4 and present method.

Figure 3 A comparison of time history of central deflection curve of beam. The time step used for this simulation is 0.05 and the parameters ( = λ = A - 1.0. 

Figure 4 A comparison of time history of central deflection curve of beam. The time step used for this simulation is 0.05 and the parameters ( = λ = A - 1.0. 

Figure 5 Time-history of central-deflection of beam with varying ( The parametric value of λ is chosen to be 1.0. 

Figure 6 Time-history of central-deflection of beam with varying λ The parametric value of ( is chosen to be 1.0. 

Table 4 A comparison of accumulated errors and accuracy for RK-4 and present method. 

Time- step size Number of time steps Accumulated error by RK-4 Accuracy (in percent) by RK-4 Accumulated error by present method Accuracy (in percent) by present method
0.05 5 7.1629284856E-04 7.838643342877E+01 4.7825078312E-04 8.556916328361E+01
10 4.1601192713E-04 9.991054140377E+01 2.8166615852E-04 9.996668334601E+01
25 8.9254489341E-05 9.910446043034E+01 6.3862965951E-05 9.935922760337E+01
50 1.8633928689E-04 9.809327174198E+01 1.4120960638E-04 9.855506398418E+01
100 3.6919714161E-04 9.584366665341E+01 2.9192878467E-04 9.671353538314E+01
250 7.2383261971E-04 7.572278868091E+01 5.9430126403E-04 8.006724623713E+01
500 3.7148542418E-04 9.975956260686E+01 1.8132325720E-04 1.000734824724E+02
0.005 5 5.04555298800E-05 9.82223873076E+01 2.39457037950E-05 9.91563623037E+01
10 2.64336613610E-05 1.00163169326E+02 1.36876932130E-05 1.00315112462E+02
25 4.97951263280E-05 1.00561378909E+02 2.28701839830E-05 1.01222287224E+02
50 1.70105932140E-05 9.98284057000E+01 3.19282474000E-06 9.99677923915E+01
100 3.59407763030E-05 9.96277324234E+01 6.42406862100E-06 9.99334607456E+01
250 1.61756733776E-04 9.35316240962E+01 2.41848699480E-05 9.90328882986E+01
500 4.74408476010E-05 9.94577701222E+01 1.21912724690E-05 1.00139341359E+02


A new numerical method is presented for nonlinear vibrational analysis of Euler-Bernoulli beams. The non-dimensional central deflection response of a clamped-clamped-Euler-Bernoulli beam (fixed from one end side) is investigated under different initial vibration amplitudes using this numerical technique. A comparison of computed results by present method is made with the results by other semi-analytical schemes in literature and the numerical values of central deflection of beam are tabulated. Relative errors between the present method and other semi-analytical approaches in literature have been shown to observe the convergence of solution obtained by present method. Moreover, a parametric study is carried out to investigate effect of different parameters on beam’s central deflection. It is worth mentioning that the obtained results by presented method are in strong agreement with results by semi-analytical approaches (for instance, Barari et al. 2011, Nikkar et al. 2014 and Johnson et al. 2014) in literature and also with classical numerical scheme RK-4. It is worthy of noting that the present scheme does not involve iterative time integration like other methods in literature and thus no associated accumulative errors. Therefore, it is capable to provide accurate results with high accuracy and convergence characteristics as can be seen by the presented results and thus can be utilized to compute dynamical response of structures like rods and shafts. We conclude that the present method has great potential in dealing with nonlinear vibrational response of beams and similar structures like rods and shafts.


The authors are grateful to the anonymous reviewers for their careful reading and constructive suggestions in order to improve the quality of this research article.


Ahmadian, M.T., Mojahedi, M., Moeenfard, H. (2009). Free vibration analysis of a nonlinear beam using homotopy and modified Lindstedt-Poincare Methods. Journal of Solid Mechanics. 29-36. [ Links ]

Bagheri, S., Nikkar, A., Ghaffarzadeh, H. (2014). Study of nonlinear vibration of Euler-Bernoulli beams by using analytical approximate techniques. Lat. Am. J. Solids Struct. 11: 157-168. [ Links ]

Barari, A., Ganjavi, B., Ghanbari, J M., Domairry, G. (2010). 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. [ Links ]

Barari, A., Kaliji, H.D., Ghadimi, M., Domairry, G. (2011). Non-linear vibration of Euler-Bernoulli beams. Lat. Am. J. Solids Struct. 8: 139-148. [ Links ]

Barari, B., Kimiaeifar, A., Domairry, G., Moghimi, M. (2010). Analytical evaluation of beam deformation problem using approximate methods. Songklanakarin Journal of Science and Technology, 32(3):207-326. [ Links ]

Bayat, M., Pakar, I., Bayat, M., (2011). Analytical study on the vibration frequencies of tapered beams, Lat. Am. J. Solids Struct. 8(2):149-162. [ Links ]

Bayat, M., Shahidi, M., Barari, A., Domairry, G. (2010). On the approximate analysis of nonlinear behavior of structure under harmonic loading. International Journal of Physical Sciences, 5(7):1074-1080. [ Links ]

Boukhalfa, A., Hadjoui, A., (2010). Free vibration analysis of an embarked rotating composite shaft using the hp-version of the FEM, Lat. Am. J. Solids Struct. 7:105-141. [ Links ]

Da Silva, J. C. R. Á., Beck, A. T., da Rosa, E., (2009). Solution of the stochastic beam bending problem by Galerkin method and the Askey-Wiener scheme, Lat. Am. J. Solids Struct. 6(1):51-72. [ Links ]

Ganji, D. D. (2012). A semi-Analytical technique for non-linear settling particle equation of Motion, Journal of Hydro-environment Research 6(4): 323-327. [ Links ]

Ganji, S.S., Barari, B., Ganji, D.D., (2011). Approximate analyses of two mass-spring systems and buckling of a column. Computers & Mathematics with Applications, 61(4):1088-1095. [ Links ]

He, J. H. (2006). Some asymptotic methods for strongly nonlinear equations, International Journal of Modern Physics B 20(10):1141-1199. [ Links ]

He., J.H. (1999). Variational iteration method - a kind of non-linear analytical technique: Some examples. International Journal of Nonlinear Mechanics, 34:699-708. [ Links ]

Jafari, S.S., Rashidi, M.M., Johnson, S. (2014). Analytical solution of the nonlinear vibration of euler-bernoulli beams via the homotopy analysis method and differential transform method. Int. J. of Appl. Math and Mech. 10(9): 96-110. [ Links ]

Javanmard, M., Bayat, M., Ardakani, A. (2013). Nonlinear vibration of Euler-Bernoulli beams resting on linear elastic foundation. Steel Comp. Struct. 15(4): [ Links ]

Jian-She Peng, J., Yan Liu, Y., Yang, J. (2010). A semi-analytical method for nonlinear vibration of Euler-Bernoulli beams with general boundary conditions. Math. Probl. Eng. [ Links ]

Khan, Y., Taghipour, R., Fallahian, M., Nikkar, A., (2012). A new approach to modified regularized long wave equation, Neural Computing and Applications doi:10.1007/s00521-012-1077-0 [ Links ]

Lai, H.Y., Hsu, J.C., Chen, C.K., (2008). An innovative eigenvalue problem solver for free vibration of Euler-Bernoulli beam by using the Adomian Decomposition Method, Computers and Mathematics with Applications, 56: 3204-3220. [ Links ]

Liu, Y., Gurram, C. S., (2009). The use of He’s Variational iteration method for obtaining the free vibration of an Euler-Beam beam, Mathematical and Computer Modeling 50:1545-1552. [ Links ]

Mirgolbabaei, H., Barari, A., Ibsen, L.B., Sfahani, M.G. (2010). Numerical solution of boundary layer flow and convection heat transfer over a flat plate. Archives of Civil and Mechanical Engineering, 10:41-51. [ Links ]

Mirzabeigy, A., Madoliat, R. (2016). Large amplitude free vibration of axially loaded beams resting on variable elastic foundation. Alexandria Engineering Journal 55, 1107-1114. [ Links ]

Pakar, I., Bayat, M. (2013). An analytical study of nonlinear vibrations of buckled Euler Bernoulli beams, Acta Physica Polonica A 123: 11-15. [ Links ]

Rafieipour, H., Tabatabaei, S.M., Abbaspour, M. (2014). A novel approximate analytical method for nonlinear vibration analysis of Euler-Bernoulli and Rayleigh Beams on the nonlinear elastic foundation. Arab J. Sci. Eng. 39: 3279-3287. [ Links ]

Sfahani, M.G., Barari, A., Omidvar, M., Ganji, S.S., Domairry, G. (2011). Dynamic response of inextensible beams by improved Energy Balance Method. In Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamic, 225(1): 66-73. [ Links ]

Received: August 29, 2016; Revised: January 20, 2017; Accepted: June 22, 2017

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License