Numerical and Soft Computing Methods for Characteristic Value Problems of ODE and ODEs Systems
View this Special IssueResearch Article  Open Access
Utkan Mutman, "Free Vibration Analysis of an Euler Beam of Variable Width on the Winkler Foundation Using Homotopy Perturbation Method", Mathematical Problems in Engineering, vol. 2013, Article ID 721294, 9 pages, 2013. https://doi.org/10.1155/2013/721294
Free Vibration Analysis of an Euler Beam of Variable Width on the Winkler Foundation Using Homotopy Perturbation Method
Abstract
Homotopy Perturbation Method (HPM) is employed to investigate the vibration of an Euler beam resting on an elastic foundation. The beam is assumed to have variable stiffness along its length. HPM is an easytouse and very efficient technique for the solution of linear or nonlinear problems. HPM produces analytical approximate expression which is continuous in the solution domain. This work assures that HPM is a promising method for the vibration analysis of the variable stiffness Euler beams on elastic foundation. Different case problems have been solved by using the technique, and solutions have been compared with those available in the literature.
1. Introduction
In geotechnical engineering, it is widely seen in application that pipelines, shallow foundations, and piles are modeled as a beam in the analysis procedure. There are also various types of foundation models such as Winkler, Pasternak, and Vlasov that have been used in the analysis of structures on elastic foundations.
The most frequently used foundation model in the analysis of beam on elastic foundation problems is the Winkler foundation model. In the Winkler model, the soil is modeled as uniformly distributed, mutually independent, and linear elastic vertical springs which produce distributed reactions in the direction of the deflection of the beam. The winkler model requires a single parameter k, representing the soil properties. However, since the model does not take into account either continuity or cohesion of the soil, it may be considered as a rather crude representation of the soil.
As stated before pipelines, shallow foundations and piles can be modeled as a beam. Hence, produced results may find various application areas. There are also different beam types in theory. The wellknown is the EulerBernoulli beam which is suitable for slender beams. For moderately short and thick beams, the Timoshenko beam model has to be used in the analysis.
Vibration of a constant stiffness Euler beam on elastic foundation was studied previously by Balkaya et al. [1] and Ozturk and Coskun [2]. Balkaya et al. [1] used Differential Transform Method while Ozturk and Coskun [2] used HPM in their studies. Avramidis and Morfidis [3] studied bending of beams on threeparameter elastic foundation. De Rosa [4] analyzed free vibration of Timoshenko beams on twoparameter elastic foundation. Matsunaga [5] studied vibration and buckling of deep beam columns on twoparameter elastic foundations. ElMously [6] determined fundamental frequencies of the Timoshenko beams mounted on the Pasternak foundation. Chen [7, 8] analyzed vibration of beam resting on an elastic foundation by the differential quadrature element method (DQEM). Coşkun [9] investigated the response of a finite beam on a tensionless Pasternak foundation subjected to a harmonic load. Chen et al. [10] developed a mixed method for bending and free vibration of beams resting on a Pasternak elastic foundation. Maheshwari et al. [11] studied the response of beams on a tensionless extensible geosyntheticreinforced earth bed subjected to moving loads. Auciello and De Rosa [12] developed different approaches to the dynamic analysis of beams on soils subjected to subtangential forces.
2. The Equations of Motion and Boundary Conditions
An Euler beam resting on the Winkler foundation shown in Figure 1 is considered in this study. The equation of motion for this problem is given as follows: where is the spring constant, is deflection, is the mass density, is the crosssectional area, EI is the beam stiffness, and is the area moment of inertia about the neutral axis. The deflection is a function of both space and time, that is, in which space variable is measured along the length of the beam and represents any particular instant of time.
Due to the support conditions at both ends of the beam, different conditions have to be imposed to the obtained solution to determine unknowns included in final approximation produced by HPM. These conditions are given as follows.(a) For clampedclamped beam the end conditions are (b) For cantilever beam (clampedfree) the end conditions are (c) For simply supported beam (pinnedpinned) the end conditions are (d) For clampedsimply supported (pinned) beam the end conditions are
Now, free vibration analysis of the beams with variable stiffness resting on elastic foundations will be formulated.
A solution is assumed as the following form to formulate the analysis of the presented problem by the separation of variables where ω is the circular frequency for the vibration. Substituting (6) into (1), equations of motion become as follows:
This equation can be rearranged as where denotes total derivative with respect to . The governing equation is now rewritten in a nondimensional form. This procedure is provided from [1] in which a constant stiffness beam was analyzed. The notation is maintained in this study for the comparison purposes. The nondimensional parameters for the Eulerbeam on the Winkler foundation are defined as [1]
Using these parameters, nondimensional form of the equations and formulation procedures are explained in the following sections.
3. Homotopy Perturbation Method
In recent years, efforts towards application of analytical approximate solution techniques for nonlinear problems have increased. In these studies, He’s Homotopy Perturbation Method [13–18] can be considered as one of the most promising methods for nonlinear problems. The HPM provides an analytical approximate expression as the solution for the problems which are continuous in the solution domain. The technique is applied to an equation of the form with boundary conditions , , where is a linear operator, N is a nonlinear operator, B is a boundary operator, Γ is the boundary of the domain , and f(r) is a known analytic function. HPM defines a homotopy as which satisfies the following inequalities:
or where and are imbedding parameters, is an initial approximation which satisifies the boundary conditions. From (10)(11), we have
The changing process of from zero to unity is that of from to . In topology, this deformation and are called homotopic. The method expresses the solution of (10)(11) as a power series in as follows:
The approximate solution of , can be obtained as
The convergence of the series in (14) has been proven in [13–18].
4. HPM Formulation and Solution Procedure for Presented Problem
In this study, a variable stiffness is assumed for the beam considered. The variation is applied by changing the width of the beam along its length. Two different variations in stiffness are assumed. These are linear and exponential variations.
In the linear variation, beam width varies in a linear manner along the beam’s length which is shown in Figure 2.
The variable width is formulated as where the dimension of α is . For the exponential variation, beam width varies exponentially along the beam’s length as shown in Figure 3.
Variable width of the beam in this case is where α has the same dimension as in linear case. By the use of variable width, both crosssectional area and moment of inertia are variable,
In the following sections, formulation procedures are given separately for each case.
4.1. Formulation for Linear Variation
Employing (15) into (17), variable crosssectional areas and variable stiffness become where and are the crosssectional area and moment of inertia of the section at the origin, respectively. Inserting (18) into (8)
This equation can be rewritten as
Employing (6)
Equation (21) can be made nondimensional in view of (9) as follows: where
By the application of HPM, the following iteration algorithm is obtained:
4.2. Formulation for Exponential Variation
Employing (16) into (17), variable crosssectional areas and variable stiffness become
Employing (6)
Equation (21) can be made nondimensional in view of (9) as follows:
By the application of HPM, the following iteration algorithm is obtained:
4.3. Solution Procedure
An initial approximation may be chosen as a cubic polynomial with four unknown coefficients. There exist four boundary conditions, that is, two at each end of the column, in the presented problem. Hence, an initial approximation of the following choice may be employed:
In the computations, twenty iterations are conducted and four boundary conditions for each case are rewritten by using the final approximation of iteration. Each boundary condition produces an equation containing four unknowns spread from the initial approximation. These nondimensional boundary conditions are as follows.
Clampedclamped beam:
Clampedfree (cantilever) beam:
Pinnedpinned (simply supported) beam:
Clampedpinned beam:
Hence, four equations in four unknowns may be written with respect to the boundary conditions of the problem. These equations can be represented in matrix form as follows: where . For a nontrivial solution, determinant of coefficient matrix must be zero. Determinant of coefficient matrix yields a characteristic equation in terms of . Positive real roots of this equation are the normalized free vibration frequencies for the case considered.
5. Numerical Results
5.1. Constant Stiffness Case
As the first example, Euler beam of constant stiffness (i.e., is constant) with different boundary conditions is investigated. For the sake of comparison, all the values are set to unity such as ; hence, , according to previous studies [7]. Both algorithms given for linear and exponential variations lead to constant stiffness when .
In Table 1, first three normalized free vibration frequencies for simply supported (pinnedpinned) beam are compared with previous available results in the literature and the exact solution. Excellent agreement is observed with the exact solution.
The first five natural frequencies for clampedclamped beam and cantilever (clampedfree) beam are presented in Tables 2 and 3, respectively. They are compared with the available solutions.
Excellent agreement is observed with previous available results for both cantilever and clampedclamped beams. Clampedpinned beam was not included in the studies used for comparison. Hence, only the result obtained from HPM is tabulated for this case in Table 4.

