# Influence of the Radial Clearance of a Squeeze Film Damper on the Vibratory Behavior of a Single Spool Gas Turbine Designed for Unmanned Aerial Vehicle Applications

Authors: Geraldo Creci; José Orlando Balastrero; Sidney Domingues; Luis Vanderlei Torres; João Carlos Menezes
Journal: Shock and Vibration
 (2017)

---

## Abstract

This study presents a numerical investigation using the finite element method on the vibratory behavior of a single spool gas turbine designed for unmanned aerial vehicle applications. The shaft of the rotor-bearing system is supported on a front bearing composed of a deep groove ball bearing with a vibration absorber element and a rear squeeze film damper bearing. Three radial clearances for the squeeze film damper were analyzed to determine the best geometric configuration for the rear bearing, considering the rotordynamic performance of the entire system. Whirl speeds and unbalanced system responses were carefully evaluated to determine the best radial clearance for the squeeze film damper. After defining the best radial clearance, a transient analysis was performed to simulate the transition of the system through resonance, and a spectral map is presented to illustrate the vibratory behavior of the system considering the influence of all related important frequencies. The rotordynamic behavior of the system is predicted, and vibration problems are avoided. Its mechanical drawings were released to manufacturing, and the first prototype is in the experimental test phase, thus indicating that the numerical results presented in this study are consistent.


---

## Body

## 1. Introduction

The demand for more powerful rotating machines has led to higher operating speeds, which has resulted in the need for an accurate prediction of the dynamic behavior of rotating machines in the design and development stages to avoid undesired vibration conditions at operating speeds. Understanding the physics of rotating machines is essential in several application fields [1–4]. For example, the safety and service life of rotating machinery in aerospace, automotive, and power plants are directly linked to the high quality of their vibratory behavior. A high level of reliability, increased operation, and decreased unscheduled maintenance is becoming increasingly desirable for the currently developed industry products [5].The rotors in rotating machines are always slightly unbalanced due to manufacturing and assembling inaccuracies. In the rotating components of turbine engines, particularly bladed discs, one of the primary sources of structural uncertainty is mistuning. Mistuning refers to variations in the characteristics of each sector of a system that is ideally expected to be cyclic symmetric [6–8]. It generally results in a slight variability in the dynamical responses and induces localization phenomena and stress magnification. In particular, if the critical speeds are exceeded, then the unbalance produces a vibration of large amplitude.A full rotordynamic analysis is mandatory to avoid vibration resonance at operating speeds. Critical speeds are defined as coincidences in the shaft rotating speed and the rotating natural frequencies of the rotor-bearing system. Because catastrophic failure occurs at critical speeds due to resonance, critical speeds are designed to be sufficiently separate from the operating speed range. Furthermore, determining both the instability threshold and unbalanced system responses is of significant concern [9–11].There are several numerical approximations for the vibration analysis of rotor-bearing systems. The most popular approach, which is particularly well-suited for modeling large scale and complicated systems, is the finite element method [12–17]. In recent years, several articles have been published addressing full dynamic analyses in numerous types of rotating machineries and have involved complete investigations of the rotordynamic behavior in the respective application and avoid vibration problems [18–30].In this study, a full rotordynamic analysis is performed for a single spool gas turbine designed for unmanned aerial vehicle applications. The referred gas turbine has been developed as an academic research project in association with industry. It is important to observe that, in a gas turbine, excitation frequencies can result from a variety of sources. The most critical source is the mass unbalance of the rotor at a frequency Ω. Coupling misalignments result in an excitation frequency of 2Ω. Blades, vanes, nozzle, diffuser, and other devices produce excitation frequencies ofsΩ, where s is the number of blades, vanes, or other devices. Bearing excitation sources can produce subharmonic excitation frequencies of 0.5Ω. In the current application, there is no coupling. Therefore, 2Ω excitation frequency can be discharged. Additionally,sΩ excitation frequencies can only exist with very low energies and at very high frequencies, being also discharged. Subharmonic excitation from the bearings is not likely to occur. However, possible vibration amplitudes at 0.5Ω frequency are considered.The rotordynamic behavior of the developed single spool gas turbine is investigated by considering the dynamic stiffness and damping coefficients of the bearings. The expected performance of the engine is to generate 5 kN of thrust at 28,150 rpm under international standard atmosphere conditions. The shaft of the rotor-bearing system is supported on a front bearing composed by a deep groove ball bearing with a vibration absorber element [31, 32] and a rear squeeze film damper bearing [33, 34]. Both bearings were designed to provide the appropriate stiffness and damping characteristics based on the entire rotor-bearing system. Three different radial clearances for the squeeze film damper were analyzed to determine the best configuration to avoid potential vibration problems.


## 2. Single Spool Gas Turbine

### 2.1. Technical Characteristics

The aeronautical gas turbine developed in this study has a single spool and is designed to be used in unmanned aerial vehicle applications. It has an air inlet duct, five-stage axial compressor, flow-through annular combustion chamber, single-stage axial turbine, and exhaustion nozzle. Figures1(a) and 1(b) illustrate the manufactured gas turbine currently in the experimental test phase.Figure 1
Manufactured single spool gas turbine developed in this study: (a) apparatus for demonstrative purposes and preliminary tests and (b) test rig for full test protocol.
(a)
(b)It should be noted that low weight and high thrust generation are highly desirable characteristics in smaller modern aircraft gas turbines. Based on user requirements, the dry weight of the developed gas turbine should be approximately 650 N, its maximum length should not exceed 1.50 m, and the maximum diameter of the circular section should not exceed 0.35 m. The operation is expected for a flight envelope from 0 to 1,200 m, maximum speeds of up to 0.9 Mach, acceleration of 10 g in the longitudinal axis and 7 g in the transverse axis, and cumulative operating time of at least 30 h. By having a single spool, the design, manufacture, and assembly of the rotor-bearing system occur more robustly, simply, and within a shorter amount of time. The operating speed range of the gas turbine is set between 80 and 100% of the nominal speed based on the operation at the design point. Therefore, its operating speed range is set between 375.33 and 469.17 Hz, which are the limits of the rotation speeds in which the gas turbine will mostly operate.Figure2 depicts a three-dimensional illustration of the cross-section of the developed gas turbine with details on the two bearings that support the rotor-bearing system. In Figure 2(a), the primary internal components of the aircraft gas turbine and its external general shape can be seen in detail, with a primary focus on the bipartite mounting of the axial compressor. The rotor-bearing system is supported by two bearings: a front or cold bearing and a rear or hot bearing. The front bearing consists of a deep groove ball bearing and a vibration absorber element. The rear bearing is an open-ended unsealed squeeze film damper with a circumferential feeding groove. The lubrication system has been carefully designed to meet the requirements of the bearings, which are continuously lubricated by aeronautical oil with specific performance characteristics.Figure 2
Three-dimensional illustration of the cross-section of the developed gas turbine: (a) bipartite mounting of the axial compressor and other inner components; (b) front or cold bearing; and (c) rear or hot bearing.
(a)
(b)
(c)


### 2.2. Front Bearing-Deep Groove Ball Bearing

