Published: 15 May 2017

Study of double-potential-well leaf spring system’s chaotic vibration

Shuyong Liu1
Jiang Jian2
Pan Su3
Jiechang Wu4
Ye Liu5
Yuan Fang6
1, 3, 4, 5, 6College of Power Engineering, Naval University of Engineering, Wuhan, China
2Institute of Acoustic, Chinese Academy of Sciences, Beijing, China
Corresponding Author:
Pan Su
Views 218
Reads 91
Downloads 1577

Abstract

Chaotic vibration has become increasingly popular in the study of acoustic and vibration engineering. Many engineering designs have taken advantage of the special characteristics of chaos, and deliberately introduced it into the system to improve efficiency. As an important component, leaf springs have long been used in the suspension system of wheeled vehicles. Recent development is considering chaotic vibration in the design of leaf springs to improve the system’s reliability. However, little experimental research has been carried out to investigate the chaos characteristics of leaf springs. Meanwhile, a preliminary study showed that some of the conventional signal processing methods may not be able to successfully identify the chaos features from a leaf spring test rig due to the complexity of the practical signal. Therefore, in this paper, a leaf spring system’s chaotic vibration and relevant signal processing strategy were investigated in theory and experiment. Firstly, the relationship between the amplitude and frequency of the double potential well system is derived with averaging method. The stability is analyzed on the Vander pol plane and the global bifurcation diagram and Lyapunov exponent spectrum are applied to determine the chaotic regime accurately. Numerical simulation was conducted using a finite element method to give an idea of the leaf spring’s natural frequencies where chaotic vibration can be potentially generated. The experimental rig was then designed based on double potential well theory to generate stable and repeatable chaotic vibration, and an experimental study was carried out to investigate the system’s response characteristics under different excitation strengths and frequencies. An improved signal processing method, Wavelet-SG-EEMD (Wavelet, Savitzky-Golay (SG) and Ensemble Empirical Mode Decomposition (EEMD)), was used to reduce noise and beneficial to identify chaotic features of the vibration signal generated by the system. The nonlinear vibration response features of the system were carefully analyzed. Sub-harmonic phenomena, periodic modes and chaotic behavior were discovered during the experiment.

1. Introduction

Chaos has become increasingly popular in resolving complicated problems in engineering fields [1]. On one hand, a lot of projects have been conducted to identify and understand chaotic behaviors in a system and try to reduce the influence from chaos. Gao et al. developed a self-adaptive tracking control system to constrain large amplitude chaotic motion occurred in a motor suspension system [2] and Han et al developed a sliding mode controlling method to reduce Lorenz chaotic vibration [3]. On the other hand, some engineering solutions/designs were trying to take advantage of a system’s chaotic behavior and deliberately introduce chaos to improve the system’s efficiency. For example, Long et al. invented a chaotic vibration roller which significantly increased the device’s working efficiency [4]. DMITRIEV et al study the speech and music signal transmission by chaos [5]. Govindan R. B. et al. applied chaos theory to diagnose the patient ECG [6].

As an important component, leaf springs have long been used in the suspension system of wheeled vehicles. Recent development is considering nonlinear vibration in the design of leaf springs to improve the system’s reliability. The nonlinear large deformation, contact and friction are usually applied to characterize the nonlinear behavior of leaf springs. OMAR et al. [7] researched a nonlinear finite-element procedure to accurately model the deformations and vibrations of the leaf springs, in which the effect of the distributed inertia and elasticity are considered. Qing Li et al studied the multileaf spring in vehicles and develop a contact finite element algorithm. Then, a piecewise contact stress pattern is approximated to the original contact state between two layered beams [8]. Yong-Jin Yum researched the frictional characteristics of automotive leaf spring. The load-strain results of the compression test are obtained, and the finite element analysis was used to calculate the friction force [9]. Hiroyuki Sugiyama et al present a nonlinear elastic model of leaf springs to simulate the multibody vehicle systems. The nonlinear dynamic coupling between the finite rotations and the leaf deformation are analyzed, and the leaf spring geometry and deformations are modelled with the distributed inertia and stiffness [10].

Because the leaf spring has the same structure dynamics as the beam, it is beneficial to design and analyze spring base on the beam theory, and the analogous complex behavior of beam will occur in the leaf spring. Rahman et al designed a parabolic leaf spring, and used both the small and large deflection theories to calculate the stress and the deflection of the beam, in which the geometric nonlinearity is analyzed [11]. Dwivedy et al. investigate the nonlinear response of a base-excited slender beam carrying an attached mass, and the fixed-point, periodic, quasiperiodic and chaotic vibration are observed [12]. Nayfeh et al. studied the nonlinear dynamics of the fixed-free in extensional beam under the principal parametric excitation. The cubic nonlinearities due to the curvature and inertial are contained in the equations of motion. For some range of parameters, Hopf bifurcation exist, and response consists of amplitude and phase modulated or chaotic motions [13]. Li Chen proposed a moment integral treatment approach and applied it to the problems of complex and varying beam properties [14]. Nallathambi et al. developed an algorithm to analyze the large deflection of curved prismatic cantilever beams with uniform curvature subjected to a follower load at tip. The novel method reducing two-point boundary value problem to an initial value problem, and a single parameter shooting integrating strategy of the cantilever beam is develop [15].

Since the leaf spring undergoes large deflections when it is subjected to external excitation, the solutions of those highly nonlinear problems become very complex. Elliptic integral formulation [16], Galerkin procedure, R-K method with shooting technique [17, 18], finite difference method [19] and finite element method(FEM) [20] are the frequently used. Nayfeh et al. [21] applied the combination of Galerkin procedure and the method of multiple scales to construct a first-order uniform expansion for the inextensional beam. The results show that the nonlinear inertial terms produce a softening effect and play a significant role in the planar response of high frequency modes. Besides, the double potential well Duffing theory is used to research the nonlinear dynamics analytically. S. Li et al. used a self-excited term of Rayleigh type and Duffing double well potential to model vehicle suspension system due to electro- or magneto-rheological fluid damping where it is causing a hysteretic effect [22]. Thompson studied the manner of metastable mechanical oscillator escaping from cubic potential well. The basin of attraction is obtained, and he Melnikov method is applied to predict the homoclinic tangle [23]. Teueba et al. studied the nonlinearly damped double-well Duffing oscillator and analytically estimates how certain nonlinear damping terms affects the dynamics of the nonlinear oscillator. The condition of Melnikov-equivalence between the nonlinear damping system and the linear damping one is obtained [24]. Grzegorz Litak et al. introduced a harmonic excitation term and damping as perturbations, the critical Melnikov amplitude of the road surface profile is found, above which the system can vibrate chaotically [25].

