Research Article  Open Access
Approximations for Large Deflection of a Cantilever Beam under a Terminal Follower Force and Nonlinear Pendulum
Abstract
In theoretical mechanics field, solution methods for nonlinear differential equations are very important because many problems are modelled using such equations. In particular, large deflection of a cantilever beam under a terminal follower force and nonlinear pendulum problem can be described by the same nonlinear differential equation. Therefore, in this work, we propose some approximate solutions for both problems using nonlinearities distribution homotopy perturbation method, homotopy perturbation method, and combinations with LaplacePadé posttreatment. We will show the high accuracy of the proposed cantilever solutions, which are in good agreement with other reported solutions. Finally, for the pendulum case, the proposed approximation was useful to predict, accurately, the period for an angle up to yielding a relative error of 0.01222747.
1. Introduction
Solving nonlinear differential equations is an important issue in science because many physical phenomena are modelled using such equations. During last century, the approximate solution of such equations was a task performed by hand using methods like the perturbation method. Nowadays, symbolic software like Maple or Mathematica allow researchers to calculate highly accurate approximate solutions using new methods like the homotopy perturbation method (HPM) [1–38]. The HPM method is one of the most famous analytic techniques for nonlinear differential equations, which is widely applied in science and engineering.
In [28, 39] was reported the nonlinearities distribution homotopy perturbation method (NDHPM) as an extended version of HPM that allows the distribution of the nonlinearities between the different iterations of HPM method. The main advantage of this process is the increase in the number of possible iterations by constructing easier to solve linear equations for HPM procedure. Additionally, several works reported that a LaplacePadé posttreatment [40–49] of power series solutions obtained by HPM can improve accuracy and convergence to the exact solution. Therefore, in this work, we propose the use of the aforementioned methods to obtain approximate solutions for two nonlinear problems, related by to same nonlinear differential equation. The first one is the large deflection of a cantilever beam under a terminal follower force [50–55]. This kind of structures arises in marine riser and other practical applications of civil engineering. The mathematical formulation of such a problem yields a nonlinear twopoint boundaryvalue problem, which usually can be solved by numerical methods [50]. The second one is the nonlinear pendulum problem [56–58], which is modelled by the same equation of the Cantilever’s case study. The pendulum or systems containing the pendulums are wellstudied complex phenomen [59]. For both problems, the results exhibited high accuracy using our proposed solutions, which are in good agreement to the results reported by other works. In particular, for the pendulum case, we succeeded to predict, accurately, the period until an initial angle up to achieving a relative error of .
This paper is organized as follows. In Section 2, the basic idea of the HPM method is provided. We give an introduction to NDHPM method in Section 3. In Section 4, the basic concept of Padé approximants is explained. The LaplacePadé coupling with NDHPM and HPM methods is presented in Section 5. In Sections 6 and 7, the cantilever and pendulum problems are solved, respectively. In addition, a discussion on the results is presented in Section 8. Finally, a brief conclusion is given in Section 9.
2. Basic Concept of HPM
It can be considered that a nonlinear differential equation can be expressed as having boundary condition as where is a general differential operator, is a known analytic function, is a boundary operator, is the boundary of domain , and denotes differentiation along the normal drawn outwards from [50]. The operator, generally, can be divided into two operators, and , which are linear and nonlinear operators, respectively. Hence, (1) can be rewritten as
Now, a possible homotopy formulation is where is the initial approximation for (3) which satisfies the boundary conditions and is known as the perturbation homotopy parameter.
For the HPM method [4–7], we assume that the solution for (4) can be written as a power series of
Considering that , results that the approximate solution of (1) is
The series (6) is convergent on most cases as reported in [4, 7, 20, 21].
3. HPM Method with Nonlinearities Distribution
Recent reports [28, 39] have introduced the NDHPM method, which eases the searching process of solutions for (3) and reduces the complexity of solving differential equations. As first step, a modified homotopy is introduced:
It can be noticed that the homotopy function (7) is essentially the same as (4), except for the nonlinear operator and the nonhomogeneous function , which embeds the homotopy parameter . The arbitrary introduction of within the differential equation is a strategy to redistribute the nonlinearities between the successive iterations of the HPM method and, thus, increase the probabilities of finding the sought solution.
Again, we establish that
Considering that turns out that the approximate solution for (1) is
The convergence of the NDHPM method is exposed in [28].
4. Padé Approximants
A rational approximation to on is the quotient of two polynomials and of degrees and , respectively. We use the notation to denote this quotient. The Padé approximations to a function are given by [49, 60]
The method of Padé requires that and its derivative should be continuous at . The polynomials used in (10) are
The polynomials in (11) are constructed so that and agree at and their derivatives up to agree at . For the case where , the approximation is just the Maclaurin expansion for . For a fixed value of the error is the smallest when and have the same degree or when has a degree higher than .
Notice that the constant coefficient of is . This is permissible because it can be noted that and are not changed when both and are divided by the same constant. Hence, the rational function has unknown coefficients. Assume that is analytic and has the Maclaurin expansion
And, from the difference
The lower index in the summation at the right side of (13) is chosen because the first derivatives of and should agree at .
When the left side of (13) is multiplied out and coefficients of the powers of are set equal to zero for , this results in a system of linear equations:
Notice that in each equation the sum of the subscripts on the factors of each product is the same, and this sum increases consecutively from to . The equations in (15) involve only the unknowns and must be solved first. Then, the equations in (14) are successively used to find [49].
5. LaplacePadé Posttreatment of Power Series Solutions
The coupling of Laplace transform and Padé approximant [40] is used to deal with the lost information of truncated power series [41–49]. The process can be summarized as follows.(1)First, Laplace transformation is applied to a power series solution obtained by HPM or NDHPM methods. (2)Next, is written in place of . (3)Afterwards, we convert the transformed series into a meromorphic function by forming its Padé approximant of order . and are arbitrarily choses, but they should be smaller than the power series order. In this step, the Padé approximant extends the domain of the truncated series solution to obtain better accuracy and convergence. (4)Then, is written in place of . (5)Finally, by using the inverse Laplace transformation, we obtain the modified approximate solution.
The LaplacePadé posttreatment of HPM and NDHPM methods will be referred throughout the rest of this paper as LPHPM and LPNDHPM, respectively.
6. Solution of Cantilever Problem
Beams subjected to follower load arise in many problems of structural engineering like lumbar spine in biomechanics, smart structures applications, and loads applied to structural components by cables in tension, among others. The governing equations of such structures involve the effects of nonlinearities due to large deformations and material properties [50–55].
In this work, we study the largedeflection problem of a cantilever beam under a tipconcentrated follower force as depicted in Figure 1. The angle of inclination of the force with respect to the deformed axis of the beam remains unchanged during deformation. We assume that the center of the Cartesian system is located at the fixed end and is the local angle of inclination of the centroidal axis beam. Therefore, the mathematical formulation of this problem yields the following nonlinear twopoint boundaryvalue problem: where prime denotes differentiation with respect to . Additionally, the related dimensionless quantities are defined as where represents the length of the originally straight elastic cantilever beam, is the flexural rigidity, is the arc length measured from the tip, and is the terminal follower force with constant angle .
In [50] is reported that (16) can be reduced to an initialvalue problem by the following transformation: resulting in the following initial value problem: with the extra condition: where the slope of can be obtained by using (18) and (20), resulting in the following:
Finally, the deformed curve of beam axis can be calculated from solution of the slope as
6.1. Solution Obtained by Using NDHPM Method
In order to circumvent difficulties due to the nonlinear term of (19), we propose to use a seventhorder Taylor series expansion, resulting in the following: and by using (7) and (23), we establish the following homotopy equation: where homotopy parameter has been arbitrarily embedded into the power series of Taylor’s expansion, and linear operator is and trial function is
Substituting (8) into (24) and reordering and equating terms with the same powers, we obtain the following linear set of differential equations: By solving (27), we obtain where .
Substituting solutions (28) into (8) and calculating the limit when results in the thirdorder approximation:
6.2. Solution Calculated by Using HPM and LPHPM Methods
By using (4) and (23), we establish the following homotopy equation: where linear operator is and trial function is
Substituting (5) into (30) and reordering and equating terms with the same powers, we obtain the following set of linear differential equations: By solving (33), we obtain where the subsequent iterations are calculated using Maple CAS software.
Substituting solutions (34) into (5) and calculating the limit when , we can obtain the 18thorder approximation:
Expression (35) is too long to be written here; nevertheless, it can be simplified if we consider a particular case. For instance, when , it results in the following:
By using (36) and (21) we can calculate, approximately, the deflection angle ; nevertheless, accuracy can be increased applying the LaplacePadé posttreatment. First, Laplace transformation is applied to (36); then, is written in place of . Afterwards, Padé approximant is applied and is written in place of . Finally, by using the inverse Laplace transformation, we obtain the modified approximate solution. For instance, for , the result of LPHPM is
On one side, for , the exact result is [51]. On the other side, by using (21), we can calculate the deflection angle using (36) and (37) approximations, resulting in and , respectively. Therefore, for this case, the LaplacePadé posttreatment was relevant in order to deal with the HPM truncated power series, resulting a notorious increase of accuracy (see Table 1).

