Next Article in Journal
An Operational Image-Based Digital Twin for Large-Scale Structures
Previous Article in Journal
Predictive Maintenance for Remanufacturing Based on Hybrid-Driven Remaining Useful Life Prediction
Previous Article in Special Issue
Stability Analysis of Rotor-Bearing Systems under the Influence of Misalignment and Parameter Uncertainty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linear and Nonlinear Performance Analysis of Hydrodynamic Journal Bearings with Different Geometries

by
Carlos Alberto Alves Viana
,
Diogo Stuani Alves
and
Tiago Henrique Machado
*
Laboratory of Rotating Machinery, School of Mechanic Engineering, UNICAMP 200, Rua Mendeleyev, Postal Box 6122, Campinas 13083-970, Brazil
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(7), 3215; https://doi.org/10.3390/app12073215
Submission received: 15 September 2021 / Revised: 8 November 2021 / Accepted: 23 November 2021 / Published: 22 March 2022
(This article belongs to the Special Issue Design and Optimization of Rotor Dynamics in Applications)

Abstract

:
In rotor dynamics, a traditional way of representing the dynamics of hydrodynamic bearings is using stiffness/damping coefficients. It is thus necessary to carry out a linearization of hydrodynamic forces around the shaft’s equilibrium position. However, hydrodynamic bearings have highly nonlinear nature, depending on operating conditions. Therefore, this paper discusses the applicability of these linear/nonlinear approaches using a computational model of the rotating system, where the finite element method is used for rotor modelling and the finite volume method for bearing calculation. The main goal is to investigate the boundaries for linear approximation of the hydrodynamic forces present in lobed hydrodynamic bearings, with the system operating under high loading conditions. Several numerical simulations were performed varying preload parameter and rotating speed. A comparison of the system’s responses, in time domain (shaft orbits) and frequency domain (full spectrum), is made for linear and nonlinear models. Results showed that trilobed bearings are more susceptible to nonlinearities, even in situations of smaller vibration amplitudes, while elliptical bearings are sensitive only under larger vibration amplitudes. These analyses are of great importance for mapping the influence of nonlinearities in different types of lobed hydrodynamic bearings with fixed geometry, varying the preload parameter to verify the influence on the system’s dynamic response. This study is important and serves as the basis for cases of monitoring and fault diagnosis (in the field of structural health monitoring) since it is crucial to distinguish what would be a fault signature or a standard nonlinear effect created by the use of hydrodynamic bearings.

1. Introduction

Hydrodynamic bearings play an important role in the dynamics of turbomachines and are widely used in industry because they allow high loads at low and high rotating speeds. The bearings properties highly influence the dynamic characteristics of a rotor, and instability phenomena tend to especially appear at high rotating speeds. The latest advances in the metal-mechanical industry have overcome most bearing machining problems. New materials and designs have brought increased bearing life, which are no longer the main reason for fault or malfunction in certain applications. Thus, the interactions of bearings with the surrounding mechanical system are drawing attention.
The dynamics of the bearing is mainly governed by support zones between its fixed and movable parts since these transmit all vibration between the fixed parts of the bearing and the moving parts of the system. Therefore, the oil film in these interfaces, mainly provided to increase the bearings’ life, completely separating the contact surfaces, has a primordial role in the dynamic interactions. This thin film will, therefore, absorb part of the transmitted energy.
Although the equation governing the behavior of the oil film inside the bearing, Reynolds’ equation, dates from the late nineteenth century, it was only after mid-twentieth century that its complete solution became viable, mainly due to the advent of computer science. The first use of computers for the solution of the Reynolds Equation, with the correct boundary conditions, was achieved in 1956 by Pinkus. Pinkus analyzed the pressure distribution in elliptical bearings [1], partial arc [2] and three-lobe bearings [3], experimentally investigating the instability of these geometries. From then, with very little time, a wide spectrum of solutions were developed for numerous types of bearings, using both gas and oil as a lubricant [4].
Regarding the oil film dynamics, Lund [5,6] published a method for linearizing the bearing forces through the calculation of coefficients of stiffness and damping, to be introduced in the equation of motion of the rotating system, which has improved the studies in rotor dynamics. From then on, the study of linear coefficients behavior and their representativeness in the dynamic response of machines spread, becoming the most widely used way for representing the influence of hydrodynamic bearings on the behavior of rotating systems, as can be seen in [7,8,9,10].
Despite the conveniences introduced by using linear coefficients of stiffness and damping for bearings, these bearings can present, in some cases, nonlinear nature that influences the equation of motion for the entire rotating system.
The instability behavior of rotors supported by bearings is widely known as a nonlinear phenomenon and continues to attract the attention of researchers [11,12]. However, nonlinear vibrations can occur in many other, and more realistic, operational conditions.
Therefore, nonlinear rotor models were developed and analyzed by several authors, as can be seen in [13,14,15,16,17,18]. In Machado et al. [17], cylindrical hydrodynamic bearings were subjected to several operating conditions in order to evaluate different nonlinear behaviors. The authors tested the influence of structural damping, the gyroscopic effect, journal eccentricity and external excitation force. Simulation of the journal orbits inside the bearings permitted the observation of differences between the linear and nonlinear models, validated through experimental procedures. It was concluded that the greatest nonlinear impact factor on hydrodynamic forces was external excitation. The results, using linear and nonlinear approaches, diverged with the increase of the unbalance mass, indicating that the linear model was unsatisfactory for system representation in these conditions.
Over the years, researchers have conducted some studies concerning bearings of different geometries, aiming at the inclusion of additional effects, such as analysis of special lubricants, thermal and/or elastohydrodynamic analyzes and misalignment effects, among many others. Some of this research can be referred to in [19,20,21,22,23,24,25,26,27,28,29].
In this context, it is possible to observe that many of the previous research takes into account cylindrical journal bearings to study the behavior of nonlinear vibrations. Additionally, the works that aim to study multi-lobed bearings are concerned in introducing new effects, neglecting their inherent nonlinear behavior. Therefore, this paper expands the studies carried out in [17] on elliptical and trilobed bearings, with its main contribution is to analyze the nonlinear conditions introduced by these kinds of bearings, in the time and frequency domain, in a rotating system subject to severe loading conditions. For this a comparison between the linear and nonlinear approaches is performed. In addition, the vast majority of studies relating to nonlinearities in hydrodynamic bearings are concerned with the phenomenon of fluid induced instability. However, the idea here is to analyze the presence of these nonlinearities in operating conditions of the rotating system below the instability threshold, which are the most common cases found in industrial machines.
The lubricant film is analyzed according to its dynamic characteristics, i.e., evaluating the pressure distribution and, consequently, the resulting hydrodynamic forces. The solution of the Reynolds equation through the finite volume method generates the pressure field inside the oil film, this method being widely used in the numerical simulation of problems for fluid dynamics and heat transfer.
Regarding the severe operating conditions, high load cases are evaluated at low and high rotating speeds, representing situations where the shaft is closer to the bearing surface (low speeds) and situations above the first critical frequency (high speeds). Two different geometries, which are elliptical and three-lobe bearings, are analyzed. In addition, for each geometry, different levels of pre-load are tested.
The analyses are made through evaluation of the temporal response for each system (shaft orbits) and of the spectrum response in the frequency domain (fullspectrum), as stated in [27]. Therefore, the results presented in this paper meet the needs of industry professionals, especially those who are knowledgeable and engaged in the mechanisms of how vibrations arise from rotor-bearing systems as well as in the reduction of these vibrations. Furthermore, the correct understanding of the vibration behavior is mandatory for fault identification in health monitoring of the rotating machinery.

