International Journal of Rotating MachineryVolume 2012 (2012), Article ID 378565, 10 pagesdoi:10.1155/2012/378565

Research Article

A Continuous Model for the Dynamical Analysis of Mistuned Bladed Rotors

Paolo Bisegna¹ and Giovanni Caruso²

¹Department of Civil Engineering, University of Rome Tor Vergata, Via Politecnico 1, 00133 Rome, Italy²Construction Technologies Institute, National Research Council, Via Salaria Km 29,300, 00016 Monterotondo Stazione, Rome, Italy

Received 26 October 2011; Revised 23 January 2012; Accepted 12 February 2012

Academic Editor: Meinhard Taher Schobeiri

Copyright © 2012 Paolo Bisegna and Giovanni Caruso. 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.

Abstract.

Ideal bladed rotors are rotationally symmetric, as a consequence they exhibit couples of degenerate eigenmodes at coinciding frequencies. When even small imperfections are present destroying the periodicity of the structure (disorder or mistuning), each couple of degenerate eigenfrequencies splits into two distinct values (frequency split) and the corresponding modal shapes exhibit vibration amplitude peaks concentrated around few blades (localization phenomenon). In this paper a continuous model describing the in-plane vibrations of a mistuned bladed rotor is derived via the homogenization theory. Imperfections are accounted for as deviations of the mass and/or stiffness of some blades from the design value; a perturbation approach is adopted in order to investigate the frequency split and localization phenomena arising in the imperfect structure. Numerical simulations show the effectiveness of the proposed model, requiring much lower computational effort than classical finite element schemes.

1. Introduction

Bladed rotors used in turbomachinery possess rotational symmetry provided that no imperfections are present, that is, if all the blades and associated disk sectors are identical to each other. Under this ideal condition, their dynamical behavior is characterized by the presence of couples of degenerate eigenmodes at coinciding frequencies, as shown in [1, 2]. That behavior can be studied in detail by using the finite element technique [1, 3, 4], and the computational cost may be significatively reduced by exploiting the periodicity properties characterizing these structures. In fact, for a given structural eigenmode, all the blades exhibit the same vibration amplitude with a constant phase shift between adjacent ones; thus only a single sector of the structure may be considered in the analysis, by enforcing suitable constraints depending on the admissible phase shifts [1].

As a matter of fact, in practice the symmetry is destroyed by the presence of unavoidable manufacturing defects, material tolerances, or damage arising during service. The loss of symmetry, called disorder or mistuning, may significantly alter the system dynamics even at small disorder level. In particular, the eigenmode degeneracy is removed, and the coinciding modal frequencies of a degenerate eigenmode are split into two distinct values (frequency split phenomenon). Moreover, the vibration localization effect may appear [5–7], consisting of a vibrational-energy confinement in small regions of the rotor, rapidly decaying far away. As a consequence, some blades may vibrate with small amplitude whereas some others with significantly larger amplitude, causing a local stress increase, possibly leading to fatigue failure.

This phenomenon has been extensively studied in the literature: recent works studied the vibration localization arising in bladed rotors by employing the Galerkin method [8, 9], the finite element method [10], or by considering simplified lumped-parameter models [11–15]. Both lumped-parameter models and finite element analysis were employed in [16] to investigate the dependence of the localization phenomenon from the internal damping and interblade coupling in bladed disks. Finite element analysis applied to the study of disordered bladed rotors with a large number of blades is quite accurate but may require a significant amount of computational effort. On the other hand, simplified lumped-parameter models supply satisfying qualitative results, being also computationally more effective than finite element schemes, but may not be able to capture all the peculiarities of the dynamical behavior of a complex structure such as a bladed rotor.