Experimental studies for these systems are particularly important. However, few experimental studies have been carried out to investigate the chaos characteristics of leaf springs carefully. Liang Shan et al investigate the chaotic vibration of a vehicle model over road excitation [26]. Moon and Holmes investigated chaotic vibration of the beam in the magnetic field [27]. Fossas et al. studied the chaos in buck converter [28]. One common problem challenging all of these experimental studies is how to keep the system under a certain chaotic regime and produce stable and repeatable chaotic vibration. Besides, the signal processing methods are required to extract the chaotic characteristic effectively.

Therefore, in this paper, an experimental study was carried out to investigate the chaotic vibration of leaf springs. The natural frequencies where the chaotic vibration can be potentially excited were obtained from static analysis using a finite element based numerical model. The dynamic characteristics of the leaf spring were investigated experimentally on a test rig which was designed based on the double potential well theory [27]. A newly developed signal processing method, Wavelet-SG-EEMD )Wavelet, Savitzky-Golay (SG) and Ensemble Empirical Mode Decomposition(EEMD)) [29-32], were used to reduce noise and help to identify chaos features of the vibration signal generated by the system. A phase space reconstruction method, based on the first minimum mutual information (FMMI) theory, was used to obtain the chaotic attractors. An improved fast space grid algorithm [33] was used to obtain the maximum Lyapunov Exponent (LE) of the signal. A modified Poincaré method [34] was used to obtain the experimental Poincaré section.

In Section 2, the nonlinear amplitude and frequency characteristic curve of the double potential well system is analyzed, and the bifurcation is researched numerically. In Section 3, natural frequencies and relevant modes of the leaf spring used in the experiment are obtained by FE methods. In Section 4 the design of the leaf spring test rig is introduced. Details of Wavelet-SG-EEMD denoising method are presented in Section 5. Results from the experimental signals are analyzed and discussed in Section 6. Conclusions are given in Section 7.

2. Mechanism of the double-potential-well

2.1. Mathematic model

The differential equation of motion for the leaf spring system with magnets is simplified to double potential well system and it is modeled as the Duffing equation [27]:

1
z¨+ξz˙+γz+βz3=fcosω~τ,

where ξ is the dimensionless damping factor, γ is dimensionless linear stiffness, β is dimensionless nonlinear stiffness, f is the dimensionless force amplitude, ω~ is dimensionless excitation frequency, z¨=d2z/dτ2, z˙=dz/dτ. For the Eq. (1), the nonlinear restoring force of the spring is:

2
Fz=-γz-βz3.

The potential function V(z) is defined as:

3
V(z)=γz2/2+βz4/4.

In Eq. (3), when the parameter γ is –1, and β is 1 respectively, the double potential well exists objectively, as shown in Fig. 1. Therefore, the mass point can move between the double potential well M and N or in one of the potential well. If the mass point oscillates between M and N periodically, as shown in Fig. 1(a), the periodic vibration with large magnitude occurs in the leaf spring system. If the mass point oscillates in M or N only, the periodic vibration with small magnitude occurs, as shown in Figs. 1(b), (c). The attractive and complicated case is that the mass point oscillates between M and N irregularly, corresponding to the chaotic motion, as shown in Fig. 1(d).

Fig. 1Double potential well diagram and different oscillation of the ball

Double potential well diagram and different oscillation of the ball

a)

Double potential well diagram and different oscillation of the ball

b)

Double potential well diagram and different oscillation of the ball

c)

Double potential well diagram and different oscillation of the ball

d)

2.2. Amplitude-frequency characteristic curve

In order to analyze the amplitude frequency curve characteristic of the system, the analytical relationship between the amplitude and frequency parameter is derived through the average method [35]. Firstly, the excitation frequency is set to ω~=ω0+εp~ in Eq. (1), where p~ is detuning factor of frequency, and perturbation parameter ε is introduced. Then Eq. (1) is rewritten as:

4
z¨+ξz˙+ω02z+βz3=fcosω0+εp~τ.

Let ξ=εδ~, β=εα~, f=εA~, and substitute into Eq. (4), we get:

5
z¨+εδ~z˙+ω02z+εα~z3=εA~cosω0+εp~τ.

To simplify analysis, substitute new time scale variables τ=ωt/(ω0+εp~),d/dτ=d/[ω/(ω0+εp~)]dt into Eq. (5), and the equivalent equation is obtained:

6
z¨+εδz˙+ω2+εpz+εαz3=εAcosωt,

where:

δ=δωω0+εp, α=αωω0+εp2, A=Aωω0+εp2, ω2+εp=ω2 ω0ω0+εp2.

In Eq. (6), when ε is 0 the equation is z¨+ω2z=0, it is the linear one. When εis small values, the system exhibits weak nonlinear characteristic. Eq. (6) is rewritten as the following state equation:

7
z˙1=z2,z˙2=-εδz2-ω2+εpz1-εαz13+εAcosωt.

To analyze the behavior on the Vander Pol plane, the coordinate transformation is conducted. And the new coordinates are:

8
z1=utcosωt+vtsinωt,z2=-utsinωt+vtcosωt.

The derivation of Eq. (8) is:

9
z˙1=-ωutsinωt+ωvtcosωt+cosωtu˙t+v˙tsinωt,z˙2=-ωutcosωt-ωvtsinωt-u˙tsinωt+v˙tcosωt.

Substitute Eqs. (8) and (9) into Eq. (7), and get:

