# Development and Experimental Validation of Auxiliary Rolling Bearing Models for Active Magnetic Bearings (AMBs) Applications

Authors: Anna Tangredi; Enrico Meli; Andrea Rindi; Alessandro Ridolfi; Pierluca D’Adamio; Amedeo Frilli; Duccio Fioravanti; Benjamin Defoy
Journal: International Journal of Rotating Machinery
 (2019)

---

## Abstract

Nowadays, the search for increasing performances in turbomachinery applications has led to a growing utilization of active magnetic bearings (AMBs), which can bring a series of advantages thanks to their features: AMBs allow the machine components to reach higher peripheral speeds; in fact there are no wear and lubrication problems as the contact between bearing surfaces is absent. Furthermore, AMBs characteristic parameters can be controlled via software, optimizing machine dynamics performances. However, active magnetic bearings present some peculiarities, as they have lower load capacity than the most commonly used rolling and hydrodynamic bearings, and they need an energy source; for these reasons, in case of AMBs overload or breakdown, an auxiliary bearing system is required to support the rotor during such landing events. During the turbomachine design process, it is fundamental to appropriately choose the auxiliary bearing type and characteristics, because such components have to resist to the rotor impact; so, a supporting design tool based on accurate and efficient models of auxiliary bearings is very useful for the design integration of the Active Magnetic Bearing System into the machine. This paper presents an innovative model to accurately describe the mechanical behavior of a complete rotor-dynamic system composed of a rotor equipped with two auxiliary rolling bearings. The model, developed and experimentally validated in collaboration withBaker Hughes a GE company (providing the test case and the experimental data), is able to reproduce the key physical phenomena experimentally observed; in particular, the most critical phenomenon noted during repeated experimental combined landing tests is the rotor forward whirl, which occurs in case of high friction conditions and greatly influences the whole system behavior. In order to carefully study some special phenomena like rotor coast down on landing bearings (which requires long period of time to evolve and involves many bodies and degrees of freedom) or other particular events like impacts (which occur in a short period of time), a compromise between accuracy of the results and numerical efficiency has been pursued. Some of the elements of the proposed model have been previously introduced in literature; however the present work proposes some new features of interest. For example, the lateral and the axial models have been properly coupled in order to correctly reproduce the effects observed during the experimental tests and a very important system element, the landing bearing compliant suspension, has been properly modelled to more accurately describe its elastic and damping effects on the system. Furthermore, the model is also useful to characterize the frequencies related to the rotor forward whirl motion.


---

## Body

## 1. Introduction

Rotors equipped with AMBs technology have the advantage to allow the machine to reach high performances, but a secondary rolling bearings system is required, because magnetic bearings are not able to support the rotor in case of overload or failure of the energy supply system. Such auxiliary bearings system, whose main elements are catcher bearings, is directly engaged during rotor landings; thus, in order to safely support the rotor, it must have very high mechanical and dynamical properties in order to guarantee the desired impact, thermal, and wear resistance. The landing system shall also be able to protect the rotor from possible damage and guarantee long service life; for these reasons, the development of a tool able to represent such rotor-dynamic systems and based on efficient rolling bearing models is useful to provide an appropriate support during the design phase.The current work presents an improved auxiliary bearing system model with a FEM rotor model able to represent real system dynamic behavior and to reproduce the key effects observed during experimental tests. The model is able to predict the behavior of the complete system during transient, steady and emergency conditions (e.g., delevitation and landing phenomena) and it is composed of different submodels: the rotor submodel, which includes a FEM model to represent the lateral dynamic of the rotor and two rigid body models to reproduce axial and torsional rotor dynamics, and the auxiliary rolling bearing submodel, characterized by a combination of multibody and contact models able to represent the dynamics and the interactions of the auxiliary system elements.The model has been developed and experimentally validated in collaboration withBaker Hughes a GE company which provided the system layout data and the results of the experimental tests performed on a dedicated scale test rig (Section 2.5).In this paper, a comparison between the most significant experimental and simulated results is presented; through the analysis of the results, it is possible to observe a satisfying reproduction of the real system behavior. Real rotor displacements, auxiliary bearing elements displacements, interaction forces and moments, whirling frequencies, rotor, and bearing speeds and acceleration are very well represented by the simulated results thanks to the optimal calibration of the model parameters.Rotating system behavior has been studied extensively in the past by many researchers. Szcygielski [1] conducted one of the initial studies concerning a high-speed rotor touching a rigid body, modelling the rotor like a gyro pendulum; the preliminary analytical results showed a good representation of the experimental trajectories of the gyro axis.Among many researchers, Muszynska [2], Beatty [3], and Ehrich [4] studied the rubbing phenomena, revealing that the rotating systems show nonlinear dynamical behavior such as sub- and super-synchronous responses and even chaotic motion.Other researchers investigated theoretically and experimentally the dry friction backward whirl motion; Bartha [5, 6] observed that, during rotor impact and rub phenomena, the rotor vibration shows backward orbiting due to the presence of a tangential friction force.Wilkes et al. [7] investigated, through experimental and numerical methods, the nature of multimode dry friction whip and whirl for a variety of rub materials and clearances; their experimental results showed multiple whirl and whip regions.In order to reduce rubbing and the related phenomena, Lahriri et al. [8] proposed new backup bearing design, using special pins to center the rotor during impacts and mitigate its lateral motion. Furthermore, they developed a mathematical model to capture phenomena arising from rotor impacts.Schmied and Pradetto [9] studied the behavior of a rotor landing onto auxiliary bearings and developed a numerical predictive model: simulation results reproduced a fully backward whirl rotor orbit for a few rotations, until the auxiliary bearing reached the rotor speed.Other researches focus on the steady-state responses of rotor-dynamic systems: H. Xie et al. [10] studied a complex rotor system supported by auxiliary bearings with clearance at each end of the rotor. The influence of rotor imbalance, clearance, support stiffness, and damping has been investigated using direct numerical integration of the governing equation of motion and the harmonic balance method.L. Della Pietra et al. [11] studied and experimentally investigated the static and the dynamic behavior of a rigid balanced rotor on journal bearings; the conducted study underlined the importance of the static behavior of the rotor bearing system for a correct interpretation of some aspect of the dynamic behavior. In fact they observed that a small dynamic imbalance can significantly modify the rotor-dynamic behavior inducing a whirl with large orbits.Other researchers designed innovative auxiliary bearings. Hawkins et al. [12] and McMullen et al. [13] have performed endurance tests, testing different variety of bearings and landing sleeves materials in order to improve the design of the auxiliary system. They observed excessive wearing of the thrust washer and concluded that a thrust surface with lower friction coefficient could provide better performance.M. Dougdag et al. [14] carried out an experimental verification of a simplified model of a ball bearing in both static and dynamic modes and test its capabilities to simulate accurately fault effects; the developed model considers the balls scrolling in the cage and the effect of the load rotating vector.Ransom et al. [15] and Masala et al. [16] performed experimental tests and developed specific numerical model on a scaled and a full scale AMB supported subsea compressor.R. G. Kirk [17] analysed AMB turbomachinery applications concerning the operation of the rotor drop onto the auxiliary bearings. His studies focused on the selection of specific auxiliary bearings able to withstand the contact forces following the speed rotor drop. The proper design of AMB machinery requires the calculation of the anticipated loading for the auxiliary bearings; so, Kirk presented results of transient response evaluation of a full-size compressor rotor to illustrate some of the important parameters in the design for rotor drop.In another research E.E. Swanson et al. [18] studied in detail the dynamic behavior of the rotor/bearing/support systems in the case of total loss of AMB support. They presented test results for the initial rotor drop transient for several backup bearing configurations. The cases reported are rigidly mounted antifriction bearing, rigidly mounted solid bushings (lubricated and unlubricated), and two soft mounted solid bushing configurations.Other researchers like J. Wilkes et al. [19] observed strong circular forward whirl following an AMB failure and developed specific models to reproduce such phenomenon. They developed an innovative model of auxiliary bearing system that represents a great step forward in the reproduction of rotor-dynamic system behavior; in fact the model allows reproducing the forward whirl rotor orbit observed during the landing tests. The model is based on a formulation where a friction force between the axial face of the rotor sleeve and the thrust face of the catcher bearing has been included.The model proposed in this research work represents an efficient and flexible tool, able to reproduce the behavior of the whole system and the physical phenomena connected to the auxiliary rolling bearings system (e.g., forward whirl). A compromise between accuracy of the results and numerical efficiency is mandatory, in order to carefully study some interesting phenomena like the rotor full coast down on landing bearings, which requires long period of time to evolve and involves many bodies and degrees of freedom, and particular events like impulsive forces which occur in a very short period of time.The present work proposes some innovative features of interest; the new formulation takes into account the reciprocal coupling between the lateral and the axial rotor-dynamic models, enabling the reproduction of the forward whirl observed in the radial dynamics and the oscillatory axial behavior of the system during the whirling phenomenon. The elastic and damping elements of the compliant axial and radial auxiliary bearing suspension have been modelled taking into account their features and their influence on the whole system behavior.The developed model can simulate the rotor dynamics in presence of different unbalances and represent the behavior of the whole system under different operating conditions: rotor constant speed, rotor decelerated speed, and different axial load profiles; it also includes the estimation of power and energy losses, so it can be used as an effective supporting tool for the machine designer to foresee the main sources of power losses. Furthermore the model is also useful to characterize the frequencies related to some characteristic phenomena of the system, like the backward and forward oscillation frequencies of the machine components.


## 2. Materials and Methods

In this Section the architecture of the rotor-dynamic system (Figure2) is described and the structure of the developed numerical model is presented; Section 2.3 is dedicated to the description of the auxiliary bearing system and the dynamic and contact models developed to represent the bearing elements behavior.


### 2.1. General Architecture

The rotor-dynamic system represented in the numerical model is composed of a rotor supported by two radial and one axial magnetic bearings and equipped with an auxiliary bearing system based on two rolling bearings (Figure1).Figure 1
General scheme of a magnetic bearings supported rotor with its auxiliary landing system.Figure 2
Model architecture.The main elements of the system, which are represented into the Matlab-Simulink model, are(1)
housing(a);(2)
rotor(b);(3)
two radial active magnetic bearings (AMBs)(c);(4)
one axial active magnetic bearing (AMB)(d);(5)
two auxiliary rolling bearings(e);(6)
radial rotor landing sleeve(f);(7)
axial and radial rotor landing sleeve(g);(8)
force actuator for the application of the axial load(h);(9)
baseplate(i);(10)
the axial support for the landing system(l);(11)
an elastic-damping element for each bearing, the ribbon damper(m).Between the rotor sleeve and the catcher bearings a radial clearance is present: in emergency conditions, the rotor passes through the radial clearance, lands on the inner rings of the catcher bearings, and continues to rotate supported by the auxiliary system.The real rotor-dynamic system is equipped with a force actuator which applies an external axial thrust load on the nondrive end of the system; at this side of the system there is also an axial clearance between the rotor sleeve and the auxiliary bearing; under the axial thrust, the rotor translates along its symmetry axis.Between the outer race and the housing specific elastic elements (Figure1) are placed: the axial support for the bearing and the elastic ribbon damper. The axial support acts on the outer race of the bearing, while the ribbon is placed between the outer ring of the auxiliary bearing and the housing; it creates both radial and tangential force components and has both elastic and damping characteristics. In the event of a full circular whirl mode, the ribbon characteristics greatly influence the frequency of the whirl motion.The whole model is able to predict the behavior of the complete rotor-dynamic system during transient, steady and emergency conditions (e.g., emergency landing phenomena) and it is composed of two main submodels (see Figure2):(1)
Rotor model, which includes a FEM model representing the lateral dynamics of the rotor and two simplified models, where the rotor is represented as a rigid body, to reproduce axial and torsional rotor dynamics; the main purpose is to have a good compromise between accuracy of the results and computational efficiency, considering the long time duration of the events to be simulated.(2)
Auxiliary rolling bearing model, characterized by a combination of multibody and contact models, able to represent the dynamics of the elements of the auxiliary system.In order to reproduce the forward whirl phenomenon observed during combined radial-axial landing experimental tests, the lateral and the axial model have been coupled (see Figure2): the axial force calculated in the rotor axial model represents an input for the catcher bearing lateral model.The effect produced by this coupling allows well reproducing the phenomena observed during the experimental tests and, especially, the forward whirling.


