Shock and Vibration
Volume 2018 (2018), Article ID 1617650, 10 pages
https://doi.org/10.1155/2018/1617650

Research Article

Bifurcation Analysis of a Rigid Impact Oscillator with Bilinear Damping

Liping Zhang,¹ Haibo Jiang,¹ and Yang Liu²

¹School of Mathematics and Statistics, Yancheng Teachers University, Yancheng 224002, China
²College of Engineering, Mathematics and Physical Sciences, University of Exeter, Rennes Drive, Exeter, EX4 4RN, UK

Correspondence should be addressed to Haibo Jiang; yctcjhb@126.com

Received 24 September 2017; Revised 28 January 2018; Accepted 11 February 2018; Published 12 March 2018

Academic Editor: Hassan Haddadpour

Copyright © 2018 Liping Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper studies a rigid impact oscillator with bilinear damping developed as the mechanical model of an impulsive switched system. The stability and the bifurcation of periodic orbits in the impact oscillator are determined by using the mapping methods. One-parameter bifurcation analyses under variation of forcing frequency and amplitude of external excitation are carried out. Coexisting attractors and various types of bifurcations, such as grazing, period-doubling, and saddle-node, are observed, which show the complex phenomena inhered in this impact oscillator.

1. Introduction

The rigid impact, known as the impulsive reactions whenever rigid bodies collide, widely exists in many engineering applications, such as rotating machinery, car suspension systems, and cutting processes. In general, such impact can be studied using the rigid impact oscillator, which is modelled using the coefficient of restitution rule assuming the instantaneous reversal of velocity for the collision body. The rigid impact oscillator is a nonsmooth dynamical system which can exhibit complex dynamical behavior, so the stability and the bifurcation of the rigid impact oscillator have received great attention; see, for example, [1–18]. In [1], Shaw and Holmes observed the chaotic and long period motions in a class of periodically forced linear oscillators with impacts. In [2], Whiston analyzed the vibroimpact response of a one-dimensional undamped linear oscillator preloaded against a stop by a constant force and subjected to harmonic excitation. Nordmark [3] studied the singularities of grazing impacts in a single-degree-of-freedom periodically forced oscillator subjected to a rigid amplitude constraint using analytical methods. In [4], Ivanov developed the linear theory for analyzing the stability and the bifurcation of a class of rigid impact oscillators. Foale [5, 6] classified the types of grazing bifurcation in a class of rigid impact oscillators and presented analytical results to show how the type of grazing bifurcation changed with control parameter. In [7], Luo and Xie observed Hopf bifurcations in a vibroimpact system at two strong resonance cases. Luo and Chen [8] investigated the mapping dynamics of periodic motions for a class of rigid impact oscillators. Dankowicz and Zhao [9, 10] studied different bifurcation scenarios associated with the switching between impacting and nonimpacting motions of impact actuators near grazing. In [11], Thota and Dankowicz studied continuous and discontinuous grazing bifurcations of impacting oscillators by using the discontinuity mapping. In particular, the concept of discontinuity geometry was introduced in [12–14] to study the grazing bifurcation in rigid impact oscillators. Later on, Mason and Piiroinen [15] focused on the intricate relationship between smooth and nonsmooth phenomena in an impacting system. Jiang et al. [16] studied the antiphase synchronization and the symmetry-breaking bifurcation of two mechanical oscillators coupled by rigid impacts. In [17], Jiang et al. investigated the differences between the grazing-induced bifurcations in impact oscillators with one-sided elastic and rigid constraints by a path-following method.

