Research Article  Open Access
Zhiyong Zhang, Thomas Sattel, Aditya Suryadi Tan, Xiaoting Rui, Shaopu Yang, Rui Yang, Yongheng Ying, Xuan Li, "Suppression of Complex Hysteretic Resonances in Varying Compliance Vibration of a Ball Bearing", Shock and Vibration, vol. 2020, Article ID 8825902, 11 pages, 2020. https://doi.org/10.1155/2020/8825902
Suppression of Complex Hysteretic Resonances in Varying Compliance Vibration of a Ball Bearing
Abstract
It is traditionally considered that, due to the Hertzian contact forcedeformation relationship, the stiffness of rolling bearings has stiffening characteristics, and gradually researchers find that the supporting characteristics of the system may stiffen, soften, and even coexist from them. The resonant hysteresis affects the stability and safety of the system, and its jumping effect can make an impact on the system. However, the ball bearing contains many nonlinearities such as the Hertzian contact between the rolling elements and raceways, bearing clearance, and timevarying compliances (VC), leading great difficulties to clarify the dynamical mechanism of resonant hysteresis of the system. With the aid of the harmonic balance and alternating frequency/time domain (HBAFT) method and Floquet theory, this paper will investigate the hysteretic characteristics of the Hertzian contact resonances of a ball bearing system under VC excitations. Moreover, the linearized dynamic bearing stiffness of the system will be presented for assessing the locations of VC resonances, and the nonlinear characteristics of bearing stiffness will also be discussed in depth. Our analysis indicates that the system possesses many types of VC resonances such as the primary, internal, superharmonic, and even combination resonances, and the evolutions of these resonances are presented. Finally, the suppression of resonances and hysteresis of the system will be proposed by adjusting the bearing clearance.
1. Introduction
Rollingelement bearings are characterized by rolling motion of balls or rollers in bearing raceway, often referred to as antifriction bearings in comparison to the sliding motion of regular sleeve bearings [1]. With the advantage of small friction and lowcost maintenance, rolling bearings are the main supporting parts of rotating machinery and machine tools, and also in many cases, the designers use a rolling bearing only because it is easier to select from a manufacturer’s catalog [2]. A major limitation of rolling bearings is that they are subjected to very high alternating and timevarying compliance (VC) oscillating Hertzian stresses at the rolling contacts. Due to the nonlinearities of Hertzian contact and bearing clearance, rolling bearing can bring complex dynamic hysteresis and jumping vibrations to its bearingrotor systems, which is a basic scientific problem in the field of bearingrotor dynamics [3–6].
There are two major effects of a health rolling bearing with respect to machine vibration [1]. The first one is the spring characteristics to support a structure, and the second effect occurs because it is an inherent VC excitation source producing timevarying forces. The VC excitation is an inevitable parametrical excitation to generate vibrations to the rotor system, when the rolling elements carry external load moving along the bearing raceway [7]. Vibrations originating in other parts of the machine can also excite and influence the rolling bearing [2], and in recent years, indepth studies have been carried out gradually on the dynamic influences on some special operation conditions such as starved lubrication [8, 9] and nonlinear instability [5, 6]. As far as VC vibrations are concerned, Sunnersjö [7] studied a linear bearing model focusing on the severity of VC vibrations for a cylindrical roller bearing considering the rotor mass and bearing cage speed, which showed that large movements occur at low speed and heavy load. Fukata et al. [10] presented beating and chaotic VC vibrations at critical speeds, which considered a twodegreeoffreedom ball bearing model including the nonlinearities of the Hertzian contact and bearing clearance. Focused on this classical model, the perioddoubling and quasiperiodic route to chaos in the VC vibrations has been found by Mevel and Guyader [11, 12], and the effects of the unbalance factors [13], bearing clearance [14, 15], and damping and radical loads [16, 17] have been studied subsequently.
Contact nonlinearity is a common nonlinear factor in scientific technology and engineering applications [18, 19]. It has been shown that contact resonances can bring soft resonant hysteresis and jump vibration to the Hertzian contact system with harmonic excitation [20–22]. Kovacic [23] found that constant load has a significant effect on the hysteretic characteristics of the resonance response. Sankaravelu et al. [24] early indicated that dynamic hysteresis and jumps exist in VC vibrations of ball bearings. Generally speaking, a high level of VC vibrations is expected whenever the exciting is close to the Hertzian contact resonant frequencies of a ball bearing system. Elsayed et al. [25] illustrated detailed orders of the linear vibration modes for a ball bearing system through finite element modeling analysis. Zhang et al. [5] found that the primary contact resonances of a balanced rigidrotor ball bearing model can exhibit soft spring behaviors for period1 VC responses, and these nonlinear vibration modes have been confirmed by experiments subsequently [17]. Besides, perioddoubling and chaotic motions can emerge in the VC resonant ranges for a ball bearingrotor system [5, 11, 14]. Dynamic hysteresis as a dangerous bifurcation [26] can bring inside impacts to the system. It has been found that some types of military aircraft had been repaired many times due to the bistability jump faults coming from there aeroengine bearings. Besides, the fatigue life of a rolling contact bearing is a function of the oscillating stress magnitude at the contact [2]. Therefore, the designers must keep in mind to prevent dynamic resonant hysteresis and jumps in bearingrotor system, which will be addressed in this paper.
As a basic parameter of rolling bearings, bearing clearance has great significance for bearing life, installation, and thermal expansion capacity [1], and it can bring remarkable influence on bearing stiffness [2, 16, 27]. Oswald et al. [28] found that bearing life declines gradually with positive clearance and rapidly with increasing negative clearance, and the life can be maximized for a small negative operating clearance. Yakout et al. [29] presented a statistical analysis to study the correlation between the internal radial clearance and the dynamic characteristics of rolling bearings. In the field of bearing industry, it is generally considered that the adverse vibration and noise behaviors can be reduced by adjusting the ball bearings to zero clearance [30]. Particularly, clearancefree operation is usually desirable for precision or highspeed machinery, such as machine tools and turbines [2].
The motivation of the present paper is to investigate the mechanism of VC contact resonances and their hysteresis characteristics in a ball bearing system, considering control of the main structural parameters of the system such as bearing clearances and damping factors. For doing this, as our previous literature [5], the HBAFT method will be applied to trace the VC periodic solution branch and its stability characteristics will be analyzed by using Floquet theory.
2. Characteristics of the Ball Bearings
2.1. Governing Equations of VC Motions
A classic 2degreeoffreedom balanced ball bearingrotor model (see Figure 1), considering the nonlinearities of Hertzian contact and bearing clearance, is verified to be effective for analysis of radial VC motions of a ball bearing [5, 10, 11, 17], and the following assumptions are required [5]:(1)There is an interference fit between the shaft and inner race of the radial ball bearings, and the outer race is fixed to a rigid support(2)The rotor shaft of the studied ball bearingrotor system is assumed to be rigid(3)It is under pure rolling condition between the raceway and balls, and the whole system is faultfree(4)The inertia of balls introduces little dynamic effect on the system
The equation of VC motion can be expressed aswith restoring forces F_{x} and F_{y}andwhere m, c, W, and C_{b} are the equivalent mass, damping factor estimated from [31], constant force, and stiffness coefficient of the bearing system, respectively; G[·] is a Heaviside function representing the contact conditions [5], obtaining the value 1 (contact case) for positive arguments and the value 0 (noncontact case) for negative or zero arguments; δ_{i} and θ_{i} are radial deformation of ith ball (total number of balls, N_{b}) and instant angular location; 2δ_{0} is the radial internal clearance of the bearing; and Ω denotes the cage velocity, which is related to the shaft speed ω_{s}, ball diameter, and pitch diameter of the ball bearing system as
Due to the parametric VC excitation, primary parametric resonance occurs as Ω_{vc} = N_{b}·Ω approaches one of the resonant frequencies of equation (8); that is,where ω_{avg_xx} and ω_{avg_yy}, discussed in the next section, are the averaged first resonant frequencies in the vertical and horizontal radical directions of the system.
2.2. Ball Bearing Nonlinear Stiffness
Bearing stiffness as supporting spring is directly related to load and vibration characteristics of the bearing moving or rotating parts. Due to Hertzian point contact between each rolling element and its raceways of a rolling bearing, the total bearing load W_{s} versus deformation δ_{Ws} along the load direction in the static condition exhibits approximately hard (i.e., stiffening) spring characteristics as [3, 4, 31]where t takes certain t_{0} in equation (4); C is constant coefficient in connection with bearing structure; n takes 3/2 or 10/9 for ball or roller bearing, representing the Hertzian point contact and line contact, respectively, between the rolling elements and raceways [31]; and then the static bearing stiffness k_{Ws} in the load direction can be expressed asand the static lateral stiffness (the bearing stiffness in orthogonal direction to the direction of W_{s}) also can be given (see [31]).
Compared with the static bearing stiffness derived in equation (8), it should be noted that the bearing stiffness is time varying and nonlinear coupling with operation conditions such as rotor speed, VC excitation, and system responses, which can make the operating bearing stiffness exhibit complex hard and soft spring characteristics [3], so it is necessary to consider the dynamic bearing stiffness of the system when specializing in hysteretic resonances.
In order to descript the timevarying bearing support characteristics fully, from equation (2), the linearized dynamic bearing stiffness in x and y directions of system (1) can be approximated as follows [11]:which is obtained by numerical integration of equation (1).
According to equation (9), dynamic resonant frequencies can be expressed asand then, the equivalent values of dynamic resonant frequencies of the system are averaged as [5]where K_{int} is the number of numerical integration steps when solving equation (1).
3. Methodology
The steadystate periodic response of a linear system is always stable, but it may get unstable in a nonlinear dynamic system. The analysis of periodic motion types and their bifurcation mechanism of a nonlinear system is beneficial to abnormal vibration (e.g., super/subharmonic, combination, chaotic, bistable, and hystertic responses) control and even machine design. The HBAFT method combined with Floquet stability analysis is an effective way to trace a VC periodic solution branch and its stability characteristics of ball bearing system [5, 6].
Introducing nondimensional time τ = Ωt and nondimensional displacements X = x/δ_{0} and Y = y/δ_{0} to equation (1), then the following nondimensional motion equation is obtained: where λ is a control parameter. Take the VC excitation period T = 2π of equation (12) as period of the response, X (τ) = X (τ +Τ). Then, X (τ) and F (·) can be represented by Fourier series as
First, inserting equation (13) into equation (12), the harmonic balance yields an algebraic relationshipwhereare the harmonic terms of X (τ) and F (·), respectively.
Second, in order to solve P from equation (14), one needs to formulate Q by P, and this can be achieved by the discrete Fourier transforms as
Herein, P_{k} = a_{k} + ib_{k}; N is the numbers of samplings in the time domain, and ; Q_{k} = c_{k} + id_{k}, and ϕ is 1 if n = 0; otherwise, ϕ is 2; F (n) = F (X^{′} (n), X (n), λ) is the value of the restoring force F (X′ (τ), X (τ), λ) at the nth discrete time.
Third, P can be obtained from equation (14) by fixedpoint iteration method. Besides, the Arclength method is adopted to bridge over the failure iteration at the turning point due to the singular iterative matrix emerging (see [5]).
Finally, Hsu’s method is applied for Floquet stability analysis of obtained periodic solutions aided as shown in our previous work [32].
4. Results and Discussion
4.1. Specifications
The JIS6306 ball bearing is adopted in this paper, which has been studied commonly during recent years [5, 6, 10–12], and the main parameters of the system studied are given in Table 1. Herein, the damping factor c is estimated by [31]where k is the linearized stiffness of the ball bearing. According to Figure 2 in our paper, we take k = 6 × 10^{7}, so the damping factor is in range of 150–1500 Ns/m. If doing experimental research, the bearing damping can be measured by experimental modal analysis referring to literature [33].

