Nonlinear Dynamics Analysis of the Semiactive Suspension System with Magneto-Rheological Damper

This paper examines dynamical behavior of a nonlinear oscillator which models a quarter-car forced by the road profile. The magneto-rheological (MR) suspension system has been established, by employing the modified Bouc-Wen force-velocity (F-v) model of magneto-rheological damper (MRD). The possibility of chaotic motions in MR suspension is discovered by employing the method of nonlinear stability analysis. With the bifurcation diagrams and corresponding Lyapunov exponent (LE) spectrum diagrams detected through numerical calculation, we can observe the complex dynamical behaviors and oscillating mechanism of alternating periodic oscillations, quasiperiodic oscillations, and chaotic oscillations with different profiles of road excitation, as well as the dynamical evolutions to chaos through period-doubling bifurcations, saddle-node bifurcations, and reverse period-doubling bifurcations.


Introduction
Magneto-rheological fluid (MRF) is a suspension of micronsized, magnetic particles in a carrier fluid.When exposed to a magnetic field, the rheology of MRF reversibly and consecutive changes from a free-flowing Newton liquid to a semisolid Bingham with controllable yield strength, which is known as magneto-rheological effect [1][2][3].Magnetorheological damper (MRD) based on MRF has significant promise for effective vibration damping in many applications, such as MRD-based semiactive suspension system, which has attracted much attention to improve the ride comfort and handling safety of the vehicle [4][5][6].Nevertheless, there are still numerous challenges in controlling a MRD to get the superior performances, because it has highly nonlinear characteristics and chaotic motion due to typical hysteretic characteristics.
In the past decades, many studies of MR suspension focused on the semiactive control algorithms.Many new ideas were proposed, such as "Skyhook" control or artificial intelligence based control, and had been implemented in a practical vehicle.However, chaotic behavior may exist due to the nonlinearities in MR suspension, which was not studied well yet.Li et al. [7] investigated the chaotic motion in nonlinear suspension system with hysteretic characteristics and verified the path from quasiperiodic to chaos by deriving Melnikov method.Siewe Siewe [8] applied the method of multiple scales to analyze local bifurcation in the quartercar system with periodically excited road profile, and a variety of nonlinear behaviors were found, such as resonance and antiresonance phenomena and saddle-node bifurcation.Litak et al. [9] used the analytical Melnikov theory and predicted the lowest critical amplitude that a single degree of freedom (DoF) vehicle model may transit to a chaotic motion, under a road surface profile consisting of harmonic and noisy components.Luo and Rajendran [10] carried out the periodic motion and stability of a single DoF semiactive suspension model by developing a mapping structure, and 2 Shock and Vibration the model of MRD was formulated with piecewise linear equations.However, all of above models are simple as single DoF model which is very different from the practical situation owing to the neglect of wheel movement, and this may explain that the vibrations appear for impossible amplitude of 0.4 m in [9].Consequently, Borowiec and Litak [11] studied a 2-DoF quarter-car and find out the transition to chaos, by applying Melnikov theory and recurrence approach.In earlier studies, the suspension system was mostly considered as single DoF system for its simplicity.What is more, calculation models of MRD were mostly made employing polynomials or piecewise linear model in [9], while the Bouc-Wen model is considered appropriate to describe well the dynamic performance of MRD [12][13][14].So far, there are no systematic nonlinear dynamics of the Bouc-Wen model based MR suspension system, owing to the complex structure of model.
In this paper, from the experimental results, such a commercial MRD is modeled using a modified Bouc-Wen model proposed by the author.A 2-DoF model is established to express the MR suspension system, by employing the identification results.The stability of the system is analyzed according to the stability criterion, and all possible motions of the MR suspension system are determined by plotting the frequency response, bifurcation diagrams, and phase portraits under different road profiles.The Lyapunov exponent (LE) is calculated to detect chaotic motion.Time series with combination of power spectrum density is used as assisted means for the special system parameters.
The paper is organized as follows.In Section 2, the dynamic model of 2-DoF MR suspension system is formulated, employing the modified Bouc-Wen calculation model of MRD.In Section 3, the stability of the system is analyzed by calculating the eigenvalue of the Jacobian matrix of fixed point.Next, the numerical calculation is conducted to determine the dynamic behavior and to confirm the stability analysis as well.Comprehensive numerical results include frequency response, bifurcation diagrams, phase plane portraits, Poincare map, and time series; thus, the process of the transition to chaotic motion is revealed.Finally, the conclusions of the research are presented in Section 4.