In this paper, we will study a rigid impact oscillator with bilinear damping through one-parameter bifurcation analysis. In general, linear damping is always considered in the physical models of mechanical systems. However, in some engineering practices, piecewise linear or viscous dampers are widely designed and used owing to their adaptability in the changing environment. For example, shock absorber with bilinear characteristics is often used for improving the driving safety and traveling comfort of vehicles [19]. Therefor, piecewise linear or viscous dampers have received considerable attention from researchers and engineers, for example, [20–29]. In [19], the nonlinear stationary oscillations of a quarter-vehicle model with two degrees of freedom subjected to a vertical road excitation and bilinear damping were studied. In [20], Natsiavas presented an appropriate stability analysis for periodic solutions of harmonically excited piecewise linear systems with dry friction and damping coefficients depending on the velocity direction. Natsiavas and Verros [21] studied the dynamics of oscillators with strongly nonlinear asymmetric damping. In [22], Verros et al. studied the dynamic response of a controlled single-degree-of-freedom quarter-car model subjected to road excitation, and the control strategy applied to the system was based on the selection of two values of the damping ratio. In [23], the effect of asymmetric damping on the sprung mass position was investigated by analyzing a two-degree-of-freedom quarter-vehicle model. In [24], Verros et al. investigated the car models involving passive damping with constant and dual-rate characteristics. In [25], the model of a single degree of system subjected to harmonic base excitation was used to study the performance of vibration isolation using four semiactive dampers, and five control algorithms were studied based on the skyhook control, the balance control, and the adaptive damping control. Papalukopoulos and Natsiavas [26] investigated biodynamic response of the lumped parameter passenger-seat models coupled with the simplified models of ground vehicles, and dual-rate suspension damping was used to model the damper between the vehicles. In [27], Luo and Rajendran investigated the periodic motions and the stability in a semiactive suspension system with the magnetorheological damping which varies with relative velocity. In addition, Waters et al. [28] considered the effect of dual-rate suspension damping on vehicle response to transient road inputs. A damping control law was proposed in [29] to adjust the damping coefficient for suppressing the chaotic behavior of an impact oscillator.

Recently, there has been an increasing interest in the analysis and synthesis of impulsive switched systems with the presence of two nonsmooth effects, impulsive disturbances and switching, because of their significance in theory and applications in physics, biology, engineering, and information science [30–37]. In [31], Xie and Wang introduced the impulsive behavior into switched linear systems and studied their controllability and observability. Guan et al. [32] studied the exponential and asymptotical stability for a class of hybrid impulsive and switching nonlinear systems and developed a new hybrid impulsive and switching control strategy. In [33], Zhong and Wu studied the robust exponential stability for a class of switched systems with impulsive effect and structural perturbations. In [34], control of a class of discrete impulsive switched systems with disturbance was studied by using the Lyapunov direct method. In [35–37], the stability analysis of impulsive switched systems with time delays was studied. In this paper, we will develop a mechanical model for impulsive switched systems and carry out bifurcation analysis for this system [38]. Firstly, the physical model of a rigid impact oscillator with bilinear damping is built, and the dual-rate damper model is used to describe the bilinear damping. Secondly, in order to describe the segments of periodic orbits of the impact oscillator, we divide the phase space into four subspaces and obtain four discontinuity boundaries. Then, five local maps are defined between these discontinuity boundaries, and periodic orbits can be described by the composite of these local maps. Thus, the stability and the bifurcation of the impact oscillator can be determined by the eigenvalues of the Jacobian matrix of these global maps. Thirdly, one-parameter bifurcation analysis is carried out to show the influence of forcing frequency and amplitude on the dynamics of the system. Finally, numerical simulations show that the system exhibits complex phenomena, including periodic and chaotic orbits.

2. Physical Model

The physical model of the rigid impact oscillator with one-sided constraint is shown in Figure 1(a). Displacement of the mass is represented by . The mass is connected to the supporting base by a linear spring with stiffness and a dual-rate damper with the damping coefficient expressed as

