Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 470927, 7 pages

http://dx.doi.org/10.1155/2013/470927

## Analytical Solution for Free Vibration Analysis of Beam on Elastic Foundation with Different Support Conditions

^{1}Department of Civil Engineering, Faculty of Engineering, Hacettepe University, 06800 Ankara, Turkey^{2}Department of Civil Engineering, Faculty of Engineering, Kocaeli University, 41380 Kocaeli, Turkey

Received 15 December 2012; Accepted 17 April 2013

Academic Editor: Mehmet Atay

Copyright © 2013 Baki Ozturk and Safa Bozkurt Coskun. 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.

#### Abstract

Analytical solutions for free vibration analyses of a beam on elastic foundation are obtained for different support conditions. The analytical solutions are applied on three different axially loaded cases, namely; (1) one end clamped, the other end simply supported; (2) both ends clamped, and (3) both ends simply supported cases. Analytical solutions and frequency factors are evaluated for different ratios of axial load, , acting on the beam to Euler buckling load, . The analytical solutions give results which are in excellent agreement with the variational iteration method (VIM) and homotopy perturbation method (HPM) results available in the literature for the particular problem considering all the cases provided in this study and the differential transform method (DTM) results available in the literature for the clamped-pinned case.

#### 1. Introduction

The free vibration equation of an axially loaded beam on elastic foundation is a fourth-order partial differential equation. For this particular engineering problem, the analytical solutions will be implemented in this study. Although the governing equation seems to be a linear one, finding the eigenvalues for the free vibration analysis is still challenging. One may not simply obtain the eigenvalues sequentially and their corresponding eigen vectors even with a software.

Free vibration equation of the beam on partially elastic foundation including only bending moment effect was analytically solved [1] while the eigenvalues for free vibration of beam-column systems on elastic foundation were obtained using a numerical approach [2]. The separation of variables technique was used to obtain the free vibration circular frequencies of piles partially embedded in soils [3]. In addition, differential transform method (DTM) has been proposed to solve eigenvalue problems for free and transverse vibration problems of a rotating twisted Timoshenko beam under axial loading [4, 5]. Furthermore, the DTM was also used to find the nondimensional natural frequencies of tapered cantilever Bernoulli-Euler beam [6]. Free vibration equations for one end clamped and other end simply supported beam on elastic foundation were solved by using the DTM for various axial loads acting on the beam [7].

Meanwhile, both the variational iteration method (VIM) and homotopy perturbation method (HPM) were used to solve the free vibration equations of beam on elastic foundation for support conditions of one end clamped, and other end simply supported, both ends clamped and both ends simply supported considering various case studies [8–11]. The beam on elastic foundation was investigated for these three different support conditions considering various values.

Recently, there have been also other studies which are helpful to better understand dynamic behavior of both infinite beams resting on elastic foundation [12] and tapered column with pinned ends embedded in Winkler-Pasternak elastic foundation [13]. In this study, the analytical solutions and analytical results for free vibration analysis of beam on elastic foundation are provided. The analytical results are in excellent agreement with the results of the particular problem solved using the VIM method and the HPM method available in the literature [8–10].

#### 2. Problem Formulation

Beam on elastic foundation and internal forces and deformations of differential segment of the beam having a length of are depicted in Figures 1(a) and 1(b), respectively. The elastic foundation is idealized by Winkler model. Hence, the relationship between displacement function of the beam on elastic foundation and the distributed force applied on the elastic foundation beneath the beam can be written by . In this equation , where is the modulus of subgrade reaction while is the beam width.

The equilibrium equations of the internal forces acting on differential beam segment, dimensionless parameter instead of position variable with and neglecting the second-order terms are used in order to define the motion equation of the beam on elastic foundation [3]:

