Effect of mistuning parameters on dynamic characteristics of mistuned bladed disk

Hongyuan Zhang1 , Haocun Dong2 , Tianyu Zhao3 , Xin Zhang4 , Huiqun Yuan5

1, 2, 4School of Automotive and Transportation, Shenyang Ligong University, Shenyang, China

1School of Mechanical Engineering and Automation, Northeastern University, Shenyang, China

3, 5School of Sciences, Northeastern University, Shenyang, China

2Corresponding author

Journal of Vibroengineering, Vol. 23, Issue 4, 2021, p. 784-798. https://doi.org/10.21595/jve.2021.21629
Received 3 August 2020; received in revised form 1 February 2021; accepted 12 February 2021; published 2 March 2021

Copyright © 2021 Hongyuan Zhang, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 80
Reads 30
Downloads 292
WoS Core Citations 0
CrossRef Citations 0
Abstract.

The aim of this paper is to study the nonlinear dynamics of the tenon and mortise of aero-engine compressor blades, based on lumped parameter model of tenon and mortise of blade disk system, considering the influence of gap and friction, The nonlinear vibration characteristics of mistuned bladed disk with different mistuning parameters are studied. The effects of tenon and mortise gap, dry friction, blade damping mistuning and tenon damping mistuning on the vibration characteristics of the mistuned bladed disk under strong coupling and weak coupling are obtained respectively. The results show that the gap and friction has little effect on the vibration response of the weakly coupled bladed disk, but the gap and friction has a great influence on the strongly coupled bladed disk, the larger the gap, the more complex the nonlinear dynamic characteristics of tenon. Under the same gap and coupling strength, the mistuning parameters have a certain influence on the vibration amplitude of the bladed disk.

Effect of mistuning parameters on dynamic characteristics of mistuned bladed disk

Highlights
  • A lumped parameter model of blade disk system considering the influence of gap and friction is proposed.
  • The nonlinear vibration characteristics of mistuned bladed disk with different mistuning parameters are studied.
  • The effects of tenon and mortise gap, dry friction, blade damping mistuning and tenon damping mistuning on the vibration characteristics of the mistuned bladed disk under strong coupling and weak coupling are obtained respectively.

Keywords: mistuned bladed disk, tenon and mortise, mistuning parameter, gap, friction, coupling.

1. Introduction

Tenon and mortise structure has the advantages of large load-bearing, small structure size and convenient maintenance and replacement. Therefore, tenon and mortise structure is widely used in aero-engine compressor and other power mechanical bladed disk system. However, due to the complexity of the working environment, the tenon and mortise structure of bladed disk system is prone to failure, The nonlinear factors such as gap and friction between tenon and mortise make the vibration characteristics of bladed disk system more complex, especially for mistuned bladed disk system.