2. Materials and Methods

The theory applied in this paper can essentially be divided into four parts. The first part introduces the basic formulation of the Reynolds equation, which is the foundation of the classical lubrication theory. The second part introduces the linear and nonlinear characteristics of the hydrodynamic forces prevenient from the bearings’ oil film. In part three, the equations of motion, which describe the rotating system, are shown. Finally, the fourth part contains the numerical methods used for solving the equations and obtaining the response of the mechanical system.

2.1. Reynolds Equation

The Reynolds equation used in this work considers a viscous fluid, and its formulation takes into account the momentum conservation and mass conservation equations. Its solution, in case of hydrodynamic bearings, provides the fluid pressure distribution, the information needed to analyze the vast majority of bearing problems. It comes from the Navier–Stokes equations, adopting some simplifying hypotheses for the model [4]. Equation (1) expresses the classical form of the isothermal Reynolds equation for a dynamically loaded journal bearing, where θ is the circumferential coordinate, z is the axial coordinate, μ is the fluid viscosity, h corresponds to the oil film thickness, P is the pressure, U the relative tangential velocity between journal and bearing surfaces and t is time.
θ h 3 P θ + z h 3 P z = 6 μ U h θ + 12 μ h t
The film thickness is shown in Equation (2), where C r is the radial clearance and n is the parameter called “eccentricity ratio” given by n = e / C r . The parameter e is the “eccentricity” and represents the distance between the center of the bearing and the center of the shaft. Figure 1 shows two types of lobed radial bearing geometries, being (Figure 1a) an elliptical bearing and (Figure 1b) a three-lobed bearing. The angle ϕ , defined between the Y direction and the line of centers O b and O j , in the same direction as the θ coordinate, is called “attitude angle”. W is the load referring to rotor weight.
h θ = C r 1 + n cos θ
For the solution of the Reynolds equation, the approach used in the work of Machado and Cavalca [28] is the basis for the numerical model for bearings applied in this paper. This approach utilizes the finite volume method, employing a uniform mesh. Regarding the boundary conditions, atmospheric pressure is considered at the axial edges of the bearing. Near the minimum film thickness location, the Swift-Steiber boundary condition (also known as Reynolds’ condition) is applied along with a relaxation process that makes the negative pressures null after each iteration, assuming no cavitation inside the bearing, as showed in Equation (3).
p = 0                               i f                 p i t e r < 0 p = p i t e r               i f                 p i t e r 0
Also, as previously seen, the oil film thickness h is a fundamental parameter to solve the Reynolds equation, and the whole analysis depends on its clearance. Its calculation differs for each bearing type. For obtaining the film thickness, it is necessary to introduce two new variables, which are called the “eccentricity ratio” and the “equilibrium angle” for each radial bearing lobe.
Two circular arcs from the elliptical bearings, being their centers located on the same line and displaced from the center of the bearing by a distance known as “ellipticity” ε . Defining m = ε / C r as the “ellipticity ratio”, the “eccentricity ratios” for the bottom (1) and top (2) lobes are given by Equation (4). With that, it is also possible to write the “equilibrium angles” for each lobe, as follows in Equation (5).
n 1 = n 2 + m 2 + 2 n m cos ϕ ;         n 2 = n 2 + m 2 2 n m cos ϕ
α 1 = arcsen n sen ϕ n 1 ;         α 2 = arcsen n sen ϕ n 2
Trilobed bearings consist in the intersection of three eccentric lobes, where the center of each lobe has the same distance from the center of the bearing (ellipticity). This configuration allows, in most cases, the formation of three hydrodynamic pressure wedges. The eccentricity ratios for the bottom (b), right (r) and left (l) lobes are given as Equation (6). Equation (7) defines the respective equilibrium angles.
n B = n 2 + m 2 + 2 n m cos ϕ ;       n R = n 2 + m 2 2 n m cos π 3 + ϕ ;
n L = n 2 + m 2 2 n m cos π 3 ϕ
α B = arcsen n sen ϕ n B ;     α R = 2 π 3 arcsen n sen π 3 + ϕ n R ;
α L = 2 π 3 arcsen n sen π 3 ϕ n L
With the eccentricity ratios of each lobe, in both geometries, it is possible to obtain the oil film thickness for each section and, consequently, calculate the pressure distribution of the lobes.