Less effort has been devoted to the derivation of continuous models, which can realistically describe the dynamics of imperfect bladed rotors with low computational effort. For this reason they are particularly useful for performing parametric analyses or designing innovative vibration control schemes. Moreover, continuous models can easily yield analytical results, thus providing a general and comprehensive understanding of the phenomena involved. A continuous model of a bladed rotor without imperfections was considered in [17]. The rotor was composed of grouped blades mounted on a flexible disk. The coupling between adjacent groups of blades due to the disk was accounted for by using lumped springs, whereas distributed springs were used to model the blade stiffness. In [18] the homogenization technique was applied to a simple lumped-parameter model, composed of pendula acted upon by rotational and linear springs, widely employed in the literature to qualitatively study the localization phenomenon in bladed rotors. A continuous model was obtained, yielding analytical closed-form expressions for the eigenfrequencies and the eigenmodes, as well as for the resonance peaks of the forced response, depending on the mistuning level. In [19] a homogenized continuous model of a bladed rotor without imperfections was developed, taking into account the flexibility of both the support disk and the blades by means of the Euler-Bernoulli beam theory. That continuous model was used for the design and optimization of a passive vibration control scheme.

In this paper the model proposed in [19] is suitably generalized to account for the presence of mistuning. In particular, the in-plane vibrations of a mistuned bladed rotor are considered. In Section 2.1 a classical model of a bladed rotor is recalled, employing the Euler-Bernoulli beam theory for the description of both the blades and the support ring. The equations governing that model are solved using the finite element technique, and the obtained numerical results are used as a benchmark for the validation of the homogenized model proposed here. In Section 2.2 the homogenization technique is applied to the Euler-Bernoulli model in order to obtain a continuous homogenized model of the bladed rotor without imperfections. In Section 2.3 the homogenized model previously derived is generalized in order to account for the presence of imperfections, consisting of deviations of the mass and/or stiffness of some blades from the design value. In Section 3.1 analytical expressions of the natural frequencies and vibration modes of a rotor without imperfections are computed using the homogenized model. In order to derive analytical expressions of the same quantities relevant to a rotor with imperfections, a linearization of the proposed homogenized model is performed in Section 3.2, by employing a perturbation technique. As a result, a linear variational formulation of the homogenized model is obtained, supplying analytical expressions for the natural frequencies and vibration mode of a mistuned bladed rotor (Section 3.3). Finally, some simulation results are presented in Section 4, showing the ability of the proposed homogenized model to describe the dynamical behavior of the mistuned bladed rotor introduced in Section 4.1. In particular the results supplied by the homogenized model, obtained by using a numerical approach, and the analytical results supplied by the linearized version of the homogenized model are compared with results supplied by the classical Euler-Bernoulli model, obtained using the finite element technique as described in Section 4.2. The comparison shows good agreement between the finite-element and the homogenized model outcomes, relevant to both the frequency split and the localized vibration modes of the mistuned rotor (Sections 4.4 and 4.5).

2. Homogenized Model

In this section a homogenized model for the analysis of the in-plane vibrations of a mistuned bladed rotor is derived. The rotor is schematically represented in Figure 1; it is made of a linearly elastic material and is composed of blades, of length , clamped on a ring of radius , representing the turbine shaft. The angular spacing between any two adjacent blades is constant and equal to .

Figure 1: Schematic representation of a bladed rotor.

2.1. Euler-Bernoulli Model

In order to develop the homogenized model, a classical model based on the Euler-Bernoulli theory is briefly recalled here, denoted in the foregoing with the acronym EB. It will be also used as a benchmark, after discretization by means of the finite-element method, in order to assess the accuracy of the homogenized model in performing a dynamical analysis of a mistuned rotor. A polar coordinate system is introduced. Let and be the tangential and radial displacement of the ring, respectively; let and be the transversal and axial displacement of the th blade, placed at an angle , respectively, as shown in Figure 1. Both the ring and the blades are assumed to be inextensible, and their flexural behavior is governed by the Euler-Bernoulli theory. The axial strain , the rotation , and the variation of curvature of the cross section of the ring are respectively given by where is the ring radius.

Assuming first that no imperfections are present, and denoting with and , respectively, the linear mass density and in-plane bending stiffness of the ring, and with and , respectively, the linear mass density and bending stiffness of each blade, possibly depending on the radial coordinate , the following Hamiltonian functional describing the dynamical behavior of the structure can be written as follows: subjected to the following constraints: where denotes the time.