Many scholars have carried out a lot of research on the non-linear factors of mortise and tenon. Cha [1] proposed an analysis method based on the equivalent linearization method to calculate the statistical characteristics of the nonlinear response, and studied the relationship between the vibration mode of the system and the performance of the friction damper. Sinclair [2-4] found that the non-linearity of the tenon system is caused by the friction caused by the deformation of the tenon. Petrov [5] developed an effective method to analyze the forced response of nonlinear periodic symmetric structures, and used this method to analyze the tuned bladed disk system with friction damping and gap. Ender [6] used the lumped mass model to calculate the effect of different dry friction on the vibration of bladed disk system. Wang [7] established a mass spring model with dry friction by harmonic balance method, taking into account the effects of coupling strength, viscosity damping ratio, excitation force amplitude to friction amplitude ratio, and the degree of dry friction disharmony on the bladed disk response. Wang [8] established a practical model of bladed disk system, which considered the random mistuning of the delayed dry friction damping, and analyzed the effect of dry friction on the vibration localization of the random mistuning bladed disk system by the incremental harmonic balance method. The results show that the nonlinear friction damping also leads to the localized phenomenon of the tuned disk. Firrone [9] proposes an overlapping method, which combines the static and dynamic displacements of blade disk and damper to calculate the forced response of blade disk system. Zhang [10] used the micro sliding model of friction damper to study the blade root dry friction, and found that the blade root dry friction can effectively reduce the vibration amplitude of mistuned blade disk system. Ma [11] established a finite element model of tenon joint structure by using 1/38 model, and comprehensively considered the effects of centrifugal force, friction coefficient, installation angle and contact area on the contact characteristics of tenon and mortise. The results show that the friction coefficient can effectively affect the maximum contact pressure and sliding distance, but the installation angle of tenon has no obvious effect on the maximum sliding distance. Peng [12], Wang [13] and Chen [14] considered the influence of different fit gap on the stress distribution of tenon and mortise. Zeng [15] established a nonlinear dynamic model of tenon and mortise with dry friction and gap by using a lumped parameter model, and analyzed the vibration response of mistuned bladed disk system by using numerical method. It was found that mistuned bladed disk system would lead to complex nonlinear dynamic behavior. Liu [16] used the anti-aliasing time-frequency fusion algorithm and considered the three-dimensional micro sliding friction contact, established the non-linear dry friction damping mistuned bladed disk system model with 18 blades in the form of density mistuned, and analyzed the influence of the positive pressure on the contact surface and the double mistuned blade mass on the vibration characteristics of the bladed disk system. The results show that the nonlinear friction damping can effectively suppress the vibration of the mistuned bladed disk system. Santosh [17] proposed a new damping method applied to the rotating condition of wind turbine blades. In this method, containers filled with spherical particles are alternately installed at four different positions on each blade. Experiments were carried out at different speeds and different positions of the container, and different damping effects were obtained, and the best installation position of the shock absorber was found. Zhao [18] based on the bladed disk model, they used intelligent optimization algorithm to optimize the vibration reduction of bladed disk system. Mashayekhi [19] proposed an efficient model order reduction technique for forced vibration analysis of a mistuned bladed disk with friction contact. Afzal [20] investigated the damping potential of multiple friction contacts in a bladed disk. Friction contacts at tip shrouds and strip dampers are considered. It is shown that friction damping effectiveness can be potentially increased by using multiple friction contact interfaces. Yang [21] established a fast prediction method for forced vibration response of bladed disks with contact interfaces. Coulomb friction model is used to simulate the contact state between tenon and mortise.

In this paper, a lumped parameter model of the mistuned blade disk system including blade, tenon and disk is established for the compressor blade disk system of the aeroengine. The effect of the coupling stiffness, gap and friction of the tenon on the vibration characteristics of the mistuned blade disk system is considered. Blade damping mistuning and tenon damping mistuning on the vibration characteristics of the mistuned bladed disk under strong coupling and weak coupling are discussed respectively.

2. Dynamic model of the bladed disk system

The compressor bladed disk system of aero-engine is connected by blade and disk through tenon and mortise structure. The blade is installed on the disk through tenon. The tenon stiffness and damping are simplified as the lumped parameter dynamic model of the blade disk system with equivalent stiffness and damping, as shown in Fig. 1. Set the number of sectors to n, mbi, mti and mdi are the equivalent mass of blade body, blade tenon and wheel disk of the ith sector represently, kbi, kti and kdi are the equivalent stiffness of blade body, tenon and disk of the ith sector represents, cbi and cti are the equivalent damping of blade body and blade tenon of the ith sector, Fi is the exciting force acting on the ith sector, kc is the coupling stiffness between sectors. Set mbi, mti and mdi displacements are xi, xi+n, xi+2n.

When there is a small gap between tenon and mortise, the coupling effect of blade vibration in adjacent sector is variable. Set the initial gap of mortise is d0, the coupling force between disks can be simplified as the elastic force of Eq. (1):

(1)
f g = - k c x i + n - x i + 1 + n - d 0 .

If the tenon of the ith sector bears the force from the sectors on the left and right sides respectively are fL and fR:

(2)
f L = - k c ( x i + n - x i - 1 + n - d 0 ) , x i + n - x i - 1 + n > d 0 , 0 , x i + n - x i - 1 + n d 0 ,
(3)
f R = - k c ( x i + n - x i + 1 + n - d 0 ) , x i + n - x i + 1 + n > d 0 , 0 , x i + n - x i + 1 + n d 0 .

Fig. 1. Lumped parameter model of bladed disk system with tenon and mortise

 Lumped parameter model of bladed disk system with tenon and mortise

The friction between the blade tenon and the disk mortise adopts the Coulomb friction model shown in Eq. (4):

(4)
g i + n = - μ N i + n s g n x ˙ i + n ,

where, μ is the coefficient of friction; Ni+n is normal pressure; sgn is a symbolic function, when relative sliding speed x˙i+n greater than zero, equal to zero and less than zero respectively, sgnreturn values are 1, 0, and –1, respectively.

The exciting force of air flow acting on the ith blade is:

(5)
F i = F   0 s i n ω t + ϕ i , i = 1,2 , , n ,