Dynamic Model of Quarter-Car Semiactive Suspension System
2.1.Mechanical Model and Formulation.In Figure 1, a classic dynamic model of a quarter MR suspension system [15] [14].Figure 2 shows the structure of typical Bouc-Wen phenomenon model, and it accurately describes the MRD inherent hysteresis nonlinear properties.However, it cannot describe nonlinear response and saturation characteristics of the magnetic field, because the linear output item in the model is just to represent the relationship between the damping force and the control current.Therefore, the author has imported the sigmoid function to improve the conventional Bouc-Wen phenomenon model [16].The issue is effectively solved by decoupling the hysteretic characteristics with the current modulation separation: where   and   express the direct current and its maximum value for driving the MRD, respectively, and 0 ≤   ≤   .(  ) denotes the saturated nonlinear direct current control function proposed by the authors and (  ) ≥ 1, (  ) = 1, for   = 0.  ℎ (V  ) denotes yielded passive damping force with hysteresis depending on the piston relative displacement velocity (V  ) of the MRD for   = 0.   is the piston travel of the MRD,  and  are inner variables without units, and  0 ,  1 ,  2 ,  0 ,  0 , , , ,  0 ,  1 , , , and  0 are constants, respectively.
As is shown in Figure 3(a), a CARRERA MagneShock MRD of the vehicle suspension is further employed from [16], which permits maximum control direct current 0.5 A at 12 V [16].On the basis of the actual measured characteristic data of MRD, the model parameters are identified as  0 = 184.1,MRD has clear strong nonlinear property from the above results.Due to the damping force dependency on driven current and excitation frequency, the -V curve varies evidently under different driven and excitation profiles [17].As a result, nonlinear behavior caused by hysteresis cannot be ignored in engineering application.

Nonlinear Dynamics Characteristics of MR Suspension
Stability is undoubtedly the key requirements in the field of vehicle suspension control.In previous studies, comprehensive performances of the MR suspension were more highly valued, by comparing the transmissibility of vibration and vibration amplitude with the traditional passive one.Here, the nonlinear dynamic response and the mechanism of the instability process are the problem to be solved.
In this section, the stability analysis of the MR suspension is conducted in view of modern control theory and nonlinear dynamics theory.Stability of the MR suspension system is examined at the equilibrium of the dimensionless nonlinear model.Further, the bifurcation analysis, combining LE assessment, is carried out to describe nonlinear dynamical evolution process and to find out the route to chaos.In addition, the phase portrait and Poincare map diagrams are plotted to give the more intuitive response of the system under different excitation conditions.
Furthermore, due to quadratic and absolute terms in the above six dimensions system, analytical solution is scarcely to obtain.Therefore, the numerical methods are usually employed to analyze such a system.Considering the application background, the range of frequencies and amplitudes of the road profile during vehicle running is certain, which are normally below 15 Hz and 10 cm, respectively [4][5][6].We use the common harmonic excitation as the road surface.It is expressed as   =   sin(Ω ⋅ ), in which   represents roughness of the road surface and Ω is angular frequency Ω = 2.The dimensionless form of road excitation is   = Ω 2  sin(Ω/⋅).The following analysis methods of nonlinear dynamics are applied to such a special complex system: (1) By obtaining the frequency band response, we find out the area that the system is sensitive to the corresponding road profiles.
(2) By drawing the bifurcation diagrams, we study the nonlinear dynamics of the system under different road parameters.Combing the calculated LE spectrum, we can further determine the chaotic motion of the system under the corresponding road profiles.
(3) Based on the above analysis, the dynamics evolution process is vividly portrayed using phase plane portraits, time series, and power spectrum of the system response at critical parameters.

