Abstract: In this research, variational iteration method and homotopy perturbation method are applied to solve a nonlinear fourth order boundary value problem. These problems used as mathematical models in viscoelastic inelastic flows and deformation of beams and plate deflection theory. Comparison is made between the exact solutions and the results of the Variational Iteration Method (VIM) and Homotopy Perturbation Method (HPM). The results reveal that these methods are very effective and simple. In this survey, it will be shown that these methods can also be used for solving nonlinear boundary value problems.
INTRODUCTION
According to the classical beam theory, if u = u(x) denotes the configuration of the deformed beam, then the bending moment satisfies the relation M = -Eiu″ where, E is the Young modulus of elasticity and I is the inertial moment. Now if the deformation is caused by a load f = f(x), one deduces, from a free body diagram, that f = -v and v = M′ = -Eiu′′′, where, v denotes the shear force. Suppose that u represents an elastic beam of length L = 1, which is clamped at its left side x = 0 and resting on a kind of elastic bearing at its right side x =1, Along its length, a load f is added to cause deformations (Fig. 1).
For simplicity it can be assumed EI = 1, then from above remarks, it can be got the following boundary value problem:
(1) |
(2) |
(3) |
where,
Fig. 1: | Beam on elastic bearing |
Chawla and Katti (1979) investigated finite difference scheme for the nonlinear differential equation of order 2n:
(4) |
Subject to the boundary conditions:
(5) |
In this study, we consider a general fourth-order boundary value problem of the form:
(6) |
With the boundary conditions:
(7) |
where, f is a continuous function on [a, b] and the parameters αi and βi, i = 1, 2 are finite real arbitrary constants. Such type of systems are not only regarded as a general boundary value problem, but also used as mathematical models in viscoelastic and inelastic flows (Momani, 1991), deformation of beams (Ma and Silva, 2004) and plate deflection theory (Chawla and Katti, 1979). With the rapid development of nonlinear science, many different methods were proposed to solve various boundary-value problems (BVPS), such as the homotopy perturbation method (Rafei and Ganji, 2006; Ganji and Sadighi, 2006; He, 1999a, 2000, 2003; Zhang and He, 2006; Choobbasti et al., 2008, Barari et al., 2008) and the variational iteration method (VIM) (He, 1999b; He and Wu, 2006; Tari et al., 2007; Momani and Abuasad, 2006; Odibat and Momani, 2006). In this letter, we apply the homotopy-perturbation method and variational iteration method to the discussed problem.
BASIC IDEA OF HOMOTOPY-PERTURBATION METHOD
Linear and Nonlinear phenomena are of fundamental importance in various fields of science and engineering. Most models of real-life problems are still very difficult to solve. Therefore, approximate analytical solutions such as Homotopy-Perturbation Method (HPM) were introduced. This method is the most effective and convenient ones for both linear and nonlinear equations.
Perturbation method is based on assuming a small parameter. The majority of nonlinear problems, especially those having strong nonlinearity, have no small parameters at all and the approximate solutions obtained by the perturbation methods, in most cases, are valid only for small values of the small parameter.
Generally, the perturbation solutions are uniformly valid as long as a scientific system parameter is small. However, we cannot rely fully on the approximations, because there is no criterion on which the small parameter should exists. Thus, it is essential to check the validity of the approximations numerically and/or experimentally. To overcome these difficulties, HPM have been proposed recently.
To explain this method, let us consider the following function:
(8) |
With the boundary conditions of:
(9) |
where, A, B, f(r) and Γ are a general differential operator, a boundary operator, a known analytical function and the boundary of the domain Ω, respectively.
Generally speaking the operator A can be divided into a linear part Land a nonlinear part N (u). Equation 8 can, therefore, be written as:
(10) |
By the homotopy technique, we construct a homotopy
(11) |
or
(12) |
where, p∈[0,1] is an embedding parameter, while u0 is an initial approximation of Eq. 8, which satisfies the boundary conditions. Obviously, from Eq. 11 and 12 we will have:
(13) |
(14) |
The changing process of p from zero to unity is just that of v (r, p) from u0 to u (r). In topology, this is called deformation, while L(v) -L(u0) and A(v)-f(r) are called homotopy.
According to the HPM, we can first use the embedding parameter p as a small parameter and assume that the solutions of Eq. 11 and 12 can be written as a power series in p:
(15) |
Setting p = 1 yields in the approximate solution of Equation 15 to:
(16) |
The combination of the perturbation method and the homotopy method is called the HPM, which eliminates the drawbacks of the traditional perturbation methods while keeping all its advantages.
The series (16) is convergent for most cases. However, the convergent rate depends on the nonlinear operator A(v). Moreover, He (1999a) made the following suggestions:
• | The second derivative of N(v) with respect to v must be small because the parameter may be relatively large, i.e., p→1. |
• | The norm of |
BASIC IDEA OF VARIATIONAL ITERATION METHOD
To clarify the basic ideas of VIM, we consider the following differential equation:
(17) |
where, L is a linear operator, N is a nonlinear operator and g (t) is a homogeneous term.
According to VIM, we can write down a correction functional as follows:
(18) |
where, λ is a general lagrangian multiplier which can be identified optimally
via the variational theory. The subscript n indicates the nth approximation
and un is considered as a restricted variation, i.e.,
APPLICATION OF VARIATIONAL ITERATION METHOD
Consider the following nonlinear boundary value problem:
(19) |
Subject to the boundary conditions:
(20) |
where,
(21) |
The exact solution for this problem is:
(22) |
According to Eq. 18, we have the following iteration formulation:
(23) |
Now we assume that an initial approximation has the form:
(24) |
where, a, b, c and d are unknown constants to be further determined.
By the iteration formula (23), we have the following first-order approximation:
(25) |
Incorporating the boundary conditions, Eq. 20, into u1(x), we obtain:
(26) |
We therefore, obtain the following first-order approximate solution:
(27) |
APPLICATION OF HOMOTOPY-PERTURBATION METHOD
To solve Eq. 19 by means of HPM, we consider the following process after separating the linear and nonlinear parts of the equation.
A homotopy can be constructed as follows:
(28) |
Substituting v = v0 + pv1+ in to Eq. 28 and rearranging the resultant equation based on powers of p-terms, one has:
(29) |
(30) |
(31) |
Table 1: | Comparison of the approximate solutions with exact solution |
With the following conditions:
(32) |
With the effective initial approximation for v0 from the conditions (32) and solutions of Eq. 29-31 may be written as follows:
(33) |
(34) |
(35) |
In the same manner, the rest of components were obtained using the maple package.
According to the HPM, we can conclude that:
(36) |
Fig. 2: | Comparison between different solutions |
Therefore, substituting the values of v0(x), v1(x) and v2 (x) from Eq. 33-35 into Eq. 36 yields:
(37) |
Comparison of the approximate solutions with exact solution is shown in Table 1 and Fig. 2 showing a remarkable agreement. Of course we can obtain even higher accurate solutions without any difficulty.
CONCLUSIONS
The homotopy perturbation method and variational iteration method are employed successfully to study a fourth order boundary value problem in structural engineering and fluid mechanic. The results revealed that The variational iteration method and homotopy perturbation method are remarkably effective for solving boundary value problems. Comparison between the approximate and exact solutions shows that the one iteration of variational iteration method is enough. These methods are very promoting method, which will be certainly found widely applications.