In the following, periodic solutions are traced by HBAFT method with the cage speed Ω as the controlled parameter, where a periodic VC motion is called periodn motion if the period of the response is n times the VC parametrically excited period T_{vc} = 2π/Ω_{vc}. Herein, the numerical verifications to the HBAFT tracing results are simulated by the classical explicit RungeKutta numerical integration process of literature [34].
4.2. Complex Resonance and Stiffness Characteristics
For δ_{0} = 6.0 μm c = 200 Ns/m, Figure 3 shows period1 frequencyresponse curves of the system, where vertical and horizontal primary resonant amplitudes are traced, respectively. It is clear that soft hysteretic resonance occurs in the vertical direction since the frequencyresponse curve bends to the left, and the system exhibits coexistence of soft and hard spring characteristics in the horizontal direction, which agrees with previous theoretical and experimental studies [16, 17]. In terms of the bifurcation and stability of the period1 solution branch considered, subcritical perioddoubling bifurcations emerge at A_{1} and A_{2}, cyclic fold bifurcations occur at turning points A_{3} and A_{4}, and secondary Hopf bifurcations include quasiperiodic motions at A_{5} and A_{6}, respectively. These bifurcation characteristics agree well with the numerical simulations as shown in Figure 4, where the blue and magenta dots illustrate the localmaxima numerical bifurcation diagram (as in [35]), and the localmaxima frequencyresponse curves of stable period1 and period2 branches traced by HBAFT method overlapped the maximum value of the numerical bifurcation dots (e.g., see the enlargement around bifurcation point A_{1} in Figure 4).
The cyclic fold bifurcation and the subcritical bifurcation all can lead to instabilities and dynamic jump to the ball bearing systems. Herein, the perioddoubling bifurcations come from onetotwo internal resonances [5, 6], and it is clear that strong couplings between the two degrees of freedom of the system arise in the period2 response range (see range BB in Figure 4). At this time, the VC motions exhibit complex hysteresis characteristics of bistability and jumped phenomenon. The most basic motion behavior is the coexistence of period1 and period2 motions in BB range, for example, the bistability characteristics when Ω = 230 rad/s (see Figures 5 and 6).
(a)
(b)
Meanwhile, as shown in Figure 2, the system possesses two types of timevarying stiffness characteristics, since the dynamic stiffness of the system is related to the response of the system [1, 3]. Therefore, compared to the dynamic stiffness of equation (10), the static stiffness of equation (8) lacks comprehensive perspectives to analyze the stiffness nonlinearity of the system. Moreover, it is worth noting that the coexistence of various stiffness characteristics makes it more difficult to predict and adjust the hysteretic resonances. For example, as shown in Figure 7(a), although the primary resonant hysteresis of period1 VC motion has been suppressed as the damping coefficient c increasing to 1400 Nm/s, the system still holds hysteretic jumping behaviors due to the nonlinear dynamic stiffness related to the period2 motions (see Figure 7(b)).
(a)
(b)
Unstable period 1 solution branch A_{5}A_{6} in Figure 4 is induced by secondary Hopf bifurcating instability and strongly coupled vibrations produced between the two degrees of freedom of the system in CC range. Moreover, the complex characteristics of bistabilities and jumps also emerge in this range. For instance, as shown in Figures 8 and 9, quasiperiodic and chaotic motions coexist at Ω = 174.3 rad/s. In terms of the frequency spectrum of the quasiperiodic response considered (see Figure 9(a)), the values of incommensurable dimensionless VC frequencies p and q are about 0.663 and 1.337, respectively, and (p + q)/2 = 1, where the value 1 is the dimensionless VC excitation frequency Ω_{vc} = N_{b}·Ω. The values p·Ω = 115.561 rad/s and q·Ω = 233.039 rad/s locate near the resonant excitation frequency ranges BB and DD, respectively, so it is reasonable to suppose that the case of combination resonance ω_{p} + ω_{q} = 2·Ω_{vc} occurs in the segment of CC. It has been found in some pieces of literature [5, 6, 10–16] that the quasiperiodic motion is a typical characteristic of VC vibrations in rolling bearing system, but so far, its physical trigger mechanism is not given. In general, the typical combination resonance of multipledegreeoffreedom system often can arise in quasiperiodic motions [36–38]. Therefore, it can be considered that the combination resonance may be one of the physical mechanisms of generation of VC quasiperiodic motions. In addition, the existence of combination resonance makes the response characteristics of the system more complicated.
(a)
(b)
(c)
(a)
(b)
4.3. Control of Resonances and Hysteresis
In certain applications, particularly in highspeed machines or high precision machines, such as machine tools, it is possible to eliminate clearance and introduce small interference (negative clearance) in the bearings, because it increases the bearing stiffness and reduces the runout and noise due to elastic deformation or clearance [2, 30]. However, the interference should be ensured not to be excessive, because the interference can result in extrarolling contact pressure to reduce the fatigue life of the bearing. Figure 10 shows the change of the resonant amplitude of the two degrees of freedom of the system as bearing clearance δ_{0} reduces from 6 μm to 0 μm. Herein, the resonant amplitude in y direction moves to high frequency as δ_{0} reduces to 0 μm, because a small bearing clearance can stiffen the support of the shaft [2], and it is worth noting that the resonant response is suppressed effectively and hard resonant hysteresis emerges when δ_{0} = 2 μm. Moreover, the combination and onetotwo internal resonances are eliminated in the condition of clearancefree operation (i.e., δ_{0} = 0 μm). However, the resonant hysteresis and jumps still exist in the two degrees of freedom of the system.
One can see in Figure 11(a) that the soft hysteresis in vertical resonant range gradually eliminates as the bearing clearance δ_{0} reduces from 0.2 μm to −0.6 μm, and then the vertical resonant amplitude gets larger as δ_{0} continues to reduce from −0.6 μm to −1.2 μm (see Figure 11)(b). This difference is due to the variation of the stiffness characteristics. As the bearing clearance decreases to −0.8 μm, −1.0 μm, and −1.2 μm, the stiffness in the horizontal direction is close to the stiffness in the vertical direction, which induces strong coupling resonance between the two degrees of freedom of the system. For example, when δ_{0} takes −1.2 μm, the two degrees of freedom are remarkably coupled in horizontal and vertical resonant ranges (see Figure 12(a)), respectively, and combination resonance again arises between these two resonant ranges (see E_{1}E_{2} branch in Figure 12(b)). In other words, the vibration amplitudes, bearing dynamic loads, and even the bearing life can be optimized for a proper small negative operating clearance, which agrees well with the results of literature [28].
(a)
(b)
(a)
(b)
In terms of hysteresis characteristics considered, due to the effect of clearance nonlinearity elimination as δ_{0} is negative, at this moment, the Hertzian contact nonlinearity of the system plays a dominant role. Therefore, as illustrated in Figure 12(a), the resonant amplitudes nearly all bend to left as the Hertzian contact deformation exhibits a stiffening nonlinearity relationship. In addition, compared with resonant responses for positive bearing clearance cases shown in Figure 10, for δ_{0} = −1.2 μm, the resonant excitation frequency is higher in the horizontal direction than in the vertical direction (see Figure 12(a)), which agrees with the equivalent dynamic resonant frequencies predicted by equation (11) (see Figure 13(a)), and Figure 13(b) illustrates that selection of proper bearing clearance can suppress the hysteresis in the vertical direction effectively resonance and make the system operate in a quasilinear stiffness interval.
(a)
(b)
Figure 14 gives the power spectra of the dynamic displacement parameter (DDP) [7]: DDP = x_{peaktopeak} + y_{peaktopeak}, where x_{peaktopeak} = x (t)_{max}x (t)_{min}, and y_{peaktopeak} = y (t)_{max}y (t)_{min}. In order to reflect the hysteresis characteristics in this figure, we get the DDP value when the controlled parameter Ω runs up and down by the way in literature [6]. In general, as shown in Figure 14, the system contains many types of VC resonant characteristics, such as primary resonances R_{1} and R_{2} in vertical and horizontal directions, respectively, and also superharmonic resonances R_{3} and R_{4} are aroused in these two directions. Moreover, the forms of ω_{p} + ω_{q} = 2·Ω_{vc} and 4·ω_{p} + ω_{q} = 8·Ω_{vc} combination resonances labeled as R_{5} and R_{6}, respectively, exist for this system. Figure 15 shows that equation (11) can provide a good prediction of typical resonance locations, where ω_{p} = ω_{avg_yy} and ω_{q} = ω_{avg_xx}, so it is a valuable way for analysis of VC resonances for the ball bearing systems.
5. Conclusions
This paper focuses on the characteristics of Hertzian contact resonances and dynamic hysteresis of a ball bearing system. The HBAFT method and Floquet theory are applied to trace the periodic response behaviors of the system. It is found that the system contains many types of resonances such as primary, combination, and superharmonic resonances, which makes the response characteristics of the system more complicated, and the linearized dynamic bearing stiffness can be used to assess the locations of these resonances. It should be emphasized that combination resonance can bring quasiperiodic motions to the system. Moreover, the coexistence of two types of timevarying stiffness characteristics is demonstrated. It is presented that bearing clearance has significant effects on the dynamic characteristics of ball bearing system, including the bearing stiffness, resonant amplitudes, and their hysteretic behaviors. Our analysis indicates that the hysteretic resonances can be suppressed effectively by adjusting the bearing clearances around the clearancefree operations, and this method may be beneficial for the engineering study and inspection of resonant vibration control and even bearing life optimization of rolling bearing systems. In the future, the intrinsic correlations between typical resonances, Hertzian contact nonlinearity, bearing clearance, and even asymmetric stiffness characteristics of the system should be further clarified.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work has been supported by the China Postdoctoral Science Foundation (Grant no. 2019M661849), Jiangsu Planned Projects for Postdoctoral Research Funds (Grant no. 2018K164C), Fundamental Research Funds for the Central Universities (Grant no. 30920021155), China Scholarship Council (ID: 201906845022), and Natural Science Foundation of Jiangsu Province (Grant no. BK20190434).
References
 T. A. Harris, Rolling Bearing Analysis, John Wiley & Sons, Hoboken, NJ, USA, 4th edition, 2001.
 A. Harnoy, Bearing Design in Machinery: Engineering Tribology and Lubrication, Marcel Dekker, New York, NY, USA, 2002.
 F. F. Ehrich, Handbook of Rotordynamics, McGrawHill, New York, NY, USA, 1992.
 M. L. Adams, Rotating Machinery Vibration: From Analysis to Troubleshooting, Taylor & Francis, Oxfordshire, UK, 2009.
 Z. Zhang, Y. Chen, and Q. Cao, “Bifurcations and hysteresis of varying compliance vibrations in the primary parametric resonance for a ball bearing,” Journal of Sound and Vibration, vol. 350, pp. 171–184, 2015. View at: Publisher Site  Google Scholar
 Z. Y. Zhang, X. T. Rui, R. Yang, and Y. S. Chen, “Control of perioddoubling and chaos in varying compliance resonances for a ball bearing,” ASME Journal of Applied Mechanics, vol. 87, Article ID 021005, p. 2020. View at: Google Scholar
 C. S. Sunnersjö, “Varying compliance vibrations of rolling bearings,” Journal of Sound and Vibration, vol. 58, no. 3, pp. 363–373, 1978. View at: Publisher Site  Google Scholar
 H. T. Shi, X. T. Bai, K. Zhang et al., “Influence of uneven loading condition on the sound radiation of starved lubricated full ceramic ball bearings,” Journal of Sound and Vibration, vol. 461, Article ID 114910, 2019. View at: Google Scholar
 H. T. Shi and X. T. Bai, “Modelbased uneven loading condition monitoring of full ceramic ball bearings in starved lubrication,” Mechanical Systems and Signal Processing, vol. 139, Article ID 106583, 2020. View at: Google Scholar
 S. Fukata, E. H. Gad, T. Kondou, T. Ayabe, and H. Tamura, “On the radial vibration of ball bearings: computer simulation,” Bulletin of JSME, vol. 28, no. 239, pp. 899–904, 1985. View at: Publisher Site  Google Scholar
 B. Mevel and J. L. Guyader, “Routes to chaos in ball bearings,” Journal of Sound and Vibration, vol. 162, no. 3, pp. 471–487, 1993. View at: Publisher Site  Google Scholar
 B. Mevel and J. L. Guyader, “Experiments on routes to chaos in ball bearings,” Journal of Sound and Vibration, vol. 318, no. 3, pp. 549–564, 2008. View at: Publisher Site  Google Scholar
 M. Tiwari, K. Gupta, and O. Prakash, “Dynamic response of an unbalanced rotor supported on ball bearings,” Journal of Sound and Vibration, vol. 238, no. 5, pp. 757–779, 2000. View at: Publisher Site  Google Scholar
 M. Tiwari, K. Gupta, and O. Prakash, “Effect of radial internal clearance of a ball bearing on the dynamics of a balanced horizontal rotor,” Journal of Sound and Vibration, vol. 238, no. 5, pp. 723–756, 2000. View at: Publisher Site  Google Scholar
 C.q. Bai, Q.y. Xu, and X.l. Zhang, “Nonlinear stability of balanced rotor due to effect of ball bearing internal clearance,” Applied Mathematics and Mechanics, vol. 27, no. 2, pp. 175–186, 2006. View at: Publisher Site  Google Scholar
 Z. Zhang, Y. Chen, and Z. Li, “Influencing factors of the dynamic hysteresis in varying compliance vibrations of a ball bearing,” Science China Technological Sciences, vol. 58, no. 5, pp. 775–782, 2015. View at: Publisher Site  Google Scholar
 Y. L. Jin, R. Yang, L. Hou, Y. S. Chen, and Z. Y. Zhang, “Experiments and numerical results for varying compliance contact resonance in a rigid rotorball bearing system,” ASME Journal of Tribology, vol. 139, Article ID 041103, 2017. View at: Publisher Site  Google Scholar
 R. A. Ibrahim, Vibroimpact Dynamics: Modelling, Mapping and Applications, SpringerVerlag, Berlin, Germany, 2009.
 C. F. Li, Z. X. Zhang, Q. Y. Yang, and P. Y. Li, “Experiments on the geometrically nonlinear vibration of thinwalled cylindrical shells with points supported boundary condition,” Journal of Sound and Vibration, vol. 473, Article ID 115226, 2020. View at: Google Scholar
 R. Nayak, “Contact vibrations,” Journal of Sound and Vibration, vol. 22, no. 3, pp. 297–322, 1972. View at: Google Scholar
 D. P. Hess and A. Soom, “Normal vibrations and friction under harmonic loads: Part Ihertzian contacts,” Journal of Tribology, vol. 113, no. 1, pp. 80–86, 1991. View at: Publisher Site  Google Scholar
 E. Rigaud and J. PerretLiaudet, “Experiments and numerical results on nonlinear vibrations of an impacting hertzian contact. Part 1: harmonic excitation,” Journal of Sound and Vibration, vol. 265, no. 2, pp. 289–307, 2003. View at: Publisher Site  Google Scholar
 I. Kovacic, “Forced vibrations of oscillators with a purely nonlinear powerform restoring force,” Journal of Sound and Vibration, vol. 330, no. 17, pp. 4313–4327, 2011. View at: Publisher Site  Google Scholar
 A. Sankaravelu, S. T. Noah, and C. P. Burger, “Bifurcation and chaos in ball bearings,” in Proceedings of ASME Nonlinear and Stochastic Dynamics, Chicago, USA, August 1994. View at: Google Scholar
 E. S. Elsayed, A. M. Elkhatib, and M. Yakout, “Vibration modal analysis of rolling element bearing,” in Proceedings of 4th International Conference on Integrity, Reliability and Failure, Funchal, Madeira, June 2013. View at: Google Scholar
 J. M. T. Thompson, H. B. Stewart, and Y. Ueda, “Safe, explosive, and dangerous bifurcations in dissipative dynamical systems,” Physical Review E, vol. 49, no. 2, pp. 1019–1027, 1994. View at: Publisher Site  Google Scholar
 S. H. Ghafari, E. M. AbdelRahman, F. Golnaraghi, and F. Ismail, “Vibrations of balanced faultfree ball bearings,” Journal of Sound and Vibration, vol. 329, no. 9, pp. 1332–1347, 2010. View at: Publisher Site  Google Scholar
 F. B. Oswald, E. V. Zaretsky, and J. V. Poplawski, “Effect of internal clearance on load distribution and life of radially loaded ball and roller bearings,” Tribology Transactions, vol. 55, no. 2, pp. 245–265, 2012. View at: Publisher Site  Google Scholar
 M. Yakout, M. G. A. Nassef, and S. Backar, “Effect of clearances in rolling element bearings on their dynamic performance, quality and operating life,” Journal of Mechanical Science and Technology, vol. 33, no. 5, pp. 2037–2042, 2019. View at: Publisher Site  Google Scholar
 Schaeffler Technologies AG & Co KG, The Design of Rolling Bearing Mountings: Design Examples Covering Machines, Vehicles and Equipment, Fag Bearing company Ltd., Schweinfurt, Germany, 2012.
 K. Erwin, Dynamics of Rotor and Foundation, SpringerVerlag, Berlin, Germany, 1993.
 Z.y. Zhang and Y.s. Chen, “Harmonic balance method with alternating frequency/time domain technique for nonlinear dynamical system with fractional exponential,” Applied Mathematics and Mechanics, vol. 35, no. 4, pp. 423–436, 2014. View at: Publisher Site  Google Scholar
 M. Yakout, A. Elkhatib, and M. G. A. Nassef, “Rolling element bearings absolute life prediction using modal analysis,” Journal of Mechanical Science and Technology, vol. 32, no. 1, pp. 91–99, 2018. View at: Publisher Site  Google Scholar
 J. R. Cash and A. H. Karp, “A variable order RungeKutta method for initial value problems with rapidly varying righthand sides,” ACM Transactions on Mathematical Software (TOMS), vol. 16, no. 3, pp. 201–222, 1990. View at: Publisher Site  Google Scholar
 G. A. Held, C. Jeffries, and E. E. Haller, “Observation of chaotic behavior in an electronhole plasma in Ge,” Physical Review Letters, vol. 52, no. 12, pp. 1037–1040, 1984. View at: Publisher Site  Google Scholar
 A. H. Nayfeh and D. T. Mook, Nonlinear Oscillations, John Wiley & Sons, Hoboken, NJ, USA, 1995.
 L. N. Panda and R. C. Kar, “Nonlinear dynamics of a pipe conveying pulsating fluid with parametric and internal resonances,” Nonlinear Dynamics, vol. 49, no. 12, pp. 9–30, 2007. View at: Publisher Site  Google Scholar
 S. Bamadev, L. N. Panda, and G. Pohit, “Combination, principal parametric and internal resonances of an accelerating beam under two frequency parametric excitation,” International Journal of Nonlinear Mechanics, vol. 78, pp. 35–44, 2016. View at: Google Scholar
Copyright
Copyright © 2020 Zhiyong 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.