6.3. Power Series Solution by Using NDHPM and LaplacePadé Posttreatment
By using (7) and (23), we establish the following homotopy equation: where homotopy parameter has been arbitrarily embedded into the power series of Taylor expansion, linear operator , and trial function are taken from (31) and (32), respectively.
Substituting (8) into (38), reordering and equating terms with the same powers, we obtain the following set of linear differential equations:
Solving (39), we obtain where the subsequent iterations are calculated using Maple CAS software.
Substituting solutions (40) into (8) and calculating the limit when , we can obtain the 18thorder approximation:
Expression (41) is too long to be written here; nevertheless, it can be simplified considering a particular case. For instance, a terminal follower normal angle , resulting
In order to increase accuracy of , (42) is transformed by the Padé approximation and Laplace transformation. First, Laplace transformation is applied to (42) and then is written in place of . Then, Padé approximant is applied and is written in place of . Finally, by using the inverse Laplace transformation, we obtain the modified approximate solution. For example, given , the resultant LPNDHPM approximation is
Now, by using (21), we can calculate for (42) and (43) approximations. Results for NDHPM and LPNDHPM are and , respectively (see Table 1). In fact, for , the exact result is ; therefore, the LaplacePadé aftertreatment was a key factor to increase accuracy of NDHPM truncated power series.
7. Solution of Nonlinear Pendulum
The pendulum [56–58] of Figure 2 can be modelled by (19). Nevertheless, we need to redefine the variables and parameters of (19), resulting in the following: considering where prime denotes differentiation with respect to time , is the length of the pendulum, is acceleration due to gravity, and is the initial angle of displacement .
Taking into account the above considerations and substituting into (29), results that the approximate solution (29) obtained by using NDHPM can be considered, also, as the solution of the nonlinear pendulum problem (44).
7.1. Solution Obtained by Using HPM and LaplacePadé Posttreatment
Although approximate solutions (35) and (41) possess good accuracy for cantilever problem, we need to increase the order of Taylor expansion for the term in (23) to obtain a highly accurate solution for pendulum problem, compared to other reported solutions [50–55]. Therefore, we propose to use a 22thorder Taylor series expansion of term, resulting in the following: where initial conditions are and .
A homotopy formulation that can generate a power series solution for this problem is where linear operator and trial function are taken from (31) and (32), respectively.
Substituting (5) into (47) and reordering and equating terms with the same powers, we obtain the following set of linear differential equations: by solving (48), we obtain where the subsequent iterations are calculated by using Maple CAS software.
Substituting solutions (49) into (5) and calculating the limit when , we can obtain the 12thorder approximation:
Expression (50) is too large to be written here; nevertheless, it can be simplified considering particular cases. For example, if we consider and , we get the following results:
Now, (51) is transformed by the LaplacePadé posttreatment. First, Laplace transformation is applied to (51) and then is written in place of . Then, Padé approximant is applied and is written in place of . Finally, by using the inverse Laplace transformation, we obtain the modified approximate LPHPM solution:
Equation (52) is more compact and computationally efficient than power series solution (51). Thus, (52) can be used when accurate and handy expressions are needed. This process can be repeated for different values of and .
We can express pendulum period in terms of where is the complete elliptic function of the first kind formulated as
For comparison purposes, we will normalize the value of to
Finally, is calculated for some values of in Tables 3 and 4.
8. Numerical Simulation and Discussion
8.1. Cantilever Problem
Structures subjected to follower forces have been reported by several researchers [50, 55]. In this work, we propose the solution for large deflections of a cantilever beam subjected to a follower force reported in [55]. Therefore, by using approximate solution (29) and (21), we show in Figure 3 the relationship between terminal slope and load parameter for some values of deflection angle in the range . Flutter instability occurs, and the maximum value of terminal slope is approximately as reported in [50, 55].
On one side, considering a fixed value for , we show in Figure 4 the deformed curves of a beam under different values of load parameter . First, the cantilever beam is a straight line, and it is deformed gradually as the amplitude of the endconcentrated follower force increases, until the curvature of the beam is notorious for . Additionally, when the load parameter is larger than the flutter load, the beam deforms towards the opposite direction as reported in [50]. On the other side, we consider a fixed value for load parameter and different values of as shown in Figure 5, resulting that the deflection and terminal slope are influenced by the angle of follower force producing a notorious bending as increases reaching the extreme case at .
Considering a normal angle , we show in Table 1 the terminal slope (using (21)) versus different values of load parameter . In such table, we compare the proposed solutions (29), (35), (41), LPHPM/LPNDHPM approximations (see (37), and (43)) to exact solution [51], a shooting method solution [53], and HPM solution reported in [50]. Our proposed approximate solutions are in good agreement with the other reported approximations. In particular, (29) and LPNDHPM approximations (see (43)) exhibited the best accuracy of all proposed solutions of this work, with a similar accuracy compared to the results reported in [50]. Furthermore, we showed in Table 2, the terminal slope and coordinates (see (22)) of cantilever beam under normal follower forces (), resulting that our proposed solution (29) exhibited an accuracy in good agreement with the values obtained by a direct method (D.M.) [55] and the results reported in [50]. For this work, the cantilever beam is considered elastic and ideal; then, Figures 4 and 5 present an ideally deformed beam. Further work should be performed in order to include nonideal effects due to material characteristics of cantilever.



