Open AccessArticle

A Prediction Method for the Damping Effect of Ring Dampers Applied to Thin-Walled Gears Based on Energy Method

Yanrong Wang

^1,2,

Hang Ye

^1,2,

Xianghua Jiang

^1,2 and

Aimei Tian

^3,*

School of Energy and Power Engineering, Beihang University, Beijing 100191, China

Collaborative-Innovation Center for Advanced Aero-Engine, Beijing 100191, China

School of Astronautics, Beihang University, Beijing 100191, China

Author to whom correspondence should be addressed.

Symmetry 2018, 10(12), 677; https://doi.org/10.3390/sym10120677

Submission received: 16 November 2018 / Revised: 26 November 2018 / Accepted: 30 November 2018 / Published: 30 November 2018

(This article belongs to the Special Issue Symmetry in Engineering Sciences)

Download

Browse Figures

Figure 1
An example of a gear with a ring damper. "> Figure 2
Finite element model: (a) The gear (fundamental sector); (b) the ring damper. "> Figure 3
Typical gear resonance failure [<a href="#B1-symmetry-10-00677" class="html-bibr">1</a>]. "> Figure 4
Mode shape of the gear with 3 nodal diameters. "> Figure 5
Local behavior in the contact region. "> Figure 6
Validation of the proposed method by [<a href="#B34-symmetry-10-00677" class="html-bibr">34</a>] results: Resonance amplitude by normalized normal pressure. "> Figure 7
The critical angle versus the normalized amplitude. "> Figure 8
Energy dissipated per cycle by the ring damper and maximum kinetic energy of the system versus normalized amplitude. "> Figure 9
Normalized frictional force and contact state: (a) B < Bc; (b) B = Bc; (c) B > Bc; (d) B ≫ Bc. "> Figure 10
Effect of the rotating speed: (a) Friction damping at various rotating speed; (b) friction damping for normalized rotating speed (for a given vibration stress). "> Figure 11
Effect of temperature. "> Figure 12
Effect of the ring damper density. "> Figure 13
Effect of the friction coefficient:(a) Friction damping at various friction coefficient; (b) friction damping for normalized friction coefficient (for a given vibration stress). "> Figure 14
Effect of the radial thickness of the ring damper. "> Figure 15
Effect of the axial thickness of the ring damper. ">

Review Reports Versions Notes

Abstract

In turbomachinery applications, thin-walled gears are cyclic symmetric structures and often subject to dynamic meshing loading which may result in high cycle fatigue (HCF) of the thin-walled gear. To avoid HCF failure, ring dampers are designed for gears to increase damping and reduce resonance amplitude. Ring dampers are installed in the groove. They are held in contact with the groove by normal pressure generated by interference or centrifugal force. Vibration energy is attenuated (converted to heat) by frictional force on the contact interface when the relative motion between ring dampers and gears takes place. In this article, a numerical method for the prediction of friction damping in thin-walled gears with ring dampers is proposed. The nonlinear damping due to the friction is expressed as equivalent mechanical damping in the form of vibration stress dependence. This method avoids the forced response analysis of nonlinear structures, thereby significantly reducing the time required for calculation. The validity of this numerical method is examined by a comparison with literature data. The method is applied to a thin-walled gear with a ring damper and the effect of design parameters on friction damping is studied. It is shown that the rotating speed, geometric size of ring dampers and friction coefficient significantly influence the damping performance.

Keywords:

thin-walled gear; ring damper; vibration; energy dissipation; friction damping

1. Introduction

Vibrations of gears are mainly caused by dynamic meshing loads. Resonance of the gear may occur if the excitation frequency is close to the resonance frequencies of the gear within its range of operating speeds. To avoid fatigue failure owing to high resonance stresses, the ideal solution is to redesign the gear to move its natural frequencies away from any potential external excitation. This method is called detuning [1]. However, for the thin-walled Gear, which is typically lightweight and operates at high rotating speed, detuning may not be feasible because each gear has multiple natural frequencies in coincidence with the mesh frequency within its operating range.

If detuning does not prevent resonance, then damping, as a passive control technique, is a feasible option to avoid high cycle fatigue failures. Friction dampers are effective approaches to provide damping in turbomachinery [2,3]. Friction dampers are substructures that remain in contact with the main structure through elastic deformation or centrifugal force. The vibration energy of the system is attenuated (converted to heat) by friction on the interface via the relative motion between the damper and primary structure [4].

