# Variational Iteration Method of Dropping Shock Response for the Suspension Spring Packaging System

Authors: Shuang Song; An-Jun Chen
Journal: Shock and Vibration
 (2015)

---

## Abstract

In accordance with dropping shock dimensionless cubic nonlinear dynamic equation of suspension spring system, by variational iteration method, a first-order approximate solution of the system was obtained. The nondimensional peak of displacement, the nondimensional peak of acceleration, and the dropping shock extended period were compared with the results of the Runge-Kutta method, at which relative errors were less than 4%. The influence of suspension angle on peaks of response were discussed. It shows that the maximum response nondimensional acceleration decreases with decrease of the suspension angle under condition of the same nondimensional dropping shock velocity. Conditions for resonance were obtained by applying the variational iteration method, which should be avoided in the packaging design. The results provide reference for suspension spring system design.


---

## Body

## 1. Introduction

The damage evaluation concept proposed by Newton [1] is the foundation for present packaging design. In procedure of transportation and storage, dropping shock might lead to serious damage to product. The suspension spring system with eight springs as cushioning components performs geometric nonlinear; the absorber effect of suspension spring system is superior to linear system and is suitable for protecting high precision instrument with low fragility. Nonlinear suspension spring natural vibration was studied by Wu et al. [2] based on foundation displacement. For suspension spring system under rectangular pulse excitation, the three-dimensional shock spectrum and the three-dimensional damage boundary were obtained by Wang and Chen [3–5].For complexity of nonlinear systems, the numerical method is mainly used to analyze dropping shock characteristic [6–8]. The variational iteration method (VIM) [9–14] has been widely applied in solving kinds of nonlinear equations. For given initial approximate solution (may contain unknown parameters), a first-order approximate analytic solution is obtained, which is not limited by small parameter. Usually, the first-order approximate solution can meet action for engineering accuracy. A dynamic model with nonlinear cubic-quintic Duffing oscillators generated from the suspension spring packaging system was established, the VIM was applied to dynamic damage evaluation [14]. However, the dynamic model with nonlinear cubic-quintic Duffing oscillators is so complicated that the expression of the frequency parameter cannot be obtained clearly, the difficulty of calculation increases, and dynamic response factors analysis is inadequate.In this paper, for the dynamic model with cubic oscillators of the suspension spring system, the VIM is applied to get first-order approximate solution of system. The analytical expressions of important parameters with more clear physical meaning such as the frequency parameter, the nondimensional peak of displacement response, the nondimensional peak of acceleration response, the dropping shock extended period, and conditions for resonance of the system were obtained.


## 2. Modeling and Equations

The dynamic model of the suspension spring packaging system is shown in Figure1. 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.The approximate dropping shock dynamic equation [3–5] and dropping shock initial conditions of suspension spring system can be expressed as(1)md2xdt2+8ka0x+b0l02x3=0,(2)x0=0,(3)dx0dt=2gH,where(4)a0=sin2⁡ϕ0,(5)b0=121-6sin2⁡ϕ0+5sin4⁡ϕ0.Here the coefficient x denotes the product displacement, k denotes the coupling spring stiffness coefficient, m is the mass of product, H is the dropping height, ϕ0 denotes the suspension angle of the system, l0 denotes the original length of springs, and g is the acceleration of gravity.To simplify dropping shock dynamic equation, by introducing the new nondimensional parameters(6)y=xl0,(7)τ=tTthe frequency parameter of the system is defined as(8)ω=8kmand the period parameter of the system is defined as(9)T=1ω.Substituting all of parameters defined above into (1), the nondimensional form of the motion equation can be written in the following form:(10)d2ydτ2+a0y+b0y3=0.Equation (10) is the nondimensional shock dynamic equation. Initial conditions can be written as(11)y0=0,(12)dy0dτ=V=2gmH8kl02,where V denotes the nondimensional dropping velocity.By (10), the nondimensional dropping velocity and the suspension angle are related to the response of the system closely, which is characterized of geometric nonlinearity.


## 3. Variational Iteration Method

VIM [9] has been widely applied in solving different kinds of nonlinear equations and is especially effective in solving nonlinear vibration problems with approximation. For the nondimensional dropping shock dynamic equation (10) and the initial conditions (11) and (12), the initial solution can be taken as below:(13)y0τ=Vαsin⁡ατ,where α is the frequency parameter. By using the view of VIM, we can construct the following iteration formula:(14)yn+1τ=ynτ+∫0τ1a0sin⁡a0s-τ·d2ynsds2+a0yns+b0yn3sds.The first-order iteration approximate solution is obtained as(15)y1τ=3b0V34α3α2-a0sin⁡ατ-b0V34α39α2-a0·sin⁡3ατ+Va0-3b0V34a0α2α2-a0+3b0V34a0α29α2-a0sin⁡a0τ.In order to ensure that no secular terms appear in the next iteration, let the coefficient of sina0τ be equal to zero, and the frequency parameter can be expressed as(16)α=5a0+16a02+54b0V23.For the nondimensional shock dynamic equation (10), the first-order approximation solutions can be written as follows:(17)y1τ=3b0V34α3α2-a0sin⁡ατ-b0V34α39α2-a0sin⁡3ατ,(18)d2y1τdτ2=-3b0V34αα2-a0sin⁡ατ+9b0V34α9α2-a0sin⁡3ατ.During product packaging design procedure, the maximum of displacement response, and the maximum of acceleration response are comparatively important for product dropping shock damage evaluation, however for general nonlinear systems getting their analytic expressions are difficult. By the first-order approximation solutions (17) and (18), when ατ=π/2, the nondimensional peak of displacement and the nondimensional peak of acceleration can be written as(19)y1τm=3b0V34α3α2-a0+b0V34α39α2-a0,(20)d2y1τdτ2m=3b0V34αα2-a0+9b0V34α9α2-a0.The nondimensional dropping shock extended period can be obtained as(21)τ=πα.From (16), (19), (20) and (21), we know the nondimensional peak of displacement, the nondimensional peak of acceleration, and the nondimensional dropping shock extended period are related to V, a0, and b0. V is related to H, and parameters a0 and b0 are just related to ϕ0 (the design parameter of the system). For given suspension spring system, by using (19), (20) and (21), the relations of the shock response peak and the shock extended period can be obtained for the different dropping height. We can also analyze the suspension angle influence expediently. The first-order approximation solutions (17) and (18) have explicit physical significance, concise forms and can be analyzed easily.


