In a dual-rotor system, speed ratio is a basic and important parameter during the operation process, and it usually has a great influence on the dynamic response of dual-rotor systems. The speed ratio
r is defined as the ratio of outer rotor speed
ωo to inner rotor speed
ωi, namely,
r =
ωo/
ωi, and
r > 0 indicates a co-rotation, while
r < 0 indicates a counter-rotation. Generally, the values of speed ratio adopted in dual-rotor systems reported in literatures [
2,
8,
11,
12] change from 1 to 2 for the typical case of co-rotation. Hence, the values of speed ratio
r fixed at 1.2, 1.3, 1.5 and 2.0 are selected to investigate their effects on MSDS responses. The rotating speed of the inner rotor is varied from 5 rad/s to 3000 rad/s at intervals of Δ
ωi = 5 rad/s. For convenient analysis,
ωi is normalized to the first order critical speed
ωc1i of MSDS excited by the inner rotor, the response displacements
x and
y are normalized to the AMB nominal air-gap length
g0, and the response frequency
ω is normalized to
ωi, namely
Ω =
ωi/
ωc1i,
X =
x/
g0 and
Y =
y/
g0, and
N =
ω/
ωi, respectively. For example, in the case of
r = 1.2, control parameters
Kp = 5 and
Kd = 0.0003, the first three order critical speeds calculated by complex Riccati transfer matrix method [
13] are provided in
Table 2. It is clear from
Table 2 that
r = 1.2 is in co-rotation with
ωc1i = 828.48 rad/s. Thus, a step of
Ω is Δ
Ω = Δ
ωi/
ωc1i ≈ 0.006. Furthermore, the simulation results show that the vibration characteristics of inner and outer rotors are similar, so the analysis is focused on the dynamic response of disk 1 on the inner rotor.
Bifurcation diagrams, dynamic trajectories, Poincaré maps and spectrum graphs are employed to exhibit the results of numerical simulation. Dynamic trajectories represent the instantaneous positions of the rotor center. Spectrum graphs exhibit the frequency contents of the rotor response and are obtained from the Fourier transformation of the time series of rotor response in the vertical direction. Poincaré maps are determined by sampling the trajectory of the rotor center at a constant interval of the forcing period of T = 2π/ωi and projecting the outcome on the X(nT) versus Y(nT) plane. Bifurcation diagrams demonstrate the trajectory of the rotor in the vertical direction against the variation of speed parameter Ω. In order to clearly demonstrate the variations of bifurcation with rotating speed, bifurcation diagrams are depicted with 3ΔΩ as the speed interval in this research. To significantly improve the computational efficiency of generating bifurcation diagrams, a parallel computing strategy is adopted in MATLAB programming and is implemented at the workstation.
3.1. Nonlinear Dynamic Responses
In the case of
r = 1.2, the types of bifurcation for the inner rotor response without rubbing are examined with the increase of speed parameter
Ω, and the corresponding bifurcation diagrams are generated in
Figure 7. After careful inspection, each speed corresponds to five points, which indicates that the MSDS always performs the sub-synchronous vibration of period 5 over the full speed range. It should be noted that for some speed ranges, the number of points displayed is less than five due to multiple close response amplitude points at the speed shown in a larger scale figure, as shown in
Figure 8a,c. In some intervals, the number of points for some rotational speeds is also less than five, while the corresponding Poincaré diagram shows that the number of points is still five. The reason is that displacements in the vertical direction that are intercepted by the bifurcation diagram at this moment are almost equal, as illustrated in
Figure 8b. It should be noted that the dynamic trajectories in
Figure 8, “○” and “△”, respectively, represent the positions of rotor center at a certain moment and the subsequent moment, according to which the change of the rotor whirling trajectory with time is demonstrated and the whirling direction of the rotor can be judged.
In the case of
r = 1.2 without rubbing, the vibration of MSDS is mainly dominated by the dual-frequency unbalanced excitation of inner and outer rotors. More specifically, it is determined by the speed ratio and rotational frequency difference of inner and outer rotors. Taking
Figure 8a as an example, there are two dimensionless frequency components of 1 and 1.2 in the spectrum, namely the dimensionless rotating frequencies of inner and outer rotors. In
Figure 8b, the dimensionless rotating frequency difference 1/5 and its multiple frequency components also appear in the spectrum. The most likely cause of this phenomenon is the operation of MSDS near the first-order critical speed. As a result, the MSDS is dominated by the sub-synchronous vibration of period 5 without rubbing when the speed ratio is
r = 1.2.
In the case of
r = 1.2 with rubbing, the bifurcation diagrams are provided in
Figure 9. For the range 0 <
Ω ≤ 0.978, the rotor response is the period 5 vibration. For
Ω = 0.984 and 0.990 ≤
Ω ≤ 1.002, the responses of period 20 and period 15 are observed, respectively. With further increase of
Ω, the sub-synchronous vibration of period 5, quasi-periodic and chaotic vibration are observed alternately. Period 5 vibration occurs in the ranges
Ω = 1.008, 1.074 ≤
Ω ≤ 1.110 and 1.123 ≤
Ω ≤ 1.129. Quasi-periodic vibration is observed to exist in the ranges 1.020 ≤
Ω ≤ 1.038,
Ω = 1.116, 1.153 ≤
Ω ≤ 1.159 and 1.183 ≤
Ω ≤ 1.195. Chaotic vibration is seen in the ranges
Ω = 1.014, 1.044 ≤
Ω ≤ 1.068, 1.165 ≤
Ω ≤ 1.177 and 1.201 ≤
Ω ≤ 1.304. Period 5 vibration is largely seen to dominate the rotor response for the range 1.310 ≤
Ω ≤ 3.6, except at
Ω = 3.253 where chaotic vibration is observed. By comparing
Figure 7 and
Figure 9, it is found that rubbing will cause stable multiple periodic motion, critical quasi-periodic motion and unstable chaotic motion.
In order to have a deep insight into the above abundant nonlinear vibrations stimulated by rubbing, the dynamic responses of period 5, 15, 20, quasi-periodic and chaotic vibration at
Ω = 1.008,
Ω = 0.996,
Ω = 0.984,
Ω = 1.032 and
Ω = 1.044 are depicted in
Figure 10. It is easy to observe that the period
n motion is characterized by
n points in Poincaré maps, and the minimum dimensionless frequency component in the spectrum is 1/
n. In addition, the response trajectory becomes more complicated with larger
n, as shown in
Figure 10a–c. When ghosting appears in the trajectory in
Figure 10d, the Poincaré map is exhibited as five closed elliptic curves, and the frequency components with clear contours can still be observed in the spectrum, which infers that MSDS performs a quasi-periodic motion. Furthermore, chaotic motion occurs when the trajectory becomes chaotic, many irregular points exist in the Poincaré diagram, and dense burrs appear in the spectrum, as depicted in
Figure 10e.
In the above nonlinear dynamic behaviors in co-rotation with
r = 1.2, the high-order sub-synchronous motion of period 15 and period 20 are multiples of five, corresponding to the rotating frequency difference 1/5 with speed ratio
r = 1.2. Under the combined actions of exponential order rubbing force, nonlinear electromagnetic force and unbalanced force, the fractional frequency components 1/15 and 1/20 are excited, which corresponds to 1/3 and 1/4 of the rotational frequency difference 1/5, respectively. Accordingly, in the spectrums in
Figure 10b,c, the frequency components can be regularly expressed as multiple frequency components of 1/15 and 1/20, respectively. Even for the quasi-periodic and chaotic motions, the main frequency components in the spectrums can be expressed as a series of combined frequency components of the rotating frequency difference 1/5 and its fractional frequency component 1/25, as shown in
Figure 10d,e. Therefore, the nonlinear vibration characteristics of MSDS are closely related to the speed ratio.
With parameters unchanged in counter-rotation
r = −1.2, the bifurcation diagrams of MSDS are shown in
Figure 11. For 0 <
Ω ≤ 0.972, the rotor response is period 5. The attractors of period 10, 20, and 15 are, respectively, observed for
Ω = 0.978,
Ω = 0.984 and
Ω = 0.996, and the rotor response is chaos at
Ω = 0.990. With further increase of
Ω, the alternation of chaotic vibration and sub-synchronous vibration of period 5 is captured, except for
Ω = 1.044 and 1.889 where period 10 vibration is observed, for
Ω = 1.919 where period 15 vibration is seen, and for
Ω = 1.183 where the response is quasi-periodic. Chaotic vibration is seen to exist in the ranges 1.002 ≤
Ω ≤ 1.038, 1.159 ≤
Ω ≤ 1.177, 1.225 ≤
Ω ≤ 1.267, 1.316 ≤
Ω ≤ 1.340,
Ω = 1.786 and 1.805 ≤
Ω ≤ 1.835. The existence of period 5 vibration is found in the ranges 1.050 ≤
Ω ≤ 1.153, 1.189 ≤
Ω ≤ 1.219, 1.273 ≤
Ω ≤ 1.310, 1.346 ≤
Ω ≤ 1.780, 1.792 ≤
Ω ≤ 1.799 and 1.841 ≤
Ω ≤ 1.925. Quasi-periodic vibration is found in the ranges 1.931 ≤
Ω ≤ 1.961 and 2.022 ≤
Ω ≤ 2.046. For the ranges 1.968 ≤
Ω ≤ 1.980, 2.004 ≤
Ω ≤ 2.016 and specific speed
Ω = 1.986, period 5, 20 and 10 vibrations are observed. The response is chaotic in the ranges 1.992 ≤
Ω ≤ 1.998 and 2.052 ≤
Ω ≤ 2.094. Period 5 vibration dominates the response for the range 2.100 ≤
Ω ≤ 3.6, except for chaos vibration in the small range 3.313 ≤
Ω ≤ 3.319.
As typical examples of nonlinear dynamic behaviors in counter-rotation with
r = −1.2, the sub-synchronous vibrations of period 10, 15, 20, quasi-periodic and chaotic vibration are illustrated in
Figure 12 at
Ω = 0.978,
Ω = 0.996,
Ω = 0.984,
Ω = 1.183 and
Ω = 2.052, respectively. In the case of rubbing, the dynamic response characteristics for counter-rotation with
r = −1.2 are almost the same as those of co-rotation with
r = 1.2. When MSDS performs sub-synchronous vibrations of period
n with
n = 10, 15 and 20, its frequency components can be expressed as multiple times of the fractional frequency 1/
n corresponding to the rotating frequency difference 1/5. For the quasi-periodic motion and chaotic motion, the main frequency components can be expressed as various combined frequency components of rotating frequency difference 1/5 and its fractional frequency. Therefore, the nonlinear dynamic characteristics for counter-rotation with rubbing are also determined by the speed ratio. In addition, it can be observed that the rotor alternately performs forward whirl and backward whirl in the dynamic trajectories in
Figure 12a–d, while the rotor only performs forward whirl in the dynamic trajectories in
Figure 10a–d. It implies that the dynamic response of MSDS in counter-rotation is more complicated than that in co-rotation.
From the above analysis, it is obvious that the sub-synchronous vibration of period 5 is the fundamental motion under fixed-point rubbing with r = ±1.2. In contrast to co-rotation r = 1.2, the dynamic responses in counter-rotation r = −1.2 are more complicated due to the opposite direction of the gyroscopic moment. More specifically, period n (n = 5m, m = 2, 3, 4,...), quasi-periodic and chaotic vibrations span wider speed ranges and occur more frequently in counter-rotation, which is also tenable for other speed ratios.
3.2. Influence of Speed Ratio
To investigate the effects of speed ratio on the bifurcation characteristics of MSDS, the bifurcation diagrams and nonlinear dynamic responses for r = 1.3, r = 1.5 and r = 2.0 in co-rotation are investigated in this section.
When the speed ratio is
r = 1.3, the bifurcation diagrams of MSDS are illustrated in
Figure 13. By careful examination, sub-synchronous motion of period 10 is found to be the fundamental motion that dominates the dynamic responses.
n = 10 seems to correspond to the denominator of the rotational frequency difference 3/10. In addition, sub-synchronous motion of period 20, quasi-periodic and chaotic motion mainly occur in the ranges 1.104 ≤
Ω ≤ 1.593 and
Ω = 3.006. The dynamic responses of various typical motions are shown in
Figure 14. Similarly, for the sub-synchronous motion of period 10 and period 20 in
Figure 14a,b, the basic frequency components are 1/10 and 1/20, respectively, and the other frequency components can be expressed as their multiple frequency components. In quasi-periodic motion with at 10 closed curves in
Figure 14c, the main frequency component can be expressed as the multiples of basic frequency component 1/10.
The bifurcation diagrams and corresponding typical dynamic responses of MSDS with speed ratio
r = 1.5 are provided in
Figure 15 and
Figure 16, respectively. In this situation, the period 2 vibration dominates the whole speed range, where
n = 2 is equal to the denominator of the rotational frequency difference ½. Quasi-periodic motion and chaotic motion exist in the ranges 1.129 ≤
Ω ≤ 1.255 and 2.384 ≤
Ω ≤ 2.619. A double period 4 motion at
Ω = 0.845, multiple period 28 vibration at
Ω = 1.213 and the quasi-period motion with two closed curves at
Ω = 1.044 are demonstrated as typical examples of dynamic responses in
Figure 16. The sub-synchronous motion of period
n is characterized by
n points in the Poincaré maps and the basic frequency component 1/
n. The main frequency components of quasi-period motion may be expressed as the combination of rotational frequency difference 1/2 and its fractional frequency 1/10, for example, 1/5 = 2 × 1/10, 3/10 = 1/2 − 2 × 1/10, and so on.
The bifurcation diagrams and corresponding typical dynamic responses of MSDS with speed ratio
r = 2 are depicted in
Figure 17 and
Figure 18. It is obvious that the simplest form of motion, sub-synchronous motion of period 1, namely synchronous motion, occupies an absolute advantage over other forms of motion. Also, the 1 in period 1 corresponds to the rotating frequency difference 1 in the case of
r = 2. The sub-synchronous motion of period 2, quasi-periodic motion and chaotic motion mainly appear in the range of 1.165 ≤
Ω ≤ 1.219, 1.986 ≤
Ω ≤ 2.456, 2.837 ≤
Ω ≤ 2.861, respectively, and their spectral characteristics are similar to those of other speed ratios.
It is shown from the above analysis that the nonlinear vibration characteristics of MSDS with fixed-point rubbing mainly depend on the speed ratio. First, the most basic sub-synchronous vibration of period
n is determined by the speed ratio
r. The relationship between
n and
r is that
n equals the greatest common divisor (GCD) of 10 times 1 and
r divided by 10, i.e.,
n = 10/GCD(10 × [1·
r]), which is listed in
Table 3. The minimum dimensionless frequency component is 1/
n for the basic motion in the spectrum and other frequency components can be expressed as integer multiples of 1/
n. These two laws are also applicable to the higher order sub-synchronous motion of period
mn (
m = 2, 3, 4,...). Second, quasi-periodic motion with speed ratio
r is characterized by
mn closed curves, and its main frequency components can be expressed as a variety of the combined frequency components of rotating frequency difference
r-1 and its fractional frequency, which is also suitable for chaotic motion with the same speed ratio. It should be pointed out that the above characteristics are also verified by the bifurcation diagrams of other speed ratios, which are neglected for simplicity. Consequently, it clearly states that speed ratio has considerable influence on the nonlinear dynamic characteristics of MSDS. Actually, the influence of speed ratio is achieved by the coupling effects of inter-shaft bearings and plentiful nonlinear dynamic behaviors reflect the strong cross-excitation phenomenon between inner and outer rotors, which is one of the most significant differences between dual-rotor systems and single-rotor systems.
Another interesting phenomenon is the speed ranges where complicated dynamic responses exist. Except for the most basic motion for each speed ratio, high-order sub-synchronous motions, quasi-periodic and chaotic motions are generally excited near the critical speeds of MSDS. Taking
r = 1.2 in co-rotation in
Figure 9 for an example, the sub-synchronous vibrations, except for period 5, mainly exist in the range 0.984 ≤
Ω ≤ 1.304, which includes the first two order critical speeds excited by the inner rotor, namely Ω
c1i = 1 and
Ωc2i = 1.238. For the chaotic vibration at
Ω = 3.253, it is close to the inner rotor speed
Ω3i =
Ωc3o/
r = 3.306, corresponding to the third order critical speed
Ωc3o = 3.968 excited by the outer rotor. It indicates that the motions near critical speeds are much more complicated than those away from them. The main reason is due to the fact that the vibration amplitude is larger and the rubbing response is more intense when the MSDS operates near the critical speeds, so the friction force during the rub impact process increases. Since the direction of the friction force is always opposite to the rotation direction of rotors, it will drive the rotor to whirl backward when its amplitude increases to a certain value, thus resulting in more complicated nonlinear vibrations.
3.3. Parameter Selection
During the operation of rotating machinery, prolonged exposure to fluctuating stresses caused by non-synchronous vibrations may result in the fatigue failure of rotor components [
20,
21,
22]. In general, the alternating frequency of fluctuating stress increases with the complexity of asynchronous motion, and high alternating frequency will accelerate fatigue failure. Therefore, non-synchronous vibrations, especially high-order sub-synchronous motion, quasi-periodic motion and chaotic motion, are undesirable.
It is found from the above analysis that the rotor motion is synchronous only when the speed ratio is an integer. However, the speed ratio is usually not an integer, so the MSDS will inevitably undergo non-synchronous motion even without rubbing. From the perspective of reducing the severity of fluctuating stress, the complexity of nonlinear dynamic response should be reduced if possible. In other words, the speed ratio corresponding to smaller n is recommended to be selected.
Another more effective method is to take full advantage of the adjustable and controllable support characteristics of MSDS. With the real time adjustment of the support characteristics by the control system, the equivalent damping of AMBs can be increased to attenuate response amplitudes, and the equivalent stiffness can be changed to keep critical speeds away from the operating speeds of MSDS [
13,
14,
15], which will bring the MSDS into fundamental sub-synchronous vibrations. For example, with a larger derivative coefficient of
Kd = 0.0005 and other parameters remaining the same, the bifurcation diagrams of rotor response for
r =1.2 is depicted in
Figure 19. Comparing
Figure 9 and
Figure 19, it is obvious that except for the basic sub-synchronous vibration of period 5, other types of sub-synchronous vibrations disappear for larger
Kd = 0.0005, especially for the range 0.984 ≤
Ω ≤ 1.304. It is exactly where the MSDS is superior to the traditional dual-rotor system. Nevertheless, the chaos vibration at
Ω = 3.254, corresponding to the first order bending modes, still exists for the reason that the derivative coefficient has little influence on this mode of the system [
13], which deserves further investigation.
In order to have an insight into the influence of AMB control parameters on MSDS nonlinear dynamics, the bifurcation diagrams of rotor response with proportional coefficient
Kp and derivative coefficient
Kd are shown in
Figure 20, where
ωi = 700rad/s,
r = 1.2, and other parameters remained the same with
Table 1. In
Figure 20a, MSDS performs the simplest sub-synchronous motion of period 5 when
Kp is in the range 1.8 ≤
Kp ≤ 39.4 and it alternately performs high-order sub-synchronous motion, quasi-periodic motion and chaotic motion when
Kp is beyond the range 1.8 ≤
Kp ≤ 39.4. The reason is that increasing
Kp will increase the equivalent stiffness of AMB, and too large or too small equivalent stiffness will make the operating speed close to a certain order critical speed of MSDS. In
Figure 20b, MSDS performs the simplest sub-synchronous motion of period 5 when
Kp is in the range 0.0002 ≤
Kd ≤ 0.00086, and it performs complicated non-synchronous motions when
Kd is small or large. The reason is that the AMB equivalent damping increases with
Kd; when equivalent damping is too small, the capacity of MSDS to attenuate vibration reduces, and when it is too large, it directly amplifies the influence of noise in differential signal on the system stability, resulting in unstable motions [
23]. Therefore, it is necessary to select appropriate control parameters to make the MSDS perform the simplest sub-synchronous motion.
In summary, there are three suggestions for parameter selection to keep the asynchronous vibrations of MSDS as simple as possible. First, the speed ratio r corresponding to smaller n is recommended, where n = 10/GCD(10 × [1·r]), and GCD represents the function of the greatest common divisor. Secondly, according to the MSDS working speed and the influence of proportional coefficient on critical speeds of the system, appropriate proportional coefficient is suggested to keep critical speeds away from working speeds. Additionally, selecting a larger derivative coefficient to increase the equivalent damping of AMBs can effectively attenuate the vibration of MSDS, but it should be noted that when equivalent damping is too large, it will increase the instability risk by directly amplifying the noise in the derivative signal. Thirdly, the designed operating speed of MSDS should be far away from critical speeds corresponding to the bending modes.