Thin-walled structural components in aircraft gas turbine engines are easily excited to high vibration level. To reduce the vibrational stress of turbomachinery blades caused by the forced response from aerodynamic exciting sources and negative aerodynamic damping, i.e., flutter [5,6,7,8], many types of friction dampers have been studied and applied in actual structures. Among them, the under-platform damper has been extensively studied in detail [9,10,11,12,13,14]. This type of damper is installed under the platform or between neighboring blades. However, gears do not have suitable positions to install the under-platform damper. Therefore, ring dampers are used as damping devices for gears. In contrast with under-platform dampers, limited work has been carried out to investigate ring dampers. Lopez [15,16] used ring dampers on the train wheels to reduce the vibration emitted by freight traffic. The results revealed that increasing the mass of the ring damper is beneficial to vibration reduction. Laxalde [17] studied the damping strategy of ring dampers by using the dynamic Lagrangian frequency-time method to derive the forced response of blisks in the presence of ring dampers. The results showed that the size of the alternating stick-slip area determines the damping effectiveness of ring dampers. A nonlinear modal analysis method is proposed by Laxalde [18], and applied to analyze the effect of design parameters of ring dampers. Zucca [19] studied the effect of the key parameters (for example, mass and friction coefficient) of ring dampers on the vibration amplitude. The authors used the contact element to link the static and dynamic differential equations and calculated the forced response of the coupling system. Tang [20] proposed a novel reduced-order modeling method to solve the forced responses of the blisk–damper systems based on Craig–Bampton component mode synthesis. The authors studied the effect of geometric parameters of ring dampers on the blisk forced responses [21] by this method.

For ring dampers to be effective, they are typically located on the rim of the gear where large vibration amplitudes occur, as shown in Figure 1. Otherwise the energy dissipation due to friction will be reduced and even equal to zero, and the ring damper will be ineffective. Ring dampers are mostly effective only for the fundamental mode shapes of the gear [22]. These modes are characterized by a large amplitude at the rim of the gear. For thin-walled gears, friction damping is produced by the relative motion caused by the different extension deformations between ring dampers and gears along the tangential direction of the contact surface [23]. However, note that the circumferential deformation is caused by radial vibration. In other applications, for example train wheels, vibration energy is attenuated by the axial component of the vibration, and friction damping is produced by relative motion in the axial direction [15]. Zucca [22] analyzed the axial and circumferential relative motion of a bevel gear with a ring damper in different response conditions. The results show that although the radial and axial components of the vibration have the same order of magnitude, the ring damper worked mainly in the circumferential direction because the relative displacement along the circumferential direction is much larger than along the axial direction. No relative motion occurs in the radial direction due to the ring damper maintaining contact with the primary structure by centrifugal force.

Although all of these papers show that vibration amplitude will decrease when ring dampers are used, limited work to investigate the nonlinear friction damping of thin-walled gears with ring dampers has been done. Most previous theoretical analyses have focused on the forced response of main structures in the presence of ring dampers. In contrast, the energy dissipation by ring dampers has been seldom studied. Niemotka [24] proposed a design method for split ring dampers to lower the vibration amplitude of annular air seals in gas turbine engines based on a quasi-static energy dissipation analysis.

The primary objective of this work is to construct a numerical model to predict the damping of ring dampers in thin-walled gears. In the model, the nonlinear friction damping is expressed as equivalent mechanical damping in the form of vibration stress dependent. Macro-slip is used in the friction model to calculate the energy dissipation. The validity of the proposed method is confirmed by a comparison with forced response analysis results. The secondary objective is to investigate the influence of rotating speed, temperature, parameters of ring dampers, and friction coefficient on the damping performance by means of method proposed in this paper.

The rest of this paper is arranged as follows. The theoretical background, including the equation of motion and modal analysis, is introduced in Section 2. Theoretical derivation of equivalent damping ratio of the ring damper is shown in Section 3. Method validation and parameter analysis are performed on a thin-walled gear in Section 4, followed by conclusions in Section 5.

2. Vibration Analysis of The Gear-Ring Damper System

2.1. The Equations of Motion

The equations of motion in time domain of the gear-ring damper system can be written as

M \ddot{X} + C \dot{X} + KX + F_{n l} (X, \dot{X}, t) = F (t)

(1)

where M, C, and K are the mass, damping, and stiffness matrices of the gear, respectively, and X is the vector of the displacements. F(t) is the vector of the external excitation force.

F_{n l} (X, \dot{X}, t)

is the vector of the nonlinear forces generated by the ring damper and depends on the vibration displacement and vibration velocity of the system.

F_{n l} (X, \dot{X}, t)

can be given by the equivalent damping and stiffness matrices as [25]

F_{n l} (X, \dot{X}, t) {= C}_{eq} \dot{X} + K_{eq} X

(2)

The equivalent damping matrix C_eq and the equivalent stiffness matrix K_eq depend on the motion of the gear.

The displacement vector X is a function of time and can be expressed as a linear combination of the natural modes of the un-damped system.

X (t) = Φ q (t)

(3)

Thus, Equation (1) can be rewritten as

M Φ \ddot{q} (t) + C Φ \dot{q} (t) + K Φ q (t) + C_{eq} Φ \dot{q} (t) + K_{eq} Φ q (t) = F (t)

(4)

where Φ is the mass-normalized eigenvector matrix of the gear.

Premultiplying Equation (4) throughout by Φ^T:

I \ddot{q} (t) + Z \dot{q} (t) + Λ q (t) + Z_{eq} \dot{q} (t) + Λ_{eq} q (t) = Q (t)

(5)

where

I = Φ^TMΦ, Z = Φ^TCΦ, Λ = Φ^TKΦ, Z_eq = Φ^TC_eqΦ, Λ_eq = Φ^TΛ_eq, Q(t) = Φ^TF(t)