### 2.2. Rotor Model

The lateral dynamic of the rotor is described by a FEM model, while the axial and the torsional dynamics are described by two simplified models (Figure2) where the rotor is represented as a rigid body.The kinematic quantities calculated by the rotor model (e.g., rotor nodes radial and axial displacementsy, z, α, and β, rotor nodes velocities y˙, z˙, α˙, and β˙, and rotor angular velocity Ω) represent the inputs of the auxiliary rolling bearing model; the interaction forces and moments (e.g., radial interaction forces between the rotor and the bearing Fr, interaction torques T, and axial forces Fax) between the bearing elements are calculated in the auxiliary bearing model and constitute the inputs of the rotor model. So, there is a reciprocal exchange of kinematic and dynamic quantities between the submodels at each time step of the simulation.During the experimental tests carried out on a dedicated test rig, in particular during combined radial-axial landings, coupling effects between the lateral and the axial dynamics have been observed. Under an axial thrust load the experiments show a forward whirl rotor lateral orbit; in order to represent such effects, a specific coupled formulation between the two dynamic models has been developed.The whole model provides a good representation of the dynamical behavior of the rotor with a low computational effort.


#### 2.2.1. Rotor Lateral Model

In order to describe the lateral dynamics, the rotor is represented through a finite element discretization, where each rotor element is considered as a 2D BEAM (Figure3). Such FEM model simulates the lateral rotor behavior, according to the classical rotor-dynamics equations of motion and accounting for magnetic and auxiliary bearings actions and external loads [20]; each BEAM element is based on the Timoshenko Theory [20, 21]. Each node has 4 degrees of freedom: two radial translations and two rotations around the translational axes.Figure 3
Discretization of the rotor in 2D BEAM elements.The inputs of the model are the forces and moments acting on the rotor in correspondence of the nodes where the magnetic and the auxiliary bearings are located and the external loads and torques (i.e., the unbalance load). The outputs of the rotor model are the kinematic quantities: positionsx and velocities x˙ for each degree of freedom of the rotor nodes.The equations of motion that represent the rotor lateral behavior are reported below [20]:(1)Mrotx¨+Crotx˙+ΩGrotx˙+Krotx=F,where x contains the two translational and the two rotational degrees of freedom of the rotor nodes (y, z, α, β), while Mrot, Crot, Grot, and Krot are, respectively, the mass, damping, stiffness, and gyroscopic matrices of the rotor; the vector F contains the external forces and moments acting on the rotor nodes.Furthermore, a modal reduction for the rotor FEM lateral model has been implemented. The modal reduction used is based on the Craig-Bampton method and generates reduced rotor matrices (mass, stiffness, damping, and gyroscopic matrix) on the base of the number of nodes and modes chosen by the user to describe the considered system [22]; by using the modal reduction, it is possible to simulate long period of time in quite short simulation time, reaching a good compromise between results accuracy and computational efficiency.


#### 2.2.2. Rotor Axial Model

In order to reproduce the axial rotor dynamics, a simple model has been developed in which the rotor is represented as a rigid body able to translate along its symmetry axis. The model describes the axial vibration of the rotor and the axial force, which are input for the lateral model and must be known at each instant of the simulation in order to have a real coupling between the lateral and the axial dynamics.Figure4 shows a scheme of the system axial model; the kax, cax parameters are the equivalent axial stiffness and damping of the system which have variable values depending on the rotor axial displacement; Fthrust is the thrust force applied by the actuator and is used to simulate the axial landing, while mr is the rotor mass. In correspondence of the initial instant of the simulation, the rotor is placed at its initial position (xr(0)=0) where the axial clearance is equal to clax; when the axial thrust force is applied to the system, the rotor moves according to the following dynamic equations:(2)x¨r=-Fthrustmr+Faxt,where Fthrust is the axial thrust force acting on the rotor and mr is the rotor mass.Figure 4
Schematic representation of the system axial model.By knowingxr(t) and the equivalent stiffness kax and damping cax of the system (Figure 4), the axial reaction force Fax can be calculated as(3)Faxt=-kaxxrt-caxx˙rt.The axial forceFax(t) depends on the axial thrust load force Fthrust and acts on the lateral surface of the rotor sleeve and on the thrust surface of the inner race of the bearing. The axial contact between the rotor sleeve and the inner race of the rolling bearing produces a tangential force that is proportional to Fax(t) and the friction coefficient; for this reason, it acts on the lateral degrees of freedom of the rotor and of the bearing inner race and it represents the coupling between the lateral and the axial dynamic model: the axial force Fax, calculated through (3), is the main input of the rotor lateral model and of the rolling bearing inner race lateral model.The equivalent axial stiffnesskax of the system can be calculated considering the stiffness elements (bearing stiffness kcb, support stiffness kws, and housing stiffness kh) and taking into account the axial displacement of the rotor: when the axial system displacement is less than the axial clearance, no stiffness is considered to work; as the system translation exceeds the axial clearance clax, the support stiffness kws, and the bearing stiffness kcb work together; so, before the axial displacement of the shaft reaches the value δm+clax corresponding to the maximum deflection of the support, the equivalent axial stiffness can be obtained through the series of kcb and kws and the axial force Fax(t) can be calculated as follows:(4)ifxr≤claxFaxt=0ifxr>clax&xr≤clax+δmFax=-1kcb+1kws-1xr-clax-caxx˙rifxr<-clax&xr≥-clax+δmFax=-1kcb+1kws-1xr+clax-caxx˙r. where xr represents the axial coordinate of the rotor. When the axial displacement of the shaft exceeds the value δm+clax (and the maximum deflection of the support is reached), the housing and the bearing stiffnesses kh and kcb work in series and the support stiffness kws does not play any role anymore:(5)ifxr>clax+δmFax=-1kcb+1kh-1xr-clax+δm-1kcb+1kws-1δm-caxx˙rifxr<-clax+δmFax=-1kcb+1kh-1xr+clax+δm+1kcb+1kws-1δm-caxx˙r.


#### 2.2.3. Rotor Spinning Model

In order to determine the rotor rotational velocity during the simulation, an approximate torsional model has been implemented, in which the rotor is described as a rigid body able to rotate around its symmetry axis. This simplified model allows obtaining a good representation of the torsional dynamics of the rotor, ensuring again a good compromise between accuracy and numerical efficiency.The torques acting on the rotor are the auxiliary bearing moments and the windage torques (see (6)): (6)IΩ˙=Taux=Taux1+Taux2+Twind,whereI represents the torsional moment of inertia and Ω˙ is the rotor angular speed; Taux1 and Taux2 are the torques exerted by the two auxiliary bearings on the rotor while Twind is the windage torque.The rotor angular speed decreases, during a simulation, because of friction effects (see Figure5).Figure 5
Scheme of the rotor motion.


### 2.3. Auxiliary Bearing Model

The catcher bearing model describes the dynamics of the bearing components and allows the calculation of the interaction forces and moments between their contact surfaces. By solving the equations of motion of each element of the auxiliary bearing system, it is possible to calculate the position of the center of gravity of each element of the catcher bearing.The catcher bearing system (diameter of the bearing equal to 60 mm) is shown in Figure6 and is composed of the following elements:(1)
housing;(2)
ribbon damper;(3)
axial support;(4)
inner race;(5)
outer races;(6)
rolling elements;(7)
landing sleeve.Figure 6
Schematic representation of the auxiliary system.From a dynamical point of view, a simplifying assumption has been made: the catcher bearing elements are able to move in the plane, so the bearing inner race has 3 degrees of freedom (i.e., it can translate along they and z directions and rotate around the x axis) while the outer race can only translate in the plane without rotating.Equations (7) describe the translational dynamics of the rolling bearing inner race; the bearing inner race degrees of freedom are two because the inner race motion is considered to move in the radial plane of the bearing OYZ (Figure 5). This simplifying hypothesis allows obtaining a good compromise between the accuracy computational efficiency:(7)Mix¨i+Cix˙i=P+F,F=Fio+Fir+Frub,where xi contains the coordinates of the inner race center of gravity, Mi is the inner race mass matrix, Ci is the inner race damping matrix, P is the weight force of the inner race, and the vector F is the sum of the radial forces acting on the bearing inner race (see also Figure 7): Fio is the interaction force on the inner race exerted by the rolling elements, Fir is the contact force exchanged between the inner race and the rotor during the contact (including a normal FrN component and a tangential one Ffr; see (15)), and Frub is the resultant force generated during the axial rub between the thrust surface of the bearing inner race and the lateral surface of the rotor sleeve (see (19)).Figure 7
Schematic representation of the interactions between the elements of the auxiliary system.Equation (8) describes the rotational dynamics of the bearing inner race:(8)IiΩ˙i=Tfr+Tcb+Trubwhere Ii is the torsional moment of inertia of the bearing inner race and Tfr is the moment acting on the inner race and due to the action of the radial friction force Ffr; Tcb is the ball bearing drag torque due to the rolling friction of the rolling elements during the motion and Trub is the resultant torque due to the axial rub between the thrust surface of the bearing inner race and the lateral surface of the rotor sleeve.The outer surface of the outer race interacts with the housing through the ribbon damper, which is an elastic-damping element; furthermore, there is another elastic element: the axial support on the outer race, placed between the external surface of the outer race and the statoric flange.Equation (9) describes the translational dynamics of the bearing outer race, which can move into the radial plane OYZ (Figure 5):(9)Mox¨o+Cox˙o=Foi+Foh,where xo contains the coordinates of the outer race center of gravity,Mo and Co are the bearing outer race mass and damping matrices, the vector Foi is the interaction force on the outer race exerted by the rolling elements, and Foh is the radial force exchanged between the outer race and the housing flange.The auxiliary bearings have spherical rolling elements, which have been modelled as force elements, without representing their dynamic behavior. The elastic force exerted by the rolling elements allows representing the interactions between the outer and the inner race of the bearings.


### 2.4. Contact Models

The catcher bearing model simulates the plane motion of the auxiliary bearing elements and calculates the possible interaction forces and moments among them; the kinematic physical quantities are determined by solving the equations of motion of each elements (see Section2.3), while the interaction forces/moments are calculated through appropriate contact models based on the Hertz theory.In this section, the contact models implemented in the proposed bearing model will be explained, showing the formulation used to calculate the interaction forces and moments between the contact surfaces.The surfaces which interact during the motion are (see Figure7):(1)
radial contact between the rotor landing sleeve and the inner surface of the bearing inner race: the contact forceFir exchanged between the inner race and the rotor during the contact includes a normal FrN component and a tangential one Ffr; furthermore, a friction moment Tfr is generated;(2)
axial contact between the thrust surface of the rolling bearing inner race and the lateral side of the landing rotor sleeve: during axial or combined axial-radial rotor landings, axial rubbing force and moment are generatedFrub, Trub;(3)
elastic interaction between the bearing outer race and the bearing inner race through the rolling elements: the balls exert an elastic forceFio on the inner race and Foi on the outer races;(4)
the interaction force between the housing and the outer race is exerted through the ribbon damper elementFoh.


#### 2.4.1. Rotor-Inner Race Radial Interaction