2.2. Hydrodynamic Forces

Once the pressure field is obtained, the oil film resultant forces can be evaluated as shown in Equation (8). This approach calculates nonlinear forces’ components in radial and tangential directions, as can also be seen in Figure 1.
F r F t = 0 2 π L 2 L 2 p cos θ sin θ d z R d θ
Through the attitude angle, it is possible to bring the forces from the reference frame to the inertial one by decomposition, as follows in Equation (9).
F x F y = F r sen ϕ + F t cos ϕ F r cos ϕ + F t sen ϕ
These forces are directly inserted into the system’s equations of motion, when considering the nonlinear model.
For characterizing the oil film behavior, the linear approach considers the lubricant inside the journal bearing as a set of springs and dampers [5,6]. Using a first-order Taylor series expansion, one can write the linearized force as presented in Equation (10), in which the bearing’s dynamic coefficients of stiffness K x x , K y y , K x y and damping C x x , C y y , C x y are the partial derivatives of the forces calculated around the equilibrium position (Equation (11)).
F x = F x 0 + K x x Δ x + K x y Δ y + C x x Δ x ˙ + + C x y Δ y ˙ F y = F y 0 + K y x Δ x + K y y Δ y + C y x Δ x ˙ + + C y y Δ y ˙
K x y = F x y 0 ;     C x y = F x y ˙ 0
F x 0 and F y 0 are the hydrodynamic forces at the journal equilibrium position and ( Δ x , Δ y ) and ( Δ x ˙ , Δ y ˙ ) disturbances around it.

2.3. Equation of Motion

Each element of the rotating system is individually modelled by the finite element method, allowing inertial coupling between translation and rotation through the mass matrix [M], stiffness matrix [K], gyroscopic matrix [G] and proportional damping matrix [C] (as showed in Equation (12)), where α and β are proportional coefficients.
C = α M + β K
The Lagrange equation is used to obtain the system’s equations of motion, being its generic form shown in Equation (13). Moreover, qi is the ith generalized coordinate and Fi the generalized force component acting in the direction of each generalized coordinate.
d d t T q ˙ i T q i + U q i + R q ˙ i = F i ;           i = 1 , 2 , , N
Defining the energies—kinetic T , deformation U and work performed by the nonconservative forces R —for each element i, the total energy of the system is given by the sum of the energies of each element, as follows in Equation (14).
T = i = 1 N T i ;           U = i = 1 N U i ;               R = i = 1 N R i
According to Lallane and Ferraris [26] the kinetic energy, deformation and the work of the global nonconservative forces are given as Equation (15), in which Ω is the shaft rotating speed and I P the polar moment of inertia.
T = 1 2 q ˙ T M q ˙ Ω q ˙ T G q + 1 2 I P Ω 2 ; U = 1 2 q T K q ; R = 1 2 q ˙ T C q ˙
Combining Equation (13) and Equation (15), the system’s complete equation of motion follows as in Equation (16). The external forces vector (Equation (17)) is composed by the sum of the exciting forces f e t (usually due to the unbalance mass force), the rotor’s weight force W and the hydrodynamic forces f h t (which can be modeled as linear (Equation (10)) or nonlinear (Equation (9))).
M q ¨ t + C + Ω G q ˙ t + K q t = f t
f t = f e t + W + f h t

2.4. Solution Methods

As said before, the bearing model takes into account the solution of the Reynolds equation (solved here through finite volume numerical method). Regarding the rotating system modeling, the finite element method is used to obtain the equation of motion, being the main elements employed: rigid disk element and Timoshenko beam element.
Considering the linear analysis, after calculating bearing’s support reactions and unbalance forces, the Reynolds equation must be statically solved to obtain the equilibrium position. Then, finite perturbations are applied to calculate the bearing dynamic coefficients, as shown in Equation (10). These dynamic coefficients are inserted into the bearing nodes aiming for the computation of the global force vector. If the analysis is nonlinear, the hydrodynamic forces (Equation (8)) are obtained directly from the dynamic Reynolds equation solution.
Thus, assembling all global matrices as well as the forces vector, the solution of the rotor dynamic problem is provided by integrating the equation of motion. The use of nonlinear hydrodynamic forces can cause severe changes to the machines’ dynamic behavior. So, a robust numeric integration scheme should be employed. A widely used method for solving time domain differential equations, especially when these equations have a certain degree of nonlinearity, is the Newmark method.
For each time step, the prediction and correction of the variables (displacement, velocity and acceleration) are calculated. Then, using Newton–Raphson iterative method, the variables are updated according to the equation of motion. Finally, the time step is incremented, repeating the procedure until the entire time interval is contemplated [29].

3. Results and Discussion