where F 0 is the amplitude of the exciting force on the ith rotor blade, ω is the excitation frequency, n is the number of rotor blades, ϕi represents the phase angle of the traveling wave excitation force on the ith blade, as shown in Eq. (6):

(6)
ϕ i = 2 π E i - 1 n , i = 1,2 , , n ,

where is the order of excitation.

The dynamic differential equation of the ith sector of the bladed disk system is obtained by analyzing the force on the ith sector:

(7)
m b i x ¨ i + c b i x ˙ i - x ˙ i + n + k b i x i - x i + n = F i , m t i x ¨ i + n + c b i ( x ˙ i + n - x ˙ i ) + c t i ( x ˙ i + n - x ˙ i + 2 n ) + k b i ( x i + n - x i )             + k t i x i + n - x i + 2 n - f R - f L = - μ N i + n s g n x ˙ i + n , m d i x ¨ i + 2 n + k d i x i + 2 n + k t i ( x i + 2 n - x i + n ) + c t i ( x ˙ i + 2 n - x ˙ i + n )             + k c ( x i + 2 n - x i + 2 n - 1 - d o ) + k c ( x i + 2 n - x i + 2 n + 1 - d o ) = 0 .

Then:

(8)
m b i x ¨ i + c b i x ˙ i - c b i x ˙ i + n + k b i x i - k b i x i + n = F i , m t i x ¨ i + n - c b i x ˙ i + ( c b i + c t i ) x ˙ i + n - c t i x ˙ i + 2 n - k b i x i             + k b i + k t i + 2 k c x i + n - k t i x i + 2 n = - μ N i + n s g n x ˙ i + n , m d i x ¨ i + 2 n - c t i x ˙ i + n + c t i x ˙ i + 2 n - k t i x i + n + ( k t i + k d i + 2 k c ) x i + 2 n               - k c x i + 2 n - 1 - k c x i + 2 n + 1 = 2 k c d o .

The dynamic equation of the whole disk system is:

(9)
M X ¨ + C X ˙ + K X = F ,

where the displacement vector is X=x1,,xi,,xi+n,,xi+2n,,x3n3n×1T.

The force vector is F=F1,F1,,Fn,g1+n,g2+n,,g2n,2kcd0,,2kcd03n×1T.

The mass matrix is:

(10)
M = M b 0 0 0 M t 0 0 0 M d 3 n × 3 n ,

where:

M b = m b 1 m b 2 m b n n × n ,               M t = m t 1 m t 2 m t n n × n ,
M d = m d 1 m d 2 m d n n × n .

The damping matrix is:

(11)
C = C b - C b 0 - C b C b + C t - C t 0 - C t C t 3 n × 3 n ,

where:

C b = c b 1 c b 2 c b n n × n ,               C t = c t 1 c t 2 c t n n × n .

The stiffness matrix is:

(12)
K = K b - K b 0 - K b K b + K t + 2 K c - K t 0 - K t K t + K d + 2 K c + K c ' 3 n × 3 n ,

where:

K b = k b 1 k b 2 k b n n × n ,             K t = k t 1 k t 2 k t n × n ,
K d = k d 1 k d 2 k d n n × n ,         K c = k c k c k c n × n ,
K c ' = 0 - k c 0 0 - k c 0 0 - k c 0 0 - k c 0 n × n .

The dimensionless parameters of disk system are defined as follows: blade body equivalent mass mb=1, blade tenon equivalent mass mt=3.75, equivalent mass of wheel disk md=3, blade equivalent stiffness kb=1, equivalent stiffness of tenon kt=9.1, equivalent stiffness of disk kd=1.1, select coupling stiffness kc=16 is weak coupling, coupling stiffness kc=71 is strong coupling, blade equivalent damping cb=0.05, equivalent damping of tenonct=0.625, order of exciting force E=6. In order to verify the accuracy of the calculation model in this paper, Using the above parameters, through comparison with the model in reference [8], the comparison results are shown in Fig. 2. The horizontal and vertical coordinates are dimensionless physical quantities without units.

Fig. 2. Model validation

 Model validation

It can be seen from Fig. 2 that the model in this paper is in good consistent with the model in reference [15], which verifies the correctness of the model in this paper.

3. Nonlinear dynamic analysis of the Tenon and mortise of mistuned bladed disk

3.1. Nonlinear vibration characteristics of tenon with dry friction

The dynamic equation of the tuned bladed disk is:

(13)
M X ¨ + C X ˙ + K X = F .

When mass, damping and stiffness are mistuning, the dynamic equation of the mistuned bladed disk is:

