Abstract and Applied Analysis

Volume 2014 (2014), Article ID 385404, 7 pages

http://dx.doi.org/10.1155/2014/385404

Research Article

Application of Variational Iteration Method for Dropping Damage Evaluation of the Suspension Spring Packaging System

Shuang Song,¹ Ning-ning Duan,¹ and An-jun Chen1,2

¹Department of Packaging Engineering, Jiangnan University, Wuxi 214122, China

²China National Light Industry Package Quality Supervising & Testing Center, Wuxi 214122, China

Received 8 April 2014; Revised 9 June 2014; Accepted 12 June 2014; Published 25 June 2014

Academic Editor: Patricia J. Y. Wong

Copyright © 2014 Shuang Song 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.

Abstract.

The dropping damage evaluation for packaging system is essential for safe transportation and storage. A dynamic model of nonlinear cubic-quintic Duffing oscillator for the suspension spring packaging system was proposed. Then, a first-order approximate solution was obtained by applying He’s variable iteration method. Based on the results, a damage evaluation equation was derived, which reveals the main controlling physical parameters for damage potential of drop to packaged products concretely. Finally, the dropping damage boundary curves and surfaces for the system were discussed. It was found that decreasing the suspension angle can improve the safe region of the system.

1. Introduction

Newton [1] proposed the concept of damage boundary for the first time in 1968, which established the foundation of present cushioning packaging dynamics. However, it can only be applied to linear packaging systems. For nonlinear packaging systems, there has been some undergoing work since then. A dropping damage evaluation for a tangent nonlinear system with a critical component was proposed by Wang et al. [2]. Wang et al. [3] proposed a three-dimensional shock spectrum for nonlinear packaging system with a critical component. They also suggested the damage boundary surfaces concept for damage evaluation of a tangent nonlinear packaging system with a critical component [4]. The fatigue damage of most packaged products is caused by dropping shock in the process of transportation. Therefore, Wang [5, 6] proposed the concept of dropping damage boundary curve with system parameter and the dimensionless dropping shock velocity as two basic evaluation quantities. These theories are all based on numerical analysis method. However, the influence of relevant parameters cannot be revealed to show their physical significance clearly. The variational iteration method (VIM) proposed by He et al. [7–9], which has been widely applied [10–19], can solve many kinds of nonlinear equations without small parameters limitation (first-order approximate analytical solutions can achieve high precision). Wang et al. [20–22] studied the dropping response of typical nonlinear packaging systems and obtained inner-resonance conditions.

The suspension spring system with eight springs as cushioning components performs geometric nonlinearity and is suitable for protecting high precision instrument with low fragility. Wu and Yang [23] studied the natural vibration characteristics of the system under the excitation of foundational displacement. They found that the shock absorption performance of the system with pendulum springs was better than the one with vertical springs. Assuming the system was excited by rectangular pulse, the three-dimensional shock spectrum and the damage boundary surface were obtained by Wang and Chen [24–26]. They achieved the conclusion that both increasing the pulse incentive amplitude and decreasing the suspension angle could expand the safe area, and the safe area would be expanded more obviously when the angle was less than 75°.

In this paper, by applying the VIM, we solve the nondimensional dynamic equation of the suspension spring system under the excitation of dropping shock to obtain a first-order approximate solution and obtain the nondimensional maximum acceleration expression. Then, a damage evaluation equation presenting the relationship between physical parameters and the damage boundary is suggested. Finally, the dropping damage boundary curves and surfaces of the system are discussed according to the damage evaluation equation.

2. Modelling and Equations

The dynamic model of the suspension spring packaging system is shown in Figure 1. A product is suspended in the middle of the container by 8 springs (four springs are on the upside, and the other four are on the downside).

Figure 1: Dynamic model of the suspension spring system.

Wang and Chen [24–26] proposed an approximate dynamic equation with a cubic oscillator of the suspension spring system. For acquiring a higher precision, we establish a more concrete dynamic model with nonlinear cubic-quintic Duffing oscillators of the system by using the Taylor series, which can be expressed as

where

Based on dropping shock, the initial conditions can be written as

Here are the coefficients: is the dropping height, denotes the product displacement, denotes the product mass, is the acceleration of gravity, denotes the coupling stiffness, denotes the original length of the springs, and denotes the suspension angle.

By introducing new nondimensional parameters, (1) can be equivalently written in the following dimensionless form:

where

The initial conditions can be written as

where denotes the dimensionless dropping velocity.

3. Variational Iteration Method

A nonlinear equation can be written as

where is a linear operator, is a nonlinear operator, and is a continuous function.

The VIM is proposed by He et al. [7–9] for the first time, and the correction functional for the nonlinear equation can be established as follows:

