Homotopy perturbation method with an auxiliary term for the optimal design of a tangent nonlinear packaging system

Abstract

The homotopy perturbation method with an auxiliary term is applied to obtain an approximate analytical solution for the tangent nonlinear packaging system. To further raise the accuracy of the solution, this work has expanded the possibility of the energy method for the homotopy perturbation method with an auxiliary term. And then, the amplitude and frequency of the tangent nonlinear packaging system can be calculated. The results show that the solution is in good consistent with the numerical solution solved by the Runge–Kutta method. It can be concluded that the proposed method is suitable for the optimal design of a tangent nonlinear packaging system.

Keywords

Tangent nonlinear packaging system homotopy perturbation method Runge–Kutta method energy method optimal design

Introduction

Shock and vibration are inevitable situations for the packaged products during transit from the manufacturing plants to the consumers, which are investigated by many scholars.^1–3 In fact it should be pointed out that the constitutive models of cushioning package materials are strong nonlinear and further effort is required to reveal the theoretical analysis of the packaging systems.

Recent developments of the analytical solutions of nonlinear problems have attracted much attention. Among the methods for analytical solutions, the perturbation method⁴ is one of the most well-known approaches based on the existence of small parameters. As we all know, most nonlinear equations have no small parameters at all, which has limited the application of the perturbation method. In order to eliminate the assumption of “small parameters”, considerable research efforts have been devoted to develop new methods, including the variational iteration method^5,6 and the homotopy perturbation method (HPM).^7–12 HPM is a coupled method of the homotopy technology and the classical perturbation method, which has been extensively used to solve many nonlinear problems.^13–17 Several excellent reviews describing the developments of the HPM are available. Li and He¹⁸ constructed a homotopy equation with a higher order coupled with the basic idea of the enhanced perturbation method. Noor¹⁹ presented a modification of the HPM with an auxiliary term. He²⁰ proposed an alternative approach by introducing an auxiliary term into the homotopy equation.

This paper investigates the applicability and validity of the HPM with an auxiliary term for the optimal design of a tangent nonlinear packaging system. Combined with the energy method, the amplitude and frequency of the tangent nonlinear packaging system are obtained, which are in accord with the results of Runge–Kutta method.

Homotopy equation with an auxiliary term

Considering a general nonlinear equation

L u + N u = 0

(1)

Where

L

and

N

are the linear operator and nonlinear operator, respectively.

The classic homotopy equation can be constructed as

\tilde{L} u + p (L u - \tilde{L} u + N u) = 0

(2)

where

\tilde{L}

is a linear operator and

\tilde{L} u = 0

can approximately describe the solution property. The embedding parameter

p

