Dynamic Analysis of Flexible Shaft and Elastic Disk Rotor System Based on the Effect of Alford Force

Taken the flexible shaft and elastic disk rotor system as the research object, vibration differential equations are established by the modal synthesis method. Based on the Runge–Kutta method of direct integration, dynamic equations are solved successfully. -en, the effects of Alford force caused by the airflow excitation are researched on the stability and dynamic response. -e results show modal vibration of the disk takes a great effect on the coupled natural frequency of shaft transverse vibration in the high rotational speed region. -e effect of shaft-disk coupling makes rotor frequency to decrease, and different rotational speeds have little effect on rotor frequency. Eccentric distance is sensitive to the amplitude of rotor response, but it has little effect on rotor stability, motion law, and relative relationship of two response frequencies. -e proportion of rotor responses caused by aerodynamic excitation force becomes larger with the increase of rotational speed and disk radius. However, the proportion decreases with the second support moving to the right, and it is helpful for the stable operation of the rotor system.


Introduction
As a kind of large rotating machinery, a gas turbine is widely applied in industrial fields such as aviation, aerospace, transportation, and electric power.It takes an important role in national economy.
A gas turbine works under the high-speed environment, and various nonlinear factors are shown on strong coupling, high dimension, and unsteady state.ere will be hidden serious problems in operation, and it is often huge and even disastrous if the accidents occur.
For a discontinuous system of the mass continuous distribution axis, it is difficult to express the dynamic characteristics of all parts in an equation.Gao and Meng [1] studied the vibration and stability characteristics of elastically supported cantilever beam by the transfer matrix method.But, this research was only on a continuous beam and shaft system, the effect of the disk was neglected on dynamic characteristics.By combining the modal synthesis method with the transfer matrix method, Han et al. [2] solved the dynamic characteristics of a complex rotor system.
With the appearance of optimization design techniques, the rotor structure is continuously improved.e wheel disk is made much thinner so that coupling vibration is more obvious between blade and disk.So elastic disk modeling is very important for the research on the gas turbine disk-shaft system.Zhang [3] established the perturbation motion equation of the flexible shaft-elastic disk rotor.Based on the finite element method, Meek et al. [4,5] analyzed the dynamic characteristics of a rigid disc-flexible shaft rotor system.Hili et al. [6] developed an elastic disk-rigid shaft model with an elastic support, and the effects of elastic deformation and elastic support are discussed on rotor responses when the disk is under axial load.
In addition, fluid-solid coupling problems affect rotor stable operation seriously.It easily leads to the instability phenomenon by airflow excitation in gas turbine working [7,8].Because there is always the relative eccentricity between rotor blade and casing, the clearance of the blade tip is unevenly distributed along the circumferential direction.
is makes the distributions of blade efficiency, pressure, and aerodynamic force different.
So except for a torque, blade circumferential aerodynamic force is also combined with Alford force, which is acting on the axis and perpendicular to eccentric shifting [9].Alford force is a cyclic force.Once it is larger than the external damping to maintain the system stable, positive work will be performed on the rotor system.
at is inputting energy to promote the development of rotor whirl motion.e self-excited vibration will occur when its amplitude increases continuously.
A steam turbine was made by GE corporation, and a strong fluid gap vibration was produced when the load was increased [10]. is is the first time that the phenomenon of aerodynamic excitation was discovered in a rotating machine.Soon after, omas also reported a similar vibration phenomenon [11].He first proposed the concept of gap excitation and suggested the calculation formula for turbine gap excitation force.However, the problem did not attract enough attention at that time.
In 1957, the first ultra-supercritical pressure units were operated at American Ohio Ferrer Power Plant.After that, the former Soviet Union, West Germany, and Japan had also developed and taken into operation the ultra-supercritical pressure units.In working operation, there are quite a few aerodynamic excitation problems on several supercritical units.Because of the higher steam parameters, aerodynamic excitation of a high-pressure rotor is increased on highpower ultra-supercritical pressure units.
In the 1960s, the problem of aerodynamic excitation occurred again.At this period, the gap excitation mainly took place in the turbine of high-power ultra-supercritical pressure units which were manufactured by the Federal Republic of Germany.In the study of aeroengines stability, American scholar J. S. Alford achieved the important progress in theory.In his research, the theory of gap-flow excitation was proposed for the first time, and the calculation formula of turbine excitation force was given.
en, the problem of gap excitation was given enough attention.
For experimental verification of Alford's force, Urlichs and Wohkab [12,13], respectively, developed the experimental research on turbine engines.Vance and Laudadio [14] performed the test of lateral force acting on the blower impeller and confirmed the presence of Alford force.In the sealing force models proposed by omas and Alford, circumferential flow velocity is ignored at the inlet, which is the primary factor of crosscut forces.
Soviet scholars studied the spiral effect of aerodynamic flow in a high-pressure turbine rotor and found that the inertial effect of aerodynamic flow is far greater than the friction effect in the sealing gap [15].At the end of the twentieth century, more researchers began to study aerodynamic excitation forces and the interaction forces between impeller and flow section [16,17].Shanghai Power Equipment Research Institute conducted an experiment on the aerodynamic oscillation of a gas turbine [18].Considering the excitation of aerodynamic force, a small-scale modeling test rig was established by the single-span flexible rotor.en, a simplified calculation method was proposed for the airflow excitation stiffness.e stability and unstable rotation speed were analyzed on the 200 MW turbine power generators.
Chai et al. [19] developed a further research on the airflow exciting force.Taken from the fluid dynamics, the Alford force formula of the uniform flow field was derived on a straight blade and a twisted blade.Chen and Sun [20] proposed an entirely new method of oscillatory fluid dynamics, and it was used for solving the airflow stability and dynamic coefficient of air seal.
Lu et al. [21] analyzed aerodynamic excitation in the cylinder of 200 MW high-pressure turbine and discussed the effect of aerodynamic excitation on the stability of the rotor system.For modifying the formula of Alford force, Yang et al. [22] decomposed rotor displacement into static displacement and dynamic displacement.
is provided a method for adjusting rotor static displacement to reduce the exciting force.
In this paper, vibration differential equations of an elastic disk-flexible shaft rotor are established by the modal synthesis method.Based on the Runge-Kutta method of direct integration, dynamic equations are solved successfully.
en, the effects of Alford force caused by the airflow excitation are researched on the stability and dynamic response.