10a
u(t)sinωt-ωu(t)sinωt-v(t)cosωt+ωv(t)cosωt+u˙(t)cosωt+v˙(t)sinωt=0,
10b
-Aεcosωt+pεu(t)cosωt-ωu(t)cos(ωt)+ω2u(t)cos(ωt)-δεu(t)sinωt +αεu3tcos3ωt+δεv(t)cosωt+pεv(t)sinωt-ωsinωtv(t)+ω2sinωtv(t) +3αεu2tvtcos2ωtsinωt+3αεu(t)v2(t)cosωtsin2ωt +αεv3(t)sin3ωt-u˙(t)sinωt+v˙(t)cosωt=0.

The u˙(t) and v˙(t) are found with Eqs.10 (a), (b):

11a
u˙(t)=αεu3(t)u(t)cos3ωtsinωt+(3/4)αεu2(t)v(t)sin22ωt -u(t)sinωt[(1-pε-ω2)cosωt+δεsinωt-3αεv2(t)cosωtsin2ωt] -(1/2)pεsin2ωt+(1/2)[1+pε-2ω+ω2 -(-1+pε+ω2)cos2ωt+δεsin2ωt]v(t)+αεv3(t)sin4ωt,
11b
v˙(t)=-αεu3(t)cos4ωt-3αεu2(t)v(t)cos3ωtsinωt -(1/4)u(t){2[1+pε+ω2-2ω+(-1+pε+ω2)cos2ωt-δεsin2ωt] +3αεsin22ωtv2(t)} +cosωt{Aεcosωt-[δεcosωt+(-1+pε+ω2)sinωt]]v(t)-αεsin3ωtv3(t)},

Let ωt=θ, and Eqs. 11(a), (b) is translated into autonomous system:

12a
u˙(t)=αεu3(t)u(t)cos3θsinθ+(3/4)αεu2(t)v(t)sin22θ -u(t)sinθ[(1-pε-ω2)cosθ+δεsinθ-3αεv2(t)cosθsin2θ]+(1/2)pεsin2θ +(1/2)[1+pε-2ω+ω2-(-1+pε+ω2)cos2θ+δεsin2θ]vt+αεv3tsin4θ,
12b
v˙(t)=-αεu3(t)cos4θ-3αεu2(t)v(t)cos3θsinθ-(1/4)u(t){2[1+pε+ω2-2ω -(1-pε-ω2)cos2θ-δεsin2θ]+3αεv2(t)sin22θ} +cosθ{Aεcosθ-[δεcosθ+(-1+pε+ω2)sinθ]]v(t)-αεsin3θv3(t)}.

Assume that the amplitude u(t) and v(t) of the response are slow variation parameters, integrate the right-hand side of the Eq. (12) in the range of [0, 2π], and then evaluate the average. The following equations are obtained:

13a
ut=-12δεut+12vt+12pεvt-ωvt+38αεu2(t)v(t)
+38αεv3t+12ω2v(t),
13b
vt=12Aε-12ut-12pεut+ωut-12ω2ut-38αεu3(t)
-12δεvt-38αεutv2t.

The, u(t)v(t) are expressed in the form of polar coordinates u(t)=rcosϕ, v(t)=rsinϕ,and substitute it into Eq. (13), we get:

14a
-r˙cosϕ+rϕ˙sinϕ-12rδεcosϕ+12rsinϕ+12pεrsinϕ,-ωrsinϕ+12ω2rsinϕ+38αεr3cos2ϕsinϕ+38αεr3sin3ϕ=0,
14b
-r˙sinϕ-rϕ˙cosϕ+12Aε-12rcosϕ-12pεrcosϕ-12ω2rcosϕ,-38αεr3cos3ϕ+ωrcosϕ-12δεrsinϕ-38αεr3cosϕsin2ϕ=0.

Therefore, the r˙ and ϕ˙ are found:

15a
r˙=12-rδε+Aεsinϕ,
15b
ϕ˙=-r[4pε+3r2αε+4(ω-1)2]+4Aεcosϕ8r.

The equilibrium of Eq. (15) on the Vander Pol plane is satisfied with:

16a
12(-rδε+Aεsinϕ)=0,
16b
-r[4pε+3r2αε+4(ω-1)2]+4Aεcosϕ8r=0.

And thus, the following equations are derived:

17a
sinϕ=rδA,
17b
cosϕ=r(4+4pε+3r2αε-8ω+4ω2)4Aε.

Remove the variable ϕ and get:

18
1Aε{16A2ε2-9r6α2ε2-24r4αε[pε+(ω2-1)2]
-16r2[p2ε2+δ2ε2+2pε(ω-1)2+(ω-1)4]}=0.

Fig. 2Nonlinear vibration system amplitude-frequency curve

Nonlinear vibration system amplitude-frequency curve

Since the system natural frequency is changed in essence when parameter p is detuned in Eq. (6), the nonlinear characteristic of the system can be revealed through the relationship between the r and p. In Eq. (18), when the parameters are ε=1, α=1, ω=1, δ=0.5, A=2, 4, 6. The amplitude-frequency curve is obtained as shown in Fig. 2. In this figure, the tongue structure is observed obviously, and the resonant frequency shift phenomenon is exhibited. The advantage for the nonlinear vibration isolation system is to avoid the large vibration amplitude at resonance regime.

2.3. Global dynamic behavior characteristic of the system

The global dynamic behavior characteristic of the system can be analyzed on the Vander Pol plane, and reveals its long-term trend. The advantage of this method is to avoid the infinite time numerical analysis to judge whether the system terminal behavior is stable or instable. Specially, when the parameters are close to the tongue area, the phase space flow moves at slow speed, but the trajectory will be increased in secular when time approaches infinity. Secondly, the system behavior can be controlled through the parameters based on the Van der Pol plane. Three cases are disused in detail as follows:

(1) The parameter γ is varied. When the parameters ξ=0.1, γ=-1, f=10, ω~=3, system approaches to a stable periodic state, as shown in Fig. 3(a). When γ is change to –6, the basic property of the dynamics is invariant except for the trajectories on the left-hand side of the plane become denser, as shown in Fig. 3(b).

Fig. 3The convergent dynamical characteristic of the system when parameter γ is varied

The convergent dynamical characteristic of the system when parameter γ is varied

a)

The convergent dynamical characteristic of the system when parameter γ is varied