where denotes the differentiation of with respect to the time . The relationship between the damping force and the velocity for the dual-rate damper model is presented in Figure 1(b). As shown in the figure, the damping coefficients in compression and extension are different depending on the sign of the velocity. Thus the dual-rate damper is easy to be designed in engineering and it can also achieve a better compromise performance for the vehicle between ride, road-holding, handling, and control, so it is used as a typical model in automotive engineering applications. However, to adapt complex engineering circumstances, we can consider the other types of dampers, such as three-linear, four-linear, and viscous. The dynamics of impact oscillators with these dampers may be complex, which needs further studies. External harmonic excitation is applied on the mass with amplitude and frequency . The gap between the mass and the stop is denoted by , and the stop is assumed to be rigid and described by the coefficient of restitution .

(a)

(b)

Figure 1: (a) Physical model and (b) the dual-rate damper model.

The governing equations of the impact oscillator can be written as

where and represent the velocities of the impact oscillator before and after impact, respectively. For simplicity, we introduce the following nondimensional variables and parameters, , , , , , , , and , where is the nondimensional time, is the nondimensionalized displacement, is the nondimensional gap, is the nondimensional forcing amplitude, and is an arbitrary reference distance. After the nondimensionalization, the equations of motion of the impact oscillator taking into account the dual-rate damper are given as

where denotes the differentiation of with respect to the nondimensional time and is the Heaviside step function.

3. Local Maps and Stability Analysis

The impact oscillator (3) is a state-dependent impulsive switched (hybrid) dynamical system which can be expressed as

where , , , , , , , , , , .

Given an initial value , , the motion of the impact oscillator is governed by the equation until , , and the state jumps according to . If , , the motion switches to continuous motion governed by the equation , and then, the procedure continues. Each segment can be denoted by the signature . Figure 2(a) shows the representative trajectory of nonimpacting period-1 orbit by using the segments and . Figure 2(b) presents the representative trajectory of period-1 orbit with one impact per period using the segments and .

(a)

(b)

(c)

Figure 2: Representative trajectories and segments of (a) the nonimpacting period-1 orbit, (b) the period-1 orbit with one impact per period, and (c) the local maps.

In the following, the mapping methods proposed in [8, 17, 27, 39–41] will be used to study the stability and the bifurcation of the impact oscillator (3). In order to describe the segments of the impact oscillator, we divide the phase space into four subspaces , , , and , where , , , and . Switches between these regimes occur on the discontinuity boundaries, , , , and , where , , , and . Here, we define five local maps as follows: : , : , : , : , and : , which are shown in Figure 2(c). These local maps project the point located on one of the discontinuity boundaries onto another. The local map can be represented by : , where and the trajectory is in the regime . The local map can be represented by : , where and the trajectory is in the regime . The local map can be represented by : , where , , and the trajectory is in the regime . Due to impacts at the boundary, represents the impact relation; that is, : , where , , . The local map can be represented by : , where , , and the trajectory is in the regime . The governing equations for the local maps , , , and can be obtained from the solutions of the linear ordinary differential equations in the corresponding regimes (see the Appendix). For example, the global map of the nonimpacting period-1 orbit shown in Figure 2(a) can be expressed as . Thus, the stability of this period-1 motion can be determined using the corresponding Jacobian matrix of the global map which can be computed by the chain rule . The global map of the period-1 orbit with one impact per period shown in Figure 2(b) can be expressed as , and its Jacobian matrix can be computed by . Therefore, for a periodic motion of mapping structure , where , , its stability and bifurcations can be determined by the eigenvalues of the Jacobian matrix . If two eigenvalues of the Jacobian matrix lie inside the unit circle, the periodic motion is stable. If one of them lies outside the unit circle, the periodic motion is unstable. If one of the two eigenvalues is +1 and the second one is inside the unit circle, the saddle-node bifurcation will occur. If one of the two eigenvalues is −1 and the second one is inside the unit circle, the period-doubling bifurcation may occur. In our numerical simulations, we will follow these bifurcations by using the numerical continuation method developed in [40, 41].

4. Bifurcation Analysis of the Impact Oscillator with Bilinear Damping