Vibration Characteristics of Disk-Shaft Coupling System
Figure 1 is the model of the elastic disk with inner diameter a and outer diameter b.For keeping the disk orthogonal to the shaft in all modes, it is assumed that the contact between disk and shaft is rigid.e differential equation of transverse elastic vibration can be derived by Hamilton's principle.e free vibration equation of the disk is expressed as In order to simulate dynamic characteristics of a flexible shaft and an elastic disk rotor system, the elastic circular plate with an inner clamp edge and outer freedom edge is adopted in the research.
Bending displacement of the elastic disk is expressed as where U i (r, θ) is the orthogonal mode of the disk and q i (t) is the time-related part of generalized coordinates.Here, it is 2 Shock and Vibration the simple harmonic vibration, so € q i (t) −ω 2 i q i (t), where ω i is the i-order natural frequency of the elastic disk.U i (r, θ) can be expressed as where U mn is indicated as the orthogonal mode with m nodal diameters and n circles.U i (r, θ) can be expressed as U mn by the reordering of natural frequencies increasing.R mn (r) is the radial distribution of mode shapes; a mn and b mn are constants which determine the initial angle of the disk combined mode.
For a full symmetric disk, theoretical analysis and experiments have shown there is no xed initial angle for the disk mode whether it is uniform or not.To simplify the following formula, it is generally assumed that the initial angle is zero.When boundary conditions are certain, R mn (r) has the inherent form.
Inertia force of the vibrating disk is expressed as e direction of inertia moment at p point is perpendicular to the axis.It is expressed as From equations ( 4) and ( 6), it can be found that the zerodiameter mode of the disk is coupled with longitudinal vibration of the shaft by the e ect of inertia force F p .In addition, one pitch diameter mode of the disk is coupled with bending vibration of the shaft by the e ect of inertia moment.However, other pitch diameter modes of the disk are not coupled with bending vibration of the shaft.So the bending coupled vibration is only discussed, and one pitch diameter mode of the disk is selected.