increases from 0 to 1, the trivial problem (

\tilde{L} u = 0

) will continuously transform to the original one (

L u + N u = 0

On the basis of He’s recent study,²⁰ the homotopy equation can be constructed with an auxiliary term as

\tilde{L} u + p (L u - \tilde{L} u + N u) + α p (1 - p) u = 0

(3)

Where

α

is an auxiliary parameter. It is obvious that when

α = 0

, equation (3) turns out to be the classical equation (2). And the parameter

α

can be determined in the following iteration procedure. When

p = 0

p = 1

, the auxiliary term will vanish completely. The solutions of equations (2) and (3) are expanded into a series in

p

u = u_{0} + p u_{1} + p^{2} u_{2} + \dots

(4)

0 = w^{2} + p a_{1} + p^{2} a_{2} + \dots

(5)

The tangent nonlinear packaging system

The tangent nonlinear packaging system commonly exists in packaging engineering, and the governing equations can be expressed as

m \ddot{x} + 2 k_{0} \frac{d_{b}}{π} tan \frac{π x}{d_{b}} = 0

(6)

x (0) = 0, \dot{x} (0) = \sqrt{2 g h}

(7)

where the coefficient

m

denotes the mass of the packaged product; the coefficient

x

denotes the displacement of the product; while

k_{0}

is the coupling stiffness of the cushioning material;

d_{b}

is the compression limit of the cushioning material, and

h

is the dropping height.

According to the Taylor series, equation (6) can be simplified as

m \ddot{x} + k_{0} x + r x^{3} + a x^{5} = 0

(8)

Introducing these parameters: $k = ε w_{0}^{2}$ , $ε = r / k_{0}$ , $b = σ w_{0}^{2}$ , $σ = a / k_{0}$ , $w_{0} = \sqrt{k_{0} / m}$ , equation (8) can be equivalently written in the following form

\ddot{x} + w_{0}^{2} x + k x^{3} + b x^{5} = 0

(9)

According to the homotopy equation with an auxiliary term, equation (9) can be calculated as

\ddot{x} + w^{2} x + p [(w_{0}^{2} - w^{2}) x + k x^{3} + b x^{5}] + α p (1 - p) x = 0

(10)

Substituting equation (4) into equation (10), the iteration equations can be obtained

{\ddot{x}}_{0} + w^{2} x_{0} = 0, x_{0} (0) = 0, {\dot{x}}_{0} (0) = V

(11)

{\ddot{x}}_{1} + w^{2} x_{1} + (w_{0}^{2} - w^{2} + α) x + k x^{3} + b x^{5} = 0

(12)

{\ddot{x}}_{2} + w^{2} x_{2} + (w_{0}^{2} - w^{2}) x + 3 k x_{1} x_{0}^{3} + 5 b x_{1} x_{0}^{5} + α (x_{1} - x_{0}) = 0

(13)

Solving equation (11), we have

x_{0} = A sin w t, A w = V

(14)

Where

A

and

w

are the displacement amplitude and the frequency, respectively.

Substituting equation (14) into equation (12) leads to

{\ddot{x}}_{1} + w^{2} x_{1} + I_{1} \sin w t + I_{2} \sin 3 w t + I_{3} \sin 5 w t = 0

(15)

where

I_{1} = (w_{0}^{2} - w^{2} + α) A + \frac{3}{4} k A^{3} + \frac{5}{8} b A^{5}

(16)

I_{2} = - \frac{1}{4} k A^{3} - \frac{5}{16} b A^{5}

(17)

I_{3} = \frac{1}{16} b A^{5}

(18)

In order to eliminate the secular term, it must be satisfied that $I_{1} = 0$ . So we have

w^{2} = (w_{0}^{2} + α) A + \frac{3}{4} k A^{3} + \frac{5}{8} b A^{5}

(19)

A special solution of equation (12) can be obtained as

x_{1} = \frac{I_{2}}{8 w^{2}} sin 3 w t + \frac{I_{3}}{24 w^{2}} sin 5 w t

(20)

The second-order iteration equation can be given by substituting equation (20) into equation (13)

{\ddot{x}}_{2} + w^{2} x_{2} + k_{1} \sin w t + k_{2} \sin 3 w t + k_{3} \sin 5 w t + k_{4} \sin 7 w t + k_{5} \sin 9 w t = 0

(21)

where

k_{1} = \frac{95 b^{2} A^{8}}{1536 w^{2}} + \frac{5 k A^{6} b}{64 w^{2}} + \frac{3 k^{2} A^{4}}{128 w^{2}} - α

(22)

k_{2} = \frac{k A^{2} (6 I_{2} - I_{3})}{32 w^{2}} + \frac{5 A^{4} b (9 I_{2} - 2 I_{3})}{192 w^{2}} + \frac{I_{2} (w_{0}^{2} - w^{2} + α)}{8 w^{2}}

(23)

k_{3} = - \frac{k A^{2} (3 I_{2} - 2 I_{3})}{32 w^{2}} - \frac{15 A^{4} b (2 I_{2} - I_{3})}{192 w^{2}} + \frac{I_{3} (w_{0}^{2} - w^{2} + α)}{24 w^{2}}

(24)

k_{4} = - \frac{k A^{2} I_{3}}{32 w^{2}} + \frac{5 A^{4} b (3 I_{2} - 4 I_{3})}{384 w^{2}}

(25)

k_{5} = \frac{5 A^{4} b I_{3}}{384 w^{2}}

(26)

No secular term in $x_{2}$ requires $k_{1} = 0$ , thus

\frac{95 b^{2} A^{8}}{1536 w^{2}} + \frac{5 k A^{6} b}{64 w^{2}} + \frac{3 k^{2} A^{4}}{128 w^{2}} - α = 0

(27)

Combining equations (19), (27) and the following parameters, $m = 10 kg$ , $h = 0.6 m$ , $k_{0} = 600 N / cm$ , $r = 72 N / {cm}^{3}$ , $a = 0.0001 N / {cm}^{5}$ , the displacement amplitude $A$ , the auxiliary term $α$ , and the frequency $w$ can be determined as 0.03186 m, 108.07, 107.64 s⁻¹, respectively.

Substituting equation (14) into equation (21) results in

x_{2} = \frac{k_{2}}{8 w^{2}} sin 3 w t + \frac{k_{3}}{24 w^{2}} sin 5 w t + \frac{k_{4}}{48 w^{2}} sin 7 w t + \frac{k_{5}}{80 w^{2}} sin 9 w t

(28)

The final solution $x$ can be obtained by $x = x_{0} + x_{1} + x_{2}$ . Thus

x = A sin w t + \frac{I_{2}}{8 w^{2}} sin 3 w t + \frac{I_{3}}{24 w^{2}} sin 5 w t + \frac{k_{2}}{8 w^{2}} sin 3 w t + \frac{k_{3}}{24 w^{2}} sin 5 w t + \frac{k_{4}}{48 w^{2}} sin 7 w t + \frac{k_{5}}{80 w^{2}} sin 9 w t

(29)

And the second-order acceleration iteration approximate expression can be written as

\ddot{x} = - A w^{2} sin w t - \frac{9 I_{2}}{8} sin 3 w t - \frac{25 I_{3}}{24} sin 5 w t - \frac{9 k_{2}}{8} sin 3 w t - \frac{25 k_{3}}{24} sin 5 w t - \frac{49 k_{4}}{48} sin 7 w t - \frac{81 k_{5}}{80} sin 9 w t

(30)

Thus the approximate solutions of the maximum displacement and acceleration can be expressed, respectively.

x_{m} = A - \frac{I_{2}}{8 w^{2}} + \frac{I_{3}}{24 w^{2}} - \frac{k_{2}}{8 w^{2}} + \frac{k_{3}}{24 w^{2}} - \frac{k_{4}}{48 w^{2}} + \frac{k_{5}}{80 w^{2}}

(31)

{\ddot{x}}_{m} = - A w^{2} + \frac{9 I_{2}}{8} - \frac{25 I_{3}}{24} + \frac{9 k_{2}}{8} - \frac{25 k_{3}}{24} + \frac{49 k_{4}}{48} - \frac{81 k_{5}}{80}

(32)

The energy method

In the analysis of dropping shock response of nonlinear packaging systems, the maximum displacement and acceleration of the products are usually obtained by the energy method. In order to raise the accuracy of the solution, our research extends the knowledge into the energy method, trying to correct the amplitude and the frequency of the tangent nonlinear packaging system. Thus

W h = \int_{0}^{x_{m}} (k_{0} x + r x^{3} + a x^{5}) d x

(33)

Where

W

and

h

are the weight of product and the drop height, respectively.

Solving equation (33), the maximum displacement $x_{m}$ can be determined as 3.401 cm.

According to equation (4), the maximum acceleration can be expressed as

{\ddot{x}}_{m} = w_{0}^{2} x + k x^{3} + b x^{5}

(34)

By solving equations (31) and (32) simultaneously, the new frequency $w$ and displacement amplitude $A$ can be determined as 111.2 s⁻¹ and 0.0333 m, respectively.

Results

The numerical solution solved by Runge–Kutta method is most equal to the exact solution, which can be plotted in Figure 1 (dotted line). The curve of the HPM with an auxiliary term is shown in Figure 1 (solid line), while the result of HPM with an auxiliary term combined with the energy method (HPM-EM) is also illustrated (dashed line). And a further comparison of the proposed method and Runge–Kutta method is performed, showing overall good agreement. To further verify the presented method, the maximum displacements and accelerations of the tangent nonlinear packaging system are listed in Table 1. On the basis of these results we concluded that good agreement is achieved between HPM-EM solution and the Runge–Kutta method solution.

Figure 1.

Comparison of solutions obtained by HPM and HPM-EM with that of Runge–Kutta method.

Table 1.

Results of the maximum displacements and accelerations for the tangent nonlinear packaging system.

	$x_{m}$ (cm)	Error (%)	${\ddot{x}}_{m}$ (g)	Error (%)
HPM	3.186	6.32	44.93	9.94
Runge–Kutta	3.401	-	49.89	-
HPM-EM	3.333	2.00	50.40	1.02

HPM-EM: homotopy perturbation method with energy method.

Conclusion

The analysis of dropping shock response for the packaged products is an unavoidable procedure in cushioning package design. The major objective of this paper is to develop the HPM with an auxiliary term for the optimal design for a tangent nonlinear packaging system. In order to raise the accuracy of HPM with an auxiliary term, our research extends the knowledge into the energy method. To gain more sight, the solutions of HPM-EM and Runge–Kutta method are compared. From the comparison results shown above, the relative error of HPM-EM can be decreased to 2%, which indicates that the solution of HPM-EM is in good agreement with the numerical solution. Our findings lead us to conclude that this study provides a scientific and effective method for the analysis of dropping shock response of a tangent nonlinear packaging system. Thus, the future looks bright for the presented method for different types of nonlinear packaging systems.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Natural Science Foundation of Jiangsu Province, China (Grant No. BK20151128), National first-class discipline program of Light Industry Technology and Engineering (LITE2018-29), the 111 Project (No. B18027), and Natural Science Foundation of Zhejiang Province (Grant number: LY16A020004).

References

Wang

ZW.

Dropping damage boundary curves for cubic and tangent package cushioning systems. Packag Technol Sci 2002; 15: 263–266.

Wang

Jiang

, et al. Dropping damage evaluation for a tangent nonlinear system with a critical component. Comput Math Appl 2011; 61: 1979–1982.

Piao

Saito

Kawaguchi

, et al. Effectiveness verification of the hybrid free-fall test. Packag Technol Sci 2017: 1–11.

Cole

JD.

Perturbation methods in applied mathematics. IEEE Trans Syst Man Cybern 1970; 23: 97.

XH.

Variational iteration method: new development and applications. Comput Math Appl 2007; 54: 881–894.

JH.

Variational iteration method – some recent results and new interpretations. J Comput Appl Math 2007; 207: 3–17.

JH.

Homotopy perturbation technique. Comput Methods Appl Mech Eng 1999; 178: 257–262.

JH.

A coupling method of a homotopy technique and a perturbation technique for non-linear problems. Int J Non-Linear Mech 2000; 35: 37–43.

JH.

Some Asymptotic methods for strongly nonlinear equations. Int J Mod Phys B 2006; 1141–1199.

10.

JH.

An elementary introduction to the homotopy perturbation method. Comput Math Appl 2009; 57: 410–412.

11.

JH.

Recent development of the homotopy perturbation method. Topological Meth Nonlinear Anal 2008; 31: 205–209.

12.

JH.

Homotopy perturbation method with two expanding parameters. Indian J Phys 2014; 88: 193–196.

13.

Seyyed

FMM.

Application of modified couple stress theory and homotopy perturbation method in investigation of electromechanical behaviors of carbon nanotubes. Adv Appl Math Mech 2017; 9: 23–42.

14.

Abou-Zeid

MY.

Homotopy perturbation method for couple stresses effect on MHD peristaltic flow of a non-Newtonian nanofluid. Microsyst Technol 2018; ( 5): 1–8.

15.

Saranya

Mohan

Kizek

, et al. Unprecedented homotopy perturbation method for solving nonlinear equations in the enzymatic reaction of glucose in a spherical matrix. Bioprocess Biosyst Eng 2018; 41: 1–14.

16.

Tripathi

Mishra

HK.

Application of homotopy perturbation method using Laplace transform intended for determining the temperature in the heterogeneous casting-mould system. Differ Equations Dyn Syst 2018; ( 1): 1–14.

17.

Fatahi-Vajari

Azimzadeh

Analysis of nonlinear axial vibration of single-walled carbon nanotubes using homotopy perturbation method. Indian J Phys 2018: 1–14.

18.

CH.

Homotopy perturbation method coupled with the enhanced perturbation method. J Low Freq Noise Vibr Active Control 2018. DOI: 10.1177/1461348418800554

19.

Noor

MA.

Some iterative methods for solving nonlinear equations using homotopy perturbation method. Int J Comput Math 2010; 87: 141–149.

20.

JH.

Homotopy perturbation method with an auxiliary term. Abstr Appl Anal 2012; ( 2): 183–194.