 Open Access
 Authors : Y. Hari Krishna , T. Pavan Kumar , M. Sridhar , Y. Bhavana, C. Sai Surya Bhagavan
 Paper ID : IJERTV8IS110405
 Volume & Issue : Volume 08, Issue 11 (November 2019)
 Published (First Online): 05122019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Dynamic Analysis and Intelligent ElectroDynamic Control of AeroEngine RotorBearing System
Y. Hari Krishna, T. Pavan Kumar, M. Sridhar, Y. Bhavana, C. Sai Surya Bhagavan
Andhra University Andhra Pradesh.
Abstract. In the present work, dynamic analysis of a typical twospool aeroengine turbocompressor rotor system mounted on ball bearings supported with squeezefilm dampers is considered. The dynamic equations are established by the finite element model and the reduced assembled system of equations are solved to obtain transient unbalance response and frequency domain plots during startup (acceleration) conditions. The nonlinear bearing forces with squeezefilm damper loads at the rotor nodes are considered during analysis. An electromagnetic actuator is employed based on a trained neural network data for amplitude control at the specific location on the rotor. Unlike, the conventional PDcontrol scheme, proposed controller can predict the required component forces based on the instantaneous radial displacements at that location. Interactive computer programs are developed to simulate the practical system by considering various forces including gyroscopic, gravity in addition to nonlinear bearing forces.
Keywords: Transient analysis, Beam elements, Electrodynamic actuator, Neural network observer, Twospool engine.
1 INTRODUCTION
Nowadays, rotor models have become more and more complex due to increasing demands of engineering applications. The simplified nonlinear dynamic models of rotors cannot reflect the practical conditions. Therefore, more complex, flexible rotorbearing models such as those using onedimensional Timoshenko beam elements, transfer matrix methods and mode synthesis techniques are needed. Most recent aspect is to control the critical response amplitudes in an online fashion. In this line, various strategies [13] have been proposed in literature to increase the damping such as semiactive magnetorheological fluids, electromagnetic bearings and viscoelastic supports. In practice, augmentation of damping is a difficult task and alternatively active control techniques became the choice. The active control capability to stabilize rotors and reduce unbalance vibration was illustrated in literature by several authors [48]. These works do not consider the timevarying plant dynamics, so that the robustness of controller may not be guaranteed. Robust stability issues were lateron shown by few other papers [911]. It was proved that the model identification and control as integral aspects in the robust controller design.
Recently, the active vibration control of rotors using electromagnetic exciters is one of the important topics in the turbomachinery analysis. Control of resonance amplitudes and reduction of local vibration levels is possible online
with electromagnetic coils. In [1213] and [1416], illustrated the local control approach using linear PD control technique for highspeed rotors mounted on hydrodynamic and rolling element bearings. A hybrid approach of intelligent PID feedback control to suppress vibrations in the critical zones presented [17]. An attenuating method [18] for lateral vibration of rotors with the aid of semiactive dampers and active control of hydrodynamic bearings illustrated. In [19] used primary optimal control gain of genetic algorithm in proportional derivative control to get stability in the rotor system. A nonlinear golden section adaptive control algorithm to reduce the vibration levels of a flexible Cartesian smart material manipulator run by a servomotor implemented [20]. In [21] employed two piezoactuators to actuate two movable bushing to control the rotor at high speeds. The output error signal coming from the proximity probe is responsible to control the piezoactuator acting force at the bushings. The change of the magnetic flux induced in the electric coils is responsible for the controlling damping effect of the dampers. An adaptive control approach [22] proposed to suppress vibrations due to imbalance of masses in the rotors and [23] proposed a method to convert the complex domain modes of a rigid rotor supported on active magnetic bearings to achieve maximum PID control performance in the system. Most of the control schemes proposed in rotors are based on linear concepts. To augment the controller performance in nonlinear rotor systems, an attempt is required to design and test a real time strategy by using some inverse approach.
In the present work, a typical twospool aero engine turbocompressor rotor system supported by ball bearings is considered. The nonlinear dynamic equations are solved using timeintegration scheme and the transient response stability issues are discussed during acceleration/deceleration condition. An electromagnetic actuator design based on a trained neural network is implemented for controlling the vibration amplitudes at a specific location of the rotor. Dynamic analysis and control simulations are illustrated at different speeds of operation. Remainder of the paper is organized as follows: section2 presents the modeling of rotor bearing system with squeeze film damper force considerations. Section3 gives the mathematical concept of electromagnetic actuator and its control schemes. Section4 outlines the results using a numerical example.
2 MATHEMATICAL MODELING
Modern aeroengine rotors have multispool configurations [2430] due to their several advantages. In multispool engine rotors, as the air flows through the compressor, its pressure and temperature rises, increasing the speed. Two or
more concentric shafts operating at different speeds are commonly used. Each spool is given by a separate turbine stage and runs at different speeds. A twospool engine rotor considered in the present work with flexible bearings is shown in Fig. 1.
Fig. 1. Schematic diagram of a twospool engine
The assembly of whole structure is consisting of high pressure (LP) turbine and low pressure (LP) compressor. Four bearings with a casing are used to support the entire structure. The HP shaft on the front end is having the casing to support the SFD bearing through a squirrel case spring and this will provide the actual axial location of parallelism of the HP rotor in the oil annulus. The dynamic analysis of
this rotor system is carried out using a finite element model. Twonode 8degrees of freedom Timoshenko beam elements are used to discretize the assembly. It is assumed that the torsional and axial deformations are negligible. Fig.

