Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation

Zhang, Hongyuan; Yuan, Huiqun; Sun, Hongyun

doi:10.21595/jve.2019.20740

Journal of Vibroengineering

Browse Journal

Submit article

Published: 30 September 2019

Check for updates

Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation

Hongyuan Zhang¹

Huiqun Yuan²

Hongyun Sun³

¹School of Automobile and Traffic, Shenyang Ligong University, Shenyang, China

²College of Sciences, Northeastern University, Shenyang, China

^{2, 3, 1}School of Mechanical Engineering and Automation, Northeastern University, Shenyang, China

Corresponding Author:

Hongyuan Zhang

Cite the article Download PDF

Downloads 1519

Abstract

This paper presents a method of mistuning identification of a random mistuned bladed parameter, to predict the forced vibration characteristics of random mistuned blade disk of aeroengine compressor, based on the prestress mode synthesis method, substructure finite element reduction model was adopted, the identification method combines Weibull distribution with Monte Carlo simulation to identify the random mistuned blade parameters. The vibration characteristics of mistuned bladed disks under different engine orders of excitation were calculated, and the vibration characteristics of mistuned bladed disk systems under different mistuning standard deviations were discussed. Vibration localization characteristics of mistuned bladed disk system were analysed by using localization factor. The correctness of the identification method was verified by calculating the vibration amplitude of mistuned bladed disk system under different mistuning standard deviations.

Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation

Highlights

A method of mistuning identification of a random mistuned bladed parameter based on Weibull distribution and Monte Carlo simulation is proposed.
The vibration characteristics of mistuned bladed disks under different engine orders of excitation are studied.
The vibration characteristics of mistuned bladed disks under different mistuning standard deviations are discussed.
Vibration localization characteristics of mistuned bladed disk system are analysed by localization factor.

1. Introduction

In theory, the blade disk system should be a circular symmetrical structure, but due to processing, installation, wear and in order to suppress flutter, the mass, stiffness and natural frequency of each blade are different, that is, mistuning. Mistuning will cause the vibration energy of the blade disk system to concentrate on a few blades, resulting in serious vibration localized, leading to fatigue failure, thereby affecting the service life of the compressor blade disk system. Therefore, it is very important to study the vibration response of the blade disk system caused by blade mistuning. For mistuned blade disk system, the correct identification of mistuned parameters is the key to the dynamic analysis of the mistuned blade disk system. The blade mistuning can be divided into mass mistuning, damping mistuning and stiffness mistuning according to the mistuning parameters, and random mistuning and deterministic mistuning according to the mistuning form. Random mistuning is based on a specific distribution function, using different mean and standard deviation of mistuning to deal with the mistuning of blade parameters. Correct identification of mistuning parameters is the key to the simulation of the blade disk system. Since SNM and FMM do not decompose the substructures, the key to the identification of mistuning parameters is to select a reasonable modal set. The corresponding identification methods are SNM identification method [1] and FMM identification method [2].

Mignolet [3] in order to accurately predict the statistical characteristics of the forced response of mistuned bladed disks, three different methods are proposed to solve the problem of insufficient data, and the three methods are validated on the lumped parameter model. Pichot [4] presented a method for identifying the mistuning of blade disk based on the system mode and the measurement of natural frequency. This method establishes the finite element model of blade disk, and presents a computational reduced-order modelling technology based on component-mode substitution method. Laxalde [5] proposed a model updating method based on the reduced-order model established by modal synthesis method. The method takes the measured modal data as input, and the modal data are extracted from the forced response measurement of the whole blade disk. In the process of mistuning identification, the measured data are regularized to reduce the impact of the final measurement error. Madden [6, 7] proposed an inverse reduced-order model suitable for solving the inverse problem of mistuning identification and a technique suitable for identifying high-order modal families. This improves the overall accuracy and robustness of mistuning identification. Wang [8] aiming at the problem of mistuning identification of integral blade disk, a method based on the classical modal method was proposed, which combines the analytical modal of the tuned blade disk system with the measurement modal of mistuned structure. Judge [9] proposed two methods for identifying single blade mistuning by measuring the response of the blisks, which is especially suitable for the blisks whose single blade cannot be tested separately. Lu [10] analysed the frequency response of a piezoelectric smart beam rotating in a plane by the Hamilton's principle.

Nyssen [11, 12] studied the experimental modal analysis of the theoretical blade disk using basic excitation to identify the mistuning of each blade. In this method, the optical measurement is used to obtain the precise geometric structure of the structure, and the geometric mistuning can be correlated to each blade. Weber [13] proposed a method to determine modal characteristics by experiments and optimizes random mistuning pattern to match these characteristics. Jamia [14] proposed a blade mistuning index (WEBMI), which was obtained by the wavelet packet transform of blade tuning and mistuning response. The proposed mistuning index was proved on the lumped parameter model. The experimental results show that the proposed method can identify fine mistuning.

