The Influence of Shaft’s Bending on the Coupling Vibration of a Flexible Blade-Rotor System

The influence of shaft bending on the coupling vibration of rotor-blades system is nonignorable. Therefore, this paper analyzed the influence of shaft bending on the coupling vibration of rotor-blades system.The vibration mode function of shaft under elastic supporting condition was also derived to ensure accuracy of the model as well. The influence of the number of blades, the position of disk, and the support stiffness of shaft on critical speed of system was analyzed. The numerical results show that there were two categories of coupling mode shapes which belong to a set where the blade’s first two modes predominate in the system: shaft-blade (SB) mode and interblade (BB) mode due to the coupling between blade and shaft. The BB mode was of repeated frequencies of (Nb−2) multiplicity for number blades, and the SBmode was of repeated frequencies of (2) multiplicity for number blades.What is more, with the increase of the number of blades, natural frequency of rotor was decreasing linearly, that of BB mode was constant, and that of SB mode was increasing linearly. Natural frequency of BB mode was not affected while that of rotor and SB mode was affected (changed symmetrically with the center of shaft) by the position of disk. In the end, vibration characteristics of coupling mode shapes were analyzed.


Introduction
In the compressor/turbine engineering equipment, the blade, the disk, and the shaft are assembled together by a certain connecting structure and forming a blade-disk-shaft coupling system which has some coupling characteristics.With the development of scientific research, many scholars in this field have carried out a lot of research on the vibration of blades rotor coupling system.However, the influence of shaft bending on the coupling vibration of rotor-blades system is nonignorable.And there are potential couplings between shaft bending and blade bending under the right condition in fact.Therefore, this present research focuses on the influence of shaft bending on the coupling vibration of rotor-blades system.
The flexible disk model is the most popular analytical model for such systems at first.Parker and Mote [1] predicted natural frequencies of stationary annular or circular plates with a perturbation solution, and they discussed the mode split of the repeated natural frequencies.Khorasany and Hutton [2] studied linear vibration behavior of annular spinning disk whose inner boundary is free to move in the axial direction while the body of the disk is constrained by a space fixed linear spring.Kim et al. [3] reported that significant changes could occur to its natural frequencies and modes when a structure deviates from axisymmetry because of circumferentially varying model features.Dopkin and Shoup [4] studied the influence of disk's flexibility on natural frequencies by a transfer matrix method.
As the rotating systems become lighter and more flexible with higher operating speed for higher productivity and economical design, the rigid shaft-disk model cannot accurately predict the vibration characteristics of the system.And the shaft-blades model is considered, considering the coupling vibration among shaft bending and blade bending.Anegawa et al. [5] successfully analyzed the coupling of the shaft with in-plane vibration of the blades with  = 1 modes in a 1-DOF mass-spring system.What is more, the existence of an unstable region due to coupling and the characteristics of coupled resonance were indicated which was confirmed by experimental data.Turhan and Bulut [6] studied coupling effects between shaft-torsional and bladebending (in-plane) vibrations in turbomachinery and other rotating bladed structures and considered an idealized model where the blades are represented by uniform Euler-Bernoulli beams.Besides, they derived a linearized mathematical model for the system by using synthetically, multiframe, and mixed (consisting in a coalescence of finite element and Galerkin methods) approach.Chiu and Chen [7] investigated the influence on coupling vibrations among shaft-torsion and blade-bending coupling vibrations of a multidisk rotor system analytically.They solved the natural frequencies and the mode shapes of the system for one-to three-disk cases as examples and drew the diagrams of the coupling mode shapes.From the results, it was found that the interblade (BB) modes were of repeated frequencies of (  −1) multiplicity for number of blades.What is more, the numerical calculations also revealed that the natural frequencies were affected by blade numbers and disk distance.Lesaffre et al. [8] presented a model of fully flexible bladed rotor developed in the rotating frame and took the gyroscopic effects as well as the spin softening effects and the centrifugal stiffening effects into consideration.Huang and Ho [9] introduced a new approach to analyzing the dynamic coupling between shaft torsion and blade bending of a rotating shaft-disk-blade unit.Numerical examples showed not only coupling between the shaft, disk, and blades, but also coupling between individual blades and the blade coupling modes occurred only in repeated frequencies.What is more, the effect of shaft speed on the modal frequencies was investigated.Plots illustrating the occurrence of critical speeds and flutter instabilities are presented.Chiu and Yang [10] investigated the influence on coupling vibrations among shaft-torsion, blade-bending, and lacing wire coupling vibrations of a multidisk rotor system with grouped blades analytically.Moreover, they also found that not only the natural frequencies in a mistuned blade but also the types of coupling modes were found changed.
Then the disk-blade model is used to predict the influence of stagger and twist on natural frequencies for stationary as well as rotating bladed disk systems [11].And the flexible disk-blades model attracts a lot of research since the early researches [12][13][14].To get a better understanding and an accurate prediction of the dynamic characteristics of shaftdisk-blade units, Chun and Lee [15] developed an analysis substructure synthesis and assumed modes method to investigate the effect of the flexibility of a bladed disk assembly on the vibrational modes of a flexible rotor system.Yang and Huang [16] studied the longitudinal motion effects on the coupled vibration in a shaft-disk-blades system.Yang and Huang [17][18][19] investigated the coupling vibrations and applications among support-longitudinal, shaft-torsion, disktransverse, and blade-bending vibrations in a shaft-diskblades unit.Chiu and Huang [20] investigated the influence on coupling vibrations among shaft torsion, disk-transverse, and blade bending of a rotor system with a mistuned blade length analytically and found that the shaft-disk-blades system appeared to have four types of coupling vibrations: shaft-blade (SB), shaft-disk-blade (SDB), disk-blade (DB), and blade-blade (BB) mode.
From the preceding research it could be seen that there were two kinds of analysis method for the typical structure.One was regarding the support of blade as clamped or elastic support and analyzing the vibration of the blade separately and the other was regarding blade-disk system as a rigid body component and analyzing the dynamic characteristics of the rotor system separately.However, the research on blade-disk-shaft coupling system had gradually become the focus of research.Although some scholars had researched the coupling vibration of shaft torsion and blade bending, the influence of shaft's bending is usually ignored.But we found that there is a new mode which is SB due to the coupling vibration of shaft's bending and blade's bending.And the present paper aim is to analyze the influence of shaft bending on the coupling vibration of shaft-disk-blades system.What is more, influences of the blades number, disk position, and support stiffness on the characteristic value of system, the Campbell diagram, and the mode shape of vibration were studied by numerical simulation.