shows the simplified model of the twinspool rotor assembly.
Fig. 2. The FE model of LP and HP rotor system.
It consists of four bearings and four disks, two each on inner and outer spools respectively. It is considered that the unbalance acts on all the four disks and bearings are rolling element bearings supported in squeezefilm dampers. The simplified configuration is different from the conventional
model in consideration of nonlinear bearing forces as compared to conventional springdamper bearing elements. As seen in the figure, the node 4 of inner spool and node 10 of outer spool is connected by the fourth bearing element.

Equations of motion
By considering rigid disc, flexible shaft model, the rotor bearing system is discretized into eight beam elements along with four lumped disk masses. Bearing and unbalance forces are accounted as external force vector. By using energy expressions and applying the Hamilton principle, the governing equation of motion for the system are given by []{ } + ([] + (0 + t)[]){ } + []{} = {} +
{} + {} +{} (1)
Where [C], [G] are the viscous and gyroscopic damping matrices. [K] and [M] are the assembled stiffness and assembled mass matrices, {Fe} is the unbalance force vector,
{Fb} bearing force vector, {Fs} is the squeeze film force vector, {Fg} is the gravity force vector and {q} is the vector of nodal displacements. Also, 0 is initial velocity, t is the time and is constant acceleration. Altogether five elements (64=24 degrees of freedom) are taken on inner spool and three elements (44=16 degrees of freedom) are on the outer spool.

