AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 424510 10.1155/2013/424510 424510 Research Article On the Convergence of the Homotopy Analysis Method for Inner-Resonance of Tangent Nonlinear Cushioning Packaging System with Critical Components Ghoreishi Mohammad 1 Ismail A. I. B. Md. 1 Rashid Abdur 2 Anderson Douglas 1 School of Mathematical Science Universiti Sains Malaysia 11800 Penang Malaysia usm.my 2 Department of Mathematics Gomal University 29050 Dera Ismail Khan Pakistan gu.edu.pk 2013 7 11 2013 2013 13 01 2013 12 08 2013 2013 Copyright © 2013 Mohammad Ghoreishi 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.

Homotopy analysis method (HAM) is applied to obtain the approximate solution of inner-resonance of tangent cushioning packaging system based on critical components. The solution is obtained in the form of infinite series with components which can be easily calculated. Using a convergence-control parameter, the HAM utilizes a simple method to adjust and control the convergence region of the infinite series solution. The obtained results show that the HAM is a very accurate technique to obtain the approximate solution.

1. Introduction

One of the most important subjects in packaging system is to investigate products damaged due to being dropped. Many researchers have investigated cushioning packaging system in this special field [1, 2]. The mechanical and electronic products are composed of large number of elements and, generally, damage at the so-called critical components. To prevent any damage, a critical component and a cushioning packaging are included in a package system . The following assumptions are made basically in the last decade by the researchers [1, 4].

The researchers considered that the packaging system is a spring-mass, single degree of freedom system.

The use of simple linear or nonlinear springs for cushioning packaging may not be appropriate.

Wang et al.  considered a linear model for this system, while the oscillation in the package system is inborn nonlinearity (see [4, 5]). Our goal of this paper is to obtain the approximate solution of inner-resonance of tangent nonlinear cushioning packaging system with critical components introduced in [1, 4] using HAM, which is one of the semilinear approximate analytical methods.

In the last two decades, many researchers have employed the approximate analytical methods such as adomian decomposition method (ADM), variational iteration method (VIM), homotopy perturbation method (HPM), and HAM to solve differential equations. These methods give the solution of the differential equations in the form of infinite series. One of the advantages of approximate analytical methods is that these methods do not produce rounding-off errors. Contrary to the implicit finite difference method (FDM), the approximate analytical methods do not require the numerical solution of the systems of differential equations.

We apply the HAM to construct the series solution for the inner-resonance of tangent nonlinear cushioning packaging system with critical components. An advantage of HAM over perturbation methods is that it is not dependent on small or large parameters. As it is well known, perturbation methods cannot be applied to all nonlinear equations because these methods are based on the existence of small or large parameters. Besides, nonperturbation methods are independent of small parameters. According to , both of the two techniques (perturbation and nonperturbative techniques) cannot provide a simple procedure to adjust or control the convergence region and rate of given approximate series. According to [7, 8], HAM also allows for fine-tuning of convergence region and rate of convergence by allowing an auxiliary parameter ħ to vary. To the best of our knowledge, this is the first attempt at solving the inner-resonance of tangent nonlinear cushioning packaging system with components approximately using the HAM.

In recent years, the HAM has been successfully applied for solving various nonlinear systems of equations in many branches of mathematics and sciences, such as strongly coupled reaction-diffusion system , fractional Lorenz system , coupled Schrodinger-KdV equation , Burgers and coupled Burgers equation , system of second-order BVPs , and HIV infection CD4+  T-cell . For more studies of HAM and its applications, the readers are also referred to see .

Our paper is organized as follows.

In Section 2, we introduce mathematical modelling of the inner-resonance of tangent nonlinear cushioning packaging system with critical components. In Section 3, we present a description of the HAM on system of equations, as expanded by previous researchers in particular [9, 12], applied to the inner-resonance of tangent nonlinear cushioning packaging system with critical components. We prove the convergence of homotopy series solution for the inner-resonance of tangent nonlinear cushioning packaging system with critical components also in this section. In Section 4, we have applied the HAM to obtain the approximate solution of inner-resonance of tangent nonlinear cushioning packaging system with critical components. Finally, in Section 5, we give the conclusion of this study.