Many typologies of imperfections can be accounted for by varying the linear mass density and/or bending stiffness of some blades with respect to the design value; for example, in [15] a crack in a single blade was modeled by reducing the bending stiffness of the damaged blade. Accordingly, if the th blade of the rotor is imperfect, the linear mass density and bending stiffness of that blade can be modified as follows: where and are the variations with respect to the corresponding nominal values.

2.2. Homogenized Model for a Perfect Rotor

The homogenized model developed in [19], valid for a rotor without imperfections, is briefly recalled here in order to illustrate the homogenization procedure. The variational formulation (2) is rewritten for a family of rotors with increasing number of blades, starting from , which is the number of blades exhibited by the reference rotor, in order to approach the homogenization limit . It is assumed that, for each chosen , the blade mass density and bending stiffness rescale as and , respectively, such that the overlined quantities, where is the angular spacing between adjacent blades, remain constant during the homogenization limit. The quantities defined in (5) represent, respectively, the homogenized linear mass density and bending stiffness. These assumptions are substituted into the functional (2) and the homogenization limit is then performed; by extending the functions and , to the annular region and assuming that they converge, together with their derivatives, to functions and , the functional (2) converges, for , to [19] subjected to the following constraints: In order to analyze the rotor dynamics, a sinusoidal time dependence of the unknown is assumed, that is, where is the imaginary unit. Enforcing this position and recasting the functional (6) into a variational form, the following is obtained: subjected to the following constraints: where indicates a test function. The homogenized model described here is denoted with the acronym HOM in the foregoing.

2.3. Homogenized Model for an Imperfect Rotor

The presence of imperfections can be accounted for in the homogenized model by letting the homogenized linear mass density and bending stiffness depend also on the angular variable . Accordingly, it is assumed that where and are suitable perturbations superimposed to the homogenized design values of linear mass density and bending stiffness , respectively. In order to choose and , it is necessary to specify how the imperfect blades behave during the homogenization limit. For the sake of example, it is assumed that the th blade of the rotor is imperfect, with increased linear mass density equal to as indicated in Figure 2. When the number of blades is increased to approach the homogenization limit, the number of imperfect blades increases as depicted in Figure 2, and their mass density and bending stiffness are rescaled according to the conditions given in Section 2.2. In particular, at the th homogenization-process step, the total number of blades becomes , and the blade linear mass density and bending stiffness are rescaled by a factor with respect to their nominal value; on the other hand, the number of imperfect blades becomes , with rescaled linear mass density variation equal to .

(a)

(b)

Figure 2: Homogenization limit for an imperfect rotor; imperfection due to an increase of linear mass density of a blade. (a) Reference rotor with blades and imperfection on the th blade, in black color. (b) Homogenization step process ; rotor with blades, with rescaled bending stiffness and linear mass density.

Under these assumptions, in the homogenization limit , the homogenized linear mass density perturbation and the homogenized bending stiffness perturbation are piecewise constant functions having the following expressions: where The expression (12) can be straightforwardly generalized to the case of multiple imperfect blades.

3. Analytical Solutions for Eigenfrequencies and Eigenmodes

In this section explicit formulas for the eigenfrequencies and eigenmodes of a bladed rotor, with or without imperfections, are reported. They are obtained by using the homogenized model developed in Sections 2.2 and 2.3.

3.1. Perfect Rotor

In order to evaluate the eigenmodes and eigenfrequencies of a perfect rotor by using the homogenized model proposed here, it is assumed for the sake of simplicity that and do not depend on the radial variable ; under this simplifying assumption, using the localization lemma, (9) and (10) yields the field equilibrium equation: and the boundary conditions: Due to the rotational symmetry of the problem, the eigenmodes depend harmonically on the angular variable . For each fixed nodal diameter, denoted by the modal number , a countable set of eigenfrequencies and eigenmodes, denoted by the modal number , can be computed as follows. Fixed , a solution of the field equation (14), is given by where is an arbitrary constant, and the modal eigenfrequency is related to by the relation In (16) the quantities , , are scalar unknowns; by substituting (16) into the four boundary conditions (15), an eigenvalue problem is obtained, exhibiting a countable set of eigenvalues and corresponding eigenvectors , . The phase orientation is left undetermined after solving (14) and (15); as a consequence, for each fixed and two independent eigenmodes exist, at the same modal frequency . These are known in the literature as degenerate eigenmodes [1, 2]. All the eigenmodes relevant to are nondegenerate.