SqueezeFilm Damper Bearing
The main advantage of Squeeze film dampers is to attenuate the unbalance synchronous response of a rotor and to unbalance and to reduce subsynchronous rotor dynamic instabilities. The required external damping for the roller element bearings used in Aircraft engines are provided with one or more SFDs to support the rotor. Normally, a lubricating film is provided in between a whirling journal and a stationary housing. The typical rolling element bearing outer race is constrained with a squirrel cage or a dowel pin (elastic) support to arrest the rotation along with the journal. Pressurized lubricant fills the squeeze film lands through the central groove or feed holes. The squeeze oil film will be pressurized due to the rotation of the inner race with the rotor and the outer race will be in stationary position within the housing. Due to the lubricant fluid film reaction forces, dynamic pressure will be generated and that will support to damp the increase rotor whirl motion amplitudes.
Lubricant viscosity, journal kinematics and damper geometry are the important factor in controlling the force generated due to fluid film reaction in a SFD bearing. Fig.3 shows the schematic of SFD bearings.
Fig. 3. Schematic of SFD bearing system.
Here, the outer race (journal) is connected to the casing through a centralizing spring. The bearing forces are due to contact between the balls and races, which are obtained using the following Hertz contact force relations [31]:
F = Nb K (x cos + y sin r
3/2
H. cos
(2a)
bx j=1 b j
j 0) j
F = Nb K
3/2
(x cos + y sin r ) H. sin
(2b)
by j=1 b j j 0 j
where r0, Kb and H are the radial clearance, Heaviside function and the Hertz contact stiffness and is depends on the contact material and shape, respectively. x & y are the coordinates. The jth ball angular location term j can be written as
= Ã— t + 2 (j 1); j=1,2,. N
(3)
j Cage Nb b
where Cage is the ball bearing cage angular velocity is written as Cage=Ã—r/(R+r) with R and r as the radii of the outer and inner races and is the rotational speed of shaft, Nb is balls number. The ball passing frequency or Variable Compliance (VC) is defined as VC=CageÃ— Number of balls (Nb). It is assumed that outer race rotates about its axis, but does not spin in the oil film. The relative displacement between the shaft (s) and outer race (o) is given by, x= xsxo, y=ysyo. The outer race (journal) is subjected to squeezefilm oil forces. The solution of the Reynolds equation for incompressible fluids written below will give the instantaneous pressure distribution p inside the oil film [32]:
2
2
1 {p (p)} + p p = 12(e sin + e cos) (4)
1
1
2 z2
Here, h= c+ecos is called squeezefilm thickness.
The simplified force components of oil film in the xy stationary coordinate system can be obtained from the following important assumptions:
(i) Oilfilm thickness when compared to radius of journal is small, (ii) pressure gradient along the radial direction is small and
(iii) laminar lubricant fluid flow.
The resulting squeezefilm forces are given from the pressure distribution as:
3
3
F = 1L [( I + I )cos ( I + I )sin] (5a)
sx c2 1 2 2 3
3
3
F = 1L [( I + I )sin ( I + I )cos] (5b)
sy c2 1 2 2 3
X2+Y2
X2+Y2
where R1, , c and L are the radius of SFD, oil dynamic viscosity, the radial clearance of squeeze film damper and the land length of squeeze film damper respectively. The eccentricity ratio is given by: = . =tan1 (Y/X). is the angular velocity of
c
SFD and (X,Y) is the displacement of SFD in the fixed coordinate system and I1, I2 and I3 are the bearing integrals and evaluated analytically in closed form and given by [33]:
1
1
+
I1 = 1
1
1
+
I2 = 1
1
1
+
I3 = 1
cos2 (1+cos)3 d
sincos (1+cos)3 d
sin2 (1+cos)3 d
(6a)
(6b)
(6c)
1
1
= tan1 ( ) (6d)
All these forces are incorporated in Eq.(1) at appropriate nodes and the resultant equations are solved. In present case, these squeezefilm damper forces are considered only on LP rotor front and rear bearings.


ELECTROMAGNETIC ACTUATOR MODEL
The electromagnetic actuator with four pole is shown in the Fig. 4. Poles 1 and 2 provide actuation force in vertical (x) direction and poles 3 and 4 provide actuation force in the horizontal (y) direction.
Fig. 4. Fourpole electromagnetic actuator.
The expressions for actuation force in vertical and horizontal direction are:
2
2
F = F F
= A N2cos [ (Ib+i1)
(Ib+i3)2
] (7a)
ax ax1
ax3 0 a
(Gxcos)2
(G+xcos)2
F = F F
= A N2cos [ (Ib+i2)
(Ib+i4)2
2
2
] (7b)
ay ay2
ay4 0 a
(Gycos)2
(G+ycos)2
Where 0 is permeability of freespace (0 = 4 Ã— 107H/m), Aa is the crosssection of the airgap which is assumed to be equal to pole face area, N is the number of turns per pole, G is a nominal air gap between stator and rotor, is half the included angle of a pole ( =22.5Â°), Ib is bias current which is kept fixed and in (n = 1, 2, 3 and 4) are small variable current of corresponding poles which are varied according to control strategy. For the sake of simplicity currents in the diametrically opposite poles are kept equal in magnitude, but opposite in direction, i.e., i1= – i3 and i2= – i4. In the actuator design, the allowable displacement of the rotor is considered about onetenth of the air gap size. At normal point of operation, i1=i2=i3=i4=0 and x=y=0. By expanding with Taylor series and neglecting higher order terms, we get:
=
402
2 1 +
40222
=k i +k x (8a)
=k i +k x (8a)
3 1a 1 2a
=
402
2 2 +
40222
= k i +k y (8b)
= k i +k y (8b)
3 1a 2 2a
Here, k1a and k2a are forcecurrent and forcedisplacement constants of electromagnetic actuator respectively. Computation of currents i1 and i2 as a function of instantaneous errors is the function of controller.