Dynamic Modeling of Flexible Shaft and
Elastic Disk Rotor System

Rotor Dynamics Model of Flexible Shaft and Elastic
Disk. e model of the exible disk-exible rotor system is shown in Figure 2. e disk is located at the leftmost end of the shaft.e distance from the rst spring support to the disc is l 1 , the distance between the two spring supports is l 2 , and the distance from shaft end to the second spring support is l 3 .
For describing motion state of the shaft, four variables v s1 , w s1 and v s2 , w s2 are taken in two planes which are relatively vertical.Among them, v s1 and w s1 indicate small lateral deformation by the e ect of bending, and v s2 and w s2 indicate small lateral deformation by the e ect of shearing.
erefore, total deformation of the shaft is expressed as In the rotating coordinate system, potential energy of the shaft due to bending and shearing can be expressed as where k is the shape factor which is a constant related to the shaft cross-sectional shape.It can characterize the distribution of uneven shear stress in the cross-sectional area.E s , G s , and A s are axis Young's modulus, shear mode, and crosssectional area, respectively.I s is the axis inertia moment related to diameter.Due to translation and rotation, kinetic energy of the axis can be expressed as where ρ s is the density of the shaft.θ sξ and θ sη denote the small-angle rotation of the shaft cross section due to bending, and they can be expressed as Shock and Vibration Kinetic energy of the disk is expressed as where m D is the mass of the disk, J Dx and J Dy are the disc inertia moments of axis and radial direction, v and w are the shaft displacements at the disc location, and θ ξ and θ η denote the small rotation angle caused by shaft bending.e tiny rotation angle is ignored which is caused by the rigid body displacements because it will be eliminated in the Lagrangian equation.
Strain energy of the disk can be expressed as For elastic support in the rotor system, it can be calculated by the principle of virtual work.
e formula is shown as follows: e elastic disk-shaft rotor is approximately analyzed by the modal synthesis method.Ritz base of each substructure is obtained with the subsystem hypothesis morphological method.Eigenfunctions of transverse vibrations from exible beam and shear beam with free ends are selected as the hypothetical function of continuous shaft transverse de ection induced by bending and shearing.Free vibration eigenfunction of the balanced disk is applied as the hypothesis function of the elastic disk in the system.
Transverse deformation of the elastic disk is expressed as where corner marker "1" indicates the rst section diameter mode, corner marker "n" indicates the number of pitch circles (not including pitch circle subjected to the inner circle clamp).J 1 and Y 1 are the Bessel functions of the rst and second types, respectively, and I 1 and K 1 are the improved Bessel functions of the rst and second types, respectively.A n , B n , C n , and D n are the constants which are determined by boundary conditions, the number of pitch circles, and the condition of orthogonality.β n is the dimensionless parameter related to natural frequency ω n , and the relation is expressed as Coe cient matrix determinant of linear equations is zero by the boundary conditions, so β n can be obtained.
According to the Lagrange equation: where L T s + T D − V s − V D , and the following is obtained: As can be seen from Figure 2, each boundary or connection condition is shown as follows.
e boundary condition of O point is expressed as Connection condition of A point is expressed as 4 Shock and Vibration e boundary condition of B point is expressed as e boundary condition of free end is expressed as By solving the equation set formed by the above equations, natural frequency and mode shape of the system can be obtained.

