Bilevel Programming, Equilibrium, and Combinatorial Problems with Applications to Engineering
View this Special IssueResearch Article  Open Access
The Solution of Fourth Order Boundary Value Problem Arising out of the BeamColumn Theory Using Adomian Decomposition Method
Abstract
Adomian decomposition method (ADM) is applied to linear nonhomogeneous boundary value problem arising from the beamcolumn theory. The obtained results are expressed in tables and graphs. We obtain rapidly converging results to exact solution by using the ADM. This situation indicates that the method is appropriate and reliable for such problems.
1. Introduction
It is possible to model many of the physical events that take place in nature using linear and nonlinear differential equations. This modeling enables us to understand and interpret the particular event in a much better manner. Thus, finding the analytical and approximate solutions of such models with initial and boundary conditions gain importance. Differential equations have had an important place in engineering since many years. Scientists and engineers generally examine systems that undergo changes.
Many methods have been developed to determine the analytical and approximate solutions of linear and nonlinear differential equations with initial and boundary value conditions and among these methods, the ADM [1–6], homotopy perturbation method [7–11], variational iteration method [12–19], and homotopy analysis method [20–25] can be listed. Package software such as Matematica, Maple, or Matlab is used to overcome tedious algebraic operations when using these methods. The determination of the analytical and approximate solutions of linear and nonlinear differential equations is an important topic for civil engineering, because these equations are the mathematical models of complex events that occur in engineering.
In this study, the approximate solution of the fourth order linear nonhomogeneous differential equation with initial and boundary conditions that arises in the beamcolumn theory will be determined via the ADM and comparisons will be made with the existing results in the literature.
The ADM was first put forth in the 1980s by an American scientist named Adomian [26]. This method is based on the decomposition of the unknown function. Using this method it is possible to determine the approximate solutions for the linear and nonlinear ordinary and partial differential equations. The ADM has been used effectively by researchers during 1990 and 2007 especially for the solution of differential and integral equations [27–32]. Using this method a nonlinear problem can be applied directly without discretization and linearization. Hence, this is a method that is preferred by researchers.
In this study, the analytical and approximate solutions of a fourth order linear boundary value problem were calculated using the ADM. In addition, the fourth order linear ordinary differential equation set used in the beamcolumn theory was solved under specific boundary conditions. The ADM used provides realistic solutions without changing the linear or nonlinear differential equation model. Numerical results close to the real solution can be found by calculating the terms of the limited number of decomposition series. It is possible to find the numerical solution of a differential equation using this method without the necessity of indexing.
2. Adomian Decomposition Method
The method has many advantages in comparison with many other traditional methods such as Finite differences, Finite elements, and Galerkin method [33–35]. The method gives converging results adapted to the problems. It is a disadvantage that convergence interval is small and the calculation of the ADM polynomials occurs in nonlinear problems in the ADM.
We now consider second order nonlinear ordinary differential equations with initial conditions. This equation can be written in an operator by where is the lowerorder derivative, which is assumed to be invertible, is a linear differential operator of order greater than , is a nonlinear differential operator, and is a source term. We next apply the inverse operator to both sides of (1) and use the given initial conditions to get where
The ADM consists in decomposing the unknown function of any equation into a sum of infinite number of components given by the decomposition series where the components , are to be determined in a recursive manner. The nonlinear term will be decomposed by the infinite series of polynomials given by where are Adomian polynomials [1, 36]. Substituting (3) and (4) into (2) gives where
The various components of the solution can be easily determined by using the recursive relation
In order to obtain the numerical solutions of closed solution function using the decomposition method; being; the term, can be calculated by taking into account the reduction formula (8). In addition, the series solution of the decomposition written as (10) generally yields results that rapidly converge for physical problems. The convergence of the decomposition series has been examined by many researchers in the literature. The convergence of the decomposition series has been examined theoretically by Cherruault [37]. In addition to these studies, Abbaoui and Cherruault have suggested a new approach in determining the convergence of the decomposition series [38]. These authors have determined the convergence of the decomposition series method by giving new conditions.
3. The Application of the ADM to the Linear Problem with Boundary Condition That Arises in BeamColumn Theory
3.1. Problem Definition
Fourth order differential equations consist of various physical problems that are related to the elastic stability theory. Differential relations should be established between the effects of various crosssection effects in order to understand beamcolumn problems better.
When a cross section of distance shown in Figure 1(b) is taken from a beamcolumn subject to both the axial load and the spread load perpendicular to the axis as shown in Figure 1(a), internal forces arise in the element. When the equilibrium equation in the direction is written with the following ordinary differential equation is found [39–42]:
(a)
(b)
The algebraic sums of the forces acting on both surfaces of the crosssection element are the same due to equilibrium. Consider
Here, is the shear force acting on the surface of the element, whereas is the bending moment that tries to bend the crosssection element.
If rotations are assumed to be small and the second order terms in terms of are neglected, then (13) becomes Since rotations are assumed to be small, if , then (14) becomes
Here, represents bending rigidity. If the derivative of both sides of (15) is taken in terms of , then the fourth order linear differential equation for the elastic curve is found as such:
3.2. Application of the ADM to the Problem
Example 1. We consider the fourth order linear nonhomogeneous differential equation [43]:
with the boundary conditions of
If (18) is written out in operator form, we obtain
Here,
are the derivative and integral operators. If the inverse operator is applied to (19) and initial conditions are taken, we find
where and . If we use (3) in (21), then we find
whereas the reduction formula given below can be written using (22):
The and constants in the reduction formula will be determined using boundary conditions (18) after finding the decomposition series. From the reduction relation, we can obtain the solution terms of the decomposition series as
If these terms are placed in (4), we obtain the approximate solution obtained via the ADM using the five terms of (17) and problem (18) can be written as
If we use boundary conditions (18) in solutions series (26), then we have and . The analytical solution of problem (17) and (18) is [43]
The numerical results and graphs obtained for five terms using the solution series in (26) have been given below.
4. Results and Discussion
In this study, the approximate solution of the fourth order boundary value problem arising in beamcolumn theory has been determined using the ADM. This method can be applied on differential equations without the need for discretization, indexing, or linearization.
Nonhomogeneous problem has been handled regarding the topic and the obtained results have been given in Table 1 and Figure 2. As can be seen from the table and figure, the method used has given results that converge rapidly to the analytical solution with the removal of a few terms from the solution series. This shows that the method is suitable and reliable for such problems. In addition, the handled example has been compared with the results obtained by Chen and Atsuta [41, 42] using nonpolynomial spline method. As seen in Table 2 values close to those obtained in this study were found only in case when the “” step was small. New terms can be added to the solution series of the ADM to obtain results that are much better than those of the spline method.