(14)
M + Δ M X ¨ + C + Δ C X ˙ + K + Δ K X = F ,

where, M, C and K are mass matrix, damping matrix and stiffness matrix respectively. ΔM is the mass mistuning matrix, ΔC is the damping mistuning matrix, ΔK is the stiffness mistuning matrix:

(15)
Δ M = Δ m b 1 Δ m b 2 Δ m b n n × n ,
(16)
Δ C = Δ c b 1 Δ c b 2 Δ c b n n × n ,
(17)
Δ K b = Δ k b 1 Δ k b 2 Δ k b n n × n ,

where, Δmbi, Δcbi and Δkbi are mass mistuning, damping mistuning and stiffness mistuning of the ith blade respectively.

In this paper, the dynamic analysis of the mistuned bladed disk is carried out when only the blade stiffness mistuning is considered, and the tenon and mortise gap and friction are considered.

3.1.1. Vibration response of tenon without friction between Tenon and mortise

Considering that the coupling stiffness between the disks is strong coupling and weak coupling, the excitation frequency ω=1 is selected, and the initial gap d0=0.008 and d0=0.8 are respectively. Fig. 3 shows the vibration response of tenon without friction.

When friction is not taken into account, 2-9 SHG appears in the response spectrum of the strong coupling small gap system, and the 2, 4, 7 SHG components are larger. When the coupling stiffness is weak coupling and small gap, the maximum amplitude is 4.165, and there are multiple frequency doubling phenomena in the system response spectrum, among which the second harmonic component is the largest. There are irregular points in the section of the tenon Poincare, and the waveform in time domain is random. When the coupling stiffness is strong coupling and large gap, the maximum amplitude is 1.041, and the system response spectrum appears 2-9 frequency doubling, among which 2, 4, 7 frequency components are larger. There are irregular points in the Poincare section of the tenon, and the position of the unsmooth waveform in time domain may be caused by the influence of frequency components of 2, 4 and 7. When the coupling stiffness is weak coupling and large gap, the maximum amplitude is 4.3453. There are multiple frequency doubling phenomena in the system response spectrum, among which the second frequency component is the largest. There are irregular points in the tongued Poincare section, and the time-domain waveform is random. These characteristics indicate that the system is in a chaotic state. Under the same coupling strength, when the friction is not considered, the vibration of tenon increases with the gap.

Fig. 3. Response characteristics of tenon without friction

 Response characteristics of tenon without friction

a)kc=71, ω=1, d0=0.008

 Response characteristics of tenon without friction

b) kc=16, ω=1, d0=0.008

 Response characteristics of tenon without friction

c)kc=71, ω=1, d0=0.8

 Response characteristics of tenon without friction

d)kc=16, ω=1, d0=0.8

3.1.2. Vibration response of tenon with friction between Tenon and mortise

When the coupling stiffness is strong coupling, the excitation frequency ω=1 and the initial gap d0=0.008 are selected.

When the coupling stiffness is strong coupling and small gap, the maximum amplitude is 1.0369, and there are 2-8 frequency doubling in the system response spectrum, among which the 2, 4, 5 frequency components are larger. There are irregular points in the section of the tenon Poincare, and the position of the unsmooth time-domain waveform may be caused by the influence of the second, fourth and fifth harmonic components. When the coupling stiffness is weak coupling and small gap, the maximum amplitude is 3.7074, and there are multiple frequency doubling phenomena in the system response spectrum, among which the second, fourth and sixth harmonic components are larger. There are irregular points in the section of the tenon Poincare, and the position of the unsmooth waveform in time domain may be caused by the influence of frequency components of 2, 4 and 6. When the coupling stiffness is strong coupling and large gap, the maximum amplitude is 1.0578, there are several frequency doubling phenomena in the system response spectrum, among which the second harmonic component is the largest. There are irregular points in the section of the tenon Poincare. When the coupling stiffness is weak coupling and large gap, the maximum amplitude is 4.3453, in which 2, 3, 4 SHG components are larger. There are irregular points in the Poincare section of Tenon and mortise, and the position of unsmooth waveform in time domain may be caused by the influence of frequency components of 2, 3 and 4. These characteristics indicate that the system is in chaotic state.

From the comparison between Fig. 3 and Fig. 4, it can be seen that when the coupling stiffness is strong coupling, the vibration amplitude of tenon is smaller than that of weak coupling, and when friction is not considered, the vibration amplitude of tenon is small regardless of whether the initial gap is small or large. When the coupling stiffness is weak coupling, the vibration amplitude of tenon is larger when the friction is considered.