As stated previously, the front bearing of the single spool gas turbine consists of a vibration absorber element and a deep groove ball bearing. Furthermore, the front bearing is constituted by lubricating channels to transport the oil from the lubrication system. Figure2(b) illustrates the configuration of the front bearing in detail. The deep groove ball bearing is continuously lubricated by aeronautic oil that can endure high temperature loads. In addition to lubrication, the oil dampens vibrations. The vibration absorber element provides a designed stiffness to the front bearing, which contributes to more acceptable vibration amplitudes. The stiffness of the vibration absorber element is lower than the stiffness of the deep groove ball bearing contacts and approximately constant for all rotational speeds of the rotor-bearing system. Thus, the characterization of the front bearing should be performed considering the influence of the stiffness of the vibration absorber element and the dynamic damping properties of the deep groove ball bearing contacts [9]. In [35], the procedures for characterization of the stiffness and damping dynamic coefficients of the deep groove ball bearing are presented. The stiffness coefficients of the vibration absorber element are determined by a simple structural finite element modeling and calculation using ANSYS. Table 1 lists the calculated stiffness and damping coefficients of the front bearing as a function of the rotational speed.Table 1
Stiffness and damping dynamics coefficients of the front bearing as a function of the rotational speed.
Rotational speed[Hz]
Stiffness[N/m]
Damping[Ns/m]
0
1.37
×
10
7
39.3
83.33
1.37
×
10
7
33.4
167.67
1.37
×
10
7
19.4
250
1.37
×
10
7
13.5
333.33
1.37
×
10
7
9.3
416.67
1.37
×
10
7
7.0
500
1.37
×
10
7
5.1


### 2.3. Rear Bearing-Squeeze Film Damper

As stated previously, the rear bearing of the gas turbine consists of an unsealed squeeze film damper with a circumferential feeding groove. The type of ball bearing of the rear bearing is an angular contact. Figure2(c) depicts the rear bearing configuration in detail. The angular contact ball bearing is continuously lubricated by aeronautic oil. The nozzle serves to improve the lubrication and cooling of the bearing. Due to the high thermal loads from the combustion chamber, the ball bearing is capable of developing small axial motions along the housing because of dilatation. Thus, there is an axial centering spring, which allows axial movement of the angular contact ball bearing without the loss of functionality. Additionally, the centering spring ensures a good contact between the balls and the pathways and prevents the rotation of the outer ring. The axial preload of the system is 150 N. This preload affects the stiffness coefficients of the bearings. However, the variation on the stiffness coefficients caused by axial preload will not impact very much on rotordynamic results, since the axial preload of 150 N can be considered small according to the characteristics of the studied system (dimensions and materials). Therefore, for simplicity, axial preload has not been considered in the analysis. A radial clearance (cr) and an oil film of the squeeze film damper are formed between the outer surface of the external ring and the housing; that is, a radial clearance exists between the external surface of the outer ring and the housing, which is filled with oil. The circumferential feeding groove acts as a power supply of the oil, which divides the total film into two adjacent films. The groove and the adjacent oil films act as vibration absorbers to address special conditions. The geometric dimensions of the feeding groove and its interactions with the two oil film regions are considered to determine the dynamic stiffness and damping coefficients of the rear bearing. The outer ring of the angular contact ball bearing acts as a nonrotating bearing. The complete modeling of the rear bearing and the procedures for calculation of the stiffness and damping dynamic coefficients can be found in [36]. Cavitation effects by vaporization and π-film model were considered. Table 2 presents the calculated stiffness and damping dynamic coefficients of the rear bearing for three different radial clearances as a function of the rotational speed. The radial clearance is one of the most important factors used to determine the damping coefficient of a squeeze film damper, which has an optimum value for each type of application. Therefore, several computer simulations in rotordynamics were conducted to determine the influence of each radial clearance in the vibratory behavior of the developed gas turbine rotor-bearing system.Table 2
Stiffness and damping dynamics coefficients of the rear bearing as a function of the rotational speed for three different radial clearances of the squeeze film damper.
Rotational speed[Hz]
c
r
=
0.08 mm
c
r
=
0.12 mm
c
r
=
0.1905 mm
Stiffness[N/m]
Damping[Ns/m]
Stiffness[N/m]
Damping[Ns/m]
Stiffness[N/m]
Damping[Ns/m]
0
9.43
×
10
6
32,668
2.94
×
10
6
9,720
8.28
×
10
5
2,449
83.33
9.43
×
10
6
32,668
2.94
×
10
6
9,720
8.28
×
10
5
2,449
167.67
1.61
×
10
7
26,418
5.30
×
10
6
7,856
1.69
×
10
6
1,978
250
2.01
×
10
7
20,352
7.16
×
10
6
6,056
2.60
×
10
6
1,529
333.33
2.13
×
10
7
14,137
8.44
×
10
6
4,216
3.56
×
10
6
1,072
416.67
2.44
×
10
7
11,204
10.5
×
10
6
3,352
4.89
×
10
6
859,6
500
2.62
×
10
7
8,111
12.5
×
10
6
2,441
6.38
×
10
6
635,9


## 2.1. Technical Characteristics

The aeronautical gas turbine developed in this study has a single spool and is designed to be used in unmanned aerial vehicle applications. It has an air inlet duct, five-stage axial compressor, flow-through annular combustion chamber, single-stage axial turbine, and exhaustion nozzle. Figures1(a) and 1(b) illustrate the manufactured gas turbine currently in the experimental test phase.Figure 1
Manufactured single spool gas turbine developed in this study: (a) apparatus for demonstrative purposes and preliminary tests and (b) test rig for full test protocol.
(a)
(b)It should be noted that low weight and high thrust generation are highly desirable characteristics in smaller modern aircraft gas turbines. Based on user requirements, the dry weight of the developed gas turbine should be approximately 650 N, its maximum length should not exceed 1.50 m, and the maximum diameter of the circular section should not exceed 0.35 m. The operation is expected for a flight envelope from 0 to 1,200 m, maximum speeds of up to 0.9 Mach, acceleration of 10 g in the longitudinal axis and 7 g in the transverse axis, and cumulative operating time of at least 30 h. By having a single spool, the design, manufacture, and assembly of the rotor-bearing system occur more robustly, simply, and within a shorter amount of time. The operating speed range of the gas turbine is set between 80 and 100% of the nominal speed based on the operation at the design point. Therefore, its operating speed range is set between 375.33 and 469.17 Hz, which are the limits of the rotation speeds in which the gas turbine will mostly operate.Figure2 depicts a three-dimensional illustration of the cross-section of the developed gas turbine with details on the two bearings that support the rotor-bearing system. In Figure 2(a), the primary internal components of the aircraft gas turbine and its external general shape can be seen in detail, with a primary focus on the bipartite mounting of the axial compressor. The rotor-bearing system is supported by two bearings: a front or cold bearing and a rear or hot bearing. The front bearing consists of a deep groove ball bearing and a vibration absorber element. The rear bearing is an open-ended unsealed squeeze film damper with a circumferential feeding groove. The lubrication system has been carefully designed to meet the requirements of the bearings, which are continuously lubricated by aeronautical oil with specific performance characteristics.Figure 2
Three-dimensional illustration of the cross-section of the developed gas turbine: (a) bipartite mounting of the axial compressor and other inner components; (b) front or cold bearing; and (c) rear or hot bearing.
(a)
(b)
(c)