As one can see, perfect agreement is obtained for constant stiffness case. This issue is mainly due to constant coefficient governing equation. In the following sections, variable stiffness cases are investigated.
5.2. Linearly Varying Stiffness
Linearly varying beam width results in a linearly varying flexural stiffness. The variation is based on parameter α. A number of case studies are conducted with respect to parameter α, and the results are given in Tables 5, 6, 7, and 8.




Variation of normalized free vibration frequencies () with respect to nondimensional variation parameter for each beam are also given in Figures 4, 5, 6, 7, and 8.
5.3. Exponentially Varying Stiffness
Exponentially varying beam width results in an exponentially varying flexural stiffness. The variation is again based on parameter α, and a number of case studies are also conducted with respect to parameter α. Results are given in Tables 9, 10, 11, and 12.




Variation of normalized free vibration frequencies () with respect to nondimensional variation parameter for each beam are also given in Figures 9, 10, 11, 12, and 13.
6. Conclusions
In this study, HPM is introduced for the free vibration analysis of variable stiffness Euler beams on elastic foundations. As a demonstration of application of the method, firstly the case studies for which previous results were available are chosen. In these studies, constant stiffness Euler beams were considered, and first analyses are conducted for constant stiffness case for comparison with the available results. HPM produced results which are in excellent agreement with the previously available solutions that encourage the application of the method for variable stiffness beams. To represent a variation in stiffness, a rectangular beam with varying width is considered, and two types of variation are taken into consideration. These are, namely, linear and exponential changes. Such varying widths lead to linearly and exponentially varying stiffnesses. The analyses are expanded for variable stiffness cases. HPM also produced reasonable results for the vibration of variable stiffness beams which show the efficiency of the method. For the variable stiffness problems, the governing equation is a differential equation with variable coefficients, and it is not easy to obtain analytical solutions for these types of problems. However, it is easy to put those variable parameters into the iteration algorithm of HPM, and the results can be obtained after performing some iterations with the method. The results obtained in this study point out that the proposed method is a powerful and reliable method in the analysis of the presented problem.
References
 M. Balkaya, M. O. Kaya, and A. Saǧlamer, “Analysis of the vibration of an elastic beam supported on elastic soil using the differential transform method,” Archive of Applied Mechanics, vol. 79, no. 2, pp. 135–146, 2009. View at: Publisher Site  Google Scholar
 B. Ozturk and S. B. Coskun, “The Homotopy Perturbation Method for free vibration analysis of beam on elastic foundation,” Structural Engineering and Mechanics, vol. 37, no. 4, pp. 415–425, 2011. View at: Google Scholar
 I. E. Avramidis and K. Morfidis, “Bending of beams on threeparameter elastic foundation,” International Journal of Solids and Structures, vol. 43, no. 2, pp. 357–375, 2006. View at: Publisher Site  Google Scholar
 M. A. De Rosa, “Free vibrations of Timoshenko beams on twoparameter elastic foundation,” Computers and Structures, vol. 57, no. 1, pp. 151–156, 1995. View at: Google Scholar
 H. Matsunaga, “VIbration and buckling of deep beamcolumns on twoparameter elastic foundations,” Journal of Sound and Vibration, vol. 228, no. 2, pp. 359–376, 1999. View at: Google Scholar
 M. ElMously, “Fundamental frequencies of Timoshenko beams mounted on Pasternak foundation,” Journal of Sound and Vibration, vol. 228, no. 2, pp. 452–457, 1999. View at: Google Scholar
 C. N. Chen, “Vibration of prismatic beam on an elastic foundation by the differential quadrature element method,” Computers and Structures, vol. 77, no. 1, pp. 1–9, 2000. View at: Publisher Site  Google Scholar
 C. N. Chen, “DQEM vibration analyses of nonprismatic shear deformable beams resting on elastic foundations,” Journal of Sound and Vibration, vol. 255, no. 5, pp. 989–999, 2002. View at: Publisher Site  Google Scholar
 I. Coşkun, “The response of a finite beam on a tensionless Pasternak foundation subjected to a harmonic load,” European Journal of Mechanics, A/Solids, vol. 22, no. 1, pp. 151–161, 2003. View at: Publisher Site  Google Scholar
 W. Q. Chen, C. F. Lü, and Z. G. Bian, “A mixed method for bending and free vibration of beams resting on a Pasternak elastic foundation,” Applied Mathematical Modelling, vol. 28, no. 10, pp. 877–890, 2004. View at: Publisher Site  Google Scholar
 P. Maheshwari, S. Chandra, and P. K. Basudhar, “Response of beams on a tensionless extensible geosyntheticreinforced earth bed subjected to moving loads,” Computers and Geotechnics, vol. 31, no. 7, pp. 537–548, 2004. View at: Publisher Site  Google Scholar
 N. M. Auciello and M. A. De Rosa, “Two approaches to the dynamic analysis of foundation beams subjected to subtangential forces,” Computers and Structures, vol. 82, no. 6, pp. 519–524, 2004. View at: Publisher Site  Google Scholar
 J. H. He, “Coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 J. H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004. View at: Publisher Site  Google Scholar
 J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons & Fractals, vol. 26, no. 3, pp. 695–700, 2005. View at: Publisher Site  Google Scholar
 J. H. He, “Asymptotology by homotopy perturbation method,” Applied Mathematics and Computation, vol. 156, no. 3, pp. 591–596, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 12, pp. 87–88, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, “Limit cycle and bifurcation of nonlinear problems,” Chaos, Solitons & Fractals, vol. 26, no. 3, pp. 827–833, 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Utkan Mutman. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.