b)

Fig. 4The convergent dynamical characteristic of the system when parameter f is varied

The convergent dynamical characteristic of the system when parameter f is varied

a)

The convergent dynamical characteristic of the system when parameter f is varied

b)

(2) The excitation amplitude parameter f is varied. When the parameters ξ=0.1, γ= –1, f= 3, ω~= 3 and the parameter f is decreased from 10 to 3, the dynamic behavior still approaches to two stable state, as shown in Fig. 4(a). But when the f=0.3, it is interesting that the system approaches to one stable regime only, as shown in Fig. 4(b), which means that the system oscillates in one potential well.

(3) The excitation frequency parameter ω~ is varied. When ξ=0.1, γ=-1, f=0.4, ω~=2 and the frequency is decreased from 3 to 2, and system works in a periodic state and approaches to one stable regime, as shown in Fig. 5(a). When the f is varied from 0.4 to 1, the system approaches to two stable domains, as shown in Fig. 5(b).

Fig. 5The convergent dynamical characteristic of the system when parameter ω~ is varied

The convergent dynamical characteristic of the system when parameter ω~ is varied

a)

The convergent dynamical characteristic of the system when parameter ω~ is varied

b)

2.4. Lyapunov exponent curve analysis

For system (1), the Lyapunov exponent spectrum combined with the bifurcation analysis is an effective way to judge whether the system works in a chaotic or not. If the Lyapunov is positive, the system works in a chaotic state. When the parameters are ξ=0.1, ω~=3, γ= –1, and the excitation force amplitude f is varied from 0.1 to 10, the bifurcation diagram is obtained in Fig. 6.

In Fig. 6, there are abundant dynamical behaviors being revealed. When the excitation force amplitude f is small, the system exhibits period-1 response, and the double-period bifurcation happens at the f is 0.8. The multiple periodic motion and chaos interleave in the range from 2 to 5. In order to identify the chaos parameter regime precisely, the Lyapunov exponent spectrum is calculated, as shown in Fig. 7.

Fig. 6The global bifurcation of the system with the variation of f

The global bifurcation of the system  with the variation of f

Fig. 7The Lyapunov exponent spectrum with the variation of the f

The Lyapunov exponent spectrum  with the variation of the f

In this figure, the chaos parameters can be determined on the abscissa axis. Some cases are selected to discover the different dynamic characteristics of the system. When f is 0.1, the system exhibits period-1 motion, the Poincaré section is one point on the phase plane, and the limit cycle is observed, as shown in Fig. 8(a).

When f is 3.4, the system works in a chaotic state, and the strange attractor is given as shown in Fig. 8(b). The corresponding Lypunov exponent is 0.035, and it is positive one, which is the index for chaos. When the f is increase to 3.5, the period-2 motion is presented, and two Poincaré section points lie on the phase plot, as shown in Fig. 8(c). When the f is increase to 4.4, the system comes back to chaos again, as shown in Fig. 8(d).

In order to verify that the calculation is accurate, the Lyapunov exponent spectrum characteristic is checked. The system (1) is rewritten as the following form:

19
z˙1=z2,z˙2=-ξz2-γz1-βz13+fcosz3,z˙3=ω~.

Fig. 8The different response phase plane diagram when the parameter f is varied

The different response phase plane diagram when the parameter f is varied

a) The limit cycle when f is 0.1

The different response phase plane diagram when the parameter f is varied

b) Chaotic response when f is 3.4

The different response phase plane diagram when the parameter f is varied

c) The period-2 response when f is 3.5

The different response phase plane diagram when the parameter f is varied

d) Chaotic response when f is 4.4

According to the phase space flow evolution, the divergence of f is get:

divf=f1z1+f2z2+f3z3=-ξ,

where f=(f1,f2,f3)T is composed of right hand side components of Eq. (19). Thus, the initial phase space volume V(0) will converge in the form of V(t)=V(0)e-ξt. In essence, it is also considered as that the volume converges in the different phase space directions with the law of V(t)=V(0)e(λ1+λ2+λ3)t. Therefore, the Lyapunov exponent spectrum λi(i=1, 2, 3) characterizes the convergent behavior of the system, and it satisfied with i=13λi=-ξ. As an example, parts of Lyapunov exponents are listed in table1. Because z˙3 keeps constant, the Lyapunov exponent λ3 is 0. The summation of the other two exponents is –0.1. It meets the elationship between Lyapunov exponent spectrum and damping coefficient.

Table 1Lyapunov spectrum when the excitation amplitude f is varied from 0.5 to 0.8

Excitation amplitude
0.5
0.6
0.7
0.8
Lyapunov exponent λ1
–0.0462
–0.0451
–0.0436
–0.0393
Lyapunov exponent λ2
–0.0538
–0.0549
–0.0564
–0.0607

When the excitation amplitude is fixed to 4.4, and the parameter ω~ is increased from 0.1 to 10, the global bifurcation diagram is obtained, as shown in Fig. 9. The chaotic vibration occurs in the low frequency regime, and the multi-periodic motion exhibits in high frequency regime. Therefore, in the experiment, the excitation frequency should be adjusted to low frequency range to observe chaos.

In order to get the global dynamics of the system, the comprehensive response characteristic diagram is computed under different excitation force and frequency, as shown in Fig. 10. The black is period-1 motion; the yellow is period-2 motion; the blue is period-3 motion, the magenta is period-4 motion; the red is chaotic motion, and the green is the other kinds of multiple-period motions. Seen from the figure, the chaotic behavior occurs close to the system’s natural frequency, 3rd sup-harmonic, and half harmonic frequencies regime.

Fig. 9the different response when the parameter ω~ is varied

the different response when  the parameter ω~ is varied

Fig. 10The global comprehensive response of the system when the parameters f and ω~ are varied

The global comprehensive response of the system when the parameters f and ω~ are varied

3. Finite element analysis of leaf spring

A finite element model (FE) was then used to numerically investigate the natural frequency and vibration mode of the leaf spring system.

Geometrical parameters and material properties of the leaf spring system are shown in Tables 2 and 3. The excitation, a longitudinal force and a transversal disturbance with 30 % strength of the longitudinal force, was applied at one end of the leaf spring (with the other end restricted). Details of excitation are shown in Table 4.