Fig. 4. Response characteristics of tenon with friction

Response characteristics of tenon with friction

a)kc=71, ω=1, d0=0.008

Response characteristics of tenon with friction

b) kc=16, ω=1, d0=0.008

Response characteristics of tenon with friction

c)kc=71, ω=1, d0=0.8

Response characteristics of tenon with friction

d)kc=16, ω=1, d0=0.8

3.2. Effects of tenon gap and dry friction on blade vibration response of mistuned bladed disk

It can be seen from Fig. 5 and Fig. 6 that the maximum vibration amplitude of the mistuned blade is smaller than that without considering the friction force when the gap between tenon and groove and friction force are considered. When the coupling stiffness between disks is weak, the maximum amplitude of blades in the mistuned bladed disk system considering the friction between tenons and mortise is small. Compared with Fig. 5 and Fig. 6, it can be seen that the mistuning will change the vibration amplitude of the bladed disk system, and the existence of the mistuning factor makes the dynamic behavior of the bladed disk system more complex. The coupling strength has a great influence on the vibration amplitude of the blade, and the strong coupling makes the energy transfer between sectors strengthen, resulting in the increase of the vibration amplitude of the blade.

Fig. 5. Blade vibration response amplitude of strongly coupled

 Blade vibration response amplitude of strongly coupled

a)d0=0.008, without friction

 Blade vibration response amplitude of strongly coupled

b)d0=0.008, with friction

 Blade vibration response amplitude of strongly coupled

c)d0=0.8, without friction

 Blade vibration response amplitude of strongly coupled

d)d0=0.8, with friction

Fig. 6. Blade vibration response amplitude of weak coupled

 Blade vibration response amplitude of weak coupled

a)d0=0.008, without friction

 Blade vibration response amplitude of weak coupled

b)d0=0.008, with friction

 Blade vibration response amplitude of weak coupled

c)d0=0.8, without friction

 Blade vibration response amplitude of weak coupled

d)d0=0.8, with friction

4. Effect of damping mistuning on vibration characteristics of mistuned bladed disk

4.1. Effect of blade equivalent damping mistuning on vibration characteristics of mistuned bladed disk

Considering the effect of blade equivalent damping mistuning on the vibration characteristics of bladed disk with different gap and coupling strength, the coupling stiffness of 16 and 71 is selected as the weak coupling and strong coupling, and the gap of 0.008 and 0.8 is selected as the small gap and large gap. Through analysis, the time domain waveform diagram, phase diagram, Poincare section and tenon frequency spectrum of tenon are obtained as shown in Fig. 7.

It can be seen from Fig. 7 that when the gap is the same, the vibration amplitude of the tenon is smaller than that of the weak coupling under strong coupling. In the case of strong coupling, when the gap is large, the vibration amplitude of the tenon is smaller than that when the gap is small. In the case of weak coupling, when the gap is large, the vibration amplitude of the tenon is larger than that of the small gap.

It can be seen from Fig. 8 and Fig. 9 that under the condition of the same gap, the vibration amplitude of the blade with strong coupling is larger; when the coupling stiffness and gap are the same, the equivalent damping mistuning of the blade body will increase the vibration amplitude of the blade, but the vibration amplitude of the blade will change slightly. The mistuning of damping leads to the enhancement of vibration energy transfer in each sector and the increase of blade vibration amplitude.

Fig. 7. Effect of blade equivalent damping mistuning on vibration response of tenon

 Effect of blade equivalent damping mistuning on vibration response of tenon

a)kc=16, d0=0.008

 Effect of blade equivalent damping mistuning on vibration response of tenon

b) kc=71, d0=0.008

 Effect of blade equivalent damping mistuning on vibration response of tenon

c)kc=16, d0=0.8

 Effect of blade equivalent damping mistuning on vibration response of tenon

d)kc=71, d0=0.8

Fig. 8. Response of blade equivalent damping mistuning to blade vibration with small gap (d0=0.008)

 Response of blade equivalent damping mistuning to blade vibration with small gap (d0=0.008)

a) Tuned, kc=16

 Response of blade equivalent damping mistuning to blade vibration with small gap (d0=0.008)

b) Mistuned, kc=16

 Response of blade equivalent damping mistuning to blade vibration with small gap (d0=0.008)

c) Tuned, kc=71

 Response of blade equivalent damping mistuning to blade vibration with small gap (d0=0.008)

d) Mistuned, kc=71