Design of PD Controller
To minimize the vibration amplitude of system appropriate magnetic force should applied by the actuator; force should act in the opposite direction of motion and should increase with displacement. As seen from above simplified equations, the magnetic force depends upon the pole winding current of an electromagnet; so by controlling pole winding currents, the system response can be controlled. ProportionalDerivative (PD) control is generally used for faster response. Proportional control increase the control input in proportion to the error e(t) within the acceptable range of error. Derivative control changes the control input in proportion to rate of change of error . Equation of control output of PD controller using proportional and derivative gains Kp, Kd can be given as
u(t) = Kpe(t) + Kde (t) (9)
The controlled outut is currents and the error signal. Hence, the control currents for pole pair in horizontal direction (i1) and in the vertical direction (i2) are given as [12]:
i1 = Kpex(t)+Kde x(t) (10 a) i2 =
Kpey(t)+Kde y(t) (10 b)
As the set point is having zero velocity and displacement, the terms ex(t) = x and ey(t) = y along with e x(t) = x and
e y(t) = y . Designer has to select the appropriate gain constants Kp and Kd and obtain the control forces Fax and Fay.

Neural network based controller
As seen from Eq.(7), it is observed that the control forces are functions of the velocities of rotor and its radial displacements. The control function in non linear in nature. So, a control concept similar to an inverse dynamic model of the system is used. A nonlinear controller based on trained neural network model with instantaneous transient response as inputs and corresponding currents as outputs is proposed in present work. Fig.5 shows the block diagram of the control methodology.
Fig. 5. Error based control scheme.
Conventional 3layer back propagation neural network is used to train the data. The inputs for the neural network are errors (vector of negative X(t) values) and error rates (vector of negative () values) over a sampling period as the values
of i1 and i2 are varied randomly from 0 to 100 milliamperes according to equations (7). A 4 input (x, y, – and – ) and two output (i1, i2) training patterns generated inversely have been used in training. Learning rate of 0.4 and momentum
factor of 0.01 are implemented. Optimum hidden layer nodes were found to be 4 and the neural network is trained using an interactive computer program to find the set of weight matrices. Further, these are employed in estimating the unknown currents corresponding to a computed nodal displacement and velocity vectors during every time step. The detailed backpropagation algorithm for 3layer multilayer perceptron network can be found in open literature [34].

RESULTS AND DISCUSSIONS
The coupled set of differential equations is evaluated in timedomain series with aid of RungeKutta equation of order four with initial value as zero. The dimensional
properties of the rotor used in the analysis shown in Table 1. A speed ratio of 1.5 is implemented between HP and LP rotors. The material for the spools is steel with modulus of elasticity value E=210 GPa and density =7.8 g/cm3. All four disks have polar moment of inertia (IP) twice the diametric mass moment of inertia (ID). Further, an eccentricity of 10 microns is considered in the disks. The resultant stiffness, mass, damping (damping factor of 0.01) and gyroscopic matrices are used. The SFD parameters used are [35]: =5103 PaS, c=2104 m, L=8.3103 m, R1=0.03

