# Analysis on Multimode Nonlinear Resonance of an Asymmetric Rolling Bearing Rotor with the Application of a Squeeze Film Damper

Authors: Chen Huizheng; Zhong Shun; Lu Zhenyong; Chen Yushu; Liu Xiyu
Journal: Shock and Vibration
 (2021)

---

## Body

## 1. Introduction

Because of the high working rotating speed, the rotor structures in aeroengines [1–3] are supported by rolling bearings [4–8]. Rolling bearings are a type of main support for aeroengines and are a key source of vibration (parametric excitation) in rotor-bearing systems, which create a demand for vibration analysis and diagnostic techniques under various operating conditions [9]. Studying the nonlinear behaviour of the system and identifying the influence of rolling element bearing parameters on the nonlinear behaviours of the rotor system will help to reasonably avoid system instability, increase the operating life, and improve the efficiency of mechanical equipment [10]. Thus, the rolling bearing supported rotor systems have drawn many research interests [11–17].The bearing introduced nonlinearity which makes the research complex. For a nonlinear system, not only the main resonances of each order could occur but also the nonlinear resonances, which are more complicated and wider in range, will appear. The early research of the nonlinear phenomena in rolling bearings could be traced back to 80s. Fukuta et al. [18] found that ball bearing systems have nonlinear dynamic behaviours such as superharmonic, subharmonic, quasi-periodic, and chaos-like solutions. Mevel and Guyader [19, 20] considered that the system entered chaos through a period-doubling bifurcation and explained the relationship between chaos motion and the contact failure. The resonance types of a multi-degree-of-freedom rolling bearing-rotor system mainly contain (1) prime resonances of each order, that is, the resonance when the external frequency is equal to the natural frequency of the system, (2) superharmonic or subharmonic resonances, which are caused by the nonlinear relationship between the excitation and response in the nonlinear system, (3) quasi-periodic resonances, and (4) combination resonances. Therefore, it is necessary to add damping structures to suppress the unexpected resonances. The squeeze film damper (SFD) is one of the dampers usually used in aeroengines.As a support structure with damping function, the SFD has been widely used in aeroengines due to its small size and light weight and has even become a standard configuration in the design. Many researchers focused on the design of SFD and its application effects to the rotor system. Zhang et al. [21] presented a multiobjective optimal design method for a squeeze film damper with centering spring. The proposed method was able to obtain the optimal design parameters for a flexible rotor system. Shaik and Dutta [22, 23] proposed a close form solution to find the nonlinear tuning parameters of symmetric or asymmetric rotor-squeeze film damper system. They also dealt with the development of a close form solution to find the stability of flexible symmetric horizontal rigid and tuned flexible rotor mounted on hydrodynamic bearing with squeeze film damper system. Zheng et al. [24] proposed a controllable clearance squeeze film damper which used hydraulic pressure to adjust the radial clearance and established a corresponding test rig. Fan and Behdinan [25] investigated the effect of a circumferential central groove on an open-ended squeeze film damper analytically. Shin et al. [26] treated the unconventional application of the SFD for the mitigation of Morton effect-induced vibration. Ma et al. [27] emphasized the comparative analysis of the influence of SFD on the nonlinear dynamic behaviour of the dual-rotor system supported by rolling bearings. Chen et al. [28] investigated dynamic characteristics of the rotor system mounted on aircraft during maneuvering flights. Iacobellis et al. [29] designed and built SFD-rotor test rig to study the effect of SFD oil supply pressure, oil temperature, oil inlet feed number/orientation, unbalance, and seals on the response of a Jeffcott rotor. Li et al. [30] manufactured and tested a pair of 3.5 inch SFD bearings to validate a new squeeze film damper bearing design.From the existing work, the complexity of the response of rolling bearings, the damping effect of squeeze film dampers, and the nonlinear stiffness imposed on the nonlinear damping system all need further study. To further investigate the dynamical behaviours and achieve better or desired performances of the rotor systems with SFD, so as to provide a theoretical basis for the corresponding vibration reduction design strategy, the design of the parameters of the SFD need focused attention. Therefore, this work will focus on the dynamic analysis of the ball bearing supported rotor system with the application of the SFD and their parameter analysis. The paper is organized as follows. In Section2, a ball bearing-rotor system model with SFD installed is established. The detailed modeling processes of the bearing and oil force are described. In Section 3, the dynamic analysis of the ball bearing-rotor system with or without SFD is conducted. In Section 4, parameters’ analysis of the ball bearing act on the dynamic behaviour is given out. Section 5 discusses the comparison results, and Section 6 concludes the whole work.