Fig. 9. Response of blade equivalent damping mistuning to blade vibration with large gap (d0=0.8)

 Response of blade equivalent damping mistuning to blade vibration with large gap (d0=0.8)

a) Tuned, kc=16

 Response of blade equivalent damping mistuning to blade vibration with large gap (d0=0.8)

b) Mistuned, kc=16

 Response of blade equivalent damping mistuning to blade vibration with large gap (d0=0.8)

c) Tuned, kc=71

 Response of blade equivalent damping mistuning to blade vibration with large gap (d0=0.8)

d) Mistuned, kc=71

4.2. Effect of tenon equivalent damping mistuning on vibration characteristics of mistuned bladed disk

In the case of different coupling strength and gap, the frequency spectrum of the coupling diagram with strong coupling stiffness kc=71 and weak coupling stiffness kc=16 are selected respectively. The gaps of 0.008 and 0.8 were selected as small and large gaps, the effect of blade equivalent damping mistuning on the vibration response of tenon is obtained by analysis.

It can be seen from Fig. 10 that when the gap is the same, the vibration amplitude of the tenon is smaller than that of the weak coupling under strong coupling; in the case of strong coupling, when the gap is large, the vibration amplitude of the tenon is smaller than that when the gap is small; in the case of weak coupling, when the gap is large, the vibration amplitude of the tenon is larger than that of the small gap. It can be seen from the comparison of Fig. 4 and Fig. 10 that under the small gap, blade equivalent damping mistuning increases the vibration amplitude of tenon.

Fig. 10. The effect of blade equivalent damping mistuning on the vibration response of tenon

 The effect of blade equivalent damping mistuning on the vibration response of tenon

a)kc=16, d0=0.008

 The effect of blade equivalent damping mistuning on the vibration response of tenon

b) kc=71, d0=0.008

 The effect of blade equivalent damping mistuning on the vibration response of tenon

c)kc=16, d0=0.8

 The effect of blade equivalent damping mistuning on the vibration response of tenon

d)kc=71, d0=0.8

It can be seen from Fig. 11 and Fig. 12 that when the gap is large, the equivalent damping mistuning of the blade body will increase the vibration amplitude of the blade; when the gap is small, the equivalent damping mistuning of the blade body has little influence on the vibration amplitude of the blade.

Fig. 11. The effect of tenon equivalent damping mistuning on the vibration response of blade with small gap (d0=0.008)

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with small gap (d0=0.008)

a) Tuned, kc=16

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with small gap (d0=0.008)

b) Mistuned, kc=16

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with small gap (d0=0.008)

c) Tuned, kc=71

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with small gap (d0=0.008)

d) Mistuned, kc=71

Fig. 12. The effect of tenon equivalent damping mistuning on the vibration response of blade with large gap (d0=0.8)

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with large gap (d0=0.8)

a) Tuned, kc=16

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with large gap (d0=0.8)

b) Mistuned, kc=16

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with large gap (d0=0.8)

c) Tuned, kc=71

 The effect of tenon equivalent damping mistuning  on the vibration response of blade with large gap (d0=0.8)

d) Mistuned, kc=71

4.3. Effect of mistuning of tenon dry friction damping force on vibration characteristics of mistuned bladed disk

In the case of different coupling strength and gap, the frequency spectrum of the coupling diagram with strong coupling stiffness kc= 71 and weak coupling stiffness kc=16 are selected respectively. The gaps of 0.008 and 0.8 were selected as small and large gaps, the effect of tenon dry friction damping mistuning on the vibration response of tenon is obtained by analysis.

Fig. 13. Effect of tenon’s dry friction damping force mistuning on vibration response of tenon

 Effect of tenon’s dry friction damping force mistuning on vibration response of tenon

a)kc=16, d0=0.008

 Effect of tenon’s dry friction damping force mistuning on vibration response of tenon

b) kc=71, d0=0.008

 Effect of tenon’s dry friction damping force mistuning on vibration response of tenon

c)kc=16, d0=0.8

 Effect of tenon’s dry friction damping force mistuning on vibration response of tenon

d)kc=71, d0=0.8

It can be seen from Fig. 13 that when the gap is the same, the vibration amplitude of the tenon is smaller than that of the weak coupling. When the coupling strength is strong and the gap is large, the vibration amplitude of the tenon is smaller than that the gap is small. In the case of weak coupling, when the gap is large, the vibration amplitude of the tenon is larger than that of the small gap. The gap has little effect on the vibration response of the weakly coupled bladed disk, but the gap has a great influence on the strongly coupled bladed disk, the larger the gap, the more complex the nonlinear dynamic characteristics of tenon. It can be seen from the comparison of Fig. 4 and Fig. 13 that under the same gap and coupling strength, the dry friction damping mistuning of tenon increases the vibration amplitude of tenon.