## 2.2. Front Bearing-Deep Groove Ball Bearing

As stated previously, the front bearing of the single spool gas turbine consists of a vibration absorber element and a deep groove ball bearing. Furthermore, the front bearing is constituted by lubricating channels to transport the oil from the lubrication system. Figure2(b) illustrates the configuration of the front bearing in detail. The deep groove ball bearing is continuously lubricated by aeronautic oil that can endure high temperature loads. In addition to lubrication, the oil dampens vibrations. The vibration absorber element provides a designed stiffness to the front bearing, which contributes to more acceptable vibration amplitudes. The stiffness of the vibration absorber element is lower than the stiffness of the deep groove ball bearing contacts and approximately constant for all rotational speeds of the rotor-bearing system. Thus, the characterization of the front bearing should be performed considering the influence of the stiffness of the vibration absorber element and the dynamic damping properties of the deep groove ball bearing contacts [9]. In [35], the procedures for characterization of the stiffness and damping dynamic coefficients of the deep groove ball bearing are presented. The stiffness coefficients of the vibration absorber element are determined by a simple structural finite element modeling and calculation using ANSYS. Table 1 lists the calculated stiffness and damping coefficients of the front bearing as a function of the rotational speed.Table 1
Stiffness and damping dynamics coefficients of the front bearing as a function of the rotational speed.
Rotational speed[Hz]
Stiffness[N/m]
Damping[Ns/m]
0
1.37
×
10
7
39.3
83.33
1.37
×
10
7
33.4
167.67
1.37
×
10
7
19.4
250
1.37
×
10
7
13.5
333.33
1.37
×
10
7
9.3
416.67
1.37
×
10
7
7.0
500
1.37
×
10
7
5.1


## 2.3. Rear Bearing-Squeeze Film Damper

As stated previously, the rear bearing of the gas turbine consists of an unsealed squeeze film damper with a circumferential feeding groove. The type of ball bearing of the rear bearing is an angular contact. Figure2(c) depicts the rear bearing configuration in detail. The angular contact ball bearing is continuously lubricated by aeronautic oil. The nozzle serves to improve the lubrication and cooling of the bearing. Due to the high thermal loads from the combustion chamber, the ball bearing is capable of developing small axial motions along the housing because of dilatation. Thus, there is an axial centering spring, which allows axial movement of the angular contact ball bearing without the loss of functionality. Additionally, the centering spring ensures a good contact between the balls and the pathways and prevents the rotation of the outer ring. The axial preload of the system is 150 N. This preload affects the stiffness coefficients of the bearings. However, the variation on the stiffness coefficients caused by axial preload will not impact very much on rotordynamic results, since the axial preload of 150 N can be considered small according to the characteristics of the studied system (dimensions and materials). Therefore, for simplicity, axial preload has not been considered in the analysis. A radial clearance (cr) and an oil film of the squeeze film damper are formed between the outer surface of the external ring and the housing; that is, a radial clearance exists between the external surface of the outer ring and the housing, which is filled with oil. The circumferential feeding groove acts as a power supply of the oil, which divides the total film into two adjacent films. The groove and the adjacent oil films act as vibration absorbers to address special conditions. The geometric dimensions of the feeding groove and its interactions with the two oil film regions are considered to determine the dynamic stiffness and damping coefficients of the rear bearing. The outer ring of the angular contact ball bearing acts as a nonrotating bearing. The complete modeling of the rear bearing and the procedures for calculation of the stiffness and damping dynamic coefficients can be found in [36]. Cavitation effects by vaporization and π-film model were considered. Table 2 presents the calculated stiffness and damping dynamic coefficients of the rear bearing for three different radial clearances as a function of the rotational speed. The radial clearance is one of the most important factors used to determine the damping coefficient of a squeeze film damper, which has an optimum value for each type of application. Therefore, several computer simulations in rotordynamics were conducted to determine the influence of each radial clearance in the vibratory behavior of the developed gas turbine rotor-bearing system.Table 2
Stiffness and damping dynamics coefficients of the rear bearing as a function of the rotational speed for three different radial clearances of the squeeze film damper.
Rotational speed[Hz]
c
r
=
0.08 mm
c
r
=
0.12 mm
c
r
=
0.1905 mm
Stiffness[N/m]
Damping[Ns/m]
Stiffness[N/m]
Damping[Ns/m]
Stiffness[N/m]
Damping[Ns/m]
0
9.43
×
10
6
32,668
2.94
×
10
6
9,720
8.28
×
10
5
2,449
83.33
9.43
×
10
6
32,668
2.94
×
10
6
9,720
8.28
×
10
5
2,449
167.67
1.61
×
10
7
26,418
5.30
×
10
6
7,856
1.69
×
10
6
1,978
250
2.01
×
10
7
20,352
7.16
×
10
6
6,056
2.60
×
10
6
1,529
333.33
2.13
×
10
7
14,137
8.44
×
10
6
4,216
3.56
×
10
6
1,072
416.67
2.44
×
10
7
11,204
10.5
×
10
6
3,352
4.89
×
10
6
859,6
500
2.62
×
10
7
8,111
12.5
×
10
6
2,441
6.38
×
10
6
635,9


## 3. Rotordynamic Model for FE Analysis