In this paper, a mistuned parameter identification method based on Weibull distribution and Monte Carlo simulation is proposed. By calculating the parameters identified by this method, only a small number of samples can be used to obtain results consistent with a large number of samples, which greatly reduces the calculation time.

Based on the reduced order model of the Prestressed Component Mode Synthesis Method in Ref. [15], this study mainly focuses on forced response characteristics of the random mistuned blade which are not involved in Ref. [15]. This paper based on the Prestressed Component Mode Synthesis Method，the finite element reduction model of mistuned blade disk system was adopted, the mistuning parameters were identified by the combination of Weibull distribution and Monte Carlo simulation technology, the vibration response of the mistuned bladed disk system were studied. The vibration response of mistuned bladed disk system under different mistuned standard deviations was discussed. The localization factor was defined and the vibration localization characteristics of mistuned bladed disk system were analysed by the localization factor. The identification method was verified by calculating the vibration amplitude under the different mistuning standard deviation.

2. Identification method of random mistuned blade parameters

For actual engineering blade disk system, mistuning destroys the cyclic symmetry and results in different stiffness and frequency of each blade. If the lumped mass model was used to study the localization characteristics of random mistuning blade disk vibration, which cannot accurately reflect the vibration characteristics of actual engineering blade disk system.

In order to overcome the problem of a large number of nodes and elements in integral blade disk, a finite element reduction model is established by using the prestressing mode synthesis method described in reference [15], refer to reference [15] for the prestressing modal synthesis method, which is not described in detail here. Fig. 1(a) is a model of integral blade disk after cutting a sector to show the cross section of the substructure sector, and Fig. 1(b) is a substructure model. The basic sector of the substructure is divided into 52163 elements, and the number of nodes is 79174. The material properties of the model are listed in Table 1.

If the integral blade disk model was used, but because of the large number of elements and nodes of the integral blade disk finite element model, it is difficult to calculate. Even if the blade disk system substructure finite element reduction model shown in Fig. 1 is used to reduce the degree of freedom of the model, the sample size will be too large, which will lead to a long calculation time. It is also difficult to obtain the statistical characteristics of the random vibration of the mistuned blade disk system.

Fig. 1Finite element model of bladed disk with one sector cut out

a) Bladed disk

b) Substructure model

Table 1Blisk material properties

Balde material property	Property value	Disk material property	Property value
Mass density	4400 kg⋅m⁻³	Mass density	4700 kg⋅m⁻³
Modulus of Elasticity	113 GPa	Modulus of Elasticity	150 GPa
Poisson’s ratio	0.3	Poisson’s ratio	0.3

Since the maximum value of the structural vibration response of a random mistuned bladed disk system is studied, the maximum value of the structural vibration response of the bladed disk obeys the Weibull distribution from the theory of extreme value statistics. In the literature [16, 17], the Weibull distribution is used to analyze the random vibration characteristics of mistuned bladed disk system. The results show that a large number of random mistuned data can be simulated by using a small number of random mistuned data satisfying Weibull distribution, and the results are consistent with a large number of data. In this paper, the random mistuning problem is studied by combining the Weibull distribution with Monte Carlo simulation.

The distribution function of the three-parameter Weibull distribution is:

1

F (x) = 1 - \exp [- {(\frac{t - γ}{η})}^{β}],

where, $γ$ , $η$ and $β$ are the position parameter, the scale parameter and the shape parameter are respectively, $γ \geq 0$ , $η > 0$ , $β > 0$ , $t$ is random variables, and here is the maximum dimensionless amplitude of the mistuned blade.

The probability density function of the Weibull distribution is:

2

f (x) = \frac{β}{η} {(\frac{t - γ}{η})}^{β - 1} \exp [- {(\frac{x - γ}{η})}^{β}] .

In order to make the sample obey the Weibull distribution, it is necessary to determine the estimated values of three-parameter $γ$ , $η$ and $β$ . Making the three-parameter Weibull distribution as close as possible to the distribution determined by sample data.

The probability density function of Weibull distribution of Eq. (2) is logarithmic on both sides:

3