Fig. 14. Response of blade vibration when tenon’s dry friction damping force mistuning under weak coupling (kc=16)

 Response of blade vibration when tenon’s dry friction damping force mistuning  under weak coupling (kc=16)

a) Tuned d0=0.008

 Response of blade vibration when tenon’s dry friction damping force mistuning  under weak coupling (kc=16)

b) Mistuned d0=0.008

 Response of blade vibration when tenon’s dry friction damping force mistuning  under weak coupling (kc=16)

c) Tuned d0=0.8

 Response of blade vibration when tenon’s dry friction damping force mistuning  under weak coupling (kc=16)

d) Mistuned d0=0.8

Fig. 15. Response of blade vibration when tenon’s dry friction damping force mistuning under strong coupling (kc=71)

Response of blade vibration when tenon’s dry friction damping force mistuning  under strong coupling (kc=71)

a) Tuned d0=0.008

Response of blade vibration when tenon’s dry friction damping force mistuning  under strong coupling (kc=71)

b) Mistuned d0=0.008

Response of blade vibration when tenon’s dry friction damping force mistuning  under strong coupling (kc=71)

c) Tuned d0=0.8

Response of blade vibration when tenon’s dry friction damping force mistuning  under strong coupling (kc=71)

d) Mistuned d0=0.8

It can be seen from Fig. 14 and Fig. 15 that under the condition of the same gap, the vibration amplitude of the blade with strong coupling is larger; when the coupling stiffness and gap are the same, the tenon’s equivalent damping mistuning will increase the vibration amplitude of the blade, but the vibration amplitude of the blade will change slightly.

5. Conclusions

Aiming at the Tenon and mortise structure of aero-engine compressor bladed disk, a dynamic model including blade, tenon and disk was established. The nonlinear dynamic characteristics of bladed disk system with damping mistuning were studied.

1) The dynamic model and nonlinear vibration equation including blade body, tenon and disk are established. The effects of tenon gap and dry friction on the nonlinear dynamics of bladed disk are considered in the model.

2) Comparative analysis of the effects of tenon gap and dry friction on the vibration response of the mistuned bladed disk. The results show that the tenon gap and dry friction have a certain influence on the vibration response of the mistuned bladed disk, which cannot be ignored in the dynamic analysis.

3) The effects of blade equivalent damping and tenon equivalent damping mistuning on the vibration characteristics of bladed disk system are obtained: When the gap is the same, the vibration amplitude of the tenon is smaller than that of the weak coupling under the condition of strong coupling; The vibration amplitude of the tenon is smaller than that of the small gap under the condition of strong coupling; In the case of weak coupling, the vibration amplitude of tenon is larger than that of small gap when the gap is large.

4) The effects of the mistuning of the tenon’s dry friction damping force on the vibration characteristics of the mistuned bladed disk is obtained: Under the same gap, the blade vibration amplitude is larger when the coupling is strong; When the coupling stiffness and gap are the same, the equivalent damping mistuning of blade will lead to the increase of blade vibration amplitude, but the change of blade vibration amplitude is small.

Acknowledgements

This research is supported by project of the National Natural Science Foundation of China (No. 51805076, No. 51775093 and No. U1708255), China Postdoctoral Science Foundation (2019M651147), Doctoral Research Initiation Fund of Liaoning Provincial (2020-BS-156), Shenyang Youth Science and Technology Innovation Talent Project (RC200006) and Scientific Research Fund Project of Liaoning Provincial Department of Education (LG202008).