The studied system is nonlinear; however, the modeling is linear. Despite being linear, the modeling approach allows a good representation of the studied system and the achievement of accurate rotordynamic results. The general rotordynamics equation in a stationary reference frame can be written as(1)Mu¨+C+Gu˙+K+Bu=F,where [M], [C], and [K] are the standard mass, damping, and stiffness matrices. {F} is the external force vector which contains all the forcing functions. [G] is the skew symmetric gyroscopic matrix which depends on the rotational velocity and captures the gyroscopic effect which couples rotational degrees of freedom perpendicular to spinning axis. [B] is the rotating damping matrix which also depends on rotational velocity and modifies apparent stiffness of the structure and can produce unstable motion. The classic linearized model with eight spring and damping coefficients is used to model the bearings. The forces at each bearing are assumed to obey the governing equations of the following form: (2)cxxcxycyxcyyu˙xu˙y+kxxkxykyxkyyuxuy=fxfy,where cij and kij are the bearing damping and stiffness coefficients, respectively.Figure3 depicts the manufactured rotor-bearing system of the developed gas turbine with its respective discs and locations of the bearings. For the front bearing, the stiffness coefficients are given by a finite element analysis of the vibration absorber element, which is implied in kyy=kzz and kyz=kzy=0. The stiffness coefficients kyy=kzz are considered the same for all rotor speeds, since the variation is very small and, thus, can be neglected. The dynamic damping properties of the deep groove ball bearing are given by [35]. In a deep groove ball bearing the damping coefficients are very small. The cross-coupled terms are neglected in this study. Therefore cyy=czz and cyz=czy=0. Regarding the rear bearing, the squeeze film damper used in our work is a short open-end squeeze film damper with a circumferential feeding groove. The mathematical model used to obtain the dynamic stiffness and damping coefficients has been implemented according to [36]. Only the principal force components fx and fy are calculated since it is a short open-end squeeze film damper with a circumferential feeding groove, whose design considerations assume very small eccentricity and circular synchronous centered motions. Therefore kyy=kzz, cyy=czz, and kyz=kzy=cyz=czy=0. Tables 1 and 2 list all of the numerical values used in the calculations.Figure 3
Manufactured rotor-bearing system of the developed gas turbine: (a) compressor discs; (b) free turbine disc; and (c) locations of the bearings.
(a)
(b)
(c)Figure4 illustrates a schematic model for the FE analysis of the developed rotor-bearing system with its primary nodes. This is a simplified model because the blades and other complex geometries of the rotor-bearing system are suppressed for finite element discretization. The centers of mass and moments of inertia of the discs are considered to derive this simplified model. The shaft is made of AISI 4340 normalized steel. Young’s modulus of AISI 4340 normalized steel is 205 GPa, Poisson’s ratio is 0.29, and the density is 7,850 kg/m3. The discs are made of different materials to withstand the high gradient of the mechanical and thermal loads. Table 3 lists the primary physical and geometric properties of the discs. All of the discs are subjected to high centrifugal loads. The first three discs of the axial compressor are made of aluminum alloys. The fourth and fifth discs are made of titanium alloy because the temperatures above these discs are significantly higher. The sixth disc is a free turbine disc, which consists of an Inconel alloy to withstand the high temperature loads from the postcombustion gases.Table 3
Physical and geometric properties of the discs.
Discs
D1
D2
D3
D4
D5
D6
Material
Al7178-T6
Al7178-T6
Al2024-T8
Ti-6Al-4V
Ti-6Al-4V
In713LC
Young’s modulus, GPa
71.7
71.7
72.4
113.8
113.8
163.3
Density, kg/m3
2,830
2,830
2,870
4,430
4,430
8,000
Poisson’s ratio
0.33
0.33
0.33
0.342
0.342
0.382
Width, m
0.02
0.02
0.02
0.02
0.02
0.03
Inner diameter, m
0.08
0.08
0.08
0.08
0.08
0.055
Outer diameter, m
0.19
0.21
0.23
0.24
0.25
0.25Figure 4
Rotor-bearing system for the FE analysis of the developed gas turbine with its primary nodes. All dimensions are in meters.The finite element model was developed using the ANSYS Parametric Design Language (APDL). The effects of the rotary inertia, gyroscopic moments, axial loads, shear deformations, and internal and external damping are included in the formulations. All of the elements in the model have 4 degrees of freedom at each node. These include translations in the nodaly- andz-directions and rotations about the nodaly- andz-axes. The shaft elements are modeled using Timoshenko beam quadratic elements of the type BEAM189. The discs are modeled using PIPE16 elements, which are elastic straight uniaxial elements with tension compression, torsion, and bending capabilities. The front and rear bearings are modeled using MATRIX27 elements, which are matrix elements that have an undefined geometry but an elastic kinematic response that can be specified by the stiffness and damping coefficients as a function of the rotational speed.Due to its importance and the high degree of dangerousness of the application, it is evident that a full rotordynamic analysis needs to be performed in the rotor-bearing system considering all possible vibration problems. Initially, it is essential to verify whether the primary system critical speeds are near or far from the operating speed range. As the rotational speed approaches a critical speed of the system, the vibration levels become higher, and the existing alternating stresses also increase, thus increasing the risk of mechanical interference of the components and the occurrence of fatigue failure. Because the number of resonance points is virtually unlimited, it is necessary to assume that certain resonance points are inevitable and eliminate those that may be potentially problematic. The vibration amplitudes tend to decrease for a given vibrational motion as the frequency increases. Therefore, it is clear that the resonance points at higher frequencies have less importance. Thus, the analyst must define which resonance points should be eliminated and which ones can be tolerated. As a general rule, the resonance points in the lowest natural frequencies derived by the stronger excitation frequencies should be avoided, especially in the operating speed range, because they tend to have the greatest vibration amplitudes.By observing Figure3 and the schematic model of the rotor-bearing system depicted in Figure 4, it is clear that the most critical component of the rotor-bearing system is the single-stage free turbine disc (D6). This fact can be noted because the disc D6 is a free disc, which is located after the rear bearing. Therefore, the existence of any significant mass unbalance over this disc or any excitation of a vibration mode over the shaft will cause large amplitudes of vibration with a large possibility of fatigue failure in the region of nodes 18, 19, and 20.


## 4. Results and Discussions

A balanced quality grade G1 from the ISO standard 1940/1 was used to manufacture the rotor-bearing system. For the rotordynamic numerical analyses, a mass unbalance of1×10-4 kg·m located at node 21 of the finite element model was used. The vibration amplitudes were analyzed in both the five-stage axial compressor and the free turbine disc because the tip blade clearances of these components were defined to be as small as possible for maximum fluid dynamic efficiency.


### 4.1. Campbell Diagrams

The groove depth and other squeeze film damper-related parameters were kept the same for all three radial clearances (cr=0.08 mm, cr=0.12 mm, and cr=0.1905 mm). The squeeze film damper with a radial clearance equal to 0.08 mm has the highest generated damping values. Because the bearings are symmetrical, only modes in Forward Whirl (FW) can potentially be induced. Figure 5 depicts the Campbell diagram for the radial clearance of the squeeze film damper equal to 0.08 mm. It can be seen that, due to the high generated damping values, the natural frequencies, notably modes 1FW and 1BW, have a strong nonlinearity. By considering the synchronous excitation line (1Ω), there is only one critical speed: 2FW @ 147.1 Hz. Other points of resonance in relation to the synchronous excitation line occur in Backward Whirl (BW) and thus cannot be induced. The critical speed of 2FW @ 147.1 Hz is far from the operating speed range. Other critical speeds of interest that can cause potential vibration problems are encountered when considering the 0.5Ω excitation line. The following critical speeds can be highlighted: 2FW @ 305.1 Hz and 1FW @ 360.9 Hz. Based on these critical speeds, 1FW @ 360.9 Hz is extremely close to the operating speed range; thus, this condition is undesirable and should be avoided.Figure 5
Campbell diagram for the radial clearance of the squeeze film damper equal to 0.08 mm.The squeeze film damper with a radial clearance equal to 0.12 mm has intermediate damping values. It can be seen in the Campbell diagram in Figure6 that the natural frequencies have less nonlinearity when compared with the natural frequencies of the diagram in Figure 5. However, there are two critical speeds when considering the synchronous excitation line: 1FW @ 50.43 Hz and 2FW @ 149.1 Hz. The critical speed of 2FW @ 149.1 Hz is far from the operating speed range and can thus be considered safe. Only one critical speed can be highlighted with respect to the 0.5Ω excitation frequency: 2FW @ 318.8 Hz. This critical speed is outside the operating speed range by a significant margin of distance; hence, it is possible to consider it safe. There are no critical speeds within the operating speed range.Figure 6
Campbell diagram for the radial clearance of the squeeze film damper equal to 0.12 mm.The squeeze film damper with a radial clearance equal to 0.1905 mm has the lowest generated damping values. It can be seen in the Campbell diagram in Figure7 that the variation in the natural frequencies for this clearance is similar to the variation presented in the Campbell diagram in Figure 6. By considering the synchronous excitation line, the existence of two critical speeds can be observed: 1FW @ 32.5 Hz and 2FW @ 149.1 Hz. The critical speed of 2FW @ 149.1 Hz is far from the operating speed range and can thus be considered safe. Only one critical speed can be highlighted with respect to the 0.5Ω excitation frequency: 2FW @ 315.3 Hz. This critical speed is outside the operating speed range by a significant margin of distance. Therefore, it can be considered safe. There are no critical speeds within the operating speed range.Figure 7
Campbell diagram for the radial clearance of the squeeze film damper equal to 0.1905 mm.By considering the Campbell diagrams presented for the three different radial clearances, it should be noted that the critical speeds that could result from the intersections of the natural frequencies with respect to the 2Ω excitation line were not considered in the analyses because there is no coupling in the rotating shaft, and this 2Ω excitation frequency is thus unlikely to occur with significant levels of energy. Furthermore, the other critical speeds with respect to thesΩ excitation lines were discharged because they only can appear at extremely high frequencies and hence have very low amplitudes of vibration. These excitation frequencies have limited possibilities of existing with significant levels of energy to cause vibration problems in the system.