2. Modeling and Equations

Generally, imagine that everything we know and have a relationship to, including things such as art, clothes, possessions, homes, gardens, trees and fields, mountains, lakes, oceans, continents, our friends, and loved ones, are instances where resonance can and does occur. Now, consider that everything we do and think is an attempt to seek out and to return to the experience of resonance, a return to the feeling of Belonging and feeling that things feel right. Even though we may not identify the motivation for making and forming certain relationships, the real attraction and value of any relationship is whether it is fulfilling and makes us feel good. We are searching for the feeling of resonance. The profound nature of these kinds of experiences is dependent on the nature and quality of the relationship we can have to something of an external nature, but we also have the potential to experience resonance within our own body and our inner being, and this I would call a “state of Inner Resonance.” When a practitioner consciously perceives resonance happening in a therapeutic relationship between themselves and patient, it is also a state of Inner Resonance. In this circumstance, the practitioners create the appropriate environment, conditions, and quality of presence so that they can receive the patients intention to be understood and received. In another way, the patient’s need to experience resonance (whether consciously or unconsciously) is met. Inner Resonance, then, is a state of receptive presence or conscious empathy between two or more people.

In Figure 1, the model of packaging system with critical component is shown to be considered as a nonlinear spring with stiffness coefficient  k2. It can also idealize the joining part between the mass of critical component  m1  and the main part of the product  m2  as a linear spring with stiffness coefficient  k1. According to Figures 1 and 2, the motion of this system can be written as [1, 4] (1)m1d2x(t)dt2+k1(x(t)-y(t))=0,m2d2y(t)dt2+2k2dbπtan(π2dby(t))-k1(x(t)-y(t))=0, with initial conditions (2)x(0)=0,x(0)=2gh,y(0)=0,y(0)=2gh.

The model of packaging system with critical component (m1m2) and  x1=x(t),  x2=y(t).

The  ħ-curves of  X(0.1)  and  Y(0.1)  obtained by the 5-order approximation of HAM.

Table 1 summarize the meanings of parameters and variables. To simplify (1), new variables are introduced as (3)X=xL,Y=yL,T=tT0, where (4)T0=m2k2,L=2dbπ. We define the frequency parameters of the critical component and main part of product as  w1=k1/m1  and  w2=k2/m2, respectively. The notations  λ1=w1/w2  and  λ2=m1/m2  are considered as parameter ratio and mass ratio, respectively. By considering all parameters defined, (1) can be equivalently written in the following system of nonlinear equations : (5)d2XdT2+w012X-w012Y=0,d2YdT2+w022Y+13Y3+215Y5+(1-w022)X=0, with initial conditions (6)X(0)=0,X(0)=T0L2gh,Y(0)=0,Y(0)=T0L2gh, where  X=X(T),  Y=Y(T)  and (7)w01=λ1,w02=1+λ12λ2.

List of variables and parameters (modified from [1, 4]).

Parameters and variables Illustration
x Displacement response of critical component
y Main body of the product
m 1 Mass of critical component
m 2 Main part of product
k 1 Stiffness coefficient
k 2 Stiffness coefficient
d b Compression limit of the cushioning pad
h Dropping height
g Gravity acceleration
2 g h Dropping shock velocity of the product
3. Homotopy Analysis Method (HAM)