Di erential Equations of Rotor System under Air ow
Excitation Force.In order to analyze the e ect of Alford force on rotor stability, the equivalent cross sti ness is applied to indicate the action of Alford force.Figure 3 illustrates the process of instability caused by equivalent cross sti ness.O point represents the axis before the deformation, and O ′ point represents the axis after deformation.R is the the radius of casing, and δ is the clearance between rotor and casing.Assuming that the rotor is precessing by the e ect of unbalance, the direction of precession is shown by the arrow in Figure 3. F c is the inertia force, F e is the elastic recovery force, F d is the external damping force, and F a is induced by the Alford force.If F a is greater than the force of outer damping F d on the turbine, it will accelerate the positive precession.So the inertial force increases, the de ection increases, the unstable force increases, and the precession is faster.With this development, vibration amplitude of self-energized motion increases until it is limited by the structure or reaches the limit circle (or has caused damage before).
By the action of rotational dragging, the uid will generate a circular ow in the sealing gap, which plays a signi cant role in rotor vibration.In a certain conditions, it will cause the rotor instable and generate self-excited lateral vibration.In 1987, Muszynska and Bently [23] proposed the Muszynska model based on a large number of experiments.
is model introduces the average circumferential speed ratio τ to characterize the overall motion of uid in the gap between the bearing and sealing.Rotational e ect of seal excitation force is the main factor of rotor instability.e reaction force caused by sealing excitation force rotates around the rotor at a xed angular velocity.In xed coordinates, this reaction force can be described as the following mathematical equation: where K, D, and m f indicate the sti ness, damping, and inertia e ects of uid's disturbing motion on the rotor.It is generally taken that cyclic force is rotor tangential driving force constituted by the cross sti ness term τω D. According to the equation, larger uid damping D is easy to make rotor instable.Reducing the average circumferential velocity ratio τ can improve the stability of the rotor-seal system.A large number of experiments [24] have proved the assertion that the application of antiswirling structure can improve the stability of the sealing system.In addition, the main sti ness term in the model, that is, rotor radial sti ness K − m f τ 2 ω 2 , can become a negative value at the e ect of uid centrifugal inertia force.Perturbation analysis and experimental data show that [25] the e ect of uid centrifugal inertia force becomes very obvious, when the ratio of seal gap and rotor radius is small.is has a great in uence on the sti ness matrix.So the radial sti ness term can easily reach a negative value, and it is also a factor that causes rotor instable.e Muszynska model correctly re ects the basic dynamic characteristics of sealing force.Equation ( 23) not only is applied in the condition of small displacement at the rotor center position but also can indicate the nonlinear sealing force at larger displacement.Experimental and numerical results con rm that [26] τ, K, and D are the nonlinear Shock and Vibration functions of disturbed displacements x and y.It can be expressed as follows: where e � ( ������ x 2 + y 2  /C r ) is rotor relative partial displacement and C r is the sealing clearance.Specific seal is approximately described by selecting the appropriate exponents n and b.Due to the effects of axial flow, surface roughness, secondary flow, and other factors, in general τ 0 < (1/2).K 0 , D 0 , and m f can be calculated by the combination of Childs seal dynamic coefficient equation and equation (24): e Muszynska sealing force model is mainly based on a large number of experiments and abstracted from the model.It has a strong generality and grasps the main mechanical characteristics of sealing force.
e correctness of the Muszynska model has been widely recognized in the application.e Muszynska sealing force model is adopted as the model of airflow excitation force, and the dynamic differential equation of the elastic disc-flexible shaft rotor is established by the modal synthesis method.
Applying the principle of virtual work, virtual work caused by the Alford force can be obtained as where According to the modal synthesis method, the dynamic equation of the rotor can be expressed as where F e is the imbalance force caused by eccentric mass and F q is the force of bearing oil film.F q is expressed as where So K q , C q , G q , and M q in equation ( 29) are expressed as Assuming that q � q/C r , _ q � (dq/dτ)(dτ/dt) � Ωq ′ , and e � e/C r , equation ( 28) is expressed as (32)

Analysis and Discussion of the Results
In order to explain the effects of different parameters on rotor vibration responses, rotor parameters are taken as follows if no special description is given.e geometric parameters of rotor structure are indicated as l 1 � 0.   4(a) and 4(b), it can be seen that the vibration of the disk is not obvious at the low-order frequency.As the frequency increases, the vibration of the disk gradually becomes prominent.
e above derivation is performed in the rotating coordinate of shaft and the rotating coordinate of the disc, so the obtained natural frequencies are also relative to the rotating coordinate system.In order to compare with the experimental data, it is necessary to convert the obtained natural frequency into an inertial coordinate system.For the transverse vibration of the shaft and disk, the conversion relationship is different between rotating coordinate system and inertial coordinate system.Hence, there are two different natural frequencies in the inertial coordinate system.