3.2. Linearized Homogenized Model for an Imperfect Rotor

In order to derive analytical expressions for the modal eigenfrequencies and eigenmodes of the rotor in the presence of imperfections, the homogenized model proposed in Section 2.3 is linearized by employing a perturbation technique. Accordingly, the eigenmodes of the imperfect rotor are computed by perturbing the eigenmodes relevant to the perfect rotor, that is, for each fixed couple of modal indices , the eigenmodes of the imperfect rotor are given by where is the eigenmode of the perfect rotor given in (16) and is an unknown perturbation. Moreover the relevant eigenfrequency is given by where is the frequency shift due to the presence of imperfections. By substituting the expansions (11), (19), and (20) into the weak formulation (9) and retaining only the first order terms, the following linearized weak formulation is obtained: under the following constraints: where it is assumed, for simplicity, that and do not depend on the radial variable . The unknowns appearing in (21) are the scalar quantities and and the function . Due to the presence of the imperfections, the indeterminacy on the phase angle is removed and, for each couple of modal indexes , two set of solutions are expected, corresponding to the two split eigenmodes belonging to the same degenerate eigenmode relevant to the perfect structure. This linearized homogenized model is denoted in the foregoing with the acronym L-HOM.

3.3. Eigenfrequencies and Eigenmodes Evaluation for an Imperfect Rotor

By considering the linearized variational formulation (21), it is observed that, for each fixed modal number , with , , and solutions of (14) and (15), are two independent functions belonging to the kernel of the self-adjoint operator at the left hand side of (21). Accordingly, by substituting first and then in (21), together with the expression of and given by (16) and (18), respectively, two independent equations in the unknowns and are obtained, whose solutions supply the two couples of frequency split and phase orientation relevant to the perturbed mode . Once the values of and have been computed, the corresponding modal shapes can be obtained by making use of the linear weak formulation (21). To this end (21) is integrated by parts with respect to the variable and, making use of (22)₂, it yields Another weak-form equation is obtained by multiplying (22)₁ by and integrating over , yielding In order to find out the unknown modal perturbation , a spectral representation is used by setting

The representation (26) of is then substituted in (24), and the test function is chosen according to the following expressions where is any positive fixed integer and is an arbitrary function of . Using the localization lemma in the first integral of (24), an explicit expression for and , , is obtained. The expression (26) for , where now and are known scalar quantities, is then substituted into (24) and (25). Linear equations in the unknown coefficients and are obtained, for each fixed positive integer , by choosing the test function and , at and , as and and using again localization lemma. Their solution provides the Fourier coefficients of the unknown modal perturbation .

4. Numerical Simulations

In this section, some numerical simulation results are presented in order to show the effectiveness of the proposed homogenized model in the investigation of the frequency split and localization phenomenon in imperfect rotors.

4.1. Case Study Bladed Rotor

A case study problem is introduced here; a bladed rotor, similar to the one schematically represented in Figure 1, is considered. It is composed of 32 elastic blades of cross-section mm and length mm; the blades are clamped to a support ring of radius mm and cross section mm. The ring and the blades are composed of steel, with Young modulus GPa and mass density Kg/m³. The degenerate modal frequencies of this bladed rotor are reported in Table 1, relevant to modal numbers and , and evaluated according to the EB model. The eigenmodes relevant to are not considered since they do not involve any deformation in the support ring, in fact for that modes [19].

Table 1: Modal frequencies in rad/s, relevant to the perfect rotor and evaluated according to the EB model. Modal frequencies relevant to modal numbers and .

27.0746
27.0746
27.9718
27.9718
28.4242
28.4242
28.6596
28.6596

4.2. Numerical Algorithms