References

  1. Cha D., Sinha A. Statistics of responses of a mistuned and frictionally damped bladed disk assembly subjected to white noise and narrow band excitations. Probabilistic Engineering Mechanics, Vol. 21, Issue 4, 2000, p. 384-396. [Publisher]
  2. Sinclair G. B., Cormier N. G. Contact stresses in dovetail attachments: physical modeling. Journal of Engineering for Gas Turbines and Power, Vol. 124, Issue 2, 2000, p. 325-331. [Publisher]
  3. Sinclair G. B., Cormier N. G., Griffin J. H., et al. Contact stresses in dovetail attachments: finite element modeling. Journal of Engineering for Gas Turbines and Power, Vol. 124, Issue 1, 2002, p. 182-189. [Publisher]
  4. Sinclair G. B., Cormier N. G. Contact stresses in dovetail attachments: alleviation via precision crowning. Journal of Engineering for Gas Turbines and Power, Vol. 125, Issue 4, 2003, p. 1033-1041. [Publisher]
  5. Petrov E. P. A Method for use of cyclic symmetry properties in analysis of nonlinear multiharmonic vibrations of bladed disks (2003-GT-38480). Journal of Turbomachinery, Vol. 126, Issue 1, 2004, p. 175-183. [Publisher]
  6. Ciğeroğlu E., Nevzat Özgüven H. Nonlinear vibration analysis of bladed disks with dry friction dampers. Journal of Sound and Vibration, Vol. 295, Issue 3, 2006, p. 1028-1043. [Publisher]
  7. Wang H. J., He E. M., Yu S. X. Forced response characteristics of bladed disks with disordered dry friction. Journal of Aerospace Power, Vol. 21, Issue 4, 2006, p. 711-715. [Search CrossRef]
  8. Wang A. L., Long Q. Forced response characteristics of bladed disks with random mistuned non-linear friction damping. Journal of Aerospace Power, Vol. 26, Issue 1, 2011, p. 178-184. [Search CrossRef]
  9. Firrone C. M., Zucca S., Gola M. M. The effect of underplatform dampers on the forced response of bladed disks by a coupled static/dynamic harmonic balance method. International Journal of Non-Linear Mechanics, Vol. 46, Issue 2, 2011, p. 363-375. [Publisher]
  10. Zhang L., Yuan H. Q., Han Q. K., et al. Vibration analysis of mistuned bladed disk system based on microslip friction model. Journal of Vibration Engineering, Vol. 25, Issue 3, 2012, p. 289-293. [Search CrossRef]
  11. Ma H., W. D., Neng H. Q., et al. Contact characteristic analysis of a compressor disc dovetail structure. Journal of Northeastern University (Natural Science), Vol. 36, Issue 1, 2015, p. 109-113. [Search CrossRef]
  12. Peng M. L., Yang Z. C., Cao Y. Y., et al. Analysis on the influences of the contact gap on the contact response of turbine disc-blade structures. Journal of Aerospace Power, Vol. 26, Issue 4, 2011, p. 779-786. [Search CrossRef]
  13. Wang S. F., Wen Z. X., Gao Y. F., et al. Effects of contact gap on contact properties of turbine blade tenon /mortise. Journal of Propulsion Technology, Vol. 36, Issue 1, 2015, p. 112-118. [Search CrossRef]
  14. Chen Z. Y., Wang C., Zhou P. Research on Robust Optimization of Fatigue Life for Blade-Disk Considering Contact Gap between Tenon and Mortise. Journal of Propulsion Technology, Vol. 39, Issue 4, 2018, p. 857-864. [Search CrossRef]
  15. Zeng H. N., Cao S. Q., Su Y. L. Vibration localization of a mistuned bladed disk system with friction and gap. Journal of Vibration and Shock, Vol. 35, Issue 2, 2016, p. 82-90. [Search CrossRef]
  16. Liu Y. L., Shangguan B., Xu Z. L. Effects of dry friction damping on forced vibration response of mistuned bladed disk. Journal of Xi’an Jiaotong University, Vol. 50, Issue 2, 2016, p. 111-117. [Search CrossRef]
  17. Santosh R. Sandanshiv, Umesh S. Chavan Vibration suppression effects on rotating wind turbine blade using a particle damping method. Vibroengineering Procedia, Vol. 29, 2019, p. 43-48. [Publisher]
  18. Zhao T. Y., Yuan H. Q., Yang W. J., Sun H. G. Genetic particle swarm parallel algorithm analysis of optimization arrangement on mistuned blades. Engineering Optimization, Vol. 49, Issue 4, 2017, p. 1-22. [Publisher]
  19. Mashayekhi F., Nobari A. S., Zucca S. Hybrid reduction of mistuned bladed disks for nonlinear forced response analysis with dry friction. International Journal of Non-Linear Mechanics, Vol. 116, 2019, p. 73-84. [Publisher]
  20. Afzal M., Arteaga I. L., Kari L. Numerical analysis of multiple friction contacts in bladed disks. International Journal of Mechanical Ences, Vol. 137, 2018, p. 224-237. [Publisher]
  21. Yang G., Zhou B., Zang C. P. Fast response prediction method for bladed disks with contact interfaces. Journal of Aerospace Power, Vol. 34, Issue 9, 2019, p. 1953-1961. [Search CrossRef]