In this section, one-parameter bifurcation analysis of the impact oscillator with bilinear damping is carried out to investigate the influence of forcing frequency and amplitude on the dynamics of the system.

4.1. Influence of the Forcing Frequency

To investigate the influence of forcing frequency on the dynamics, we chose , , , , , and the forcing frequency was considered as a branching parameter. Bifurcation diagrams were constructed by starting from an initial frequency, plotting the solution once it converged to a periodic response (or plotting the next 500 iterates if there was no convergence after 2000 periods), and using the final conditions as the initial conditions of the next value of frequency. Bifurcation diagram constructed by the mass displacement under variation of forcing frequency is shown in Figure 3. As can be seen from Figure 3, an interesting bifurcation structure was recorded. As the forcing frequency exceeds 0.8, the nonimpacting period-1 orbit is unchanged until the orbit grazes the discontinuity boundary; that is, grazing bifurcation occurs at . Then the period-1 response changes to a stable period-3 orbit with one impact via a grazing bifurcation. As the forcing frequency is further increased, the trajectory of the period-3 orbit with one impact per period contacts the discontinuity boundary with zero velocity again; that is, the second grazing bifurcation is observed at which leads to a small regime of chaotic motions for the impact oscillator. Thereafter, a reverse period-doubling cascade occurs and the response of the impact oscillator changes from chaotic motion to a stable period-4 orbit with two impacts per period followed by a stable period-2 orbit with one impact per period , which bifurcates into a stable period-1 orbit with one impact per period eventually. With further increase of the forcing frequency, the period-1 orbit with one impact undergoes a saddle-node bifurcation at resulting in an unstable period-1 orbit with one impact per period . As the forcing frequency becomes less than 1.4, there is a nonimpacting period-1 orbit followed by a grazing bifurcation occurring at , and the state of the impact oscillator jumps to the coexisting stable period-1 orbit with one impact per period.

Figure 3: (a) Bifurcation diagrams constructed for the mass displacement under varying forcing frequency calculated at ,  ,  ,  ,  . Black and blue dots mark the mass displacement which is swept in the increasing and decreasing directions, respectively. The magenta dots show the unstable period-1 orbit with one impact. Additional windows demonstrate the trajectories of the attractors on the phase plane obtained for ,  ,  ,  , ,  ,  , and (stable impacting period-1 orbit in black, unstable impacting period-1 orbit in magenta, stable nonimpacting period-1 orbit in blue). Green solid lines indicate the location of the discontinuity boundaries and , and the red points show the Poincaré sections. (b) The corresponding amplitude frequency curve showing the maximum of the absolute values of the mass displacement as a function of the forcing frequency . Blue, green, and red dots represent one-parameter grazing (GR), period-doubling (PD), and saddle-node (SN) bifurcations, respectively. Solid and dotted lines denote stable and unstable solutions, respectively. The dashed horizontal line denotes .

Figure 3(b) shows the corresponding amplitude frequency curve of the period-1 solution. As the frequency increases, a grazing bifurcation () occurs at the critical value , leading a stable nonimpacting period-1 orbit to an unstable period-1 orbit with one impact per period. With the further increase in the frequency, a period-doubling bifurcation (PD) takes place and the unstable period-1 orbit with one impact per period regains stability. Later on, the stable period-1 impact orbit with one impact per period loses stability via a saddle-node bifurcation (SN). On the other hand, as the frequency decreases, the second grazing bifurcation (GR) occurs at the critical value . Hence, the stable nonimpacting period-1 orbit becomes unstable, and an unstable period-1 orbit with one impact per period emerges. Then, the lower branch of unstable solutions coalesces with the upper stable branch and they annihilate each other through a saddle-node bifurcation (SN). Here, the hysteresis phenomenon is observed.