In order to find a stationary solution for the Hamiltonian functional (2), describing the dynamical behavior of the case study bladed rotor according to the EB model, a finite-element scheme is adopted here. Accordingly, (2) is discretized by using two-node beam elements and Hermite polynomials as interpolation scheme, both for the ring and the blades. The axial inextensibility constraint is enforced by using a penalization method. A discrete formulation for the dynamical problem is finally obtained by assuming as unknowns the nodal displacements and rotations. The results provided by the latter discretized model are used in the foregoing as benchmark to assess the accuracy of the homogenized model.

The homogenized weak formulation in (9), defined in the annular region , is also discretized in the angular direction with the Ritz-Rayleigh method assuming, as shape functions, whereas a finite element discretization is assumed in the radial direction, employing two-node beam elements with Hermite polynomials as interpolation scheme. Accordingly, the nodal unknowns are the real and imaginary part of the complex Fourier coefficients of the tangential displacement and its radial derivative.

4.3. Validation of the Homogenization Limit

In this section, the homogenization limit is numerically studied, and a convergence analysis is presented. The modal eigenfrequencies of the perfect rotor described in Section 4.1 are computed increasing the number of blades, starting from the reference configuration with 32 blades and rescaling at the same time the linear mass density and the bending stiffness of the blades as described in Section 2, in order to approach the homogenization limit relevant to the HOM model. In particular, , , , and have been considered in the computations. In Figure 3 the relative error between degenerate eigenfrequencies evaluated according to the reference EB model and the HOM model is reported, relevant to modal numbers and . The figure shows that a quadratic convergence, indicated by the triangle, is achieved. The relative error slightly increases with the increase of . The homogenized model turns out to be quite accurate in describing the rotor dynamics, as the relative error is very small even in the reference case of 32 blades.

Figure 3: Convergence analysis for the homogenization limit: relative error between modal frequencies evaluated according to the EB model and the HOM model, as a function of the number of blades. Rotor without imperfections. Modal frequencies relevant to modal numbers and — □  ; – ; – –   ; . The triangle in the figure indicates the quadratic convergence slope.

4.4. Modal Frequency Analysis

In this section the modal frequencies relevant to an imperfect bladed rotor are analyzed. The imperfection is introduced by increasing the linear mass density of the fourth blade of the perfect rotor described in Section 4.1 by a factor . The modal frequencies , relevant to modal numbers and , are evaluated according to the reference EB model and the HOM model, considering different values of . In Table 2 the modal frequencies evaluated according to the EB model are reported.

Table 2: Modal frequencies in rad/s, relevant to the imperfect rotor and evaluated according to the EB model. The imperfection is introduced by increasing the mass of the fourth blade of the perfect rotor described in Section 4.1 by , , , and 1%, with respect to its nominal value. Modal frequencies relevant to modal numbers and .

27.0742
27.0746
27.9713
27.9718
28.4237
28.4242
28.6591
28.6596

27.0738
27.0746
27.9709
27.9718
28.4233
28.4242
28.6587
28.6596

27.0703
27.0746
27.9671
27.9718
28.4193
28.4242
28.6544
28.6596

27.0656
27.0746
27.9619
27.9718
28.4134
28.4242
28.6476
28.6596

From Table 2 it can be seen that, due to the presence of imperfections, a frequency split occurs between frequencies belonging to the same modal number , which were coinciding in the perfect case as shown in Table 1. The difference between eigenfrequencies relevant to the imperfect rotor in Table 2 and the corresponding ones relevant to the perfect rotor in Table 1 is the frequency shift , which occurs due to the presence of imperfections. It can be noticed that, for each couple of modal numbers , one of the two frequency shifts is vanishing; in fact, for each fixed , one of the two eigenfrequencies of the imperfect rotor in Table 2 coincides with the corresponding frequency of the perfect rotor in Table 1. In particular, this eigenfrequency corresponds to a modal shape antisymmetric with respect to the location of the imperfect blade [18]. Table 3 contains the relative error between the frequency shifts evaluated according to the EB model and to the HOM model. For each couple of modal numbers , only the value relevant to the nonvanishing frequency shift is reported. The errors are almost independent from the imperfection level , and they slightly increase with the increase of , ranging from when to when . Finally, in Table 4 the relative errors between frequency shifts evaluated according to the HOM model and the L-HOM model are reported. It is remarked here that while the results supplied by the HOM model are computed by using the numerical procedure described in Section 4.2, the results relevant to the L-HOM model are readily obtained using the analytical expressions developed in Section 3.3. For each couple of modal numbers , only one value is reported, relevant to the nonvanishing frequency shift. Results in Table 4 show that the linearized model is suitable for accurately evaluating frequency shifts of an imperfect rotor when imperfections are sufficiently small; for the relative errors range from 0.3% to 2.5%, increasing with mode number . The relative difference between frequency shifts becomes larger for larger values of together with higher values of mode number ; for an increased accuracy a refined theory may be used, which can be derived by retaining also the higher-order terms in the perturbation expansion performed in Section 3.2. Accordingly, a nonlinear theory would be obtained, whose solution would require a numerical iterative procedure.