Table 2Geometrical parameters of leaf spring

Parameters
Symbol
Value
Length
L
180 mm
Width
w
20 mm
Thickness
d
1.5 mm

Table 3Material properties of leaf spring

Parameters
Symbol
Value
Young’s modulus
E
1.1e11 Pa
Possion ratio
ν
0.35
Density
ρ
8700 kg/m3

Table 4Excitation force details

Parameters
Symbol
Value
Longitudinal force
Fz
0 to 30 N
Transversal force
Fy
0 to 9 N (30 % of longitudinal force)
Normalized longitudinal force
NFC
0 to 1

A four-node rectangular element was used with a maximum mesh size of 0.054 mm. The geometry and mesh of the system are shown in Fig. 11. Displacements on y and z directions at the free end of the leaf spring are shown in Fig. 12, against the value of normalized longitudinal force.

From Fig. 12, the displacement in the y direction is a concave curve, which indicates a nonlinear buckling and two possible positions on the curve nearby the point P. This means the system can potentially cause chaotic states.

When the leaf spring is installed as a cantilever, the first 6 order natural frequencies obtained from the FE model are shown in Table 5, and the relevant vibration modes are shown in Fig. 13.

Fig. 11Geometry and mesh of the leaf spring system

Geometry and mesh  of the leaf spring system
Geometry and mesh  of the leaf spring system

Fig. 12Displacement of the leaf spring in y and z direction

Displacement of the leaf spring  in y and z direction

To drive the system entering chaotic vibration regime, the excitation frequency should be close to the system’s natural frequency or its harmonic/half harmonic frequencies. Since the system’s first order natural frequency is 26.664 Hz, the excitation frequency ranges in later experimental study will be focused on the range from 5 Hz to 27 Hz which covers the first order natural frequency and its half harmonics.

Table 5The first 6 order natural frequencies of leaf spring

Order
Natural frequency (Hz)
Order
Natural frequency (Hz)
1st
26.664
4th
935.86
2nd
167.75
5th
1567.4
3th
472.85
6th
2379.5

4. Test rig set-up

The test rig was designed based on the double potential well principle. As seen in Fig. 14(a).

It has 8 major components, including:

1) Exciter (JZK-50) which provides single frequency excitation to the external frame. Frequency of the excitation was controlled by signal generator (YE1311). Strength of the excitation was controlled by power amplifier (YE5872A);

2) Support base which is used to keep the exciter in a horizontal state;

3) Leaf spring which has the same geometry and material properties as shown in Tables 2 and 3;

4) Rectangular magnet which is used to produce external force;

5) Iron tip which is a small piece of iron. It was placed on the tip of the leaf spring and can be attracted by the magnets;

6) External frame which is used to install magnets and connect exciter;

7) Clamp which is used to fix the end of leaf spring;

8) Laser vibrometer (OPTEX CD33-30NV) which is used to measure the displacement of the leaf spring. It can measure maximum displacement of ±4 mm, with a maximum error of ±0.1 % F.S. Distance between the vibrometer and the leaf spring is 30 mm. The laser vibrometer is connected to a data acquisition card (NI4472).

Fig. 13Different vibration modes of the leaf spring

Different vibration modes of the leaf spring

a) The first-order vibration shape

Different vibration modes of the leaf spring

b) The second-order vibration shape

Different vibration modes of the leaf spring

c) The third-order vibration shape

Different vibration modes of the leaf spring

d) The fourth-order vibration shape

Different vibration modes of the leaf spring

e) The fifth-order vibration shape

Different vibration modes of the leaf spring

f) The sixth-order vibration shape

LabVIEW was used to control the DAQ card and online monitor the time history and frequency spectrum. During the experiment, the vibration signal was recorded by the DAQ card with a sampling rate of 2 kHz and a total record time of 5 s for each record. The excitation frequency was adjusted from 5 Hz to 27 Hz with a step of 0.1 Hz. Measurement was repeated under each frequency. Recorded signals were analysed in both time domain and frequency domain. Characteristic frequencies were identified from the frequency spectrum. After phase space reconstruction, the characteristic exponents of the measured signals were obtained with the nonlinear time series analysis method. The chaos identification code was programmed using newly developed signal processing techniques, including improved LE exponent [33] Poincaré section [34], and Wavelet-SG-EEMD, which will be discussed in Section 5.

5. Wavelet-SG-EEMD denoising model

In experiment, recorded chaotic vibration signals are always contaminated by noise from background, electronics and local disturbances. Since the chaotic vibration signals are non-stationary and nonlinear, traditional linear filtering methods are not effective. Therefore, a new combined denoising model is used here. This model uses a Wavelet-SG algorithm which was developed based on the wavelet transform [30] and Savitzky-Golay (SG) theory [32] as pre-filter, and used Ensemble Empirical Mode Decomposition(EEMD) [31] method to reduce the noise and local disturbances.

Fig. 14Leaf spring chaotic vibration experimental rig. 1. Shaker; 2. Support base; 3. Leaf spring; 4. Magnet; 5. Iron tip; 6. External frame; 7. Clamp; 8. Laser vibrometer

Leaf spring chaotic vibration experimental rig. 1. Shaker; 2. Support base; 3. Leaf spring; 4. Magnet; 5. Iron tip; 6. External frame; 7. Clamp; 8. Laser vibrometer

a) Sketch of the experimental rig

Leaf spring chaotic vibration experimental rig. 1. Shaker; 2. Support base; 3. Leaf spring; 4. Magnet; 5. Iron tip; 6. External frame; 7. Clamp; 8. Laser vibrometer

b) Photo of the experimental rig

Fig. 15Flow chart of denoising procedure of Wavelet-SG-EEMD method

Flow chart of denoising procedure of Wavelet-SG-EEMD method

The signal processing flow chart of using this Wavelet-SG-EEMD denoising model is shown in Fig. 15, and the detailed procedures are as follows.

Obtain the signal’s wavelet coefficients of each scale by applying wavelet decomposition:

20
cj,k=ncj-1,nhn-2k,dj,k=ndj-1,ngn-2k, k=0,,N-1,