In this work, the computational simulations are presented, and discussed, essentially in two parts, namely, the numerical results of orbits, and the spectrum of the system response (by Discrete Fourier Transform—full spectrum). Two types of radial lobed bearings were used in the system for comparation, which are elliptical and three-lobed bearings. The preload parameter is varied, being 0.25 and 0.8, in order to simulate higher and lower eccentricity conditions, respectively. Thus, combining all cases, it was possible to verify the influence of four different bearing configurations on the simulated system response. Both three-lobe (Equation (6)) and elliptical (Equation (4)) bearing models were considered to calculate the hydrodynamic forces. The orbits were plotted simultaneously and compared for each case, facilitating the discussion.
The rotating system used is shown in Figure 2. It consists of 18 beam elements (19 nodes) and 1 disk element. The SAE 1020 steel shaft has elastic modulus of 2 × 1011 N/m2, poison coefficient of 0.3, 7800 kg/m3 density, 700 mm length, 15 mm diameter. Moreover, it is supported by two hydrodynamic bearings, placed approximately at 580 mm from each other, positioned at nodes 3 and 17, respectively, having a length of 20 mm, an inner diameter of 31 mm with 90 μm of radial clearance. The lubricant is the ISO VG 32 oil at a constant average temperature of 40 °C ( μ = 2.7622 × 10−2 Pa·s). The disk was placed at node 7, positioned 200 mm from bearing 1, with 94.7 mm diameter and 47 mm length. In this case, the unbalanced mass was m 0 = 3 g, located 37 mm away from the center of the disk. In addition to the unbalance, a 200 N external force was applied to the disk node to simulate a severe loading condition. As shown in Equation (12), the structural damping can be proportional to the mass and stiffness matrices, but for steel shafts, the recommendation is to only use a matrix proportional to the stiffness. Therefore, the coefficients used here were α = 0   and   β = 1.5 × 10 5 (this value was experimentally estimated in [17]).
Is worth noting that the procedures used for the simulations—the finite element modelling and time integration scheme, as well as the solution of the Reynolds equation—have been experimentally validated in [17] for a rotor system supported by cylindrical journal bearings. Although the bearing geometries presented are different from the cylindrical one, the procedure used to solve the Reynolds equation is the same of [17].
Figure 3 shows the pressure field distribution in each bearing (elliptical and trilobed) in order to exemplify its behavior in each analyzed geometry. Preload 0.8 was used, since the pressure peaks can be better visualized in this situation. When considering the elliptical bearing, it is possible to observe the presence of two pressure peaks due to the two lobes. Similarly, for the three-lobe bearing, three pressure peaks are seen due to the three lobes.
Figure 4 depicts the behavior of the stiffness (Figure 4a) and damping (Figure 4b) coefficients for the elliptical and three-lobe bearings, which are a function of the rotating speed. These coefficients are used in the linear bearing approach. Analyzing the dynamic coefficients, it is possible to note that the cross-coupled coefficients (xy, yx) for the three-lobe case are smaller than for the elliptical case, confirming that three-lobe bearings present greater stability in comparison to the elliptical one, as expected.

3.1. Numerical Results of Orbits

The numerical analyses made here compare the orbits of shaft inside the bearings. Since the rotor is asymmetric to the presence of a decentralized disk, as seen in Figure 2, bearing 1 was chosen for viewing the orbits, since gyroscopic effects are more evident there. The first critical speed of the system is about 37 Hz. Thus, five rotating speeds were selected to evaluate the orbit behavior in conditions of high eccentricity (below first critical speed), low eccentricity (above first critical speed) and high amplitudes (around first critical speed) being 15, 30, 35, 45 and 60 Hz. Additionally, the Sommerfeld number range goes from 0.049 to 0.2.
Figure 5 shows the orbits results for the elliptical bearing with preload of 0.25. A total simulated time of 20 s is used, in which only the last second of the converged orbit is plotted. It is possible to visualize the “self-centering” effect, in which with increasing rotating speed the shaft tends to be centered in the radial clearance of the bearing. Additionally, the major discrepancies between linear and nonlinear models occur near the first critical speed of the system (35 Hz, green lines), where the vibration amplitude is higher. In this case, either shape or size can be different when considering the nonlinear model. For the other rotating frequencies, just slight differences can be observed between the models, maintaining similar orbit shapes and positions.
By increasing the preload level to 0.8, the shapes of the orbits remained similar, except for 35 Hz, which has differences due to the proximity with the critical speed (Figure 6). However, the orbits are smaller than the previous case, since the higher the preload, the smaller the radial clearance in the vertical direction, which creates a stiffer oil film and, consequently, reduces the vibration amplitudes. Additionally, it is interesting to note that Figure 6 shows a zoom in the whole region containing the orbits.
Changing the bearing type, Figure 7 shows orbits for the three-lobed bearing with preload 0.25. In this case, differences in size can be observed at all rotating speeds, indicating a greater influence of nonlinearities in this bearing geometry. At 35 Hz, near the first critical speed, the discrepancies in size are ever higher, and it is possible to observe differences in shape and position as well.
For the three-lobed bearing with preload 0.8, the central position of the orbits remained the same for all rotating speeds (comparing linear and nonlinear approaches), as shown in Figure 8. However, the orbits’ size and shape change at all rotations in an evident way. In the proximity of the critical region (35 Hz), the amplitude of the orbit increases, but as stated, maintaining the same position. This shows that the nonlinearity substantially affects three-lobed bearings and has an even greater effect in the case of high preload.
According to results obtained by the calculation of the orbits, the conclusion is that for elliptical bearings the nonlinearities are only influential in cases of higher orbit amplitudes, i.e., in the vicinity of the critical speed. This behavior is similar to that observed in cylindrical bearings, as shown in the work of reference [17]. On the other hand, in the case of three-lobed bearings, the influence of nonlinearities is much more evident, being practically observed in all rotating speeds, with greater sensitivity the higher the value of the bearing preload.

