Issue 
Mechanics & Industry
Volume 22, 2021



Article Number  35  
Number of page(s)  13  
DOI  https://doi.org/10.1051/meca/2021034  
Published online  28 May 2021 
Regular Article
Dynamics of onboard rotors on finitelength journal bearings subject to multiaxial and multifrequency excitations: numerical and experimental investigations
^{1}
Univ Lyon, INSA Lyon, CNRS, LaMCoS, UMR5259,
69621
Villeurbanne,
France
^{2}
AVNIR Engineering,
62 Boulevard Niels Bohr,
69603
Villeurbanne Cedex,
France
^{*} email: yvon.briend@insalyon.fr
Received:
15
January
2021
Accepted:
30
April
2021
Onboard rotating machinery subject to multiaxial excitations is encountered in a wide variety of hightechnology applications. Such excitations combined with mass unbalance forces play a considerable role in their integrity because they can cause parametric instability and rotor–stator interactions. Consequently, predicting the rotordynamics of such machines is crucial to avoid triggering undesirable phenomena or at least limiting their impacts. In this context, the present paper proposes an experimental validation of a numerical model of a rotorshafthydrodynamic bearings system mounted on a moving base. The model is based on a finite element approach with Timoshenko beam elements having six degrees of freedom (DOF) per node to account for the bending, torsion and axial motions. Classical 2D rectangular finite elements are also employed to obtain the pressure field acting inside the hydrodynamic bearing. The finite element formulation is based on a variational inequality approach leading to the Reynolds boundary conditions. The experimental validation of the model is carried out with a rotor test rig, designed, built, instrumented and mounted on a 6DOF hydraulic shaker. The rotor’s dynamic behavior in bending, torsion and axial motions is assessed with base motions consisting of mono and multiaxial translations and rotations with harmonic, random and chirp sine profiles. The comparison of the predicted and measured results achieved in terms of shaft orbits, full spectrums, transient history responses and power spectral densities is very satisfactory, permitting the experimental validation of the model proposed.
Key words: Onboard rotor / experimental validation / hydrodynamic bearings / 6DOF shaker / multiaxial excitation / random motion
© Y. Briend et al., published by EDP Sciences 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Rotating machinery subject to base motions, socalled onboard rotor, is a common feature in many industrial fields. Examples include a helicopter turboengine subject to airflight maneuvers and spectral lines associated with rotor blades, a turbopropeller suspended from an aircraft wing excited by the aerodynamic forces of a broad frequency range, and a spatial turbopump undergoing pyrotechnic shocks. These systems fulfill vital functions that must be continuously ensured despite the internal and external excitations they undergo during their lifetime. Internal excitations include the unavoidable residual mass unbalances that generate forces synchronous with the rotation speed and which are responsible for resonance phenomena when passing through critical speeds. The external forces generated from the base motions are generally asynchronous and may create parametric excitation. In the presence of nonlinearities, these forces are able to trigger dynamic behaviors such as bifurcations, and quasiperiodic and chaotic motions. These phenomena may be sources of failures by fatigue [1,2] or may lead to rotor–stator contacts, with the destruction of labyrinth seals or turbine blades [3]. In order to ensure the integrity of onboard rotors, predicting the emergence of such phenomena is crucial to ensure they can be avoided or that their impacts can at least be limited. This goal can be achieved by the development of reliable onboard rotor modeling, rigorously validated by a complete series of experimental comparisons.
In the last decade, contributions devoted to the topic of onboard rotors have been relatively abundant. In particular, Dakel et al. [4] analyzed the bending vibrations of a simply supported rotor composed of a shaft with one disk, excited by both its base and the shaft asymmetry. Their numerical results in terms of stability charts and shaft orbits evidenced the interactions between the two types of parametric excitations. Afterwards, the same authors [5] focused on the same shaftdisk system, this time symmetrical, mounted on short journal bearings with the great advantage of having an analytical solution to the Reynolds equation. The following year, the lateral and torsional dynamics of a multishaft onboard rotor was examined by Han and Chu [6]. Their timehistory and spectral analysis revealed that the roll rotation of the base was responsible for a noticeable increase in system response and the occurrence of new resonances in the vicinity of the natural frequencies. Later, they determined in [7] the stability of a singleshaft rotor in bending only, parametrically excited by one or several harmonic rotations of its base. The influence of the phase between the combined base rotations acting with the same frequency was shown in stability charts. Thereafter, the nonstationary shaft’s speed of rotation was introduced numerically in onboard rotor applications by Bouziani et al. [8], who performed a numerical analysis of the bending behavior in terms of spectrograms. Saimi et al. [9] used a new discretization method with “hp” finite elements to model an onboard rotor with a cantilevered disk, undergoing constant base rotations which altered the Campbell diagram. Following the work done by Dakel et al. [4], onboard rotors supported by journal bearings have attracted much attention. Wang et al. [10] considered the case of a centrifugal pump modeled by a simplified six degrees of freedom (DOF) system and subject to the El Centro earthquake, Reddy et al. [11] introduced viscoelastic damping in the journal bearings of a rotor representing a turbocharger, and Liu et al. [12] compared two journal bearing models in terms of shaft orbits for a rotor with one disk excited by a harmonic base translation.
Several authors addressed support flexibility in combination with noninertial support motion. Vicencio and Cruz [13] examined the dynamics of rotating machinery used in the nuclear field, such as gas and steam turbines supported by hydrodynamic bearings mounted on flexible foundations. By imposing base translations defined with accelerograms corresponding to shock response spectrums from real seismic excitations, they deduced the impact of the nonlinearities stemming from the journal bearings and the importance of knowing the foundation characteristics of such rotors. Zhu et al. [14] analyzed an onboard hollow rotor mounted on hydrodynamic bearings included within rubber damping rings showing viscoelastic behavior. They first demonstrated the apparent passive vibration control related to this additional damping, in particular in the vicinity of the critical speeds, and secondly the influence of the base translations on the whole rotorbearingfoundationsystem. Sales et al. [15] modeled the bending, torsional and longitudinal effects of a rotor mounted on a flexible moving support represented by a Timoshenko beam. The impact of support flexibility was exhibited on the bending effect through the disk orbits and the Campbell diagram for harmonic base rotations composed of roll and yaw.
Lastly, a novel control law intended to master the parametric instabilities of an onboard rotor was proposed by Soni et al., for a shaft mounted on active magnetic bearings [16] and in the presence of an electromagnetic actuator [17]. Other contributions were devoted to onboard rotors with geometric nonlinearities caused by large shaft displacements. This was the case of Phadatare et al. [18] with their highly flexible shaft mounted on a base subject to sinusoidal vertical translations and of Shahgholi et al. [19] with their asymmetric shaft whose support was subject to sinusoidal transverse rotations. In both these works, the shaft dynamics was restricted to the first mode of the linear system to facilitate the application of perturbation methods. As performed previously by Dakel et al. [4], Yi et al. [20] introduced the shaft and disk asymmetries of an onboard rotor and afterwards demonstrated the frequency combinations stemming from the rotational base motions and the asymmetry through the spectral displacement response of the disk. Qiu et al. [21] focused on the dynamics of a planetary gear excited by harmonic yaw motions of its support. Sousa et al. [22] and Stanica et al. [23] carried out several numerical investigations such as the analysis of the influence of a constant support rotation on the variation of critical speeds. Other simplified models have been developed recently to assess the impact of the wave motions on the engine rotors of maritime vessels [24,25].
All the research of this last decade mentioned previously has dealt exclusively with numerical results. The experimental investigations performed after those of Lee et al. [26], Duchemin et al. [27] and Driot et al. [28] are however much rarer. Among the few cases was the work done by Sousa et al. [29] who built an onboard rotor test bench composed of a slender shaft with a constant circular cross section and a disk supported by ball bearings. This rotor was excited by its base with a monoaxial electrodynamic shaker reproducing shock and sinusoidal base translations. Their numerical model was validated in terms of timehistory disk displacements in the case of a shock profile for the rotor at rest and of a sinusoidal profile for the rotor operating at both sub and supercritical speeds. Chen et al. [30] designed a support for combining simultaneous rotations around two orthogonal axes. Nonetheless, they excited their rotor with only one roll rotation of the base since the shaker they used was monoaxial. The validation of their model was carried out on a harmonic profile from the disk displacements in the time domain. More recently, Jarroux et al. [31] mounted a rotor equipped with active magnetic bearings on a 6DOF hydraulic shaker, the bearings being included within touchdown bearings with elastic dampers. This rotor was subjected to a harmonic vertical translation of its support for several acceleration levels until the occurrence of rotor–stator contacts within the active magnetic bearings. Their model was validated in terms of shaft orbits using numericalexperimental comparisons.
Therefore, strong interest in onboard investigations can be found in the literature. Nevertheless, most of the works remain numerical, with assumptions that are often far from those occurring in real industrial systems. This is the case for instance with the papers devoted to journal bearings [5,10–13,24,25] which exploited the short bearing hypothesis so as to use analytical solutions of the Reynolds equation and thus facilitate the numerical investigations. Furthermore, the base excitations were still too restricted to monoaxial and monofrequency cases as well as limited to the shaft bending dynamics without accounting for the axial and torsional effects [5,7–11,14]. This also holds for the experimental works, mostly owing to test benches using monoaxial shakers that moreover do not permit a straightforward change of position of the instantaneous axis of the base rotation. In addition, the experimental validations have been exclusively performed with shafts mounted on ball bearings rather than on hydrodynamic bearings for which numerical results are however abundant. The purpose of the present paper is thus to overcome all these issues by proposing a novel and complete experimental study in order to validate an onboard rotor model developed previously in [32]. The latter has six DOFs per node to account for all bending, torsion and axial shaft motions. The rotor is mounted on finitelength hydrodynamic bearings. Simultaneous translational and rotational base motions are properly defined to be able to correctly represent the multiaxial experimental excitations. Experimental validations of the linearized model are then carried out afterwards for one subcritical shaft speed of rotation, not only with monoaxial base excitations but also with multiaxial and multifrequency base excitations. For the sake of completeness, several time profiles of base motions are imposed such as harmonic, random and chirp sine. The bending, axial and torsional shaft dynamics are all addressed successively. In order to obtain an experimental validation of the model proposed, a test rig of an onboard rotor has been designed, built and instrumented accordingly. A 6DOF hydraulic shaker is used to subject the rotor to base excitations.
2 Modeling and base motion definition
A typical onboard rotor is sketched in Figure 1a. Two frames of reference are necessary to represent its dynamics. The first one, , is considered as Galilean. The second one, , is fixed to the rotor base and thus independent of the shaft’s rotation, with its origin O located at the left end of the shaft’s neutral axis and y⃗ along the shaft’s axis of rotation at rest. A complete description and model of the rotor dynamics by a finite element method was proposed in [32]. By wayof a brief recall, the shaft is modeled with 1D Timoshenko beam elements having two nodes and six DOFs per node which are (u, w) the bending transverse displacements along (x⃗, z⃗) respectively,v the axial displacement in y⃗, (ψ, θ) the bending transverse rotations around (z⃗, x⃗) respectively, and β the torsion angle of rotation around y⃗. All these displacements depend on only the axial coordinate y, and the total shaft speed of rotation is defined as where is the constantnominal shaft speed of rotation provided by the motor. Cubic shape functions are employed to approximate (u, w) while linear shape functions are used for v and β. The following system of equations governing the rotor dynamics can be reached in a matrix form by applying the Lagrange’s equations to all the energies and virtual works of the system [32]: (1)
with M, G and K being the classical mass, gyroscopic and structural stiffness matrices, respectively, F_{u} the classical mass unbalance force vector and δ the global DOF vector expressed with respect to the static equilibrium position. C_{ba}, K_{ba} and F_{ba} are the timedependent damping and stiffness matrices and the force vector, respectively, related to the base motion. F_{gl} is a force vector related to the gravity load oriented along X⃗_{0}, which is nonzero and timedependent only in the case of a rotational base motion [32] (pitch or roll only, given here the gravity is along the yaw axis). It is obtained as F_{gl} = F_{gd} −F_{gs} where F_{gd} is the full dynamic gravity load vector and F_{gs} is the static (without base motion) gravity load vector. The remaining terms are related to the linearization of the journal bearing restoring forces in thevicinity of the static equilibrium position at a given constant shaft speed of rotation characterized by the static solution vector δ_{s}. In this context,C_{jb} and K_{jb} are respectively the damping and stiffness matrices derived from these restoring forces with respect to and δ. The computation of the pressure fields providing the restoring forces of the journal bearings are obtained using a 2D FEM with 4node rectangular elements, taking into account the cavity problem by applying the Reynolds boundary conditions as in [33]. The only damping is provided by either the hydrodynamic bearings or the motorshaft coupling. Unlike in [32], nonstationary effects are not accounted for. Furthermore, the mass unbalance is considered small enough so that the additional matrices related to mass unbalance are neglected.
Given the aim of the present paper to deal with multiaxial excitations, a complete description of the rotor base kinematics is essential, especially in the case of combined translational and rotational motions where the rotations are defined around a fixed axis that may not necessarily pass through O. Let the base be rigid but mobile. To begin with the rotations only, the orientation of the moving frame with respect to is done with a classical approach for rigid solids using the Euler angles. Thus, is obtained from by three successive rotations: a pitch rotation of angle α_{1} around Z_{0} that creates a first intermediary frame with Z⃗_{1} =Z⃗_{0}, a yaw rotation of angle α_{2} around X⃗_{1} that creates a second intermediary frame with X⃗_{2} =X⃗_{1}, and a third roll rotation of angle α_{3} around Y⃗_{2} that creates with Y⃗_{2} =y⃗. The translation motion of with respect to is defined through the motion of the origin O of in by the vector .
In order to now impose base motions composed of simultaneous translations and rotations, one possible approach decomposes the variables (X_{0}, Y_{0}, Z_{0}) such that: (2a)
where (X_{0,rot}, Y_{0,rot}, Z_{0,rot}) and (X_{0,tr}, Y_{0,tr}, Z_{0,tr}) are the variations of (X_{0}, Y_{0}, Z_{0}) related to pure rotational and translational motions, respectively. In this way, it is possible to interpret the whole motion as a combination of a translation in of a fictitious point I and a rotation around a combination of the axes (Z⃗_{0}, X⃗_{1}, Y⃗_{2}) passing through this point I. This is illustrated in Figure 1b with an inplane example composed of two translations along X⃗_{0} and Y⃗_{0} and a rotation around Z⃗_{0}. The previous variables (X_{0,rot}, Y_{0,rot}, Z_{0,rot}) are then related to the constant coordinates (x_{I}, y_{I}, z_{I}) of I in the moving frame by: (3)
while the variables (X_{0,tr}, Y_{0,tr}, Z_{0,tr}) represent the coordinates of I in with the orthogonal transformation matrix from to in which appear cosine and sine functions of the angles (α_{1}, α_{2}, α_{3}). If (X_{0,tr}, Y_{0,tr}, Z_{0,tr}) are nil or constant, then the actual instantaneous axis of rotation of with respect to passes through I.
Therefore, according to the previous approach, the complete motion of is defined by three sets of parameters: (1) (X_{0,tr}(t), Y_{0,tr}(t), Z_{0,tr}(t)) the variables that produce the translational motion along the fixed axes of the Galilean frame , referred to as pure translations in the sequel; (2) (α_{1}, α_{2}, α_{3}) the Euler angles that produce the rotational motion of and all fixed points in around an axis passing through point I; (3) (x_{I}, y_{I}, z_{I}) the constant coordinates of I in . The matrices C_{ba}, K_{ba} and vector F_{ba} of equation (1) will then depend on these variables, as can be found in [32].
Fig. 1
(a) Sketch of an onboard rotor and (b) Example of an inplane motion of the moving frame , composed of two translations along X⃗_{0} and Y⃗_{0} and a rotation of angle α_{1} and radius C. 
3 Experimental validation
3.1 Presentation of the onboard rotor test bench
As mentioned in Section 1, experimental investigations of onboard rotors are not abundant in the literature, even less so when the shaft is supported by hydrodynamic bearings and subject to multiaxial motions. Nonetheless, this type of bearing is widely used in industry and the excitations experienced by onboard rotating machines are often naturally multiaxial. These reasons thus motivated the design and fabrication of the rotor test bench presented in Figure 2. This rotor is composed of a steel shaft, two steel disks, two hydrodynamic bearings, a flexible coupling and a DC electric motor providing a maximum torque of 0.177 N.m. All the numerical data are available in Table 1. The whole system is mounted on a thick plate, itself fixed to the upper face of the 6DOF hydraulic shaker driven in closed loop by means of accelerations. The latter is a cubeshaped longstroke shaker, with an edge of 0.8128 m and offering accelerations up to 10 g within a frequency range of [0–250] Hz. The motorshaft coupling (model 321.12.2222 of Huco^{®}) characteristics are available in Table 2. These characteristics come either from the manufacturer or from an optimization performed through preliminary modal analysis with a fixed base.
Before any motion occurs, the Galilean frame always has the same orientation as the moving frame so that the angles (α_{1}, α_{2}, α_{3}) do not present any nonzero 0Hz static components. However, the location of I in may depend on the test. In the sequel, unless otherwise specified, it is set by x_{I} = y_{I} = z_{I} = 0. The shaft is discretized into 27 beam finite elements, as shown in Figure 3. The shaft radius is 6.35 mm for the FEs 4 to 24 and 6 mm for the FEs (1, 2, 3, 25, 26, 27). In order to model the shaft local stiffening due to the tightening of the disks on the shaft, angular stiffness coefficients of 100 N.m/rad are added to the coefficients related to the DOFs ψ and θ of Nodes 9 and 21 in the structural stiffness matrix K of equation (1). These new coefficients were found by a preliminary modal analysis, the shaft having freefree boundary conditions, and proved to provide better agreements with the experimental results even with the supported shaft. The hydrodynamic bearings have finite lengths with a length/diameter ratio of 0.4851. They are fed by two circular feeding holes located in the axial middle at 0 ° and 180 ° with respect to z⃗. The hydrodynamic bearings are discretized using a 2D rectangular mesh with 200 nodes and 15 nodes in the circumferential and axial directions, respectively. The linearized stiffness and damping matrices at the operating speed rpm are summedin Table 3.
In order to assess the whole shaft dynamics, the rotor is equipped with several types of sensor. Firstly, four radial proximity probes (inductive sensors with eddy current) are distributed in two probe supports to monitor the shaft bending motion. These supports are localized at Nodes 15 and 19, i.e. relatively far from the bearings, because of constraints regarding the available space and to maximize the measurement amplitudes. Each probe was rigorously calibrated with respect to the specific circular shape and material of the shaft. These sensors are directed along either x⃗ or z⃗ of . Secondly, anaxial sensor of the same model is set on one shaft end (Node 28) in order to measure the axial rigidbody motion. Thirdly, two optical encoders with a resolution of 4096 line counts are installed at Nodes 2 and 27 to access the shaft torsional motion. The stator part of these sensors is not fixed to the stator nor any support so as to avoid additional localized stiffness of the shaft. Their weight of 215 g each is then fully supported by the shaft itself and they can follow its deflection so that their presence affects the original bending dynamics. However, the rotation of the stator part is prevented with an axial screw. These encoders are only used and mounted on the shaft for the test in Section 3.4. Finally, another proximity probe directed along x⃗ is used as a keyphasor probe (tachometer) which detects an axial groove located on the shaft near Node 4. Regarding the 6DOF shaker, it is equipped with six monoaxial accelerometers for its driving process and with four triaxial accelerometers to retrieve the real shaker motion to be introduced into the onboard rotor MEF for validation purposes. The process to obtain this real 6DOF motion is described indepth in [34].
The Campbell diagram of the rotor test bench on a fixed base and without encoders is computed and presented in Figure 4. The tracking of the mode shapes is done by applying the NC^{2}O criterion [35]. The solid, dashed and dashdotted curves refer to the bending, axial and torsional modes, respectively. The first four critical speeds are [2084, 2281, 4749, 5402] rpm. The natural frequencies and the modal damping factors of each mode appearing in the Campbell diagram at rpm are listed in Table 4. The first two modes of the latter, i.e. the axial and torsion ones, are almost only rigidbody modes with a main localized deformation in the coupling. The higher modes, after the seventh one, are not displayed since their natural frequencies are higher than 250 Hz which is above the frequency range under study for the following test. As could be expected, the additional mass from the encoders tends to lower the natural frequencies related to bending and axial dynamics.
Regarding the further experimental validations, the shaft is rotating in the direct sense, i.e. from z⃗ to x⃗, at a constant nominal speed of rpm. The keyphasor and proximity probe signals have the same sampling frequency, fixed to 4096 Hz, in order to accurately determine the shaft angle of rotation at any time. Thus, the shaft angular position as well as the phase between mass unbalance and base motion can be introduced in the numerical model. The initial phases of the mass unbalance and the base motions are given with respect to the keyphasor pulse. Finally, except for the torsion measurements in Section 3.4, the measurements presented below are not filtered or processed with any data processing tools (including windowing).
Fig. 2
Onboard rotor test bench used for the experimental validations: (a) CAD view and (b) Photo view. 
Main rotor data.
Motorshaft coupling characteristics of the rotor test bench.
Fig. 3
Onboard rotor FE mesh. 
Nonzero coefficients of the linearized stiffness and damping matrices of the hydrodynamic bearings at rpm.
Fig. 4
Campbell diagram of the rotor test bench. 
3.2 Multiaxial harmonic excitation
The first case of base excitation focuses on multiaxial harmonic translations and on bending dynamics. The motion comprises two translations: the first one along X⃗_{0} defined such that with A_{X} = 15 μm, f_{X} = 42.5 Hz and φ_{X} = 57 ° and the second one along Z⃗_{0} defined such the with A_{Z} = 15 μm, f_{Z} = 56.7 Hz and φ_{Z} = 91 °. The frequencies were chosen intentionally as harmonics of the shaft speed of rotation frequency X (f_{X}≡ 1.5X and f_{Z} ≡ 2X) to obtain a periodic shaft dynamics and periodic orbits. In fact, the shaft speed of rotation delivered by the motor is slightly lower than rpm, inducing a slow variation of the experimental orbits with time. This may be interpreted as a variation of the phase between the mass unbalance and the base motion.
The raw (without filter) experimental and numerical orbits of Nodes 15 and 19 are shown in Figure 5a and Figure 5b, respectively, within a time interval of 2 s (≈ 57 full shaft rotations). In addition, the shaft position at the experimental and numerical keyphasor pulses is superimposed on the respective shaft orbits. Excellent agreement can be noticed between the two sets of results in terms of orbit shape, amplitude and phase. Furthermore, the shaft occupies the same position at the keyphasor pulses (two pulses are seen in the orbits since a period is obtained for two full shaft rotations owing to the 1X frequency of excitation). However, some small deviations may be noted in amplitude in both Nodes 15 and 19. This can be attributed mainly to the presence of a residual mass unbalance since initial shaft balancing is never perfect, but also to some likely weaknesses in the modeling of the hydrodynamic bearings. It is noteworthy that the two orbits of Nodes 15 and 19 are relatively similar to one another. This is due to the fact that the frequencies of the base translations involved tend to excite the first two bending modes more (Modes 3 and 4 of Tab. 4 which are characterized by a maximum displacement in the middle of the shaft). Given the position of the radial probes supports, these modes do not induce any phase shift between these two orbits.
The full spectrum [36] of these experimental and numerical orbits are shown in Figure 6. They are obtained from the time signal on 12 s once the steady state is reached and normalized by the signals length (=12 s × 4096 Hz). On looking at these new results, it can be seen that the accuracy of the numerical model is also established in the frequency domain. Besides, the amplitude discrepancy of the orbits of Nodes 15 and 19 mostly occurs at the shaft speed rotation frequency (± 28.3 Hz) rather than on the base motion frequencies. This might mean that this discrepancy is due more to the residual mass unbalance than to weaknesses in the modeling of the hydrodynamic bearings. Other deviations may be noticed in some harmonics of the shaft speed of rotation, mostly present experimentally, for instance at 14.2 Hz (0.5X) or 85 Hz (3X). These deviations may be explained by certain asymmetric properties either in the shaft or disks, nonlinearities in the hydrodynamic bearings, cyclic phenomena in the motor (ball bearings), etc. Lastly, the whirl of all the orbits is mainly forward since the peaks at +28.3 Hz (1X), +42.5 Hz (1.5X) and +56.7 Hz (2X) are higher than those at the corresponding negative frequencies.
Fig. 5
Experimentalnumerical comparison of shaft orbits in response to mass unbalance and two harmonic translations along X⃗_{0} and Z⃗_{0} within a time interval of 2 s (≈ 57 full shaft rotations): (a) Node 15 and (b) Node 19. 
Seven first numerical modes of the rotor test bench at rpm (FW: forwardwhirl, BW: backward whirl, WE: with encoders, ∅E: without encoders).
Fig. 6
Full spectrum of the numerical and experimental orbits of Figure 5: (a) Node 15 and (b) Node 19. 
3.3 Multiaxial random excitation
In this section, the numerical model is tested with more general motions of the random type. Both bending and axial dynamics are assessed here. The shaker’s motion is inplane with only three DOFs involved: the translation along X⃗_{0}, the translation along Y⃗_{0} and the pitch α_{1} around Z⃗_{0}. The same power spectral density (PSD) is built for the only two monoaxial accelerometers involved in the driving process of the 6DOF shaker, with a theoretically nil crosscorrelation between each other. This PSD is defined by two plateaus at 2.51× 10^{−4} g^{2}/Hz in [2031.6] Hz and [4560] Hz and a third one at 2.51 × 10^{−6} g^{2}/Hz in [3540] Hz. This specific profile permits avoiding the resonance of the axial mode and the first two bending modes (rows 1, 3 and 4 of Tab. 4) so that the rotor response remains small enough and thus linear. The PSDs in acceleration of the corresponding six DOFs of the shaker, obtained by the process described in [34] and computed on the random part of the test between 5 s and 19 s (for a full test occurring between 0 s and 22 s), are presented in dB in Figure 7a for (Ẍ_{0,tr}, Ÿ_{0,tr}, ) and Figure 7b for , , . The real base motion is mainly composed of the three targeted DOFs. The peak around 28.3 Hz present in the other three DOFs (, , ) corresponds to the shaft frequency of rotation which spreads in the rotor base, but remains low enough.
The numerical and experimental timehistory responses in (u, w) for Node 15 and in − v (the minus sign is because the axial sensor is directed along −y⃗) for Node 28 are presented in Figure 8 in a short time interval for the sake of clarity. The corresponding PSDs, computed in the random part of the test in t ∈ [5;19] s, are depicted in Figure 9. Firstly, as expected, the amplitude of u in Figure 8 shows random variations while the amplitude of w is quite constant and unaffected by the base excitation, which is in agreement with the level of the corresponding translational base accelerations Ẍ_{0,tr} and in Figure 7a. Secondly, as for the harmonic case of Section 3.2, the amplitudes of u and especially w are slightly overestimated numerically. This overestimation indeed mostly occurs at the shaft frequency of rotation (see Fig. 9) in response to the residual mass unbalance, as in the previous harmonic case. A small peak around 10.4 Hz, which is related to the first torsion mode, is predicted numerically because of a coupling between bending and torsion in the damping matrix C_{jb}. Nonetheless,since this peak does not appear experimentally, this coupling is likely to be lower than expected. Another peak around 56.7 Hz (2X) is found in the measurements but is not predicted as the model considers only an axisymmetric shaft and disks. Regarding the axial displacement, remarkable agreement is seen in the whole frequency range of the study. In this frequency range, the shaft responds on its first axial rigidbody mode in the same way as for a 1DOF massdamperspring system (all the deformation acts on the motorshaft coupling).
Fig. 7
PSD in acceleration of the 6 DOFs of the rotor base for the random excitation: (a) pure translational DOFs and (b) angular DOFs. 
Fig. 8
Experimental (Exp.) and numerical (Num.) comparisons of displacements (u, w) of Node 15 and − v of Node 28 for a random base excitation. 
3.4 Monoaxial chirp sine excitation
So far, only the bending and axial shaft dynamics have been evaluated. The encoders were not installed in the previous cases since they may alter the original shaft dynamics and introduce additional sources of uncertainties. For instance, the balancing is perturbed, the axial rigidbody mode is deeply damped and several harmonics of the shaft speed of rotation appear due in particular to the ball bearings inside the encoders. In this context, a test focusing only on the torsion dynamics was conducted. Consequently, a monoaxial roll (around y⃗) of the rotor base was created, with a targeted instantaneous axis of rotation defined to be coaxial with the shaft axis of rotation (i.e. with x_{I} = y_{I} = z_{I} = 0) in order to minimize the bending and axial responses. However, maintaining this coaxiality efficiently with a random excitation profile through the shaker driving process is quite difficult. Moreover, a monofrequency harmonic excitation does not permit distinguishing the torsional shaft response from the bending one precisely. For these reasons, a chirp sine profile of the base rotation waschosen.
Since the first torsional mode is almost a pure coupling mode with no shaft deformation, it can hardly be detected by the encoders both mounted on the shaft. Therefore, the second torsional mode was targeted in this study. To this end, the roll rotation is defined by an angle with , f_{i} = 10 Hz the initial frequency, f_{f} = 125 Hz the final frequency and T = 16 s. In this way, the instantaneous frequency f_{3}(t) = f_{i} k^{t} of the roll rotation increases exponentially from f_{i} at t = 0 to f_{f} at t = T. In order to maintain a relatively constant amplitude of angular acceleration, A_{3} is defined so that it decreases with time such that with A_{0} = 1200° / s^{2}. The time variations of the targeted and real (experimental) α_{3} are presented in Figure 10. Despite a real amplitude slightly higher than expected, it can be seen that the chirp sine profile is globally satisfactory. The other real DOFs of the shaker, not shown here, are also present in practice with maximum accelerations around 0.15 g, 0.25 g, 0.2 g, 250° / s^{2}, 60° / s^{2} for Ẍ_{0,tr}, Ÿ_{0,tr}, , and , respectively. Although they are not nil, they are not high enough to perturb the monoaxial motion.
To remove the rigidbody rotation ϕ^{*} of the shaft and emphasize the torsional dynamics, the experimental response is established in terms of an angle B obtained as the difference of the total shaft angle of rotation ϕ between Nodes 27 and 2 (namely B = ϕ_{27} − ϕ_{2} = β_{27} − β_{2}) where the two encoders are located. Moreover, a highpfilter with a cutoff frequency of 5 Hz is applied to the experimental angle B. The corresponding result is presented in Figure 11 in terms of its PSD, computed between 2.3 s and 17.4 s. The response of angle B is also plotted without base motion for two different shaft speeds of rotation to emphasize the impact of the base on the torsional response. The two corresponding PSDs are computed from the same time duration of 15 s as for the PSD with base motion.
Firstly, it is noteworthy that phenomena related to bending dynamics were observed experimentally around 28.25 Hz and 56.75 Hz, which correspond to the harmonics 1X and 2X of the mass unbalance forces, and around 35.75 Hz and 81.25 Hz, which correspond to Modes 4 and 5 of Table 4 (with the encoders). These frequencies, that should, theoretically, not be detected by the encoders, are present since the shaft deflections inherently induce a rotational motion of the stator part of the encoders (see Sect. 3.1 for more details). Secondly, a peak at 93.75 Hz related to the second torsional mode appears. This peak is not related to any bending phenomena, because it is missing from the measured radial deflection, as seen for instance in Figure 6, and it is insensitive to the speed of rotation as seen with the dashed lines of Figure 11 in the case of no base excitation. Thus, it canbe stated that it is related to the second torsion mode, well predicted by the model, see Mode 7 at 92.50 Hz in Table 4. On the other hand, erasing the base excitation leads to an attenuation of 22 dB on the 93.75 Hz torsion peak: the base roll excitation magnifies the torsional behavior of the rotorbearing system.
Fig. 10
Comparison between the targeted and real roll angle α_{3} of the rotor base with a chirp sine profile. 
Fig. 11
Experimental comparison of angle B = β_{27} − β_{2} with and without base roll excitation in chirp sine. 
4 Conclusion
This work addressed the experimental validation of an onboard rotor model supported by finitelength hydrodynamic bearings and subject to a wide variety of excitations, such as translations and/or rotations, mono or multiaxial, harmonic or random or in chirp. To this end, a rotor test bench was designed, built, instrumented and mounted on a 6DOF hydraulic shaker. The real excitations produced by the latter were introduced in the FEM by solving an inverse problem to obtain better concordance between the calculated and measured results. All the shaft dynamics, i.e. bending, torsion and axial motions, were assessed. Excellent agreements were obtained in terms of shaft orbits and full spectrums for the bending dynamics with the multiaxial harmonic base excitation. However, some slight discrepancies in the amplitudes were noted and could mainly be assigned to the presence of residual mass unbalance in the experiment and possibly to certain limitations in the hydrodynamic bearing modeling (temperature and viscosity variations in the fluid, angular misalignment, bearing roughness, etc.). The axial dynamics was predicted with great accuracy in both the time and frequency domains, even in a complex case of random excitation. Regarding torsion, a roll base rotation in chirp sine permitted demonstrating experimentally the response of the second mode, whose frequency matched that of the model.
Conflict of interest
None.
Acknowledgements
This work was supported by the French National Research Agency(ANR) in the framework of ANRLabcomSME AdViTAM, 16LCV10006 and Equipex PHARE 10EQPX0043.
References
 J. Zhao, I. Linnett, L. McLean, Stability and bifurcation of unbalanced response of a squeeze film damped flexible rotor, J. Tribol. 116, 361–368 (1994) [Google Scholar]
 J. InayatHussain, H. Kanki, N. Mureithi, On the bifurcations of a rigid rotor response in squeezefilm dampers, J. Fluids Struct. 17, 433–459 (2003) [Google Scholar]
 H. Ma, F. Yin, Y. Guo, X. Tai, B. Wen, A review on dynamic characteristics of blade–casing rubbing, Nonlinear Dyn. 84, 437–472 (2016) [Google Scholar]
 M. Dakel, S. Baguet, R. Dufour, Steadystate dynamic behavior of an onboard rotor under combined base motions, J. Vib. Control 20, 2254–2287 (2014) [CrossRef] [Google Scholar]
 M. Dakel, S. Baguet, R. Dufour, Nonlinear dynamics of a supportexcited flexible rotor with hydrodynamic journal bearings, J. Sound Vib. 333, 2774–2799 (2014) [CrossRef] [Google Scholar]
 Q. Han, F. Chu, Dynamic behaviors of a geared rotor system under timeperiodic base angular motions, Mech. Mach. Theory 78, 1–14 (2014) [CrossRef] [Google Scholar]
 Q. Han, F. Chu, Parametric instability of flexible rotorbearing system under timeperiodic base angular motions, Appl. Math. Model. 39, 4511–4522 (2015) [CrossRef] [MathSciNet] [Google Scholar]
 R. Bouziani, N. Ouelaa, Simulation of the dynamic behavior of a rotor subject to base motion under variable rotational speed, Mech. Ind. 18, 308 (2017) [CrossRef] [Google Scholar]
 A. Saimi, A. Hadjoui, An engineering application of the hp version of the finite elements method to the dynamics analysis of a symmetrical onboard rotor, Eur. J. Comput. Mech. 25, 388–416 (2016) [CrossRef] [Google Scholar]
 R. Wang, X. Guo, Y. Wang, Nonlinear analysis of rotor system supported by oil lubricated bearings subjected to base movements, Proc. Inst. Mech. Eng. C 230, 543–558 (2016) [CrossRef] [Google Scholar]
 M.R. Reddy, J. Srinivas, Vibration analysis of a support excited rotor system with hydrodynamic journal bearings, Proc. Eng. 144, 825–832 (2016) [CrossRef] [Google Scholar]
 Z. Liu, Z. Liu, Y. Li, G. Zhang, Dynamics response of an onboard rotor supported on modified oilfilm force considering base motion, Proc. Inst. Mech. Eng. C 232, 245–259 (2018) [CrossRef] [Google Scholar]
 F. Vicencio, E.F. Cruz, A high order nonlinear study to evaluate the seismic response of rotating machines–structure–soil foundation systems, J. Earthquake Eng. 0, 1–33 (2019) [Google Scholar]
 H. Zhu, W. Chen, R. Zhu, J. Gao, M. Liao, Study on the dynamic characteristics of a rotor bearing system with damping rings subjected to base vibration, J. Vib. Eng. Technolog. 8, 121–132 (2020) [CrossRef] [Google Scholar]
 T.d.P. Sales, E. Spuldaro, L.F. Damy, D.A. Rade, Dynamic modeling of flexible rotors mounted on an elastic base undergoing arbitrary attitude motion, in: International Conference on Rotor Dynamics. (Springer, 2018), pp. 562–576 [Google Scholar]
 T. Soni, J.K. Dutt, A. Das, Parametric stability analysis of active magnetic bearingsupported rotor system with a novel control law subject to periodic base motion, IEEE Trans. Ind. Electr. 67, 1160–1170 (2019) [Google Scholar]
 T. Soni, A. Das, J. Dutt, Active vibration control of ship mounted flexible rotorshaftbearing system during seakeeping, J. Sound Vib. 467, 115046 (2020) [Google Scholar]
 H. Phadatare, B. Choudhary, B. Pratiher, Evaluation of nonlinear responses and bifurcation of a rotorbearing system mounted on moving platform, Nonlinear Dyn. 90, 493–511 (2017) [Google Scholar]
 M. Shahgholi, G. Payganeh, Forced vibrations of nonlinear symmetrical and asymmetrical rotating shafts mounted on a moving base, ZAMM J. Appl. Math. Mech. 99, e201700097 (2019) [Google Scholar]
 Y. Yi, Z. Qiu, Q. Han, The effect of timeperiodic base angular motions upon dynamic response of asymmetric rotor systems, Adv. Mech. Eng. 10, 1687814018767172 (2018) [Google Scholar]
 X. Qiu, Q. Han, F. Chu, Dynamic modeling and analysis of the planetary gear under pitching base motion, Int. J. Mech. Sci. 141, 31–45 (2018) [Google Scholar]
 M. Sousa, V. Del Claro, A. Cavalini, V. Steffen, Numerical investigation on the dynamic behavior of an onboard rotor system by using the fem approach, J. Braz. Soc. Mech. Sci. Eng. 39, 2447–2458 (2017) [Google Scholar]
 C.M. Stanica, M.V. Predoi, I. Stroe, Study of rotating machineries in a noninertial reference frame subjected to rotations, Romanian J. Acoustics Vib. 16, 125–136 (2019) [Google Scholar]
 Y. Han, M. Li, Nonlinear dynamic characteristics of marine rotorbearing system under heaving motion, Shock Vib. 2019 (2019) [Google Scholar]
 W. Zhao, M. Li, Y. Liu, Nonlinear dynamics of marine rotor system coupled with air bagfloating raft subjected to the basement excitations in lateral directions, Shock Vib. 2020 (2020) [Google Scholar]
 A.S. Lee, B.O. Kim, Y.C. Kim, A finite element transient response analysis method of a rotorbearing system to base shock excitations using the statespace newmark scheme and comparisons with experiments, J. Sound Vib. 297, 595–615 (2006) [CrossRef] [Google Scholar]
 M. Duchemin, A. Berlioz, G. Ferraris, Dynamic behavior and stability of a rotor under base excitation, J. Vib. Acoustics 128, 576–585 (2006) [CrossRef] [Google Scholar]
 N. Driot, C.H. Lamarque, A. Berlioz, Theoretical and experimental analysis of a baseexcited rotor, J. Comput. Nonlinear Dyn. 1, 257–263 (2006) [Google Scholar]
 M. Sousa, V. Del Claro, A. Cavalini, V. Steffen, Experimental validation of an onboard rotor fe model, in: Proceedings of the 24th ABCM International Congress of Mechanical Engineering, COBEM, Curitiba, Parana, Brazil (2017) [Google Scholar]
 L. Chen, J. Wang, Q. Han, F. Chu, Nonlinear dynamic modeling of a simple flexible rotor system subjected to timevariable base motions, J. Sound Vib. 404, 58–83 (2017) [Google Scholar]
 C. Jarroux, J. Mahfoud, R. Dufour, F. Legrand, B. Defoy, T. Alban, Dynamic behavior of a rotorAMB system due to strong base motions, in International Conference on Rotor Dynamics (Springer, 2018), pp. 340–349 [Google Scholar]
 Y. Briend, M. Dakel, E. Chatelet, M.A. Andrianoely, R. Dufour, S. Baudin, Effect of multifrequency parametric excitations on the dynamics of onboard rotorbearing systems, Mech. Mach. Theory 145, 103660 (2020) [Google Scholar]
 T. Zheng, N. Hasebe, Calculation of equilibrium position and dynamic coefficients of a journal bearing using free boundary theory, ASME J. Tribol. 122, 616–621 (2000) [Google Scholar]
 Y. Briend, E. Chatelet, R. Dufour, F. Legrand, S. Baudin, Identification of real translational and rotational displacements of sixaxial shakers with only six measured linear accelerations, Mech. Syst. Signal Process. 154, 107584 (2021) [Google Scholar]
 G. Mogenier, T. Baranger, G. Ferraris, R. Dufour, L. Durantay, A criterion for mode shape tracking: application to Campbell diagrams, J. Vibrat. Control 20, 179–190 (2014) [Google Scholar]
 P. Goldman, A. Muszynska, Application of full spectrum to rotating machinery diagnostics, Orbit 20, 17–21 (1999) [Google Scholar]