Concerning the development of this part of the model, see also the paper [19]. During a radial landing, the inner surface of the bearing inner ring and the rotor landing sleeve get in contact; the rotor bounces, after the impact, on the inner ring and continues to rotate supported by the auxiliary bearings. As a result of the interaction with the rotor, the inner ring accelerates and tends to reach the rotor speed.The contact forcesFrN exchanged between the rotor and the bearing inner ring is made up of two components: FrN is the normal contact force and Ffr is the friction force due to the radial rubbing.The mathematical steps required to calculate such forces are briefly reported below.It is very useful to define a radial-tangential (R(i1),T(j1)) reference system and express the quantities of interest (e.g., bodies centers of gravity coordinates; velocities) in this specific system (Figure 8).Figure 8
Interaction normal and friction forces between the rotor and the bearing inner ring.Furthermore, it is convenient to refer the analysis to the relative displacements and velocities between the rotor and the inner race in radial and tangential coordinates. The relative displacementsdri between the two bodies can be obtained from the difference between the coordinates of the rotor node Or and the coordinates of the inner race center of gravity Oi:(10)dri=Or-Oi.Through the use of the relative coordinates, it is possible to define a relative angle of rotationσrel and, consequently, to determine the rotational matrix R needed to convert the quantities expressed in the Cartesian reference system (Z(i0),Y(j0)) into the radial-tangential reference system (R(i1),T(j1)).Equation (11) defines the relative angle of rotation σrel:(11)σrel=atan⁡2yri,zri,while (12) defines the transformation matrix R:(12)R=cos⁡σrelsin⁡σrelsin⁡σrel-cos⁡σrel.The relative displacements dri and velocities vri expressed in the Cartesian reference system can be expressed into the radial-tangential reference system using the rotational matrix R:(13)drt=RTdri,vrt=RTvri.The contact forces and moments are calculated through the contact model if there is a interpenetrationδri between the two bodies. By knowing the physical value of the radial clearance, the penetration between the two bodies along the radial direction can be calculated as follows:(14)δri=yri-cr.Using the interpenetrationδri, it is possible to calculate the normal contact force FrN and the friction force Ffr due to the radial rubbing between the rotor sleeve radial surface and the inner surface of the inner ring:(15)FrN=Khertzδriγ+Kviscδ˙riδrii1,Ffr=μradFrNSj1where Khertz is the Hertzian contact constant, Kvisc is the viscous contact constant, μrad is the radial friction coefficient, and S is the sign function of the tangential component of the relative velocity vt:(16)vt=vrtj1+rrΩrk1-riΩik1.The friction force Ffr due to the radial rubbing between the rotor sleeve and the inner ring generates a torque Tfr responsible for the acceleration of the inner race:(17)Tfr=ri+rr2Ffrk1.During a pure radial landing, the rotor falls and bounces on to the inner ring of the bearing, which begins to rotate and quickly reaches the rotor angular speed; in case of high radial friction, the rotor performs a backward whirl orbit because of the friction effects between the rotor and the bearing inner ring.On the other hand, during a combined axial-radial rotor landing due to AMBs failure (both radial and axial magnetic bearings), the rotor gets in contact with the inner surface of the inner ring and the thrust face of the catcher bearing. In case of high friction effects, a circular forward whirl rotor orbit has been observed in some experimental tests; in fact, as a consequence of repeated rotor landings, the bearing elements begin to be worn and the friction coefficient increases; the tangential actions due to that axial friction result to be stronger than the radial rubbing force.


#### 2.4.2. Rotor Sleeve-Inner Race Axial Interaction

As regards the development of this part of the model, see also the paper [19]. Concerning the lateral dynamics, the rotor can perform small rotations around Y and Z axes; as a consequence of the rotor dynamics, the interpenetration between the lateral face of the rotor sleeve and the thrust surface of the inner ring cannot have a constant value; by knowing the inclination angles α and β of the rotor at any instant of time, it is possible to calculate the interpenetration value δri at each point of the thrust surface and, consequently, calculate the axial pressure Paxial at any point on the contact surface.In order to calculate the resultant force and the resultant torque acting on the rotor due to the friction at the thrust surface, it is necessary to evaluate the relative slip velocity between the same pointA(r,θ) on the rotor and on the bearing.The vectorv˙A,r represents the velocity of the point A(r,θ) on the rotor, while the vector v˙A,i represents the velocity of the point A(r,θ) on the bearing inner ring:(18)vA,rrA,θA=r˙r-ω˙rrAsin⁡θAi1+t˙r-ω˙rrri-rAcos⁡θAj1,vA,irA,θA=r˙i-ω˙irAsin⁡θAi1+t˙i-ω˙irAcos⁡θAj1.The relative velocity vA,ri between the point A on the rotor and the same point A on the bearing can be simply calculated by subtracting the inner ring velocity vA,r from the rotor velocity vA,i.By knowing the relative velocity between the rotor and the bearing, the resultant forceFrub, acting on the rotor and due to friction at the contact interface, can be evaluated by integrating the pressure field over the contact area:(19)Frub,r=-μax∫r1r2∫02πPA,rr˙ri1+t˙rj1vA,rirAdθAdRA,where μax is the friction coefficient between the lateral face of the rotor sleeve and the thrust face of the bearing.The torque due to the axial contact between the rotor and the bearing inner race can be determined as follows:(20)Trub,r=-μax∫r1r2∫02πPA,rt˙rcos⁡θA+r˙rsin⁡θAvA,rirA2dθAdRAk1.The torque on the rotor and on the bearing is considered to be equal and opposite:(21)Trub,r=-Trub,i.


#### 2.4.3. Outer Race-Inner Race Radial Interaction

The interaction forceFio between the bearing inner and outer races is the radial reaction due to the relative motion of the races.The races can interact through the rolling elements, which are modelled as force elements; through the load-deflection curve of the catcher bearing it is possible to extrapolate the stiffness characteristics of the rolling elements (kb1 and kb2) and define the interaction force Fio with a second order polynomial function:(22)Fio=kb2rio2+kb1rioi1,where rio is the relative displacement between the inner and the outer races centers of gravity.The rolling elements exert also another important action on the catcher bearing: the rolling friction due to the balls generates a drag torque onto the inner ring, which decelerates under this effect. The approach to calculate the resistant torqueTCB is based on the Harris-Palmgren formulation [23], but, instead of considering the applied static load, the dynamic radial and axial forces acting on the bearing have been taken into account:(23)TCB=f1Fβdmk1,where f1 is a factor depending on the bearing design and relative bearing load, Fβ depends on the magnitude and direction of the applied load, and dm is the pitch diameter of the bearing.For a ball bearing, Harris [23] defines f1 as(24)f1=zFsyCs,where Fs represents the equivalent dynamic load and Cs is the basic static load rating usually provided in the manufacturer catalog, while z and y represent empirical coefficients.For a single row angular contact bearing the equivalent loadFs can be defined as follows:(25)Fs=XsFradial+YsFax,where Xs=0.5 and Ys varies with the bearing contact angle α:(26)Ys=c1α2+c2α+c3.For a radial ball bearing Fβ is given by (27)Fβ=c4Faxcot⁡α+c5Fradial.The numeric parametersc1, c2, c3, c4, and c5 of (26) and (27) have been determined on the base of the bearing geometry and characteristics: c1=7.142910-5, c2=-1.244310-2, c3=0.64129, c4=0.9, and c5=-0.1.


#### 2.4.4. Outer Race-Housing Interaction

The ribbon damper is located between the bearing outer ring and the housing; it has both elastic and damping characteristics (i.e., structural and viscous damping) and generates radial and tangential force components. The characteristics of the ribbon greatly influence the frequency of the rotor whirl motion.The interaction force between the outer race and the housingFoh can be expressed through the force exerted by the ribbon:(28)Foh=-Cribx˙o-KribQxo-KribIxo,where I is the identity matrix, xo is the vector of the coordinates of the outer race center of gravity, Crib is the ribbon viscous damping, Krib=6.3∗107N/m is the ribbon stiffness, and Q is a matrix defined as follows:(29)Q=0-ηSvηSv0.The parameter η is a cross coupled stiffness factor, while Sv is the sign of the relative tangential velocity between the outer bearing and the housing.For combined axial-radial rotor landings, because of the axial thrust force applied to the bearing inner race, a reaction force exerted by the housing on the outer race through the axial support is present. The force contribution due to this interaction is the damping forceFoh:(30)Foh=-Faxμdx˙ox˙o,where μd is the Coulomb friction coefficient between the outer race and the housing and must be added to (28).


### 2.5. Test Rig Architecture

Figure9 shows an overview of the scaled test rig considered for the development and the validation of the proposed model. The main components of the test rig are highlighted in the following:(1)
electric motor(1);(2)
power transmission coupling(2);(3)
force sensor(3);(4)
radial AMBs and auxiliary bearing(4); they are concentric;(5)
axial centering AMB(5);(6)
rotor(6);(7)
axial actuator(7).Figure 9
Experimental apparatus; from left to right: axial actuator, NDE radial magnetic bearing, speed sensor, DE radial magnetic bearing, magnetic thrust bearing, and power transmission coupling.Figure10 shows a scheme of the test rig that is equipped with a sensors system able to estimate components temperature, forces, acceleration, vibrations, and speeds:(1)
load cells(LC);(2)
proximity sensors(PS);(3)
encoders(EN);(4)
accelerometers(A);(5)
thermocouples(TC).Figure 10
Scheme of the experimental apparatus.The tests performed on the test rig allow the identification of the system behavior in landing conditions. In particular the tests carried out on the test rig belong to three different types (with different rotor speeds): radial landings (deactivating just radial AMBs), axial landings (deactivating just axial AMB), and combined radial-axial landings (deactivating both the AMBs).The experimental tests taken as a reference in this paper for the model validation (Section3) are no whirling rotor radial landing, a pure radial landing characterized by a normal friction condition; whirling rotor radial landing, a pure radial landing with high friction coefficient; whirling rotor combining axial-radial landing with high axial thrust force value.


## 2.1. General Architecture

The rotor-dynamic system represented in the numerical model is composed of a rotor supported by two radial and one axial magnetic bearings and equipped with an auxiliary bearing system based on two rolling bearings (Figure1).Figure 1
General scheme of a magnetic bearings supported rotor with its auxiliary landing system.Figure 2
Model architecture.The main elements of the system, which are represented into the Matlab-Simulink model, are(1)
housing(a);(2)
rotor(b);(3)
two radial active magnetic bearings (AMBs)(c);(4)
one axial active magnetic bearing (AMB)(d);(5)
two auxiliary rolling bearings(e);(6)
radial rotor landing sleeve(f);(7)
axial and radial rotor landing sleeve(g);(8)
force actuator for the application of the axial load(h);(9)
baseplate(i);(10)
the axial support for the landing system(l);(11)
an elastic-damping element for each bearing, the ribbon damper(m).Between the rotor sleeve and the catcher bearings a radial clearance is present: in emergency conditions, the rotor passes through the radial clearance, lands on the inner rings of the catcher bearings, and continues to rotate supported by the auxiliary system.The real rotor-dynamic system is equipped with a force actuator which applies an external axial thrust load on the nondrive end of the system; at this side of the system there is also an axial clearance between the rotor sleeve and the auxiliary bearing; under the axial thrust, the rotor translates along its symmetry axis.Between the outer race and the housing specific elastic elements (Figure1) are placed: the axial support for the bearing and the elastic ribbon damper. The axial support acts on the outer race of the bearing, while the ribbon is placed between the outer ring of the auxiliary bearing and the housing; it creates both radial and tangential force components and has both elastic and damping characteristics. In the event of a full circular whirl mode, the ribbon characteristics greatly influence the frequency of the whirl motion.The whole model is able to predict the behavior of the complete rotor-dynamic system during transient, steady and emergency conditions (e.g., emergency landing phenomena) and it is composed of two main submodels (see Figure2):(1)
Rotor model, which includes a FEM model representing the lateral dynamics of the rotor and two simplified models, where the rotor is represented as a rigid body, to reproduce axial and torsional rotor dynamics; the main purpose is to have a good compromise between accuracy of the results and computational efficiency, considering the long time duration of the events to be simulated.(2)
Auxiliary rolling bearing model, characterized by a combination of multibody and contact models, able to represent the dynamics of the elements of the auxiliary system.In order to reproduce the forward whirl phenomenon observed during combined radial-axial landing experimental tests, the lateral and the axial model have been coupled (see Figure2): the axial force calculated in the rotor axial model represents an input for the catcher bearing lateral model.The effect produced by this coupling allows well reproducing the phenomena observed during the experimental tests and, especially, the forward whirling.