To apply the HAM, the nonlinear system (5) is considered. We make initial gusses on  X(T)  and  Y(T)  such that they satisfy the initial conditions (6) that are defined as (8)X(0)=X0=T0L2ghT,Y(0)=Y0=T0L2ghT. The auxiliary linear operators X and Y are selected as (9)X=d2XdT2,Y=d2YdT2, satisfying the following properties: (10)X(c1,XT+c2,X)=0,Y(c1,YT+c2,Y)=0, where c1,X, c2,X, c1,Y, and c2,Y are integral constants. Define the homotopy maps (11)X(X^(T;q),Y^(T;q))0000=(1-q)X[X^(T;q)-X0(T)]000000-qħHX(t)NX[X^(T;q),Y^(T;q)],Y(X^(T;q),Y^(T;q))0000=(1-q)Y[Y^(T;q)-Y0(T)]000000-qħHY(t)NY[X^(T;q),Y^(T;q)], where q[0,1] is an embedding parameter, ħ is nonzero auxiliary parameter, HX and HY are auxiliary functions, and NX and NY are nonlinear operators that are defined as (12)NX[X^(T;q),Y^(T;q)]=2X^(T;q)T2+w012X^(T;q)-w012Y^(T;q),NY[X^(T;q),Y^(T;q)]=2Y^(T;q)T2+w022Y^(T;q)+(Y^(T;q))33+215(Y^(T;q))5+(1-w022)X^(T;q).   Clearly, when  q=0, we have the homotopy maps (13)X(X^(T;0),Y^(T;0))=X[X^(T;0)-X0(T)],Y(X^(T;0),Y^(T;0))=Y[Y^(T;0)-Y0(T)]. And when  q=1, we have (14)X(X^(T;1),Y^(T;1))0000=-ħHX(T)NX[X^(T;1),Y^(T;1)],Y(X^(T;1),Y^(T;1))0000=-ħHY(T)NY[X^(T;1),Y^(T;1)]. Thus, by requiring (15)X(X^(T;q),Y^(T;q))=Y(X^(T;q),Y^(T;q))=0, we can obtain (16)(1-q)X[X^(T;q)-X0(T)]0000=qħHX(T)NX[X^(T;q),Y^(T;q)],(1-q)Y[Y^(T;q)-Y0(T)]0000=qħHY(T)NY[X^(T;q),Y^(T;q)]. If  q=0  and  q=1, the homotopy equations are as follows: (17)X^(T;0)=X0,X^(T;1)=X(T),Y^(T;0)=Y0,Y^(T;1)=Y(T). As  q  varies from  0  to  1, the solution of the nonlinear system (5) will vary from the initial guesses  X0(T)  and  Y0(T)  to the exact solutions  X(T)  and  Y(T)  of the nonlinear system (5). Expanding  X^(T;q)  and  Y^(T;q)  as a Taylor series with respect to  q  yields (18)X^(T;q)=X0+m=1Xmqm,Y^(T;q)=Y0+m=1Ymqm, where (19)Xm=1m!mX^(T;q)qm|q=0,Ym=1m!mY^(T;q)qm|q=0. According to , the convergence of the series (18) strongly depends on the auxiliary parameter ħ. Note that if q=1, then (20)X^(T;1)=X=X0+m=1Xm,Y^(T;1)=Y=Y0+m=1Ym. According to definitions (18), the governing equations for the unknowns can be deduced from the zeroth-deformation equations (16). For further analysis, the vectors are defined as (21)X^n={X0,X1,,Xn},Y^n={Y0,Y1,,Yn}. Differentiating (16) m-times with respect to q, dividing by m!, and setting q=0 give the linear equations (22)X[Xm-χmXm-1]=ħHX(T)Rm,X(Xm-1,Ym-1),Y[Ym-χmYm-1]=ħHY(T)Rm,Y(Xm-1,Ym-1), with initial conditions (23)Xm(0)=0,Xm(0)=0,Ym(0)=0,Ym(0)=0, where (24)Rm,X(Xm-1,Ym-1)=Xm-1′′+w012Xm-1-w012Ym-1,(25)Rm,Y(Xm-1,Ym-1)=Ym-1′′+w022Ym-1+13i1=0m-1Yi1i2=0m-1-i1Yi2Ym-1-i1-i2+215i1=0m-1Yi1i2=0m-1-i1Yi2i3=0m-1-i1-i2Yi3×i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i4+(1-w022)Xm-1,(26)χm={0m11m>1. Using  HX(T)=HY(T)=1, the solution of the  m-order deformation equations (22) for  m1  becomes (27)Xm=χmXm-1+ħRm,X(Xm-1(τ),Ym-1(τ))dτdτ+c1,XT+c2,X,Ym=χmYm-1+ħRm,Y(Xm-1(τ),Ym-1(τ))dτdτ+c1,YT+c2,Y.   The coefficients c1,X, c2,X, c1,Y, and c2,Y are determined using initial conditions (23).

3.1. Convergence Theorem Theorem 1.

The series  X(T)=X0+m=1Xm  and  Y(T)=Y0+m=1Ym  converge where  Xm  and  Ym  are governed by (22) under definitions (24)–(26); X and Y must be the solutions of system of (5).

Proof.

If the series  m=0Xm  and  m=0Ym  are convergent, we can write (28)SX=m=0Xm,SY=m=0Ym. And it holds that (29)limmXm=limmYm=0. From (22) and using (9), we have (30)  ħm=1Rm,X(Xm-1,Ym-1)=m=1X[Xm-χmXm-1]=limmm=0nX[Xm-χmXm-1]=X[limmm=0n(Xm-χmXm-1)]=X[limmXn]=0,ħm=1Rm,Y(Xm-1,Ym-1)=m=1Y[Ym-χmYm-1]=limmm=0nY[Ym-χmYm-1]=Y[limmm=0n(Ym-χmYm-1)]=Y[limmYn]=0. Since  ħ0, then (31)m=1Rm,X(Xm-1,Ym-1)=0,(32)m=1Rm,Y(Xm-1,Ym-1)=0. Substituting (24) into (31) and simplifying it, we obtain (33)m=1Rm,X(Xm-1,Ym-1)0000=m=1(Xm-1′′+w012Xm-1-w012Ym-1)0000=m=1Xm-1′′+w012m=1Xm-1-w012m=1Ym-10000=d2dT2m=0Xm+w012m=0Xm-w012m=0Ym0000=X′′+w012X-w012Y=0. We repeat this process and substitute (25) into (32), and simplifying it, we obtain (34)m=1Rm,Y(Xm-1,Ym-1)  0000=m=1Ym-1′′  +w022m=1Ym-1  000000+m=1(1-w022)Xm-1000000+13m=1[i1=0m-1Yi1  i2=0m-1-i1Yi2Ym-1-i1-i2]  000000+215m=1[i1=0m-1Yi1i2=0m-1-i1Yi2i3=0m-1-i1-i2Yi3000000000000000×i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i4].

For the first three terms of (34), we can easily conclude that (35)m=1Ym-1′′+w022m=1Ym-1+m=1(1-w022)Xm-10000=Y′′(T)+w022Y(T)+(1-w022)X(T). For the fourth and fifth terms in (34), we have (36)  13m=1[i1=0m-1Yi1i2=0m-1-i1Yi2Ym-1-i1-i2]00=13m=1i1=0m-1Yi1i2=0m-1-i1Yi2Ym-1-i1-i200=13i1=0m=i1+1Yi1i2=0m-1-i1Yi2Ym-1-i1-i200=13i1=0j=0Yi1i2=0jYi2Yj-i200=13i1=0Yi1j=0i2=0jYi2Yj-i200  =13i1=0Yi1i2=0j=i2Yi2Yj-i200  =13i1=0Yi1i2=0Yi2j=i2Yj-i200  =13(m=0Ym)3=13Y3(T),215m=1[i1=0m-1Yi1i2=0m-1-i1Yi2i3=0m-1-i1-i2Yi3000000×i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i4]00=215m=1i1=0m-1Yi1i2=0m-1-i1Yi20000×i3=0m-1-i1-i2Yi3i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i400=215i1=0Yi1m=i1+1i2=0m-1-i1Yi20000×i3=0m-1-i1-i2Yi3i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i400=215i1=0Yi1j=0i2=0jYi2i3=0m-1-i1-i2Yi30000×i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i4  00=215i1=0Yi1i2=0Yi2m=1+i1+i2i3=0m-1-i1-i2Yi3  0000×i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i400=215i1=0Yi1i2=0Yi2l=0i3=0lYi30000×i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i400=215i1=0Yi1i2=0Yi2i3=0Yi30000×m=1+i1+i2+i3i4=0m-1-i1-i2-i3Yi4Ym-1-i1-i2-i3-i400=215i1=0Yi1i2=0Yi2i3=0Yi3k=0i4=0kYi4Yk-i400=215i1=0Yi1i2=0Yi2i3=0Yi3i4=0k=i4kYi4Yk-i400=215i1=0Yi1i2=0Yi2  i3=0Yi3  i4=0Yi4  k=i4kYk-i4  00=215Y5(T). Thus, (37)m=1Rm,Y(Xm-1,Ym-1)=Y′′+w022Y+13Y3+215Y5+(1-w022)X=0. From (2) and (23), it holds that (38)X0=m=0Xm(0)=X0(0)+m=1Xm(0)=0,X0=m=0Xm(0)=X0(0)+m=1Xm(0)=T0L2gh,Y0=m=0Ym(0)=Y0(0)+m=1Ym(0)=0,Y0=m=0Ym(0)=Y0(0)+m=1Ym(0)=T0L2gh. Thus,  X  and  Y  satisfy the system (5) and it must be the exact solution for (5) with the initial conditions (6).

4. Example

In this section, the HAM is applied to obtain the approximate solutions of the system (5) with the initial conditions (6). We have also used symbolic software Mathematica to solve the system of linear equations (22) with the initial conditions (23). Few components of the series solutions of (20) are given as follows: (39)X0=T0ghπT2db,Y0=T0ghπT2db,X1=0,Y1=ħπT3ghT062db(1+gπ2T2hm220db2k20+g2π4T4h2m22210db4k22),X2=ħ2πT5ghT0k1m21202dbk2m1(1+gπ2T2hm242db2k20+g2π4T4h2m22756db4k22),  Y2=ħπT3ghT062db(1+ħ+ħT22+ħT2k120k20  +għπ2T4hm2k1840db2k22+gπ2T2hm220db2k20  +11għπ2T4hm2840db2k2+), It is clear that the HAM series solutions (20) on depend on the convergence-control parameter  ħ  which provides a simple way to adjust and control the convergence of the series solutions. In fact, it is very important to ensure that the series (20) are convergent. To this end, we have plotted  ħ-curve of X(0.1) and Y(0.1) by fifth- and sixth-order approximation of the HAM in Figures 2 and 3, respectively, for values m1=0.01, m2=1, k1=0.1, k2=0.02, g=0.8, h=0.01, and  db=0.9. According to these ħ-curves, it is easy to discover the valid region of convergence-control parameter  ħ  which corresponds to the line segment nearly parallel to the horizontal axis. For clearer presentation, these valid regions have been listed in Table 2. Furthermore, these valid regions ensure us the convergence of the obtained series. Liao  has pointed out that when ħ=-1, the solution obtained by the HAM is the same as the series solution obtained using HPM.

The admissible value of ħ derived from Figures 2 and 3.

m 5 6
T ( t ) - 1.8 ħ - 0.7 - 1.3 ħ - 0.8
I ( t ) - 1.3 ħ - 0.75 - 1.25 ħ - 0.75

The  ħ-curves of  X(0.1)  and  Y(0.1)  obtained by the 6-order approximation of HAM.

Now, an error analysis is introduced to obtain the optimal value of convergence-control parameter ħ. Toward this end, we define φX(T;ħ) and φY(T;ħ) to be m-order approximation HAM solution as follows: (40)φX(T;ħ)=j=0m-1Xj,φY(T;ħ)=j=0m-1Yj. We substitute (40) into nonlinear system (5) and obtain the residual error functions ERX(X,Y,ħ1) and ERY(X,Y,ħ2) as follows: (41)ERX(X,Y;ħ1)=d2φX(T;ħ1)dT2+w012φX(T;ħ1)-w012φY(T;ħ1),(42)ERY(X,Y;ħ2)=d2φY(T;ħ2)dT2+w022φY(T;ħ2)+13(φY(T;ħ2))3+215(φY(T;ħ2))5+(1-w022)φX(T;ħ2).   Following , we define the square residual error for the m-order approximation to be (43)RX(ħ1)=01(ERX(X,Y,;ħ1))2dT,RY(ħ2)=01(ERY(X,Y,;ħ2))2dT.

We can obtain the values of ħ1 and ħ2 for which the RX(ħ1) and RY(ħ2) are minimum. The optimal values of convergence-control parameters ħ1 and ħ2 are determined by solving the system of equations as (44)dRX(ħ1*)dħ1=0,dRY(ħ2*)dħ2=0.

In , several methods have been introduced to find the optimal value of ħ. In Table 3, the minimum values of RX(ħ1) and RY(ħ2) have been given with optimal values of ħ1* and ħ2* for 4-, 5-, and 6-order approximations.

The minimum values of RX(ħ*) and RY(ħ*) for various orders of approximations.

m    X ( T ) Y ( T )
ħ * Minimum RX(ħ*) ħ * Minimum RY(ħ*)
4 −0.950213 1.60718 × 1 0 - 12 −0.950432 2.54142 × 1 0 - 8
5 −0.973770 3.80411 × 1 0 - 15 −0.971307 3.62018 × 1 0 - 11
6 −0.983485 3.84892 × 1 0 - 18 −0.993411 4.43599 × 1 0 - 14

In Table 4, the absolute errors ERX and ERY have been calculated for various T(0,1) when 5- and 6-order approximation HAM solutions are considered. From the table, it can be seen that the HAM provides us with the accurate approximate solution for the inner-resonance of tangent cushioning packaging system based on critical components (5).

The residual errors ERX and ERY for various T(0,1).

T    5-order 6-order
ER X ER Y ER X ER Y
0.1 7.65813 × 1 0 - 11 1.37420 × 1 0 - 9 3.18762 × 1 0 - 13 4.47395 × 1 0 - 11
0.2 3.06638 × 1 0 - 10 2.43550 × 1 0 - 8 1.76197 × 1 0 - 14 4.84994 × 1 0 - 11
0.3 1.58864 × 1 0 - 10 8.87158 × 1 0 - 8 3.49043 × 1 0 - 12 1.21465 × 1 0 - 9
0.4 2.05543 × 1 0 - 9 5.62311 × 1 0 - 7 9.40057 × 1 0 - 12 7.36418 × 1 0 - 10
0.5 2.20192 × 1 0 - 9 1.18024 × 1 0 - 6 7.17227 × 1 0 - 11 1.31996 × 1 0 - 8
0.6 1.04361 × 1 0 - 8 6.66260 × 1 0 - 7 5.64075 × 1 0 - 11 4.10080 × 1 0 - 8
0.7 5.15236 × 1 0 - 8 3.16150 × 1 0 - 7 6.52645 × 1 0 - 10 1.44031 × 1 0 - 8
0.8 1.17256 × 1 0 - 7 1.07335 × 1 0 - 5 2.96157 × 1 0 - 9 2.22440 × 1 0 - 7
0.9 1.17456 × 1 0 - 7 1.27292 × 1 0 - 5 4.99250 × 1 0 - 9 5.66406 × 1 0 - 7

The residual errors ERX and ERY have been plotted in Figures 4 and 5 for T(0,1) and various convergence-control parameters  ħ. By considering these Figures, it is to be noted that the solution obtained using HAM gives an analytical solution with high order of accuracy with few iterations.

The errors of residual equation (41) using the sixth-order approximate solution for various  ħ  and  T(0,1).

The errors of residual equation (42) using the sixth-order approximate solution for various  ħ  and  T(0,1).

5. Conclusion

The homotopy analysis method (HAM) is applied to obtain approximate solution of inner-resonance of tangent cushioning packaging system based on critical components. It is shown that the HAM solution contains the convergence-control parameter ħ, which provides a simple way to adjust and control the convergence region of the resulting infinite series. The convergence of HAM is also proved for inner-resonance of tangent cushioning packaging system based on critical components. The obtained results show that HAM is an accurate and effective technique for obtaining the approximate solution of inner-resonance of tangent cushioning packaging system based on critical components.

Wang J. Khan Y. Yang R. H. Lu L. X. Wang Z. W. Faraz N. A mathematical modelling of inner-resonance of tangent nonlinear cushioning packaging system with critical components Mathematical and Computer Modelling 2011 54 11-12 2573 2576 Wang J. Yang R. H. Li Z. B. Inner-resonance in a cushioning packaging system International Journal of Nonlinear Sciences and Numerical Simulation 2010 11 351 352 Wang J. Wang Z. Lu L. X. Zhu Y. Wang Y. G. Three-dimensional shock spectrum of critical component for nonlinear packaging system Shock and Vibration 2011 18 3 437 445 Wang J. Jiang J. H. Lua L. X. Wang Z. W. Dropping damage evaluation for a tangent nonlinear system with a critical component Computers and Mathematics with Applications 2011 61 1979 1982 Wang J. Wang Z. Damage boundary surface of a tangent nonlinear packaging system with critical components Journal of Vibration and Shock 2008 27 2 166 167 Rashidi M. M. pour S. A. M. Abbasbandy S. Analytic approximate solutions for heat transfer of a micropolar fluid through a porous medium with radiation Communications in Nonlinear Science and Numerical Simulation 2011 16 1874 1889 Alomari A. K. Noorani M. S. M. Nazar R. Li C. P. Homotopy analysis method for solving fractional Lorenz system Communications in Nonlinear Science and Numerical Simulation 2010 15 1864 1872 Bataineh A. S. Noorani M. S. M. Hashim I. Modified homotopy analysis method for solving systems of second-order BVPs Communications in Nonlinear Science and Numerical Simulation 2009 14 2 430 442 10.1016/j.cnsns.2007.09.012 MR2458820 ZBL1221.65196 Ghoreishi M. Ismail A. I. B. Md. Rashid A. Solution of a strongly coupled reaction-diffusion system by the homotopy analysis method Bulletin of the Belgian Mathematical Society 2011 18 3 471 481 MR2883141 ZBL1227.65097 Alomari A. K. Noorani M. S. M. Nazar R. Comparison between the homotopy analysis method and homotopy perturbation method to solve coupled Schrodinger-KdV equation Journal of Applied Mathematics and Computing 2009 31 1-2 1 12 10.1007/s12190-008-0187-4 MR2545707 ZBL1177.65152 Alomari A. K. Noorani M. S. M. Nazar R. The homotopy analysis method for the exact solutions of the K(2,2), Burgers and coupled Burgers equations Applied Mathematical Sciences 2008 2 40 1963 1977 MR2443893 Ghoreishi M. Ismail A. I. B. Md. Alomari A. K. Application of the homotopy analysis method for solving a model for HIV infection of CD4+ T-cells Mathematical and Computer Modelling 2011 54 11-12 3007 3015 10.1016/j.mcm.2011.07.029 MR2841844 ZBL1235.65095 Ghoreishi M. Ismail A. I. B. Md. Alomari A. K. Comparison between homotopy analysis method and optimal homotopy asymptotic method for nth-order integro-differential equation Mathematical Methods in the Applied Sciences 2011 34 15 1833 1842 10.1002/mma.1483 MR2838272 ZBL1226.65109 Ghoreishi M. Ismail A. I. B. Md. Alomari A. K. Bataineh A. S. The comparison between homotopy analysis method and optimal homotopy asymptotic method for nonlinear age-structured population models Communications in Nonlinear Science and Numerical Simulation 2012 17 3 1163 1177 10.1016/j.cnsns.2011.08.003 MR2843782 ZBL1239.92075 Ghoreishi M. Ismail A. I. B. Md. Rashid A. The solution of coupled modified KdV system by the homotopy analysis method TWMS Journal of Pure and Applied Mathematics 2012 3 1 122 134 MR2962059 ZBL1254.65114 Hosseini K. Daneshian B. Amanifard N. Ansari R. Homotopy analysis method for a fin with temperature dependent internal heat generation and thermal conductivity International Journal of Nonlinear Science 2012 14 2 201 210 MR2995565 ZBL06174868 Liao S. J. Beyond Perturbation: Introduction to the homotopy Analysis Method 2003 Boca Raton, Fla, USA Chapman and Hall/CRC Press Liao S. J. Comparison between the homotopy analysis method and homotopy perturbation method Applied Mathematics and Computation 2005 169 2 1186 1194 10.1016/j.amc.2004.10.058 MR2174713 ZBL1082.65534 Liao S. J. An optimal homotopy-analysis approach for strongly nonlinear differential equations Communications in Nonlinear Science and Numerical Simulation 2010 15 8 2003 2016 10.1016/j.cnsns.2009.09.002 MR2592613 ZBL1222.65088 Niu Z. Wang C. A one-step optimal homotopy analysis method for nonlinear differential equations Communications in Nonlinear Science and Numerical Simulation 2010 15 8 2026 2036 10.1016/j.cnsns.2009.08.014 MR2592615 ZBL1222.65093