MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 649471 10.1155/2014/649471 649471 Research Article The Solution of Fourth Order Boundary Value Problem Arising out of the Beam-Column Theory Using Adomian Decomposition Method Kelesoglu Omer Kavun Sergii V. Department of Civil Engineering Technology Faculty Firat University 23119 Elazig Turkey firat.edu.tr 2014 3032014 2014 23 11 2013 02 03 2014 30 3 2014 2014 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.

Adomian decomposition method (ADM) is applied to linear nonhomogeneous boundary value problem arising from the beam-column 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 , homotopy perturbation method , variational iteration method , and homotopy analysis method  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 beam-column 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 . 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 . 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 beam-column 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 . 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 (1)Lu+Ru+Nu=f, where L is the lower-order derivative, which is assumed to be invertible, R is a linear differential operator of order greater than L, N is a nonlinear differential operator, and f is a source term. We next apply the inverse operator L-1 to both sides of (1) and use the given initial conditions to get (2)u(x)=u(0)+xu(0)+L-1(g)-L-1(Ru)-L-1(Nu), where (3)L=d2dx2,L-1=0x0xdxdx.

The ADM consists in decomposing the unknown function u(x) of any equation into a sum of infinite number of components given by the decomposition series (4)u(x)=n=0un(x), where the components un(x), n0 are to be determined in a recursive manner. The nonlinear term Nu will be decomposed by the infinite series of polynomials given by (5)Nu=n=0An, where An are Adomian polynomials [1, 36]. Substituting (3) and (4) into (2) gives (6)n=0un=u0(x)-L-1R(n=0un(x)-L-1n=0An),   where (7)u0(x)=u(0)+xu(0)+L-1(g).

The various components un of the solution u can be easily determined by using the recursive relation (8)u0(x)=u(0)+xu(0)+L-1g(x),uk+1(x)=-L-1(Ruk)-L-1(Ak),k0.

In order to obtain the numerical solutions of u(x) closed solution function using the decomposition method; (9)φn=k=0nuk(x)k0, being; the term, (10)u(x)=limφnn, 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 . In addition to these studies, Abbaoui and Cherruault have suggested a new approach in determining the convergence of the decomposition series . 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 Beam-Column 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 cross-section effects in order to understand beam-column problems better.

When a cross section of distance dx shown in Figure 1(b) is taken from a beam-column subject to both the P axial load and the q spread load perpendicular to the axis as shown in Figure 1(a), internal forces arise in the element. When the equilibrium equation in the y direction is written with (11)q=-dVdx, the following ordinary differential equation is found : (12)-V+qdx+(V+dV)=0.

Cross-sectional analysis of thecolumn-beam element.

The algebraic sums of the forces acting on both surfaces of the cross-section element are the same due to equilibrium. Consider (13)M+qdxdx2+(V+dV)dx-(M+dM)+Pdydx=0.

Here, V is the shear force acting on the surface of the element, whereas M is the bending moment that tries to bend the cross-section element.

If rotations are assumed to be small and the second order terms in terms of dx are neglected, then (13) becomes (14)V=dMdx-Pdydx. Since rotations are assumed to be small, if d2/dx2=-M/EI, then (14) becomes (15)-V=EId3ydx3+Pdydx.

Here, EI represents bending rigidity. If the derivative of both sides of (15) is taken in terms of x, then the fourth order linear differential equation for the elastic curve is found as such: (16)EId4ydx4+Pd2ydx2=q(x).

3.2. Application of the ADM to the Problem Example 1.

We consider the fourth order linear nonhomogeneous differential equation : (17)d4ydx4-2d2ydx2+y=-8ex,x[0,1], with the boundary conditions of (18)y(0)=y′′(0)=1,y(1)=y′′(1)=-e.

If (18) is written out in operator form, we obtain (19)Ly=2y′′-y-8ex. Here, (20)L=d4dx4,L-1(·)=0x0x0x0x(·)dxdxdxdx are the derivative and integral operators. If the L-1 inverse operator is applied to (19) and initial conditions are taken, we find (21)y(x)=Ax+13!Bx3-8L-1(ex)+2L-1(y′′)-L-1(y), where A=y(0) and B=y′′′(0). If we use (3) in (21), then we find (22)n=0yn(x)=Ax+13!Bx3-8L-1ex+2L-1(n=0yn′′(x))-L-1(n=0yn(x)), whereas the reduction formula given below can be written using (22): (23)y0(x)=Ax+13!Bx3-8L-1(ex)yn+1(x)=2L-1(yn′′)-L-1(yn)aaaaaaaaaaaaaaan0.