## 2. Modeling of the Rotor System with SFD

Figure1 shows a schematic diagram of a Jeffcott rotor system supported by ball bearings. The left and right ends have different ball bearing support schemes, respectively, as the left end is connected to the foundation with a spring and the right end is directly connected to the foundation without elastic structure. So, the degree of freedom of the outer ring of the left bearing must be considered while the right one can be ignored.Figure 1
Schematic diagram of a Jeffcott rotor.Considering the gravity, unbalanced force, gyroscopic effect, and bearing force, the dynamic model of the rotor support structure shown in Figure1 can be established into form(1)MX..+C+ΩGX.+KX+Fn=Fu+Fg,where X=x,y,θx,θy,xa,ya,xb,yb,xo,yoT is displacement coordinate vector, in which x and y are disc displacements along the x-axis and y-axis, respectively, θx and θy are disc rotating angles along the x-axis and y-axis, respectively, xa and ya represent the displacements of the equivalent mass at the left end journal, xb and yb represent the displacements of the equivalent mass at the right end, and xo and yo represent the displacements of the outer ring of the bearing mounted at the left. M, C, G, and K are corresponding mass matrix, damping matrix, gyro matrix, and stiffness matrix, whose expressions can be found in appendix. Fn=0,0,0,0,Fbx1,Fby1,Fbx2,Fby2,Fcx−Fbx1,Fcy−Fby1T is the bearing force vector, Fu=mδΩ2cosΩt,mδΩ2sinΩt,0,0,0,0,0,0,0,0T is unbalanced force vector, and Fg=0,−mg,0,0,0,−mag,0,−mbg,0,−mogT is gravity vector. Ω is rotating frequency which determines the frequency of the unbalanced force.The bearing force model adopted in this work is a two-degree-of-freedom ball bearing model, including nonlinear Hertzian contact force, clearance, and variable stiffness characteristics. Assuming that the initial positions of the bearing balls at both ends are the same, the position of the balli at time t can be expressed as(2)θi=2πi−1Nb+Ωct,where Nb is the number of the rolling balls and Ωc=ri/ri+ro·Ω, in which, ri is the radius of the inner ring and ro is the radius of the outer ring.Then, the bearing force of the left end can be expressed as(3)Fbx1Fby1=Cb∑i=1Nbδ1iHδ1i3/2cosθ1isinθ1i,where δ1i=xa−xocosθi+ya−yosinθi−δ0, Cb is the Hertz contact stiffness, H· is Heaviside function, and δ0 is the radial clearance of bearing.Meanwhile, the bearing force of the right end can be expressed as(4)Fbx2Fby2=Cb∑i=1Nbδ2iHδ2i3/2cosθ2isinθ2i,where δ2i=xbcosθi+ybsinθi−δ0.The oil film force adopts a short bearing approximate model, and the expression is(5)Fcx=Frxor−Fτyor,Fcy=Fryor+Fτxor,where(6)Fr=μRL3cl3I302r˙+I311rψ˙,Fτ=μRL3cl3I311r˙+I320rψ˙,where r is the radial displacement of the journal center at each support, and its expression is(7)r=xo2+yo2,ψ is the whirl angle of each journal, and the expression is(8)ψ=arctanyoxo,μ is the viscosity coefficient of the lubricating oil at each journal, R is the radius of each journal bearing, L is the length of each damper.(9)Inlm=∫θ1θ2sinlθcosmθ1+rcosθ/cl3dθ,is the Sommerfeld coefficient.Defining the dimensionless time asτ=Ωt and dimensionless displacement vector Q=EX, in which E=1/cdiag1,1,l,l,1,1,1,1,1,1; equation (1) can be transformed into(10)Q″+C¯+G¯Q′+KQ¯+F¯n=F¯u+F¯g,where C=¯1/ΩEM−1CE−1 is dimensionless damping matrix, G=EM¯−1GE−1 is the dimensionless gyro matrix, K=¯1/Ω2EM−1KE−1 is the dimensionless stiffness matrix, F¯n=1/Ω2EM−1Fn is the dimensionless bearing force vector, F¯u=1/Ω2EM−1Fu is the dimensionless unbalanced force vector, and F¯g=1/Ω2EM−1Fg is the dimensionless gravity vector. The superscript prime denotes the derivative with respect of dimensionless time. The values of the system parameters are listed in Table 1. It should be mentioned that, as the stiffness is time varying, the mean of the response valves are used here as(11)kxx=1N∑i=1N∂Fx∂xiΔτ,kyy=1N∑i=1N∂Fy∂yiΔτ,where N=T/Δτ is the number of discrete points.Table 1
The values of the parameters.
ParametersValuesRotor parametersm120 kgma12 kgmb18 kgkrr7.4×107 N/mkrφ and kφr2.77×107 N/mkφφ6.15×107 N/mka1.0×107 N/mc1266 Ns/mc2 and c399.6 Ns/mc4221 Ns/ml10.8 ml21.25 mJd2.5 kg·m2Jp2.5 kg·m2δ1.67×10−6 mBearing parametersri39.6 mmro70.4 mmmo0.965 kgCb1.0×109 N/m3/2