## 2.2. Rotor Model

The lateral dynamic of the rotor is described by a FEM model, while the axial and the torsional dynamics are described by two simplified models (Figure2) where the rotor is represented as a rigid body.The kinematic quantities calculated by the rotor model (e.g., rotor nodes radial and axial displacementsy, z, α, and β, rotor nodes velocities y˙, z˙, α˙, and β˙, and rotor angular velocity Ω) represent the inputs of the auxiliary rolling bearing model; the interaction forces and moments (e.g., radial interaction forces between the rotor and the bearing Fr, interaction torques T, and axial forces Fax) between the bearing elements are calculated in the auxiliary bearing model and constitute the inputs of the rotor model. So, there is a reciprocal exchange of kinematic and dynamic quantities between the submodels at each time step of the simulation.During the experimental tests carried out on a dedicated test rig, in particular during combined radial-axial landings, coupling effects between the lateral and the axial dynamics have been observed. Under an axial thrust load the experiments show a forward whirl rotor lateral orbit; in order to represent such effects, a specific coupled formulation between the two dynamic models has been developed.The whole model provides a good representation of the dynamical behavior of the rotor with a low computational effort.


### 2.2.1. Rotor Lateral Model

In order to describe the lateral dynamics, the rotor is represented through a finite element discretization, where each rotor element is considered as a 2D BEAM (Figure3). Such FEM model simulates the lateral rotor behavior, according to the classical rotor-dynamics equations of motion and accounting for magnetic and auxiliary bearings actions and external loads [20]; each BEAM element is based on the Timoshenko Theory [20, 21]. Each node has 4 degrees of freedom: two radial translations and two rotations around the translational axes.Figure 3
Discretization of the rotor in 2D BEAM elements.The inputs of the model are the forces and moments acting on the rotor in correspondence of the nodes where the magnetic and the auxiliary bearings are located and the external loads and torques (i.e., the unbalance load). The outputs of the rotor model are the kinematic quantities: positionsx and velocities x˙ for each degree of freedom of the rotor nodes.The equations of motion that represent the rotor lateral behavior are reported below [20]:(1)Mrotx¨+Crotx˙+ΩGrotx˙+Krotx=F,where x contains the two translational and the two rotational degrees of freedom of the rotor nodes (y, z, α, β), while Mrot, Crot, Grot, and Krot are, respectively, the mass, damping, stiffness, and gyroscopic matrices of the rotor; the vector F contains the external forces and moments acting on the rotor nodes.Furthermore, a modal reduction for the rotor FEM lateral model has been implemented. The modal reduction used is based on the Craig-Bampton method and generates reduced rotor matrices (mass, stiffness, damping, and gyroscopic matrix) on the base of the number of nodes and modes chosen by the user to describe the considered system [22]; by using the modal reduction, it is possible to simulate long period of time in quite short simulation time, reaching a good compromise between results accuracy and computational efficiency.


### 2.2.2. Rotor Axial Model

In order to reproduce the axial rotor dynamics, a simple model has been developed in which the rotor is represented as a rigid body able to translate along its symmetry axis. The model describes the axial vibration of the rotor and the axial force, which are input for the lateral model and must be known at each instant of the simulation in order to have a real coupling between the lateral and the axial dynamics.Figure4 shows a scheme of the system axial model; the kax, cax parameters are the equivalent axial stiffness and damping of the system which have variable values depending on the rotor axial displacement; Fthrust is the thrust force applied by the actuator and is used to simulate the axial landing, while mr is the rotor mass. In correspondence of the initial instant of the simulation, the rotor is placed at its initial position (xr(0)=0) where the axial clearance is equal to clax; when the axial thrust force is applied to the system, the rotor moves according to the following dynamic equations:(2)x¨r=-Fthrustmr+Faxt,where Fthrust is the axial thrust force acting on the rotor and mr is the rotor mass.Figure 4
Schematic representation of the system axial model.By knowingxr(t) and the equivalent stiffness kax and damping cax of the system (Figure 4), the axial reaction force Fax can be calculated as(3)Faxt=-kaxxrt-caxx˙rt.The axial forceFax(t) depends on the axial thrust load force Fthrust and acts on the lateral surface of the rotor sleeve and on the thrust surface of the inner race of the bearing. The axial contact between the rotor sleeve and the inner race of the rolling bearing produces a tangential force that is proportional to Fax(t) and the friction coefficient; for this reason, it acts on the lateral degrees of freedom of the rotor and of the bearing inner race and it represents the coupling between the lateral and the axial dynamic model: the axial force Fax, calculated through (3), is the main input of the rotor lateral model and of the rolling bearing inner race lateral model.The equivalent axial stiffnesskax of the system can be calculated considering the stiffness elements (bearing stiffness kcb, support stiffness kws, and housing stiffness kh) and taking into account the axial displacement of the rotor: when the axial system displacement is less than the axial clearance, no stiffness is considered to work; as the system translation exceeds the axial clearance clax, the support stiffness kws, and the bearing stiffness kcb work together; so, before the axial displacement of the shaft reaches the value δm+clax corresponding to the maximum deflection of the support, the equivalent axial stiffness can be obtained through the series of kcb and kws and the axial force Fax(t) can be calculated as follows:(4)ifxr≤claxFaxt=0ifxr>clax&xr≤clax+δmFax=-1kcb+1kws-1xr-clax-caxx˙rifxr<-clax&xr≥-clax+δmFax=-1kcb+1kws-1xr+clax-caxx˙r. where xr represents the axial coordinate of the rotor. When the axial displacement of the shaft exceeds the value δm+clax (and the maximum deflection of the support is reached), the housing and the bearing stiffnesses kh and kcb work in series and the support stiffness kws does not play any role anymore:(5)ifxr>clax+δmFax=-1kcb+1kh-1xr-clax+δm-1kcb+1kws-1δm-caxx˙rifxr<-clax+δmFax=-1kcb+1kh-1xr+clax+δm+1kcb+1kws-1δm-caxx˙r.


### 2.2.3. Rotor Spinning Model

In order to determine the rotor rotational velocity during the simulation, an approximate torsional model has been implemented, in which the rotor is described as a rigid body able to rotate around its symmetry axis. This simplified model allows obtaining a good representation of the torsional dynamics of the rotor, ensuring again a good compromise between accuracy and numerical efficiency.The torques acting on the rotor are the auxiliary bearing moments and the windage torques (see (6)): (6)IΩ˙=Taux=Taux1+Taux2+Twind,whereI represents the torsional moment of inertia and Ω˙ is the rotor angular speed; Taux1 and Taux2 are the torques exerted by the two auxiliary bearings on the rotor while Twind is the windage torque.The rotor angular speed decreases, during a simulation, because of friction effects (see Figure5).Figure 5
Scheme of the rotor motion.


## 2.2.1. Rotor Lateral Model

In order to describe the lateral dynamics, the rotor is represented through a finite element discretization, where each rotor element is considered as a 2D BEAM (Figure3). Such FEM model simulates the lateral rotor behavior, according to the classical rotor-dynamics equations of motion and accounting for magnetic and auxiliary bearings actions and external loads [20]; each BEAM element is based on the Timoshenko Theory [20, 21]. Each node has 4 degrees of freedom: two radial translations and two rotations around the translational axes.Figure 3
Discretization of the rotor in 2D BEAM elements.The inputs of the model are the forces and moments acting on the rotor in correspondence of the nodes where the magnetic and the auxiliary bearings are located and the external loads and torques (i.e., the unbalance load). The outputs of the rotor model are the kinematic quantities: positionsx and velocities x˙ for each degree of freedom of the rotor nodes.The equations of motion that represent the rotor lateral behavior are reported below [20]:(1)Mrotx¨+Crotx˙+ΩGrotx˙+Krotx=F,where x contains the two translational and the two rotational degrees of freedom of the rotor nodes (y, z, α, β), while Mrot, Crot, Grot, and Krot are, respectively, the mass, damping, stiffness, and gyroscopic matrices of the rotor; the vector F contains the external forces and moments acting on the rotor nodes.Furthermore, a modal reduction for the rotor FEM lateral model has been implemented. The modal reduction used is based on the Craig-Bampton method and generates reduced rotor matrices (mass, stiffness, damping, and gyroscopic matrix) on the base of the number of nodes and modes chosen by the user to describe the considered system [22]; by using the modal reduction, it is possible to simulate long period of time in quite short simulation time, reaching a good compromise between results accuracy and computational efficiency.


## 2.2.2. Rotor Axial Model

In order to reproduce the axial rotor dynamics, a simple model has been developed in which the rotor is represented as a rigid body able to translate along its symmetry axis. The model describes the axial vibration of the rotor and the axial force, which are input for the lateral model and must be known at each instant of the simulation in order to have a real coupling between the lateral and the axial dynamics.Figure4 shows a scheme of the system axial model; the kax, cax parameters are the equivalent axial stiffness and damping of the system which have variable values depending on the rotor axial displacement; Fthrust is the thrust force applied by the actuator and is used to simulate the axial landing, while mr is the rotor mass. In correspondence of the initial instant of the simulation, the rotor is placed at its initial position (xr(0)=0) where the axial clearance is equal to clax; when the axial thrust force is applied to the system, the rotor moves according to the following dynamic equations:(2)x¨r=-Fthrustmr+Faxt,where Fthrust is the axial thrust force acting on the rotor and mr is the rotor mass.Figure 4
Schematic representation of the system axial model.By knowingxr(t) and the equivalent stiffness kax and damping cax of the system (Figure 4), the axial reaction force Fax can be calculated as(3)Faxt=-kaxxrt-caxx˙rt.The axial forceFax(t) depends on the axial thrust load force Fthrust and acts on the lateral surface of the rotor sleeve and on the thrust surface of the inner race of the bearing. The axial contact between the rotor sleeve and the inner race of the rolling bearing produces a tangential force that is proportional to Fax(t) and the friction coefficient; for this reason, it acts on the lateral degrees of freedom of the rotor and of the bearing inner race and it represents the coupling between the lateral and the axial dynamic model: the axial force Fax, calculated through (3), is the main input of the rotor lateral model and of the rolling bearing inner race lateral model.The equivalent axial stiffnesskax of the system can be calculated considering the stiffness elements (bearing stiffness kcb, support stiffness kws, and housing stiffness kh) and taking into account the axial displacement of the rotor: when the axial system displacement is less than the axial clearance, no stiffness is considered to work; as the system translation exceeds the axial clearance clax, the support stiffness kws, and the bearing stiffness kcb work together; so, before the axial displacement of the shaft reaches the value δm+clax corresponding to the maximum deflection of the support, the equivalent axial stiffness can be obtained through the series of kcb and kws and the axial force Fax(t) can be calculated as follows:(4)ifxr≤claxFaxt=0ifxr>clax&xr≤clax+δmFax=-1kcb+1kws-1xr-clax-caxx˙rifxr<-clax&xr≥-clax+δmFax=-1kcb+1kws-1xr+clax-caxx˙r. where xr represents the axial coordinate of the rotor. When the axial displacement of the shaft exceeds the value δm+clax (and the maximum deflection of the support is reached), the housing and the bearing stiffnesses kh and kcb work in series and the support stiffness kws does not play any role anymore:(5)ifxr>clax+δmFax=-1kcb+1kh-1xr-clax+δm-1kcb+1kws-1δm-caxx˙rifxr<-clax+δmFax=-1kcb+1kh-1xr+clax+δm+1kcb+1kws-1δm-caxx˙r.