Analysis of Coupled Natural Frequency in Transverse
Vibration.For transverse vibration of the shaft, conversion relation is expressed as where ω is the natural frequency in the rotating coordinate system, Ω is the angular velocity of the shaft, and ω 0 is the natural frequency in the inertial coordinate system.For transverse vibration of the disk, the transformation is difficult from the rotating coordinate system to the inertial coordinate system.Compared with the rotation angle Ωt, declination angles θ ξ and θ η of transverse vibration can be ignored in the transformation of the coordinate system.So, the relationship between two angular coordinate systems can be expressed as Excitation force of the disc is F 0 cos ω 0 t, and vibration response at r 0 � ρ and θ 0 � 0 is expressed as where ω is the natural frequency in the rotating coordinate system and m is the number of mode nodal diameters in the disk.It can be seen from the above equation that there are two resonant frequencies in the inertial coordinate system: Plotting the chart of frequency-rotational angular velocity, ω − mΩ and −(ω − mΩ) are the continuous curves.
erefore, natural frequencies are generally divided into two parts when rotating.at is, at a certain rotational angular velocity, two different excitation frequencies will excite the same response mode in the inertial space.
For the specific problem, natural frequency of diskshaft coupling is indicated as ω. e antiprecession mode of the rotating system may be negative, and only the first nodal diameter mode of the disk will couple with the vibration of shaft, so m � 1. Equation (36) can be expressed as According to equation (28), the coupled natural frequency of shaft transverse vibration is shown as Figure 5.
Here, 5(a) is the diagram of coupled natural frequency in the flexible shaft and rigid disk, and 5(b) is the diagram of coupled natural frequency in the flexible shaft and elastic disk.
In Figure 5, the difference between Figures 5(a) and 5(b) is obvious when the rotational speed is large, and the difference between two models is smaller in the region with lower rotational speed.It is illustrated that modal vibration of the shaft occupies a dominant position in the low rotational speed region.Modal vibration of the disk has little effect on the coupled natural frequency of shaft transverse vibration, while it takes a great effect in the high rotational speed region.Hence, the difference between two models is quite obvious.
According to equation (37), the coupled natural frequency of elastic disk transverse vibration is shown as Figure 6.Here, 6(a) is the diagram of forward natural frequency, and 6(b) is the diagram of backward natural frequency.
Mass matrix and stiffness matrix of the rotor system will be affected by the coupling effect of shaft and disk, but the gyroscopic matrix and damping matrix are not affected.Without considering the physical background, the disk-shaft coupling term can be directly removed in the mass matrix and stiffness matrix.So, it is to analyze the effect of the diskshaft coupling on rotor natural frequency.
From Figure 7, it can be seen that the effect of shaft-disk coupling makes rotor frequency decrease, and different rotational speeds have little effect on rotor frequency.e effects of different structural parameters are discussed on rotor frequency.Figure 8(a) shows the frequency of the rotor system decreases with the increase of outer diameter.From Figure 8(b), the first-order frequency of the rotor system decreases when the first support moves to the right, but the third-order frequency increases.8

Shock and Vibration
In Figure 8(c), the second-order frequency of the rotor system decreases obviously when the second support moves to the right, and it has little e ect on the frequencies of other orders.

Stability Analysis of Elastic Disk-Flexible Shaft at Aerodynamic
Force. e e ects of di erent parameters are discussed on rotor time-domain chart, axis orbit, and spectrum graph.Based on the variation of rotor responses, the stability of rotor system is analyzed.Here, the dimensionless displacement in the y direction is denoted as Y, the dimensionless displacement in the z direction is denoted as Z, the ratio of average circumferential speed is denoted as τ, the frequency of rotor response is denoted as f, and the amplitude of rotor response is denoted as A.