From above, we can conclude that the deflection and terminal slope are governed by the deflection angle and magnitude of the follower force. This means that the bending of the beam becomes more notorious as the deflection angle or the magnitude of the follower force increases [50].
8.2. Pendulum Problem
On one hand, in [56, 58] the authors achieve approximate solutions for pendulum period, with relative error <2% for , and <1% for , respectively. On the other hand, in Table 3 and Table 4 we show a comparison between , calculated using (29) and (50) and approximations calculated by LPHPM procedure explained in Section 7.1 (considering ). The value is obtained, applying NewtonRaphson method to equation , that is, calculating the first crossing of by zero since the initial angle . In Table 4, we show the relative error (RE) for the range of , resulting that our proposed solutions are in good agreement with exact results. Furthermore, LPHPM solution can predict the for exhibiting a RE of , which is highly accurate compared to the results obtained in [56, 58].
The coupling of LaplacePadé and HPM method was a key factor to increase the accuracy range until value is really close to . LaplacePadé posttreatment was able to deal with the truncate power series, in order to obtain better accuracy for angles close to although the behaviour was different for small values of . However, in general terms, the LaplacePadé posttreatment generates handy and computationally more efficient expressions (see (51) versus (52)). Furthermore, we can appreciate from Tables 3 and 4 that HPM and NDHPM methods do not present useful values for and , due to numerical noise of the proposed approximations. Fortunately, LPHPM approximations achieve a good accuracy within the range . From the above discussion we can conclude that, as the initial angle of displacement approximates to the vicinity , the pendulum period increases notoriously.
Finally, further research should be done in order to obtain approximate solutions of the pendulum problem involving additional nonlinear effects like friction and mass, among others.
9. Conclusions
In this paper we presented the solution for two nonlinear problems related by the same differential equation. The first one was the large deflection of a cantilever beam under a terminal follower force. The second one was the nonlinear pendulum problem. We proposed a series of solutions obtained by NDHPM, HPM, and combinations with LaplacePadé posttreatment. For both problems, results exhibited high accuracy, which are in good agreement with the results reported by other works. In particular, for the pendulum case, we succeeded to predict, accurately, the period for an initial angle up to . By means of two case studies, we confirmed that NDHPM, HPM, LPNDHPM, and LPHPM are powerful and useful methods; all of them capable of generating highly accurate approximations for nonlinear problems.
Acknowledgments
The authors gratefully acknowledge the financial support provided by the National Council for Science and Technology of Mexico (CONACyT) through Grant CB201001 no. 157024. The authors would like to express their gratitude to RogelioAlejandro CallejasMolina and Roberto RuizGomez for their contribution to this paper.
References
 J.H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012. View at: Publisher Site  Google Scholar
 J.H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “An elementary introduction to the homotopy perturbation method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 410–412, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “The homotopy perturbation method nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “A coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005. View at: Publisher Site  Google Scholar
 J.H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 12, pp. 87–88, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. D. Ganji, “A semiAnalytical technique for nonlinear settling particle equation of Motion,” Journal of HydroEnvironment Research, vol. 6, no. 4, pp. 323–327, 2012. View at: Publisher Site  Google Scholar
 M. Sheikholeslami and D. D. Ganji, “Heat transfer of Cuwater nanofluid flow between parallel plates,” Powder Technology, vol. 235, pp. 873–879, 2013. View at: Publisher Site  Google Scholar
 S. S. Ganji, D. D. Ganji, M. G. Sfahani, and S. Karimpour, “Application of AFF and HPM to the systems of strongly nonlinear oscillation,” Current Applied Physics, vol. 10, no. 5, pp. 1317–1325, 2010. View at: Publisher Site  Google Scholar
 D. D. Ganji, H. B. Rokni, M. G. Sfahani, and S. S. Ganji, “Approximate traveling wave solutions for coupled WhithamBroerKaup shallow water,” Advances in Engineering Software, vol. 41, no. 78, pp. 956–961, 2010. View at: Publisher Site  Google Scholar
 A. Barari, M. Omidvar, A. R. Ghotbi, and D. D. Ganji, “Application of homotopy perturbation method and variational iteration method to nonlinear oscillator differential equations,” Acta Applicandae Mathematicae, vol. 104, no. 2, pp. 161–171, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Koçak and A. Yildirim, “Numerical solution of 3D Green's function for the dynamic system of anisotropic elasticity,” Physics Letters A, vol. 373, no. 35, pp. 3145–3150, 2009. View at: Publisher Site  Google Scholar
 Y. Khan, H. VazquezLeal, and Q. Wu, “An efficient iterated method for mathematical biology model,” Neural Computing and Applications, pp. 1–6, 2012. View at: Publisher Site  Google Scholar
 Y. Khan, Q. Wu, N. Faraz, A. Yildirim, and M. Madani, “A new fractional analytical approach via a modified RiemannLiouville derivative,” Applied Mathematics Letters, vol. 25, no. 10, pp. 1340–1346, 2012. View at: Publisher Site  Google Scholar
 N. Faraz and Y. Khan, “Analytical solution of electrically conducted rotating flow of a second grade fluid over a shrinking surface,” Ain Shams Engineering Journal, vol. 2, no. 34, pp. 221–226, 2011. View at: Google Scholar
 Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers & Mathematics with Applications, vol. 61, no. 11, pp. 3391–3399, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Biazar and H. Aminikhah, “Study of convergence of homotopy perturbation method for systems of partial differential equations,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2221–2230, 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, vol. 10, no. 5, pp. 2633–2640, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Biazar, F. Badpeima, and F. Azimi, “Application of the homotopy perturbation method to ZakharovKuznetsov equations,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2391–2394, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Biazar and H. Ghazvini, “Homotopy perturbation method for solving hyperbolic partial differential equations,” Computers & Mathematics with Applications, vol. 56, no. 2, pp. 453–458, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Biazar and H. Ghazvini, “Exact solutions for nonlinear Schrödinger equations by He's homotopy perturbation method,” Physics Letters A, vol. 366, no. 12, pp. 79–84, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Fathizadeh, M. Madani, Y. Khan, N. Faraz, A. Yildirim, and S. Tutkun, “An effective modification of the homotopy perturbation method for mhd viscous flow over a stretching sheet,” Journal of King Saud University—Science, 2011. View at: Publisher Site  Google Scholar
 M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and Laplace transformation,” Mathematical and Computer Modelling, vol. 53, no. 910, pp. 1937–1945, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. VázquezLeal, U. FilobelloNiño, R. CastañedaSheissa, A. SarmientoReyes, and J. Sanchez Orea, “High accurate simple approximation of normal distribution integral,” Mathematical Problems in Engineering, vol. 2012, Article ID 124029, 22 pages, 2012. View at: Publisher Site  Google Scholar
 H. VázquezLeal, U. FilobelloNiño, R. CastañedaSheissa, L. HernándezMartínez, and A. SarmientoReyes, “Modified HPMs inspired by homotopy continuation methods,” Mathematical Problems in Engineering, vol. 2012, Article ID 309123, 19 pages, 2012. View at: Google Scholar  MathSciNet
 U. FilobelloNiño, H. VázquezLeal, R. CastañedaSheissa et al., “An approximate solution of blasius equation by using hpm method,” Asian Journal of Mathematics and Statistics, vol. 5, no. 2, pp. 50–59, 2012. View at: Google Scholar
 D. Ağırseven and T. Öziş, “An analytical study for Fisher type equations by using homotopy perturbation method,” Computers & Mathematics with Applications, vol. 60, no. 3, pp. 602–609, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Öziş and D. Ağırseven, “He's homotopy perturbation method for solving heatlike and wavelike equations with variable coefficients,” Physics Letters A, vol. 372, no. 38, pp. 5944–5950, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Öziş and A. Yildirim, “A note on He's homotopy perturbation method for van der Pol oscillator with very strong nonlinearity,” Chaos, Solitons and Fractals, vol. 34, no. 3, pp. 989–991, 2007. View at: Publisher Site  Google Scholar
 D. D. Ganji and A. Rajabi, “Assessment of homotopyperturbation and perturbation methods in heat radiation equations,” International Communications in Heat and Mass Transfer, vol. 33, no. 3, pp. 391–400, 2006. View at: Publisher Site  Google Scholar
 M. Rafei, D. D. Ganji, and H. Daniali, “Solution of the epidemic model by homotopy perturbation method,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 1056–1062, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. D. Ganji, “The application of He's homotopy perturbation method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 355, no. 45, pp. 337–341, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 D. D. Ganji and A. Sadighi, “Application of homotopyperturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24–34, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Sheikholeslami, H. R. Ashorynejad, D. D. Ganji, and A. Yildirim., “Homotopy perturbation method for threedimensional problem of condensation film on inclined rotating disk,” Scientia Iranica, vol. 19, no. 3, pp. 437–442, 2012. View at: Google Scholar
 D. D. Ganji, H. Tari, and M. B. Jooybari, “Variational iteration method and homotopy perturbation method for nonlinear evolution equations,” Computers & Mathematics with Applications, vol. 54, no. 78, pp. 1018–1027, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Homotopy perturbation method with an auxiliary term,” Abstract and Applied Analysis, vol. 2012, Article ID 857612, 7 pages, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Yamamoto, C. Dang, Y. Hao, and Y. C. Jiao, “An aftertreatment technique for improving the accuracy of Adomian's decomposition method,” Computers & Mathematics with Applications, vol. 43, no. 67, pp. 783–798, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. H. Sweilam and M. M. Khader, “Exact solutions of some coupled nonlinear partial differential equations using the homotopy perturbation method,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2134–2141, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Momani, G. H. Erjaee, and M. H. Alnasr, “The modified homotopy perturbation method for solving strongly nonlinear oscillators,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2209–2220, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Khan and N. Faraz, “Application of modified Laplace decomposition method for solving boundary layer equation,” Journal of King Saud University—Science, vol. 23, no. 1, pp. 115–119, 2011. View at: Publisher Site  Google Scholar
 D. Bahuguna, A. Ujlayan, and D. N. Pandey, “A comparative study of numerical methods for solving an integrodifferential equation,” Computers & Mathematics with Applications, vol. 57, no. 9, pp. 1485–1493, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Momani and V. S. Ertürk, “Solutions of nonlinear oscillators by the modified differential transform method,” Computers & Mathematics with Applications, vol. 55, no. 4, pp. 833–842, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 A. Gökdoğan, M. Merdan, and A. Yildirim, “The modified algorithm for the differential transform method to solution of Genesio systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 1, pp. 45–51, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 P.Y. Tsai and C.K. Chen, “An approximate analytic solution of the nonlinear Riccati differential equation,” Journal of the Franklin Institute, vol. 347, no. 10, pp. 1850–1862, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. E. Ebaid, “A reliable aftertreatment for improving the differential transformation method and its application to nonlinear oscillators with fractional nonlinearities,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 1, pp. 528–536, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Merdan, A. Gökdoğan, and A. Yildirim, “On the numerical solution of the model for HIV infection of ${\mathrm{CD}4}^{+}\text{T}$ cells,” Computers & Mathematics with Applications, vol. 62, no. 1, pp. 118–123, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y.G. Wang, W.H. Lin, and N. Liu, “A homotopy perturbationbased method for large deflection of a cantilever beam under a terminal follower force,” International Journal for Computational Methods in Engineering Science and Mechanics, vol. 13, no. 3, pp. 197–201, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 B. N. Rao and G. V. Rao, “On the large deflections of cantilever beams with end rotational load,” Zeitschrift Fur Angewandte Mathematik Und Mechanik, vol. 6, pp. 507–508, 1986. View at: Google Scholar
 B. N. Rao and G. V. Rao, “Large deflections of a nonuniform cantilever beam with end rotational load,” Forschung im Ingenieurwesen, vol. 54, no. 1, pp. 24–26, 1988. View at: Publisher Site  Google Scholar
 B. N. Rao and G. V. Rao, “Large deflections of a cantilever beam subjected to a rotational distributed loading,” Forschung im Ingenieurwesen, vol. 55, no. 4, pp. 116–120, 1989. View at: Publisher Site  Google Scholar
 J. H. Argyris and S Symeonidis, “Nonlinear finite element analysis of elastic systems under nonconservative loading: natural formulation. I. Quasistatic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 26, no. 1, pp. 75–123, 1981. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. S. Shvartsman, “Large deflections of a cantilever beam subjected to a follower force,” Journal of Sound and Vibration, vol. 304, no. 3–5, pp. 969–973, 2007. View at: Publisher Site  Google Scholar
 A. Beléndez, A. Hernández, T. Beléndez, C. Neipp, and A. Márquez, “Application of the homotopy perturbation method to the nonlinear pendulum,” European Journal of Physics, vol. 28, no. 1, article 010, pp. 93–104, 2007. View at: Publisher Site  Google Scholar
 A. Beléndez, E. Arribas, M. Ortuño, S. Gallego, A. Márquez, and I. Pascual, “Approximate solutions for the nonlinear pendulum equation using a rational harmonic representation,” Computers & Mathematics with Applications, vol. 64, no. 6, pp. 1602–1611, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. Beléndez, C. Pascual, M. L. Lvarez, D. I. Méndez, M. S. Yebra, and A. Hernndez, “Higher order analytical approximate solutions to the nonlinear pendulum by He's homotopy method,” Physica Scripta, vol. 79, no. 1, Article ID 015009, 2009. View at: Publisher Site  Google Scholar
 P. Brzeski, P. Perlikowski, S. Yanchuk, and T. Kapitaniak, “The dynamics of the pendulum suspended on the forced Duffing oscillator,” Journal of Sound and Vibration, vol. 331, no. 24, pp. 5347–5357, 2012. View at: Publisher Site  Google Scholar
 G. A. Baker, Essentials of Padé Approximants, Academic Express, London, UK, 1975.
Copyright
Copyright © 2013 H. VázquezLeal et al. 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.