Figure 4 presents the evolution of basins of attraction of the impact oscillator with bilinear damping. Here, the basin of the period-1 orbit with one impact per period (red) is given in cyan, and the nonimpacting period-1 orbit (blue) is shown in magenta. As can be seen from this series of graphs, the basin for the period-1 orbit with one impact per period gradually shrinks as the frequency increases and disappears completely at , while the nonimpacting period-1 orbit gradually grows as the frequency increases.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: Evolution of the basins of attraction for the impact oscillator with bilinear damping obtained for ,, ,,, (a) , (b) , (c) , (d) , (e) , and (f) . The period-1 orbit with one impact per period is shown in red with a cyan basin, and the nonimpacting period-1 orbit is presented in blue with a magenta basin.

4.2. Influence of the Forcing Amplitude

To investigate the influence of forcing amplitude on the dynamics, we fixed , , , , , and used the forcing amplitude as a branching parameter. The constructed bifurcation diagram is shown in Figure 5. As can be seen from the figure, an interesting bifurcation structure was recorded. As the forcing amplitude exceeds 0.3, the nonimpacting period-1 orbit is unchanged until the orbit grazes the discontinuity boundary , so the grazing bifurcation (GR) occurs at . Then, a small regime of chaotic motions can be observed which is followed by a period-2 orbit with one impact per period . As the forcing amplitude increases, the trajectory of the period-2 orbit with one impact contacts the discontinuity boundary with zero velocity, so the second grazing bifurcation (GR) occurs at leading to a period-2 orbit with two impacts per period (). Shortly after, a reverse period-doubling bifurcation (PD) occurs and the response of the impact oscillator changes to a stable period-1 orbit with one impact per period . With further increase in the forcing amplitude, a saddle-node bifurcation (SN) is observed at , and the state changes to a coexisting stable period-1 orbit with two impacts per period . Now we consider the variation of the branching parameter in a reversed direction. As the forcing amplitude becomes less than 2.4, there is a stable period-1 orbit with two impacts per period which is followed by a saddle-node bifurcation (SN) recorded at . Since then, the state of the impact oscillator jumps to the coexisting stable period-1 orbit with one impact per period . Again, the hysteresis phenomenon is observed.

Figure 5: Bifurcation diagrams constructed for the mass displacement under varying forcing amplitude (a) and (b) calculated at ,,,, and . Black and blue dots show the mass displacements swept in increasing and decreasing directions, respectively. The magenta dots show the unstable period-1 orbit with one impact. Additional windows demonstrate the trajectories of the attractors on the phase plane obtained for ,,,,,,, and , respectively. Stable period-1 orbit with two impacts per period, unstable period-1 orbit with two impacts per period, and stable period-1 orbit with one impact per period are shown in blue, magenta, and black. Green lines indicate the location of the discontinuity boundaries and , and the red points show the Poincaré sections.

Figure 6 presents the evolution of basins of attraction for the impact oscillator with bilinear damping. Here the basin of the period-1 orbit with one impact per period (red) is given in cyan and the period-1 orbit with two impacts per period (blue) is in magenta. As can be seen from these graphs, the basin of the period-1 orbit with one impact per period shrinks gradually as the frequency increases and disappears completely at . The basin for the period-1 orbit with two impacts per period grows gradually as the amplitude increases.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: Evolution of the basins of attraction for the impact oscillator with bilinear damping obtained for ,,,,, (a) , (b) , (c) , (d) , (e) , and (f) . Here the period-1 orbit with one impact per period shown in red has a cyan basin, and the period-1 orbit with two impacts per period shown in blue has a magenta basin.

5. Conclusions

The rigid impact oscillator with bilinear damping constructed as the mechanical model of an impulsive switched system was investigated through one-parameter bifurcation analysis in this paper. The stability and the bifurcation of the periodic orbits in the impact oscillator were studied by using the mapping methods. In order to describe the periodic orbit, we divided the phase space of the system into four subspaces and defined four discontinuous boundaries. Local maps were built to describe the motion of the system, and the eigenvalues of the Jacobian matrix of these global maps were used to determine the stability and bifurcation of the system.

