In this paper, the stability of vehicle concerning the slow-varying sprung mass is analyzed based on two degrees of freedom quarter-car model. A mathematical model of vehicle is established, the nonlinear vibration caused by sprung mass vibration is solved, and frequency curve is obtained. The characteristics of a stable solution and the parameters affecting the stability are analyzed. The numeric solution shows that a slow-varying sprung mass is equivalent to adding a negative damping coefficient to the suspension system, making the effective damping coefficient change from negative to positive. Such changing parameters lead to Hopf bifurcation and a shrinking limit cycle. The simulation results indicate the existence of static as well as dynamic bifurcation and the result is a change in the final stable vibration of the suspension. Even the tiny vibration of the sprung mass will lead to amplitude mutation, leading to the sprung mass instability.
1. Introduction
The stability and bifurcation of a nonlinear system are closely related. Bifurcations that include both static and dynamic bifurcations affect the stability of the system. Static bifurcation such as saddle node bifurcation occurs mainly due to the nonlinear stiffness, discussed in vibration absorbers [1, 2]. Dynamic bifurcation such as Hopf bifurcation occurs due to the parameter variations, reported in induction motor drive system [3] and circuit system [4]. They separated static bifurcation from dynamic bifurcation.
Research on static bifurcation varies. Static bifurcation caused by maneuvering load of an aircraft is analyzed in [5]. The bifurcation of a bearing system is studied in [6] and it was found that when the speed exceeds a critical value, the saddle node bifurcation occurs. The static bifurcation of an oscillator and its control is studied in [7, 8], illustrating the connection between nonlinear stiffness and saddle node bifurcation. Vibration analysis of a drill string using FEM [9] enables easy reconfiguration of the drill string for different boundary conditions. The model of these works is based on a smooth vibration state without any restriction. Compared with that, a vibro-impact system [10–13] is repeatedly reported. In these literatures, bifurcation of a piecewise-linear impact oscillator with drift is analyzed in [10]. Works in [11–13] concentrated on slip and contact motion of an oscillator excited by moving base. These researches go further to find more parameters which could result in bifurcation. However, the bifurcation is still static. The systems of all the works mentioned above can be treated as static models since their parameters remain still. But the parameters of any practical system always change with time.
Considering the case of high-speed rail CRH3 in China as the object, the critical speed leading to Hopf bifurcation between the wheel and the rail, which is critical for the safety of high-speed trains, is analyzed [14]. The stability of the rotational motion is studied in [15–17], of which the work presented in [16] analyzed an asymmetric rotation movement and found that both the Hopf bifurcation and the saddle node bifurcation will occur when the damping coefficient or the speed reaches a certain value. The kinematic stability of a symmetrical vane wheel of an aircraft is studied in [15] and it was found that any variation of the damping between the blade and the bearing, as well as the rotational speed, may lead to the generation of Hamilton Hopf bifurcation. Literatures [18–22] concern the effects of Coulomb damping on the system. Hopf bifurcation occurs as the effect of nonsmooth Coulomb friction [18]. After analyzing the vibration-impact dynamics of SDOF oscillator under the effect of Coulomb damping, it has been found out that there exist abundant dynamic phenomena, such as period-doubling routes to chaos [19, 20] and jump between multiperiodic orbits [21, 22]. A 1 : 2 resonant Hopf-Hopf bifurcation [23] and cellular shock instability are also studied concerning Hopf bifurcation [24]. Compared with the literature mentioned above, these references focused on the effect of parameter variations on the stability of the system. From these studies, it can be seen that the variations in system parameters can lead to the bifurcation. However, these studies are limited to the time delay and any changes in the characteristics of the system, such as variations in mass, stiffness, and damping coefficient, are neglected. The effect of the slowly varying parameters leading to bifurcation is rarely studied, and this kind of changes is widespread. Hence, it is worth paying attention to.
Slowly varying parameter defines the parameter which demands one-tenth time or less compared with other parameters with same self-variation in one system. Vehicle suspension is a system where the sprung mass is slowly varying. The fast-varying parameters like velocity or acceleration vary back and forth quickly. They are considered as variables in system equations. Different from that, the sprung mass is usually considered as constant. Actually, it decreases slowly, adding a slowly varying parameter to a suspension system, which changes the form of the equations. This is always ignored by researchers. For example, an adaptive control of nonlinear uncertain active suspension systems is analyzed in [25]. A vehicle seat with negative stiffness structures [26] and nonlinear asymmetrical shock absorber for vehicles [27] are presented. A controller considering actuator time delay is presented in [28, 29]. They focus only on the control of nonlinear suspension and time delay, ignoring the slowly changing parameters, especially the variations in sprung mass. Enlighted by bifurcation research and the inadequate nonlinear suspension study mentioned above, we decide to figure out the effect of varying sprung mass on the stability of the nonlinear suspension.
This paper is organized as follows. In Section 2, a mathematical model for 1/4 nonlinear suspension is established and is solved using the method of multiple scales. Section 3 analyzes the stability of the solution and possible bifurcation. Section 4 verifies the effectiveness of the proposed theoretical analysis. Section 5 concludes the paper.
2. Governing Equation for Nonlinear Suspension with Slowly Varying Sprung Mass under Forced Vibration
Leaving the turn of the vehicle out of account, the suspension model can be described in Figure 1, which includes eccentric mass me, the sprung mass ms(τ), and the unsprung mass mu, connected by a nonlinear spring ks(δs) and a shock absorber cs(ξs). Since the sprung mass is a slowly varying parameter, it is ms(τ) instead of ms(t) based on the definition of slowly varying parameter mentioned in Introduction. The wheel can be simplified as a nonlinear spring ku(δu) and a damper cu(ξu). When the automobile travels in the freeway, the construction vehicle crawls in the construction site, or the military vehicle slows down preparing for shooting, the unevenness of the road is small, which can be treated as a plane. And the vertical vibration is meecos(Ωt).
Nonlinear quarter-car model.
According to Newton’s second law, the governing equation is(1a)ddtmsτdzsdt+fcsdzsdt-dzudt+fkszs-zu=meecosΩt,(1b)mud2zsdt2-fcsdzsdt-dzudt-fkszs-zu+fcudzudt+fkuzu=0,where fcs(dzs/dt-dzu/dt)=cs(dzs/dt-dzu/dt)+ξs(dzs/dt-dzu/dt)3, fks(zs-zu)=ks(zs-zu)+δs(zs-zu)3, fcu(dzs/dt)=cudzs/dt+ξu(dzs/dt)3, fku(zu)=kuzu+δuzu3. The damping term and the spring term share the same form ax+bx3 in accordance with literatures [29, 30].
Owing to the fact that the slow time satisfies τ=εt, so we can get d/dt=d/dτdτ/dt=εd/dτ. And (1a) and (1b) can be simplified as follows:(2a)msτzs¨+εm′τzs.+cszs.-zu.+ξszs.-zu.3+kszs-zu+δszs-zu3=meecosΩt,(2b)muzu¨-cszs.-zu.-ξszs.-zu.3-kszs-zu-δszs-zu3+cuzu.+ξuzu.3+kuzu+δuzu3=0,where m′(τ)=dm(τ)/dτ,zs..=d2zs/dt2,zs.=dzs/dt,zu..=d2zu/dt2,zu.=dzu/dt.
Comparing to suspension system with constant sprung mass, this system contains an additional damping term εm′(τ)zs., and the damping coefficient εm′(τ) of this term is the function of the slowly varying sprung mass. When the sprung mass varies from big to small, this slowly varying damping coefficient is negative, resulting in self-excited vibration. Meanwhile, the natural frequency is slowly varying, leading to slowly varying amplitude and phase angle.
The sprung mass vibration makes the whole suspension in a forced vibration state. Equations (2a) and (2b) can be changed as(3a)x¨+εμx.+2ςxwxx.-y.+2ζxwxx.-y.3+wx2x-y+βxx-y3=fcosΩt,(3b)y¨-2ςxywxyx.-y.-2ζxywxyx.-y.3-wxy2x-βxyx-y3+2ςywyy.+wy2y+βyy3=0,where zs=x, zu=y, μ=ms′(τ)/ms(τ), ςx=cs/2wxms(τ), ζx=ξs/2wxms(τ), wx2=ks/ms(τ), βx=δs/ms(τ), ςxy=cs/2wxymu, ζxy=ξs/2wxymu, wxy2=ks/mu, βxy=δs/mu, ςy=cu/2wymu, wy2=ks+ku/mu, βy=δu/mu, f=mee/ms(τ).
The sprung mass varies due to the change of load. When ms≈9muks/(ks+ku), the 3 : 1 internal resonance occurs.
Because the exact solution of (3a) and (3b) cannot be found, this paper applies the method of multiple scales to solve the equation. A small perturbation parameter ε is introduced, and the scale transformation is carried out.(4a)ςx⟶εςx,ζx⟶εζx,ςxy⟶εςxy,ζxy⟶εζxy,βx⟶εβx,βxy⟶εβxy,wxy2x⟶εwxy2x,f⟶εf,(4b)ςy⟶ε2ςy,βy⟶ε2βy.
Because the stiffness of the tire ku is much bigger, the nonlinearity term δu is much smaller than that of the spring, and the equivalent damping coefficient cu is much smaller than the damper, the perturbation parameter is ε2 (4b) when rescaling stiffness and damping coefficient of the tire in the same equation.
Substituting (4a) and (4b) into (3a) and (3b) and retaining ε0 and ε1 terms yield the following equations:(5a)x¨+εμx.+2εςxwxx.-y.+2εζxwxx.-y.3+wx2x-y+εβxx-y3=εfcosΩt,(5b)y¨-2εςxywxyx.-y.-2εζxywxyx.-y.3-εwxy2x-εβxyx-y3+wy2y=0.
The term wxy2x in (3b) is much smaller in value than the term wy2y, as such this term can be considered as a small perturbation term. Because of that, it can be rescaled as εwxy2x in (5b).
According to the method of multiscale, the approximate solution can be expressed as(6a)xt,ε=x0T0,T1+εx1T0,T1+Oε2,(6b)yt,ε=y0T0,T1+εy1T0,T1+Oε2,where T0=t can be seen as faster time scale and T1=εt can be seen as slower time scale.
Under 3 : 1 internal resonance condition,(7a)Ω=wx+εσ,(7b)wy=3wx+εσ2,where σ is a detuning parameter expressing the closeness of Ω to wx and σ2 is a detuning parameter expressing the closeness of wy to 3wx.
There are differential operators like these (8a)ddt·=D0+εD1·,d2dt2·=D02+2εD0D1·,(8b)D0=∂∂T0,D1=∂∂T1.
Substituting (6a), (6b), (7a), (7b), (8a), and (8b) into (5a) and (5b) and listing equations with the same order of ε result in the following equations:(9a)ε0:D02x0+wx2x0=wy2y0,(9b)D02y0+wy2y0=0,(10a)ε1:D02x1+wx2x1=-2D0D1x0-μD0x0-2ςxwxD0x0-y0-2ζxwxD0x0-y03+wx2y1-βxx0-y03+fcosΩT0,(10b)D02y1+wy2y1=-2D0D1y0+2ςxywxyD0x0-y0+2ζxywxyD0x0-y03+wxy2x0+βxyx0-y03.
The general solution of (9a) and (9b) can be expressed as(11a)x0=GBeiwyT0+GB¯e-iwyT0+AeiwxT0+A¯e-iwxT0,(11b)y0=BeiwyT0+B¯e-iwyT0,where amplitudes A and B are functions of slower time T1. A¯ and B¯ are the complex conjugates of A and B and G=wx2/(wx2-wy2).
Substituting (11a) and (11b) into (10a) and (10b) and eliminating the secular terms result in the following equations:(12a)-2A′iwx-μAiwx-2ςxwxAiwx-2ςxwxi3A2A¯wx3+6ABB¯G-12wxwy2+3A¯2B1-Gwx2wyeiσ2T1-βx3A2A¯+6ABB¯G-12+3A¯2BG-1eiσ2T1+12feiσT1=0,(12b)-2B′iwy-2ςxywxyBiwy+2ςxywxyGBiwy+2ζxywxyi6AA¯BG-1wx2wy+3B2B¯G-13wy3-A3wx3e-iσ2T1+wxy2GB+βxy6AA¯BG-1+3B2B¯G-13+A3e-iσ2T1=0.
Change A and B into polar form:(13a)A=12aeiβ,(13b)B=12beiθ,where a,b,β, and θ are real functions of slower time T1.
Substituting (13a) and (13b) into (12a) and (12b) and separating real and imaginary parts yield the following equations:(14a)a′=-l1a-h1a3-h2ab2-h3a2bcosϕ-3γ+h4a2bsinϕ-3γ-h5sinγ,(14b)aγ′=-g1a+h6a3+h7ab2-h3a2bsinϕ-3γ-h4a2bcosϕ-3γ-h5cosγ,(14c)b′=-l2b-s1a2b-s2b3-s3a3cosϕ-3γ-s4a3sinϕ-3γ,(14d)bϕ′=-g2b+s5a2b+s6b3+s3a3sinϕ-3γ-s4a3cosϕ-3γ,where l1=μ/2+wxςx, h1=3/4ζxwx3, h2=3/2ζxwxwy2(G-1)2, h3=3/4ζxwx2wy(1-G), h4=-3/8βx/wx(G-1), h5=1/2f/wx, g1=σ, h6=3/8βx/wx, h7=3/4βx/wx(G-1)2, l2=ςxywxy(1-G), s1=-3/2ζxywxywx2(G-1), s2=-3/4ζxywxywy2(G-1)3, s3=1/4ζxywxywx3/wy, s4=-1/8βxy/wy, g2=3σ-σ2+1/2wxy2/wyG, s5=-3/4βxy/wy(G-1), s6=-3βxy/8wy(G-1)3, γ=β-σT1, ϕ=θ+σ2T1-3σT1.
3. Stability Analysis of Solutions3.1. Frequency Response and Stable Region of Solution
Since sprung mass m(τ) is function of slow time scale τ, it can be considered constant in fast time scale t. The steady-state response of system can be sought by letting a′=b′=0 and γ′=ϕ′=0 in (14a) and (14b). There are two kinds of solutions, uncoupled case a≠0, b=0 and coupled case a≠0,b≠0. As for uncoupled case, (14a) and (14b) can be simplified into(15)h52-l1a+h1a32-g1a-h6a32=0.
The number of solutions concerning (15) can be 1 or 3. When the intensity of excitation from sprung mass vibration is beyond a certain value, the number of solutions is 3, resulting in saddle node bifurcation and jump phenomenon in amplitude-frequency curve. When the intensity is within this value, only one solution can be sought. And this critical value can be solved when there is only one solution in (15) throughout the whole frequency band.(16)s5_stable=827h12+h622l1h6h1+3h6h6-3h1-l1h13.
As can be seen from (16), this critical value is related to damping coefficient and nonlinearity stiffness of the spring.
Combining (7a), (7b), (11a), and (11b) and (13a) and (13b), then we get(17a)x=Gbcos3Ωt+ϕ+acosΩt+γ+Oε,(17b)y=bcos3Ωt+ϕ+Oε.
From (17a) and (17b) we can see that, in the coupled case, there exist two vibration frequencies, that is, forced vibration frequency Ω and free vibration frequency 3Ω, due to the existence of internal resonance. The nonlinear stiffness of the spring makes the free vibration frequency exactly three times the forcing frequency. As for uncoupled case, the vehicle vibration is mainly reflected by the vibration of sprung mass, and there is only one vibration frequency in the vehicle.
When the suspension is in the uncoupled case, the steady statement of sprung mass vibration is that of autonomous system in the singularity (a¯,r¯). Linearize the system.
Linearizing (14a) and (14b) in the singularity (a¯,r¯) yields the following equation about disturbance of Δa¯ and Δr¯:(18a)dΔa¯dT1=-l1Δa¯-3h1a¯2Δa¯-h5cosγ¯Δγ¯,(18b)dΔγ¯dT1=2h6a¯Δa¯+h5a¯2cosγ¯Δa¯+h5a¯sinγ¯Δγ¯.
Use a′=γ′=0 to eliminate γ¯ in the former equation and secular equation can be got:(19)det-l1-3h1a¯2-λNM-l1-h1a¯2-λ=0,where M=1/a¯-g1+3h6a¯2,N=-a¯-g1+h6a¯2.
As for 2l1+4h1a¯2>0, the condition of losing stability for steady-state solution is(21)l12+g12-4l1h1-h6g1a¯2+3h12+h62a¯4<0.The stable region is solved from (21).
3.2. Stability of Solutions Accompanied with Bifurcation
In the slow time scale τ, the excitation only changes the position of the equilibrium point of the dynamic bifurcation. It does not change the nature of the bifurcation. Therefore, the excitation is out of consideration when analyzing the dynamic bifurcation of the suspension. We can make h5=0. As for uncoupled case, (14a) and (14b) can be transformed into(22a)a′=-l1a-h1a3,(22b)γ′=-g1+h6a2.
Focus on (22a) and let(23)za=def-l1a-h1a3=0.
Then we get two solution curves intersecting at (0, 0). That is,(24a)a=0,(24b)a=-l1h1.
They are point and circle in polar coordinates, which correspond to the equilibrium point and limit cycle of the two-dimensional systems (22a) and (22b). Focus on its stability and define(25)z′a=dzda=-l1-3h1a2.
For trivial solution a=0:z′(0)=-l1. At the time l1>0, a=0 is asymptotic stable. At the time l1<0, (0, 0) is not stable.
For nontrivial solution a=-l1/h1:z′(-l1/h1)=2l1. At the time l1<0, the system has asymptotically stable solution. And at the time l1>0, the system has unstable solutions. The former corresponds to a supercritical pitchfork bifurcation, while the latter corresponds to subcritical pitchfork bifurcation.
For the nonlinear suspension system, the sprung mass m(τ) reduced slowly from large to small. In the beginning when l1<0, the system has nontrivial solution a=-l1/h1. As the sprung mass continues decreasing, the limit cycle is reduced until the trivial solution a=0 emerges. Then, the sprung mass continues to decrease, l1>0, and the system solution is asymptotically stable. In the range l1<0, the initial amplitude inside the limit cycle increases and that outside the limit cycle reduces, all converging in the stable limit cycle in the end. This process leads to change in the static bifurcation nature of some initial values after they go through the dynamic bifurcation, destroying the system stability.
4. Numerical Validation4.1. Frequency Response and Stable Region of Solution
Let mu=45kg,ks=2000N/m,δs=660N/m,ku=16000N/m. When the damper is out of work, the damping decreases sharply to cs=20Ns/m,ξs=0.3Ns/m. At the same time, the slowly varying sprung mass satisfies(26)msτ=80e-5/9ττ≤0.5245τ>0.52.When τ>0.52, the sprung mass remains constant, and wy=20rad/s,wx=6.6667rad/s, which means wy≈3wx. These parameters are simplified from real trucks like SY5500THB, SY5411THB, and XZJ5470THB53.
According to (16), the maximum excitation acceleration in which the bifurcation does not occur in the steady-state is fstable=2wxs5_stable=2.8677m/s2. When the excitation acceleration increases, the amplitude of the sprung mass vibration will show bifurcation caused by the saddle node, that is, jumping phenomenon in the amplitude-frequency curve.
Figure 2 shows that when the excitation acceleration is less than fstable, the amplitude-frequency curve is single-valued and remains the same in the forward frequency sweeping or backward frequency sweeping. When the excitation acceleration is much bigger than fstable, jump phenomenon occurs. When the frequency is swept up, the amplitude is getting higher to the top and then jumps down to a much lower level. As the frequency continues to increase, the amplitude declines. When the frequency is swept down, the opposite jump-up phenomenon will occur. In other words, there exists jump phenomenon caused by bifurcation in amplitude-frequency curve. The curve is different in the forward frequency sweeping or backward frequency sweeping. With the increase of excitation acceleration, the vibration amplitude of the sprung mass under all frequency bands increases, and the unstable frequency band increases.
The effect of varying acceleration f on frequency response of vibration a. f=2m/s2 (black squares) and f=6m/s2 (red triangles), f=10m/s2 (blue crosses).
Figures 3 and 4 are frequency response curves under the conditions of δs=660N/m and ξs=0.3Ns/m, respectively, and fstable=2.8677m/s2 in common. As can be seen in Figure 3, when the nonlinear damping coefficient of the damper is bigger than 0.3Ns/m, under which the solution is stable, the steady-state solution of the sprung mass vibration does not bifurcate. On the contrary, bifurcation occurs in the steady-state solution of the sprung mass vibration. And the smaller the damping coefficient is, the wider the unstable bands become. Different from Figure 2, only the amplitude of which the frequency is close to the resonant frequency will get larger as the nonlinear damping coefficient decreases. Those beyond the resonant frequency change later. This indicates that the decreasing nonlinear damping coefficient has effect on the amplitude of which the frequency is close to the resonant frequency, enlarging both the amplitude and the unstable frequency bands. In the rest of the frequency band, the amplitude has nothing to do with the damping coefficient. Figure 4 shows that the nonlinearity of the spring does not increase the vibration amplitude of the sprung mass but causes the resonance point to shift to right and generates unstable bands. Concluding from Figures 2, 3, and 4, the increasing acceleration of the excitation, the decreasing nonlinear damping coefficient of the damper, and the increasing nonlinearity of the spring result in wider and wider unstable frequency bands.
The effect of varying nonlinear damping coefficient ξs on frequency response of vibration a. ξs=1Ns/m(black squares), ξs=0.5Ns/m (red triangles), and ξs=0.05Ns/m (blue crosses).
The effect of varying nonlinear stiffness δs on frequency response of vibration a. δs=60N/m(black squares), δs=460N/m (red triangles), and δs=860N/m (blue crosses).
From Figure 2, it can be found that when the amplitude-frequency curve shows an unstable frequency band, the excitation acceleration and force frequency will show multiple values. This is reflected more clearly in Figure 5. Let ξs=0.3Ns/m and δs=660N/m; then we get Figure 5. In a smaller force frequency like Ω=6.67rad/s, the amplitude of the sprung mass vibration is a single-valued curve; when Ω=6.77rad/s and Ω=6.87rad/s, the curve changes to multivalue. Similar to Figure 2, in the multivalue curve, there is a jumping phenomenon in the process of increasing or decreasing the excitation acceleration and the acceleration response curve is different. The force frequency gets bigger; the unstable interval of the excitation acceleration gets wider and gradually shifts to right. Different from Figure 2, with small excitation acceleration, the amplitude decreases as the force frequency increases; with big excitation acceleration, the opposite will occur.
The effect of varying force frequency Ω on force acceleration response of vibration a. Ω=6.67rad/s(black squares), Ω=6.77rad/s (red triangles), and Ω=6.87rad/s (blue crosses).
Figures 6 and 7 are response curves under the conditions of δs=660N/m and ξs=0.3Ns/m, respectively, and fstable=6m/s2 in common. As can be seen in Figure 6, in a small force frequency like Ω=6.67rad/s, the amplitude of the sprung mass vibration is a single-valued curve, which decreases as the damping coefficient increases. When Ω=6.77rad/s and Ω=6.87rad/s, the curve changes to multivalue. In the process of increasing in the nonlinear damping coefficient, the unstable damping coefficient interval under small force frequency gradually moves to right and decreases, and curves converge in the end. Figure 7 shows that when the nonlinear stiffness of the spring increases, the bigger the force frequency is, the larger the vibration amplitude is, and the unstable nonlinear stiffness interval increases and shifts to right. But the maximum amplitude of the sprung mass vibration in the whole nonlinear stiffness range is equal.
The effect of varying force frequency Ω on nonlinear damping coefficient response of vibration a. Ω=6.67rad/s(black squares), Ω=6.77rad/s (red triangles), and Ω=6.87rad/s (blue crosses).
The effect of varying force frequency Ω on nonlinear stiffness response of vibration a. Ω=6.69rad/s(black squares), Ω=6.73rad/s (red triangles), and Ω=6.77rad/s (blue crosses).
Let the nonlinear stiffness δs=660N/m and the nonlinear damping coefficient ξs=0.3Ns/m; Figures 8 and 9 are obtained. As can be seen in Figure 8, with the increase of the damping coefficient, the unstable region reduces along the direction of the saddle node. Similar to Figure 8, as the nonlinear stiffness increases, the whole unstable region in Figure 10 enlarges and gradually moves down.
The effect of varying nonlinear damping coefficient ξs on unstable area of vibration a. ξs=1Ns/m(black squares), ξs=0.5Ns/m (red triangles), and ξs=0.05Ns/m (blue crosses).
The effect of varying nonlinear stiffness δs on unstable area of vibration a. δs=260N/m(black squares), δs=460N/m (red triangles), and δs=860N/m (blue crosses).
The effect of varying force frequency Ω on unstable area of vibration a. Ω=6.77rad/s (black squares), Ω=6.87rad/s (red triangles), and Ω=6.97rad/s (blue crosses).
Let the nonlinear damping coefficient ξs=0.3Ns/m and the nonlinear stiffness δs=660N/m; Figures 10 and 11 are obtained. Both Figures 10 and 11 indicate that, in the whole range of nonlinear damping coefficient, the force frequency increases, and the unstable region increases and moves to upper right.
The effect of varying force frequency Ω on unstable area of vibration a. Ω=6.72rad/s (black squares), Ω=6.77rad/s (red triangles), and Ω=6.82rad/s (blue crosses).
4.2. Stability of Solutions Accompanied with Bifurcation
When τ>0.52, the sprung mass remains constant and the saddle node bifurcation occurs.
The critical lower limit frequency is wdown; the critical upper limit frequency is wupper. When the force frequency Ω gets close to wdown, all vibrations under different initial conditions are attracted to the stable focus P1. Once the force frequency Ω satisfies Ω>wdown, another stable focus P3 is generated, which is shown in Figures 12(a) and 12(b). It can be seen in Figure 12(c), when the force frequency Ω satisfies wdown<Ω<wupper, there exist two stable focuses P1,P3 and a saddle node P2. Only certain initial leads to trajectory which can get to the saddle node P2. With a small disturbance, the system state will be attracted to focus P1 or P3. The dashed line in Figure 12(c) is the dividing line of the two regions. The state outside the dashed line is attracted to the stable focus P3 and those inside the dashed line is attracted to the stable focus P1. As the force frequency continues to increase, at the time Ω>wupper, the focus P1 disappears, and all states are attracted to the focus P3. The region covered by red dash dotted line varies from large to small and ultimately disappears. These phenomena indicate that the upper and the lower two solutions are asymptotically stable in the amplitude-frequency curve, and the intermediate solution is unstable. Because of the fact that only asymptotically stable motion can be achieved in actual physical world, the jumping phenomenon can be watched.
Evolution of the phase plane with the force frequency: (a) 6.81, (b) 6.82, (c) 6.85, and (d) 6.88 rad/s.
When τ≤0.52, the sprung mass changes slowly, and dynamic Hopf bifurcation occurs. The limit cycle shrinks as the sprung mass reduces. The initial amplitude inside the limit cycle increases and that outside the limit cycle reduces, all converging in the stable limit cycle in the end. At that time τ=0.52, due to the impact of Hopf bifurcation, different initials converge in the limit cycle a=0.7, as shown in Figure 13.
Phase portrait and time-domain vibration amplitude of Hopf bifurcation.
This process changes the initial value. Since the nonlinear vibration is closely related to initial state, it will change the nature of the vibration.
In Figure 14, blue line represents vibration beginning with the initial value (2, 0.7) while the red line represents vibration beginning with the initial value (0.7, 0.7). The initial (0.7, 0.7) can be obtained from initial (2, 0.7) after being affected by Hopf bifurcation. As can be seen in Figure 14(b), when the excitation frequency satisfies wdown<Ω<wupper, the final stable state is different from these two initial values. At the time τ=0.52, the effect of the Hopf bifurcation changes the initial value of the blue line to the initial value of the red line. Once the sprung mass is constant, the Hopf bifurcation disappears and the new initial value is affected by the saddle node bifurcation. These indicate that the stability of the system is extremely fragile.
Different solution induced by Hopf bifurcation.
5. Conclusion
In this paper, the stability of vehicle concerning the slow-varying sprung mass is analyzed based on two-degree-of-freedom quarter-car model. The vibration amplitude solution is solved using multiscale method, followed by numerical validation. After analyzing amplitude-frequency curve, phase plane of steady-state response and phase portrait, and time-domain vibration amplitude of Hopf bifurcation, the main points of our work can be concluded as follows:
(a) There exist rich mechanical phenomena, which include internal resonance, jumping phenomenon, and Hopf bifurcation.
When the unsprung mass mu, sprung mass ms, the linear stiffness of tire ku, and the linear stiffness of the suspension spring ks satisfy the condition ms≈9muks/(ks+ku), the vibration frequency ratio of unsprung mass and sprung mass is 3 : 1, and internal resonance will occur. Under the condition of internal resonance, when the sprung mass and the unsprung mass are not coupled, the vibration amplitude of the sprung mass creates a jumping phenomenon. Meanwhile, since the sprung mass is slowly varying, dynamic Hopf bifurcation occurs.
(b) The essential of the effect of slowly varying sprung mass is changing the damping coefficient.
The numeric solution shows that a slow-varying sprung mass is equivalent to adding a negative damping coefficient to the suspension system, making the effective damping coefficient change from negative to positive. Such changing parameters lead to Hopf bifurcation and a shrinking limit cycle.
(c) The stability of nonlinear suspension is closely related to damping coefficient, nonlinear stiffness, and a slowly varying sprung mass.
Because of the slowly varying process, vibration with different initial states converges in the stable limit cycle in the end. They cannot result in final vibration with different states. Since the nonlinear vibrations are closely related to the initial states, Hopf bifurcation can change the final stable state of the vibration.
All of the above observations indicate the existence of static as well as dynamic bifurcation and the result is a change in the final stable vibration of the suspension. Even the tiny vibration of the sprung mass will lead to amplitude mutation, leading to the sprung mass instability.
Conflicts of Interest
The authors declare that no conflicts of interest exist regarding the publication of this paper.
FebboM.MacHadoS. P.Nonlinear dynamic vibration absorbers with a saturation201333261465148310.1016/j.jsv.2012.11.0252-s2.0-84871712963DjemalF.ChaariF.DionJ.-L.RenaudF.TawfiqI.HaddarM.Performance of a non linear dynamic vibration absorbers201531334535310.1017/jmech.2014.762-s2.0-84930765028DeivasundariP. S.UmaG.SanthiR.Experimental verification of Hopf bifurcation in pulse-width modulated inverter fed cage induction motor drive system20147234034910.1049/iet-pel.2013.00632-s2.0-84894271763HanX. J.XiaF. B.JiP.BiQ. S.KurthsJ.Hopf-bifurcation-delay-induced bursting patterns in a modified circuit system201636517527HouL.ChenY. S.LuZ. Y.LiZ. G.Bifurcation analysis for 2:1 and 3:1 super-harmonic resonances of an aircraft cracked rotor system due to maneuver load2015811-253154710.1007/s11071-015-2009-12-s2.0-84930758970SaeedN. A.EissaM.El-GaniniW. A.Nonlinear oscillations of rotor active magnetic bearings system2013741-212010.1007/s11071-013-0967-8MR31051712-s2.0-84883812480AcarM. A.YilmazC.Design of an adaptive–passive dynamic vibration absorber composed of a string–mass system equipped with negative stiffness tension adjusting mechanism2013332223124510.1016/j.jsv.2012.09.0072-s2.0-84867671464ElnaggarA. M.KhalilK. M.Control of the nonlinear oscillator bifurcation under superharmonic resonance2013543443GhasemlooniaA.RideoutD. G.ButtS. D.Vibration analysis of a drillstring in vibration-assisted rotary drilling: finite element modeling with analytical validation201313531803290210.1115/1.40233332-s2.0-84878523770ChávezJ. P.PavlovskaiaE.WiercigrochM.Bifurcation analysis of a piecewise-linear impact oscillator with drift2014771-221322710.1007/s11071-014-1285-5MR32151412-s2.0-84902260737AndreausU.NisticòN.An analytical-numerical model for contact-impact problems: theory and implementation in a two-dimensional distinct element algorithm19983298110AndreausU.CasiniP.Friction oscillator excited by moving base and colliding with a rigid or deformable obstacle200237111713310.1016/S0020-7462(00)00101-3ZBL1117.743032-s2.0-0036028167AndreausU.CasiniP.Dynamics of friction oscillators excited by a moving base and/or driving force2001245468569910.1006/jsvi.2000.35552-s2.0-0035939863ZhangT. T.DaiH. Y.Bifurcation analysis of high-speed railway wheel-set20168331511152810.1007/s11071-015-2425-2MR34494892-s2.0-84955213370NajafiA.GhazaviM.-R.JafariA.-A.Stability and Hamiltonian HOPf bifurcation for a nonlinear symmetric bladed rotor20147821049106410.1007/s11071-014-1495-xMR32676872-s2.0-84911007981ShahgholiM.KhademS. E.Hopf bifurcation analysis of asymmetrical rotating shafts20147741141115510.1007/s11071-014-1367-4MR32474662-s2.0-84897359583HuangJ.LuoA. C. J.Analytical periodic motions and bifurcations in a nonlinear rotor system20142342545910.1007/s40435-013-0055-42-s2.0-84975229472HetzlerH.On the effect of nonsmooth Coulomb friction on Hopf bifurcations in a 1-DoF oscillator with self-excitation due to negative damping2012691-260161410.1007/s11071-011-0290-1MR29298972-s2.0-84861745483ZhangH.ZhangN.MinF.YanW.WangE.Bifurcations and chaos of a vibration isolation system with magneto-rheological damper20166303531010.1063/1.49445562-s2.0-84962268407AndreausU.CasiniP.Dynamics of SDOF oscillators with hysteretic motion-limiting stop200022214516410.1023/A:1008354220584AndreausU.CasiniP.Forced motion of friction oscillators limited by a rigid or deformable obstacle200129217719810.1081/SME-1001044792-s2.0-0035327970AndreausU.De AngelisM.Nonlinear dynamic response of a base-excited SDOF oscillator with double-side unilateral constraints20168431447146710.1007/s11071-015-2581-42-s2.0-84953402455RevelG.AlonsoD. M.MoiolaJ. L.Numerical semi-global analysis of a 1:2 resonant Hopf-Hopf bifurcation20132471405310.1016/j.physd.2012.12.007MR30199412-s2.0-84873132677YaoJ. H.O(2)-Hopf bifurcation for a model of cellular shock instability2014269637510.1016/j.physd.2013.11.005MR31513052-s2.0-84890536146HuangY.NaJ.WuX.LiuX.GuoY.Adaptive control of nonlinear uncertain active suspension systems with prescribed performance20155414515510.1016/j.isatra.2014.05.0252-s2.0-84923702691DanhL. T.AhnK. K.Active pneumatic vibration isolation system using negative stiffness structures for a vehicle seat201433351245126810.1016/j.jsv.2013.10.0272-s2.0-84890554873SilveiraM.PontesB. R.Jr.BalthazarJ. M.Use of nonlinear asymmetrical shock absorber to improve comfort on passenger vehicles201433372114212910.1016/j.jsv.2013.12.0012-s2.0-84892621422KongY.ZhaoD.YangB.HanC.HanK.Robust non-fragile H∞/L2-L∞ control of uncertain linear system with time-delay and application to vehicle active suspension201525132122214110.1002/rnc.3196ZhuY.ZhuS.Nonlinear time-delay suspension adaptive neural network active control20142014676587110.1155/2014/765871MR32562582-s2.0-84934924599YazidI. I. M.MazlanS. A.KikuchiT.ZamzuriH.ImaduddinF.Design of magnetorheological damper with a combination of shear and squeeze modes201454879510.1016/j.matdes.2013.07.0902-s2.0-84883521016