## 2.2.3. Rotor Spinning Model

In order to determine the rotor rotational velocity during the simulation, an approximate torsional model has been implemented, in which the rotor is described as a rigid body able to rotate around its symmetry axis. This simplified model allows obtaining a good representation of the torsional dynamics of the rotor, ensuring again a good compromise between accuracy and numerical efficiency.The torques acting on the rotor are the auxiliary bearing moments and the windage torques (see (6)): (6)IΩ˙=Taux=Taux1+Taux2+Twind,whereI represents the torsional moment of inertia and Ω˙ is the rotor angular speed; Taux1 and Taux2 are the torques exerted by the two auxiliary bearings on the rotor while Twind is the windage torque.The rotor angular speed decreases, during a simulation, because of friction effects (see Figure5).Figure 5
Scheme of the rotor motion.


## 2.3. Auxiliary Bearing Model

The catcher bearing model describes the dynamics of the bearing components and allows the calculation of the interaction forces and moments between their contact surfaces. By solving the equations of motion of each element of the auxiliary bearing system, it is possible to calculate the position of the center of gravity of each element of the catcher bearing.The catcher bearing system (diameter of the bearing equal to 60 mm) is shown in Figure6 and is composed of the following elements:(1)
housing;(2)
ribbon damper;(3)
axial support;(4)
inner race;(5)
outer races;(6)
rolling elements;(7)
landing sleeve.Figure 6
Schematic representation of the auxiliary system.From a dynamical point of view, a simplifying assumption has been made: the catcher bearing elements are able to move in the plane, so the bearing inner race has 3 degrees of freedom (i.e., it can translate along they and z directions and rotate around the x axis) while the outer race can only translate in the plane without rotating.Equations (7) describe the translational dynamics of the rolling bearing inner race; the bearing inner race degrees of freedom are two because the inner race motion is considered to move in the radial plane of the bearing OYZ (Figure 5). This simplifying hypothesis allows obtaining a good compromise between the accuracy computational efficiency:(7)Mix¨i+Cix˙i=P+F,F=Fio+Fir+Frub,where xi contains the coordinates of the inner race center of gravity, Mi is the inner race mass matrix, Ci is the inner race damping matrix, P is the weight force of the inner race, and the vector F is the sum of the radial forces acting on the bearing inner race (see also Figure 7): Fio is the interaction force on the inner race exerted by the rolling elements, Fir is the contact force exchanged between the inner race and the rotor during the contact (including a normal FrN component and a tangential one Ffr; see (15)), and Frub is the resultant force generated during the axial rub between the thrust surface of the bearing inner race and the lateral surface of the rotor sleeve (see (19)).Figure 7
Schematic representation of the interactions between the elements of the auxiliary system.Equation (8) describes the rotational dynamics of the bearing inner race:(8)IiΩ˙i=Tfr+Tcb+Trubwhere Ii is the torsional moment of inertia of the bearing inner race and Tfr is the moment acting on the inner race and due to the action of the radial friction force Ffr; Tcb is the ball bearing drag torque due to the rolling friction of the rolling elements during the motion and Trub is the resultant torque due to the axial rub between the thrust surface of the bearing inner race and the lateral surface of the rotor sleeve.The outer surface of the outer race interacts with the housing through the ribbon damper, which is an elastic-damping element; furthermore, there is another elastic element: the axial support on the outer race, placed between the external surface of the outer race and the statoric flange.Equation (9) describes the translational dynamics of the bearing outer race, which can move into the radial plane OYZ (Figure 5):(9)Mox¨o+Cox˙o=Foi+Foh,where xo contains the coordinates of the outer race center of gravity,Mo and Co are the bearing outer race mass and damping matrices, the vector Foi is the interaction force on the outer race exerted by the rolling elements, and Foh is the radial force exchanged between the outer race and the housing flange.The auxiliary bearings have spherical rolling elements, which have been modelled as force elements, without representing their dynamic behavior. The elastic force exerted by the rolling elements allows representing the interactions between the outer and the inner race of the bearings.


## 2.4. Contact Models

The catcher bearing model simulates the plane motion of the auxiliary bearing elements and calculates the possible interaction forces and moments among them; the kinematic physical quantities are determined by solving the equations of motion of each elements (see Section2.3), while the interaction forces/moments are calculated through appropriate contact models based on the Hertz theory.In this section, the contact models implemented in the proposed bearing model will be explained, showing the formulation used to calculate the interaction forces and moments between the contact surfaces.The surfaces which interact during the motion are (see Figure7):(1)
radial contact between the rotor landing sleeve and the inner surface of the bearing inner race: the contact forceFir exchanged between the inner race and the rotor during the contact includes a normal FrN component and a tangential one Ffr; furthermore, a friction moment Tfr is generated;(2)
axial contact between the thrust surface of the rolling bearing inner race and the lateral side of the landing rotor sleeve: during axial or combined axial-radial rotor landings, axial rubbing force and moment are generatedFrub, Trub;(3)
elastic interaction between the bearing outer race and the bearing inner race through the rolling elements: the balls exert an elastic forceFio on the inner race and Foi on the outer races;(4)
the interaction force between the housing and the outer race is exerted through the ribbon damper elementFoh.


### 2.4.1. Rotor-Inner Race Radial Interaction

Concerning the development of this part of the model, see also the paper [19]. During a radial landing, the inner surface of the bearing inner ring and the rotor landing sleeve get in contact; the rotor bounces, after the impact, on the inner ring and continues to rotate supported by the auxiliary bearings. As a result of the interaction with the rotor, the inner ring accelerates and tends to reach the rotor speed.The contact forcesFrN exchanged between the rotor and the bearing inner ring is made up of two components: FrN is the normal contact force and Ffr is the friction force due to the radial rubbing.The mathematical steps required to calculate such forces are briefly reported below.It is very useful to define a radial-tangential (R(i1),T(j1)) reference system and express the quantities of interest (e.g., bodies centers of gravity coordinates; velocities) in this specific system (Figure 8).Figure 8
Interaction normal and friction forces between the rotor and the bearing inner ring.Furthermore, it is convenient to refer the analysis to the relative displacements and velocities between the rotor and the inner race in radial and tangential coordinates. The relative displacementsdri between the two bodies can be obtained from the difference between the coordinates of the rotor node Or and the coordinates of the inner race center of gravity Oi:(10)dri=Or-Oi.Through the use of the relative coordinates, it is possible to define a relative angle of rotationσrel and, consequently, to determine the rotational matrix R needed to convert the quantities expressed in the Cartesian reference system (Z(i0),Y(j0)) into the radial-tangential reference system (R(i1),T(j1)).Equation (11) defines the relative angle of rotation σrel:(11)σrel=atan⁡2yri,zri,while (12) defines the transformation matrix R:(12)R=cos⁡σrelsin⁡σrelsin⁡σrel-cos⁡σrel.The relative displacements dri and velocities vri expressed in the Cartesian reference system can be expressed into the radial-tangential reference system using the rotational matrix R:(13)drt=RTdri,vrt=RTvri.The contact forces and moments are calculated through the contact model if there is a interpenetrationδri between the two bodies. By knowing the physical value of the radial clearance, the penetration between the two bodies along the radial direction can be calculated as follows:(14)δri=yri-cr.Using the interpenetrationδri, it is possible to calculate the normal contact force FrN and the friction force Ffr due to the radial rubbing between the rotor sleeve radial surface and the inner surface of the inner ring:(15)FrN=Khertzδriγ+Kviscδ˙riδrii1,Ffr=μradFrNSj1where Khertz is the Hertzian contact constant, Kvisc is the viscous contact constant, μrad is the radial friction coefficient, and S is the sign function of the tangential component of the relative velocity vt:(16)vt=vrtj1+rrΩrk1-riΩik1.The friction force Ffr due to the radial rubbing between the rotor sleeve and the inner ring generates a torque Tfr responsible for the acceleration of the inner race:(17)Tfr=ri+rr2Ffrk1.During a pure radial landing, the rotor falls and bounces on to the inner ring of the bearing, which begins to rotate and quickly reaches the rotor angular speed; in case of high radial friction, the rotor performs a backward whirl orbit because of the friction effects between the rotor and the bearing inner ring.On the other hand, during a combined axial-radial rotor landing due to AMBs failure (both radial and axial magnetic bearings), the rotor gets in contact with the inner surface of the inner ring and the thrust face of the catcher bearing. In case of high friction effects, a circular forward whirl rotor orbit has been observed in some experimental tests; in fact, as a consequence of repeated rotor landings, the bearing elements begin to be worn and the friction coefficient increases; the tangential actions due to that axial friction result to be stronger than the radial rubbing force.


### 2.4.2. Rotor Sleeve-Inner Race Axial Interaction

As regards the development of this part of the model, see also the paper [19]. Concerning the lateral dynamics, the rotor can perform small rotations around Y and Z axes; as a consequence of the rotor dynamics, the interpenetration between the lateral face of the rotor sleeve and the thrust surface of the inner ring cannot have a constant value; by knowing the inclination angles α and β of the rotor at any instant of time, it is possible to calculate the interpenetration value δri at each point of the thrust surface and, consequently, calculate the axial pressure Paxial at any point on the contact surface.In order to calculate the resultant force and the resultant torque acting on the rotor due to the friction at the thrust surface, it is necessary to evaluate the relative slip velocity between the same pointA(r,θ) on the rotor and on the bearing.The vectorv˙A,r represents the velocity of the point A(r,θ) on the rotor, while the vector v˙A,i represents the velocity of the point A(r,θ) on the bearing inner ring:(18)vA,rrA,θA=r˙r-ω˙rrAsin⁡θAi1+t˙r-ω˙rrri-rAcos⁡θAj1,vA,irA,θA=r˙i-ω˙irAsin⁡θAi1+t˙i-ω˙irAcos⁡θAj1.The relative velocity vA,ri between the point A on the rotor and the same point A on the bearing can be simply calculated by subtracting the inner ring velocity vA,r from the rotor velocity vA,i.By knowing the relative velocity between the rotor and the bearing, the resultant forceFrub, acting on the rotor and due to friction at the contact interface, can be evaluated by integrating the pressure field over the contact area:(19)Frub,r=-μax∫r1r2∫02πPA,rr˙ri1+t˙rj1vA,rirAdθAdRA,where μax is the friction coefficient between the lateral face of the rotor sleeve and the thrust face of the bearing.The torque due to the axial contact between the rotor and the bearing inner race can be determined as follows:(20)Trub,r=-μax∫r1r2∫02πPA,rt˙rcos⁡θA+r˙rsin⁡θAvA,rirA2dθAdRAk1.The torque on the rotor and on the bearing is considered to be equal and opposite:(21)Trub,r=-Trub,i.


### 2.4.3. Outer Race-Inner Race Radial Interaction

The interaction forceFio between the bearing inner and outer races is the radial reaction due to the relative motion of the races.The races can interact through the rolling elements, which are modelled as force elements; through the load-deflection curve of the catcher bearing it is possible to extrapolate the stiffness characteristics of the rolling elements (kb1 and kb2) and define the interaction force Fio with a second order polynomial function:(22)Fio=kb2rio2+kb1rioi1,where rio is the relative displacement between the inner and the outer races centers of gravity.The rolling elements exert also another important action on the catcher bearing: the rolling friction due to the balls generates a drag torque onto the inner ring, which decelerates under this effect. The approach to calculate the resistant torqueTCB is based on the Harris-Palmgren formulation [23], but, instead of considering the applied static load, the dynamic radial and axial forces acting on the bearing have been taken into account:(23)TCB=f1Fβdmk1,where f1 is a factor depending on the bearing design and relative bearing load, Fβ depends on the magnitude and direction of the applied load, and dm is the pitch diameter of the bearing.For a ball bearing, Harris [23] defines f1 as(24)f1=zFsyCs,where Fs represents the equivalent dynamic load and Cs is the basic static load rating usually provided in the manufacturer catalog, while z and y represent empirical coefficients.For a single row angular contact bearing the equivalent loadFs can be defined as follows:(25)Fs=XsFradial+YsFax,where Xs=0.5 and Ys varies with the bearing contact angle α:(26)Ys=c1α2+c2α+c3.For a radial ball bearing Fβ is given by (27)Fβ=c4Faxcot⁡α+c5Fradial.The numeric parametersc1, c2, c3, c4, and c5 of (26) and (27) have been determined on the base of the bearing geometry and characteristics: c1=7.142910-5, c2=-1.244310-2, c3=0.64129, c4=0.9, and c5=-0.1.