The Establishment of Dynamic Model
A dynamic model of rotor-blade coupling system with elastic restraints is shown in Figure 1.The rotor system is composed of continuous flexible shaft and rigid disk.The Timoshenko beam model is used to derive the energy equation of shaft and the bending and torsional vibration are also taken into consideration.The disk is regarded as a mass point.Considering the radial and transverse vibration of blade, the energy equation of blade is derived by using the Euler-Bernoulli beam model.- is the global coordinate system.  -      is the local coordinate system of disk where the geometric center of disk is the origin.  -      is the local coordinate system of blade where the geometric center of blade root is the origin, and  point is an arbitrary point on the blade.As shown in the figure, both ends of the rotor are elastic supported, where the elastic support stiffness and damping at both ends are assumed to be equal.

The Establishment of the Energy Equation of Flexible
Rotor Model.According to the relevant knowledge of elastic mechanics, the kinetic energy of the shaft can be expressed as where   and   are cross-sectional moment of inertia (transverse moment of inertia) and polar moment of inertia, respectively.
Regarding the disk as a mass point and ignoring the influence of the disk on the vibration mode of beam, the kinetic energy of disk can be represented in the following form: where   is the mass of disk and   is the distance between the position of disk and left end point of shaft.The total kinetic energy of rotor system can be obtained as follows: Analyzing the model of shaft and considering the torsion, bending, and shearing, the total potential energy of shaft can be given as follows: where   ,   ,   ,   , and   denote the shaft's shear elastic modulus, elastic modulus, polar rotary inertia, torsional rigidity, and cross-sectional area, respectively.

