Abstract

The dynamic response and its stability of a blade rotor with coupled rubbing in the labyrinth seal and tip seal are investigated. The dynamic equations are established based on the Hertz contact rubbing force at the labyrinth seal and the tip rubbing force considering both the contact deformation of the tip seal and the bending deformation of the blade. Numerical simulations show that the coupled rubbing response includes periodic motions, almost periodic motions, and chaotic motions. Compared with the single rubbing fault, coupled rubbing increases the range of rotation velocity of contact. A new continuation shooting method is used in the solution and stability analysis of the periodic response to give the bifurcation diagrams. The paths of the system for entering and exiting chaos are analyzed.

1. Introduction

To improve the operational efficiency of the turbo machinery, the clearance between the rotor and stator is designed to be smaller and smaller, which also increases the possibility of rubbing. Severe rubbing may cause breakage of the blade and instability of the rotor, which eventually lead to serious operational accidents [1].

There are usually two kinds of seals for turbo machines: labyrinth seals and tip seals. The rubbing of labyrinth seals is usually described as a rubbing problem between two cylindrical surfaces. The contact forces may be expressed as a linear or nonlinear function of the deformation of interaction. The friction forces are given by the Coulomb law of friction. Many researchers focused on related problems, such as forward periodic motion [24], reverse precession [5, 6], chaos [7, 8], and transient response [9]. Other interesting rubbing problems have also attracted people's attention. Patel [10, 11] investigated a rotor with coupled faults of fatigue crack and rotor-stator rub. Chang-Jian et al. [12] focused on the response of a flexible rub-impact rotor supported by oil film bearings. Hou et al. [13, 14] presented the effects of maneuver flight on the rotor/seal rubbing response. Cao et al. [15] investigated a rub-impact rotor system with fractional order damping. Yang et al. [16] studied the dynamic behaviors of a rotor-casing system subjected to axial load and radial rub. Sun et al. [17] considered a dual-rotor system induced by impact rub. Tai et al. [18] studied the periodic response of the rotor system with single-point impact rub. Roques et al. [19] investigated speed transients with rotor-to-stator rubbing caused by an accidental blade-off imbalance. Jin et al. [20] studied a dual-rotor-bearing-coupling misalignment system with blade-casing rubbing numerically and experimentally to reveal the nonlinear vibration characteristics of an actual aero-engine rotor system. Wang et al. [21] considered the stick-slip transitions, which are related to the dynamical scenarios of sliding bifurcations, of the self-excited dry friction backward whirls of a general rotor/stator rubbing system. Hildebrandt et al. [22] discussed the heat flux distribution problem between bristle package and rotor during brush seal rubbing event based on a 3D CFD model.

When rubbing happens at the tip seal, the number of blades in contact and the deformation of each blade is different at each moment, which is more likely to lead to a complex dynamic response. For the tip rubbing, by assuming the blade as a cantilever beam, Choy [9] and Padovan [23] derived the relationship between the normal contact force and the radial deformation of the blade and investigated the nonlinear dynamics of rotor/casing rub interactions. Jiang et al. [24] proposed a rubbing contact force model in the form of periodic pulse force and derived the expression of the maximum contact force through theoretical calculation. Lesaffre et al. [25] studied the blade-machine impact friction stability, and the rubbing was achieved by using an additional stiffness matrix and damping matrix. Padova [26] conducted an experimental study to analyze the interaction between the blade and casing with different incursion depths. Kim [27] studied the rubbing of the asymmetric blade rotor and the anisotropic casing. Sinha [28] investigated the dynamics of an asymmetric turbofan rotor with fan blade-loss event. The fan blades are modeled as pretwisted thin shallow shells. Chen [29] carried out numerical simulation and experimental research on the rotor dynamics of the blade tip rubbing, in which the frictional force of the blade and the casing adopts the linear rubbing model. Ma et al. [30] studied a rotor system with fixed-point rubbing fault. They established a point-point contact model to simulate the rotor-stator rubbing. The augmented Lagrangian method is applied to deal with contact constraint conditions. Ma et al. [31] investigated blade disk dovetail structure under blade tip rubbing based on ANSYS software and discussed the effects of rubbing under different rotating speeds and penetration depths. The local rubbing phenomenon was modeled by a pulse force. The conditions of no rubbing, single-blade rubbing, and four-blade rubbing are considered in the numerical simulations. Yang et al. [32] analyzed the vibration of a dual-rotor-bearing-double casing system with pedestal looseness and multistage turbine blade-casing rub considering the distribution of soft and hard coatings.