In (1), is dimensionless displacement function of the beam considering axial and shear forces; is the time variable; is the phase angle; is the ratio of axial load acting on the beam to Euler buckling load; is the distributed mass of the beam; is the beam circular frequency; is the shape factor for the shape of the beam section considered; is the beam length; , , , and are cross-section area, shear modulus, elastic modulus, and moment of inertia of the beam; and , , respectively. If the axial and shear force effects are neglected, the dimensionless equation of motion for the beam on elastic foundation becomes as follows [14]:

In (2), is the dimensionless displacement function of the beam neglecting axial and shear forces. Division of both sides of (1) and (2) by gives the following equations:

#### 3. Analytical Solutions of the Problem

The equation of motion in (3) can be rearranged as follows: where is relative stiffness, , and is frequency factor, .

In order to simplify the symbolic representation of solution process, (5) is rewritten as where

Assuming a solution of the form of , characteristic equation for (6) becomes

Equation (8) is a quadratic equation in terms of and its solution is

By taking the square root of (9), roots of (8) would be obtained as follows:

Hence the solution finally takes the following form:

Coefficients , , , and would be obtained by inserting the boundary conditions for the beam. In this study three different support configurations for the beam will be considered as shown in Figure 2. These are SS for *Simply Supported*, CC for *Clamped*-*Clamped,* and CS for *Clamped*-*Simply Supported* beams. The boundary conditions related to these supports are given below.

For a simple support:

For a clamped support:

Inserting the boundary conditions given in (12)-(13) into (11) for SS, CC, and CS beams, a nonlinear system of equations in the following form would be produced: where and includes the coefficients of (11). Coefficient matrix is dependent on frequency factor and equating the determinant to zero gives the frequency factors for vibration modes of the corresponding beam:

Positive real roots of this equation are free vibration frequencies of the beams on elastic foundations shown in Figure 2.

#### 4. Numerical Study

Analytical solution for free vibration of a one end clamped, and one end pinned beam on elastic foundation was previously given by Çatal [7] with additional analyses by DTM. However, Çatal [7] gives analytical solutions for only CS beam, and in this study it is observed that these results fail to coincide with the exact results.

Previously, both the variational iteration method (VIM) and homotopy perturbation method (HPM) were used to solve the free vibration equations of beam on elastic foundation for support conditions of one end clamped, and other end simply supported, both ends clamped and both ends simply supported considering various case studies [8–10]. In this study, analytical solutions for three cases are provided which are shown in Figure 2. Each case was previously analyzed by the use of both VIM and HPM with the previously explained procedure [8–10]. Numerical values in this study are chosen as the same used in Çatal [7]. Hence, an IPB 500 steel profile resting on a Winkler foundation having a modulus of subgrade reaction of 50,000 kN/m^{2} is considered. Other numerical values are as follows: ; ; ; ; ; .

Frequency factors are calculated taking bending moment, shear, and axial effects into consideration for , , and due to circular frequencies of the beam.

Variation of frequency factor, is tabulated in Tables 1–3. Table 1 includes comparison of analytical results of this study for the CS beam with analytical results obtained by Çatal [7]. Tables 2 and 3 show the analytical results for both CC and SS beams which are not solved by Çatal [7].

The graphical presentations of analytical solutions are provided in Figures 3, 4, and 5.

#### 5. Conclusion

In this study, the analytical results for free vibration analysis of beam on elastic foundation are provided for three different axially loaded cases which are namely one end clamped, the other end simply supported (CS beam); both ends clamped (CC beam), and both ends simply supported (SS beam) cases. In addition, the analytical solution and the frequency factors evaluated for different ratios of axial load acting on the beam to Euler buckling load, are obtained. In Figures 3, 4, and 5, it is observed that for the same values of , the largest values of are obtained for both ends clamped (CC) beam while the lowest values of are obtained for both ends pinned (SS) beam.