## 3. Dynamical Analysis with the Application of SFD

### 3.1. Main Resonance and Quasi-Periodic Motion

In this section, the amplitude-frequency characteristics of the structure are going to be investigated. A dimensionless excitation frequency defined asλ=Ω/300 is introduced for further analysis. Figure 2 gives the amplitude-frequency curves of the responses of the disc in horizontal direction with different parameters described by the governing equations of 1 with or without SFD. Obviously, the SFD plays a significant role in reducing the amplitude in the first-order main resonance region of the system. Accordingly, the resonance band as well as jump-bistable region reduces with the application of the SFD.Figure 2
The amplitude-frequency curves in first-order primary resonance region.
(a)(b)(c)Figure3 gives the amplitude-frequency curves in the frequency region which is from λ=1.5 to λ=7.0. The blue dashed line in Figure 3 indicates the amplitude-frequency curve when the SFD is not applied. Quasi-periodic solutions exist in the third-order resonance region of λ=5.96 and the region from λ=2.0 to λ=4.0, which covers the second-order main resonance region. The thick solid line indicates the amplitude-frequency response after adding the SFD. Although being suppressed, the quasi-periodic solutions which are located on the third-order resonance region still have large amplitudes. The spectrum is shown in Figures 3(b) and 3(c), where the quasi-periodic response of λ=6.09 not only corresponds to the first-order resonance λf=0.93 but also to the second-order resonance of λf=2.68, by introducing a definition of an absolute quasi-periodic frequency λf=λ·ω, where ω is the frequency of the response.Figure 3
The amplitude-frequency curves and spectrum of responses of disc in horizontal direction whenNb=8. (a) Amplitude-frequency curves, (b) spectrum at λ=6.09, and (c) spectrum at λ=6.38 (blue dot line means no SFD installed, solid line means the SFD is mounted, blue means run-up process, and brown means run-down process).In the quasi-periodic motion, fromλ=2.0 to λ=4.0, the damping effect of the SFD is significant. The large and continuous quasiperiodic motion is reduced to a relatively independent area with a small amplitude. For the horizontal direction of the disc, there is a relatively smooth area between the first-order and third-order resonance regions. In addition, in the second-order resonance region (λ=2.6 to 2.76), a cross structure with soft and hard characteristics coexisting appears.Figure4 gives amplitude-frequency curves of response of the left-end journal in horizontal direction in the same frequency region as Figure 3(a). There are obvious resonance peaks in this region, as well as some jump phenomena. The frequency spectrum analysis is performed around the second-order main resonance (λ=2.67). The leftmost resonant peak in Figure 4 is near λ=1.88, and the corresponding frequency component is shown in Figure 5(a). The dominant frequency is 1/2 VC frequency. Figure 5(b) corresponds to the second-order resonance frequency at λ=2.67. It can be seen that the second-order resonance area is periodic motion and there is jump behaviour. In addition, there are obvious jump phenomena in the region from λ=5.25 to 5.6. Taking the two solutions at λ=5.42, the spectrums are shown in Figures 5(c) and 5(d). The run-down curve is a periodic motion, while the run-up curve has subharmonic frequency component. Therefore, the resonance peak is the subharmonic resonance peak of the second-order natural frequency. Compared with the results under the condition of no SFD, subharmonic resonance of the second-order frequency also occurs, indicating that this is a phenomenon caused by the nonlinearity of the bearing.Figure 4
The amplitude-frequency curves of responses of left-end journal in horizontal direction whenNb=8 (blue dot line means no SFD installed, solid line means the SFD is mounted, blue means run-up process, and brown means run-down process).
(a)(b)(c)(d)Figure 5
Spectrums of the responses at (a)λ=1.88, (b) λ=2.67, (c) λ=5.42 (run-up), and (d) λ=5.42 (run-down) when Nb=8.
(a)(b)(c)(d)Corresponding to the peaks of the amplitude-frequency curves atλ=2.05, λ=2.86, λ=2.97, and λ=3.58 in Figure 4, the spectrums are shown in Figure 6. Besides the excitation frequency and the VC frequency, there are unconventional frequency components. These components are quasi-periodic motions, and the corresponding absolute quasi-periodic frequencies are (a) λf=2.677, (b) λf=2.688, (c) λf=2.682, and (d) λf=0.918. The first three relate to the second-order main resonance.Figure 6
Spectrums of the responses at (a)λ=2.05, (b) λ=2.86, (c) λ=2.97, and (d) λ=3.58 when Nb=8.
(a)(b)(c)(d)(e)(f)