where cj,k and dj,k are approximated coefficients and detailed coefficients respectively, h and g are impulse responses of the filter, and j is the decomposition level of the response.

Reduce the noise using SG method which is presented in Eq. (21). When si is one of the IMF (Intrinsic Mode Function) coefficients, a polynomial psi with M the order is fitted by nl+nr+1 points in the vicinity of si, among which nl is the number of the left point of si, and nr is the number of the right point of si. The smooth value qiis:

21
qi=k=0Mbks-siΔspsi,

where Δs=si+1-si; bk is the coefficient of psi. Suppose the measurement data is yi, in order to fit the testing data with psi, the coefficient bkis determined by optimized function:

22
mini=1i+nr[pi(sj)-yj]2.

Keep the wavelet coefficients in low frequency band unchanged, process the detailed wavelet coefficients with Eq. (23):

23
dj'(i)=sgndjidji-Tj, djiTj,0, dj(i)<Tj,

where the Tj is the threshold at the jth level of wavelet decomposition. Because the variance and the amplitude of the noise will decrease with the increase of the decomposition level, a variable threshold is set as Tj=σ2logn/j+1 depending on the variance σ of noise, data length n and the decomposition level j, and thus details of the useful signals are preserved as much as possible.

Reconstruct signal with the approximated and detailed coefficients with wavelet reconstruction algorithm Eq. (24):

24
cj-1,n=ncj,nhk-2n+ndj,ngk-2n.

Add random white noise nit to reconstructed signal x'(t):

25
xit=x't+nit,

where i=1, 2, 3,…, M, M>1.

Decompose the signal xit into IMFs aijt (j= 1,…, N) through the EMD (Empirical Mode Decomposition), where aijt denotes the ith IMF of the jth sample, and N is the number of IMFs. If j<M then go to the step (5) with j=j+1. Repeat steps (5) and (6) with different white noise time series at each time.

Finally, calculate the ensemble average ajt of the M samples for IMF, and reconstruct the original signal:

26
ajt=1Mi=1Maijt.

A noisy Lorenz time series is used to test performance of the proposed denoising model. The SNR of Lorenz is set to 5 dB, and the length is 1000. As shown in Figs. 16(a), (b), the proposed denoising model can significantly reduce the noise from chaotic signal. The measured signal of experiment is processed in Figs. 16(c), (d).

Fig. 16The denoising result of Wavelet-SG-EEMD method

The denoising result of Wavelet-SG-EEMD method

a) The noisy Lorenz time series

The denoising result of Wavelet-SG-EEMD method

b) The denoising result

The denoising result of Wavelet-SG-EEMD method

c) Reconstructed attractor on phase plane with original measured signal

The denoising result of Wavelet-SG-EEMD method

d) Reconstructed attractor on phase plane with denoised signal

6. Experimental results analysis and discussion

6.1. Vibration responses under different strength of excitation.

The excitation frequencies studied in the experiment are from 5 Hz to 27 Hz and the output voltage of power amplifier (OVA)is up to1.5 V. The frequency was adjusted with a step of 0.1 Hz. The output voltage of power amplifier was then adjusted from 0 to 1.5 V under each frequency step.

The iron tip falls into the right potential well and oscillates in a period-1 mode (a limit cycle in phase space), when the excitation frequency is 15.6 Hz and output voltage of amplifiers 0.1 V. As shown in Fig. 17(a), the amplitude of the vibration is 0.056 mm. Fundamental frequency (15.6 Hz) and third harmonics (46.8 Hz) are observed, in Fig. 17(b). The reconstruction attractor is a limit cycle, as seen in Fig. 17(c). If a minor disturbance is added to the system, it can be observed from the experiment that the iron tip jumps out of the well and vibrates between two wells. It falls into the left potential well and works in a periodic state after a while, which indicated that the external excitation is not enough to push it out of the well barrier and it is constrained in the potential well fully.

When the excitation is 15.6 Hz and the output of power amplifier is increased to 0.5 V, the vibration amplitude increased to 0.1721 mm. The system shows similar responds, as shown in Fig. 18.

Fig. 17Vibration response when system OVA is 0.1 V

Vibration response when system OVA is 0.1 V

a) Time history of the response

Vibration response when system OVA is 0.1 V

b) Frequency spectrum of the response

Vibration response when system OVA is 0.1 V

c) Reconstructed periodic-1 attractor

Vibration response when system OVA is 0.1 V

d) Poincaré section on limit cycle

Fig. 18Vibration response when system OVA is 0.5 V

Vibration response when system OVA is 0.5 V

a) Time history of the response

Vibration response when system OVA is 0.5 V

b) Frequency spectrum of the response

Vibration response when system OVA is 0.5 V

c) Reconstructed period-1 attractor

Vibration response when system OVA is 0.5 V

d) Poincaré section of cycle

The limit cycle is reconstructed and an improved Poincaré section algorithm is applied to analyze the signal as shown in Fig. 18(c). In order to identify the chaos and periodic mode accurately, at least one point must be located on the limit cycle. However due to the uncertainty from experimentation, the location of the point in the limit cycle may vary slightly. As a result, a number of points are searched in the potential regime as shown in Fig. 18(d) using modified Poincaré method [34].

Increasing the output of power amplifier to 0.8 V, the system reaches a critical state where after the same minor disturbance is applied, the system enters chaotic vibration mode which shows that the iron tip jumps between left and right wells randomly. However, the system can still return back to period-1 where the iron tip stays in one of the wells without any disturbance.

When the output voltage of power amplifier exceeds the critical level, reaches 0.95 V, half harmonics (7.83 Hz) can be recognized Fig. 19(b). This indicates that the doubling periodic is happening in the system.

The system enters into a stable chaotic state automatically even without any external disturbance when the output of power amplifier reaches 0.98 V. As shown in Fig. 20(a), the time history is random-like motion and the power spectrum is broadband. The reconstructed phase space parameters are calculated with FNN (false nearest neighbors) and FMMI (first minimum mutual information) methods. The embedding dimension and delay parameters are 4 and 16 respectively, as shown in Figs. 20(c), (d). The chaotic attractor is obtained in Fig. 20(e), and the maximal LE [33] is 0.313.