(6)

because I denotes the unity matrix and Z, Λ, Z_eq, and Λ_eq are all diagonal. In the vicinity of the jth natural frequency, Equation (5) can be rewritten as

{\ddot{q}}_{j} (t) + 2 (ζ_{j} + ζ_{j, eq}) ω_{j} {\dot{q}}_{j} (t) + (k_{j} + k_{j, eq}) q_{j} (t) = Q_{j} (t), with j = 1, 2, \dots, n

(7)

where ζ_j and ζ_j,_eq are the modal damping ratio and the equivalent damping ratio caused by the ring damper for the jth mode, respectively; k_j and k_j,_eq are the modal stiffness and equivalent stiffness for the jth mode, respectively; and

k_{j} = ω_{j}^{2}

; ω_j is the jth natural frequency of the undamped system.

The n equations represented in Equation (7) can be uncoupled from all other equations. Therefore, the forced response of the jth mode can be calculated if the relationship between the equivalent damping and the equivalent stiffness and response amplitude can be pre-calculated.

In general, the mass of the ring damper is much smaller than the mass of the main structure. Let the weight penalty be defined as

β = \frac{mass of the ring damper}{mass of the gear}

(8)

In this study, the weight penalty is less than 5%. Note that the magnitudes of M and K are much larger than the magnitude of

F_{n l} (X, \dot{X}, t)

, thus k_j,_eq is much smaller than k_j. Generally, for the ring damper, k_j,_eq is two orders of magnitude lower than k_j. In other words, the ring damper does not affect the shape of the vibration mode; rather, it affects only the vibration amplitude. Moreover, the influence of the damper on the resonance frequency of the primary structure can be neglected. However, the equivalent damping matrix is of the same order of magnitude or even larger with respect to the damping matrix because the structural damping is usually small (For steel, the damping ratio is 1~5×10⁻⁴). The results of other scholars [2,3,20,25,26,27] also showed that the influence of the ring damper on the frequency is negligible. With or without ring dampers, the frequency variation is less than 1%. Thus, the damper ring reduces the resonant amplitude of the gear, primarily by providing damping, rather than changing the stiffness of the gear system.

2.2. Modal Analysis

Modal analysis was performed with the FEM software ANSYS 14.5. The gear and ring damper finite element models are shown in Figure 2. The gear is a cyclic symmetry structure, comprising z fundamental sectors (Figure 2a). The ring damper is machined to be C-shaped for ease of installation. There is a split in the axial direction, as shown in Figure 2b.

Typical gear resonance failure in practice [1] is shown in Figure 3. The mode shapes (Figure 4) that lead to gears failure have the following features:

The modal amplitude has an integer number of harmonic distributions along the circumferential direction.
The nodal line passes through the center of rotation, and the vibration amplitude of the nodal line is zero.
For thin-walled gears, the gear rim vibrates mainly in the radial direction.

Therefore, in this study, we focused on nodal diameter vibration. For N nodal diameters (ND), the radial displacement of the groove of the gear can be assumed as

w = B \cos (N θ)

(9)

where B is the maximum amplitude of the groove of the gear, N is the number of nodal diameters, θ is circumferential angle.

3. Theoretical Model of Equivalent Damping Ratio of The Ring Damper

3.1. Energy Dissipated by Frictional Force

In this paper, the motion of the gear is assumed to be small amplitude vibrations, i.e., only elastic deformation is considered and in the same mode shape, and the vibration stress is proportional to the vibration amplitude. The following energy dissipation analysis is based on the method proposed by Alford [28,29,30,31] and Niemotka [24].

Generally, deflections of the structure at a resonance are very small compared to its size; otherwise, the structure will suffer fatigue failure in a short time. For small deformations, the strain-curvature relation is

ε = κ y

(10)

where ε, κ, and y are strain, curvature, and the distance from the neutral line, respectively.

The preceding equation shows that the circumferential strains are proportional to the curvature and are linearly related with the distance y from the neutral line. Here tensile strain is defined as positive and compressive strain is defined as negative.

The curvature can be expressed by the bending moment:

κ = \frac{M}{E I}

(11)

where M, E, and I are bending moment, Young’s modulus, and moment of inertia of ring dampers, respectively. Equation (11) is known as the moment-curvature equation. When the radius of curvature of a ring is sufficiently large compared to its radial height, the relationship between the bending moment M and radial displacement w can be expressed as [32]

\frac{M}{E I} = \frac{1}{R^{2}} [w + \frac{d^{2} w}{d θ^{2}}]

(12)

By substituting Equation (9) into Equation (12), the following relationship is obtained:

\frac{M}{E I} = \frac{1}{R^{2}} B (1 - N^{2}) \cos (N θ)

(13)

At a distance y from the mean radius R, the bending strain is:

ε_{y} = \frac{y}{R^{2}} B (1 - N^{2}) \cos (N θ)

(14)

For the gear, the strain on the contact surface of the gear is tensile on the groove interface; in contrast, it is compressive for the ring damper on the contact surface and vice versa, as shown in Figure 5.