3.2. Numerical Results of DFT Full Spectrum

Spectral analysis is an alternative way of identifying and analyzing signals, which becomes complementary to time domain analysis. For this, the signal is decomposed in the frequency domain, allowing the identification of the components of that signal. A known tool for the treatment of discrete and aperiodic signals is the DFT (“Discrete Fourier Transform”). A conventional DFT is called a “half spectrum”, as it relates only the positive part of the frequency domain, commonly in harmonics, e.g., 1×, 2×, 3× etc. For cases in which information about the relative phase is important, such as in the diagnosis of unbalance and misalignment faults, among others, the application of the “DFT full spectrum” is indicated [27]. In these situations, the full spectrum enables the determination of the rotor’s behavior under “forward” and “backward” precession conditions.
Each orbit has, in general, an elliptical shape, which can be represented as the sum of two circular orbits, one in forward rotation and the other in backward rotation, both rotating at the same frequency, which is the frequency of the filtered orbit. Thus, the orbit will be more or less circular, depending on the amplitude of the signal components in the frequency domain.
Therefore, Figure 9 shows the first results for the elliptical bearing 1, with preload 0.25, using the DFT full spectrum decomposition. To facilitate the interpretation of the graphs, Table 1 brings a comparison of the magnitudes of the most relevant amplitude peaks for all analyzed cases. The intention is to compare the influence of linear and nonlinear models on the frequency spectrum for the different bearing types.
It is possible to observe the presence of the −1× component in all rotations, for both models, on the left-hand side of the spectrum, i.e., frequency related to backward mode, due to the fact that hydrodynamic bearings inherently introduce some anisotropy into the system. In linear cases, there are only the 1× and −1× components, indicating the purely elliptical shape of the orbit, as can be attested in Figure 5. For nonlinear cases, there are significant values related to high order harmonics, such as 2× and −2×. For 35 Hz, the components −3× and 3× appear, indicating the proximity of the system’s critical speed region, which is an indication that the system is more nonlinear under these conditions. These additional harmonics influences the orbit size and shape, which can again be clearly observed in Figure 5.
By increasing the preload to 0.8, for the elliptical bearing 1, it is possible to observe similar behavior to the case of smaller preload, as shown in Figure 10 and Table 2. The amplitudes are at least 1 order of magnitude lower, and there is no presence of −3× and 3× components for the rotation of 35 Hz, due to the smaller orbit originated by the low journal eccentricity.
For the three-lobed bearing 1, with preload 0.25, it is possible to see in Figure 11 and Table 3 that amplitudes 1× and −1× are greater for the nonlinear case (mainly the forward component 1×). Higher vibration orbits are more prone to nonlinearities coming from the hydrodynamic bearings. So, it is possible to state that the three-lobed bearing is more susceptible to present nonlinear behavior, as can also be seen by the presence of additional harmonics in the signal. Comparing Table 1 and Table 3, it is perceivable that these differences are higher for the three-lobed bearing. Regarding the other harmonic components for all rotations, they have similar behavior to those observed for the elliptical bearing cases.
Finally, Figure 12 and Table 4 show the results for the three-lobed bearing 1 with preload 0.8. This case has somewhat different behavior. The only frequency components observed are 1× and −1×, indicating a purely elliptical shape of the orbit for both linear and nonlinear models, even though the orbits obtained with the two models are quite different, as shown in Figure 8. By observing the rotating speeds and comparing the models, it can be seen that the 1× amplitudes have very close values, with just slight differences. However, the nonlinear −1× amplitudes have higher values when compared to the linear ones, indicating the deformations in the orbits (Figure 8), which means more anisotropy in the system, since in this case, the backward component is greater.
Analyzing the decomposition of the signals by DFTs, it is possible to conclude that, for low preloaded bearings, for both elliptical and three-lobed cases, the nonlinearities are more detectable through the appearance of higher order harmonics, such as 2×, 3× (forward) and −2×, −3× (backward). On the other hand, for the case of three-lobed bearings with high preload, the nonlinearities are evidenced by the increased growth in the backward component −1×.

4. Conclusions

The paper presents comparative results of linear and nonlinear models for elliptical and three-lobe hydrodynamic bearings, operating under severe loading conditions. The analysis included the comparison of the orbits (by the numerical integration of the system in time-domain) and signal amplitudes in the frequency domain (by DFT full spectrum).
Regarding the presented orbits, the results showed that for elliptical bearings, the situations with the greatest presence of nonlinearities occur with high vibration amplitudes, near the critical speed region, whereas for three-lobed bearings, the nonlinearities have greater influence even in smaller vibration amplitude situations.
An important point was seen in the analysis of the full spectrum, which showed that nonlinearities are not always evidenced by the appearance of higher order harmonic components but can also be presented through a significant increase in the −1× component.
Finally, the combined analyses presented here are of great importance for mapping the influence of nonlinearities in different types of multi-lobed hydrodynamic bearings with fixed geometry. As the number of lobes, and/or the preload value, increases, the vibration levels decrease, generally pointing to a more linear bearing. Consequently, the traditional linear analysis could be employed without major drawbacks. This study is important for monitoring and fault diagnosis, mainly those that are model based since representative mathematical models are mandatory for a reliable diagnosis.

Author Contributions