### 2.4.4. Outer Race-Housing Interaction

The ribbon damper is located between the bearing outer ring and the housing; it has both elastic and damping characteristics (i.e., structural and viscous damping) and generates radial and tangential force components. The characteristics of the ribbon greatly influence the frequency of the rotor whirl motion.The interaction force between the outer race and the housingFoh can be expressed through the force exerted by the ribbon:(28)Foh=-Cribx˙o-KribQxo-KribIxo,where I is the identity matrix, xo is the vector of the coordinates of the outer race center of gravity, Crib is the ribbon viscous damping, Krib=6.3∗107N/m is the ribbon stiffness, and Q is a matrix defined as follows:(29)Q=0-ηSvηSv0.The parameter η is a cross coupled stiffness factor, while Sv is the sign of the relative tangential velocity between the outer bearing and the housing.For combined axial-radial rotor landings, because of the axial thrust force applied to the bearing inner race, a reaction force exerted by the housing on the outer race through the axial support is present. The force contribution due to this interaction is the damping forceFoh:(30)Foh=-Faxμdx˙ox˙o,where μd is the Coulomb friction coefficient between the outer race and the housing and must be added to (28).


## 2.4.1. Rotor-Inner Race Radial Interaction

Concerning the development of this part of the model, see also the paper [19]. During a radial landing, the inner surface of the bearing inner ring and the rotor landing sleeve get in contact; the rotor bounces, after the impact, on the inner ring and continues to rotate supported by the auxiliary bearings. As a result of the interaction with the rotor, the inner ring accelerates and tends to reach the rotor speed.The contact forcesFrN exchanged between the rotor and the bearing inner ring is made up of two components: FrN is the normal contact force and Ffr is the friction force due to the radial rubbing.The mathematical steps required to calculate such forces are briefly reported below.It is very useful to define a radial-tangential (R(i1),T(j1)) reference system and express the quantities of interest (e.g., bodies centers of gravity coordinates; velocities) in this specific system (Figure 8).Figure 8
Interaction normal and friction forces between the rotor and the bearing inner ring.Furthermore, it is convenient to refer the analysis to the relative displacements and velocities between the rotor and the inner race in radial and tangential coordinates. The relative displacementsdri between the two bodies can be obtained from the difference between the coordinates of the rotor node Or and the coordinates of the inner race center of gravity Oi:(10)dri=Or-Oi.Through the use of the relative coordinates, it is possible to define a relative angle of rotationσrel and, consequently, to determine the rotational matrix R needed to convert the quantities expressed in the Cartesian reference system (Z(i0),Y(j0)) into the radial-tangential reference system (R(i1),T(j1)).Equation (11) defines the relative angle of rotation σrel:(11)σrel=atan⁡2yri,zri,while (12) defines the transformation matrix R:(12)R=cos⁡σrelsin⁡σrelsin⁡σrel-cos⁡σrel.The relative displacements dri and velocities vri expressed in the Cartesian reference system can be expressed into the radial-tangential reference system using the rotational matrix R:(13)drt=RTdri,vrt=RTvri.The contact forces and moments are calculated through the contact model if there is a interpenetrationδri between the two bodies. By knowing the physical value of the radial clearance, the penetration between the two bodies along the radial direction can be calculated as follows:(14)δri=yri-cr.Using the interpenetrationδri, it is possible to calculate the normal contact force FrN and the friction force Ffr due to the radial rubbing between the rotor sleeve radial surface and the inner surface of the inner ring:(15)FrN=Khertzδriγ+Kviscδ˙riδrii1,Ffr=μradFrNSj1where Khertz is the Hertzian contact constant, Kvisc is the viscous contact constant, μrad is the radial friction coefficient, and S is the sign function of the tangential component of the relative velocity vt:(16)vt=vrtj1+rrΩrk1-riΩik1.The friction force Ffr due to the radial rubbing between the rotor sleeve and the inner ring generates a torque Tfr responsible for the acceleration of the inner race:(17)Tfr=ri+rr2Ffrk1.During a pure radial landing, the rotor falls and bounces on to the inner ring of the bearing, which begins to rotate and quickly reaches the rotor angular speed; in case of high radial friction, the rotor performs a backward whirl orbit because of the friction effects between the rotor and the bearing inner ring.On the other hand, during a combined axial-radial rotor landing due to AMBs failure (both radial and axial magnetic bearings), the rotor gets in contact with the inner surface of the inner ring and the thrust face of the catcher bearing. In case of high friction effects, a circular forward whirl rotor orbit has been observed in some experimental tests; in fact, as a consequence of repeated rotor landings, the bearing elements begin to be worn and the friction coefficient increases; the tangential actions due to that axial friction result to be stronger than the radial rubbing force.


## 2.4.2. Rotor Sleeve-Inner Race Axial Interaction

As regards the development of this part of the model, see also the paper [19]. Concerning the lateral dynamics, the rotor can perform small rotations around Y and Z axes; as a consequence of the rotor dynamics, the interpenetration between the lateral face of the rotor sleeve and the thrust surface of the inner ring cannot have a constant value; by knowing the inclination angles α and β of the rotor at any instant of time, it is possible to calculate the interpenetration value δri at each point of the thrust surface and, consequently, calculate the axial pressure Paxial at any point on the contact surface.In order to calculate the resultant force and the resultant torque acting on the rotor due to the friction at the thrust surface, it is necessary to evaluate the relative slip velocity between the same pointA(r,θ) on the rotor and on the bearing.The vectorv˙A,r represents the velocity of the point A(r,θ) on the rotor, while the vector v˙A,i represents the velocity of the point A(r,θ) on the bearing inner ring:(18)vA,rrA,θA=r˙r-ω˙rrAsin⁡θAi1+t˙r-ω˙rrri-rAcos⁡θAj1,vA,irA,θA=r˙i-ω˙irAsin⁡θAi1+t˙i-ω˙irAcos⁡θAj1.The relative velocity vA,ri between the point A on the rotor and the same point A on the bearing can be simply calculated by subtracting the inner ring velocity vA,r from the rotor velocity vA,i.By knowing the relative velocity between the rotor and the bearing, the resultant forceFrub, acting on the rotor and due to friction at the contact interface, can be evaluated by integrating the pressure field over the contact area:(19)Frub,r=-μax∫r1r2∫02πPA,rr˙ri1+t˙rj1vA,rirAdθAdRA,where μax is the friction coefficient between the lateral face of the rotor sleeve and the thrust face of the bearing.The torque due to the axial contact between the rotor and the bearing inner race can be determined as follows:(20)Trub,r=-μax∫r1r2∫02πPA,rt˙rcos⁡θA+r˙rsin⁡θAvA,rirA2dθAdRAk1.The torque on the rotor and on the bearing is considered to be equal and opposite:(21)Trub,r=-Trub,i.


## 2.4.3. Outer Race-Inner Race Radial Interaction

The interaction forceFio between the bearing inner and outer races is the radial reaction due to the relative motion of the races.The races can interact through the rolling elements, which are modelled as force elements; through the load-deflection curve of the catcher bearing it is possible to extrapolate the stiffness characteristics of the rolling elements (kb1 and kb2) and define the interaction force Fio with a second order polynomial function:(22)Fio=kb2rio2+kb1rioi1,where rio is the relative displacement between the inner and the outer races centers of gravity.The rolling elements exert also another important action on the catcher bearing: the rolling friction due to the balls generates a drag torque onto the inner ring, which decelerates under this effect. The approach to calculate the resistant torqueTCB is based on the Harris-Palmgren formulation [23], but, instead of considering the applied static load, the dynamic radial and axial forces acting on the bearing have been taken into account:(23)TCB=f1Fβdmk1,where f1 is a factor depending on the bearing design and relative bearing load, Fβ depends on the magnitude and direction of the applied load, and dm is the pitch diameter of the bearing.For a ball bearing, Harris [23] defines f1 as(24)f1=zFsyCs,where Fs represents the equivalent dynamic load and Cs is the basic static load rating usually provided in the manufacturer catalog, while z and y represent empirical coefficients.For a single row angular contact bearing the equivalent loadFs can be defined as follows:(25)Fs=XsFradial+YsFax,where Xs=0.5 and Ys varies with the bearing contact angle α:(26)Ys=c1α2+c2α+c3.For a radial ball bearing Fβ is given by (27)Fβ=c4Faxcot⁡α+c5Fradial.The numeric parametersc1, c2, c3, c4, and c5 of (26) and (27) have been determined on the base of the bearing geometry and characteristics: c1=7.142910-5, c2=-1.244310-2, c3=0.64129, c4=0.9, and c5=-0.1.


## 2.4.4. Outer Race-Housing Interaction

The ribbon damper is located between the bearing outer ring and the housing; it has both elastic and damping characteristics (i.e., structural and viscous damping) and generates radial and tangential force components. The characteristics of the ribbon greatly influence the frequency of the rotor whirl motion.The interaction force between the outer race and the housingFoh can be expressed through the force exerted by the ribbon:(28)Foh=-Cribx˙o-KribQxo-KribIxo,where I is the identity matrix, xo is the vector of the coordinates of the outer race center of gravity, Crib is the ribbon viscous damping, Krib=6.3∗107N/m is the ribbon stiffness, and Q is a matrix defined as follows:(29)Q=0-ηSvηSv0.The parameter η is a cross coupled stiffness factor, while Sv is the sign of the relative tangential velocity between the outer bearing and the housing.For combined axial-radial rotor landings, because of the axial thrust force applied to the bearing inner race, a reaction force exerted by the housing on the outer race through the axial support is present. The force contribution due to this interaction is the damping forceFoh:(30)Foh=-Faxμdx˙ox˙o,where μd is the Coulomb friction coefficient between the outer race and the housing and must be added to (28).


## 2.5. Test Rig Architecture

Figure9 shows an overview of the scaled test rig considered for the development and the validation of the proposed model. The main components of the test rig are highlighted in the following:(1)
electric motor(1);(2)
power transmission coupling(2);(3)
force sensor(3);(4)
radial AMBs and auxiliary bearing(4); they are concentric;(5)
axial centering AMB(5);(6)
rotor(6);(7)
axial actuator(7).Figure 9
Experimental apparatus; from left to right: axial actuator, NDE radial magnetic bearing, speed sensor, DE radial magnetic bearing, magnetic thrust bearing, and power transmission coupling.Figure10 shows a scheme of the test rig that is equipped with a sensors system able to estimate components temperature, forces, acceleration, vibrations, and speeds:(1)
load cells(LC);(2)
proximity sensors(PS);(3)
encoders(EN);(4)
accelerometers(A);(5)
thermocouples(TC).Figure 10
Scheme of the experimental apparatus.The tests performed on the test rig allow the identification of the system behavior in landing conditions. In particular the tests carried out on the test rig belong to three different types (with different rotor speeds): radial landings (deactivating just radial AMBs), axial landings (deactivating just axial AMB), and combined radial-axial landings (deactivating both the AMBs).The experimental tests taken as a reference in this paper for the model validation (Section3) are no whirling rotor radial landing, a pure radial landing characterized by a normal friction condition; whirling rotor radial landing, a pure radial landing with high friction coefficient; whirling rotor combining axial-radial landing with high axial thrust force value.


## 3. Results