where is the Laplace multiplier which can be obtained by the variational theory and is the restriction on varying which is equal to 0 solving the Laplace multiplier.

Applying the VIM, construct the following correction functional of the system:

According to the principle of stationary, (9) can be turned into the following form:

The conditions of stationary can be written as

The Laplace multiplier is obtained as

The following iteration formula can be constructed as

For the nondimensional dynamic equation (4) and the initial conditions equation (6), we can take the initial solution below:

where is the frequency parameter and the amplitude . The first-order iteration approximate solution can be obtained as

Let the coefficient of be equal to zero so that there is no secular term appearing in the next iteration; namely,

As the dropping shock pulse is a half-sine pulse, the nondimensional dropping shock extended period can be obtained as

The first-order nondimensional displacement iteration approximate expression, namely, the first-order approximate solution, can be written as

The first-order nondimensional acceleration iteration approximate expression can be written as

Substitute into (18) and (19); the nondimensional maximum displacement can be written as

and the nondimensional maximum acceleration can be written as

For the following amounts, ,  N/m, and  m, we choose the dropping height , and the suspension angle . According to (16), we gain the frequency parameter . By applying the VIM the nondimensional maximum acceleration and the dropping shock extended period are obtained as and . Applying the Runge-Kutta (R-K) method, results are and . Comparing with the R-K method, the relative errors by using the VIM are and , which can meet the requirement of the packaging design.

The impact energy under dropping shock is related not only to the maximum acceleration, but also to the whole waveform. Hence, it is necessary to prove that the overall precision of the waveform can meet the requirement. The nondimensional acceleration response curves can be obtained just as shown in Figure 2. It indicates that the curve by applying the VIM can meet very well with the one by the R-K method.

Figure 2: Comparison of the nondimensional acceleration -time response of the system by the VIM with the one by the R-K method when the suspension angle .

4. Dropping Damage Evaluation

Wang [5, 6] proposed the concept of dropping damaging boundary, which can provide the theoretical foundation for the dropping damage evaluation of the suspension spring system.

The dropping shock acceleration of the system can be written as follows:

where the system parameter

Set as the product fragility. The relational expression about the nondimensional maximum acceleration , the product fragility , and the system parameter can be obtained just as follows:

By combining (21) and (24), the system parameter can be written as

Equation (25) includes more than one variable such as , , (related to and ), and . Therefore, according to (25), we can evaluate the dropping shock characteristics of the system.

Respectively, we choose the suspension angle , equal to 60°, 65°, 70°, and 75°, and select the product fragility , equal to 10 and 15. The system parameter and the dimensionless dropping shock velocity are selected as two basic evaluation quantities [5, 6]. According to (25), dropping damage boundary curves of the system can be obtained just as Figure 3 indicated. The safe area is under the dropping damage boundary curve, and the product is safe when the coordinate point enters the area.

(a)

(b)

Figure 3: Dropping damage boundary curves of the system when the suspension angle , , , and and the product fragility (a) and (b) . : the dimensionless dropping shock velocity; : the system parameter.

By selecting the suspension angle as the third evaluation quantity, set product fragility , , , and , respectively, and then dropping damage boundary surfaces can be obtained by applying the VIM, just like Figure 4 indicated. The safe area is under the dropping damage boundary surfaces, and the product is safe when the coordinate point enters the area.

(a)

(b)

(c)

(d)

Figure 4: Dropping damage boundary surfaces of the system when the product fragility (a) , (b) , (c) , and (d) . : the dimensionless dropping shock velocity; : the system parameter; : the suspension angle.

According to Figures 3 and 4, it is shown that decreasing the suspension angle , the security performance of the system will improve, and increasing the product fragility , the damage boundary curves and surface will move up obviously. So the characteristics of suspension geometry nonlinear () to protect products are superior to the linear system ().

5. Conclusion

A dynamic model with nonlinear cubic-quintic Duffing oscillators is proposed for thesuspension spring packaging system, and the first-order approximate solution of the equation is obtained by He’s variable iteration method. Based on the results, the damage evaluation equation for the packaging system is derived, revealing the main controlling parameters of the damage potential of dropping shock to packaged products. Finally, the damage boundary curves and surfaces for the system are discussed. It is found that decreasing the suspension angle can help to protect the packaged product.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

R. E. Newton, Fragility Assessment Theory and Practice, Monterey Research Laboratory, Monterey, Calif, USA, 1968.

J. Wang, J. Jiang, L. Lu, and Z. Wang, “Dropping damage evaluation for a tangent nonlinear system with a critical component,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 1979–1982, 2011.

J. Wang, Z.-W. Wang, L.-X. Lu, Y. Zhu, and Y.-G. Wang, “Three-dimensional shock spectrum of critical component for nonlinear packaging system,” Shock and Vibration, vol. 18, no. 3, pp. 437–445, 2011.