Since the clearance of the labyrinth seal and the tip seal is very small, when the blade rotor is excessively vibrating, rubbing may occur at both two kinds of seals. Coupled rubbing may result in a more complex dynamic response. Thus, it is necessary to investigate the dynamic response of the blade rotor with coupled rubbing. Furthermore, the existing research mainly adopted numerical simulation and experimental methods. Few studies focused on the stability and bifurcation of periodic motion.

In this paper, we propose a study on the coupled rubbing response of the blade rotor and the bifurcation behavior of the rubbing periodic solutions. The structure of this paper is as follows. In Section 2, dynamic equations of coupling rubbing of the blade rotor are established. The complex dynamic behaviors appearing at the rotor's rising up and slowing down are analyzed by numerical simulation. In Section 3, a new continuation shooting method is introduced. And it is employed to solve the periodic response of the coupled rubbing and analyze the stability of the solutions. The influence of rotation velocity on the dynamic response of the coupled rubbing is analyzed. Then, the paths of entering and exiting chaos of the system are discussed. Finally, the conclusion of this paper is summarized.

2. Numerical Simulation of the Blade Rotor System with Coupled Rubbing

2.1. Dynamical Equations of the Blade Rotor

Figure 1 is the model of the blade rotor with coupled rubbing. The rotor is rigidly supported, and the mass of the rotor is concentrated at the disks. The shaft of the rotor is a mass less Euler–Bernoulli beam, and the bending stiffness is . and are the lumped masses of the labyrinth disk and the blade disk, respectively. Since the diameter of the labyrinth disk is small, the moment of inertia can be negligible; thus, . and are the equivalent moment of inertia and the polar moment of inertia of the blade disk, respectively. The equation of motion of the system can be given aswhere are the lateral displacements of the labyrinth disk, and are the lateral displacements and the deflection angles of the blade disk, respectively. We assume that the eccentricity of the rotor only exists at the blade disk. is the eccentricity, is the angular velocity of the rotor, and is the rotation velocity. are the stiffness coefficients, and their expressions can be given as

The damping is supposed to be proportional to the stiffness; there are , where is the proportional damping coefficient. are horizontal and vertical components of rubbing forces of the labyrinth seal, respectively. And their expressions can be given as follows:where the contact force is given by Hertz contact theory (see Figure 2), where , , are Young's modulus and the Poisson ratio of the rotor and the seal, respectively, is the radius of the labyrinth seal, and is the clearance between the seal and disk, and the friction force is .

2.2. The Model of Tip Rubbing Forces

For the tip rubbing, we consider two aspects of deformation: first is the contact deformation of the tip seal , which is caused by the contact of blade tip and casing and can be described by the Hertz contact model, and the second is the radial deformation of the blade , which is caused by the bending of the blade (shown in Figure 3(a)) under the action of rubbing forces. The blade is regarded as an Euler beam, and is the contact force of a single blade. The friction force is set to obey the Coulomb friction law, and there is . are the span length of the blade and the radius of the disk, respectively. is the span coordinate of the blade. The bending moment equation of the blade can be given aswhere is the deflection of the blade. are the density, cross-sectional area, and bending stiffness of the blade, respectively. Considering the boundary conditions , the deflection function is assumed to be a high-order polynomial aswhere are undetermined coefficients. Substituting (5) into (3), the expression of can be obtained. Figure 4 shows the deflection curves of the beam with .. The calculation results of are essentially the same. Hence, we set . The radial deformation of the blade tip can be obtained as

And based on Hertz contact theory, the deformation of the seal can be given aswherewhere are Young's modulus and the Poisson ratio of the blade and seal at the tip seal, respectively, is the radius of the seal, is the tip clearance, and there is . We can find the relationship between the total deformation and the contact force as

Hence, the contact force can be written as

Figure 3(b) shows the diagrammatic sketch of tip rubbing of the blade rotor. This paper ignores the influence of the initial position angle of the blade. Assuming that the position angle of the first blade is 0, then the coordinates of the first blade tip arewhere is the position angle of th blade and is the number of blades. Then, we can find the radius of the tip of th blade as

If , the blade will contact with the seal. Considering (9), the rubbing forces of th blade can be written asand if ,