Numerical Results.
It is known that the dynamics of a vibration system may be analyzed through the frequency response diagram [18].Therefore, for the studied system, the frequency response is obtained by plotting the vibration amplitude of sprung mass   expressed by V  .In addition, the LE spectrum [19] is used to reveal the detail of the system frequency response.Figure 4(a) shows the frequency response of the model, which covers pass-band of the road, and Figure 4(b) presents the LE spectrum as well.The road frequency  is slowly increased by increment of 0.001 Hz.As is illustrated in Figure 4(a) there exists a critical jump of the system response for  = 1.752Hz [20] near the resonance point.The phenomenon of the jump causes the motion to change, but the system remains stable, as is shown in gray shadow.However, the diagram exhibits a more complicated and different behavior while  is increased to 2.71 Hz, with that restoration of stability in a short time.With  being increased, it is shown that the system falls into instability area as 2.75 Hz-4 Hz, which indicates that the chaotic motion may appear when  is within or near this area, as is shown in red shadow.This is confirmed by the LE spectrum showed in Figure 4(b).The diagram illustrates that there exist positive LE among the LE spectra in instable area shaded in red, which indicates the existence of the chaos.Next, the frequency response goes back to normal and all the LEs are less than zero, which indicates that the system remains stable up to 15 Hz.The bifurcation and max LE diagrams under parameter variations are efficient methods for analyzing the nonlinear dynamic behavior.Figure 5 shows the bifurcation diagram and the corresponding LE diagram, with  varying in the above-mentioned instable area.The bifurcation diagram is obtained by plotting the stroboscopic point of the displacement   , because   is important for the safety of the vehicle.As is shown in Figure 5(a), when  = 2.08 Hz, the system loses period-1 stability, bifurcating from period-1 into period-2 motion.For  ∈ (2.08-2.74)Hz, the period-2 motion of the system develops into period-8 through the perioddoubling bifurcation.Moreover, it is interesting that complex motion occurs at  ∈ (2.439-2.488)from the bifurcation diagram, and the max LEs are positive.At the same time, note that the points in the bifurcation diagram are limited in enveloping curve, not wide distribution.Consequently, the exact dynamic behavior in this area will be identified later.Next, after short stay in the strange area, the system gets back to periodic motion until  = 2.719 Hz, through series of inverse period-doubling bifurcation from period-8 to period-2.Nevertheless, the system directly enters the chaotic state after period-2 motion, as is shown in Figure 5(b), and the max LE turns positive.When  ∈ (2.78-4) Hz, periodic motion, period-doubling bifurcation, and chaotic motions appear alternately and the system suffers the saddle-node bifurcation.Beyond the threshold for the onset of chaotic motion, there are some "periodic windows, " which could be the feature of the transient chaos.The max LEs before and after the saddle-node bifurcation point are opposite in sign.
Then, for larger frequency over 3.73 Hz, there are no longer "periodic windows" presented up to the cutoff value of  in simulation.In Figure 5(c), we can observe that the system escapes from chaos to periodic motions through reverse period-doubling bifurcation.
In order to give a rather clear presentation of the dynamical behavior, we depict phase portraits (with lines) and Poincare maps (with points) [21].Figures 6(a The displacement   and velocity V  are adopted.A periodic dynamic response is exhibited, including period-1, period-2, period-4, and period-8.In Figure 6(a), the period-1 motion is presented, according to the one loop circle as well as a single point of the Poincare map.The period-1 motion is replaced by period-2, period-4, and period-8 in turn, which confirms the analysis of process of bifurcating.Figure 6(e) shows the phase plane at  = 2.471 Hz in the above-mentioned strange area.Note that the phase plane is keeping regular even consist of amount of closed curves, and the Poincare map contains limited points.In order to verify the dynamic behavior, in the above strange area, we plot the timing diagram and calculate the power spectrum density (PSD) diagram of   for  = 2.4617 Hz.It is indicated that periodic and chaotic coexisting state exists when the system falls in the strange area [22], as is shown in Figure 7. Figure 6(f) reveals the chaotic attractor, because both of the phase plane and Poincare map distribute irregularly throughout the phase space.In the chaotic area, the saddle-node bifurcation is confirmed by phase plane portraits and Poincare maps, as is shown in Figures 6(g), 6(h), and 6(i).Note, as shown in Figure 6(i), a new chaotic attractor exists as  is increased, which covers a larger area compared with Figure 6(f).It means that the vibration amplitude rapidly increases with a small increase of road excitation frequency, which is harmful to driving safety.
The road surface amplitude is also essential in vehicle dynamics analysis.Figure 8 illustrates the influence of road amplitude on the system dynamics, of the global bifurcation graph and LE spectrum under different amplitude of road excitation.The amplitude varies from 0.005 m to the max 0.1 m, containing the normal road condition, and frequency  is fixed at 5 Hz.It is observed that the displacement   and velocity V  exhibit a periodic response, including period-1, period-2, and period-4 motion, until the road amplitude is increased to   = 0.0315 m. Figure 9 expresses the phase portraits and Poincare map corresponding to   = 0.01 m, 0.02316 m, 0.02720 m, and 0.06 m, respectively.The phase plane and Poincare map is plotted in Figures 9(a), 9(b), and 9(c), which confirms the process of perioddoubling bifurcation shown in Figure 8(a).Note, for   = 0.02720 m, a limit circle occurs, which corresponds to the phenomenon of jump in the process of period-2 shown in Figure 8(a).When the amplitude of road excitation is increased larger than   = 0.0315 m, the system enters the chaotic region from Figure 8 because the max LE turns positive.In addition, just as is shown in Figure 9 The system maintains chaotic motion and vibration amplitude of suspension is getting far larger than the sustainable limit, which seriously threatens the handling safety of the vehicle.
The results of analysis above are complete descriptions of the set of parameters where chaos occurs under road profiles.The proposed analytical method is useful for estimating the suspension parameters so that chaos does not occur as desired, and thus avoiding the unexpected dangers in the running process of the vehicle.

Conclusions
This paper studies the nonlinear dynamics problems in the practical application of MRD.The 2-DoF MR vehicle suspension was set up based on the identified data of a commercial MRD.By calculating the eigenvalues of the Jacobian matrix at fixed point, the possibility of chaotic movement in the system was discovered.The theoretical analysis was then confirmed by numerical simulations.Under the single frequency harmonic excitation, the nonlinear dynamic evolution process was analyzed with frequency varying.The main conclusion would be that the loss of stability of the system appears near the midfrequency band and high amplitude of road surface, through the dynamical evolution to chaos by period-doubling bifurcations, saddlenode bifurcations, and reverse period-doubling bifurcations.These show the importance of parametric excitation in the control of vibrations in a MR suspension system.The research

Figure 3 :
Figure 3: (a) A pictorial view of the candidate MRD.(b) Comparison of tested and calculated results.

Figure 4 :
Figure 4: Frequency response diagram and LE spectrum diagram.