Cite this article as: Y. Briend, E. Chatelet, R. Dufour, M.A. Andrianoely, F. Legrand, S. Baudin, Dynamics of onboard rotors on finitelength journal bearings subject to multiaxial and multifrequency excitations: numerical and experimental investigations, Mechanics & Industry 22, 35 (2021)
All Tables
Nonzero coefficients of the linearized stiffness and damping matrices of the hydrodynamic bearings at rpm.
Seven first numerical modes of the rotor test bench at rpm (FW: forwardwhirl, BW: backward whirl, WE: with encoders, ∅E: without encoders).
All Figures
Fig. 1
(a) Sketch of an onboard rotor and (b) Example of an inplane motion of the moving frame , composed of two translations along X⃗_{0} and Y⃗_{0} and a rotation of angle α_{1} and radius C. 

In the text 
Fig. 2
Onboard rotor test bench used for the experimental validations: (a) CAD view and (b) Photo view. 

In the text 
Fig. 3
Onboard rotor FE mesh. 

In the text 
Fig. 4
Campbell diagram of the rotor test bench. 

In the text 
Fig. 5
Experimentalnumerical comparison of shaft orbits in response to mass unbalance and two harmonic translations along X⃗_{0} and Z⃗_{0} within a time interval of 2 s (≈ 57 full shaft rotations): (a) Node 15 and (b) Node 19. 

In the text 
Fig. 6
Full spectrum of the numerical and experimental orbits of Figure 5: (a) Node 15 and (b) Node 19. 

In the text 
Fig. 7
PSD in acceleration of the 6 DOFs of the rotor base for the random excitation: (a) pure translational DOFs and (b) angular DOFs. 

In the text 
Fig. 8
Experimental (Exp.) and numerical (Num.) comparisons of displacements (u, w) of Node 15 and − v of Node 28 for a random base excitation. 

In the text 
Fig. 9
PSD of the timehistory responses of Figure 8 computed for t ∈ [5; 19] s. 

In the text 
Fig. 10
Comparison between the targeted and real roll angle α_{3} of the rotor base with a chirp sine profile. 

In the text 
Fig. 11
Experimental comparison of angle B = β_{27} − β_{2} with and without base roll excitation in chirp sine. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.