e E ect of Eccentric Distance and Rotational Speed.
Figures 9 and 10 show the responses of the rotor system at di erent imbalance and rotational speed.According to the spectrum graph, the responses of the rotor system include the frequencies caused by aerodynamic excitation force and rotor imbalance.Summarizing the above gures, rotor frequency caused by aerodynamic excitation force can be obtained, and it is about 50 Hz.
Comparing with rotor responses at di erent eccentric distances, there is no clear change in rotor motion law and relative relationship between two response frequencies.It illustrates eccentric distance has little e ect on rotor stability, and only the amplitude of rotor response is changed.
With the increase of rotational speed, the value of vibration frequency caused by the imbalance increases, and the proportion of vibration response caused by aerodynamic excitation force becomes larger.Ultimately, the rotor system is caused unstable when rotational speed increases to be a certain value.

e E ect of Disk Radius.
Figure 11 shows the responses of the rotor system at di erent disk radii.According to spectrum graphs, the proportion of rotor responses caused by aerodynamic excitation force becomes larger with the increase in the disk radius and gradually covers the response caused by the imbalance.Comparing with the

Shock and Vibration 9
charts of the time domain, vibration amplitude of the rotor system increases rapidly with the increase in the disk radius.e rotor system becomes unstable when the disk radius increases to be a certain value.It is because that support center of the whole rotor system is also moving to the right; the mass of disk and Alford force can be balanced.e rotor system would be stable when the second support moves to the right at a certain position.

Conclusions
Vibration di erential equations of the elastic disk-exible shaft rotor are established by the modal synthesis method.
Based on the Runge-Kutta method of direct integration, vibration characteristics of the elastic disk and the exible shaft rotor were analyzed successfully.en, the e ects of Alford force caused by the air ow excitation are further researched on the stability and dynamic response.e results show the following:

Shock and Vibration 11
increase of rotational speed and disk radius.e rotor system becomes unstable when rotational speed or disk radius increases to a certain value.However, the proportion decreases with the second support moving to the right, and it is helpful for the stable operation of the rotor system.12 Shock and Vibration (r)cos(mθ) r dr dθ.

Figure 1 :
Figure 1: Schematic diagram of the disk.

Figure 2 :
Figure 2: Rotor model of the exible shaft and the elastic disk.

Figure 3 :
Figure 3: e model of rotor damage-instability caused by aerodynamic excitation force.
35, l 2 � 0.95 m, and l 3 � 0.65 m. e inner and outer radii, thickness, density, elastic modulus, and Poisson's ratio of the elastic disk are, respectively, indicated as a � 0.11 m, b � 0.445 m, h D � 0.088 m, ρ D � 4640 kg/m 3 , E D � 1.13 × 10 11 N/m 2 , and v D � 0.28.e radius, density, Poisson's ratio, and elastic modulus of the flexible shaft are indicated as r s � 0.11 m, ρ s � 7900 kg/m 3 , v s � 0.28, and E s � 2.06 × 10 11 N/m 2 .e stiffness coefficient of elastic support is indicated as k 11 � k 21 � 6.0 × 10 8 N/m. e damping coefficient of the support is indicated as c 11 � c 21 � 3000 N • s/m.In the model of aerodynamic excitation force, the circumferential loss coefficient of sealed inlet airflow is indicated as ξ � 0.1.e seal length is indicated as l � 0.088.e radial seal clearance is indicated as δ � 2 mm.Fluid average axial velocity in the sealing cavity is indicated as v � 5 m/s.e sealing radius is indicated as R � 0.445 m. e coefficient of fluid viscosity is indicated as μ � 1.8 × 10 −5 pa • s. e empirical coefficients are indicated 6 Shock and Vibration as m 0 � −0.25, n 0 � 0.079, τ 0 � 0.

Figure 4 :Figure 5 :Figure 6 :
Figure 4: Mode shape comparison of the di erent rotor models: rst eight mode shapes of (a) exible shaft and rigid disk rotor system and (b) exible shaft and elastic disk rotor system.

Figure 9 :Figure 10 :
Figure 9: Rotor responses of the time domain, axis orbit, and spectrum at eccentric distances e 0.001 (a) and e 0.01 (b) and rotational speed Ω 800 rad/s.

Table 1 :
First eight natural frequencies of the rigid disk and elastic disk rotor models.