The Vibration Mode Function of Shaft with Elastic
Support.In order to simulate the boundary condition of shaft more accurately, the shaft is simplified as an elastic support beam with springs at both ends and simplified model is shown in Figure 2, the span and mass of uniform vertical elastic support beam  are  and , respectively, and the stiffness of spring at  and  ends is  1 and  2 , respectively.Assuming the central axis of each section of beam has bending vibration at the same plane, applicating plane assumption in the process of vibration, excluding the effects of rotational inertia and shear deformation, and ignoring the rotation around the center axis of the cross-section and damping, the free vibration differential equation of beam has the following expression: where (, ) donate the transverse displacement of the beam, and it is the function of two variables including the position of the section  and the time .The solution of ( 5) is separated from space and time.Assumed modes method is used to discretize shaft deformation, so the beam transverse displacement (, ) can be written as follows: where () is vibration mode function.Substituting (6) into (5), the differential equation of beam can be gotten: where Equation ( 7) is a four-order linear ordinary differential equation with constant coefficients; it can be solved by using Euler transformation  () =  1 sin  +  2 cos  +  3 sh  +  4 ch , (9) where  is the natural frequency of beam and  1 ,  2 ,  3 , and  4 are integral constants.
As the elastic support beam is adopted, boundary conditions can be expressed as Assuming  =  and substituting (10) into (9), we can get the equation about  2 and  4 : Substituting ( 11) and ( 12) into (9), then If there are nonzero solutions for  1 ,  2 , and  3 , the determinant of the coefficients must be zero in (13).The root is obtained by the method of numerical solution, and then substituting the root and ( 12) into (9), the vibration mode function of elastic support beam can be written as where  and  are Young modulus and moment of inertia of shaft, assuming  1 =  2 .

Discretization of the Energy Equation of the Elastic Rotor
Model.Vibration model functions of the direction of , , and torsion are   (),   (), and Φ  , respectively.Introducing regular coordinates   (),   (), and   .The assumed mode method is adopted to discretize the continuous system: that is, Assuming that the shear force  is constant, the bending angle obtained as follows: Substituting ( 15)∼( 17) into (3), the discrete kinetic energy of shaft is shown in Appendix A.
Ignoring the member does not appear in the equation; it is written as follows: Discretizing the kinetic energy of disk, the discrete kinetic energy of disk is shown in Appendix A.
The total kinetic energy of the rotor system is In the same way, discretizing the potential energy of shaft, the discrete potential energy of shaft is shown in Appendix A.
In order to facilitate the representation and computation, it is written as Substituting the kinetic energy equation ( 20) and potential energy equation (21) into the Lagrange equations yields the discretized equations of motion in matrix notation as where matrices M  , C  , G  , and K  are mass, damping, gyroscopic, and stiffness matrices of the rotor system, and the meaning of specific parameters is shown in Appendix A, where matrices M  , C  , G  , and K  are mass, damping, gyroscopic, and stiffness matrices of the rotor system.

The Establishment of Energy Equation of Rotating Blade.
Assuming there are   blades evenly distributed on the rigid disk, so the disk-blade system is a circular symmetrical structure.When the th blade is deformed, the schematic diagram of disk-blade system is shown in Figure 3.
The displacement of arbitrary point Q on the blade in the fixed coordinate system can be written as where   ,   , and   are the displacement of , , and  direction of the disk in the global coordinate system.Α 0 and A 1 are the transformation matrix from local coordinate system to the rotational coordinate system of the blade.A 2 is the transformation matrix from rotating coordinate system to the swinging coordinate system.A 3 is the transformation matrix of the coordinate system caused by the swing.The specific meaning is given in Appendix B.
When   > 3, the circular symmetric structure has good geometrical and mechanical properties.Without taking into account the axial motion of rotor, the kinetic energy of the th blade is Specific expression is as follows: where   ,   , and   are density, cross-sectional area, and length of the blade, respectively;  is radius of disk;  = Ω + ( − 1)2/  ;   is the twist angle at disk hub.
Considering the effects of bending, circumferential compression, and centrifugal and normal-force-generating potential energy of blade, total potential energy of th blade can be expressed as where   and   are Young's modulus and cross-section area moment of inertia of blade.
Because the movement of end of the rigid disk and the blade is consistent, the modal function of cantilever beam is chosen for blade: Vibration model functions of the displacement of blades  and V are   () and   ().The assumed mode method is adopted to discretize the continuous system: that is, Substituting (29) into ( 26) and ( 27), the discrete energy equation of th blade is given by Substituting kinetic energy equation (32) and potential energy equation (30) into simplified Lagrange equation, the differential equation of vibration of th blade is obtained as follows: where M, C, G, and K are mass, damping, gyroscopic, and stiffness matrices of the blade.