The proposed model allows performing the simulation for all the classical operating conditions and to set up parameters like the rotor angular velocity, the friction coefficients, and the landing time.In this section the most interesting results of the main simulations are reported and compared to the related experimental results. In particular, it is easy to observe the different dynamic behavior of the rotor in case of radial landing (with normal friction conditions and high friction conditions) and of combined axial-radial landing. During a pure radial landing characterized by normal friction conditions, the rotor bounces on to the bearing inner rings and, after the first bounces, the rotor and the inner ring rotate at the same angular speed. In case of a pure radial landing characterized by high radial friction coefficient, the rotor realizes a backward orbit; while in a combined axial-radial landing with a high axial thrust force value acting on the rotor, the tangential forces generated in the axial contact cause a rotor forward whirl. The backward and forward whirl phenomena occur in particular conditions of high axial or radial friction; in fact, after a sequence of radial or combined axial-radial landings, some bearing elements begin to be worn (e.g., bearing inner race) as a consequence of the high impact and friction forces and consequently, the friction coefficient increases causing the whirling phenomenon.All the simulation results reported in the paper have been obtained by using both variable step ODE algorithms and fixed step ODE algorithms [24]. The choice of the constant step for the fixed step algorithms and of the maximum step size for the variable step algorithms (together with the related tolerances) has been made to have a characteristic integration frequency (i.e., the inverse of the time step) always greater than the most significant frequencies of the considered system.


### 3.1. No Whirling Landing

In case of overloads or failure of the energy supply, the AMBs are disabled and the rotor falls and impacts on the inner rings of the auxiliary bearings.The reported results show the displacements, forces, and orbits of specific rotor nodes; in particular, the nodes where the auxiliary bearings are located at the drive end and at the nondrive end of the machine (for brevity denoted DE node and NDE node) have been considered.The reported simulation and experimental test reproduces a rotor landing, in which the rotor has an angular speed equal to 7000 RPM and the AMBS are deactivated; so, as a consequence, the rotor lands on the auxiliary bearing system bouncing on to the bearing inner rings. During this first radial landing test (low radial friction without backward whirling), the following values were set for the speedΩ, the radial friction μrad, the axial friction μax, and the axial load Fthrust: Ω=7000rpm, μrad=0.15 (low radial friction), μax=0.3, and Fthrust=0N (no axial load).In Figures11 and 14 it is possible to observe the vertical and the horizontal displacements of the DE rotor node; by comparing the simulated plots with the experimental displacements reported in Figures 13 and 16, it is possible to note a good representation of dynamic rotor behavior. The vertical displacements show some initial peak values, which correspond to the first rotor bounces on the inner rings characterized by a high energy content; after that, the rotor and the bearing inner ring rotates at the same angular speed. In Figures 12 and 15 the vertical and the horizontal forces acting on DE node show peak values in correspondence of the first rotor bounces.Figure 11
Predicted rotor vertical displacements during a landing test (DE node).Figure 12
Predicted rotor vertical forces during a landing test (DE node).Figure 13
Experimental rotor vertical displacements during a landing test at DE (blue line) and NDE (red line) rotor nodes.Figure 14
Predicted rotor horizontal displacements during a landing test (DE node).Figure 15
Predicted rotor horizontal forces during a landing test (DE node).Figure 16
Experimental rotor horizontal displacements during a landing test at DE (blue line) and NDE (red line) rotor nodes.Figure17 shows the orbit described by the DE rotor node during the simulation; the blue circle is the nominal gap between the landing bearing bore and the landing sleeve, while in Figure 18 the experimental rotor orbit is reported.Figure 17
Predicted rotor orbit during a landing test (DE node).Figure 18
Experimental rotor orbit during a landing test.


### 3.2. Whirling Radial Landing

During a radial landing the rotor falls and impacts on the inner ring of the auxiliary bearing, it bounces on the ring, and in case of high friction and any conditions of wear, after some second tenths the rotor begins to whirl; in case of radial landing the contact force between the rotor sleeve and the inner ring of the bearing causes a backward whirl rotor orbit (Figure27), considering the rotor angular speed convention shown in Figure 19. The reported results show the displacements, forces, and orbits of specific rotor nodes; in particular, the nodes where the auxiliary bearings are located at the drive end and at the nondrive end of the machine (for brevity denoted DE node and NDE node) have been considered. The rotor backward whirl is the predominant phenomenon in case of pure radial landings because of the radial friction effects.Figure 19
Rotor angular speed convention: rotor positive angular speed is clockwise.The reference simulation reproduces a pure radial landing where the rotor rotates with an angular speed equal to 7000 RPM, the AMBS are deactivated, and the rotor lands on the auxiliary bearing system; between the inner race of the catcher bearing and the rotor landing sleeve surface, friction radial forces arise: such forces are responsible for the rotor backward whirl orbit.In order to emphasize such behavior, the radial friction coefficient has been set to a quite high value with respect to the combined axial-radial landings, where the axial friction coefficient is predominant on the radial one. The high value of the friction coefficient simulates the wear condition of the bearing elements due to the repeated landings.During this second radial landing test (high radial friction with backward whirling), the following values were set for the speedΩ, the radial friction μrad, the axial friction μax, and the axial load Fthrust: Ω=7000rpm, μrad=0.3 (high radial friction), μax=0.3, and Fthrust=0N (no axial load).In Figures20 and 24 it is possible to observe the vertical and the horizontal displacements of the DE rotor node (the time zooms of such quantities are reported in Figure 21 to better highlight amplitudes, phases and frequencies). By comparing the simulated plots with the experimental displacements reported in Figures 23 and 26, it is possible to note a good representation of the backward rotor orbit. The vertical displacements show some initial peak values, which represent the rotor bounces on the inner ring; after 0,5 s the vertical displacements switch from maximum negative values to maximum positive values: such oscillatory trend produces the backward whirl rotor orbit. Figure 27 show the orbits described by the DE rotor node during the simulation; the blue circle represents the nominal gap between the landing bearing bore and the landing sleeve. In Figures 22 and 25 the vertical and the horizontal forces on the rotor node in correspondence of the auxiliary bearing located at the DE side of the machine are reported; it is possible to observe the high dynamical values of such forces during the backward whirl.Figure 20
Predicted rotor vertical displacements during a radial landing test (DE node).Figure 21
Time zoom of the predicted rotor vertical (in blue) and horizontal (in red) displacements during a radial landing test (DE node).Figure 22
Predicted rotor vertical forces during a radial landing test (DE node).Figure 23
Rotor vertical displacements during a radial landing experimental test at DE (blu line) and NDE (red line) rotor nodes.Figure 24
Predicted rotor horizontal displacements during a radial landing test (DE node).Figure 25
Predicted rotor horizontal forces during a radial landing test (DE node).Figure 26
Rotor horizontal displacements during a radial landing experimental test at DE (blue line) and NDE (red line) rotor nodes.Figure 27
Predicted rotor backward whirl orbit at the DE rotor nodes.In Table1 the comparison between some significant results of the experimental tests and the related numerical simulations is reported: there is a good alignment between the displacements, and the activation time of the backward whirl has the same values.Table 1
Radial landing results comparison.
Simulated Experimental V.Disp.[mm] Act.Time[s] V. Disp.[mm] Act.Time[s]  0.25 -0.28 0.4 0.2 -0.4 0.4  Backw. Freq.   Backw. Freq. [Hz]   [Hz]  51   49


### 3.3. Whirling Combined Axial-Radial Landing

In this section the most significant results of the axial-radial combined landing have been reported. The results reported in this paper are aimed at highlighting the forward whirl, the most important phenomenon observed during the combined landings; it occurs when the tangential forces between the lateral surface of the inner ring and the thrust surface of the rotor sleeve, caused by the axial friction, have a high value. In this condition such tangential forces prevail on the radial friction effects, which are responsible for the rotor backward whirl, and the forward whirl may occur. As in the previous case, the forward whirl is a particular phenomenon which may take place in case of high friction and wear conditions.The considered simulation corresponds to a test in which the rotor rotates with an angular speed equal to 7000 RPM supported by the AMBs; the AMBs are deactivated, the rotor falls, and lands on to the auxiliary bearings. At the same time also the axial AMB is deactivated and an axial thrust force acts on the system determining the axial displacement of the rotor and the contact between the inner race and the landing sleeve.In case of a combined radial-axial landing the rotor performs a forward whirl orbit: in fact, during the axial landing the thrust surface of the inner ring gets in contact with the lateral surface of the rotor sleeve; because of this interaction, friction tangential forces are generated and the rotor begins to perform a forward whirl orbit. During the axial-radial landing test (with forward whirling), the following values were set for the speedΩ, the radial friction μrad, the axial friction μax, and the axial load Fthrust: Ω=7000rpm, μrad=0.15 (low radial friction), μax=0.3, and Fthrust=4000N.In Figures28, 30, 31, and 32 it is possible to observe vertical and horizontal displacements and the corresponding force trends during the whole simulation. The time zooms of vertical and horizontal displacements of the DE rotor node are reported in Figure 29 to better highlight amplitudes, phases, and frequencies. Figures 33 and 34 show the results obtained from the experimental test corresponding to the simulation of the combined landing.Figure 28
Predicted rotor vertical displacements at DE node during a combined landing test.Figure 29
Time zoom of the predicted rotor vertical (in blue) and horizontal (in red) displacements at DE node during a combined landing test.Figure 30
Predicted rotor vertical forces at DE node during a combined landing test.Figure 31
Predicted rotor horizontal displacements at DE node during a combined landing test.Figure 32
Predicted rotor horizontal forces at DE node during a combined landing test.Figure 33
Experimental vertical displacements and forces at DE (blue line) and NDE (red line) rotor nodes.Figure 34
Experimental horizontal displacements and forces at DE (blu line) and NDE (red line) rotor nodes.By making a comparison between the resultant vertical displacements and forces of the simulation and the experimental results, it is possible to find a good agreement between experimental and numerical physical quantities.The orbits of the DE and NDE rotor nodes resulting from the numerical simulation are a bit different from the experimental orbits: in the experimental test the rotor begins to whirl slower and the orbits are fuller as a consequence of the whirling dynamics observed on the test rig: the whirling activation and deactivation on the real test are slower compared to the results obtained through the numerical simulation (see Figures35, 36, and 37).Figure 35
Predicted rotor DE node orbit.Figure 36
Predicted rotor NDE node orbit.Figure 37
Experimental rotor orbit of the DE node (blue line) and NDE node (red line) rotor nodes.In order to extrapolate the frequency content of the vertical displacement of the DE rotor node signal, the Fourier Transform Method has been applied; through this method the signal has been described in the frequency domain to identify the characteristic forward whirl frequency (Figure38).Figure 38
Simulated whirling frequency (rotor angular speed 7000 RPM).In Table2, the comparison between some numerical results of the simulation and the experimental data obtained from the experimental tests is reported for a combined axial-radial landing: there is a good alignment between the dynamic forces and the related vertical displacements. Furthermore, the forward whirl frequencies, reported for three different speeds Ω=3000rpm - 7000 rpm - 10000 rpm, are well reproduced (see Figure 39).Table 2
Combined landing results comparison.
Simulated Experimental V.Disp. V.Forces V. Disp. V.Forces [mm] [KN] [mm] [KN]  0.2 -0.28 4.5 -1.5 0.15 -0.3 4.5 -0.5  Forw. Freq.   Forw. Freq. [Hz]   [Hz]  50   48 67   65 80   77Figure 39
Comparison between experimental and simulated forward whirl frequency.


## 3.1. No Whirling Landing