l n [- l n (1 - F (t)] = β l n (t - γ) - β l n (η) .

The $t$ value of the formula is obtained from Monte Carlo simulation data, where the $t$ value is the dimensionless of vibration amplitude, and the position parameter $γ$ is the maximum value of $t$ . As can be seen from the reference [18], the position parameter $γ$ value can be approximated by the following equation for mistuned blade disk vibration:

4

γ \approx (\sqrt{N + 1}) / 2 .

For Eq. (3), make Weibull transformation, $x$ and $y$ are given by:

5

\{\begin{array}{l} x = l n t, \\ y = l n [- l n (1 - F (t)] . \end{array}

When $t - γ = 0$ , $y \to - \infty$ , when $t \to + \infty$ , $y = β x - β l n η$ .

The $t$ data obtained by Monte Carlo simulation are arranged from small to large and can be obtained by using the median rank formula [19]:

6

\begin{array}{l} F (t) = \frac{i - 0.3}{n + 0.4}, & (i = 1,2, \dots, n) \end{array} .

For $t_{1}, y_{1}$ ; $t_{2}, y_{2}$ ;…; $t_{n}, y_{n}$ calculated $n$ times, if there is a linear relationship between $y$ and $t$ , the data calculated above satisfies the requirement:

7

\begin{array}{l} y_{i} = B_{0} + B_{1} x_{i} + ε_{i}, & (i = 1,2, \dots, n), \end{array}

where $B_{0}$ and $B_{1}$ are regression coefficients, $ε_{i}$ independent of each other and subject to $N (0, σ^{2})$ distribution, the expression becomes:

y = [\begin{matrix} y_{1} \\ y_{2} \\ ⋮ \\ y_{n} \end{matrix}], β = [\begin{array}{l} β_{1} \\ β_{2} \end{array}], ε = [\begin{matrix} ε_{1} \\ ε_{2} \\ ⋮ \\ ε_{n} \end{matrix}], Z = [\begin{matrix} 1 & x_{1} \\ 1 & x_{2} \\ ⋮ & ⋮ \\ 1 & x_{n} \end{matrix}] .

Eq. (7) can be written as:

8

y = Z β + ε,

i.e.

9

\hat{y} = B_{0} + B_{1} x .

The Eq. (9) is called empirical regression equation.

In order to obtain estimates ${\hat{B}}_{0}$ and ${\hat{B}}_{1}$ of $B_{0}$ and $B_{1}$ , $Q (B_{0}, B_{1})$ is given by:

Q (B_{0}, B_{1}) = \sum_{i = 1}^{n} (y_{i} - B_{0} - B_{1} x_{i})^{2},

represents all computed data and:

10

y = B_{0} + B_{1} x,

sum of squares of deviations. The least square method is to determine ${\hat{B}}_{0}$ and ${\hat{B}}_{1}$ so that the square of the deviation and $Q$ are minimized, $Q ({\hat{B}}_{0}, {\hat{B}}_{1})$ is given by $Q ({\hat{B}}_{0}, {\hat{B}}_{1}) = m i n Q (B_{0}, B_{1})$ .

According to the extremum principle of multivariate functions, $B_{0}$ and $B_{1}$ must satisfy the equations:

\{\begin{array}{l} \frac{\partial Q}{\partial B_{0}} = - 2 \sum_{i = 1}^{n} (y_{i} - B_{0} - B_{1} x_{i}) = 0, \\ \frac{\partial Q}{\partial B_{1}} = - 2 \sum_{i = 1}^{n} (y_{i} - B_{0} - B_{1} x_{i}) x_{i} = 0 . \end{array}

Upper equation can be rearranged as follows:

11

\{\begin{array}{l} n B_{0} + \sum_{i = 1}^{n} B_{1} x_{i} = \sum_{i = 1}^{n} y_{i}, \\ \sum_{i = 1}^{n} B_{0} x_{i} + \sum_{i = 1}^{n} B_{1} x_{i}^{2} = \sum_{i = 1}^{n} y_{i} x_{i} . \end{array}

Divide the two sides of the first equation by $n$ and rewrite it as follows:

12

{\hat{B}}_{0} = \bar{y} - {\hat{B}}_{1} \bar{x},

where:

\bar{y} = \frac{1}{n} \sum_{i = 1}^{n} y_{i}, \bar{x} = \frac{1}{n} \sum_{i = 1}^{n} x_{i} .

Substituting Eq. (12) into the Second Equation in Eq. (11), and rewrite it as follows:

13

{\hat{B}}_{1} \sum_{i = 1}^{n} (x_{i} - \bar{x}) x_{i} = \sum_{i = 1}^{n} (y_{i} - \bar{y}) x_{i},

i.e.

14

{\hat{B}}_{1} = \frac{\sum_{i = 1}^{n} (y_{i} - \bar{y}) x_{i}}{\sum_{i = 1}^{n} (x_{i} - \bar{x}) x_{i}} .

The data and $γ$ obtained by Monte Carlo simulation are substituted into Eq. (12) and Eq. (14), and the estimates ${\hat{B}}_{0}$ and ${\hat{B}}_{1}$ of $B_{0}$ and $B_{1}$ are obtained. Finally, the scale parameter $η$ and shape parameter $β$ are obtained by Eq. (3).

Whether the calculated data have good linearity should be tested, that is, the fitting degree test of a regression equation. Generally, the size of the correlation coefficient $R$ is used to describe the degree of correlation between two variables, and its mathematical expression is as follows:

15

R = \frac{\sum_{i = 1}^{n} (x_{i} - \bar{x}) (y_{i} - \bar{y})}{\sqrt{\sum_{i = 1}^{n} (x_{i} - \bar{x})^{2} \sum_{i = 1}^{n} (y_{i} - \bar{y})^{2}}} .

The closer the absolute value of $R$ is to 1, the better the regression line fits the calculated data points. When $R =$ 1, the two variables were positively correlated, when $R =$ –1, the two variables were negatively correlated; when $R \approx 0$ at that time, it was meaningless to calculate the scattering of data points along both sides of the regression line, that is, the regression line. Therefore, after finding the linear equation by regression method, the value of correlation coefficient $R$ must be calculated, then the significant value of $R$ can be found in the test table of correlation coefficient $R$ according to the number of calculation groups, and then the degree of fitting can be judged.

For the statistical characteristics of vibration response of random mistuned bladed disk system, this paper adopts the method of combining Weibull distribution with Monte Carlo simulation. The Monte Carlo method is a very important numerical calculation method guided by probability and statistics theory. It refers to the use of random numbers to solve many computational problems.

Fig. 2Identification flow chart of random mistuned blade parameters

Fig. 2 is an identification flow chart of random mistuned blade parameters, the analysis steps of this method are as follows:

(1) According to Weibull distribution, a group of mistuning values of elastic modulus of blades is generated randomly.

(2) The mistuning of the elastic modulus of the blade is written into the command stream file of the finite

element analysis, and the finite element calculation are carried out.

(3) The amplitude of blade vibration response is extracted through the analysis results of finite element method.

(4) Dimensionalization of blade vibration response amplitude, and the three-parameter estimation and correlation test of Weibull distribution are carried out.

(5) Repeat steps (1)-(4) to obtain the amplitude of the vibration response of the random mistuned bladed disk system.

(6) Statistical analysis of the obtained amplitude is carried out to obtain the statistical law that satisfies the Weibull distribution.

The mistuning parameters of the blade are analysed by introducing random mistuning of blade stiffness when the mistuning standard deviation is 1 %, 3 %, 5 % and 8 %. The maximum amplitude of blade vibration under each mistuning standard deviation is calculated by combining the Weibull distribution with Monte Carlo simulation technology.

3. Forced vibration response of mistuned blade disk system

When the frequency of the excitation force is equal to the dynamic frequency of the compressor blade disk system and the order of the excitation force is equal to the nodal diameter, the resonance will occur. Because the actual force situation of the blade disk system is more complex, the aerodynamic excitation is usually simplified to a single point excitation acting on the blade tip node in the dynamic analysis.

The dynamic equation of blade disk vibration is as follows:

16

[M] \{\ddot{X}\} + [C] \{\dot{X}\} + [K] \{X\} = \{F\},

where, [ $M$ ] is the mass matrix, [ $C$ ] is the damping matrix, [ $K$ ] is the stiffness matrix, [ $F$ ] is the excitation force vector, and the excitation force component acting on the $i$ th blade is:

17

\begin{array}{l} F_{i} = F_{i}^{0} s i n (ω t + ϕ_{i}), & i = 1,2, \dots, N \end{array},

F_{i} = F_{i}^{0} e^{j (ω t + ϕ_{i})} = \{F_{i}^{0} e^{j ϕ_{i}}\} e^{j ω t} = \{F_{i}^{0} c o s ϕ_{i} + j F_{i}^{0} s i n ϕ_{i}\} e^{j ω t},

where, $F_{i}^{0}$ , $ω$ , $N$ , $ϕ_{i}$ are the magnitude, frequency, number of blades and phase angle of the excitation force on the $i$ th blade are respectively, $ϕ_{i}$ is given as:

18

\begin{array}{l} ϕ_{i} = \frac{2 π E (i - 1)}{N}, & i = 1,2, \dots, N, \end{array}

where, $E$ is engine orders of excitation.

The real and imaginary values applied in finite element analysis can be obtained by Eq. (19) and Eq. (20).

Real part value:

19

F_{r} = F_{i}^{0} c o s [\frac{2 π E (i - 1)}{N}] .

Imaginary part value:

20

F_{i} = F_{i}^{0} s i n [\frac{2 π E (i - 1)}{N}],

where, $F_{i}^{0} =$ 8; $i$ is the number of blades, $i =$ 1, 2, …, $N$ . In this study, the damping ratio is selected according to the actual working conditions of the blade disk, $ε$ . According to Eq. (18), the magnitude of the phase angle between adjacent blades is determined by the order of excitation.

The Stator blade number and rotor blade number are $N_{s}$ and $N$ respectively. Then the excitation phase difference $φ_{e}$ is:

21

φ_{e} = \frac{2 π}{N} N_{s},

and the mode phase difference $φ_{v}$ can be expressed as:

22

φ_{v} = \frac{2 π}{N} E .

When the order of excitation is equal to the number of nodal diameter, the blade disk system resonates, i.e.:

23

φ_{e} = 2 k π \pm φ_{v},

where, $k =$ 0, 1, 2, 3….

The above equation can be rewritten as:

24

φ_{e} = φ_{v} + \frac{2 π}{N} (N_{s} - E) .

When the number of stator blades is larger than that of rotor blades, the following relationship is satisfied for the order of excitation:

25

E = N_{s} - N,

so:

26

φ_{e} = φ_{v} + 2 π .

The condition described in Eq. (23) is satisfied and resonance can be induced. When the standard deviation of elastic modulus mistuning is 1 %, 3 %, 5 % and 8 %, respectively, the random mistuning of the blade is introduced to analyse the influence of integer order excitation on the blade disk system. The maximum excitation amplitude of blade mistuning system is 8, and the frequency range of excitation frequency tuning system is 0.96-1.16 near the maximum dimensionless resonance frequency of 1.062. Table 2 and Fig. 3 are a set of random sampling blade mistuning values. The calculation results using the random mistuning values shown in Fig. 4.

From Fig. 4, it can be seen that the frequency response curves of each blade in the blade disk system under different excitation orders were obtained when the blade is mistuned. Under the condition of blade mistuning, the vibration response of each blade in the blade disk system shows great difference. The maximum resonance peak value of mistuning system is much higher than the minimum one. All these indicate that the vibration response of each blade is very different in mistuned blade disk system, and there is obvious localization of vibration response.

When the standard deviation of mistuning is 3 %, the frequency response curves of each blade in the blade disk system under different excitation orders are obtained. It can be clearly seen from the figure that under the condition of blade mistuning, the vibration responses of each blade in the blade disk system show great differences. The maximum resonance peak of the mistuning system appears at the engine orders of excitation of 1 and the maximum value is 1.9067.

Table 2Mistuning value of each blade (Pa)

Blade number	The standard deviation is 1 %	The standard deviation is 3 %	The standard deviation is 5 %	The standard deviation is 8 %
1	–2.96×10⁸	3.27×10⁹	–1.11×10⁹	–3.16×10⁸
2	5.03×10⁸	4.22×10⁸	–1.23×10⁹	–7.25×10⁹
3	4.45×10⁸	4.89×10⁹	–1.72×10⁹	9.25×10⁹
4	–1.42×10⁹	–6.68×10⁹	1.31×10⁸	–1.21×10⁹
5	–1.08×10⁹	–6.73×10⁸	2.91×10⁸	–6.49×10⁹
6	–8.41×10⁸	–4.11×10⁹	4.69×10⁹	1.23×1010
7	–5.76×10⁸	9.90×10⁹	8.67×10⁹	–2.04×10⁹
8	–3.64×10⁸	2.81×10⁹	2.65×10⁹	–5.35×10⁹
9	1.42×10⁷	4.70×10⁹	–1.19×10⁹	–2.67×10⁹
10	–3.44×10⁹	–3.60×10⁹	3.55×10⁹	–7.70×10⁹
11	–5.19×10⁸	–1.60×10⁹	1.04×10⁹	–1.02×10¹⁰
12	1.41×10⁹	–9.28×10⁸	–5.84×10⁹	2.29×10¹⁰
13	–1.21×10⁹	3.74×10⁹	5.39×10⁹	1.50×10¹⁰
14	1.06×10⁹	–9.46×10⁸	1.74×10⁹	2.79×10⁹
15	3.98×10⁸	2.39×10⁹	7.67×10⁸	–1.14×10¹⁰
16	–3.29×10⁷	–6.99×10⁹	2.92×10⁹	–7.86×10⁹
17	2.07×10⁸	–1.20×10⁹	1.48×10⁹	–1.60×10⁹
18	–1.78×10⁹	–2.80×10⁹	–5.34×10⁹	7.19×10⁹
19	–9.60×10⁷	–5.37×10⁹	–9.21×10⁸	–1.21×10¹⁰
20	1.82×10⁹	1.73×10⁹	–8.29×10⁸	–2.12×10¹⁰
21	1.12×10⁸	9.60×10⁸	–3.02×10⁹	–1.32×10¹⁰
22	4.70×10⁷	1.14×10⁸	9.55×10⁹	3.03×10⁹
23	–8.33×10⁸	–4.54×10⁹	–4.97×10⁹	3.55×10⁹
24	–3.50×10⁷	3.84×10⁹	–2.75×10⁹	4.10×10⁹
25	2.64×10⁸	1.19×10⁹	–4.04×10⁹	–1.18×10⁹
26	4.84×10⁸	–1.02×10⁹	–6.66×10⁹	1.67×10⁹
27	–4.23×10⁸	7.79×10⁷	–1.09×10⁹	–4.32×10⁹
28	–2.68×10⁸	–8.92×10⁸	–1.56×10⁹	7.83×10⁹
29	2.30×10⁹	–5.96×10⁹	8.68×10⁹	–1.24×10¹⁰
30	–2.56×10⁹	–9.73×10⁸	–1.41×10⁹	4.13×10⁹
31	2.53×10⁹	–2.83×10⁹	–6.04×10⁹	–7.71×10⁹
32	3.83×10⁸	–3.33×10⁹	9.10×10⁹	–3.04×10⁹
33	1.14×10⁹	–3.94×10⁹	7.01×10⁹	5.02×10⁹
34	–1.89×10⁹	–1.82×10⁹	–1.30×10⁹	9.43×10⁹
35	–6.70×10⁸	–6.82×10⁹	–8.55×10⁹	–1.01×10¹⁰
36	–3.16×10⁸	3.28×10⁹	–2.52×10⁹	1.14×10¹⁰
37	4.80×10⁸	1.77×10⁹	–8.85×10⁸	5.99×10⁹
38	–1.90×10⁹	–6.82×10⁷	1.57×10⁹	–6.16×10⁸

From Fig. 4, it can be seen that the frequency response curves of each blade in the blade disk system under different excitation orders were obtained when the mistuning standard deviation is increased to 5 %. It can be seen from the figure that with the increase of the standard deviation of mistuning, the resonance band of the blade-disk system widens to 0.05, and there are many peaks. Under the condition of mistuning, the vibration response of each blade in the blade-disk system shows great difference.

When the mistuning standard deviation increases to 8 %, the frequency response curves of each blade in the blade disk system under different excitation orders are obtained. It can be clearly seen from the figure that with the increase of mistuning standard deviation, the resonance band of the blade disk system widens further. The maximum bandwidth reaches 0.1. The bandwidth is twice that of the standard deviation of 5 %. There are many peaks. Under the condition of blade mistuning, the vibration response of each blade in the blade disk system shows great difference.

Fig. 3Mistuning value of each blade

a) The standard deviation is 1 %

b) The standard deviation is 3 %

c) The standard deviation is 5 %

d) The standard deviation is 8 %

Fig. 4Amplitude-frequency characteristics of mistuned bladed disk system

a) The standard deviation is 1 %