### 4.2. Unbalanced Responses

Figure8 illustrates the unbalanced responses of the rotor-bearing system when the radial clearance of the squeeze film damper equals 0.08 mm. It is possible to see that the critical speed of 2FW @ 147.1 Hz is the most troublesome. The vibration amplitudes were analyzed at the five-stage axial compressor and the free turbine disc. The vibration amplitudes reached 5.4 μm at the axial compressor region and 11.2 μm at the single-stage free turbine disc for an operating speed range of 22,520 to 28,150 rpm. The vibration amplitudes due to the critical speed of 2FW @ 147.1 Hz reached 117.6 μm over the axial compressor region and 58.1 μm at the single-stage turbine disc.Figure 8
Unbalanced responses of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.08 mm.
(a)
(b)Figure9 depicts the unbalanced responses of the rotor-bearing system when the radial clearance of the squeeze film damper equals 0.12 mm. It is possible to see that the critical speed of 1FW @ 50.43 Hz is largely attenuated by the damping performance of the squeeze film damper. The critical speed of 2FW @ 149.1 Hz is the most troublesome. The vibration amplitudes reached 5.9 μm at the five-stage axial compressor region and 8.6 mm at the single-stage free turbine disc for an operating speed range of 22,520 to 28,150 rpm. The vibration amplitudes due to the critical speed of 2FW @ 149.1 Hz reached 189.1 mm at the five-stage axial compressor region and 113.3 mm at the free turbine disc region.Figure 9
Unbalanced responses of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.12 mm.
(a)
(b)Figure10 depicts the unbalanced responses of the rotor-bearing system when the radial clearance of the squeeze film damper equals 0.1905 mm. It is possible to observe the presence of an amplitude peak related to the critical speed of 1FW @ 32.5 Hz. This peak occurs due to the lower generated damping values of the squeeze film damper with a radial clearance of 0.1905 mm. The critical speed of 2FW @ 149.1 Hz is the most troublesome. The vibration amplitudes reached 6.2 μm at the five-stage axial compressor region and 8.0 μm at the single-stage free turbine disc region for an operating speed range of 22,520 to 28,150 rpm. The vibration amplitudes due to the critical speed of 1FW @ 32.5 Hz reached 5.0 μm for the five-stage axial compressor region and 11.5 μm for the free turbine disc. The vibration amplitudes due to the critical speed of 2FW @ 149.1 Hz reached 329.5 μm at the axial compressor region and 190.7 μm at the single-stage free turbine disc.Figure 10
Unbalanced responses of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.1905 mm.
(a)
(b)


### 4.3. Transient Responses

Transient analyses were performed to simulate the transition of the rotor-bearing system through the primary resonance points. A linear acceleration model is considered for the systemω=ω0+αt. The angular acceleration during the acceleration phase of the system is set to 314.15 rad/s2, since a time of about 10 seconds for the aeronautical gas turbine to reach maximum rotational speed is expected. Figure 11 depicts the total displacement at node 21 of the rotor-bearing system when the radial clearance is 0.08 mm for the squeeze film damper. The total displacement is obtained by the sum of the horizontal displacement (z) and the vertical displacement (y). A slight salience appears in the solution approximately 6.45 seconds after starting the gas turbine. This salience is related to the critical speed of 2FW @ 305.1 Hz for the 0.5Ω excitation frequency. This result can be observed by the Campbell diagram in Figure 5. This amplitude of vibration at the 0.5Ω frequency occurs smoothly due to the high damping performance of the squeeze film damper with a radial clearance of 0.08 mm. The vibration amplitudes during the passage through the resonance point of 2FW @ 147.1 Hz reached 26.6 μm, thus demonstrating the strong attenuation by the rear bearing with 0.08 mm of radial clearance.Figure 11
Total displacement at node 21 of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.08 mm.Figure12 depicts the transient response of the rotor-bearing system for a radial clearance of 0.12 mm for the total displacements at node 21. A small peak appears in the solution approximately 6.45 seconds after the gas turbine starts up. This small peak is related to the critical speed of 2FW @ 318.8 Hz for the 0.5Ω excitation frequency. This result can be observed by the Campbell diagram in Figure 6. This vibration at the 0.5Ω frequency occurs smoothly due to the damping performance of the squeeze film damper for a radial clearance of 0.12 mm. The vibration amplitudes during the passage through the resonance point of 2FW @ 149.1 Hz reached 44.2 μm.Figure 12
Total displacement at node 21 of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.12 mm.Figure13 depicts the transient response of the rotor-bearing system for a radial clearance of 0.1905 mm for the total displacements at node 21. A significant peak appears in the solution approximately 6.45 seconds after the gas turbine starts up. This peak is related to the critical speed of 2FW @ 315.3 Hz for the 0.5Ω excitation frequency. This result can be observed by the Campbell diagram in Figure 7. This amplitude of vibration at the 0.5Ω frequency occurs in a significant manner due to the low damping performance of the squeeze film damper with a radial clearance of 0.1905 mm. The vibration amplitudes during the passage through the resonance point of 1FW @ 32.5 Hz reached 11.3 μm. The vibration amplitudes reached 47.7 μm for the resonance point of 2PD @ 149.1 Hz.Figure 13
Total displacement at node 21 of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.1905 mm.


### 4.4. Selection of the Squeeze Film Damper

From the analysis presented above, it is possible to select the best radial clearance for the squeeze film damper and, consequently, the best combination of bearings for the rotor-bearing system. It is clear that, for a radial clearance of 0.08 mm, the greatest attenuation occurs in the amplitudes of the response. Nevertheless, for a clearance of 0.08 mm, the 1FW at 360.9 Hz resonance point is extremely close to the operating range. Comparing the numerical analyses for the clearances of 0.12 and 0.1905 mm, it is possible to see that the Campbell diagrams for these two conditions are similar. However, because the 0.12 mm radial clearance provides higher damping values, the first critical speed with respect to the synchronous excitation line is considerably reduced, thus making it unnoticeable. Additionally, the unbalanced response amplitudes indicated lower values compared to the 0.1905 mm radial clearance. Thus, although the response amplitudes for the radial clearance of 0.1905 mm is also acceptable, considering the tip clearances, the radial clearance of 0.12 mm was selected to be used in the design of the squeeze film damper of the rear bearing of the developed single spool gas turbine to ensure a greater safety margin for the system.