Matrices Assembly of the Rotor-Blade Coupling System.
Discretize the differential equations of motion of the th blade and the rotor system obtained and write in matrix forms.
Because each blade has a corresponding mass matrix , damping matrix , gyroscopic matrix , and stiffness matrix  and blade-disk system is cyclic symmetric structure, the matrices of each blade and rotor system can be assembled to the global matrix of rotor-blade coupling system.The mass, gyroscopic, stiffness, and damping matrices have been assembled as shown in Figure 4.
Assembled differential equations of motion of a rotorblade coupled system are written as where the meaning of the parameters in the formula is shown in Appendix C.

Dynamic Eigenvalue Analysis of a Rotor-Blade Coupling System
In order to further study the vibration characteristics of the coupling system, it is necessary to analyze the critical speed.Specific parameters are shown in "Geometric and Material Properties of Rotor-Blade System." In order to reduce the length of article, only these effects on critical speed of rotor-blade coupling system are studied in this section: the number of blades   from three to seven, the position of disk, and the support stiffness of shaft.

The Influence of Blades Number on Critical Speed
3.1.1.When the Rotational Speed  = 0. Figure 5 shows how the first two-order frequencies of blade change due to the number of blades in a three-to-seven-blade rotor when the position of disk is different.1 1 and 1 1 belong to a set where the blade's first mode predominates in the rotorblade coupling system.1 1 and 1 1 belong to a set where the blade's second mode predominates in the shaft-blade coupling system.Note that there are two reference marks at 163.15 and 1022.5, respectively, denoting the cantilevered blade's first-and second-order frequency.Because of the coupling between shaft and blade, there are two coupling modes dominated by blade's first two-order modes which Mathematical Problems in Engineering    could be found: shaft-blade (SB) and blade-blade (BB).What is more, it also reveals a significant phenomenon that the BB modes are of repeated frequencies of (  −2) multiplicity and the SB modes are of repeated frequencies of (2) multiplicity as two transverse vibration directions of shaft are taken into consideration.And with the increase of the number of blades, the frequency of SB mode increases linearly but that of BB mode is constant.In addition, it also shows frequency changes of both two modes, respectively, when the disk is located in 1/3 and 1/2 of the shaft.It could be found that the frequency of BB mode does not change but the frequency of SB mode changes with the change of disk position, and the changing trend of blades number between first-and secondorder modes is different.
Figure 6 demonstrates how the first two-order frequencies of rotor change due to the number of blades in a three-to-seven-blade rotor when the position of disk is different.It could be found that 1 1 and 1 2 belong to a set where the rotor's first mode predominates in the shaftblade coupling system and 1 1 and 1 2 belong to a set where the rotor's second mode predominates in the shaft-blade coupling system.Note that 93.75 Hz, 95.48 Hz, 322 Hz, and 329.22 Hz, respectively, denoting the rotor's first-and secondorder frequency when the disk is located in 1/3 and 1/2 of the shaft.Furthermore, it also shows that the frequency of rotor decreases with the increase of the number of blades and the trend of decrease is mainly linear.In addition, it is also found that the position of disk can also have a certain influence on the vibration of shaft-blade coupling system.Therefore, it will be investigated in next section.first mode predominates in the rotor-blade coupling system.1 2 , 1 3 , and 1 4 belong to a set where the blade's second mode predominates in the rotor-blade coupling system.It could be found that the Campbell diagram of blade is affected by the number of blades, and with the increase of the number of blades, frequency and corresponding Campbell diagram of BB mode are constant, but frequency of SB mode increases and corresponding Campbell diagram of SB mode shows the trend of increasing.These phenomena can be verified with Figure 5.

When the Rotational
Campbell diagram of rotor changes due to the number of blades are shown in Figure 8. 1 1 and 1 1 belong to a set where the rotor's first mode predominates in the rotorblade coupling system and 1 2 and 1 2 belong to a set where the rotor's second mode predominates in the rotor-blade coupling system.It could be found that the Campbell diagram of rotor is affected by the number of blades, and with the increase of the number of blades, frequency of rotor decreases and corresponding Campbell diagram of it shows the trend of decreasing.What is more, the trend is more obvious in the first-order Campbell diagram.These phenomena can be verified with Figure 6.