Our one-parameter bifurcation analyses show that, as the forcing frequency varies for , the system experiences grazing, period-doubling, and saddle-node bifurcations. The coexistence of impacting period-1 and nonimpacting period-1 orbits was observed for . A series of basins of attraction of the impact oscillator were presented to show the evolution of these two coexisting attractors, where the basin for the impacting period-1 orbit gradually shrinks as the forcing frequency increases and disappears completely at .

Our one-parameter bifurcation analyses using the forcing amplitude as the branching parameter reveal that the impact oscillator encounters two grazing, one period-doubling, and two saddle-node bifurcations when the forcing amplitude varies for . When , a period-1 orbit with one impact per period and a period-1 orbit with two impacts per period coexist. The basins of attraction show that the period-1 orbit with one impact disappears completely at , and the basin for the period-1 orbit with two impacts grows gradually as the forcing amplitude increases.

Compared to the research on the rigid impact oscillators with linear damping in [1–18], the bilinear damping was considered and the mapping method was adopted to investigate the dynamics of the rigid impact oscillator with bilinear damping focusing on grazing bifurcations, which can be utilized for analyzing this type of nonsmooth dynamical systems. Compared to the research on piecewise linear dampers in [19–29], numerical methods including bifurcation diagrams, phase-portraits, and basin of attraction were carried out to show the influence of the forcing frequency and amplitude on the dynamics of the system. Furthermore, coexisting attractors and various types of bifurcations, such as grazing, period-doubling, and saddle-node, are observed, which show the complex phenomena inhered in the impact oscillator. Future works include two-parameter bifurcation analysis and investigation of the bifurcation mechanism inhered in this impact oscillator.

Appendix

The Analytical Solution of the Equation in the Smooth Regimes and

Consider the equation

with the initial conditions , . The solution is

and its directive is

where , , , , .

Disclosure

This paper was presented at the International Conference on Engineering Vibration (Sofia, Bulgaria, September 4–7, 2017).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (Grants nos. 11402224, 11672257, and 11202180), the Natural Science Foundation of Jiangsu Province of China (Grants nos. BK20161314 and BK20151295), the Jiangsu Overseas Research and Training Program for University Prominent Young and Middle-Aged Teachers and Presidents, the Excellent Scientific and Technological Innovation Team of Jiangsu University, and Jiangsu Key Laboratory for Big Data of Psychology and Cognitive Science. The authors also thank Professor Marian Wiercigroch for helpful comments.

References

S. W. Shaw and P. Holmes, “Periodically forced linear oscillator with impacts: chaos and long-period motions,” Physical Review Letters, vol. 51, no. 8, pp. 623–626, 1983.

G. S. Whiston, “The vibro-impact response of a harmonically excited and preloaded one-dimensional linear oscillator,” Journal of Sound and Vibration, vol. 115, no. 2, pp. 303–319, 1987.

A. B. Nordmark, “Non-periodic motion caused by grazing incidence in an impact oscillator,” Journal of Sound and Vibration, vol. 145, no. 2, pp. 279–297, 1991.

A. P. Ivanov, “Impact oscillations: linear theory of stability and bifurcations,” Journal of Sound and Vibration, vol. 178, no. 3, pp. 361–378, 1994.

S. Foale, “Analytical Determination of Bifurcations in an Impact Oscillator,” Philosophical Transactions of the Royal Society A: Mathematical, Physical & Engineering Sciences, vol. 347, no. 1683, pp. 353–364, 1994.

S. Foale and S. R. Bishop, “Bifurcations in impact oscillations,” Nonlinear Dynamics, vol. 6, no. 3, pp. 285–299, 1994.

G. Luo and J. Xie, “Bifurcations and chaos in a system with impacts,” Physica D: Nonlinear Phenomena, vol. 148, no. 3-4, pp. 183–200, 2001.