### 4.5. Stability Analysis of the System

It is of fundamental importance to predict the stability of the rotor-bearing system from rest to maximum speed. One way to check the dynamic stability of a rotor-bearing system is to calculate the logarithmic decrements in their respective natural frequencies. During the Campbell diagram calculations, the ANSYS software provides the real and imaginary parts of the natural frequencies and the respective damping rate. The calculation of the logarithm decrementξ is done as a function of the damping rate ζ, according to the following equation: (3)ξ=4π21/ζ2-1.It is well known that internal damping can destabilize a rotor-bearing system. Furthermore, numerical formulations, which include hysteretic damping effects, can occasionally lead to incorrect prediction system instability ranges. Therefore, in this work we calculate the values of the logarithmic decrements throughout the rotational speed. The instabilities are determined when the values of the logarithmic decrements are less than zero. Table 4 presents the logarithmic decrements calculated for the first six natural frequencies of the system. As seen in Table 4, there is no occurrence of instabilities because all of the logarithmic decrements in the system are positive.Table 4
Logarithmic decrements calculated for the six first natural frequencies of the rotor-bearing system.
Rotational speed[Hz]
1FWmode
1BWmode
2FWmode
2BWmode
3FWmode
3BWmode
55.56
5.0437
5.2713
0.2024
0.0216
0.8447
0.8903
111.11
3.3526
3.8727
0.4609
0.0234
0.7276
0.7742
166.67
1.0670
2.4758
1.3331
0.0453
0.6332
0.5997
222.22
0.3705
1.7818
1.3666
0.0743
0.5428
0.4473
277.78
0.1605
1.2929
1.0972
0.1012
0.4522
0.3222
333.33
0.0779
0.9097
0.8126
0.1171
0.3563
0.2216
388.89
0.0419
0.6765
0.6255
0.1460
0.3253
0.1705
444.44
0.0256
0.4636
0.4822
0.1671
0.2886
0.1283
500
0.0171
0.3288
0.3251
0.1720
0.2445
0.0934


### 4.6. Alternating Bending Stresses

The alternating bending stresses were evaluated at several specific points of the rotor-bearing system. From Figure4, it can be seen that node 20 is located in one of the most critical locations; that is, it is located after the rear bearing and before the free turbine disc. Additionally, node 20 is situated in a region with a significant transition in the cross-section area. The smallest cross-sectional area of the shaft is located at this node. Figure 14 depicts the alternating bending stress calculated at node 20 during the transition through the primary resonance points of the rotor-bearing system. As shown, extremely low stress values were found throughout the transient acceleration phase of the gas turbine. The greatest stress values on node 20 occur during the passage through the second critical speed, which is 4.4 MPa. For the critical speed with respect to the 0.5Ω excitation line, the magnitude of the alternating bending stresses can reach 1.7 MPa. By assuming that the gas turbine acceleration phase has been completed, it can be observed that the alternating bending stresses within the operating speed range of the gas turbine will not exceed 1.8 MPa. These stress values are extremely low considering the material and dimensions of the shaft cross-section. Thus, the service life of the shaft will not be compromised. Therefore, an infinite life of the shaft to fatigue failure can be expected. The alternating bending stresses calculated for other nodes of the shaft indicated lower values and were less significant.Figure 14
Alternating bending stresses evaluated at node 20 during the system transition through resonance.


### 4.7. Spectral Map of the Rotor-Bearing System

In addition to the analyses presented above, it is also important to plot the spectral map of the rotor-bearing system considering a broader range of excitation frequencies. In this spectral map, the vibration amplitudes, which are given by unbalance responses, are plotted as a function of the rotational speed and of the various excitation frequencies that may exist in the system. The number of excitation frequencies and intervals between the values of the rotational speeds are defined by the amount of substeps. Thus, the regions where the amplitudes of vibration are undesirable can be defined, considering the interaction between the rotational speeds and excitation frequencies. Figure15 illustrates the spectral map of the rotor-bearing system when the radial clearance of the squeeze film damper is equal to 0.12 mm. It was considered that the rotational speed of the rotor-bearing system ranges from 0 to 500 Hz in intervals of 50 substeps. The excitation frequencies range from 0 to 500 Hz in 250 substep intervals. The vibration amplitudes are calculated on node 21 using clusters with 5 Hz of bandwidth and a minimum frequency of 0.5 Hz. The higher vibration amplitudes occur due to the second critical speed. The first critical speed can be observed in the rotational speed of the rotor-bearing system at 500 Hz, with an excitation frequency of 125 Hz. The 0.5Ω subharmonic vibrations may be observed at a rotational speed of approximately 330 Hz with an excitation frequency of 500 Hz. It can be seen in the spectral map that the operating speed range of the rotor-bearing system is free from undesirable vibration amplitudes in all excitation frequencies and rotational speeds of interest.Figure 15
Spectral map of the developed rotor-bearing system.


## 4.1. Campbell Diagrams

The groove depth and other squeeze film damper-related parameters were kept the same for all three radial clearances (cr=0.08 mm, cr=0.12 mm, and cr=0.1905 mm). The squeeze film damper with a radial clearance equal to 0.08 mm has the highest generated damping values. Because the bearings are symmetrical, only modes in Forward Whirl (FW) can potentially be induced. Figure 5 depicts the Campbell diagram for the radial clearance of the squeeze film damper equal to 0.08 mm. It can be seen that, due to the high generated damping values, the natural frequencies, notably modes 1FW and 1BW, have a strong nonlinearity. By considering the synchronous excitation line (1Ω), there is only one critical speed: 2FW @ 147.1 Hz. Other points of resonance in relation to the synchronous excitation line occur in Backward Whirl (BW) and thus cannot be induced. The critical speed of 2FW @ 147.1 Hz is far from the operating speed range. Other critical speeds of interest that can cause potential vibration problems are encountered when considering the 0.5Ω excitation line. The following critical speeds can be highlighted: 2FW @ 305.1 Hz and 1FW @ 360.9 Hz. Based on these critical speeds, 1FW @ 360.9 Hz is extremely close to the operating speed range; thus, this condition is undesirable and should be avoided.Figure 5
Campbell diagram for the radial clearance of the squeeze film damper equal to 0.08 mm.The squeeze film damper with a radial clearance equal to 0.12 mm has intermediate damping values. It can be seen in the Campbell diagram in Figure6 that the natural frequencies have less nonlinearity when compared with the natural frequencies of the diagram in Figure 5. However, there are two critical speeds when considering the synchronous excitation line: 1FW @ 50.43 Hz and 2FW @ 149.1 Hz. The critical speed of 2FW @ 149.1 Hz is far from the operating speed range and can thus be considered safe. Only one critical speed can be highlighted with respect to the 0.5Ω excitation frequency: 2FW @ 318.8 Hz. This critical speed is outside the operating speed range by a significant margin of distance; hence, it is possible to consider it safe. There are no critical speeds within the operating speed range.Figure 6
Campbell diagram for the radial clearance of the squeeze film damper equal to 0.12 mm.The squeeze film damper with a radial clearance equal to 0.1905 mm has the lowest generated damping values. It can be seen in the Campbell diagram in Figure7 that the variation in the natural frequencies for this clearance is similar to the variation presented in the Campbell diagram in Figure 6. By considering the synchronous excitation line, the existence of two critical speeds can be observed: 1FW @ 32.5 Hz and 2FW @ 149.1 Hz. The critical speed of 2FW @ 149.1 Hz is far from the operating speed range and can thus be considered safe. Only one critical speed can be highlighted with respect to the 0.5Ω excitation frequency: 2FW @ 315.3 Hz. This critical speed is outside the operating speed range by a significant margin of distance. Therefore, it can be considered safe. There are no critical speeds within the operating speed range.Figure 7
Campbell diagram for the radial clearance of the squeeze film damper equal to 0.1905 mm.By considering the Campbell diagrams presented for the three different radial clearances, it should be noted that the critical speeds that could result from the intersections of the natural frequencies with respect to the 2Ω excitation line were not considered in the analyses because there is no coupling in the rotating shaft, and this 2Ω excitation frequency is thus unlikely to occur with significant levels of energy. Furthermore, the other critical speeds with respect to thesΩ excitation lines were discharged because they only can appear at extremely high frequencies and hence have very low amplitudes of vibration. These excitation frequencies have limited possibilities of existing with significant levels of energy to cause vibration problems in the system.