Conceptualization, C.A.A.V., D.S.A. and T.H.M.; methodology, C.A.A.V., D.S.A. and T.H.M.; software, C.A.A.V. and D.S.A.; validation, C.A.A.V., D.S.A. and T.H.M.; formal analysis, C.A.A.V., D.S.A. and T.H.M.; investigation, C.A.A.V., D.S.A. and T.H.M.; resources, T.H.M.; data curation C.A.A.V., D.S.A. and T.H.M.; writing—original draft preparation, C.A.A.V.; writing—review and editing, D.S.A. and T.H.M.; visualization, C.A.A.V. and D.S.A.; supervision, D.S.A. and T.H.M.; project administration, D.S.A. and T.H.M.; funding acquisition, C.A.A.V., T.H.M. All authors have read and agreed to the published version of the manuscript.

Funding

Coordination of Superior Level Staff Improvement—CAPES, National Council for Scientific and Technological Development—CNPq, grant # 424899/2018-3, and grants # 2017/07454-8 and # 2018/21581-5 from the São Paulo Research Foundation (FAPESP).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pinkus, O. Analysis of Elliptical Bearings. Trans. ASME 1956, 78, 965–973. [Google Scholar]
  2. Pinkus, O. Solution of Reynolds Equation for Finite Journal Bearings. Trans. ASME 1958, 80, 858–864. [Google Scholar]
  3. Pinkus, O. Analysis and Characteristics of the Three-Lobe Bearing. J. Basic Eng. 1959, 81, 49–55. [Google Scholar] [CrossRef]
  4. Pinkus, O.; Sternlicht, S.A. Theory of Hydrodynamic Lubrication; McGraw-Hill: New York, NY, USA, 1961. [Google Scholar]
  5. Lund, J.W.; Thomsen, K.K. A Calculation Method and Data for the Dynamic Coefficients of Oil-Lubricated Journal Bearings. Topics in Fluid Film Bearing and Rotor Bearing System Design and Optimization 1978. pp. 11–28. Available online: https://dyrobes.com/wp-content/uploads/2015/09/A-Calculation-Method_Oil-Lubricated_ASME-1982_linked.pdf (accessed on 22 November 2021).
  6. Lund, J.W. Review of the Concept of Dynamic Coefficients for Fluid Film Journal Bearings. J. Tribol. 1987, 109, 37–41. [Google Scholar] [CrossRef]
  7. Qiu, Z.L.; Tieu, A.K. The Effect of Perturbation Amplitudes on Eight Force Coefficients of Journal Bearings. Tribol. Trans. 1996, 39, 469–475. [Google Scholar] [CrossRef]
  8. Müller-Karger, C.M.; Granados, A.L. Derivation of Hydrodynamic Bearing Coefficients Using the Minimum Square Method. J. Tribol. 1997, 119, 802–807. [Google Scholar] [CrossRef]
  9. Zhao, S.; Zhou, H.; Meng, G.; Zhu, J. Experimental identification of linear oil-film coefficients using least-mean-square method in time domain. J. Sound Vib. 2005, 287, 809–825. [Google Scholar] [CrossRef]
  10. Li, K.; Liu, J.; Han, X.; Jiang, C.; Qin, H. Identification of oil-film coefficients for a rotor-journal bearing system based on equivalent load reconstruction. Tribol. Int. 2016, 104, 285–293. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, S.; Palazzolo, A.B. Effects of thermo hydrodynamic (THD) floating ring bearing model on rotordynamic bifurcation. Int. J. Non-Linear Mech. 2017, 95, 30–41. [Google Scholar] [CrossRef]
  12. Tehrani, G.G.; Gastaldi, C.; Berruti, T.M. Stability analysis of a parametrically excited ball bearing system. Int. J. Non-Linear Mech. 2020, 120, 103350. [Google Scholar] [CrossRef]
  13. Hattori, H. Dynamic Analysis of a Rotor-Journal Bearing System with Large Dynamic Loads: Stiffness and Damping Coefficient Variations in Bearing Oil Films. JSME Int. J. Ser. C Dyn. Control Robot. Des. Manuf. 1993, 36, 251–257. [Google Scholar] [CrossRef]
  14. Hu, A.; Hou, L.; Xiang, L. Dynamic Simulation and Experimental Study of an Asymmetric Double-Disk Ro-tor-Bearing System with Rub-Impact and Oil-Film Instability. Nonlinear Dyn. 2016, 84, 641–659. [Google Scholar] [CrossRef]
  15. Dakel, M.; Baguet, S.; Dufour, R. Dynamic analysis of a harmonically excited on-board rotor-bearing system. In Proceedings of the 10th International Conference on Vibrations in Rotating Machinery, London, UK, 11–13 September 2012; Elsevier: Amsterdam, The Netherlands, 2012; pp. 57–67. [Google Scholar]
  16. Smolík, L.; Rendl, J.; Dyk, Š.; Polach, P.; Hajžman, M. Threshold stability curves for a nonlinear rotor-bearing system. J. Sound Vib. 2019, 442, 698–713. [Google Scholar] [CrossRef]
  17. Machado, T.H.; Alves, D.S.; Cavalca, K.L. Discussion about nonlinear boundaries for hydrodynamic forces in journal bearing. Nonlinear Dyn. 2018, 92, 2005–2022. [Google Scholar] [CrossRef]
  18. Li, C.; She, H.; Tang, Q.; Wen, B. The effect of blade vibration on the nonlinear characteristics of rotor–Bearing system supported by nonlinear suspension. Nonlinear Dyn. 2017, 89, 987–1010. [Google Scholar] [CrossRef]
  19. Arumugam, P.; Swarnamani, S.; Prabhu, B.S. Effects of Journal Misalignment on the Performance Characteristics of Three-Lobe Bearings. Wear 1997, 206, 122–129. [Google Scholar] [CrossRef]
  20. Nair, K.P.; Nair, V.P.S.; Jayadas, N.H. Static and Dynamic Analysis of Elastohydrody-namic Elliptical Journal Bearing with Micropolar Lubricant. Tribol. Int. 2007, 40, 297–305. [Google Scholar] [CrossRef]
  21. Müller-Karger, C.M.; Barrett, L.E.; Flack, R.D. Influence of Fluid Film Nonlinearity on the Experimental Determination of Dynamic Stiffness and Damping Coefficients for Three-Lobe Journal Bearings. Tribol. Trans. 1997, 40, 49–56. [Google Scholar] [CrossRef]
  22. Rahmatabadi, A.; Nekoeimehr, M.; Rashidi, R. Micropolar lubricant effects on the performance of noncircular lobed bearings. Tribol. Int. 2010, 43, 404–413. [Google Scholar] [CrossRef]
  23. Chauhan, A.; Sehgal, R.; Sharma, R.K. Thermohydrodynamic Analysis of Elliptical Journal Bearing with Different Grade Oils. Tribol. Int. 2010, 43, 1970–1977. [Google Scholar] [CrossRef]
  24. Zhang, S.; Xing, Y.; Xu, H.; Pei, S.; Zhang, L. An experimental study on vibration suppression of adjustable elliptical journal bearing-rotor system in various vibration states. Mech. Syst. Signal Process. 2020, 141, 106477. [Google Scholar] [CrossRef]
  25. Zhang, S.; Xu, H.; Zhang, L.; Xing, Y.; Guo, Y. Vibration suppression mechanism research of adjustable elliptical journal bearing under syn-chronous unbalance load. Tribol. Int. 2019, 132, 185–198. [Google Scholar] [CrossRef]
  26. Lalanne, M.; Ferraris, G. Rotordynamics Prediction in Engineering, 2nd ed.; John Wiley & Sons: Chichester, UK, 1998. [Google Scholar]
  27. Goldman, P.; Muszynska, A. Application of full spectrum to rotating machinery diagnostics. Orbit 1999, 20, 17–21. [Google Scholar]
  28. Machado, T.H.; Cavalca, K.L. Modeling of hydrodynamic bearing wear in rotor-bearing systems. Mech. Res. Commun. 2015, 69, 15–23. [Google Scholar] [CrossRef]
  29. Bathe, K.J. Finite Element Procedures in Engineering Analysis; Prentice-Hall: Hoboken, NJ, USA, 1982. [Google Scholar]