A. C. J. Luo and L. Chen, “Periodic motions and grazing in a harmonically forced, piecewise, linear oscillator with impacts,” Chaos, Solitons & Fractals, vol. 24, no. 2, pp. 567–578, 2005.

H. Dankowicz and X. Zhao, “Local analysis of co-dimension-one and co-dimension-two grazing bifurcations in impact microactuators,” Physica D: Nonlinear Phenomena, vol. 202, no. 3-4, pp. 238–257, 2005.

X. Zhao and H. Dankowicz, “Unfolding degenerate grazing dynamics in impact actuators,” Nonlinearity, vol. 19, no. 2, pp. 399–418, 2006.

P. Thota and H. Dankowicz, “Continuous and discontinuous grazing bifurcations in impacting oscillators,” Physica D: Nonlinear Phenomena, vol. 214, no. 2, pp. 187–197, 2006.

D. R. Chillingworth, “Discontinuity geometry for an impact oscillator,” Dynamical Systems. An International Journal, vol. 17, no. 4, pp. 389–420, 2002.

D. R. Chillingworth, “Dynamics of an impact oscillator near a degenerate graze,” Nonlinearity, vol. 23, no. 11, pp. 2723–2748, 2010.

N. Humphries and P. T. Piiroinen, “A discontinuity-geometry view of the relationship between saddle-node and grazing bifurcations,” Physica D: Nonlinear Phenomena, vol. 241, no. 22, pp. 1911–1918, 2012.

J. F. Mason and P. T. Piiroinen, “Saddle-point solutions and grazing bifurcations in an impacting system,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 22, no. 1, 013106, 6 pages, 2012.

H. Jiang, Y. Liu, L. Zhang, and J. Yu, “Anti-phase synchronization and symmetry-breaking bifurcation of impulsively coupled oscillators,” Communications in Nonlinear Science and Numerical Simulation, vol. 39, pp. 199–208, 2016.

H. Jiang, A. S. E. Chong, Y. Ueda, and M. Wiercigroch, “Grazing-induced bifurcations in impact oscillators with elastic and rigid constraints,” International Journal of Mechanical Sciences, vol. 127, pp. 204–214, 2017.

M. di Bernardo, C. J. Budd, and A. R. Champneys, Piecewise-Smooth Dynamical Systems: Theory and Applications, vol. 163 of Applied Mathematical Sciences, Springer, London, UK, 2008.

Z. Szabó and A. Lukács, “Numerical stability analysis of a forced two-D.O.F. oscillator with bilinear damping,” Journal of Computational and Nonlinear Dynamics, vol. 2, no. 3, pp. 211–217, 2007.

S. Natsiavas, “Stability of piecewise linear oscillators with viscous and dry friction damping,” Journal of Sound and Vibration, vol. 217, no. 3, pp. 507–522, 1998.

S. Natsiavas and G. Verros, “Dynamics of oscillators with strongly nonlinear asymmetric damping,” Nonlinear Dynamics, vol. 20, no. 3, pp. 221–246, 1999.

G. Verros, S. Natsiavas, and G. Stepan, “Control and dynamics of quarter-car models with dual-rate damping,” Journal of Vibration and Control, vol. 6, no. 7, pp. 1045–1063, 2000.

C. Rajalingham and S. Rakheja, “Influence of suspension damper asymmetry on vehicle vibration response to ground excitation,” Journal of Sound and Vibration, vol. 266, no. 5, pp. 1117–1129, 2003.

G. Verros, S. Natsiavas, and C. Papadimitriou, “Design optimization of quarter-car models with passive and semi-active suspensions under random road excitation,” Journal of Vibration and Control, vol. 11, no. 5, pp. 581–606, 2005.

Y. Liu, T. P. Waters, and M. J. Brennan, “A comparison of semi-active damping control strategies for vibration isolation of harmonic disturbances,” Journal of Sound and Vibration, vol. 280, no. 1-2, pp. 21–39, 2005.