Table 3: Relative error between frequency shifts relevant to the imperfect rotor, evaluated according to the HOM model and the EB model. The imperfection is introduced by increasing the mass of the fourth blade of the perfect rotor described in Section 4.1 by , , , and 1%, with respect to its nominal value. Modal frequencies relevant to modal numbers and .

0.0114
0.0281
0.0513
0.0824

0.0116
0.0283
0.0519
0.0834

0.0125
0.0299
0.0556
0.0906

0.0103
0.0261
0.0500
0.0828

Table 4: Relative error between frequency shifts relevant to the imperfect rotor, evaluated according to the HOM model and the L-HOM model. The imperfection is introduced by increasing the mass of the fourth blade of the perfect rotor described in Section 4.1 by , , , and 1%, with respect to its nominal value. Modal frequencies relevant to modal numbers and .

0.0032
0.0054
0.0087
0.0125

0.0064
0.0108
0.0173
0.0250

0.0319
0.0539
0.0860
0.1233

0.0665
0.1122
0.1777
0.2507

4.5. Eigenmode Analysis

In this section the eigenmodes of the rotor described in Section 4.1, made imperfect by increasing the linear mass density of its 4th blade by , are evaluated. For the sake of comparison, both the HOM model and the L-HOM model are employed in the analysis, and their results are compared with the eigenmodes computed by using the finite-element formulation of the EB model, used as a benchmark. In Figure 4 the eigenmodes relevant to mode numbers and are reported; in particular, for each couple of modal number only the eigenmode corresponding to a nonvanishing frequency shift exhibits a vibration localization and is reported in the figure. All the modal shapes have been normalized using the same procedure: in particular, the modulus of the complex Fourier coefficient relevant to the dominant mode (i.e., the th + 1 Fourier coefficient for a eigenmode of mode number ) has been set to 1. Finally, the eigenmodes have been all rescaled by a factor of 0.2 to be clearly shown in the figure.

(a)

(b)

(c)

(d)

Figure 4: Eigenmodes of the imperfect rotor, evaluated according to the HOM model, the L-HOM model, and the EB model. The imperfection is introduced by increasing the mass of the fourth blade of the perfect rotor described in Section 4.1 by , with respect to its nominal value. – – green: perfect rotor, — red: imperfect rotor using HOM model, blue: imperfect rotor using L-HOM model, – · – black: imperfect rotor using the EB model.

The homogenized model turns out to be suitable for studying the localization phenomenon in imperfect rotors: in fact, as shown in Figure 4, the localization effect due to the imperfection placed on the fourth blade clearly appears, increasing with the increase of the modal number . It turns out that the HOM model and the L-HOM model are both in good agreement with the reference EB model, with the L-HOM model slightly underestimating the modal amplitude around the imperfect blade with respect to the HOM model.

5. Conclusions

In this paper a homogenized model has been proposed, suitable for the dynamical analysis of mistuned bladed rotors. The model seems to be a good compromise between accuracy and simplicity; indeed, it considers the distributed elasticity of the blades and the coupling between a single blade and the other ones, due to the presence of a flexible support disk. On the other end, it is a continuous model, and thus it is computationally more effective than classical finite element models. Numerical simulations have been reported, showing the capability of the model in evaluating the frequency split and vibration localization phenomenon arising in the presence of mistuning. The proposed model may be an effective tool for the design and the parametric analysis of bladed rotors.

