Abstract and Applied Analysis

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

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

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

^{1}Department of Packaging Engineering, Jiangnan University, Wuxi 214122, China^{2}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).

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.

#### 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.

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.

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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - Z. Wang, “On evaluation of product dropping damage,”
*Packaging Technology and Science*, vol. 15, no. 3, pp. 115–120, 2002. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Google Scholar - E. Yusufoglu, “The variational iteration method for studying the Klein-Gordon equation,”
*Applied Mathematics Letters*, vol. 21, no. 7, pp. 669–674, 2008. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Google Scholar - 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. View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Google Scholar