Fig. 19Vibration response when system OVA is 0.8 V

Vibration response when system OVA is 0.8 V

a) Time history of the response

Vibration response when system OVA is 0.8 V

b) Time history of the response

Fig. 20Vibration response when system OVA is 0.98 V

Vibration response when system OVA is 0.98 V

a) Time history of the response

Vibration response when system OVA is 0.98 V

b) Frequency spectrum of the response

Vibration response when system OVA is 0.98 V

c) FNN embedding dimension

Vibration response when system OVA is 0.98 V

d) FMI reconstruction delay

Vibration response when system OVA is 0.98 V

e) Reconstructed chaotic attractor of system

6.2. Vibration responses under different excitation frequencies

When the output of the power amplifier was less than 0.5 V, changing the excitation frequency from 10 Hz to 17 Hz, the system responses in all the cases stay at period-1 mode. However, when the output voltage of the power amplifier was changed to 1.0 V, the system shows much more complicated responses when frequency varies.

It works in chaotic pattern when the excitation frequencies are in the range of 10 Hz-16 Hz and 17 Hz-19 Hz. When the frequency is 13.4 Hz and system OVA is 1.0 V, the chaotic response is shown in Fig. 21.

However, there is a certain range, between 16.1 and 16.9 Hz, (e.g. 16.6 Hz), where it works in period-2 mode as shown in Fig. 22, and the LE is –0.0085.

Fig. 21Chaotic vibration response of the system

Chaotic vibration response of the system

a) Time history of the response

Chaotic vibration response of the system

b) Frequency spectrum of the response

Chaotic vibration response of the system

c) FNN embedding dimension

Chaotic vibration response of the system

d) FMI reconstruction delay

Chaotic vibration response of the system

e) Reconstructed chaotic attractor of system

Fig. 22Period-2 vibration response of the system

Period-2 vibration response of the system

a) Time history of the response

Period-2 vibration response of the system

b) Frequency spectrum of the response

Period-2 vibration response of the system

c) FNN embedding dimension

Period-2 vibration response of the system

d) FMI reconstruction delay time

Period-2 vibration response of the system

e) Reconstructed period-2 attractor

Period-2 vibration response of the system

f) Poincaré section on limit cycle

If the frequency is further increased continuously, the system state is changed from chaos to periodic motion. At 19 Hz, the period-3 motion occurs Fig. 23.

Furthermore, when the excitation frequency varied from 19 Hz to 23 Hz, large amplitude multiple periodical vibration between two magnets were observed. This means the iron tip can jump over the potential well peak under these frequencies.

At 23.0 Hz, the multi-periodic motion between double potential wells is shown in Fig. 24 and the vibration amplitude is increased to 1.257 mm.

Fig. 23Period-3 vibration response of the system

Period-3 vibration response of the system

a) Time history of the response

Period-3 vibration response of the system

b) Frequency spectrum of the response

Period-3 vibration response of the system

c) FNN embedding dimension

Period-3 vibration response of the system

d) FMI reconstruction delay time

Period-3 vibration response of the system

e) Reconstructed period-3 attractor

Fig. 24Multi-periodic vibration response between double potential wells

Multi-periodic vibration response between double potential wells

a) Time history of the response

Multi-periodic vibration response between double potential wells

b) Reconstructed multi-periodic attractor

7. Conclusions

A theoretical and experimental study was carried out to investigate the chaotic vibration characteristics of a double potential well leaf spring and magnet system. The nonlinear amplitude frequency characteristic curve is obtained with averaging method [35], and the Lyapunov exponent spectrum combined with the global bifurcation diagram is applied to determine the chaotic regime of the system. A test rig was designed based on the double potential well theory to generate stable and repeatable chaos. The relationship between excitation (strength and frequency) and the leaf system’s nonlinear status was carefully studied. During the experiment, sub-harmonic phenomenon which is an indication of chaos was identified. Some interesting dynamic behaviors occurred, including the two patterns of period-1 motion: one is oscillation in any of the unilateral potential wells; the other is vibration between two wells. The critical state seems to be very sensitive to the excitation. The range of excitation amplitude and frequency where stable and repeatable chaos can be generated were obtained. The signal processing methods, Wavelet-SG-EEMD denoising algorithm and improved Poincaré section method were tested and proved to be effective in identifying chaotic behavior from the system.