### 3.2. VC Frequency Resonance and Its 1/2 Subresonance

Figure7 shows the amplitude-frequency curves during the VC frequency resonance region with different numbers of balls when the SFD is implemented. With the same damper configuration, the vertical resonance peak no longer shows the jump phenomenon of the soft stiffness characteristics as indicated in [10], but the horizontal direction still maintains the cross structure of the coexisting of the soft and the hard stiffness characteristics. This means that the strong coupling between two directions is not easy to dissipate the vibration energy. The SFD has no obvious effect on the amplitude and resonance zone. This is because, under the small amplitude vibration case, the nonlinear damping characteristics is not easy to exert its vibration. Meanwhile, different numbers of balls does not lead the vibration characteristics change, compared to no SFD cases.Figure 7
The amplitude-frequency curves in VC frequency resonance region, where° indicates run-up and × indicates run-down.
(a)(b)(c)(d)(e)(f)In Figure8, the influence of the SFD on the 1/2 subharmonic VC frequency resonance is illustrated by the reduction or disappearance of the hysteresis areas. When Nb=13, the 1/2 subharmonic VC frequency resonance disappears. Basically, the SFD does not change the subharmonic resonance range of the original system. In addition, compared to [10], when Nb=8, the quasi-period with large amplitude near λ=6.78 disappears after adding the SFD. That means the 1/2 subharmonic resonance of the VC frequency can be eliminated by adjusting the number of balls. In case that the quasi-period exists, the SFD could improve the jump phenomenon in the subharmonic resonance region of the VC frequency, but its amplitude suppression effect is limited.Figure 8
The amplitude-frequency curves in subharmonic VC frequency resonance region, where° indicates run-up and × indicates run-down.


## 3.1. Main Resonance and Quasi-Periodic Motion