## 4.2. Unbalanced Responses

Figure8 illustrates the unbalanced responses of the rotor-bearing system when the radial clearance of the squeeze film damper equals 0.08 mm. It is possible to see that the critical speed of 2FW @ 147.1 Hz is the most troublesome. The vibration amplitudes were analyzed at the five-stage axial compressor and the free turbine disc. The vibration amplitudes reached 5.4 μm at the axial compressor region and 11.2 μm at the single-stage free turbine disc for an operating speed range of 22,520 to 28,150 rpm. The vibration amplitudes due to the critical speed of 2FW @ 147.1 Hz reached 117.6 μm over the axial compressor region and 58.1 μm at the single-stage turbine disc.Figure 8
Unbalanced responses of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.08 mm.
(a)
(b)Figure9 depicts the unbalanced responses of the rotor-bearing system when the radial clearance of the squeeze film damper equals 0.12 mm. It is possible to see that the critical speed of 1FW @ 50.43 Hz is largely attenuated by the damping performance of the squeeze film damper. The critical speed of 2FW @ 149.1 Hz is the most troublesome. The vibration amplitudes reached 5.9 μm at the five-stage axial compressor region and 8.6 mm at the single-stage free turbine disc for an operating speed range of 22,520 to 28,150 rpm. The vibration amplitudes due to the critical speed of 2FW @ 149.1 Hz reached 189.1 mm at the five-stage axial compressor region and 113.3 mm at the free turbine disc region.Figure 9
Unbalanced responses of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.12 mm.
(a)
(b)Figure10 depicts the unbalanced responses of the rotor-bearing system when the radial clearance of the squeeze film damper equals 0.1905 mm. It is possible to observe the presence of an amplitude peak related to the critical speed of 1FW @ 32.5 Hz. This peak occurs due to the lower generated damping values of the squeeze film damper with a radial clearance of 0.1905 mm. The critical speed of 2FW @ 149.1 Hz is the most troublesome. The vibration amplitudes reached 6.2 μm at the five-stage axial compressor region and 8.0 μm at the single-stage free turbine disc region for an operating speed range of 22,520 to 28,150 rpm. The vibration amplitudes due to the critical speed of 1FW @ 32.5 Hz reached 5.0 μm for the five-stage axial compressor region and 11.5 μm for the free turbine disc. The vibration amplitudes due to the critical speed of 2FW @ 149.1 Hz reached 329.5 μm at the axial compressor region and 190.7 μm at the single-stage free turbine disc.Figure 10
Unbalanced responses of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.1905 mm.
(a)
(b)


## 4.3. Transient Responses

Transient analyses were performed to simulate the transition of the rotor-bearing system through the primary resonance points. A linear acceleration model is considered for the systemω=ω0+αt. The angular acceleration during the acceleration phase of the system is set to 314.15 rad/s2, since a time of about 10 seconds for the aeronautical gas turbine to reach maximum rotational speed is expected. Figure 11 depicts the total displacement at node 21 of the rotor-bearing system when the radial clearance is 0.08 mm for the squeeze film damper. The total displacement is obtained by the sum of the horizontal displacement (z) and the vertical displacement (y). A slight salience appears in the solution approximately 6.45 seconds after starting the gas turbine. This salience is related to the critical speed of 2FW @ 305.1 Hz for the 0.5Ω excitation frequency. This result can be observed by the Campbell diagram in Figure 5. This amplitude of vibration at the 0.5Ω frequency occurs smoothly due to the high damping performance of the squeeze film damper with a radial clearance of 0.08 mm. The vibration amplitudes during the passage through the resonance point of 2FW @ 147.1 Hz reached 26.6 μm, thus demonstrating the strong attenuation by the rear bearing with 0.08 mm of radial clearance.Figure 11
Total displacement at node 21 of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.08 mm.Figure12 depicts the transient response of the rotor-bearing system for a radial clearance of 0.12 mm for the total displacements at node 21. A small peak appears in the solution approximately 6.45 seconds after the gas turbine starts up. This small peak is related to the critical speed of 2FW @ 318.8 Hz for the 0.5Ω excitation frequency. This result can be observed by the Campbell diagram in Figure 6. This vibration at the 0.5Ω frequency occurs smoothly due to the damping performance of the squeeze film damper for a radial clearance of 0.12 mm. The vibration amplitudes during the passage through the resonance point of 2FW @ 149.1 Hz reached 44.2 μm.Figure 12
Total displacement at node 21 of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.12 mm.Figure13 depicts the transient response of the rotor-bearing system for a radial clearance of 0.1905 mm for the total displacements at node 21. A significant peak appears in the solution approximately 6.45 seconds after the gas turbine starts up. This peak is related to the critical speed of 2FW @ 315.3 Hz for the 0.5Ω excitation frequency. This result can be observed by the Campbell diagram in Figure 7. This amplitude of vibration at the 0.5Ω frequency occurs in a significant manner due to the low damping performance of the squeeze film damper with a radial clearance of 0.1905 mm. The vibration amplitudes during the passage through the resonance point of 1FW @ 32.5 Hz reached 11.3 μm. The vibration amplitudes reached 47.7 μm for the resonance point of 2PD @ 149.1 Hz.Figure 13
Total displacement at node 21 of the rotor-bearing system for the radial clearance of the squeeze film damper equal to 0.1905 mm.


## 4.4. Selection of the Squeeze Film Damper