ε_{g} = - \frac{c_{g}}{R_{g}^{2}} B (1 - N^{2}) \cos (N θ)

(15)

ε_{d} = \frac{c_{d}}{R_{d}^{2}} B (1 - N^{2}) \cos (N θ)

(16)

where c and R are the half of the radial thickness and the radius. Subscript g and d represent gear and damper respectively.

When there is no relative motion on the contact surface, the contact state is the stick state. The relationship between strain caused by friction and bending strain is

ε_{f} = ε_{g} - ε_{d} = - (\frac{c_{d}}{R_{d}^{2}} + \frac{c_{g}}{R_{g}^{2}}) B (1 - N^{2}) \cos (N θ)

(17)

The strain caused by friction in the ring damper also can be calculated by dividing the frictional force by the product of the damper cross-sectional area and its Young’s modulus. F_f is defined as the frictional force per unit length, where F_f is a function of circumferential angle θ.

ε_{f} = \frac{R_{d}}{A_{d} E} \int F_{f} d θ

(18)

By substituting Equation (17) into Equation (18), F_f can be written as

F_{f} = \frac{A_{d} E}{R_{d}} \frac{d ε_{f}}{d θ} = - \frac{B A_{d} E}{R_{d}} (\frac{c_{d}}{R_{d}^{2}} + \frac{c_{g}}{R_{g}^{2}}) N (1 - N^{2}) \sin (N θ)

(19)

When no slipping occurs on the entire contact surface, F_{f max} will appears at θ = π/2N. And

F_{f \max} = \frac{B A_{d} E}{R_{d}} (\frac{c_{d}}{R_{d}^{2}} + \frac{c_{g}}{R_{g}^{2}}) N (1 - N^{2})

. When tangential force is greater than the maximum static friction, slipping occurs at θ < π/2N, and over the zone θ₀ < θ < π/2N, F_{f max} = μP. Where μ is friction coefficient, and P is normal pressure on the contact surface.

At θ_{0} = θ, F_{f \max} = μ P = - \frac{B A_{d} E}{R_{d}} (\frac{c_{d}}{R_{d}^{2}} + \frac{c_{g}}{R_{g}^{2}}) N (1 - N^{2}) \sin (N θ_{0})

(20)

Thus,

N θ_{0} = \arcsin \frac{μ P}{\frac{B E A_{d}}{R_{d}} (\frac{c_{g}}{R_{g}^{2}} + \frac{c_{d}}{R_{d}^{2}}) N (N^{2} - 1)}

(21)

where θ₀ represents the angle where slippage starts, which is called the critical slip angle.

When normal pressure P is constant, over a vibration cycle, the condition that no slipping occurs on the entire contact surface is the maximum vibration amplitude of the gear B is less than the critical vibration amplitude B_c.

B_{c} = \frac{μ P R_{d}}{E (\frac{c_{g}}{R_{g}^{2}} + \frac{c_{d}}{R_{d}^{2}})} \frac{1}{N (N^{2} - 1)}

(22)

In the sliding zone, the frictional force is equal to sliding frictional force μP. Therefore, the strain caused by friction can be written as

ε_{f} = - \frac{μ P R_{d}}{E A_{d}} (θ - \frac{π}{2 N})

(23)

where R_d and A_d are respectively radius and the cross-sectional area of the ring damper.

The relative displacement on the contact surface can be obtained by integrating the strain. Note that displacement is 0 at the beginning of the sliding zone.