The A and B 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 (24)y0(x)=8-8ex+(8+A)x+4x2+16(8+B)x3,y1(x)=8-8ex+8x+4x2+4x33+x43-1120(-8+A-2B)x5-x690-(8+B)x75040,y2(x)=8-8ex+8x+4x2+4x33+x43-x515+x690-(A-2(2+B))x72520-x81680+(A-4(6+B))x9362880+x10453600+(8+B)x1139916800,y3(x)=8-8ex+8x+4x2+4x33+x43+x515+x690+x7630+x85040-(-2+A-2B)x990720-x1064800+(-14+A-3B)x119979200+x1211975040-(-40+A-6B)x136227020800-x1410897286400-(8+B)x151307674368000,y4(x)=8-8ex+8x+4x2+4x33+x43+x515+x690+x7630+x85040+x945360+x10453600-(-1+A-2B)x114989600-x123991680+(-30+3A-8B)x131556755200+17x1410897286400-(-6B+3A-12B)x15653837184000-x16373621248000+(A-8(7+B))x17355687428096000+x18800296713216000+(8+B)x19121645100408832000.

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 (25)y(x)=i=04yi(x)=y0(x)+y1(x)+y2(x)+y3(x)+y4(x),(26)y(x)=40-40ex+(40+A)x+20x2+16(40+B)x3+4x43-1120(A-2(16+B))x5+x645-(-16+2A-3B)x75040-x85040-(8+3A-4B)x9362880-x1090720-(40+4A-5B)x1139916800-x125987520+(-80+11A-26B)x136227020800+x14681080400-(-128+6A-23B)x151307674368000-x16373621248000+(A-8(7+B))x17355687428096000+x18800296713216000+(8+B)x19121645100408832000.

If we use boundary conditions (18) in solutions series (26), then we have A=1 and B=-3. The analytical solution of problem (17) and (18) is  (27)y(x)=x(1-x)ex.

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 beam-column 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 “h” 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.

Numerical results obtained using the five terms of the decomposition series.

x Analytical solution y(x) Approximate solution 5 Error |y(x)-5|
0 0. 0. 0.
0.1 0.099465 0.099465 2.92821 × 10−15
0.2 0.195424 0.195424 −2.16493 × 10−15
0.3 0.283470 0.283470 6.10623 × 10−16
0.4 0.358038 0.358038 −1.30451 × 10−13
0.5 0.412180 0.412180 −2.37055 × 10−12
0.6 0.437309 0.437309 −2.58711 × 10−11
0.7 0.422888 0.422888 −1.94957 × 10−10
0.8 0.356087 0.356087 −1.12373 × 10−9
0.9 0.221364 0.221364 −5.27758 × 10−9
1.0 0. 2.10884 × 10−8 2.10884 × 10−8

Comparison of the absolute error obtained using the spline method (SM) and the ADM.

x SM h=1/5 SM h=1/10 ADM 5 CPU time
0 0 0 0 0.542
0.2 1.960 E - 5 1.228 E - 6 - 2.164 E - 15 0.846
0.4 3.211 E - 5 2.012 E - 6 - 1.304 E - 13 0.870
0.6 3.568 E - 5 2.235 E - 6 - 2.587 E - 11 0.875
0.8 2.683 E - 5 1.681 E - 6 - 1.123 E - 9 0.891
1 0 0 - 2.108 E - 8 0.902

Graph showing the ADM and analytical solution obtained using the five terms of the decomposition series.

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 beam-column 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.