Figure 1. Basic journal bearing geometry: (a) elliptical bearing; (b) three-lobed bearing.
Figure 1. Basic journal bearing geometry: (a) elliptical bearing; (b) three-lobed bearing.
Applsci 12 03215 g001
Figure 2. Sketch of the rotor used in the simulations (out of scale, for better visualization).
Figure 2. Sketch of the rotor used in the simulations (out of scale, for better visualization).
Applsci 12 03215 g002
Figure 3. Pressure field distribution—(a) Elliptical bearing; (b) Three-lobe bearing.
Figure 3. Pressure field distribution—(a) Elliptical bearing; (b) Three-lobe bearing.
Applsci 12 03215 g003
Figure 4. Dynamic Coefficients for elliptical and three-lobe bearings—(a) Stiffness; (b) Damping.
Figure 4. Dynamic Coefficients for elliptical and three-lobe bearings—(a) Stiffness; (b) Damping.
Applsci 12 03215 g004
Figure 5. Orbits for all rotations simulated in elliptical bearing 1 (preload 0.25), with zoom in rotation 15 Hz.
Figure 5. Orbits for all rotations simulated in elliptical bearing 1 (preload 0.25), with zoom in rotation 15 Hz.
Applsci 12 03215 g005
Figure 6. Orbits for all rotations simulated in elliptical bearing 1 (preload 0.8), with zoom in rotation 15 Hz.
Figure 6. Orbits for all rotations simulated in elliptical bearing 1 (preload 0.8), with zoom in rotation 15 Hz.
Applsci 12 03215 g006
Figure 7. Orbits for all rotations simulated in three-lobed bearing 1 (preload 0.25), with zoom in rotation 15 Hz.
Figure 7. Orbits for all rotations simulated in three-lobed bearing 1 (preload 0.25), with zoom in rotation 15 Hz.
Applsci 12 03215 g007
Figure 8. Orbits for all rotations simulated in three-lobed bearing 1 (preload 0.8), with zoom in rotation 15 Hz.
Figure 8. Orbits for all rotations simulated in three-lobed bearing 1 (preload 0.8), with zoom in rotation 15 Hz.
Applsci 12 03215 g008
Figure 9. DFT full spectrum of simulated results in elliptical bearing 1 (preload 0.25): (a) linear, (b) nonlinear.
Figure 9. DFT full spectrum of simulated results in elliptical bearing 1 (preload 0.25): (a) linear, (b) nonlinear.
Applsci 12 03215 g009
Figure 10. DFT full spectrum of simulated results in elliptical bearing 1 (preload 0.8): (a) linear, (b) nonlinear.
Figure 10. DFT full spectrum of simulated results in elliptical bearing 1 (preload 0.8): (a) linear, (b) nonlinear.
Applsci 12 03215 g010
Figure 11. DFT full spectrum of simulated results in three-lobed bearing 1 (preload 0.25): (a) linear, (b) nonlinear.
Figure 11. DFT full spectrum of simulated results in three-lobed bearing 1 (preload 0.25): (a) linear, (b) nonlinear.
Applsci 12 03215 g011
Figure 12. DFT full spectrum of simulated results in three-lobed bearing 1 (preload 0.8): (a) linear, (b) nonlinear.
Figure 12. DFT full spectrum of simulated results in three-lobed bearing 1 (preload 0.8): (a) linear, (b) nonlinear.
Applsci 12 03215 g012
Table 1. Magnitudes of the most relevant peaks of Figure 9—elliptical bearing (preload 0.25).
Table 1. Magnitudes of the most relevant peaks of Figure 9—elliptical bearing (preload 0.25).
Rot. SpeedAmplitudes (m)
LinearNonlinear
−3×−2×−1×−3×−2×−1×
15 Hz--1.27 × 10−72.04 × 10−7----1.46 × 10−72.36 × 10−7--
30 Hz--1.31 × 10−62.44 × 10−6---6.94 × 10−81.39 × 10−62.62 × 10−67.41 × 10−8-
35 Hz--6.67 × 10−61.27 × 10−5--2.30 × 10−71.11 × 10−65.32 × 10−61.11 × 10−51.24 × 10−61.78 × 10−7
45 Hz--1.89 × 10−63.97 × 10−6---1.26 × 10−71.85 × 10−63.93 × 10−61.26 × 10−7-
60 Hz--6.07 × 10−71.39 × 10−6---2.06 × 10−87.19 × 10−71.65 × 10−62.30 × 10−8-
Table 2. Magnitudes of the most relevant peaks of Figure 10—elliptical bearing (preload 0.8).
Table 2. Magnitudes of the most relevant peaks of Figure 10—elliptical bearing (preload 0.8).
Rot. SpeedAmplitudes (m)
LinearNonlinear
−3×−2×−1×−3×−2×−1×
15 Hz--5.42 × 10−89.07 × 10−8----5.46 × 10−89.10 × 10−8--
30 Hz--4.42 × 10−77.38 × 10−7---1.00 × 10−83.93 × 10−76.56 × 10−71.34 × 10−8-
35 Hz--1.77 × 10−62.94 × 10−6---7.58 × 10−81.85 × 10−63.08 × 10−61.25 × 10−7-
45 Hz--6.28 × 10−71.05 × 10−6---7.08 × 10−96.32 × 10−71.06 × 10−61.14 × 10−8-
60 Hz--1.74 × 10−72.89 × 10−7----2.53 × 10−74.19 × 10−7--
Table 3. Magnitudes of the most relevant peaks of Figure 11—three-lobed bearing (preload 0.25).
Table 3. Magnitudes of the most relevant peaks of Figure 11—three-lobed bearing (preload 0.25).
Rot. SpeedAmplitudes (m)
LinearNonlinear
−3×−2×−1×−3×−2×−1×
15 Hz--1.10 × 10−71.74 × 10−7----1.42 × 10−72.17 × 10−7--
30 Hz--1.07 × 10−61.96 × 10−6---4.12 × 10−81.44 × 10−62.55 × 10−63.97 × 10−8-
35 Hz--4.94 × 10−69.31 × 10−6--2.43 × 10−78.27 × 10−76.97 × 10−64.32 × 10−55.23 × 10−72.54 × 10−7
45 Hz--1.64 × 10−63.36 × 10−6---8.38 × 10−81.87 × 10−63.80 × 10−67.54 × 10−8-
60 Hz--5.83 × 10−71.31 × 10−6----7.84 × 10−71.73 × 10−6--
Table 4. Magnitudes of the most relevant peaks of Figure 12—three-lobed bearing (preload 0.8).
Table 4. Magnitudes of the most relevant peaks of Figure 12—three-lobed bearing (preload 0.8).
Rot. SpeedAmplitudes (m)
LinearNonlinear
−3×−2×−1×−3×−2×−1×
15 Hz--9.82 × 10−93.08 × 10−8----1.96 × 10−83.84 × 10−8--
30 Hz--3.73 × 10−81.95 × 10−7----9.50 × 10−82.24 × 10−7--
35 Hz--1.20 × 10−77.04 × 10−7----3.29 × 10−78.06 × 10−7--
45 Hz--3.57 × 10−82.71 × 10−7----9.36 × 10−82.52 × 10−7--
60 Hz--6.83 × 10−96.99 × 10−8----3.56 × 10−81.02 × 10−7--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Viana, C.A.A.; Alves, D.S.; Machado, T.H. Linear and Nonlinear Performance Analysis of Hydrodynamic Journal Bearings with Different Geometries. Appl. Sci. 2022, 12, 3215. https://doi.org/10.3390/app12073215

AMA Style

Viana CAA, Alves DS, Machado TH. Linear and Nonlinear Performance Analysis of Hydrodynamic Journal Bearings with Different Geometries. Applied Sciences. 2022; 12(7):3215. https://doi.org/10.3390/app12073215

Chicago/Turabian Style

Viana, Carlos Alberto Alves, Diogo Stuani Alves, and Tiago Henrique Machado. 2022. "Linear and Nonlinear Performance Analysis of Hydrodynamic Journal Bearings with Different Geometries" Applied Sciences 12, no. 7: 3215. https://doi.org/10.3390/app12073215

APA Style

Viana, C. A. A., Alves, D. S., & Machado, T. H. (2022). Linear and Nonlinear Performance Analysis of Hydrodynamic Journal Bearings with Different Geometries. Applied Sciences, 12(7), 3215. https://doi.org/10.3390/app12073215

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop