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

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, N, acting on the beam to Euler buckling load, N r . 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.


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][9][10][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][9][10].

Problem Formulation
Beam on elastic foundation and internal forces and deformations of differential segment of the beam having a length of dx 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 0 ≤  ≤ / and neglecting the second-order terms are used in order to define the motion equation of the beam on elastic foundation [3]: In ( 1),  1 (z) is dimensionless displacement function of the beam considering axial and shear forces; t is the time variable;  is the phase angle;   =  2 /( 2 ) 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; L is the beam length; A, G, E, and I are cross-section area, shear modulus, elastic modulus, and moment of inertia of the beam; and 4 , 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),  2 () is the dimensionless displacement function of the beam neglecting axial and shear forces.Division of both sides of ( 1) and ( 2) by sin(+) gives the following equations:

Analytical Solutions of the Problem
The equation of motion in (3) can be rearranged as follows: where  is relative stiffness,  =    4 /, and  is frequency factor,  = 4 √  2  4 /.
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  2 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  1 ,  2 ,  3 , and  4 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 {}  = { 1  2  3  4 } 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.

Numerical Study
Analytical solution for free vibration of a one end clamped, and one end pinned beam on elastic foundation was previously given by C ¸atal [7] with additional analyses by DTM.However, C ¸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][9][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][9][10].Numerical values in this study are chosen as the same used in C ¸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:  = 107.2* 10 −5 m 4 ;  = 2.39 * 10 −2 m 2 ;  = 0.19 kN s 2 /m;  = 3.705;  = 2.1 * 10 8 kN/m 2 ;  = 8.1 * 10 7 kN/m 2 .
Frequency factors  = 4  √  2  4 / are calculated taking bending moment, shear, and axial effects into consideration for   = 0.25,   = 0.5, and   = 0.75 due to circular frequencies of the beam.for the CS beam with analytical results obtained by C ¸atal [7].Tables 2 and 3 show the analytical results for both CC and SS beams which are not solved by C ¸atal [7].
The graphical presentations of analytical solutions are provided in Figures 3, 4, and 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 C ¸atal [7] fail to coincide with the exact results provided in this study.Previous studies [8][9][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 C ¸atal [7], the DTM converges the analytical results provided in the same study.Presumably, the difference is due to the fact that C ¸atal [7] used the length of the beam, L, in three-digit accuracy while in this study the exact values of length, L, 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.

Figure 1 :
Figure 1: (a) beam on elastic foundation and (b) internal forces and deformations of the beam on elastic foundation.

Table 1 :
[7]iation of frequency factor, , with relative stiffness, , for a CS beam on elastic foundation.(Theanalyticalresults of C ¸atal[7]are provided below in bold characters.)

Table 2 :
Variation of frequency factor, , with relative stiffness, , for a CC beam on elastic foundation.
Variation of frequency factor,  is tabulated in Tables1-3.Table1includes comparison of analytical results of this study

Table 3 :
Variation of frequency factor, , with relative stiffness, , for an SS beam on elastic foundation.