The Influence of Disk Position on Critical Speed
3.2.1.When the Rotational Speed  = 0. Figure 9 shows how the first two-order frequencies of blade change due to the position of disk when the number of blades is different.1 1 , 1 1 , and 1 2 belong to a set where the blade's first mode predominates in the rotor-blade coupling system.1 2 , 1 3 , and 1 4 belong to a set where the blade's second mode predominates in the rotor-blade coupling system.Note that Mathematical Problems in Engineering  there are two reference marks at 163.15 Hz and 1022.5 Hz, respectively, denoting the cantilevered blade's first-and second-order frequency.It could be found that the natural frequency of the SB mode is affected by the position of the disk.And the influence is symmetrical about the middle of shaft.When the disk moves from the end to the center of the shaft, the first-order frequency increases gradually and reaches the maximum when the disk is in the center of shaft, but the second-order frequency increases first and then decreases and reaches the minimum when the disk is in the center of shaft.What is more, the trend of the change is more obvious with the increase of the number of blades.
Figure 10 illustrates how first two-order frequencies of rotor change due to the position of disk in a three-to-sevenblade rotor when the number of blades is different.1 1 and 1 2 , respectively, belong to a set where the rotor's first mode predominates in the three-blade and seven-blade rotor-blade coupling system.1 1 and 1 2 , respectively, belong to a set where the rotor's second mode predominates in the threeblade and seven-blade rotor-blade coupling system.It could be found that the natural frequency of rotor is affected by the position of the disk and changes symmetrically with the center of shaft.When the disk moves from the end to the center of the shaft, the first-order frequency decreases gradually and reaches the minimum when the disk is in the center of shaft, but the second-order frequency decreases first and then increases and reaches the maximum when the disk is in the center of shaft.What is more, the trend of change is not affected by the number of blades.to the position of disk.1 1 , 1 1 , and 1 2 belong to a set where the blade's first mode predominates in the rotorblade coupling system.1 2 , 1 3 , and 1 4 belong to a set where the blade's second mode predominates in the rotorblade coupling system.It could be found that the Campbell diagram of blade is affected by the position of disk.When the disk moves from 1/3 to the center of the shaft, firstorder frequency and corresponding Campbell diagram of BB mode are constant, and first-order frequency of SB mode increases and corresponding Campbell diagram of SB mode shows the trend of increasing but second-order frequency and corresponding Campbell diagram behave oppositely.These phenomena can be verified with Figure 9.

When the Rotational
Campbell diagram of rotor changes due to the position of disk is shown in Figure 12. 1 1 and 1 1 belong to a set where the rotor's first mode predominates in the rotorblade coupling system.1 2 and 1 2 belong to a set where the rotor's second mode predominates in the rotor-blade coupling system.It could be found that the Campbell diagram of rotor is affected by the position of disk.When the disk moves from 1/3 to the center of the shaft, firstorder frequency of rotor decreases and corresponding Campbell diagram of rotor shows the trend of decreasing but second-order frequency and corresponding Campbell diagram behaves oppositely.These phenomena can be verified with Figure 10.

The Influence of Support Stiffness on Critical Speed.
In order to make the effect of support stiffness on critical speed more obvious, in this section   = 7, D_loc = 1/2.

3.3.1.
When the Rotational Speed  = 0. Figure 13 shows how the first two-order frequencies of blade change due to the support stiffness of shaft.1 1 and 1 1 belong to a set where the blade's first mode predominates in the rotor-blade coupling system.1 2 , 1 1 , and 1 1 belong to a set where the blade's second mode predominates in the rotor-blade coupling system.It could be found that the natural frequency of BB mode is constant but that of SB mode is constant at first and then increases with the increase of support stiffness of shaft.What is more, the increasing trend of first-order frequency is more obvious.That is, because the first-order frequency of blade is closer to that of shaft, the coupling phenomenon is more obvious.
Figure 14 demonstrates how the first two-order frequencies of rotor change due to the support stiffness of shaft.1 1 , 1 2 , 1 3 , and 1 4 belong to a set where the rotor's first mode predominates in the rotor-blade coupling system.1 1 , 1 2 , 1 3 , and 1 4 belong to a set where the blade's second mode predominates in the rotor-blade coupling system.It could be found that the natural frequency of rotor is affected by the support stiffness quite obviously.When the support stiffness is small, the frequency of rotor is nearly zero.However, when the support stiffness reaches a critical point, the frequency of rotor increases rapidly with the increase of support stiffness.When the support stiffness is large enough, the frequency of rotor is constant.support stiffness, the frequency and corresponding Campbell diagram of BB mode are constant, but the frequency of SB mode increases and corresponding Campbell diagram of SB mode shows the trend of increasing.These phenomena can be verified with Figure 9.