In case of overloads or failure of the energy supply, the AMBs are disabled and the rotor falls and impacts on the inner rings of the auxiliary bearings.The reported results show the displacements, forces, and orbits of specific rotor nodes; in particular, the nodes where the auxiliary bearings are located at the drive end and at the nondrive end of the machine (for brevity denoted DE node and NDE node) have been considered.The reported simulation and experimental test reproduces a rotor landing, in which the rotor has an angular speed equal to 7000 RPM and the AMBS are deactivated; so, as a consequence, the rotor lands on the auxiliary bearing system bouncing on to the bearing inner rings. During this first radial landing test (low radial friction without backward whirling), the following values were set for the speedΩ, the radial friction μrad, the axial friction μax, and the axial load Fthrust: Ω=7000rpm, μrad=0.15 (low radial friction), μax=0.3, and Fthrust=0N (no axial load).In Figures11 and 14 it is possible to observe the vertical and the horizontal displacements of the DE rotor node; by comparing the simulated plots with the experimental displacements reported in Figures 13 and 16, it is possible to note a good representation of dynamic rotor behavior. The vertical displacements show some initial peak values, which correspond to the first rotor bounces on the inner rings characterized by a high energy content; after that, the rotor and the bearing inner ring rotates at the same angular speed. In Figures 12 and 15 the vertical and the horizontal forces acting on DE node show peak values in correspondence of the first rotor bounces.Figure 11
Predicted rotor vertical displacements during a landing test (DE node).Figure 12
Predicted rotor vertical forces during a landing test (DE node).Figure 13
Experimental rotor vertical displacements during a landing test at DE (blue line) and NDE (red line) rotor nodes.Figure 14
Predicted rotor horizontal displacements during a landing test (DE node).Figure 15
Predicted rotor horizontal forces during a landing test (DE node).Figure 16
Experimental rotor horizontal displacements during a landing test at DE (blue line) and NDE (red line) rotor nodes.Figure17 shows the orbit described by the DE rotor node during the simulation; the blue circle is the nominal gap between the landing bearing bore and the landing sleeve, while in Figure 18 the experimental rotor orbit is reported.Figure 17
Predicted rotor orbit during a landing test (DE node).Figure 18
Experimental rotor orbit during a landing test.


## 3.2. Whirling Radial Landing

During a radial landing the rotor falls and impacts on the inner ring of the auxiliary bearing, it bounces on the ring, and in case of high friction and any conditions of wear, after some second tenths the rotor begins to whirl; in case of radial landing the contact force between the rotor sleeve and the inner ring of the bearing causes a backward whirl rotor orbit (Figure27), considering the rotor angular speed convention shown in Figure 19. The reported results show the displacements, forces, and orbits of specific rotor nodes; in particular, the nodes where the auxiliary bearings are located at the drive end and at the nondrive end of the machine (for brevity denoted DE node and NDE node) have been considered. The rotor backward whirl is the predominant phenomenon in case of pure radial landings because of the radial friction effects.Figure 19
Rotor angular speed convention: rotor positive angular speed is clockwise.The reference simulation reproduces a pure radial landing where the rotor rotates with an angular speed equal to 7000 RPM, the AMBS are deactivated, and the rotor lands on the auxiliary bearing system; between the inner race of the catcher bearing and the rotor landing sleeve surface, friction radial forces arise: such forces are responsible for the rotor backward whirl orbit.In order to emphasize such behavior, the radial friction coefficient has been set to a quite high value with respect to the combined axial-radial landings, where the axial friction coefficient is predominant on the radial one. The high value of the friction coefficient simulates the wear condition of the bearing elements due to the repeated landings.During this second radial landing test (high radial friction with backward whirling), the following values were set for the speedΩ, the radial friction μrad, the axial friction μax, and the axial load Fthrust: Ω=7000rpm, μrad=0.3 (high radial friction), μax=0.3, and Fthrust=0N (no axial load).In Figures20 and 24 it is possible to observe the vertical and the horizontal displacements of the DE rotor node (the time zooms of such quantities are reported in Figure 21 to better highlight amplitudes, phases and frequencies). By comparing the simulated plots with the experimental displacements reported in Figures 23 and 26, it is possible to note a good representation of the backward rotor orbit. The vertical displacements show some initial peak values, which represent the rotor bounces on the inner ring; after 0,5 s the vertical displacements switch from maximum negative values to maximum positive values: such oscillatory trend produces the backward whirl rotor orbit. Figure 27 show the orbits described by the DE rotor node during the simulation; the blue circle represents the nominal gap between the landing bearing bore and the landing sleeve. In Figures 22 and 25 the vertical and the horizontal forces on the rotor node in correspondence of the auxiliary bearing located at the DE side of the machine are reported; it is possible to observe the high dynamical values of such forces during the backward whirl.Figure 20
Predicted rotor vertical displacements during a radial landing test (DE node).Figure 21
Time zoom of the predicted rotor vertical (in blue) and horizontal (in red) displacements during a radial landing test (DE node).Figure 22
Predicted rotor vertical forces during a radial landing test (DE node).Figure 23
Rotor vertical displacements during a radial landing experimental test at DE (blu line) and NDE (red line) rotor nodes.Figure 24
Predicted rotor horizontal displacements during a radial landing test (DE node).Figure 25
Predicted rotor horizontal forces during a radial landing test (DE node).Figure 26
Rotor horizontal displacements during a radial landing experimental test at DE (blue line) and NDE (red line) rotor nodes.Figure 27
Predicted rotor backward whirl orbit at the DE rotor nodes.In Table1 the comparison between some significant results of the experimental tests and the related numerical simulations is reported: there is a good alignment between the displacements, and the activation time of the backward whirl has the same values.Table 1
Radial landing results comparison.
Simulated Experimental V.Disp.[mm] Act.Time[s] V. Disp.[mm] Act.Time[s]  0.25 -0.28 0.4 0.2 -0.4 0.4  Backw. Freq.   Backw. Freq. [Hz]   [Hz]  51   49


## 3.3. Whirling Combined Axial-Radial Landing

In this section the most significant results of the axial-radial combined landing have been reported. The results reported in this paper are aimed at highlighting the forward whirl, the most important phenomenon observed during the combined landings; it occurs when the tangential forces between the lateral surface of the inner ring and the thrust surface of the rotor sleeve, caused by the axial friction, have a high value. In this condition such tangential forces prevail on the radial friction effects, which are responsible for the rotor backward whirl, and the forward whirl may occur. As in the previous case, the forward whirl is a particular phenomenon which may take place in case of high friction and wear conditions.The considered simulation corresponds to a test in which the rotor rotates with an angular speed equal to 7000 RPM supported by the AMBs; the AMBs are deactivated, the rotor falls, and lands on to the auxiliary bearings. At the same time also the axial AMB is deactivated and an axial thrust force acts on the system determining the axial displacement of the rotor and the contact between the inner race and the landing sleeve.In case of a combined radial-axial landing the rotor performs a forward whirl orbit: in fact, during the axial landing the thrust surface of the inner ring gets in contact with the lateral surface of the rotor sleeve; because of this interaction, friction tangential forces are generated and the rotor begins to perform a forward whirl orbit. During the axial-radial landing test (with forward whirling), the following values were set for the speedΩ, the radial friction μrad, the axial friction μax, and the axial load Fthrust: Ω=7000rpm, μrad=0.15 (low radial friction), μax=0.3, and Fthrust=4000N.In Figures28, 30, 31, and 32 it is possible to observe vertical and horizontal displacements and the corresponding force trends during the whole simulation. The time zooms of vertical and horizontal displacements of the DE rotor node are reported in Figure 29 to better highlight amplitudes, phases, and frequencies. Figures 33 and 34 show the results obtained from the experimental test corresponding to the simulation of the combined landing.Figure 28
Predicted rotor vertical displacements at DE node during a combined landing test.Figure 29
Time zoom of the predicted rotor vertical (in blue) and horizontal (in red) displacements at DE node during a combined landing test.Figure 30
Predicted rotor vertical forces at DE node during a combined landing test.Figure 31
Predicted rotor horizontal displacements at DE node during a combined landing test.Figure 32
Predicted rotor horizontal forces at DE node during a combined landing test.Figure 33
Experimental vertical displacements and forces at DE (blue line) and NDE (red line) rotor nodes.Figure 34
Experimental horizontal displacements and forces at DE (blu line) and NDE (red line) rotor nodes.By making a comparison between the resultant vertical displacements and forces of the simulation and the experimental results, it is possible to find a good agreement between experimental and numerical physical quantities.The orbits of the DE and NDE rotor nodes resulting from the numerical simulation are a bit different from the experimental orbits: in the experimental test the rotor begins to whirl slower and the orbits are fuller as a consequence of the whirling dynamics observed on the test rig: the whirling activation and deactivation on the real test are slower compared to the results obtained through the numerical simulation (see Figures35, 36, and 37).Figure 35
Predicted rotor DE node orbit.Figure 36
Predicted rotor NDE node orbit.Figure 37
Experimental rotor orbit of the DE node (blue line) and NDE node (red line) rotor nodes.In order to extrapolate the frequency content of the vertical displacement of the DE rotor node signal, the Fourier Transform Method has been applied; through this method the signal has been described in the frequency domain to identify the characteristic forward whirl frequency (Figure38).Figure 38
Simulated whirling frequency (rotor angular speed 7000 RPM).In Table2, the comparison between some numerical results of the simulation and the experimental data obtained from the experimental tests is reported for a combined axial-radial landing: there is a good alignment between the dynamic forces and the related vertical displacements. Furthermore, the forward whirl frequencies, reported for three different speeds Ω=3000rpm - 7000 rpm - 10000 rpm, are well reproduced (see Figure 39).Table 2
Combined landing results comparison.
Simulated Experimental V.Disp. V.Forces V. Disp. V.Forces [mm] [KN] [mm] [KN]  0.2 -0.28 4.5 -1.5 0.15 -0.3 4.5 -0.5  Forw. Freq.   Forw. Freq. [Hz]   [Hz]  50   48 67   65 80   77Figure 39
Comparison between experimental and simulated forward whirl frequency.


## 4. Conclusions

The model presented in this paper is able to predict the behavior of a complete auxiliary bearing system under transient, steady, and emergency conditions (e.g., delevitation and landing phenomena) and it is composed of different submodels: the rotor submodel, which includes a FEM model to describe the lateral dynamics of the rotor and two rigid body models to reproduce axial and torsional rotor dynamics, and the auxiliary rolling bearing submodel, characterized by a combination of multibody and contact models able to describe the dynamics and the interactions of the auxiliary system components.The model is able to reproduce the key effects observed during the experimental tests and, in particular, the forward whirl phenomenon. In order to carefully study some special phenomena like rotor landings, which require long period of time to evolve and involve many bodies and degrees of freedom, or other particular events like impacts, which occur in a short period of time, a compromise between accuracy and numerical efficiency has been pursued.To validate the model, experimental tests have been performed on a scaled test rig in collaboration with Baker Hughes a GE company. A radial and a combined axial-radial landing of the rotor on the auxiliary bearings under specific operating conditions have been taken as reference in order to validate the numerical model. The physical and the contact parameters have been correctly set thanks to the many experimental data.The results of the simulations well match the data coming from the experimental tests, confirming the reliability of the model and its robustness against the uncertainties affecting the system parameters.Furthermore, thanks to the model modularity, different rotor geometries and different auxiliary bearings sizes and locations can be easily taken into account.Concerning the future developments of the work, the authors will improve the description of the coupling between the axial and the lateral dynamics of the system. Up to now, the authors considered the coupling between axial and lateral dynamics by studying the effect of axial dynamics on the lateral one. In the future, also the effect of lateral dynamics on the axial one will be investigated, especially in terms of axial forces and related vibrations. This way, a complete understanding of the phenomenon will be obtained.Moreover, it could be interesting to develop a thermal and a wear model of the system. The thermal model could allow the monitoring of the temperature trend of the elements of the system, while the wear model could lead to an estimation of the wear conditions of the auxiliary bearing system after the emergency landings.Finally, detailed three-dimensional (3D) multibody and flexible multibody models of the system will be developed to reach a higher level of accuracy, especially concerning the contacts description.



---

*Source: 4675286-2019-04-01.xml*