## 4. Numerical Results

For the following amounts:l0=0.075m, k=2×105N/m, H=0.5m and ϕ0=60°. As shown in Figures 2 and 3, the nondimensional dropping shock response displacement and acceleration of the system are calculated and compared with the numerical integration solutions using the Runge-Kutta (ODE45) method, showing good agreement. From (19), (20) and (21), the nondimensional peak of displacement y1m=0.6908, the nondimensional peak of acceleration y¨1m=0.3912, and the nondimensional dropping shock extended period τ1=4.0109. With numerical integration solutions using ordinary differential equation-solver in MATLAB, the results are ym=0.6871, y¨m=0.4038 and τ=3.9666. Respectively, Compared with numerical integration solutions (ODE45), the VIM results relative errors are 0.54%, 3.12% and 1.12%. The first-order approximation solutions can meet requirement of product packaging design.Figure 2
Comparison of the nondimensional displacementy-time τ response of the system by the VIM with the one by the Runge-Kutta (ODE45) method when the suspension angle ϕ0=60°.Figure 3
Comparison of the nondimensional accelerationy′′-time τ response of the system by the VIM with the one by the Runge-Kutta (ODE45) method when the suspension angle ϕ0=60°.In the engineering application, analyzing the main factors affecting the maximum response of the system is very significant. According to (16), (19) and (20), we can find that the nondimensional dropping shock velocity and the suspension angle are the key factors affecting the maximum response. Respectively, we choose ϕ0 equal to 60°, 70° and 80°. Figure 4 indicates that the maximum response nondimensional displacement of the system increases with the increase of nondimensional dropping shock velocity, and increases with the decrease of the suspension angle under condition of the same nondimensional dropping shock velocity. Figure 5 indicates that the maximum response nondimensional acceleration of the system increases with the increase of the nondimensional dropping shock velocity, and decreases with the decrease of the suspension angle under condition of the same nondimensional dropping shock velocity.Figure 4
Maximum response nondimensional displacement of the system when the suspension angleϕ0=60°, 70°, 80°.Figure 5
Maximum response nondimensional acceleration of the system when the suspension angleϕ0=60°, 70°, 80°.


## 5. Resonance

In a cushioning packaging system, any small vibration might lead to serious damage due to inner-resonance. The inner-resonance is the key problem to optimal design [12]. By (17) and formula (16), the resonance can be expected when one of the following conditions is met:(22)α=5a0+16a02+54b0V23,(23)α=a0,(24)α=a03.To verify the inner-resonance phenomenon, for the nondimensional form of the motion equation (10) and the initial conditions (11) and (12), under the condition of harmonic excitation, the dynamic equation of the system can be established as follows:(25)d2ydτ2+a0y+b0y3=βsin⁡ωτ,where β is the nondimensional excitation amplitude and ω is the frequency parameter of the harmonic excitation.We chooseϕ0=65°, β=0.1, and V=0.1. When an external excitation frequency ω is equal to type (22), (23) and (24), the system acceleration response, respectively, as shown in Figures 6(a), 6(b), and 6(c). Figure 6 shows that the frequency of the exciting force to the system natural frequency, the system acceleration response amplitude regularly and alternately increased and reduced, known as the clap-frequency vibration phenomenon. Peak acceleration increase may cause product damage, which is dangerous for products. So inner-resonance should be avoided during cushioning packaging design procedure.Figure 6
The nondimensional acceleration response of the system under conditions for resonance.
(a)
ω
=
5
a
0
+
16
a
0
2
+
54
b
0
V
2
/
3
(b)
ω
=
a
0
(c)
ω
=
a
0
/
3


## 6. Conclusion

For the nondimensional dropping shock dynamic equation of the suspension spring system, first-order approximate analytical solutions of the nondimensional displacement and acceleration are obtained by He’s variational iteration method. Expressions of the nondimensional peak of displacement, the nondimensional peak of acceleration, and the nondimensional dropping shock extended period are obtained. Compared with the results of the Runge-Kutta method, relative errors are less than 4%. The influence of dropping height and suspension angle on the nondimensional displacement and the acceleration peak values was discussed. It shows that decreasing suspension angle appropriately the nondimensional acceleration peak value would decrease, while the nondimensional displacement peak value would increase. Conditions for resonance, which should be avoided in the product packaging design procedure, can be obtained. The results provide reference for suspension spring system design.



---

*Source: 408674-2015-07-26.xml*