When the Rotational Speed
In order to further analyze the influence of the change of support stiffness, it is analyzed that the critical speed of blade changes due to the support stiffness of shaft when  = 30 rad/s.
First two-order critical speeds of blade change due to the support stiffness of shaft are shown in Figure 16.1 1 and 1 1 belong to a set where the blade's first mode predominates in the rotor-blade coupling system.1 2 , 1 2 , and 1 2 belong to a set where the blade's second mode predominates in the rotor-blade coupling system.It could be found that the critical speed is affected by support stiffness.The critical speed of BB mode is constant but that of SB mode is constant at first and then increases with the increase of support stiffness of shaft.What is more, the trend of increasing of first-order frequency is more obvious.
Figure 17 shows how the Campbell diagrams of rotor change due to the support stiffness of shaft.1 1 , 1 1 , and 1 1 belong to a set where the rotor's first mode predominates in the rotor-blade coupling system.1 2 , 1 2 , and 1 2 belong to a set where the blade's second mode predominates in the rotor-blade coupling system.It could be found that the Campbell diagram of rotor is affected by support stiffness and with the increase of support stiffness, the natural frequency  of rotor increases and corresponding Campbell diagram of rotor shows the trend of increasing.These phenomena can be verified with Figure 10.
In order to further analyze the influence of the change of support stiffness and make the change more obvious, it is analyzed that the critical speed of rotor changes due to the support stiffness of shaft when  = 1000 rad/s.
Figure 18 demonstrates how the first two-order critical speeds of rotor change due to the support stiffness of shaft.1ℎ 1 and 1 1 belong to a set where the rotor's first mode predominates in the rotor-blade coupling system.1ℎ 2 and 1 2 belong to a set where the rotor's second mode predominates in the rotor-blade coupling system.It could be found that the critical speed is affected by support stiffness.The critical speed of rotor increases with the increase of support stiffness of shaft.What is more, the trend of increasing of second-order frequency is more obvious.

Analysis of Mode Shapes
In order to further analyze the vibration characteristics of rotor-blade coupling system and verify the above phenomena, the influence of different factors on mode shapes of system is analyzed in detail.In order to make the analysis more reasonable, the mode shape of first-order natural frequency of the blade is analyzed prior to the formal analysis.
Here, the number of blades   = 4 and the support stiffness of shaft  = 1 × 10 7 N/m. Figure 19 demonstrates the mode shape of 4-blade rotorblade coupling system.It could be found that there is no obvious difference between the mode shape of SB mode and BB mode of blade; that is to say it is difficult to distinguish these two modes through their mode shapes.Furthermore, it is also revealed the shaft does not bend obviously in the mode of blade but the blade bends quite obviously in the mode of rotor.Therefore, a conclusion is drawn that the mode of blade has no effect on the deformation of rotor but the mode of rotor has an obvious effect on the deformation of blade.
Based on the conclusion that has been drawn, when discussing the influence of different factors on the mode shape of coupling system, the mode shape of blade is represented by the first-order mode shape of SB mode.Similarly, the mode shape of rotor is represented by the first-order mode shape of rotor.

The Influence of Blades Number on Mode Shapes.
Mode shapes of rotor-blade coupling system changes due to the number of blades are shown in Figure 20. Figure 20(a) is blade mode shape of rotor-blade coupling system with 4 blades; Figure 20(b) is blade mode shape of rotor-blade coupling system with 4 blades.Figure 20(c) is rotor mode shape of rotor-blade coupling system with 4 blades; and Figure 20(d) is rotor mode shape of rotor-blade coupling system with 7 blades.It could be found that there are no obvious changes of blade and shaft bending deformations with the increase of the number of blades.This shows that the number of blades has no significant effect on the mode shape of blade and rotor.

The Influence of Disk Position on Mode Shapes. Figure 21
illustrates how the mode shape of rotor-blade coupling system changes due to the position of disk. Figure 21(a) is blade mode shape of rotor-blade coupling system when the position of disk is 1/6; Figure 21(b) is blade mode shape of rotor-blade coupling system when the position of disk is 1/2.Figure 21(c) is rotor mode shape of rotor-blade coupling system when the position of disk is 1/6; Figure 21(d) is rotor mode shape of rotor-blade coupling system when the position of disk is 1/2.It could be found that when the disk moves from the end to the center of shaft, the bending deformation of blade remains almost unchanged but that of shaft increases obviously.This shows that the position of disk has an obvious effect on the mode shape of rotor.