b) The standard deviation is 3 %

c) The standard deviation is 5 %

d) The standard deviation is 8 %

In order to evaluate the effect of different mistuning standard deviations on the localization of the structural vibration of the mistuned bladed disk system under different excitation orders, the maximum displacement and strain energy of the structure of the tuned and mistuned bladed disk system under different excitation orders are extracted respectively, and calculated by the following formula of vibration localization factors.

Displacement localization factor:

27

L_{u} = \frac{U_{m, m a x E i} - U_{t, m a x E i}}{U_{t, m a x E i}} \times 100 %,

where, $U_{m, m a x E i}$ is the maximum amplitude of mistuned bladed disk system for engine orders of excitation is $i$ , $U_{t, m a x E i}$ is the maximum amplitude of the tuned bladed disk system for engine orders of excitation is $i$ .

Strain energy localization factor:

28

L_{S} = \frac{S_{m, m a x E i} - S_{t, m a x E i}}{S_{t, m a x E i}} \times 100 %,

where, $S_{m, m a x E i}$ is the maximum strain energy of the mistuned bladed disk system for engine orders of excitation is $i$ , $S_{t, m a x E i}$ is the maximum strain energy of the tuned bladed disk system for engine orders of excitation is $i$ .

Fig. 5Vibration localization characteristics of mistuned bladed disk system