From the analysis presented above, it is possible to select the best radial clearance for the squeeze film damper and, consequently, the best combination of bearings for the rotor-bearing system. It is clear that, for a radial clearance of 0.08 mm, the greatest attenuation occurs in the amplitudes of the response. Nevertheless, for a clearance of 0.08 mm, the 1FW at 360.9 Hz resonance point is extremely close to the operating range. Comparing the numerical analyses for the clearances of 0.12 and 0.1905 mm, it is possible to see that the Campbell diagrams for these two conditions are similar. However, because the 0.12 mm radial clearance provides higher damping values, the first critical speed with respect to the synchronous excitation line is considerably reduced, thus making it unnoticeable. Additionally, the unbalanced response amplitudes indicated lower values compared to the 0.1905 mm radial clearance. Thus, although the response amplitudes for the radial clearance of 0.1905 mm is also acceptable, considering the tip clearances, the radial clearance of 0.12 mm was selected to be used in the design of the squeeze film damper of the rear bearing of the developed single spool gas turbine to ensure a greater safety margin for the system.


## 4.5. Stability Analysis of the System

It is of fundamental importance to predict the stability of the rotor-bearing system from rest to maximum speed. One way to check the dynamic stability of a rotor-bearing system is to calculate the logarithmic decrements in their respective natural frequencies. During the Campbell diagram calculations, the ANSYS software provides the real and imaginary parts of the natural frequencies and the respective damping rate. The calculation of the logarithm decrementξ is done as a function of the damping rate ζ, according to the following equation: (3)ξ=4π21/ζ2-1.It is well known that internal damping can destabilize a rotor-bearing system. Furthermore, numerical formulations, which include hysteretic damping effects, can occasionally lead to incorrect prediction system instability ranges. Therefore, in this work we calculate the values of the logarithmic decrements throughout the rotational speed. The instabilities are determined when the values of the logarithmic decrements are less than zero. Table 4 presents the logarithmic decrements calculated for the first six natural frequencies of the system. As seen in Table 4, there is no occurrence of instabilities because all of the logarithmic decrements in the system are positive.Table 4
Logarithmic decrements calculated for the six first natural frequencies of the rotor-bearing system.
Rotational speed[Hz]
1FWmode
1BWmode
2FWmode
2BWmode
3FWmode
3BWmode
55.56
5.0437
5.2713
0.2024
0.0216
0.8447
0.8903
111.11
3.3526
3.8727
0.4609
0.0234
0.7276
0.7742
166.67
1.0670
2.4758
1.3331
0.0453
0.6332
0.5997
222.22
0.3705
1.7818
1.3666
0.0743
0.5428
0.4473
277.78
0.1605
1.2929
1.0972
0.1012
0.4522
0.3222
333.33
0.0779
0.9097
0.8126
0.1171
0.3563
0.2216
388.89
0.0419
0.6765
0.6255
0.1460
0.3253
0.1705
444.44
0.0256
0.4636
0.4822
0.1671
0.2886
0.1283
500
0.0171
0.3288
0.3251
0.1720
0.2445
0.0934


## 4.6. Alternating Bending Stresses

The alternating bending stresses were evaluated at several specific points of the rotor-bearing system. From Figure4, it can be seen that node 20 is located in one of the most critical locations; that is, it is located after the rear bearing and before the free turbine disc. Additionally, node 20 is situated in a region with a significant transition in the cross-section area. The smallest cross-sectional area of the shaft is located at this node. Figure 14 depicts the alternating bending stress calculated at node 20 during the transition through the primary resonance points of the rotor-bearing system. As shown, extremely low stress values were found throughout the transient acceleration phase of the gas turbine. The greatest stress values on node 20 occur during the passage through the second critical speed, which is 4.4 MPa. For the critical speed with respect to the 0.5Ω excitation line, the magnitude of the alternating bending stresses can reach 1.7 MPa. By assuming that the gas turbine acceleration phase has been completed, it can be observed that the alternating bending stresses within the operating speed range of the gas turbine will not exceed 1.8 MPa. These stress values are extremely low considering the material and dimensions of the shaft cross-section. Thus, the service life of the shaft will not be compromised. Therefore, an infinite life of the shaft to fatigue failure can be expected. The alternating bending stresses calculated for other nodes of the shaft indicated lower values and were less significant.Figure 14
Alternating bending stresses evaluated at node 20 during the system transition through resonance.


## 4.7. Spectral Map of the Rotor-Bearing System

In addition to the analyses presented above, it is also important to plot the spectral map of the rotor-bearing system considering a broader range of excitation frequencies. In this spectral map, the vibration amplitudes, which are given by unbalance responses, are plotted as a function of the rotational speed and of the various excitation frequencies that may exist in the system. The number of excitation frequencies and intervals between the values of the rotational speeds are defined by the amount of substeps. Thus, the regions where the amplitudes of vibration are undesirable can be defined, considering the interaction between the rotational speeds and excitation frequencies. Figure15 illustrates the spectral map of the rotor-bearing system when the radial clearance of the squeeze film damper is equal to 0.12 mm. It was considered that the rotational speed of the rotor-bearing system ranges from 0 to 500 Hz in intervals of 50 substeps. The excitation frequencies range from 0 to 500 Hz in 250 substep intervals. The vibration amplitudes are calculated on node 21 using clusters with 5 Hz of bandwidth and a minimum frequency of 0.5 Hz. The higher vibration amplitudes occur due to the second critical speed. The first critical speed can be observed in the rotational speed of the rotor-bearing system at 500 Hz, with an excitation frequency of 125 Hz. The 0.5Ω subharmonic vibrations may be observed at a rotational speed of approximately 330 Hz with an excitation frequency of 500 Hz. It can be seen in the spectral map that the operating speed range of the rotor-bearing system is free from undesirable vibration amplitudes in all excitation frequencies and rotational speeds of interest.Figure 15
Spectral map of the developed rotor-bearing system.


## 5. Conclusions

It was observed that the radial clearance of the squeeze film damper strongly influence the damping performance of the rear bearing. For smaller values of radial clearance, the damping generated in the bearing is higher, as described in literature. Regarding the oil viscosity, the bearing damping capacity decreases with the oil temperature because the lubricating oil loses viscosity as the temperature increases. With respect to the bearing stiffness, as the viscosity of the lubricant decreases, the bearing stiffness increases. From the Campbell diagrams, unbalanced responses, and transient analyses, the radial clearance of 0.12 mm proved to be the most suitable for the developed single spool gas turbine squeeze film damper and was thus used for the final design of the rear bearing. With respect to the system stability analysis, the logarithmic decrements were calculated, and there was no occurrence of instabilities. Furthermore, the shaft alternating bending stresses were evaluated during the acceleration phase of the rotor-bearing system. Extremely low values for the bending stresses were observed during the passage through the resonance points. Thus, the fatigue life of the shaft is expected to be infinite. With respect to the spectral map of the rotor-bearing system, it can be concluded that the operating speed range is free from undesirable vibration amplitudes considering all excitation frequencies of interest. A good agreement among the results could be verified in all numerical analyses, the rotordynamic behavior of the developed single spool gas turbine was predicted, and potential vibration problems were avoided. The finite element analyses were strictly evaluated and the mechanical drawings of the aeronautical gas turbine were released to manufacturing. The first prototype is in its final experimental test phase and indicates that the performance on the test rig is consistent with the behavior indicated by the rotordynamic numerical analysis presented in this paper.



---

*Source: 4312943-2017-02-19.xml*