The Influence of Rotational Speed on Mode Shapes.
Figure 22 shows the mode shape of rotor-blade coupling system change due to the rotational speed.Figure 22 is blade mode shape of rotor-blade coupling system when the rational speed is 100 rad/s; Figure 22(b) is blade mode shape of rotorblade coupling system when the rational speed is 1000 rad/s.Figure 22(c) is rotor mode shape of rotor-blade coupling system when the rational speed is 100 rad/s; Figure 22(d) is rotor mode shape of rotor-blade coupling system when the rational speed is 1000 rad/s.It could be found that, with the increase of rotational speed, the bending deformation of blade remains almost unchanged but that of shaft increases obviously.This shows that the position of disk has an obvious effect on the mode shape of rotor.

Conclusion
The model of rotor-blade coupling system model is established, the mode function of elastic support beam is deduced, the differential equation of system is deduced by the semianalytical method, and the influences of blades number, disk position, and support stiffness on system characteristic value, Campbell diagram, and mode shape of vibration are studied by numerical simulation.
(1) Because of the coupling between shaft and blade, there are SB and BB modes dominated by blade's first two-order modes which could be found and the BB modes are of repeated frequencies of (  − 2) multiplicity and the SB modes are of repeated frequencies of (2) multiplicity.Frequencies of both modes could be found in the Campbell diagram of blade all the time.
(2) With the increase of the number of blades, the frequency of SB mode increases linearly, that of BB mode is constant, and that of rotor decreases linearly.
The changing trend of Campbell diagram is the same as natural frequency.
(3) Frequencies of SB mode of blade and rotor mode are affected by the position of disk and change symmetrically with the center of shaft.What is more, the changing trend of first-and second-order frequency is different.But that of Campbell diagram is the same as natural frequency.
(4) It could be found that there is no obvious difference between the mode shape of SB mode and BB mode of blade.Furthermore, the mode of blade has no effect on the deformation of rotor but the mode of rotor has an obvious effect on the deformation of blade.The number of blades has no significant effect on the mode shape of blade and rotor but the position of disk and rotational speed have an obvious effect on the mode shape of rotor.

A. Vectors and Matrices Related to the Rotor System
(1) Expression of discrete kinetic energy of shaft can be shown as follows: (2) Expression of discrete kinetic energy of disk can be expressed as (3) Expression of discrete potential energy of shaft can be written as

C. Vectors and Matrices Related to the Rotor-Blade Coupling System
(1)   is generalized coordinate vector of rotor-blade coupling system.
where   and   are vectors of translational and torsional degrees of freedom of rotor;   is the vector of degrees of freedom of blade.
(2)   is mass matrix of rotor-blade coupling system.
where   and   are vectors of translational and torsional mass matrix of rotor;   is the vector of mass matrix of blade;   and   are coupling mass matrix of system.
(3)   is damping matrix of rotor-blade system (including proportional damping and gyro matrix) where   and   are vectors of translational and torsional stiffness matrix of rotor;   is the vector of stiffness matrix of blade.

Figure 1 :
Figure 1: Schematic diagram of the model of rotor-blade coupling system.

Figure 2 :
Figure 2: Schematic diagram of the simplified elastic support beam model.

Figure 4 :
Figure 4: Schematic diagram of assembly of matrices of the rotor-blade coupling system.
First-order frequency changing diagram Second-order frequency changing diagram

Figure 5 :Figure 6 :
Figure 5: Frequency of blade changes due to the number of blades in a three-to-seven-blade rotor.

Figure 7 :Figure 8 :
Figure 7: Campbell diagram of blade changes due to the number of blades.
Second-order frequency changing diagram

Figure 9 :Figure 10 :
Figure 9: Frequency of blade changes due to the position of disk.

1 Figure 11 :Figure 12 :
Figure 11: Campbell diagram of blade changes due to the position of disk.

Figure 13 :Figure 14 :
Figure 13: Frequency of blade changes due to the support stiffness of shaft.

MathematicalFigure 15 :Figure 16 :
Figure 15: Campbell diagram of blade changes due to support stiffness of shaft.

Figure 17 :
Figure 17: Campbell diagram of rotor changes due to support stiffness of shaft.

Figure 18 :
Figure 18: Critical speed of rotor changes due to the support stiffness of shaft.

Figure 22 :
Figure 22: Mode shapes of rotor-blade coupling system due to the rotational speed.