a) Displacement localization factor

b) Strain energy localization factor

Fig. 6Amplitude for different mistuning standard deviation

The Weibull distribution and mote Carlo simulation methods proposed in this paper are used for data sampling calculation, and the statistical results as shown in Fig. 5 and Fig. 6 are obtained by analysing the sampling data. As can be seen from Fig. 5, with the increase of excitation order, the displacement localization factor increases first, and then decreases, reaching its maximum at engine orders of excitation 3, with a higher degree of localization, and then decreases with the increase of excitation order. With the increase of excitation order, the strain energy localization factor increases first and then decreases when the mistuning standard deviation is 1 % and 5 %, reaching the maximum at engine orders of excitation 3, and then decreases with the increase of excitation order.

Fig. 6 shows the variation trend of the maximum dimensionless amplitude of mistuned bladed disk system with the standard deviation of mistuning, which is similar to the conclusions obtained from the analysis in reference [20], thus verifying the accuracy of the identification method.

4. Conclusions

A method for mistuning identification of mistuned bladed disk has been presented, which is combined Weibull distribution with Monte Carlo Simulation. Based on the Prestressed Component Mode Synthesis Method, the finite element reduce model of bladed disk system is used to identify the mistuning parameters for random mistuning blades. The forced vibration response characteristics of mistuned bladed disk system under different engine orders of excitation are calculated, and the vibration response of the mistuned bladed disk system under different mistuning standard deviations are discussed. The localization factor is defined and the vibration localization characteristics of mistuned bladed disk system are analysed by the localization factor. The following conclusions are obtained:

