In recent years, the ADM has been applied to solve many problems in science and engineering due to the rapid convergence and reduced computational cost. The advantage of this method is that it solves the problems directly without linearization, perturbation or other transforms [1, 2]. In the solution of linear as well as nonlinear deterministic and stochastic equations, this method introduces a rapidly convergent series solution [3, 4]. The method is modified by researchers to solve many problems, which include linear and nonlinear ordinary and partial differential equations, integral equations and their systems. Having exact solutions in some types of problems as a sum of infinite series has encouraged scientists to conduct research by applying the ADM to a wide range of problems and functional equations. Many well known problems in science and engineering are examined by the ADM; the results, which are similar or exactly the same, are compared with either analytical or numerical solutions of these problems .
The important part of the ADM is representing the nonlinear terms in the equation by Adomian polynomials built on different techniques. The most widely used Adomian polynomials are generated by George Adomian based on a Taylor series expansion [1, 2]. Using a Taylor series expansion an around initial iteration point, the nonlinear terms in the equation can be expressed as a power series of unknown functions. Then, setting the same degrees of unknown values, the Adomian polynomials can be obtained. Another way to calculate Adomian polynomials is built on a Neumann series expansion, where the nonlinear terms can be represented by a Neumann series depending on a parameter, such as λ, and the Adomian polynomials can be represented as the coefficients of λ in the series. Finally, the parametrization method is also suitable while working with nonlinear cases. Basically, the Adomian polynomials can be obtained by comparison with the expansion of a Neumann series and the expansion of nonlinear terms.
Although the method gives rapidly convergent results compared with other semi-analytic methods, there are some studies about the treatment of the method for convergence and computing time issues for the calculation of Adomian polynomials [6, 7]. The most important modification on the ADM, called Modified ADM, is to eliminate the repeating process in the recurrence formulation by choosing the initial term as two parts. Thus, the acceleration in the convergence of the solution can be increased and the number of iterations is minimized to approach the solution more rapidly .
Applications of the ADM to engineering problems have been carried out by researchers. The simple solution of the Cauchy problem for the wave equation is examined by Lesnic  using the ADM. As a conclusion, it is shown that the ADM can be applied to similar problems with various types of boundary conditions. Since the method can be easily applied to linear as well as nonlinear problems, many engineering problems can be solved with it. The study of bending vibration of beams is very important in a wide variety of areas such as bridges, tall buildings, trains, aircraft, the dynamics and control of rockets and producing machine tools. The Euler-Bernoulli beam equation, which is governed by the fourth-order differential equation, was solved by both numerical and analytical methods by many researchers. The vibration of beams was examined numerically by Reddy  using the finite element method. Also, in Haddadpour , the vibration of an Euler-Bernoulli beam was studied using the ADM and an exact closed-form solution was given for a simply supported beam. These results were compared with the solution of the same problem by the modal analysis method. Many researchers have examined the vibration of composite materials made of Functionally Graded Materials (FGMs) which are essentially two-phase particulate composites synthesized in a way such that the volume fractions of the constituents vary constantly from one point to another. The concept of FGMs could provide great flexibility in material design by controlling both the composition profile and the microstructure [12–14]. As an application of the ADM on FGMs, the one dimensional diffusion equation with variable coefficients is examined by Sahin and Karatay . In the problem, the material composition, which is named as a nonhomogeneity parameter, is modeled as both an exponential function and a power function. Finally, results were introduced in terms of a varying nonhomogeneity parameter of the composition profile. For more details about the beams and iterative methods, one can see [16–21].
In this study, the clamped simply supported composite Euler-Bernoulli beam that is governed by the fourth-order differential equation with variable coefficients was solved using the ADM. The composition profile was modeled with an exponential function of the spatial coordinate with a parameter named as the nonhomogeneity parameter. The results in terms of the natural frequencies of the beam were compared with closed-form solutions for clamped simply supported Euler-Bernoulli beam made of homogeneous material, in which the nonhomogeneity parameter was assumed as zero. These results were also introduced in figures for different values of the nonhomogeneity parameter for a composite beam. In the Euler-Bernoulli beam theory, any point on the beam can only move through the vertical direction and this movement is called deflection by w(x, t) of that point. The transverse deflection of the composite beam is governed by the fourth-order partial differential equation with variable coefficients: (1)
where μ(x) is the mass per unit length, EI(x) = E(x)I(x), E(x) is the elastic modulus, I(x) is the moment of inertia and q(x, t) is the load per unit length. The composition profile for a functionally graded Euler-Bernoulli beam can be represented mathematically by an exponential function such that the engineering parameters are given as (2) where the parameter α is named as a nonhomogeneity parameter of FGM and the values E0I0 and μ0 are constants. The ADM will be applied to solve the problem in which the corresponding eigenfunctions of the eigenvalues will be held in terms of generalized Fourier series.
2 Adomian Decomposition Method (ADM)
In this section, we examine a differential equation of the form (3) where L is an easily invertible operator, R is the remainder of the linear operator, N represents the nonlinear operator and the function g represents the nonhomogeneous part of a differential equation. Since L is ans easily invertible n-th order differential operator, its inverse L−1 to both sides of (3) along with the initial conditions, we obtain (4) where
The ADM defines a solution with an infinite series of the same form as in , (5) and the nonlinear term Nu can be decomposed by an infinite series of polynomials as (6) where A0, A1, A2, ...,An are Adomian polynomials of u0, u1, u2, ..., un given by (7) and explicitly formulated by
If the series converges in a suitable way, the approximated general solution of a differential equation can be expressed as (10)
3 Vibration of functionally graded beam
The Equation (1) can be rewritten in an operator form, with definitions given in (2), as (11) where and (12) subject to the initial conditions which also satisfy the boundary conditions (13) related to the clamped simply-supported beam like.
In the solution process, the problem will be separated into two parts: (14) where u(x, t) and v(x, t) are the respective solutions of the homogeneous and nonhomogeneous equations, (15) (16)
3.1 Solution of homogeneous part
Applying the inverse operator, to both sides of (15), the solution of the homogeneous equation subject to initial conditions can be obtained as (17)
Consequently, the solution is given by (18) where . It is clear that a deficient solution may be obtained due to the iteration process, along with the solution ui(x, t) that is found by applying the differential operator . In the solution process, it is easy to obtain a zero solution or a solution which does not satisfy the boundary conditions of the problem. From now on, to accomplish the deficiency, both functions f(x) and g(x) will be expressed in the form of a generalized Fourier series in terms of the Sturm-Liouville boundary value problem. Thus, the coefficients of the Fourier series will be chosen such that these functions satisfy the boundary conditions of the problem. Assume that f(x) and g(x) are given functions in a continuous function space Cp(a, b). When an orthonormal set of functions ϕj (x), (j = 1, 2, . . .) in Cp(a, b) is specified, it may be possible to represent f(x) and g(x) by a linear combination of those functions. This may be generalized to an infinite series that converges to f(x) and g(x) at all but possibly a finite number of points on the fundamental interval 0 < x < l as follows: (19) and, due to orthogonality it can be shown that (20)
The appropriate boundary value problem to find suitable eigenvalues and eigenfunctions for the graded beam problem can be expressed as (21) where and ϕj(x) represent eigenvalues and eigenfunctions, respectively. Then, the equation (21) can be written explicitly as the fourth order differential equation (22) and the solution is obtained as (23)
The unknowns (c1, c2, c3, c4) are determined by the end conditions of the beam such as (24)
The application of boundary conditions gives us the zero solution for ϕ(x) that is not acceptable for general vibrating systems. The characteristic equation of the determinant obtained by the end conditions for fixed values of the nonhomogeneity parameter α infinitely many solutions. For fixed values of l = 1, ϵ = 1 and α = 1, the first four roots of the characteristic equation which are named as natural frequencies of the vibration can be found as (25)
For each value of βj, it is easy to obtain all unknown coefficients assuming that c1 = 1 and then, the normalized shape functions (eigenfunctions) ϕj(x) in Equation (23) can be easily found. On the other hand, the solution of the homogeneous equation in (15) can be obtained by substituting initial conditions, along with Equation (20), into (18) as (26)
Finally, using the identities for series expansions of sine and cosine functions, the solution of the homogeneous equation can be simplified to (27) where the normalized shape functions for the first four eigenvalues can be obtained as (28) (29) (30) (31) and, for initial conditions associated with (17) as (32)
So, the corresponding solutions are given as follows: (33) (34) (35) (36)
3.2 Solution of nonhomogeneous part
Now, let us consider the nonhomogeneous problem given by (16). Expressing the function v(x, t) as an infinite series, the problem can be rewritten using the ADM as (37)
The components v0(x, t) are identified as and, using the following recursive scheme, it is possible to obtain: (38)
Using the eigenfunction expansion method, the function Q(x, t) can be expressed as a multiplication of unknown functions hj(t) and orthogonal functions ϕj(x), namely, (39) where the unknown functions hj(t) can be represented using orthogonality as (40)
Finally, using the identities for series expansion of sine functions, the solution of the nonhomogeneous equation in (16) can be simplified as (42)
Let us choose the source term in (1) as (43) and then, for fixed values of the nonhomogeneity parameter α = 1 and μ0 = 1, it can be obtained that (44)
Using these values in (42), the solutions for the nonhomogeneous equation are obtained as
Finally, the general solution of the functionally graded clamped simply-supported beam problem can be written as a superposition of homogeneous and nonhomogeneous solutions as
4 Results and discussion
This study considered a specific functionally graded composite beam which has a clamped end and a simply supported end as shown in Figure 1. The solution of the problem is obtained as an expansion of convergent series by the Adomian decomposition method. The deflection w(x, t) was represented by the first six terms of the series along with the corresponding eigenvalues and eigenfunctions. Natural frequencies for different nonhomogeneity parameters are shown in Table 1. Results for sufficiently small nonhomogeneity parameters (like α = 0:001) are compared with the exact solution of the homogeneous beam problem subject to the same boundary conditions given in Appendix A. It is shown that eigenvalues of corresponding eigenfunctions agreed with the exact solution for the beam made of homogeneous material. The figures show how the deflection is changing on the composite beam in terms of the effect of the nonhomogeneity parameter α. For different values of nonhomogeneity parameters – deflections u(x, t) and v(x, t) – solutions of the homogeneous and nonhomogeneous parts of the beam problem are shown in Figure 2 and Figure 3, respectively. The resultant of the deflections w(x, t) is given in Figure 4. In both cases, as the nonhomogeneity parameter α increases, the deflection decreases through the simply supported end for fixed values of time.
The aim of the present work is to examine the behavior of a functionally graded composite beam. Using the Adomian decomposition method, deflections of the composite beam for various values of the nonhomogeneity parameter α are obtained and the results for small values closely match the results for a homogeneous beam.
Appendix: An exact solution for a uniform beam
An exact solution of the vibration problem can be obtained by using the eigenfunction expansion method. In the solution process, the problem will be separated into two parts: where u(x, t) and v(x, t) are the solutions of homogeneous and nonhomogeneous equations, respectively. Setting q(x, t) = 0 in the solution of (46), the function u(x, t) can be defined as (48)
For a uniform beam, μ(x) and EI(x) can be expressed as constants. So, equation (48) becomes (49)
To solve this PDE, we use separation of variables. This method is based on the assumption that the solution u(x, t) can be expressed as the product of ϕ(x) and ψ(t) such that (50)
Eigenfunctions related to the homogeneous problem are given by (51) where ω2 is the separation constant of the system and λ is defined by (52)
The solution of equation (51) will provide ϕ(x) of the beam. The solution of (51) can be found as (53) and the unknown values of C1, C2, C3, C4, can be determined by the end conditions of the beam. An infinite number of nontrivial solutions can be found, corresponding to an infinite number of natural frequencies λi, i = 1, 2, 3, ..., such that (54) where the first four roots of the equation are given as, (55)
Assuming C1 = 1 and using natural frequencies λi, the corresponding nodal shapes ϕi(x) can be expressed. Let us apply the method of eigenfunction expansion by defining ψn(t) as an unknown function such as (56)
Using (51) we write,
Using orthogonality properties of eigenfunctions, (58)
Thus, we obtain the differential equation (59) where
The solution of this ODE, given in (59), can be written as
So the general solution is
Now, we should determine the initial conditions to find Cn and Dn. Using initial conditions, we have (60)
Finally, the nonhomogeneous PDE can be determined as
Adomian G., Nonlinear stochastic systems and applications to physics, Kluwer, Boston, 1989.
Adomian G., Solving frontier problems of physics: The decomposition method, Kluwer, Boston, 1994.
Cherruault Y., Convergence of Adomian’s method, Kybernetes, 1989, 18, 31-38.
Lesnic D., Convergence of ADM: Periodic temperature Comput. Math. Appl., 2002, 44, 13–24.
Wazwaz A. M., Partial differential equations methods and applications, A.A. Belkema Publishers, Lisse, 2002.
Pourdavish A., A reliable symbolic implementation of algorithm for calculating Adomian polynomials, Appl. Math. And Comput., 2006, 172, 545-550.
Zhu Y., Chang Q., and Wu S., A new algorithm for calculating Adomian polynomials, Appl. Math. And Comput., 2005, 169, 402-416.
Wazwaz A. M., A reliable modification of Adomian decomposition method, Appl. Math. Comput., 1999, 102, 72-86.
Lesnic D., The Cauchy problem for the wave equation using the decomposition method, Appl. Math. Letters, 2002, 15, 697-701.
Reddy J. N., An introduction to finite element method, Mc-Graw Hill, New York, 2006.
Haddadpour H., An exact solution for variable coefficients fourth-order wave equation using the Adomian method, Math. And Comput. Modelling, 2006, 44, 1144-1152.
Miyamoto M., Kaysser W. A, Rabin B. H., Kawasaki A. and Ford R. G., Functionally graded materials: Design, processing and applications, Kluwer Academic Publishers, Norwell, Massachusetts, 1999.
Pan W., Gong J., Zhong L. and Chen L., Functionally graded materials VII, Trans Tech Publications LTD, Switzerland, 2003.
Suresh S. and Mortensens A., Fundamentals of functionally graded materials: Processing and thermomechanical behaviour of graded metals and metal-ceramic composites, Maney Publishing, 1998.
Sahin A., Karatay I., Diffusion in graded materials by decomposition method, Scientific Research and Essays, 2011, 6 (14), 2988-2997.
Atangana A., Koca I., Chaos in a simple nonlinear system with Atangana-Baleanu derivatives with fractional order, Chaos, Solitons and Fractals, 2016, 89, 447-454.
Atangana A., On the new fractional derivative and application to nonlinear Fisher’s reaction-diffusion equation, Appl Math Comput, 2016, 1 (273),948-956.
Alkahtani B. S. T., Atangana A., Controlling the wave movement on the surface of shallow water with the Caputo-Fabrizio derivative with fractional order, Chaos, Solitons and Fractals, 2016, 89, 539-546.
Alkahtani B. S. T., Atangana A., Analysis of non-homogeneous heat model with new trend of derivative with fractional order, Chaos, Solitons and Fractals, 2016, 89, 566-51.
Yildirim K., Kucuk I., Optimal Piezoelectric Vibration Control of a Timoshenko Beam, Journal of Franklin Institute,2016, 353(1), 95-107.
Yildirim K., Optimal Forced Vibrations Control of a Smart Plate, Applied Mathematical Modelling, 2016, 40, 6424-6436.
Wazwaz A. M., Linear and Nonlinear Integral Equations: Methods and Applications, HEP and Springer, Beijing, Berlin, 2011.
About the article
Published Online: 2016-09-30
Published in Print: 2016-01-01
© 2016 Yiğit et al., published by De Gruyter Open. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License. (CC BY-NC-ND 3.0)