{\begin{array}{l} s (θ) = 0, & 0 \leq θ \leq θ_{0} \\ s (θ) = \int_{θ_{0}}^{θ} (ε_{g} - ε_{d} - ε_{f}) R_{f} d θ, & θ_{0} \leq θ \leq \frac{π}{2 N} \end{array}

(24)

Therefore,

\begin{array}{l} s (θ) & = \int_{θ_{0}}^{θ} (ε_{g} - ε_{d} - ε_{f}) R_{d} d θ \\ = \frac{μ P R_{d} R_{f}}{A_{d} E} [- \frac{1}{N^{2}} \frac{R_{d}}{R_{g}} (\frac{\sin (N θ)}{\sin (N θ_{0})} - 1) + \frac{1}{2} [{(\frac{π}{2 N} - θ)}^{2} - {(\frac{π}{2 N} - θ_{0})}^{2}]] \end{array}

(25)

The energy dissipated by the ring damper in a complete vibration cycle, ΔW, can be obtained by integrating the product of the frictional force F_f and the relative displacement s(θ) in the slip region.

\begin{array}{l} Δ W = 16 N \int_{θ_{0}}^{\frac{π}{2 N}} F_{f} \cdot Δ s (θ) R_{f} d θ = 16 \frac{{(μ P)}^{2} R_{f}^{3}}{N^{2} E A_{r}} {[\cot (N θ_{0}) + N θ_{0} - \frac{π}{2}] - \frac{1}{3} {(\frac{π}{2} - N θ_{0})}^{3}} \end{array}

(26)

Note that ΔW depends on the critical slip angle θ₀. According to Equation (21), θ₀ is a nonlinear function of B. Therefore, ΔW is a function of B.

3.2. Equivalent Damping Ratio

The loss coefficient η or damping ratio ζ is commonly used to indicate the damping capacity of engineering structures. The loss coefficient η is defined as the ratio of the energy dissipated per radian and the total vibration energy [33]:

η = \frac{Δ W / 2 π}{W} ≃ 2 ζ

(27)

For small damping, the total vibration energy of the system W approximately equal to the maximum kinetic energy [33]. Thus, the total vibration energy for the jth normal mode can be expressed as

W = \frac{1}{2} {[\dot{X}]}^{T} [M] [\dot{X}] = \frac{1}{2} ω_{j}^{2} q_{j}^{2}

(28)

Thus, the equivalent structural damping ratio k_j,_eq in Equation (7) can be rewritten as

ζ_{j, eq} = \frac{Δ W}{W} = \frac{16 \frac{{(μ P)}^{2} R_{f}^{3}}{N^{2} E A_{r}} {[\cot (N θ_{0}) + N θ_{0} - \frac{π}{2}] - \frac{1}{3} {(\frac{π}{2} - N θ_{0})}^{3}}}{ω_{j}^{2} q_{j}^{2} / 2}

(29)

4. Application and Discussion

To validate the method shown in this article, the numerical simulation is applied to a real thin-walled gear made of 4310 steel (Young’s modulus E = 207 GPa and density ρ = 7.84 × 10³ kg/m³). The mass of the gear is 425 g. Figure 4 shows the mode shape of the model with 3 ND. The corresponding natural frequency is 3758 Hz. For reasons of confidentiality, some of results are given in a normalized form.

4.1. Method Validation

The influence of the normal pressure on damping effect is compared with the results from the forced response analysis based on the harmonic balance method in [34], as shown in Figure 6. The results obtained by the two analysis methods are highly consistent. However, the method shown in this article does not need to calculate the equation of motion in the frequency domain or time domain, so it has faster calculation speed. Since the numerical method shown in this article is independent of excitation and inherent mechanical damping, the excitation and mechanical damping are given in accordance with [34]. Since the normal pressure is not directly given in [34], the normal pressure in this section is a relative value (defined as normalized normal pressure P’).

At P’ = 0, the frictional force at the contact surfaces is 0, and the ring damper can freely slide relative to the gear. The energy dissipated by frictional is 0, and the ring dampers is ineffective. An increment of P’ leads to the vibration to decrease down to a minimum value, corresponding to the optimum normalized normal pressure (about 0.45). A further increment of P’ causes the vibration to increase again. When P’ is large enough (about 1.65), the vibration amplitude increases to the amplitude at P’= 0. In this case, no relative motion takes place on the contact surface of the two structures Thus, the ring damper ceased to be effective. It is worth mentioning that two different analysis methods show that when the normal pressure is greater than about 3.7 times of the optimal normal pressure, the ring damper ceased to be effective, which will be further explained in the following parameter sensitivity analyses.

For a given normal pressure (or rotating speed), when the vibration amplitude B is small, the ring damper is full-stick, and there is no slip, as shown in Figure 7. When B increases to the critical vibration amplitude B_c, sliding appears in θ₀ = π/2N. When B increases, the critical slip angle decreases and the slip area increases. When the vibration amplitude is large enough, the critical slip angle approaches 0, and the ring damper is approximately full-slip. In this case, the energy dissipation caused by the ring damper is approximately linear with the vibration amplitude, as shown in Figure 8.

The normalized frictional force and the contact state between the gear and the ring damper along the circumferential direction are shown in Figure 9. In Figure 9a, when the vibration amplitude B is less than B_c, the contact state is stick. Thus, no relative motion occurs on the contact surface. Frictional force is a function of θ, and the maximum frictional force appears at the position of the nodal line. When B = B_c, slip appears at the position of the nodal line, as shown in Figure 9b. When B > B_c, the slip region expands to both sides as B increases, as shown in Figure 9c. When B ≫ B_c, the slip region increases slowly as the vibration amplitude increases. In this case, the contact status is approximately full-slip, as shown in Figure 9d. This observation is highly consistent with other studies [9,20], according to those studies, when the excitation frequency is far from the natural frequency, the response amplitude is small and the contact status is stick. When the excitation frequency gradually approaches the natural frequency, relative slip appears on the contact surface and the slip region gradually increases.

4.2. Effect of Ring Damper Parameters

4.2.1. Effect of Rotating Speed or Normal Pressure

The normal pressure on the contact surface depends on the rotating speed of the system, and the normal pressure is proportional to the square of the rotational speed. Thus, in this article, only the effect of rotating speed is shown.

Figure 10 shows the effect of rotating speed on the equivalent damping performance. In Figure 10a, a decrement of the rotating speed causes the contact surface to slide more easily at a given resonance stress, resulting in a lower critical vibration stress. Also, the vibration stress corresponding to the maximum damping ratio decreases as the rotating speed decreases. In Figure 10b, for a given resonance stress, when the rotating speed is greater than about 1.9 times of the optimal rotating speed, the contact surface is full-stick, where the optimum rotating speed is defined as the rotating speed corresponding to the maximum damping ratio.

4.2.2. Effect of Temperature

Figure 11 shows the effect of temperature on the damping performance. The effect of temperature is negligible. This, of course, is because the change in Young’s modulus E is small during the operating temperature range. The stiffness of the ring damper is almost unchanged. This also indicates that the ring damper can work at high temperatures and with good temperature adaptability.

4.2.3. Effect of the Ring Damper Density

The effect of the density is investigated according to its effect on the normal pressure acting on the contact surface. The normal direction is defined as along the radial direction of the gear.

The effect of the ring damper density on the damping performance is shown in Figure 12. The critical vibration stress increases with an increase of the ring damper density. If the density is too large, then the ring damper ceases to be effective due to the contact surface tends to stick. In this case no energy is dissipated by frictional force.

4.2.4. Effect of the Friction Coefficient

As shown in Figure 13a, the effect of the friction coefficient μ on the damping performance is similar to density. Increasing μ results in an increase in critical vibration stress. Moreover, in this case, the contact surface tends to be full-stick due to an increase in μ. In contrast, a decrease in μ results in the contact surface tending to be full-slip. However, due to F_{f max} = μP, the maximum frictional force on the contact surface F_{f max} decreased with a decrease in μ. For a given vibration stress, there is an optimum density that maximizes frictional damping. When the density is greater than 3.7 times the optimal density, the ring damper will cease to be effective again, as shown in Figure 13b.

4.2.5. Effect of the Cross-Sectional Area of the Ring Damper

The cross-sectional area is equal to the product of the radial thickness and the axial thickness of the ring damper. The effect of the radial thickness is shown in Figure 14. The critical vibration stress decreases and the peak damping ratio increased with an increase in the radial thickness. Increasing the radial thickness can significantly improve the damping performance. The effect of axial thickness is shown in Figure 15. The critical vibration stress is not affected by the axial thickness. However, the peak damping ratios increase with an increase in the axial thickness. In the premise that the mass of ring dampers is much smaller than the mass of gears, the equivalent damping ratio is approximately linear with the axial thickness.

Therefore, for a given cross-sectional area, the ring damper with a large ratio of radial thickness to axial thickness has a better damping effect.

5. Conclusions

In this article, a theoretical study of ring dampers for thin-walled gears has been shown. A numerical method to predict the damping performance of ring dampers is proposed. In the proposed method, the energy dissipated by the ring damper is calculated through a quasi-static process then it is expressed as the equivalent mechanical damping function that depends on vibration stress. The validity of the model is confirmed by a comparison with forced response analysis results. Compared with forced response analysis, the method shown in this article only needs once modal analysis of the primary structure. The proposed method avoids computation of the periodical response of the nonlinear structure. Therefore, minimal computation is required to obtain the damping performance, which greatly improves the efficiency of ring dampers design.

The damping performance of the ring damper depends on the vibration amplitude of the gear B and the damper parameters. When B is less than the critical vibration amplitude B_c, the ring damper is ineffective. When B is greater than B_c, the ring damper can provide friction damping. By increasing B, slip first appears at position of the nodal line, and the slip region expands to both sides as B increases. At approximately 3.7 times the critical vibration amplitude, the efficiency of the damper is theoretically maximized.

For a given amplitude, there is optimum speed, density and friction coefficient to maximize damping. Excessively increasing or decreasing the rotating speed, the ring damper density and the friction coefficient will cause the contact surface to be full-stick or full-slide. In both cases, the ring damper does not provide frictional damping. For a given mass of ring dampers, different damping performances may be observed if the density and the ratio of radial thickness to axial thickness are different.

The proposed method works well when the mass of the ring damper is significantly less than the mass of the primary structure. The ring damper can provide substantial damping and only weakly affects the mode shape of the system. This methodology is suitable for specific applications such as gears or blisks with ring dampers.

Author Contributions

Conceptualization and Methodology, Y.W. and H.Y.; Software development, parameter analysis, and Writing—Original Draft Preparation, H.Y.; Review and Editing, X.J.; Funding acquisition, A.M.

Funding

This work was supported by the National Nature Science Foundation of China (No. 51475022).

Acknowledgments

The authors would like to thank all of the reviewers for their constructive comments. In addition, Hang Ye especially wishes to thank Yiheng Zhang for her continuous support.

Conflicts of Interest

The authors declare no conflict of interest.

Notation

B	vibration amplitude	Subscript g	gear
B_c	critical vibration amplitude	Subscript d	ring damper
C	damping matrices of the gear	Subscript eq	equivalent
c	half-width of the gear rim or the ring damper	W	total energy of the system
E	Young’s modulus	w	radial displacement of the groove
F(t)	external periodic force	X	displacement vector
$F_{n l} (X, \dot{X}, t)$	nonlinear frictional force	z	number of teeth of the gear
F_f	frictional force per unit length	ε	strain
I	sectional moment of inertia	η	loss coefficient
K	stiffness matrices of the gear	κ	curvature
M	mass matrices of the gear	μ	friction coefficient
M	bending moment	θ	circumferential angle
N	number of nodal diameters	θ₀	critical slip angle
P	normal pressure	ρ	density of the ring damper
P’	normalized normal pressure	ζ	damping ratio
R	radius	ζ_eq	equivalent damping ratio provided by the ring damper
s	relative displacement	ΔW	energy dissipated per cycle by the ring damper

References

Drago, R.J.; Brown, F.W. The analytical and experimental evaluation of resonant response in high-speed, lightweight, highly loaded gearing. J. Mech. Des. 1981, 103, 346–356. [Google Scholar] [CrossRef]
Griffin, J.H. Friction damping of resonant stresses in gas turbine engine airfoils. J. Eng. Power 1980, 102, 329–333. [Google Scholar] [CrossRef]
Cameron, T.M.; Griffin, J.H.; Kielb, R.E. An integrated approach for friction damper design. J. Vib. Acoust. 1990, 112, 175–182. [Google Scholar] [CrossRef]
Ferri, A.A. Friction damping and isolation systems. J. Vib. Acoust. 1995, 117, 196–206. [Google Scholar] [CrossRef]
Wang, Y.; Fu, Z.; Jiang, X.; Tian, A. Mistuning effects on aero-elastic stability of axial compressor rotor blades. J. Eng. Gas Turbines Power 2015, 137, 102504. [Google Scholar] [CrossRef]
Fu, Z.; Wang, Y.; Jiang, X.; Wei, D. Tip Clearance effects on aero-elastic stability of axial compressor blades. J. Eng. Gas Turbines Power 2015, 137, 012501. [Google Scholar] [CrossRef]
Zhang, X.; Wang, Y.; Xu, K. Mechanisms and key parameters for compressor blade stall flutter. J. Turbomach. 2013, 135, 024501. [Google Scholar] [CrossRef]
Zhang, X.; Wang, Y.; Xu, K. Flutter prediction in turbomachinery with energy method. J. Aerosp. Eng. 2011, 225, 995–1002. [Google Scholar] [CrossRef]
Koh, K.H.; Griffin, J.H.; Filippi, S.; Akay, A. Characterization of turbine blade friction dampers. J. Eng. Gas Turbines Power 2005, 127, 856–862. [Google Scholar] [CrossRef]
Petrov, E.P.; Ewins, D.J. Advanced modeling of underplatform friction dampers for analysis of bladed disk vibration. J. Turbomach. 2007, 129, 143–150. [Google Scholar] [CrossRef]
Botto, D.; Gastadi, C.; Gola, M.M.; Umer, M. An experimental investigation of the dynamics of a blade with two under-platform dampers. J. Eng. Gas Turbines Power 2017, 140, 032504. [Google Scholar] [CrossRef]
Botto, D.; Umer, M. A novel test rig to investigate under-platform damper dynamics. Mech. Syst. Signal. Pr. 2018, 100. [Google Scholar] [CrossRef]
Pesaresia, L.; Sallesa, L.; Jonesb, A.; Greenb, J.S.; Schwingshackla, C.W. Modelling the nonlinear behaviour of an underplatform damper test rig for turbine applications. Mech. Syst. Signal. Pr. 2017. [Google Scholar] [CrossRef]
Kaneko, Y. Vibration response analysis of mistuned bladed disk with under-platform damper: Effect of variation of contact condition on vibration characteristics. In Proceedings of the ASME Turbo Expo 2017: Turbomachinery Technical Conference and Exposition, Charlotte, NC, USA, 26–30 June 2017. [Google Scholar]
Lopez, I.; Busturia, J.M.; Nijmeijer, H. Energy dissipation of a friction damper. J. Sound Vib. 2004, 278, 539–561. [Google Scholar] [CrossRef]
Lopez, I.; Nijmeijer, H. Prediction and validation of the energy dissipation of a friction damper. J. Sound Vib. 2009, 328, 396–410. [Google Scholar] [CrossRef]
Laxalde, D.; Thouverez, F.; Lombard, J.P. Forced response analysis of integrally bladed disks with friction ring dampers. J. Vib. Acoust. 2010, 132, 011013. [Google Scholar] [CrossRef]
Laxalde, D.; Salles, L.; Blanc, L.; Thouverez, F. Non-linear modal analysis for bladed disks with friction contact interfaces. In Proceedings of the ASME Turbo Expo 2008: Power for Land, Sea, and Air. American Society of Mechanical Engineers, Berlin, Germany, 9–13 June 2008. [Google Scholar]
Zucca, S.; Firrone, C.M.; Facchini, M. A method for the design of ring dampers for gears in aeronautical applications. J. Mech. Des. 2012, 134, 091003. [Google Scholar] [CrossRef]
Tang, W.; Epureanu, B.I. Geometric optimization of dry friction ring dampers. Int. J. Nonlin. Mech. 2018. [Google Scholar] [CrossRef]
Tang, W.; Epureanu, B.I. Nonlinear dynamics of mistuned bladed disks with ring dampers, International. Int. J. Nonlin. Mech. 2017, 97, 30–40. [Google Scholar] [CrossRef]
Firrone, C.M.; Zucca, S. Passive control of vibration of thin-walled gears: Advanced modelling of ring dampers. Nonlinear Dyn. 2014, 76, 263–280. [Google Scholar] [CrossRef]
Buyukataman, K. A theoretical study on the vibration damping of aircraft gearbox gears. In Proceedings of the 27th Joint Propulsion Conference, Sacramento, CA, USA, 24–26 June 1991. [Google Scholar]
Niemotka, M.A.; Ziegert, J.C. Optimal design of split ring dampers for gas turbine engines. J. Eng. Gas Turbines Power 1995, 117, 569–575. [Google Scholar] [CrossRef]
Baek, S.; Epureanu, B. Reduced-order modeling of bladed disks with friction ring dampers. J. Vib. Acoust. 2017, 139, 061011. [Google Scholar] [CrossRef]
Tangpong, X.W.; Wickert, J.A.; Akay, A. Finite element model for hysteretic friction damping of traveling wave vibration in axisymmetric structures. J. Vib. Acoust. 2008, 130, 011005. [Google Scholar] [CrossRef]
Tangpong, X.W.; Wickert, J.A.; Akay, A. Distributed friction damping of travelling wave vibration in rods. Philos. Trans. R. Soc. A 2008, 366, 811–827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Alford, J.S. Protection of labyrinth seals from flexural vibration. J. Eng. Power 1964, 86, 141–147. [Google Scholar] [CrossRef]
Alford, J.S. Protecting turbomachinery from self-excited rotor whirl. J. Eng. Power 1965, 87, 333–343. [Google Scholar] [CrossRef]
Alford, J.S. Protecting turbomachinery from unstable and oscillatory flows. J. Eng. Power 1967, 89, 513–527. [Google Scholar] [CrossRef]
Alford, J.S. Nature, causes, and prevention of labyrinth air seal failures. J. Aircraft 1975, 12, 313–318. [Google Scholar] [CrossRef]
Weaver, J.W.; Timoshenko, S.P.; Young, D.H. Vibration Problems in Engineering, 5th ed.; Wiley-Interscience: Hoboken, NJ, USA, 1990; ISBN 9780471632283. [Google Scholar]
Rao, S.S. Mechanical Vibrations, 5th ed.; Prentice Hall: New York, NY, USA, 1995; ISBN 978-0-13-212819-3. [Google Scholar]
Tang, W.; Baek, S.; Epureanu, B.I. Reduced-order models for blisks with small and large mistuning and friction dampers. J. Eng. Gas Turbines Power 2017, 139, 012507. [Google Scholar] [CrossRef]

Figure 1. An example of a gear with a ring damper.

Figure 2. Finite element model: (a) The gear (fundamental sector); (b) the ring damper.

Figure 3. Typical gear resonance failure [1].

Figure 4. Mode shape of the gear with 3 nodal diameters.

Figure 5. Local behavior in the contact region.

Figure 6. Validation of the proposed method by [34] results: Resonance amplitude by normalized normal pressure.

Figure 7. The critical angle versus the normalized amplitude.

Figure 8. Energy dissipated per cycle by the ring damper and maximum kinetic energy of the system versus normalized amplitude.

Figure 9. Normalized frictional force and contact state: (a) B < B_c; (b) B = B_c; (c) B > B_c; (d) B ≫ B_c.

Figure 10. Effect of the rotating speed: (a) Friction damping at various rotating speed; (b) friction damping for normalized rotating speed (for a given vibration stress).

Figure 11. Effect of temperature.

Figure 12. Effect of the ring damper density.

Figure 13. Effect of the friction coefficient:(a) Friction damping at various friction coefficient; (b) friction damping for normalized friction coefficient (for a given vibration stress).

Figure 14. Effect of the radial thickness of the ring damper.

Figure 15. Effect of the axial thickness of the ring damper.

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Wang, Y.; Ye, H.; Jiang, X.; Tian, A. A Prediction Method for the Damping Effect of Ring Dampers Applied to Thin-Walled Gears Based on Energy Method. Symmetry 2018, 10, 677. https://doi.org/10.3390/sym10120677

AMA Style

Wang Y, Ye H, Jiang X, Tian A. A Prediction Method for the Damping Effect of Ring Dampers Applied to Thin-Walled Gears Based on Energy Method. Symmetry. 2018; 10(12):677. https://doi.org/10.3390/sym10120677

Chicago/Turabian Style

Wang, Yanrong, Hang Ye, Xianghua Jiang, and Aimei Tian. 2018. "A Prediction Method for the Damping Effect of Ring Dampers Applied to Thin-Walled Gears Based on Energy Method" Symmetry 10, no. 12: 677. https://doi.org/10.3390/sym10120677

APA Style

Wang, Y., Ye, H., Jiang, X., & Tian, A. (2018). A Prediction Method for the Damping Effect of Ring Dampers Applied to Thin-Walled Gears Based on Energy Method. Symmetry, 10(12), 677. https://doi.org/10.3390/sym10120677

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

Article Menu