1) The mistuning identification method of mistuned bladed disk has been presented, which is combined Weibull distribution with Monte Carlo Simulation. This method can use smaller samples to obtain results consistent with larger samples. Improving computational efficiency greatly.

2) The frequency response curves of each blade in the blade disk system under different excitation orders are obtained when the blade is mistuned. When the condition of blade mistuning, the vibration responses of each blade in the blade disk system show great differences.

3) With the increase of excitation order, the displacement localization factor increases first, and then decreases, reaching its maximum at engine orders of excitation 3, with a higher degree of localization, and then decreases with the increase of excitation order. With the increase of excitation order, the strain energy localization factor increases first and then decreases.

4) The method is validated by comparing simulation results with relevant literature.

References

Yang M. T., Griffin J. H. A Normalized modal eigen value approach for resolving modal interaction. Journal of Engineering for Gas Turbines and Power, Vol. 119, Issue 3, 1996, p. 491-499.

Publisher
Feiner D. M., Griffin J. H. Mistuning Identification of bladed disks using a fundamental mistuning model: part 1-theory. Journal of Turbomachinery, Vol. 126, Issue 1, 2004, p. 150-158.

Publisher
Mignolet M. P., Rivas-Guerra A. J., Delor J. P. Identification of mistuning characteristics of bladed disks from free response data – part I. Journal of Engineering for Gas Turbines and Power, Vol. 123, Issue 2, 2001, p. 395.

