Comparison of normalized buckling loads() for the
In most of the real world engineering applications, stability analysis of compressed members is very crucial. There have been many researches dedicated tothe buckling behavior of axially compressed members. On the other hand, obtaining analytical solutions for the buckling behavior of columns with variable cross-section subjected to complicated load configurations are almost impossible in most of the cases. Some of the works related to obtaining analytical or analytical approximate solutions for the column buckling problem are provided below.
The problems of buckling of columns under variable distributed axial loads were solved in detail by Vaziri and Xie  and others. Some analytical closed-form solutions are given by Dinnik , Karman and Biot , Morley, Timoshenko and Gere  and others. One of the detailed references related to the structural stability topic is written by Simitses and Hodges  with detailed discussions. Iyengar  made some analysis on buckling of uniform with several elastic supports. Wang et al.  have given exact mathematical solutions for buckling of structural members for various cases of columns, beams, arches, rings, plates and shells. Ermopoulos  found the solution for buckling of tapered bars axially compressed by concentrated loads applied at various locations along their axes. Li  gave the exact solution for buckling of non-uniform columns under axially concentrated and distributed loading. Lee and Kuo  established an analytical procedure to investigate the elastic stability of a column with elastic supports at the ends under uniformly distributed follower forces. Furthermore, Gere and Carter  investigated and established the exact analytical solutions for buckling of several special types of tapered columns with simple boundary conditions. Solution of the problem of buckling of elastic columns with step varying thickness is established by Arbabei and Li . Stability problems of a uniform bar with several elastic supports using the moment-distribution method were analyzed by Kerekes . The research of Siginer  was about the stability of a column whose flexural stiffness has a continuous linear variation along the column. Moreover, the analytical solutions of a multi-step bar with varying cross section were obtained by Li et al. [16-18]. The energy method was used by Sampaio et al.  to find the solution for the problem of buckling behavior of inclined beam-column. Some of the important researchers who studied the mechanical behavior of beam-columns are Keller , Tadjbakhsh and Keller  and Taylor . Later on, analytical approximate techniques were used for the stability analysis of elastic columns. Coşkun and Atay  and Atay and Coşkun  studied column buckling problems for the columns with variable flexural stiffness and for the columns with continuous elastic restraints by using the variational iteration method which produces analytical approximations. Coşkun [25, 26] used the homotopy perturbation method for buckling of Euler columns on elastic foundations and tilt-buckling of variable stiffness columns. Pınarbaşı  also analyzed the stability of nonuniform rectangular beams using homotopy perturbation method. These techniques were also used successfully in the vibration analysis of Euler-Bernoulli beams and in the vibration of beams on elastic foundations. [28-29]
Recently, by the emergence of new and innovative semi analytical approximation methods, research on this subject has gained momentum. Analytical approximate solution techniques are used widely to solve nonlinear ordinary or partial differential equations, integrodifferential equations, delay equations, etc. The main advantage of employing such techniques is that the problems are considered in a more realistic manner, and the solution obtained is a continuous function which is not the case for the solutions obtained by discretized solution techniques.
The methods that will be used throughout this study are, Adomian Decomposition Method (ADM), Variational Iteration Method (VIM) and Homotopy Perturbation Method (HPM). Each technique will be explained first, and then all will be applied to a selected case study related to the topic of the article.
2. Problem formulation
Derivation of governing equations related to stability analysis is given in detail in Timoshenko and Gere , Simitses and Hodges , and Wang et al. . The reader can also refer to any textbook related to the subject. In this section, only the governing equation will be given for the related cases.
Consider the elastic columns given in Fig.1. The governing equation for the buckling of such columns is
In the case of constant flexural rigidity (
The governing equation given in Eq.(2) is a fourth order differential equation with constant coefficients which makes it possible to obtain analytical solutions easily. However, Eq.(1) includes variable coefficients due to variable flexural rigidity. For this type of differential equations, analytical solutions are limited for the special cases of
In some problems, obtaining analytical solutions is very difficult even for a constant coefficient governing equation. Consider the buckling of a column on an elastic foundation shown in Fig.2.
The governing equation for the column in Fig.2 is
which, for the constant
Although Eq.(8) is a linear equation with constant coefficients, obtaining a solution from Eq.(9) is not that easy. It is very interesting that, even with a software, one can not easily produce the buckling loads in a sequential order from Eq.(9). In view of this experience, an analytical solution for Eq.(7) is almost impossible to obtain except very limited
Hence, analytical approximate techniques are efficient alternatives for solving these problems. By the use of these techniques, a solution which is continuous in the problem domain is possible for any variation in flexural rigidity. These techniques produce the buckling loads in a sequential order, and it is also very easy to obtain the buckling mode shapes from the solution provided by the method used. These are great advantages in the solution of such problems.
3. The methods used in the elastic stability analysis of Eulercolumns
3.1. Adomian Decomposition Method (ADM)
In the ADM a differential equation of the following form is considered
Applying the inverse operator
After integrating source term and combining it with the terms arising from given conditions of the problem, a function
The nonlinear operator is represented by an infinite series of specially generated (Adomian) polynomials forthe specific nonlinearity. Assuming
where, the components of the series are determined recursively as follows:
3.2. Variational Iteration Method (VIM)
According to VIM, the following differential equation may be considered:
3.3. Homotopy Perturbation Method (HPM)
HPM provides an analytical approximate solution for problems at hand as the other previously explained techniques. Brief theoretical steps for the equation of following type can be given as
with boundary conditions. In Eq.(24)
whichsatisfies the following inequalities:
Theapproximate solution of can be obtained as :
4. Case studies
4.1. Buckling of a clamped-pinned column
The governing equation for this case was previously given in Eq.(1). ADM, VIM and HPM will be applied to this equation in order to compute the buckling loads for the clamped-pinned column with constant flexural stiffness,
4.2. Formulation of the algorithms for uniform column
The linear operator and its inverse operator for Eq.(2) is
To keep the formulation a general one for all configurations to be considered, the boundary conditions are chosen as , , and . Suitable values should be replaced in the formulation with these constants. In this case, and should be inserted for the
Finally, the solution is defined by
Based on the formulation given previously, Lagrange multiplier,
An iterative algorithm can be constructed inserting Lagrange multiplier and governing equation into the formulation given in Eq.(31) as
Based on the formulation, Eq.(2) can be divided into two parts as
The solution can be expressed as a power series in
Hence, an approximate solution would be obtained as
Initial guess is very important for the convergence of solution in HPM. A cubic polynomial with four unknown coefficients can be chosen as an initial guess which was shown previously to be an effective one in problems related to Euler beams and columns [23-29].
4.3. Computation of buckling loads
By the use of described techniques, an iterative procedure is constructed and a polynomial including the unknown coefficients resulting from the initial guess is produced as the solution to the governing equation. Besides two unknowns from the initial guess, an additional unknown
where . For a nontrivial solution, determinant of coefficient matrix must be zero. Determinant of matrix yields a characteristic equation in terms of
4.4. Determination of bucklingmode shapes
Buckling mode shapes for the column can also be obtained from the polynomial approximations by the methods considered in this study. Introducing, the buckling loads into the solution, normalized polynomial eigen functions for the mode shapes are obtained from
The same approach can also be employed to predict mode shapes for the cases including variable flexural stiffness.
4.5. Analysis of a uniform column
After applying the procedures explained in the text, the following results are obtained for the buckling loads. Comparison with the exact solutions is also provided in order that one can observe an excellent agreement between the exact results and the computed results.
Twenty iterations are conducted for each method, and the computed values are compared with the corresponding exact values for the first four modes of buckling in the following table.
From the table it can be seen that the computed values are highly accurate which show that the techniques used in the analysis are very effective. Only a few iterations are enough to obtain the critical buckling load which is Mode 1. Additional modes require additional iterations. The table also shows that additional two or three iterations will produce an excellent agreement for Mode 4. Even with twenty iterations, the error is less than 0.014% for all the methods used in the analyses.
The buckling mode shapes of uniform column for the first four modes are depicted in Fig.4. To prevent a possible confusion to the reader, the exact mode shapes and the computed ones are not shown separately in the figure since the obtained mode shapes coincide with the exact ones.
4.6. Buckling of a rod with variable cross-section
A circular rod having a radius changing exponentially is considered in this case. Such a rod is shown below in Fig.5. The function representing the radius would be as
Employing Eq.(49), cross-sectional area and moment of inertia for a section at an arbitrary point
Governing equation for the rod was previously given in Eq.(1) as
4.6.1. Formulation of the algorithms
Application of ADM leads to the following
ADM gives the following formulation with the previously defined fourth order linear operator.
Lagrange multiplier is the same as used in the uniform column case due to the fourth order derivative in Eq.(38). Hence, an algorithm by using VIM can be constructed as
Application of HPM produces the following set of recursive equations as the solution algorithm.
4.6.2. Results of the analyses
The proposed formulations are applied for two different variations,
The buckling mode shapes of the rod for the first four modes are depicted in between Figs.5-9. To demonstrate the effect of variable cross-section in the results, a comparison is made with normalized mode shapes for a uniform rod which are given in Fig.4. Constant flexural rigidity is defined as in these figures.
In this article, some analytical approximation techniques were employed in the elastic stability analysis of Euler columns. In a variety of such methods, ADM, VIM and HPM are widely used, and hence chosen for use in the computations. Firstly, a brief theoretical knowledge was given in the text, and then all of the methods were applied to the selected cases. Since the exact values for the buckling of a uniform rod were available, the analyses were initially conducted for that case. Results showed an excellent agreement with the exact ones that all three methods were highly effective in the computation of buckling loads and corresponding mode shapes. Finally, ADM, VIM and HPM were applied to the buckling of a rod having variable cross section. To this aim, a rod with exponentially varying radius was chosen and buckling loads with their corresponding mode shapes were obtained easily.
This study has shown that ADM, VIM and HPM can be used effectively in the analysis of elastic stability problems. It is possible to construct easy-to-use algorithms with these methods which are highly accurate and computationally efficient.
Struct., 45(3) ( Vaziri H. H. Xie J. Buckling of. columns under. variably distributed. axial loads. Comput 1992
Dinnik, A.N., Design of columns of varying cross-section, Trans.ASME 51, 1929
Karman T. R. Biot M. A. Mathematical Methods. in Engineering. Mc Graw Hill. New York. 1940
Morley A., Critical loads for long tapering struts, Engineering, 104 295 1917
Timoshenko S. P. Gere J. M. Theory of. Elastic Stability. Mc Graw Hill. New York. 1961
Simitses G. J. Hodges D. H. Fundamentals of. structural Stability. Elsevier-Butterworth-Heinemann Publishing. 2005
Iyengar NGR, Structural Stability of Columns and Plates.New York, John Wiley and Sons, 1988
Wang C. M. Wang C. Y. Reddy J. N. Exact Solutions. for Buckling. of Structural. CRC Members Press. L. L. C. Florida 2005
structural Eng. ASCE, 112(6) ( Ermopoulos J. C. Buckling of. tapered bars. under stepped. axial loading. J. 1986
J. Mech. A / Solids, 20( Li Q. S. exact solutions. for buckling. of non-uniform. columns under. axial concentrated. distributed loading. Eur 2001
Sound Vibration, 148(1) ( Lee S. Y. Kuo Y. H. Elastic stability. of non-uniform. columns J. 1991
Struct. Eng.- ASCE, 88( Gere J. M. Carter W. O. Critical buckling. loads for. tapered columns. J. 1962
Struct. Eng.-ASCE, 117(8) ( Arbabi F. Li F. Buckling of. variable cross-section. columns-integral equation. approach J. 1991
Kerekes F., Hulsbos CL., Elastic stability of the top chord of a three-span continuous pony truss bridge, Iowa Eng.Expt. Sta. Bull., ( 1954
Eng. Mech.-ASCE, 118( Siginer A. Buckling of. columns of. variable flexural. stiffness J. 1992
Li, Q.S., Cao, H., Li, G., Stability analysis of bars with multi-segments of varying cross section, Comput.Struct., 53 ( 1994
J. Solids Struuct., 32 ( Li Q. S. Cao H. Li G. Stability analysis. of bars. with varying. cross section. Int 1995
Li, Q.S., Cao, H., Li, G.,Static and dynamic analysis of straight bars with variable cross-section, Comput.Struct., 59 ( 1996
Sampaio Jr., J.H.B, Hundhausen, J.R.,A mathematical model and analytical solution for buckling of inclined beam columns, Appl.Math. Model., 22 ( 1998
Keller, J.B.,The shape of the strongest column, Archive for Rational Mechanics and Analysis, 5 1960
Appl. Mech. ASME, 29 ( Tadjbakhsh I. Keller J. B. Strongest columns. isoperimetric inequalities. for eigenvalues. J. 1962
- 22. Taylor, J.E.,The strongest column- an energy approach, J. Appl. Mech. ASME, 34 (1967), 486-487.
Coskun, S.B., Atay, M.T., Determination of critical buckling load for elastic columns of constant and variable cross-sections using variational iteration method,, Computers and Mathematics with Applications, 58 2009
Atay, M.T., Coskun, S.B., Elastic stability of Euler columns with a continuous elastic restraint using variational iteration method, Computers and Mathematics with Applications, 58 2009
Coskun, S.B., Determination of critical buckling loads for Euler columns of variable flexural stiffness with a continuous elastic restraint using Homotopy Perturbation Method, Int.Journal Nonlinear Sci. and Numer. Simulation, 10(2) ( 2009
Coskun, S.B., Analysis of Tilt-Buckling of Euler Columns with Varying Flexural Stiffness Using Homotopy Perturbation Method, Mathematical Modelling and Analysis, 15(3) 2010
Pinarbasi S. Stability analysis. of nonuniform. rectangular beams. using homotopy. perturbation method. Mathematical Problems. in Engineering. 2012Article ID.197483.
Safa Bozkurt Coşkun, Mehmet Tarik Atay and Baki Öztürk, Transverse Vibration Analysis of Euler-Bernoulli Beams Using Analytical Approximate Techniques, Advances in Vibration Analysis Research, Dr.Farzad Ebrahimi (Ed.), 978-9-53307-209-8InTech, 2011
Ozturk, B., Coskun, S.B., The Homotopy Perturbation Method for free vibration analysis of beam on elastic foundation, Structural Engineering and Mechanics, 37(4) 2011
- 30. Adomian, G., Solving Frontier Problems of Physics: The Decomposition Method, Kluwer, Boston, MA, 1994.
Comput. Modell., 13(7) ( Adomian G. review A. of the. decomposition method. some recent. results for. nonlinear equation. Math 1992
He,J.H., Variational iteration method: a kind of nonlinear analytical technique, Int. J. Nonlin. Mech., 34 ( 1999
J.H., A coupling method of a homotopy technique and a perturbation technique for non-linear problems, Int. J. Nonlin. Mech., 35 ( He J. H. coupling A. method of. a. homotopy technique. perturbation a. technique for. non-linear problems. Int 2000
He, J.H., An elemantary introduction to the homotopy perturbation method, Computers and Mathematics with Applications, 57 2009
He J. H. New interpretation. of homotopy. perturbation method. International Journal. of Modern. Physics B. . 2006
He, J.H., The homotopy perturbation method for solving boundary problems, Phys.Lett. A, 350 ( 2006