The value of is set as 2000 rpm and angular acceleration is varied.
Table I Geometric and material data for rotor (Referring to Fig.2)
Node
Axial dist.(mm)
douter (mm)
dinner (mm)
Disc mass(kg)
IP (kgmm2)
1
0
30.4
0
–
–
2
76.2
30.4
0
4.904
0.02712
3
323.85
30.4
0
–
–
4
406.4
30.4
0
–
–
5
457.2
30.4
0
4.203
0.02034
6
508
30.4
0
–
–
7
152.4
50.8
38.1
–
–
8
203.2
50.8
38.1
3.327
0.01469
9
355.6
50.8
38.1
2.227
0.0972
10
406.4
50.8
38.1
–
–
Fig.6 depicts the frequency response plot of a system considering ball bearing forces and not considering the forces. Fig.6(a) coincides well with the available results. It is seen that there are multiple ball passing frequencies (VC) along with harmonics in the system, when ball bearing forces are accounted.
Fig. 6. Frequency plots of rotor at speed of 2000 rpm.
Fig.7 shows the transient analysis results at different acceleration conditions in terms of frequency responses. It is obvious that the multiplicities of ball passing frequencies increase with acceleration of rotor. This may be due to a randomly changing bearing excitation.
Fig. 7. Effect of acceleration of rotor at LP compressor disk.
In order to minimize the amplitude levels, the following electromagnetic actuator parameters are considered: N=106, G=9 cm
,Ag=2.4 cm2, =22.5o, and Ib=3 A. The instantaneous current components i1 and i2 depend on the displacements x and y at the disk node1. The software code is modified to incorporate an extra component of magnetic force on the force vector at right side located at disk1 Fig.8 shows the frequency response plots of rotor computed when controlled with bearing force and when uncontrolled with force. It indicates that the amplitudes from second mode have reduced drastically and the ball passing frequency amplitudes have been dropped. The methodology can be extended for other nodes of interest.
Fig. 8. Amplitude reduction at disk1 node with EM controller.


CONCLUSIONS
In current study, a process to reduce vibration with electromagnetic actuators has been briefly presented for a twinspool rotor case study. The considered system has unbalance and gravity loading along under accelerating conditions. Gyroscopic effects arising due to combined effect of the shaft and disks were accounted, and for all disks unbalance was taken into consideration. Frequency response diagrams were used to distinguish the dynamic response at different angular accelerations of the rotor. It was observed that the squeezefilm damper at the front and rear ends of LP rotor have drastically reduced the amplitudes at critical speeds of operation. Further, the neuralnetwork based electromagnetic actuator designed to provide the nonlinear control forces has considerably reduced the response. Unlike the conventional PD control scheme, the present approach, do not search for optimum set of controller gains. However, it requires huge set of inputoutput training data before its implementation in the plant. Although the results were illustrated for a selected speed of LP rotor, the focus may be given at the critical speeds of operation in practice so as to identify the effective operational speed range for amplitude attenuation using EM actuator.
REFERENCES

Zhu C, A disktype magnetorheological fluid damper for rotor system vibration control. J Sou Vib 2005;283:10511069.

Fan YH, Lee AC. Design of a permanent/electromagnetic magnetic bearing controlled rotor system. J Fran Inst 1997;334:337356.

Verichev NN, Verichev SN, Erofeyev VI. Damping lateral vibrations in rotary machinery using motor speed modulation. J Sou Vib 2010;329:1320.

Zhong J, Li L. Fractionalorder system identification and proportionalderivative control of a solidcore magnetic bearing. ISA Trans 2013;

Burrows C, Sahinkaya M, Clements S. Active vibration control of flexible rotors: an experimental and theoretical study. Proc Roy Soc Lon A 1989;422:123146.

Lin YH, Yu HC. Active modal control of flexible rotor. Mechanical Systems and Sinal Processing 2004;18:11171131.

Jang MJ, Chen CL, Tsao YM. Sliding mode control for active magnetic bearing system with flexible rotor. J Fran Inst 2005; 342:401419.

Dimitri A, ElShafei A. Instability control and unbalance compensation of flexible rotors supported on journal bearings using magnetic bearings. Proceedings of the Eighth IFToMM International Conference on Rotor Dynamics, Seoul, Korea; Sep, 2006. p. 657664.

Fittro R, Knospe C. Rotor compliance minimization via control of active magnetic bearings. IEEE Trans Con Sys Tech 2002;10(2):238249.

Riemann B, Perini EA, Cavalca KL, De Castro HF, Rinderknecht
S. Oil whip instability control using synthesis technique on a magnetic actuator. J Sou Vib 2013;332:654673.

Xiaochun G, Dengqing C. Fuzzy proportionalintegralderivative control of an overhang rotor with double discs based on the active tilting pad journal bearing. J Vib Con 2015;19(10):14871498.

Fan CC, Pan MC. Active elimination of oil and dry whip in a rotating machine with an electromagnetic actuator. Int J Mech Sci 2011:53:126134.