Adomian G. Solving Frontier Problems of Physics: The Decomposition Method 1994 Boston, Mass, USA Kluwer Academic Publishers MR1282283 Adomian G. A review of the decomposition method and some recent results for nonlinear equations Mathematical and Computer Modelling 1990 13 7 17 43 10.1016/0895-7177(90)90125-7 MR1071436 ZBL0713.65051 Adomian G. Rach R. Equality of partial solutions in the decomposition method for linear or nonlinear partial differential equations Computers & Mathematics with Applications 1990 19 12 9 12 10.1016/0898-1221(90)90246-G MR1044436 ZBL0702.35058 Inc M. On numerical solutions of partial differential equations by the decomposition method Kragujevac Journal of Mathematics 2004 26 153 164 MR2126008 ZBL1080.65098 Inc M. Decomposition method for solving parabolic equations in finite domains Journal of Zhejiang University SCIENCE A 2005 6 10 1058 1064 2-s2.0-27944463199 10.1631/jzus.2005.A1058 Inc M. Cherruault Y. Abbaoui K. A computational approach to the wave equations: an application of the decomposition method Kybernetes 2004 33 1 80 97 10.1108/03684920410514535 MR2110374 ZBL1058.65112 Odibat Z. M. A new modification of the homotopy perturbation method for linear and nonlinear operators Applied Mathematics and Computation 2007 189 1 746 753 10.1016/j.amc.2006.11.188 MR2330252 ZBL1122.65092 Alizadeh S. R. S. Domairry G. G. Karimpour S. An approximation of the analytical solution of the linear and nonlinear integro-differential equations by homotopy perturbation method Acta Applicandae Mathematicae 2008 104 3 355 366 10.1007/s10440-008-9261-z MR2453431 ZBL1162.65419 Wang Y.-G. Song H.-F. Li D. Solving two-point boundary value problems using combined homotopy perturbation method and Green's function method Applied Mathematics and Computation 2009 212 2 366 376 10.1016/j.amc.2009.02.036 MR2531244 ZBL1166.65362 Biazar J. Ghazvini H. Convergence of the homotopy perturbation method for partial differential equations Nonlinear Analysis: Real World Applications 2009 10 5 2633 2640 10.1016/j.nonrwa.2008.07.002 MR2523226 ZBL1173.35395 Liu C.-S. The essence of the homotopy analysis method Applied Mathematics and Computation 2010 216 4 1299 1303 10.1016/j.amc.2010.02.022 MR2607239 ZBL1194.34010 Momani S. Abuasad S. Application of He's variational iteration method to Helmholtz equation Chaos, Solitons & Fractals 2006 27 5 1119 1123 10.1016/j.chaos.2005.04.113 MR2164845 ZBL1086.65113 Abulwafa E. M. Abdou M. A. Mahmoud A. A. Nonlinear fluid flows in pipe-like domain problem using variational-iteration method Chaos, Solitons & Fractals 2007 32 4 1384 1397 10.1016/j.chaos.2005.11.050 MR2286305 ZBL1128.76019 Sweilam N. H. Khader M. M. Variational iteration method for one dimensional nonlinear thermoelasticity Chaos, Solitons & Fractals 2007 32 1 145 149 10.1016/j.chaos.2005.11.028 MR2271107 ZBL1131.74018 Xu L. Variational iteration method for solving integral equations Computers & Mathematics with Applications 2007 54 7-8 1071 1078 10.1016/j.camwa.2006.12.053 MR2398128 ZBL1141.65400 He J.-H. Wazwaz A.-M. Xu L. The variational iteration method: reliable, efficient, and promising Computers & Mathematics with Applications 2007 54 7-8 879 880 10.1016/j.camwa.2006.12.056 MR2395624 ZBL1138.65301 Xu L. He J.-H. Wazwaz A.-M. Preface variational iteration method—reality, potential, and challenges Journal of Computational and Applied Mathematics 2007 207 1 1 2 10.1016/j.cam.2006.07.021 MR2332940 Coşkun S. B. Atay M. T. Analysis of convective straight and radial fins with temperature-dependent thermal conductivity using variational iteration method with comparison with respect to finite element analysis Mathematical Problems in Engineering 2007 2007 15 42072 10.1155/2007/42072 MR2365727 ZBL05311292 Atay M. T. Coşkun S. B. Effects of nonlinearity on the variational iteration solutions of nonlinear two-point boundary value problems with comparison with respect to finite element analysis Mathematical Problems in Engineering 2008 2008 10 857296 10.1155/2008/857296 MR2402980 ZBL1151.76030 Inc M. Uǧurlu Y. Numerical simulation of the regularized long wave equation by He's homotopy perturbation method Physics Letters A: General, Atomic and Solid State Physics 2007 369 3 173 179 2-s2.0-34548509923 10.1016/j.physleta.2007.04.074 Bataineh A. S. Noorani M. S. M. Hashim I. Homotopy analysis method for singular IVPs of Emden-Fowler type Communications in Nonlinear Science and Numerical Simulation 2009 14 4 1121 1131 10.1016/j.cnsns.2008.02.004 MR2468944 ZBL1221.65197 Abbasbandy S. Babolian E. Ashtiani M. Numerical solution of the generalized Zakharov equation by homotopy analysis method Communications in Nonlinear Science and Numerical Simulation 2009 14 12 4114 4121 10.1016/j.cnsns.2009.03.001 MR2537572 ZBL1221.65269 Rashidi M. M. Dinarvand S. Purely analytic approximate solutions for steady three-dimensional problem of condensation film on inclined rotating disk by homotopy analysis method Nonlinear Analysis: Real World Applications 2009 10 4 2346 2356 10.1016/j.nonrwa.2008.04.018 MR2508442 ZBL1163.34307 Song L. Zhang H. Solving the fractional BBM-Burgers equation using the homotopy analysis method Chaos, Solitons & Fractals 2009 40 4 1616 1622 10.1016/j.chaos.2007.09.042 MR2531252 ZBL1198.65205 Liao S. On the relationship between the homotopy analysis method and Euler transform Communications in Nonlinear Science and Numerical Simulation 2010 15 6 1421 1431 10.1016/j.cnsns.2009.06.008 MR2576768 ZBL1221.65206 Adomian G. Nonlinear Stochastic Operator Equations 1986 San Diego, Calif, USA Academic Press MR872695 Inc M. Işık M. Adomian decomposition method for three-dimensional parabolic equation with non-classic boundary conditions Journal of Analysis 2003 11 43 51 MR2217980 ZBL1132.35385 Abdelwahid F. A mathematical model of Adomian polynomials Applied Mathematics and Computation 2003 141 2-3 447 453 10.1016/S0096-3003(02)00266-7 MR1972923 ZBL1027.65072 Babolian E. Javadi S. New method for calculating Adomian polynomials Applied Mathematics and Computation 2004 153 1 253 259 10.1016/S0096-3003(03)00629-5 MR2063159 ZBL1055.65068 Cherruault Y. Inc M. Abbaoui K. On the solution of the non-linear Korteweg-de Vries equation by the decomposition method Kybernetes 2002 31 5 766 772 10.1108/03684920210428263 MR1926711 El-Sayed A. M. A. Gaber M. The Adomian decomposition method for solving partial differential equations of fractal order in finite domains Physics Letters A 2006 359 3 175 182 10.1016/j.physleta.2006.06.024 MR2253516 ZBL1236.35003 Momoniat E. Selway T. A. Jina K. Analysis of Adomian decomposition applied to a third-order ordinary differential equation from thin film flow Nonlinear Analysis: Theory, Methods & Applications 2007 66 10 2315 2324 10.1016/j.na.2006.03.021 MR2311034 ZBL1127.34302 Thompson L. L. Pinsky P. M. A Galerkin least-squares finite element method for the two-dimensional Helmholtz equation International Journal for Numerical Methods in Engineering 1995 38 3 371 397 10.1002/nme.1620380303 MR1313407 ZBL0844.76060 Dolbow J. Belytschko T. Numerical integration of the Galerkin weak form in meshfree methods Computational Mechanics 1999 23 3 219 230 10.1007/s004660050403 MR1691704 ZBL0963.74076 Atluri S. N. Zhu T. A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics Computational Mechanics 1998 22 2 117 127 10.1007/s004660050346 MR1649420 ZBL0932.76067 Wazwaz A. M. Partial Differential Equations: Methods and Applications 2002 Rotterdam, The Netherlands Balkema Publishers Cherruault Y. Convergence of Adomian's method Kybernetes 1989 18 2 31 38 10.1108/eb005812 MR1009979 ZBL0697.65051 Abbaoui K. Cherruault Y. New ideas for proving convergence of decomposition methods Computers & Mathematics with Applications 1995 29 7 103 108 10.1016/0898-1221(95)00022-Q MR1321262 ZBL0832.47051 Timoshenko S. P. Gere J.M. Theory of Elastic Stability 1961 2nd New York, NY, USA McGraw-Hill MR0134026 Timoshenko S. P. Gere J.M. Theory of Elastic Stability 1985 New York, NY, USA McGraw-Hill MR0134026 Chen W. F. Atsuta T. Theory of Beam-Columns 1976 1 New York, NY, USA McGraw-Hill In-plane Behavior and Design Chen W. F. Atsuta T. Theory of Beam-Columns 1977 2 New York, NY, USA McGraw-Hill Space Behavior and Design Taiwo O. A. Ogunlaran O. M. A non-polynomial spline method for solving linear fourth-order boundary-value problems International Journal of Physical Sciences 2011 6 13 3246 3254 2-s2.0-79960865963