Publisher
Pichot F., Laxalde D., Sinou J. J., Thouverez F., Lombard J. P. Mistuning identification for industrial blisks based on the best achievable eigenvector. Computers and Structures, Vol. 84, Issues 29-30, 2006, p. 2033-2049.

Publisher
Laxalde D., Thouverez F., Sinou J. J., Baumhauer S., Lombard J. P., et al. Mistuning identification and model updating of an industrial blisk. ASME Turbo Expo: Power for Land, Sea, and Air, 2006.

Publisher
Madden A. C., Castanier M. P., Epureanu B. I. Reduced-order model construction procedure for robust mistuning identification of blisks. AIAA Journal, Vol. 46, Issue 11, 2008, p. 2890-2898.

Publisher
Madden A. C., Castanier M. P., Epureanu B. I. Mistuning identification of blisks at higher frequencies. AIAA Journal, Vol. 49, Issue 6, 2011, p. 1299-1302.

Publisher
Wang S., Wang J. J., Li Q. H. Mistuning identification of integrally bladed disk based on the modal reduced technique. Journal of Aerospace Power, Vol. 24, Issue 3, 2009, p. 662-669.

Search CrossRef
Judge J. A., Pierre C., Ceccio S. L. Experimental mistuning identification in bladed disks using a component-mode-based reduced-order model. AIAA Journal, Vol. 47, Issue 5, 2009, p. 1277-1287.

Publisher
Lu En, Li Wei, Yang Xuefeng, Wang Yuqiao, Liu Yufei Dynamic modeling and analysis of a rotating piezoelectric smart beam. International Journal of Structural Stability and Dynamics, Vol. 18, Issue 1, 2018, https://doi.org/10.1142/S0219455418500037.

Publisher
Nyssen F., Golinval J. C. Identification of mistuning and model updating of an academic blisk based on geometry and vibration measurements. Mechanical Systems and Signal Processing, Vol. 68, Issue 69, 2016, p. 252-264.

Publisher
Nyssen F., Epureanu B., Golinval J. C. Experimental modal identification of mistuning in an academic two-stage drum. Mechanical Systems and Signal Processing, Vol. 88, 2017, p. 428-444.

Publisher
Weber R., Kühhorn A. Mistuning identification approach with focus on high-speed centrifugal compressors. Journal of Engineering for Gas Turbines and Power, Vol. 141, Issue 3, 2018, p. 1-10.

Publisher
Jamia N., Rajendran P., El-Borgi S., Friswell M. I. Mistuning identification in a bladed disk using wavelet packet transform. Acta Mechanica, Vol. 229, Issue 3, 2018, p. 1275-1295.

Publisher
Zhang H. Y., Yuan H. Q., Yang W. J., Zhao T. Y. Study on localization influences of frequency veering on vibration of mistuned bladed disk. Journal of Mechanical Science and Technology, Vol. 31, Issue 11, 2017, p. 5173-5184.