For the CS Beam, it is observed that analytical solutions given by Çatal [7] fail to coincide with the exact results provided in this study. Previous studies [8–10] employing both the variational iteration method (VIM) and homotopy perturbation method (HPM) to solve the free vibration equations of beam on elastic foundation in this study converge to the same results obtained in this study. The reason for that is the governing equation is a linear equation in nature; and as expected analytical approximate methods converge to the exact solution. In the study by Çatal [7], the DTM converges the analytical results provided in the same study. Presumably, the difference is due to the fact that Çatal [7] used the length of the beam, , in three-digit accuracy while in this study the exact values of length, , are used in all calculations conducted.

Accordingly, this study provides the analytical results for free vibration analysis of beam on elastic foundation for three different axially loaded cases which are namely one end clamped, the other end simply supported (CS beam); both ends clamped (CC beam) and both ends simply supported (SS beam) cases.

#### References

- P. F. Doyle and M. N. Pavlovic, “Vibration of beams on partial elastic foundations,”
*Earthquake Engineering and Structural Dynamics*, vol. 10, no. 5, pp. 663–674, 1982. View at Google Scholar · View at Scopus - H. H. West and M. Mafi, “Eigenvalues for beam-columns on elastic supports,”
*Journal of Structural Engineering*, vol. 110, no. 6, pp. 1305–1320, 1984. View at Publisher · View at Google Scholar · View at Scopus - H. H. Çatal, “Free vibration of partially supported piles with the effects of bending moment, axial and shear force,”
*Engineering Structures*, vol. 24, no. 12, pp. 1615–1622, 2002. View at Publisher · View at Google Scholar · View at Scopus - C.-K. Chen and S.-H. Ho, “Application of differential transformation to eigenvalue problems,”
*Applied Mathematics and Computation*, vol. 79, no. 2-3, pp. 173–188, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. K. Chen and S. H. Ho, “Transverse vibration of a rotating twisted Timoshenko beams under axial loading using differential transform,”
*International Journal of Mechanical Sciences*, vol. 41, no. 11, pp. 1339–1356, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - O. Özdemir and M. O. Kaya, “Flapwise bending vibration analysis of a rotating tapered cantilever Bernoulli-Euler beam by differential transform method,”
*Journal of Sound and Vibration*, vol. 289, no. 1-2, pp. 413–420, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Çatal, “Solution of free vibration equations of beam on elastic soil by using differential transform method,”
*Applied Mathematical Modelling*, vol. 32, no. 9, pp. 1744–1757, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - B. Öztürk, “Free vibration analysis of beam on elastic foundation by the variational iteration method,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 10, no. 10, pp. 1255–1262, 2009. View at Google Scholar · View at Scopus - B. Öztürk and S. B. Coşkun, “The homotopy perturbation method for free vibration analysis of beam on elastic foundation Structural Engineering and Mechanics,”
*An International Journal*, vol. 37, no. 4, 2011. View at Google Scholar - B. Öztürk, S. B. Coşkun, M. Z. Koç, and M. T. Atay, “Homotopy perturbation method for free vibration analysis of beams on elastic foundations,”
*IOP Conference Series*, vol. 10, no. 1, 2010. View at Google Scholar - S. B. Coşkun, M. T. Atay, and B. Öztürk, “Transverse vibration analysis of euler-bernoulli beams using analytical approximate techniques,” in
*Advances in Vibration Analysis Research*, chapter 1, pp. 1–22, InTech, Vienna, Austria, 2011. View at Google Scholar - I. G. Raftoyiannis, T. P. Avraam, and G. T. Michaltsos, “Dynamic behavior of infinite beams resting on elastic foundation under the action of moving loads. Structural Engineering and Mechanics,”
*An International Journal*, vol. 35, no. 3, 2010. View at Google Scholar - Ö. Civalek and B. Öztürk, “Free vibration analysis of tapered beam-column with pinned ends embedded in Winkler-Pasternak elastic foundation,”
*Geomechanics and Engineering*, vol. 2, no. 1, pp. 45–56, 2010. View at Google Scholar · View at Scopus - J. Tuma and F. Cheng,
*Theory and Problems of Dynamic Structural Analysis*, Schaum’s Outline Series, McGraw-Hill, New York, NY, USA, 1983.