In this section, the amplitude-frequency characteristics of the structure are going to be investigated. A dimensionless excitation frequency defined asλ=Ω/300 is introduced for further analysis. Figure 2 gives the amplitude-frequency curves of the responses of the disc in horizontal direction with different parameters described by the governing equations of 1 with or without SFD. Obviously, the SFD plays a significant role in reducing the amplitude in the first-order main resonance region of the system. Accordingly, the resonance band as well as jump-bistable region reduces with the application of the SFD.Figure 2
The amplitude-frequency curves in first-order primary resonance region.
(a)(b)(c)Figure3 gives the amplitude-frequency curves in the frequency region which is from λ=1.5 to λ=7.0. The blue dashed line in Figure 3 indicates the amplitude-frequency curve when the SFD is not applied. Quasi-periodic solutions exist in the third-order resonance region of λ=5.96 and the region from λ=2.0 to λ=4.0, which covers the second-order main resonance region. The thick solid line indicates the amplitude-frequency response after adding the SFD. Although being suppressed, the quasi-periodic solutions which are located on the third-order resonance region still have large amplitudes. The spectrum is shown in Figures 3(b) and 3(c), where the quasi-periodic response of λ=6.09 not only corresponds to the first-order resonance λf=0.93 but also to the second-order resonance of λf=2.68, by introducing a definition of an absolute quasi-periodic frequency λf=λ·ω, where ω is the frequency of the response.Figure 3
The amplitude-frequency curves and spectrum of responses of disc in horizontal direction whenNb=8. (a) Amplitude-frequency curves, (b) spectrum at λ=6.09, and (c) spectrum at λ=6.38 (blue dot line means no SFD installed, solid line means the SFD is mounted, blue means run-up process, and brown means run-down process).In the quasi-periodic motion, fromλ=2.0 to λ=4.0, the damping effect of the SFD is significant. The large and continuous quasiperiodic motion is reduced to a relatively independent area with a small amplitude. For the horizontal direction of the disc, there is a relatively smooth area between the first-order and third-order resonance regions. In addition, in the second-order resonance region (λ=2.6 to 2.76), a cross structure with soft and hard characteristics coexisting appears.Figure4 gives amplitude-frequency curves of response of the left-end journal in horizontal direction in the same frequency region as Figure 3(a). There are obvious resonance peaks in this region, as well as some jump phenomena. The frequency spectrum analysis is performed around the second-order main resonance (λ=2.67). The leftmost resonant peak in Figure 4 is near λ=1.88, and the corresponding frequency component is shown in Figure 5(a). The dominant frequency is 1/2 VC frequency. Figure 5(b) corresponds to the second-order resonance frequency at λ=2.67. It can be seen that the second-order resonance area is periodic motion and there is jump behaviour. In addition, there are obvious jump phenomena in the region from λ=5.25 to 5.6. Taking the two solutions at λ=5.42, the spectrums are shown in Figures 5(c) and 5(d). The run-down curve is a periodic motion, while the run-up curve has subharmonic frequency component. Therefore, the resonance peak is the subharmonic resonance peak of the second-order natural frequency. Compared with the results under the condition of no SFD, subharmonic resonance of the second-order frequency also occurs, indicating that this is a phenomenon caused by the nonlinearity of the bearing.Figure 4
The amplitude-frequency curves of responses of left-end journal in horizontal direction whenNb=8 (blue dot line means no SFD installed, solid line means the SFD is mounted, blue means run-up process, and brown means run-down process).
(a)(b)(c)(d)Figure 5
Spectrums of the responses at (a)λ=1.88, (b) λ=2.67, (c) λ=5.42 (run-up), and (d) λ=5.42 (run-down) when Nb=8.
(a)(b)(c)(d)Corresponding to the peaks of the amplitude-frequency curves atλ=2.05, λ=2.86, λ=2.97, and λ=3.58 in Figure 4, the spectrums are shown in Figure 6. Besides the excitation frequency and the VC frequency, there are unconventional frequency components. These components are quasi-periodic motions, and the corresponding absolute quasi-periodic frequencies are (a) λf=2.677, (b) λf=2.688, (c) λf=2.682, and (d) λf=0.918. The first three relate to the second-order main resonance.Figure 6
Spectrums of the responses at (a)λ=2.05, (b) λ=2.86, (c) λ=2.97, and (d) λ=3.58 when Nb=8.
(a)(b)(c)(d)(e)(f)


## 3.2. VC Frequency Resonance and Its 1/2 Subresonance

Figure7 shows the amplitude-frequency curves during the VC frequency resonance region with different numbers of balls when the SFD is implemented. With the same damper configuration, the vertical resonance peak no longer shows the jump phenomenon of the soft stiffness characteristics as indicated in [10], but the horizontal direction still maintains the cross structure of the coexisting of the soft and the hard stiffness characteristics. This means that the strong coupling between two directions is not easy to dissipate the vibration energy. The SFD has no obvious effect on the amplitude and resonance zone. This is because, under the small amplitude vibration case, the nonlinear damping characteristics is not easy to exert its vibration. Meanwhile, different numbers of balls does not lead the vibration characteristics change, compared to no SFD cases.Figure 7
The amplitude-frequency curves in VC frequency resonance region, where° indicates run-up and × indicates run-down.
(a)(b)(c)(d)(e)(f)In Figure8, the influence of the SFD on the 1/2 subharmonic VC frequency resonance is illustrated by the reduction or disappearance of the hysteresis areas. When Nb=13, the 1/2 subharmonic VC frequency resonance disappears. Basically, the SFD does not change the subharmonic resonance range of the original system. In addition, compared to [10], when Nb=8, the quasi-period with large amplitude near λ=6.78 disappears after adding the SFD. That means the 1/2 subharmonic resonance of the VC frequency can be eliminated by adjusting the number of balls. In case that the quasi-period exists, the SFD could improve the jump phenomenon in the subharmonic resonance region of the VC frequency, but its amplitude suppression effect is limited.Figure 8
The amplitude-frequency curves in subharmonic VC frequency resonance region, where° indicates run-up and × indicates run-down.