Projecting the rubbing forces of ith blade in the and directions, then there are

By summing up the rubbing forces of all blades, the tip rubbing forces of the rotor can be given as

2.3. Numerical Simulation of the Rubbing Response of Blade Rotor System

By setting and , (1) can be written in a nondimensional form. Then, the Runge–Kutta method is used to obtain the dynamic response. The parameters of the system are given as follows (Table 1).

Figure 5 shows the bifurcation diagrams and Lyapunov exponents of the rubbing response, where Figures 5(a) and 5(b) are the bifurcation diagrams considering only the labyrinth seal rubbing and the tip rubbing, respectively. The rubbing response of the labyrinth seal is dominated by the period 1 motion. The tip rubbing response is dominated by a chaotic response. And the coupled rubbing response includes both chaotic response and complex periodic response. At the rotor's rising up, the range of rotation velocity of coupled rubbing, in which contact occurs, is greater than single rubbing fault.

Furthermore, for the coupled rubbing, the contact of the labyrinth seal can maintain at the whole range of rotation velocity of rubbing, and the contact of tip seal ends at . Figure 6 shows the partial enlargement of the bifurcation diagram and Lyapunov exponents of coupled rubbing. Figure 7 shows the typical periodic coupled rubbing responses, and Figure 8 gives the aperiodic coupled rubbing responses. It can be seen that the response of the system includes period 1–8, period 12, and period 16 motions, almost periodic motions, and chaotic motions. The transformation between various types of response is also very complicated. Thus, it is necessary to study the change of response with rotation velocity, especially the path for the periodic response to chaos.

3. Stability Analysis of the Periodic Solutions of the Coupled Rubbing

In this section, the stability of the coupled rubbing response of the blade rotor is studied by the continuation shooting method. The principle of the continuation shooting method is shown in Figure 9. From a regular point of the solution curve on the Poincare map, we predict the next point in the tangent direction of the solution curve and then make the correction by the Newton–Raphson method on the hyperplane, which is orthogonal to the tangent direction, to get the next regular point of the solution curve. The calculation process is as follows.

Consider the following nonlinear differential equations:

It can be rewritten in the form of nonlinear algebraic equations as follows:where , is the Poincare map, is the intersection of the periodic solution and Poincare map corresponding to parameter , and is the period of the system. Deriving as a regular point of the solution curve, the next prediction point iswhere is the step. is the direction of prediction, and it can be determined bywhere and is dimensional zero vector. However, it is difficult to solve the value of . The Secant line can be used instead of the tangent as the prediction direction in practical calculations. Assuming and are two adjacent regular points, then

The correction process is carried out by the following formula:where and . can be solved by the following initial value problem:

When , . can be given aswhere is a small perturbation. Iterate from until it is satisfied , and set . For this method, obtained in the last iteration is the Floquet multiplier matrix of the periodic solution, and the stability of the periodic solution can be determined by its eigenvalues.

Combining the numerical simulation with the results of periodic response calculation and the stability analysis, Figure 10 shows the bifurcation diagram of the periodic response of blade rotor with coupled rubbing. The response is period 1 motion at the beginning of rubbing. As the rotation velocity increases, period-doubling bifurcation occurs at point a and the response becomes period 2 motion. Then, the response enters almost periodic motion through Hopf bifurcation at point b. The response breaks into chaos through almost periodic torus. The response returns to the period 2 motion at the c point through Hopf bifurcation. Further, the period 2 motion enters chaos through the period-doubling bifurcation near the point d, and the chaotic motion changes into the periodic motion at about . Then, the response goes back to period 2 motion (shown in Figures 10(a) and 10(b)). The period 2 motion enters chaos again through the Hopf bifurcation at point e and then leaves chaos by the saddle node bifurcation at point f. The response of the system becomes period 6 motion, and the period 6 motion jumps into period 2 motion through a saddle node bifurcation at point . Then, period 2 motion enters chaos through the saddle node bifurcation at point h (shown in Figure 10(c)).