References

D. L. Thomas, “Dynamics of rotationally periodic structures,” International Journal for Numerical Methods in Engineering, vol. 14, no. 1, pp. 81–102, 1979.

I. Y. Shen, “Vibration of rotationally periodic structures,” Journal of Sound and Vibration, vol. 172, no. 4, pp. 459–470, 1994.

A. J. Fricker and S. Potter, “Transient forced vibration of rotationally periodic structure,” International Journal for Numerical Methods in Engineering, vol. 17, no. 7, pp. 957–974, 1981.

C. W. Cai, Y. K. Cheung, and H. C. Chan, “Uncoupling of dynamic equations for periodic structures,” Journal of Sound and Vibration, vol. 139, no. 2, pp. 253–263, 1990.

D. J. Ewins, “The effects of detuning upon the forced vibrations of bladed disks,” Journal of Sound and Vibration, vol. 9, no. 1, pp. 65–79, 1969.

A. Sinha, “Calculating the statistics of forced response of a mistuned bladed disk assembly,” AIAA Journal, vol. 24, no. 11, pp. 1797–1801, 1986.

C. Pierre and E. H. Dowell, “Localization of vibrations by structural irregularity,” Journal of Sound and Vibration, vol. 114, no. 3, pp. 549–564, 1987.

B. W. Huang, “Effect of number of blades and distribution of cracks on vibration localization in a cracked pre-twisted blade system,” International Journal of Mechanical Sciences, vol. 48, no. 1, pp. 1–10, 2006.

Y. J. Chiu and S. C. Huang, “The influence on coupling vibration of a rotor system due to a mistuned blade length,” International Journal of Mechanical Sciences, vol. 49, no. 4, pp. 522–532, 2007.

E. P. Petrov and D. J. Ewins, “Analysis of the worst mistuning patterns in bladed disk assemblies,” Journal of Turbomachinery, vol. 125, no. 4, pp. 623–631, 2003.

M. P. Mignolet, W. Hu, and I. Jadic, “On the forced response of harmonically and partially mistuned bladed disks. Part I: harmonic mistuning,” International Journal of Rotating Machinery, vol. 6, no. 1, pp. 29–41, 2000.

M. P. Mignolet, W. Hu, and I. Jadic, “On the forced response of harmonically and partially mistuned bladed disks. Part II: partial mistuning and applications,” International Journal of Rotating Machinery, vol. 6, no. 1, pp. 43–56, 2000.

J. Tang and K. W. Wang, “Vibration delocalization of nearly periodic structures using coupled piezoelectric networks,” Journal of Vibration and Acoustics, vol. 125, no. 1, pp. 95–108, 2003.

H. H. Yoo, J. Y. Kim, and D. J. Inman, “Vibration localization of simplified mistuned cyclic structures undertaking external harmonic force,” Journal of Sound and Vibration, vol. 261, no. 5, pp. 859–870, 2003.

X. Fang, J. Tang, E. Jordan, and K. D. Murphy, “Crack induced vibration localization in simplified bladed-disk structures,” Journal of Sound and Vibration, vol. 291, no. 1-2, pp. 395–418, 2006.

Y. J. Chan and D. J. Ewins, “Management of the variability of vibration response levels in mistuned bladed discs using robust design concepts. Part 1: parameter design,” Mechanical Systems and Signal Processing, vol. 24, no. 8, pp. 2777–2791, 2010.

L. F. Wagner and J. H. Griffin, “A continuous analog model for grouped-blade vibration,” Journal of Sound and Vibration, vol. 165, no. 3, pp. 421–438, 1993.

P. Bisegna and G. Caruso, “Dynamical behavior of disordered rotationally periodic structures: a homogenization approach,” Journal of Sound and Vibration, vol. 330, no. 11, pp. 2608–2627, 2011.

P. Bisegna and G. Caruso, “Optimization of a passive vibration control scheme acting on a bladed rotor using an homogenized model,” Structural and Multidisciplinary Optimization, vol. 39, no. 6, pp. 625–636, 2009.