Publisher
Castanier M. P., Pierre C. Modeling and analysis of mistuned bladed disk vibration: status and emerging directions. Journal of Propulsion and Power, Vol. 22, Issue 2, 2006, p. 354-396.

Publisher
Lan H. Q., Zang C. P. Effect of blade stagger angle on forced response of bladed disk assemblies. Journal of Aerospace Power, Vol. 27, Issue 11, 2012, p. 2547-2552.

Search CrossRef
Whitehead D. S. Effect of mistuning on the vibration of turbo-machine blades induced by wakes. Journal of Mechanical Engineering Science, Vol. 8, Issue 1, 1996, p. 15-21.

Publisher
Nelson W. Applied Life Data Analysis. Wiley-Interscience, 2004.

Search CrossRef
Yao J. Y., Wang J. J., Li Q. H. Dynamic characteristics assessment of mistuned bladed disk using nodal diameter spectrum of mode shapes. Journal of Propulsion Technology, Vol. 32, Issue 5, 2011, p. 645-653.

Search CrossRef

About this article

Received

21 April 2019

Accepted

02 September 2019

Published

30 September 2019

SUBJECTS

Mechanical vibrations and applications

DOI

https://doi.org/10.21595/jve.2019.20740

Keywords

mistuned bladed disk

identification method

prestressed mode synthesis method

Weibull distribution

Monte Carlo simulation

engine orders of excitation

Acknowledgements

This research is supported by project of the National Natural Science Foundation of China (No. 51775093), China Postdoctoral Science Foundation (2019M651147), Natural Science Foundation of Liaoning Province (cn) (No. 20170540782) and Basic scientific research projects in Colleges and Universities of Liaoning Province (cn) (No. LG201713).

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.

Previous article in issue Previous Next article in issue Next

Research article

2022 07 07

Vibration characteristics of mistuned multistage bladed disks of the aero-engine compressor

Hongyuan Zhang, Haiyu Feng, Xiaomei You, Tianyu Zhao

Research article

2021 03 02

Effect of mistuning parameters on dynamic characteristics of mistuned bladed disk

Hongyuan Zhang, Haocun Dong, Tianyu Zhao, Xin Zhang, Huiqun Yuan

Research article

2020 12 04

Vibration reduction optimization for mistuned bladed disk based on reduced order modeling technique

Hongyuan Zhang, Huiqun Yuan, Hongyun Sun

Research article

2020 08 15

Effect of Coriolis force on vibration localization characteristics of mistuned bladed disk

Hongyuan Zhang, Huiqun Yuan, Hongyun Sun

H. Zhang, H. Yuan, and H. Sun, “Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation,” Journal of Vibroengineering, Vol. 21, No. 6, pp. 1600–1612, Sep. 2019, https://doi.org/10.21595/jve.2019.20740

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/jve.2019.20740
UR  - https://doi.org/10.21595/jve.2019.20740
TI  - Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation
T2  - Journal of Vibroengineering
AU  - Zhang, Hongyuan
AU  - Yuan, Huiqun
AU  - Sun, Hongyun
PY  - 2019
DA  - 2019/09/30
PB  - JVE International Ltd.
SP  - 1600-1612
IS  - 6
VL  - 21
SN  - 1392-8716
SN  - 2538-8460
ER  - 

Copy Ris

Copied to clipboard!

@article{Zhang_2019,
	doi = {10.21595/jve.2019.20740},
	url = {https://doi.org/10.21595/jve.2019.20740},
	year = 2019,
	month = {sep},
	publisher = {{JVE} International Ltd.},
	volume = {21},
	number = {6},
	pages = {1600--1612},
	author = {Hongyuan Zhang and Huiqun Yuan and Hongyun Sun},
	title = {Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation},
	journal = {Journal of Vibroengineering}
}

Copy Bibtex

Copied to clipboard!

[1]H. Zhang, H. Yuan, and H. Sun, “Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation,” Journal of Vibroengineering, vol. 21, no. 6, pp. 1600–1612, Sep. 2019, doi: 10.21595/jve.2019.20740.

Copy IEEE

Copied to clipboard!

Zhang, Hongyuan, Huiqun Yuan, and Hongyun Sun. “Forced Response Characteristics of Random Mistuned Blade Disk Based on Weibull Distribution and Monte Carlo Simulation.” Journal of Vibroengineering 21, no. 6 (September 30, 2019): 1600–1612. https://doi.org/10.21595/jve.2019.20740.

Copy Chicago

Copied to clipboard!

Forced response characteristics of random mistuned blade disk based on weibull distribution and monte carlo simulation

Abstract

Highlights

1. Introduction

2. Identification method of random mistuned blade parameters

3. Forced vibration response of mistuned blade disk system

4. Conclusions

References

About this article

Related Articles