Fan CC, Pan MC. Experimental study on the whip elimination of rotorbearing systems with electromagnetic exciters. Mechanism Machine theo 2011;46:290304.

Das AS, Nighil MC, Dutt JK, Irretier H. Vibration control and stability analysis of rotorshaft systems with electromagnetic exciters. Mechanism Machine theo 2008;43:12951316.

Das AS, Dutt JK, Ray K. Active vibration control of unbalanced flexible rotor shaft systems parametrically excited due to base motion. Appli Math Mod 2010; 34:23532369.

Das AS, Dutt JK, Ray K. Active control of coupled flexible torsional vibration in a flexible rotorbearing system using electromagnetic actuator. Int J Nonli Mechan 2011;46:10931109.

Toha SF, Tokhi MO. A Hybrid Control Scheme for a Twin Rotor System with Multi Objective Genetic Algorithm. 12th International Conference on Computer Modelling and Simulation, Emmanuel College Cambridge, United Kingdom; Mar, 2010. p. 2028.

Ferfecki P, Zapomel J. Investigation of vibration mitigation of flexibly support rigid rotors equipped with controlled elements. Proc Eng 2012; 48:135142.

Lin J, Zheng YB. Vibration suppression control of smart piezoelectric rotating truss structure by parallel neurofuzzy control with genetic algorithm tuning. J Sou Vib 2012;331:36773694.

Qiu ZC. Adaptive nonlinear vibration control of a Cartesian flexible manipulator driven by a ballscrew mechanism. Mech Sys Sig Proc 2012;30:248266.

Tuma J, Simek J, Skuta J, Los J. Active vibration control of journal bearings with the use of piezoactuators. Mech Sys Sig Proc 2013;36:618629.

Fang J, Xu X, Tang X, Liu H. Adaptive complete suppression of imbalance vibration in AMB systems using gain phase modifier. J Sou Vib 2013;332:62036215.

Han WJ, Chong WL. Proportionalintegralderivative control of rigid rotoractive magnetic bearing system via eigenvalue assignment for decoupled translational and conical modes. J Vib Con 2015; 21(12):23722393.

Gunter EJ, Li DF, Barrett LE. Unbalance response of a two spool gas turbine engine with squeeze film bearings. ASME Gas turbine conference and products show, Houston, Texas; Mar, 1981.

Shanmugam A, Padmanabhan C. A fixedfree interface component mode synthesis method for rotordynamic analysis. J Sou Vib 2006; 297:664679.

Sun G, Palazzolo A, Provenza A, Lawrence C, Carney K. Long duration blade loss simulations, including thermal growths for dualrotor gas turbine engine. J Sou Vib 2008; 316:147163.

Hai PM, Bonello P. A computational parametric analysis of the vibration of a threespool aeroengine under multifrequency unbalance excitation. J Eng Gas Turb Pow 2011;133:0725041.

Hai PM, Bonello P. A computational parametric analysis of the vibration of a threespool aeroengine under multifrequency unbalance excitation. Journal of Engineering for Gas Turbines and Power 2011; 133:07250419.

Lin FJ, Chen SY, Huang MS. Adaptive complementary sliding mode control for thrust active magnetic bearing system. Con Eng Prac 2011; 19:711722.

Groves KH, Bonello P, Hai PM. Efficient dynamic analysis of a whole aero engine using identified nonlinear bearing models. Proc Inst Mech Eng, Part C: J Mech Eng Sci 2012;226:6681.

Chen G. Study of nonlinear dynamic response of an unbalanced rotor supported on ball bearing. J Vib Acou, Trans. ASME 2009;131:06100119.

Inayat HJI, Kanki H, Mureithi NW. On the bifurcations of a rigid rotor response in squeezefilm dampers. J Flui Struc 2003;17:433 439.

Liang MA, JunHong Z, JiewWei L, Jun W, Xin L. Dynamic characteristic analysis of a misaligned rotorbearing system with squeezefilm dampers. J Zhej UniverSci 2016;17:614631.

Freeman JA, Skapura BM. Neural networks: Algorithms, Applications and Programming Techniques, AddisonWesley, NY. 1990.

Groves KH, Bonello P. Empirical identification of SFD bearings using neural networks. Mech Sys Sig Proc 2013;35:307323.