## 4. Parameter Analysis of the Ball Bearing with the Application of SFD

WhenNb=9, the amplitude-frequency curve of the disc in horizontal direction is shown in Figure 9(a). The second-order resonance shows the cross structure which means the coexistence of soft and hard stiffnesses. The spectrum of the run-up curve contains double frequency component, while the run-down curve has only periodic-1 response, as shown in Figures 9(b) and 9(c). The results are same as those when no SFD is mounted. Thus, the nonlinear characteristics are introduced by ball bearing support structure and do not been suppressed by the SFD.Figure 9
The amplitude-frequency curves and spectrums of responses of disc in horizontal direction whenNb=9. (a) The amplitude-frequency curves, (b) λ=2.87 (run-up), and (c) λ=2.87 (run-down) (blue line means run-up process and brown line means run-down process).
(a)(b)(c)It is noting that, in Figure9(a), there is a jump region between λ=5.51and 5.64. In the region, corresponding to λ=5.6, the frequency spectrum of run-up and run-down curves are shown in Figures 10(a) and 10(b). The adjacent region of λ=5.7 to 6.0 also merges jump behaviour. The spectrums of run-up and run-down curves at λ=5.85 are shown in Figures 10(c) and 10(d). Furthermore, taking the response at λ=5.69, the frequency spectrum shows periodic motion (Figure 10(e)). Therefore, it can be concluded that the two regions are independent to each other. The former is the 1/2 subharmonic vibration caused by the second-order natural frequency region. And, the latter is the 1/2 subharmonic vibration caused by the third-order main resonance of the system. From λ=6.0 to 6.3, period-1 and its 1/2 subharmonic components coexist. Then, the quasiperiodic motion merges, which is shown in Figure 10(f). Corresponding to the frequency λf=0.910, it locates on the first-order main resonance frequency region. Meanwhile, other quasi-periodic responses’ frequencies could be calculated as λf=2.754, λf=0.927, and λf=0.943, which correspond to the first and second natural frequencies of the system, respectively.Figure 10
Spectrums of the responses at (a)λ=5.6 (run-up), (b) λ=5.6 (run-down), (c) λ=5.85 (run-up), (d) λ=5.85 (run-down), (e) λ=5.69, and (f) λ=6.48 when Nb=9.
(a)(b)(c)(d)(e)(f)Figure11 shows the amplitude-frequency curve of the disc in horizontal direction when Nb=10. Its main characteristics, including crossover characteristic in the second-order resonance, the subharmonic resonance of the VC frequency, the subharmonic response of the third-order main frequency, and the quasi-periodic motion corresponding to the third-order main resonance, are the same as those of the case of Nb=9. However, in this case, the subharmonic resonance region corresponding to the second-order resonance has been merged into the subharmonic response region of the third-order main resonance. In addition, near λ=5.18, the frequencies of run-up and run-down curves are all quasi-periodic, which is not consistent with other cases. It can be seen that this is due to the coupling effect of the SFD and the bearing. The stability condition at this frequency is destroyed, which leads the instability interval expanding. The quasi-periodic regions and frequencies are λ=1.74 and λf=2.84, λ=3.39 and λf=0.952, λ=3.78 and λf=0.946, and λ=5.18 and λf=0.948.Figure 11
The amplitude-frequency curves of responses of disc in horizontal direction whenNb=10 (blue line means run-up process and brown line means run-down process).Figure12 shows the amplitude-frequency curves of the disc in horizontal direction with Nb varies from 11to 13. As the number of balls increases, the jump phenomenon of the second-order main resonance weakens and disappears, and the corresponding quasi-periodic area also decreases and disappears. In general, compared with the system with even number of balls, the quasi-period areas are less than those when the number of balls is odd. The quasiperiodic areas and corresponding frequencies of the system are listed in Table 2.Figure 12
The amplitude-frequency curves of responses of disc in horizontal direction when (a)Nb=11, (b) Nb=12, and (c) Nb=13 (blue line means run-up process and brown line means run-down process).
(a)(b)(c)Table 2
The quasi-periodic frequency with different numbers of balls.
Number of ballsQuasi-periodic frequencyNb=11(1)λ=3.48, λf=0.935;(2) λ=3.87, λf=0.968;(3) λ=6.44, λf=0.943;(4) λ=6.9, λf=0.927Nb=12(1)λ=3.06, λf=0.971;(2) λ=3.57, λf=0.959;(3) λ=3.95, λf=0.964;(4) λ=4.4, λf=0.967;(5) λ=6.56, λf=3.043;(6) λ=6.66, λf=0.976;(7) λ=6.91, λf=0.928Nb=13(1)λ=3.08, λf=0.978;(2) λ=3.54, λf=0.972;(3) λ=4.07, λf=0.993;(4) λ=6.66, λf=0.976;(5) λ=7.15, λf=0.96For horizontal comparisons, Figure13 plots the amplitude-frequency curve of the left end journal in horizontal direction when Nb=9 to 13. The change process of the quasi-periodic motion frequencies with the increase of the number of balls can be clearly seen.Figure 13
The amplitude-frequency curves of the left end journal in horizontal direction when (a)Nb=9, (b) Nb=10, (c) Nb=11, (d) Nb=12, and (e) Nb=13 (blue line means run-up process and brown line means run-down process).
(a)(b)(c)(d)(e)