J. Wang, F. Duan, J. Jiang, L. Lu, and Z. Wang, “Dropping damage evaluation for a hyperbolic tangent cushioning system with a critical component,” Journal of Vibration and Control, vol. 18, no. 10, pp. 1417–1421, 2012.

Z. Wang, “Dropping damage boundary curves for cubic and tangent package cushioning systems,” Packaging Technology and Science, vol. 15, no. 5, pp. 263–266, 2002.

Z. Wang, “On evaluation of product dropping damage,” Packaging Technology and Science, vol. 15, no. 3, pp. 115–120, 2002.

J. H. He and X. H. Wu, “Variational iteration method: new development and applications,” Computer & Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.

J. He, “Variational iteration method—some recent results and new interpretations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 3–17, 2007.

J. H. He, G. C. Wu, and F. Austin, “The variational iteration method which should be followed,” Nonlinear Science Letters A, vol. 1, no. 1, pp. 1–30, 2010.

E. Yusufoglu, “The variational iteration method for studying the Klein-Gordon equation,” Applied Mathematics Letters, vol. 21, no. 7, pp. 669–674, 2008.

Z. Odibat and S. Momani, “The variational iteration method: an efficient scheme for handling fractional partial differential equations in fluid mechanics,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2199–2208, 2009.

D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation 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.

M. Rafei, H. Daniali, and D. D. Ganji, “Variational interation method for solving the epidemic model and the prey and predator problem,” Applied Mathematics and Computation, vol. 186, no. 2, pp. 1701–1709, 2007.

H. Shahmohamadi and M. M. Rashidi, “A novel solution for the glauert-jet problem by variational iteration method-padé approximant,” Mathematical Problems in Engineering, vol. 2010, Article ID 501476, 7 pages, 2010.

I. Kucuk, “Active optimal control of the KdV equation using the variational iteration method,” Mathematical Problems in Engineering, vol. 2010, Article ID 929103, 10 pages, 2010.

Y. Khan, H. Vázquez-Leal, and L. Hernandez-Martínez, “Removal of noise oscillation term appearing in the nonlinear equation solution,” Journal of Applied Mathematics, vol. 2012, Article ID 387365, 9 pages, 2012.

Y. Khan, H. Vázquez-Leal, and N. Faraz, “An efficient new iterative method for oscillator differential equation,” Scientia Iranica, vol. 19, no. 6, pp. 1473–1477, 2012.

Y. Khan, M. Madani, A. Yildirim, M. A. Abdou, and N. Faraz, “A new approach to Van der Pol's oscillator problem,” Zeitschrift fur Naturforschung A: Journal of Physical Sciences, vol. 66, no. 10-11, pp. 620–624, 2011.

Y. Khan, N. Faraz, A. Yildirim, and Q. B. Wu, “Fractional variational iteration method for fractional initial-boundary value problems arising in the application of nonlinear science,” Computers & Mathematics with Applications, vol. 62, no. 5, pp. 2273–2278, 2011.

J. Wang, Y. Khan, R. Yang, L. Lu, Z. Wang, and N. Faraz, “A mathematical modelling of inner-resonance of tangent nonlinear cushioning packaging system with critical components,” Mathematical and Computer Modelling, vol. 54, no. 11-12, pp. 2573–2576, 2011.

J. Wang, Y. Khan, L. X. Lu, and Z. W. Wang, “Inner resonance of a coupled hyperbolic tangent nonlinear oscillator arising in a packaging system,” Applied Mathematics and Computation, vol. 218, no. 15, pp. 7876–7879, 2012.

J. Wang, Z. W. Wang, F. Duan, L. Lu, D. Gao, and A. Chen, “Dropping shock response of corrugated paperboard cushioning packaging system,” Journal of Vibration and Control, vol. 19, no. 3, pp. 336–340, 2013.

X. Wu and L. J. Yang, “Natural vibration characteristics of suspension spring geometric nonlinear damping system,” Journal of Vibration and Shock, vol. 27, no. 11, pp. 71–72, 2008.

L. Wang and A. J. Chen, “Shock characteristics of a suspension spring system under action of a rectangular pulse,” Journal of Vibration and Shock, vol. 31, no. 11, pp. 142–144, 2012.

L. Wang and A. Chen, “The damage boundary curve of the suspension packaging system under rectangular pulse,” Applied Mechanics and Materials, vol. 105–107, pp. 70–73, 2012.

L. Wang and A. J. Chen, “The shock response spectra of the suspension packaging system under rectangular pulse,” Journal of Applied Packaging Research, vol. 5, no. 4, pp. 237–246, 2011.