Application of ADM Using Laplace Transform to Approximate Solutions of Nonlinear Deformation for Cantilever Beam

We investigate semianalytical solutions of Euler-Bernoulli beam equation by using Laplace transform and Adomian decomposition method (LADM). The deformation of a uniform flexible cantilever beam is formulated to initial value problems. We separate the problems into 2 cases: integer order for small deformation and fractional order for large deformation. The numerical results show the approximated solutions of deflection curve, moment diagram, and shear diagram of the presented method.


Introduction
The Euler-Bernoulli beam theory states that the action load produces the bending moment () ∈ ([0, ]) which is proportional to deflection characteristics of the beam.The equation of this law can be written as follows: where  ∈  2 ([0, ]) is deflection curve of a uniform beam, the modulus of elasticity , and the moment of inertia .We note that  and  are both constant and the product of  and  is called beam stiffness.In a case of small deformation, we assume that   () is infinitesimal.Equation ( 1) is reduced to the well-known fourth-order linear differential equation: In this study we consider the uniform flexible of cantilever beam; see Figure 1.The parameters  and Δ are undeformed length and horizontal displacement, respectively.The deformed length of beam is verified by the integral ∫ where () = .If the slope is very small, the linear Euler-Bernoulli beam theory [2] governs the problem The solution of the problem (4) is Our method proposes to reform problem (3) in a sense of fractional calculus without linearization.Because fractional calculus is the great idea to describe behavior in nature, for example, force response of viscoelastic material [3,4], fluid flow [5], and fitting experimental data [6].For convenience, we set  = −/ and  = (0)/[1 −  2 (0)] 1/2 .We develop (1) to deal with fractional order  ∈ (1, 2]; system (3) becomes We use Adomian polynomial to approximate a nonlinear term and derive a semianalytical solution by use of Laplace transform for the initial value problem (6).

Preliminaries
In this section we introduce some definitions and theory of fractional calculus and Laplace transform which are used in our method.See [7] for more details.
Definition 1.The Riemann-Liouville fractional integral operator of order  > 0 is defined as Definition 2. The Caputo fractional derivative of order  > 0 is defined as Theorem 3. If  > 0 and  > 0, then (

Theorem 4. The Laplace transform of the Caputo fractional derivative is given by
where  = ⌊⌋ and  is a nonnegative integer.The inverse Laplace transform of power function is given by

Semianalytical Solution via LADM
To formulate the general solution of problem ( 6), we replace   with  in the equation of problem ( 6) and applied Theorem 3. We obtain Taking Laplace transform of (11) gives Using Theorem 4 and replacing   (0) by , we can rewrite as We take inverse Laplace transform of (13), and the following operator equation is obtained: We apply Adomian polynomial [8,9] to approximate a nonlinear term of (14), by setting Equation ( 14) becomes Assume that the solution of (11) can be written as  = ∑ ∞ =1   .From the combination of Laplace transform and Adomian decomposition method, we obtain the formulation of LADM for fractional order of cantilever beam deflection equation: In this work, we apply the first four terms of LADM which are provided from following recurrence: where The components of adomian polynomial are given by The iteration results of recurrence (18) are shown as follows: The first four terms' approximate solution is  =  0 + 1 + 2 +  3 : Since  =   and (0) = 0, thus () = ∫  0 ().This yields (24) 3.1.Integer Order for Small Deformation.We simulate some deflection curve of the results in (24) and ( 5) which are solution of LADM and classical method, respectively.The parameters are  = 100 in., Δ = 0 in.,  = 1.8 × 10 5 kip in. 2 ,  = 2, and  = 1kip.Figure 2 shows the deflection of cantilever beam, obtained by using classical method and LADM.LADM has accurate slope around the free end of the beam as well as classical method.Moreover LADM shows that the deflection curve remains nearly straight line in the case of small deformation.The effects of various loads influencing maximum deflection at the free end and length of the deformed beam   1.We observe that the calculation of nonlinear integer order deformation via LADM fails at  = 50kips and  = 100 kips.Since (0) > 1 thus integral of bending moment is greater than stiffness of the beam.This yields large deformation, large rotation, and large strain, namely, some concept of large deformation.Not only does LADM fails, but also classical method does, because it contradicts the assumption that the lengths of deformed and undeformed beam are identical.This is the motivation for reform in fractional order model.

Fractional Order for Large Deformation.
There are many methods for describing a large deformation of a beam, for example, the Laplace-Padé coupling with NDHPM and HPM [10], VIM [11], and pseudolinear system [1].In this section, we introduce the process of investigating semianalytical solution via LADM and some numerical results of pseudolinear system in [1].
To illustrate the numerical results we suppose that horizontal displacement Δ = 40.The beam parameters are  = 1000 in.and  = 1.8 × 10 5 kip in. 2 .The fact from our assumption is In this case we fixed the concentrate load  = 0.4 kip.We get the deformed length  = 840 and initial condition Table 2 shows that order  = 1.54 is given,   (840) = −0.0238and   (829.711)= 0, which is the most accurate for the fact in (25).We design order  = 1.54 for this problem.
Using the first four terms' approximation, the deflection is shown in Figure 3(a).The slope deflection of psudolinear system and LADM are presented in Table 3. Moment diagram (  ) can be computed from the deflection curve as follows:

Conclusion
In this study, we use Euler-Bernoulli beam equation for describing the uniform flexible cantilever beam with a concentrated load.The initial conditions are given by calculating slope of the beam.We use LADM to determine the semianalytical solution.LADM with integer order system can approximate solution without cancellation a nonlinear term in the case of small deformation.For the large deformation, we reform the problem to fractional order system and estimate fractional order  which conserves the fact from the assumption.Finally, we show that LADM gives the solution as a polynomial expression which is the advantage for analyzing moment and shear diagrams.LADM may be a powerful and successive method for solving nonlinear science and engineering problem.

Figure 1 :
Figure 1: Cantilever beam with a concentrated load at the free end.
2 and () =    () for shear diagram.The results of both moment and shear diagram are shown in Figures3(b) and 3(c), respectively.

Table 1 :
Numerical results under various load testing.

Table 2 :
Values of slope at 840 in.and maximum vertical displacement for various fractional orders.

Table 3 :
Numerical results of slope deflection.
=   (0) = 1.2626.Fractional order is determined by simulating the results of slope at 840 and deformed length.