The chaotic response of the system changes into period 6 motion by an inverse period-doubling bifurcation near point i and then changes to period 2 motion by saddle node bifurcation at point j. Period 2 motion becomes period 1 motion again through saddle node bifurcation at point k (shown in Figure 10(d)). Point l and m are period-doubling bifurcation points, and the response between the two points is period 2 motion (shown in Figure 10(e)). With further increase in rotating velocity, the motion of period 1 enters chaos by saddle node bifurcation at point n. The chaotic response turns into period 6 motion through saddle node bifurcation at point o, and then period 6 motion becomes unstable through period-doubling bifurcation at p point and becomes stable again through period-doubling bifurcation at point q. After period-doubling inverse bifurcation at point r, the response of the system is period 3 motion (shown in Figures 10(f) and 10(g)). Period 3 motion enters chaos through intermittent at point s, then exits chaos at point t through Hopf bifurcation, and enters chaos again through intermittent at point u (shown in Figure 10(h)).

The response of the system moves into period 5 motion at the point through the saddle node bifurcation. The period 5 motion enters chaos through period-doubling bifurcation at point, exits chaos through period-doubling bifurcation at point x, and then enters chaos again through Hopf bifurcation at y point. It is worth noting that the path to chaos here is not almost periodic torus. Figure 11 is the Poincare map of the chaotic response at . As the fixed points of period 5 motion become unstable, almost periodic torus disappears, and finally the response is separated from the position near the fixed points and moves towards the chaotic attractor.

With further increase in rotation velocity, the chaotic response returns to period 5 motion through the period-doubling bifurcation (near z point) and then enters chaos through the saddle node bifurcation at point . The response of the system becomes period 1 motion through the Hopf bifurcation point and then jumps to period 3 motion through the saddle node bifurcation at the point . The motion of period 3 enters chaos at the point through the period-doubling bifurcation, then exits chaos through the period-doubling bifurcation at point , and enters chaos again through the saddle node bifurcation at point . At the point , the response of the system returns to period 1 motion through the Hopf bifurcation and then jumps through saddle node bifurcation at point , and the rotor and the seal are out of contact.

When the rotor slows down, the response of the system jumps up near and rubbing occurs. The following changes of the response are mostly the same as the process of rotor’s rising up, and the difference is in the range of rotational velocity of 5100–5140 rpm. When the rotation velocity is reduced, the response changes into period 2 motion through a period-doubling bifurcation at the point . Then, saddle node bifurcation occurs at the point , the period 2 response jumps to the upper and lower solution branches, and then the periodic response becomes unstable through the Hopf bifurcation at the point .

In addition, it is shown that there are coexistence of multiple solutions corresponding to some rotation velocities, such as the coexistence of periodic solutions (see Figures 10(d) and 10(k)), coexistence of periodic solution and almost periodic solution (see Figure 10(j)), and coexistence of periodic solution and chaos (see Figure 10(i)). If the system is disturbed in the process of work, the response may change into other motions.

4. Conclusions

In this paper, the dynamic response and its stability of a blade rotor with coupled rubbing of the labyrinth seal and the tip seal are investigated. Dynamic equations of the blade rotor with coupled rubbing are established. The rubbing force of the labyrinth seal adopts the Hertz contact model. The tip rubbing model considers the Hertz contact deformation of the tip seal and the bending deformation of the blade. The bending deformation of the blade is described more accurately by using a high-order polynomial. The dynamic equations of the blade rotor with coupled rubbing are numerically simulated. It is found that the labyrinth seal rubbing response is dominated by periodic motion, the tip rubbing response is dominated by chaos, and the coupled rubbing response includes both chaotic motion and complex periodic motions. Moreover, the range of rotational speed at which the coupled rubbing rotor contacts is greater than a single rubbing fault.

A new continuation shooting method is employed to calculate the periodic solutions and analyze their stability. Considering the results of numerical simulation, the bifurcation diagrams of the periodic response of coupled rubbing are given, and the paths of the system entering and exiting chaos are analyzed. It shows that the paths of the system entering chaos include period-doubling bifurcation and intermittent, and the paths exiting chaos include Hopf bifurcation, intermittent, and period-doubling bifurcation. In addition, the investigation finds the coexistence of multiple solutions corresponding to some rotational velocities, such as coexistence of periodic solutions, coexistence of periodic solution and almost periodic solution, and coexistence of periodic solution and chaos.

Data Availability

The data that support the findings of this study are available from the corresponding author (Huabiao Zhang) upon reasonable request.

Conflicts of Interest

The author(s) declare that there are no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Authors’ Contributions

Huabiao Zhang and Xinye Li contributed equally to this work.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (Grant nos. 11972145 and 11302223) and Tianjin Excellent Science and Technology Commissioner Support Project (Grant no. 20YDTPJC00080).