C. Papalukopoulos and S. Natsiavas, “Nonlinear biodynamics of passengers coupled with quarter car models,” Journal of Sound and Vibration, vol. 304, no. 1-2, pp. 50–71, 2007.

A. C. Luo and A. Rajendran, “Periodic motions and stability in a semi-active suspension system with {MR} damping,” Journal of Vibration and Control, vol. 13, no. 5, pp. 687–709, 2007.

T. P. Waters, Y. Hyun, and M. J. Brennan, “The effect of dual-rate suspension damping on vehicle response to transient road inputs,” Journal of Vibration and Acoustics, vol. 131, no. 1, 8 pages, 2009.

S. L. T. de Souza, I. L. Caldas, and R. L. Viana, “Damping control law for a chaotic impact oscillator,” Chaos, Solitons & Fractals, vol. 32, no. 2, pp. 745–750, 2007.

Z. Sun and S. S. Ge, “Switched and impulsive systems: analysis, design and applications, Zhengguo Li, Yengchai Soh, Changyun Wen; Springer, Berlin. Heidelberg, 2005, ISBN: 3-540-23952-9.,” Automatica, vol. 42, no. 9, pp. 1617–1620, 2006.

G. Xie and L. Wang, “Necessary and sufficient conditions for controllability and observability of switched impulsive control systems,” Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 49, no. 6, pp. 960–966, 2004.

Z.-H. Guan, D. J. Hill, and X. Shen, “On hybrid impulsive and switching systems and application to nonlinear control,” Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 50, no. 7, pp. 1058–1062, 2005.

G. Zong and Y. Wu, “Exponential stability of switched systems with impulsive effect,” Control Theory and Technology, vol. 3, no. 1, pp. 60–66, 2005.

Z. Zhang, “Robust H∞ control of a class of discrete impulsive switched systems,” Nonlinear Analysis, vol. 71, no. 12, pp. e2790–e2796, 2009.

H. Xu, X. Liu, and K. L. Teo, “Delay independent stability criteria of impulsive switched systems with time-invariant delays,” Mathematical and Computer Modelling, vol. 47, no. 3-4, pp. 372–379, 2008.

H. Xu, X. Liu, and K. L. Teo, “A {LMI} approach to stability analysis and synthesis of impulsive switched systems with time delays,” Nonlinear Analysis: Hybrid Systems, vol. 2, no. 1, pp. 38–50, 2008.

X. Liu, S. Zhong, and X. Ding, “Robust exponential stability of nonlinear impulsive switched systems with time-varying delays,” Lithuanian Association of Nonlinear Analysts. Nonlinear Analysis: Modelling and Control, vol. 17, no. 2, pp. 210–222, 2012.

L. P. Zhang, H. B. Jiang, and Y. Liu, “Bifurcation analysis of a rigid impact oscillator with bilinear damping,” in Proceedings of the Abstract, the International Conference on Engineering Vibration, Sofia, Bulgaria, 2017, http://www.icoev2017.org/ocs/index.php/icoev/ICOEV2017/paper/view/242.

J. Ing, E. Pavlovskaia, M. Wiercigroch, and S. Banerjee, “Experimental study of impact oscillator with one-sided elastic constraint,” Philosophical Transactions of the Royal Society A: Mathematical, Physical & Engineering Sciences, vol. 366, no. 1866, pp. 679–704, 2008.

P. Thota and H. Dankowicz, “TC-HAT (TC): A novel toolbox for the continuation of periodic trajectories in hybrid dynamical systems,” SIAM Journal on Applied Dynamical Systems, vol. 7, no. 4, pp. 1283–1322, 2008.

H. Dankowicz and F. Schilder, “An extended continuation problem for bifurcation analysis in the presence of constraints,” Journal of Computational and Nonlinear Dynamics, vol. 6, no. 3, Article ID 031003, 2011.