In Figure 2, five terms of the decomposition method have been taken from the analytical solution and decomposition series of example have been drawn in two dimensional graphs. As seen in these graphs, the analytical and approximate solutions cannot be distinguished.
In conclusion, it has been determined that the ADM can be applied to the linear homogeneous and nonhomogeneous boundary value problems that arise in civil engineering in the beamcolumn theory. Solutions that converge rapidly to the analytical solution can be found without changing the nature of the physical phenomenon. In addition, the calculations for this method can be carried out using software such as Mathematica, Maple, and Matlab.
Conflict of Interests
The author declares that there is no conflict of interests.
References
 G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic Publishers, Boston, Mass, USA, 1994. View at: MathSciNet
 G. Adomian, “A review of the decomposition method and some recent results for nonlinear equations,” Mathematical and Computer Modelling, vol. 13, no. 7, pp. 17–43, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian and R. Rach, “Equality of partial solutions in the decomposition method for linear or nonlinear partial differential equations,” Computers & Mathematics with Applications, vol. 19, no. 12, pp. 9–12, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Inc, “On numerical solutions of partial differential equations by the decomposition method,” Kragujevac Journal of Mathematics, vol. 26, pp. 153–164, 2004. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 M. Inc, “Decomposition method for solving parabolic equations in finite domains,” Journal of Zhejiang University SCIENCE A, vol. 6, no. 10, pp. 1058–1064, 2005. View at: Publisher Site  Google Scholar
 M. Inc, Y. Cherruault, and K. Abbaoui, “A computational approach to the wave equations: an application of the decomposition method,” Kybernetes, vol. 33, no. 1, pp. 80–97, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. M. Odibat, “A new modification of the homotopy perturbation method for linear and nonlinear operators,” Applied Mathematics and Computation, vol. 189, no. 1, pp. 746–753, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. R. S. Alizadeh, G. G. Domairry, and S. Karimpour, “An approximation of the analytical solution of the linear and nonlinear integrodifferential equations by homotopy perturbation method,” Acta Applicandae Mathematicae, vol. 104, no. 3, pp. 355–366, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y.G. Wang, H.F. Song, and D. Li, “Solving twopoint boundary value problems using combined homotopy perturbation method and Green's function method,” Applied Mathematics and Computation, vol. 212, no. 2, pp. 366–376, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Biazar and H. Ghazvini, “Convergence of the homotopy perturbation method for partial differential equations,” Nonlinear Analysis: Real World Applications, vol. 10, no. 5, pp. 2633–2640, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C.S. Liu, “The essence of the homotopy analysis method,” Applied Mathematics and Computation, vol. 216, no. 4, pp. 1299–1303, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Momani and S. Abuasad, “Application of He's variational iteration method to Helmholtz equation,” Chaos, Solitons & Fractals, vol. 27, no. 5, pp. 1119–1123, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. M. Abulwafa, M. A. Abdou, and A. A. Mahmoud, “Nonlinear fluid flows in pipelike domain problem using variationaliteration method,” Chaos, Solitons & Fractals, vol. 32, no. 4, pp. 1384–1397, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. H. Sweilam and M. M. Khader, “Variational iteration method for one dimensional nonlinear thermoelasticity,” Chaos, Solitons & Fractals, vol. 32, no. 1, pp. 145–149, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Xu, “Variational iteration method for solving integral equations,” Computers & Mathematics with Applications, vol. 54, no. 78, pp. 1071–1078, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, A.M. Wazwaz, and L. Xu, “The variational iteration method: reliable, efficient, and promising,” Computers & Mathematics with Applications, vol. 54, no. 78, pp. 879–880, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Xu, J.H. He, and A.M. Wazwaz, “Preface variational iteration method—reality, potential, and challenges,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 1–2, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 S. B. Coşkun and M. T. Atay, “Analysis of convective straight and radial fins with temperaturedependent thermal conductivity using variational iteration method with comparison with respect to finite element analysis,” Mathematical Problems in Engineering, vol. 2007, Article ID 42072, 15 pages, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. T. Atay and S. B. Coşkun, “Effects of nonlinearity on the variational iteration solutions of nonlinear twopoint boundary value problems with comparison with respect to finite element analysis,” Mathematical Problems in Engineering, vol. 2008, Article ID 857296, 10 pages, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Inc and Y. Uǧurlu, “Numerical simulation of the regularized long wave equation by He's homotopy perturbation method,” Physics Letters A: General, Atomic and Solid State Physics, vol. 369, no. 3, pp. 173–179, 2007. View at: Publisher Site  Google Scholar
 A. S. Bataineh, M. S. M. Noorani, and I. Hashim, “Homotopy analysis method for singular IVPs of EmdenFowler type,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1121–1131, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Abbasbandy, E. Babolian, and M. Ashtiani, “Numerical solution of the generalized Zakharov equation by homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 12, pp. 4114–4121, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Rashidi and S. Dinarvand, “Purely analytic approximate solutions for steady threedimensional problem of condensation film on inclined rotating disk by homotopy analysis method,” Nonlinear Analysis: Real World Applications, vol. 10, no. 4, pp. 2346–2356, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Song and H. Zhang, “Solving the fractional BBMBurgers equation using the homotopy analysis method,” Chaos, Solitons & Fractals, vol. 40, no. 4, pp. 1616–1622, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Liao, “On the relationship between the homotopy analysis method and Euler transform,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 6, pp. 1421–1431, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian, Nonlinear Stochastic Operator Equations, Academic Press, San Diego, Calif, USA, 1986. View at: MathSciNet
 M. Inc and M. Işık, “Adomian decomposition method for threedimensional parabolic equation with nonclassic boundary conditions,” Journal of Analysis, vol. 11, pp. 43–51, 2003. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 F. Abdelwahid, “A mathematical model of Adomian polynomials,” Applied Mathematics and Computation, vol. 141, no. 23, pp. 447–453, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Babolian and S. Javadi, “New method for calculating Adomian polynomials,” Applied Mathematics and Computation, vol. 153, no. 1, pp. 253–259, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Cherruault, M. Inc, and K. Abbaoui, “On the solution of the nonlinear Kortewegde Vries equation by the decomposition method,” Kybernetes, vol. 31, no. 5, pp. 766–772, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 A. M. A. ElSayed and M. Gaber, “The Adomian decomposition method for solving partial differential equations of fractal order in finite domains,” Physics Letters A, vol. 359, no. 3, pp. 175–182, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Momoniat, T. A. Selway, and K. Jina, “Analysis of Adomian decomposition applied to a thirdorder ordinary differential equation from thin film flow,” Nonlinear Analysis: Theory, Methods & Applications, vol. 66, no. 10, pp. 2315–2324, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. L. Thompson and P. M. Pinsky, “A Galerkin leastsquares finite element method for the twodimensional Helmholtz equation,” International Journal for Numerical Methods in Engineering, vol. 38, no. 3, pp. 371–397, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Dolbow and T. Belytschko, “Numerical integration of the Galerkin weak form in meshfree methods,” Computational Mechanics, vol. 23, no. 3, pp. 219–230, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. N. Atluri and T. Zhu, “A new meshless local PetrovGalerkin (MLPG) approach in computational mechanics,” Computational Mechanics, vol. 22, no. 2, pp. 117–127, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. M. Wazwaz, Partial Differential Equations: Methods and Applications, Balkema Publishers, Rotterdam, The Netherlands, 2002.
 Y. Cherruault, “Convergence of Adomian's method,” Kybernetes, vol. 18, no. 2, pp. 31–38, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. Abbaoui and Y. Cherruault, “New ideas for proving convergence of decomposition methods,” Computers & Mathematics with Applications, vol. 29, no. 7, pp. 103–108, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. P. Timoshenko and J.M. Gere, Theory of Elastic Stability, McGrawHill, New York, NY, USA, 2nd edition, 1961. View at: MathSciNet
 S. P. Timoshenko and J.M. Gere, Theory of Elastic Stability, McGrawHill, New York, NY, USA, 1985. View at: MathSciNet
 W. F. Chen and T. Atsuta, Theory of BeamColumns, vol. 1 of Inplane Behavior and Design, McGrawHill, New York, NY, USA, 1976.
 W. F. Chen and T. Atsuta, Theory of BeamColumns, vol. 2 of Space Behavior and Design, McGrawHill, New York, NY, USA, 1977.
 O. A. Taiwo and O. M. Ogunlaran, “A nonpolynomial spline method for solving linear fourthorder boundaryvalue problems,” International Journal of Physical Sciences, vol. 6, no. 13, pp. 3246–3254, 2011. View at: Google Scholar
Copyright
Copyright © 2014 Omer Kelesoglu. 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.