## 5. Discussion

Through the previous analysis, we can see that, in the area whose corresponding frequency is smaller than first-order natural frequency which is normalized to 1, as the number of balls increases, the VC frequency peak generally decreases, and the stability of the system with odd-numbered balls is slightly better than that with even-numbered balls. The subresonance of VC frequency corresponds to the first-order natural frequency which is suppressed. In this way, with the application of the SFD, the main resonance peaks of each order can be reduced. The peak and area of the quasi-period motion can also be reduced. In this process, the SFD and the bearing have a coupling effect in certain intervals, which will increase the instability region. Therefore, it is necessary to adjust the damper parameters to further control the vibration of the system because the outer ring of the bearing could offset against the effects of gravity and bearing forces. Meanwhile, to increase the damping by reducing the gap is likely to cause rubbing at the journal. Therefore, it is much more safe to adjust the length or viscosity.Figure14 shows the amplitude-frequency curves of the disc and the left journal in horizontal direction when Nb=13 and L=16 mm. The first-order main resonance peak no longer has the phenomenon of hysteresis jump. Comparing to Figure 12(c), it can be seen that the quasi-periodic motion is suppressed, and the amplitude-frequency response forms a wide range between the first-order main resonance region and the third-order main resonance region. For the journal, only in the second-order main resonance region can there be an obvious peak. At the same time, it should also be noted that the quasi-periodic motion of the third-order main resonance has not been eliminated but relatively suppressed. In general, the SFD has a more obvious suppression effect on the main resonance peak and quasi-periodic motion than other frequency components and can be used together with the bearing as means to reduce vibration and extend the life of the bearing.Figure 14
The amplitude-frequency curves of responses of (a) disc and (b) left-end journal in horizontal direction whenNb=13 (blue line means run-up process and brown line means run-down process).
(a)(b)


## 6. Conclusions

This paper studies the nonlinear response characteristics of an asymmetric rolling bearing supported rotor system with application of SFD in different resonance regions. The following conclusions can be drawn:(1)
In first three orders’ prime resonance regions, the SFD suppresses the vibration significantly. Due to the nonlinearity introduced by the ball bearing, quasi-periodic solutions and jump behaviours still exist but with small amplitude.(2)
In VC frequency resonance region, the SFD almost has no damping effect because the amplitude of the response in VC frequency resonance case is small when no SFD is implemented. However, with the SFD, the soft and hard springs’ coexisting characteristic would disappear when the number of balls increases.(3)
When the number of balls increases, the nonlinear responses, such as quasiperiodic motion and double period motion, still exist. The SFD suppresses the vibration of nonlinear components generally but very little. However, with the SFD, the nonlinear behaviours and system stability would be affected by the number of balls.All these characteristics can help us understand the dynamic behaviours of the ball bearing supported rotor systems with the application of the SFD, based on which the design basis of such kind of system can be given.



---

*Source: 4370981-2021-09-21.xml*