References

  • Samuel Zambrano, Juan Sabuco How to minimize the control frequency to sustain transient chaos using partial control. Communications in Nonlinear Science and Numerical Simulation, Vol. 19, Issue 3, 2014, p. 726-737.
  • Gao Yuan, Geng Zhaoyun, Fan Jianwen Control chaos in automobile suspension system via adaptive tracking control method. Machinery Design and Manufacture, Vol. 12, 2013, p. 198-201.
  • Han Jianqun, Zheng Ping Active control isolation for the Lorenz chaotic vibration with a simple method. Systems Engineering and Electronics, Vol. 128, Issue 10, 2006, p. 1566-1568.
  • Long Yunjia, Yang Yong, Wang Congling Road roller engineering based on chaotic vibration mechanics. Engineering Science, Vol. 2, Issue 9, 2000, p. 76-79.
  • Dmitriev Panas Starkov Experiments on speech and music signals transmission using chaos. International Journal of Bifurcation and Chaos, Vol. 5, Issue 4, 1995, p. 1249-1254.
  • Govindan R. B., Narayanan K., Gopinathan M. S. On the evidence of deterministic chaos in ECG: Surrogate and predictability analysis. Chaos, Vol. 8, Issue 2, 1998, p. 495-502.
  • Omar Mohamed A., Shabana Ahmed A., Mikkola Aki Multibody system modeling of leaf springs. Journal of Vibration and Control, Vol. 10, Issue 11, 2004, p. 1601-1638.
  • Qing Li, Wei Li A contact finite element algorithm for the multileaf spring of vehicle suspension systems. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, Vol. 128, Issue 3, 2004, p. 305-314.
  • Yum Yong-Jin Frictional behavior of automotive leaf spring. Science and Technology, Proceedings of the 4th Korea-Russia International Symposium, Vol. 3, 2000, p. 5-10.
  • Sugiyama Hiroyuki, Shabana Ahmed A., Omar Mohamed A., Loh Wei-Yi Development of nonlinear elastic leaf spring modelfor multibody vehicle systems. Computer Methods in Applied Mechanics and Engineering, Vol. 195, 2006, p. 6925-6941.
  • Rahman Muhammad Ashiqur, Siddiqui Muhammad Tareq, Kowser Muhammad Arefin Design and non-linear analysis of a parabolic leaf spring. Journal of Mechanical Engineering, Vol. 37, 2007, p. 47-51.
  • Dwivedy S. K., Kar R. C. Nonlinear dynamics of a cantilever beam carrying an attached mass with 1:3:9 internal resonances. Nonlinear Dynamics, Vol. 31, 2003, p. 49-72.
  • Nayfeh Ali H., Pai Perngjin F. Non-linear non-planar parametric responses of an inextensional beam int. Journal of Non-Linear Mechanics, Vol. 24, Issue 2, 1989, p. 139-158.
  • Li Chen An integral approach for large deflection cantilever beams. International. Journal of Non-Linear Mechanics, Vol. 45, 2010, p. 301-305.
  • Nallathambi Ashok Kumar, Rao Lakshmana C., Srinivasan Sivakumar M. Large deflection of constant curvature cantilever beam under follower load. International Journal of Mechanical Sciences, Vol. 52, 2010, p. 440-445.
  • Doyle J. F. Non-Linear Analysis of Thin-Walled Structures-Statics Dynamics, and Stability. Springer, New York, 2001.
  • Dado M., Al-Sadder S. A new technique for large deflection analysis of non-prismatic cantilever beams. Mechanics Research Communications, Vol. 32, 2005, p. 692-703.
  • Baker G. On the large deflections of non-prismatic cantilevers with a finite depth. Computers and Structures Vol. 46, 1993, p. 365-370.
  • Wang C. M., Kitipornchai S. Shooting optimization technique for large deflection analysis of structural members. Engineering Structures, Vol. 14, 1992, p. 231-240.
  • Srpcic S., Saje M. Large deformations of thin curved plane beam of constant initial curvature. International Journal of Mechanical Sciences, Vol. 28, 1986, p. 275-287.
  • Argyris J. H. Non-linear finite element analysis of elastic systems under nonconservative loading-natural formulation. Part 1: Quasi static problems. Computer Methods in Applied Mechanics and Engineering, Vol. 26, 1981, p. 75-123.
  • Nayfeh Ali H., Pai Perngjin F. Non-linear non-planar parametric responses of an inextensional beam. International Journal of Non-Linear Mechanics, Vol. 24, Issue 2, 1989, p. 139-158.
  • Li S., Yang S., Guo W. Investigation on chaotic motion in hysteretic non-linear suspension system with multi-frequency excitations. Mechanics Research Communications. Vol. 31, 2004, p. 229-236.
  • Thompson J. M. T. Chaotic phenomena triggering the escape from a potential well. Proceedings of the Royal Society of London A, Vol. 421, 1989, p. 195-225.
  • Trueba Jos E. L., Rams Joaquin, Sanjuan Miguel A. F. Analytical estimates of the effect of nonlinear damping in some nonlinear oscillators. International Journal of Bifurcation and Chaos, Vol. 10, Issue 9, 2000, p. 2257-2267.
  • Litak Grzegorz, Borowiec Marek, Friswell Michael I., Szabelski Kazimierz Chaotic vibration of a quarter-car model excited by the road surface profile. Communications in Nonlinear Science and Numerical Simulation, Vol. 13, 2008, p. 1373-1383.
  • Liang Shan, Zheng Jian, Zhu Qin, Liu Fei Numerical and experimental investigations on chaotic vibration of a nonlinear vehicle model over road excitation. Journal of Mechanical Strength, Vol. 34, Issue 1, 2012, p. 6-12.
  • Moon Francis C. Chaotic Vibrations an Introduction for Applied Scientists and Engineers. New York, 1987.
  • Fossas Enric, Olivar Gerard Study of chaos in the buck converter. IEEE Transactions on Circuits and Systems, Vol. 43, Issue 1, 1996, p. 13-25.
  • Chiementin X., Kilundu B., Rasolofondraibe L., et al. Performance of wavelet denoising in vibration analysis: highlighting. Journal of Vibration and Control, Vol. 18, Issue 6, 2012, p. 850-858.
  • Chen R. X., Tang B. P., Ma J. H. Adaptive de-noising method based on ensemble empirical mode decomposition for vibration signal. Journal of Vibration and Shock, Vol. 31, Issue 15, 2012, p. 82-86.
  • Li K., Yang S. Q. Image smooth denoising based on Savitaky-Golay algorithm. Journal of Data Acquisition and Processing, Vol. 25, 2010, p. 72-74.
  • Yang Aibo, Wang Ji, Liu Shuyong, Wei Xiulei An algorithm for computing the largest Lyapunov exponent based on space grid method. Acta Electronica Sinica, Vol. 40, Issue 9, 2012, p. 1871-1875.
  • Liu Shuyong, Zhu Shijian, Yang Qingchao Study on the application of poincaré section algorithm to chaotic vibration identification. Applied Physics Vol. 1, 2011, p. 108-115.
  • Sanders J. A., Verhulst F. Averaging Methods in Nonlinear Dynamical Systems. Springer-Verlag. New York, 1985.

About this article

Received
18 November 2016
Accepted
27 December 2016
Published
15 May 2017
SUBJECTS
Chaos, nonlinear dynamics and applications
Keywords
chaotic vibration
double potential well
nonlinear system
nonlinear time series
Acknowledgements

This research is supported by the National Natural Science Foundation China (51179197, 51579242, 51509253) and the Naval University of Engineering Foundation (425517K143).

Author Contributions

Shuyong Liu contribution is mathematic model analysis. Jian Jiang – mathematic model simulation. Pan Su – experimental design. JieChang Wu – experimental analysis and signal measuring. Ye Liu – experimental design and signal analysis. Yuan Fang – signal analysis.