Next Article in Journal
Adopting New Machine Learning Approaches on Cox’s Partial Likelihood Parameter Estimation for Predictive Maintenance Decisions
Previous Article in Journal
Analysis of the Control Characteristics of the Electro-Hydraulic Vibration System Based on the Single-Neuron Control Algorithm
Previous Article in Special Issue
Vibro-Acoustic Performance of a Fluid-Loaded Periodic Locally Resonant Plate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Dynamics and Combination Resonance of a Flexible Turbine Blade with Contact and Friction of Shrouds

1
China Oil & Gas Pipeline Network Corporation, Beijing 100020, China
2
National Key Laboratory of Helicopter Aeromechanics, China Helicopter Research and Development Institute, Jingdezhen 333001, China
3
Department of Engineering Mechanics, Dalian University of Technology, Dalian 116024, China
4
State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment, Dalian 116023, China
5
Department of Theoretical and Applied Mechanics, Belarus State University, 220030 Minsk, Belarus
*
Author to whom correspondence should be addressed.
Machines 2024, 12(1), 59; https://doi.org/10.3390/machines12010059
Submission received: 18 November 2023 / Revised: 7 January 2024 / Accepted: 8 January 2024 / Published: 12 January 2024
(This article belongs to the Special Issue Advanced Dynamic Analysis and Vibro-Acoustic Control Methods)

Abstract

:
Flexible shrouded blades are commonly adopted in the last stages of steam turbines where complicated dynamical behavior can be induced by dry friction force generated on contacting interfaces between adjacent shrouds and the geometric nonlinearity due to the structural flexibility of the blades. In this paper, combination resonance caused by contact and friction forces generated on shroud interfaces is investigated, which is a concurrence of 1:3 internal resonance involving the first and second modes in the flapwise direction and the primary resonance of the first flapwise mode. The stiffness and damping at the contact interface are obtained by linearizing the contact and friction forces between shrouds through the harmonic balance method. The vibrating blade is modeled as a beam with a concentrated mass of which the responses under the combination resonance are solved through the multiple-scale method. Sensitivities of response with respect to the angle of shrouds, contact stiffness and rotation speed are illustrated, and the influences of these parameters on the periodicity and amplitudes of steady responses are demonstrated. The parametric regions where the combination resonance occurs are pointed out. Finally, parametric analyses are presented to show how the amplitude–frequency relation of the multiple-scale solutions under the combination resonance vary with detuning and design parameters. The present research provides a designing basis for improving the dynamic performance of flexible shrouded blades and suppressing vibrations of blades by adjusting structural parameters in practical engineering.

1. Introduction

Shrouds have been widely used by developers of multi-stage steam turbines to dissipate vibratory energy and reduce the stress of rotating bladed blisks for the purpose of improving their structural reliability [1]. To ensure the constant flow rate of the steam, the blades mounted on the last-stage blisk are generally designed to be the longest with considerable structural flexibility, which causes an effect of geometric nonlinearity that cannot be ignored. On the other hand, there exists contact and dry friction on the contacting interfaces between adjacent shrouds, which intensifies the nonlinear characteristics as far as the dynamic of the shrouded blades is concerned.
Various studies were published previously devoted to modeling the shrouded blade of turbines, covering the topics of modeling of blades as well as the contact force between adjacent shrouds. Single-blade modeling was adopted to understand the dynamical properties of flexible turbine blades [2,3]. In some publications, the interaction between the reference blade and the adjacent blades [4,5,6] was considered, while a few studies were aimed at the vibration of the whole blade disk [7]. Zhou et al. [8] investigated the transverse vibration of composite blades of a wind turbine using a generalized Timoshenko beam. Slender blades that resembled the beam model and a rigid disk were integrated into a continuous model by Shadmani [9] in modal analysis of a turbine blisk. As for the contacts between adjacent shrouds, various models were published for evaluating the inter-shroud action of friction and contact. The effects of macro- and micro-slip between the neighboring blades were established by Griffin [10] and later publications such as [11,12]. Later, Yang and Menq established a special model to predict the resonance of structural response possessing a three-dimensional frictional constraint [13]. Nan et al. [14] translated the contact and friction that acted on the shroud interfaces into an equivalent stiffness as well as damping by means of the harmonic balance method.
In essence, the contact between the shrouds is physically nonlinear and has been regularly simplified as a piecewise linear function of the blade motion. In the view of dynamics, such contact force is attributable to complicated behaviors, e.g., self-excited vibrations, bifurcation of motion, chaos and switch of stability [15]. He et al. [16] proposed a full set of blades subjected to friction and impact forces due to shroud contacting and demonstrated various types of periodic, quasi-periodic and even chaotic motions of the blades. He et al. [17] studied the interaction of bending and torsional vibration of blades, considering the rub and impact between shrouds. As for the nonlinear dynamics of full blisks, achievements were made in mistuned disks, nonlinear resonance frequencies and localization phenomena [18,19,20,21,22,23].
For the combination resonance of the turbine blade, most of the existing research attended to the internal resonance [24,25,26], primary resonance [27,28] or combination resonance with both internal and sub- or super-harmonic responses [29,30,31]. A few publications were dedicated to the combination resonance consisting of the internal and primary resonances of the flexible turbine blade. Li et al. [32] investigated the nonlinear vibration of a blade simultaneously undergoing two resonances (i.e., internal as well as primary), considering aerodynamic force, damping and structural nonlinearity. Yuan and Wang studied internal resonance, primary resonance and the combined multiple modal resonances of a blade [33], where they demonstrated the contribution of primary and internal resonances to the combination resonance in terms of frequencies and energy. It is worth pointing out that, despite the existent works on the nonlinear vibrations of bladed blisks, the combination resonance induced by contact and friction between neighboring blades has yet to be thoroughly understood. Such understanding will provide a better understanding of the dynamical behaviors of flexible shrouded blades and how shroud design influences the blade vibration under the situation of combination resonance.
In the present study, the dynamics of a flexible blade mounted with a shroud are investigated, focusing particularly on the combination resonance of blade vibration that is induced by the contact as well as friction forces applied on the shroud. Firstly, the equivalent stiffness and damping of the shroud due to the normal contact as well as the friction are obtained through the method of harmonic balance, and the blade is simplified as one single continuous beam attached by a concentrated mass. A steady response under the combination resonance is presented through the multiple-scale method. Afterward, influences of design parameters on periodicities and amplitudes are demonstrated for the steady responses. Jumping of amplitude due to the primary resonance and the transfer of vibration energy between resonance modes excited by the internal resonance are pointed out. Finally, amplitude–frequency curves of the multiple-scale solutions under the combination resonance are derived. Parametric analyses are presented to show the effect of detuning and key design parameters on the first and second flapwise modes.

2. Governing Equations of Flexible Shrouded Blade

Figure 1 depicts the sketch of a perfectly tuned, rotating shrouded disk that is commonly used in steam turbines. For such a blisk, all of the shrouded blades are mounted circumferentially and separated in equal space [34]. Hence, the blisk can be considered to have a cyclic symmetry. In addition, the blades of the final stage of the turbine generally have a large span size to ensure the constant flow rate of the working medium, which will enhance the flexibility of the blade.
To focus on the overall dynamics of the blade, a macro-slip model is adopted in the present work. The dimension of the shroud mounted at the tip of the blade is small compared with the blade length and is generally rigid as opposed to the long and flexible blade, which can be further treated as a single mass point. The normal contact force between the adjacent shrouds is modeled by a linear spring [16]. Accordingly, the flexible shrouded blade is modeled as a beam with an attached mass that is placed at the tip of the blade, as shown in Figure 2. A Cartesian frame can be set up for the blade with the origin placed at the center of gravity of the root and three axes of x (axial direction), y (flapwise direction) and z (edgewise direction). α represents the tilting angle of the shroud, and is the rotation speed. L is the length of the blade.  k s  and  δ s  are contact (normal) stiffness and the circumferential gap between shrouds, respectively.

2.1. Dynamic Model of Contact between Adjacent Shrouds

Prior to deriving the governing equation of the blade motion, it is necessary to describe the contact action between adjacent shrouds, which is a contribution to the nonlinear dynamics and combination resonance. Hereafter, the equivalent stiffness and damping of the normal and friction between adjacent shrouds are established through the method of harmonic balance.
In this research, the contact force is presumed to change as a piecewise linear function of the relative displacement of two adjacent shrouds, whereas friction carries typical nonlinear characteristics. The states of contact (e.g., separation, stick or slip) depend upon relative displacements in both normal and tangential directions.
In Figure 3, the non-conservative forces on the shroud are aerodynamic force (Q), normal contact force ( N s ) and friction ( F s ), where
Q = Q 0 sin l Ω t ,
with  Q 0  being the amplitude and l the number of guided vanes of the upstream inlet [27].
The flapwise displacement of the blade, v, can be assumed as [14]
v = B cos l Ω t ϕ = B cos ϑ ,
where B and ϕ stand for the amplitude and phase angle of displacement v, respectively. Then, the contact force between two adjacent shrouds is modeled as
N s = k s v δ s sin α , v > δ s 0 , δ s v δ s k s v + δ s sin α v < δ s
As for the friction force between two neighboring shrouds, there are two possible situations during one cycle of the vibration based on the change in the normal force when the macro-slip model is adopted: (1)  N s  is not zero and time varies during shroud contacting, and (2)  N s  is zero when two shrouds separate. The hysteretic loop of the friction force  F s  changes along with the tangential relative displacement are depicted in Figure 4. Different from the hysteretic loop of the friction in [14], the influence of the normal force on sliding friction is considered in this paper.
As depicted in Figure 4, the friction force between shrouds starts from point A ( ϑ = 0 ) and changes clockwisely according to the hysteresis loop. Section AB corresponds to the phase of static friction. Then, from point B ( ϑ = ϑ 1 ) to point C ( ϑ = ϑ 2 ), the blade begins to slide towards the opposite direction. The segment of BC corresponds to the phase of sliding friction. The separation of adjacent shrouds takes place at point C, and the friction declines to zero at D ( ϑ = ϑ 2 ). Subsequently,  v t  continues to decrease, and the state of friction moves from D to D′ ( ϑ = π ϑ 2 ) where the blade starts to contact the blade from the other side. Section D′E′ also corresponds to the phase of static friction, where the blade begins to slide at point E′ ( ϑ = ϑ 3 ). After that, the negative  v t  reaches the maximum at A′ ( ϑ = π ). Notably, the section A′A in the loop is symmetric with section AA′.
Based on the aforementioned hysteretic behavior, the friction force can be modeled in terms of the state of contact between adjacent blades:
F s = k t v t B cos α + μ k s v δ s sin α , 0 ϑ < ϑ 1 , ( A B ) μ k s v δ s sin α , ϑ 1 ϑ < ϑ 2 , ( B C ) 0 , ϑ 2 ϑ < π ϑ 2 , ( D D ) k t v t + δ s cos α , π ϑ 2 ϑ < ϑ 3 , ( D E ) μ k s v + δ s sin α , ϑ 3 ϑ < π , ( E A ) k t v t + B cos α + μ k s v + δ s sin α , π ϑ < π + ϑ 1 , ( A B ) μ k s v + δ s sin α , π + ϑ 1 ϑ < π + ϑ 2 , ( B C ) 0 , π + ϑ 2 ϑ < 2 π ϑ 2 , ( D D ) k t v t δ s cos α , 2 π ϑ 2 ϑ < π + ϑ 3 , ( D E ) μ k s v δ s sin α , π + ϑ 3 ϑ < 2 π , ( E A ) ,
where  k t  and μ are shear stiffness and the friction coefficient of the contact surface of the shroud.  ϑ 1 ϑ 2  and  ϑ 3  are the phase angles of the friction force, defined as
k t = 2 k s 1 ν 2 ν
ϑ 1 = cos 1 k t B + 2 μ k s δ s tan α k t B + 2 μ k s B tan α , ϑ 2 = cos 1 δ s B ,   ϑ 3 = cos 1 k t δ s + μ k s δ s tan α k t B + μ k s B tan α
To facilitate the dynamical analysis in the present study, the method of harmonic balance (HBM) [35] is used to simplify contact (normal) and friction forces to obtain the stiffness and damping at the shroud surface. For this purpose, one linearizes  N s  as
N s = N k s cos ϑ + N cs sin ϑ = k n v + c n d v d t ,
where
k n = k s sin α B π 2 B ϑ 2 + B sin 2 ϑ 2 4 δ s sin ϑ 2 .
Based on Equation (3), the contribution of the blade velocity to the contact force is negligible; hence, one assumes
c n = 0
The friction force between adjacent shrouds can be expressed as
F s = F k s cos ϑ + F cs sin ϑ = k f v + c f d v d t ,
where  k f  and  c f  represent the stiffness and damping of  F s , respectively.
k f = 1 B π k t B cos α ϑ 3 + ϑ 2 + ϑ 1 π + sin 2 ϑ 3 + sin 2 ϑ 2 + sin 2 ϑ 1 2 2 sin ϑ 1 + 2 k t δ s cos α sin ϑ 3 sin ϑ 2 + μ k s B sin α ϑ 1 + B sin 2 ϑ 1 2 2 μ k s δ s sin α sin ϑ 1
c f = 1 l Ω B π 2 k t B cos α cos ϑ 1 1 + 1 cos 2 ϑ 1 + cos 2 ϑ 2 cos 2 ϑ 3 4 2 k t δ s cos α cos ϑ 3 + cos ϑ 2 + 2 μ k s δ s sin α cos ϑ 1 cos ϑ 3 1 μ k s B sin α cos 2 ϑ 1 1 2

2.2. Governing Equations of Blade Vibration

The blade in the present study is modeled as a beam attached by a concentrated shroud mass at the tip of the blade. For blades carrying small installation angles, their torsional rigidities are sufficient; hence, the twisting movement of the rotating blade can be effectively neglected. Thus, the edgewise deformation of the blade caused by the torsion should be neglected. Additionally, the dynamic spanwise (axial) displacement can be described as minor as opposed to the pre-stretching caused by centrifugal load on the blade [32]. Therefore, the interaction between the flapwise and axial motions is fulfilled by means of the centrifugal stiffening upon lateral vibration frequencies. Based on Figure 2, the equation of the blade flapwise vibration is established, which reads
m l + m s δ x L v ¨ + c v v ˙ m l v + m s δ x L Ω 2 v + E I v E A u v 1 2 E A v 3 = Q + N s sin α + F s cos α δ x L ,
where v = v(x,t) is the dynamic displacement of the blade in the flapwise. u = u(x) is the spanwise (axial) displacement generated by centrifugal force of the blade. The EI and EA are bending and tensile stiffnesses of the blade, respectively.  m l  represents the linear mass density of the blade.  m s  is the mass of shroud. The term E A v 3 is the outcome of geometric nonlinearity due to the structural flexibility of the blade, and  δ  is the Dirac function of x. The overhead dot and prime represent the partial derivatives of a physical quantity with respect to time t and x, respectively, i.e.,
( · ) ˙ ( · ) / t ,   a n d ( · ) ( · ) / x
Subsequently, the flapwise displacement of the blade is decomposed in the modal space
v x , t = k γ k x q k t ,
where  γ k  stands for the mode shape of order k of the beam under given boundary conditions in the flapwise direction, and  q k  are the generalized (modal) coordinates of  γ k . Thus, Equation (13) can be simplified using the Galerkin method considering the first- and second-order modal coordinates, as
q ¨ k + c k q ˙ k + ω k 2 q k + g k , 1 q 1 3 + 3 g k , 2 q 1 2 q 2 + 3 g k , 3 q 1 q 2 2 + g k , 4 q 2 3 = f k sin l Ω t , k = 1,2 ,
where  c k ω k g k , j  and  f k  are defined in Equations (A1)–(A8), Appendix A. Further, Equation (16) can be rewritten with the introductions of non-dimensional notations:
q ¯ i = q i δ s   a n d   τ = Ω t ,
such that one has the compact governing equations
q ¯ ¨ k + c ¯ k , k q ¯ ˙ k + ω ¯ k , k 2 q ¯ k + g ¯ k , 1 q ¯ 1 3 + 3 g ¯ k , 2 q ¯ 1 2 q ¯ 2 + 3 g ¯ k , 3 q ¯ 1 q ¯ 2 2 + g ¯ k , 4 q ¯ 2 3 = f ¯ k sin l τ , k = 1,2 ,
where  c ¯ k ω ¯ k g ¯ k , j  and  f ¯ k  are defined in Equations (A9)–(A15), Appendix A.

3. Multiple-Scale Solutions under Combination Resonance

The multiple-scale method will be adopted in this section for solving Equation (18). To construct the analytical solution of the problem, a small, artificial bookkeeping parameter is introduced to order the quantities of dynamic displacements, damping and aerodynamic forces before they are inserted into Equation (18). Thus, one has
q ¯ k = ε q ~ k , c ¯ k , k = ε c ~ k , k , f ¯ k = ε 3 f ~ k , k = 1,2
Herein, the regular scale of time is replaced with three temporal scales, i.e.,
T 0 = t , T 1 = ε t , T 2 = ε 2 t
Then,  q ˜ k  is expanded using the power series of  ε  and the new scales
q ~ k = j = 0 2 ε j q ~ k , j T 0 , T 1 , T 2 , k = 1,2
Let  D j = / T j  and substitute Equations (19) and (21) back to Equation (18), then equate their same orders of  ε . In doing so, one obtains
order   ε   D 0 2 q ~ k , 0 + ω ¯ k 2 q ~ k , 0 = 0
order   ε 2   D 0 2 q ~ k , 1 + ω ¯ k 2 q ~ k , 1 = 2 D 0 D 1 q ~ k , 0 c ~ k , k D 0 q ~ k , 0
order   ε 3   D 0 2 q ~ k , 2 + ω ¯ k 2 q ~ k , 2 = F k q ~ j , 0 , q ~ j , 1 , l ,
where  k , j = 1,2 .
For Equation (22), let us assume that
q ~ k , 0 = A k e i ω ¯ k T 0 + c o m p l e x   c o n j u g a t e , k = 1,2
and then substitute these back to Equation (23) to obtain
2 D 1 A k + c ~ k , k A k = 0 , k = 1,2 ,
which further leads to
q ~ k , 1 = A k e i ω ¯ k T 0 + c o m p l e x   c o n j u g a t e , k = 1,2
The highest  q ~ k  can be evaluated through the perpetual term of Equation (24).
A particular type of combination resonance is of present interest, i.e., when a particular internal resonance takes places with  ω ¯ 2 = 3 ω ¯ 1  (“1:3”). To be specific, the frequencies of the first-order flap mode are triple the second-order flapwise frequency, and alongside, there is a primary resonance of the first flapwise mode (i.e.,  ω ¯ 1 = l ). Under this situation, one can assume that  ε A k = 1 / 2 α k e i β k . By introducing new phases  γ 1 = σ 1 T 1 β 1 γ 2 = β 2 3 β 1 σ 2 T 1 , one obtains the equations of motion for the amplitudes of  ε A k
α ˙ 1 = c 1,1 2 α 1 3 g 1,2 8 ω ¯ 1 α 1 2 α 2 sin γ 2 f 1 2 ω ¯ 1 cos γ 1
α ˙ 2 = c 2,2 2 α 2 + g 2,1 8 ω ¯ 2 α 1 3 sin γ 2
and also for the phases of the motion
γ ˙ 1 = ε σ 2 + c 1,1 2 8 ω ¯ 1 3 g 1,1 8 ω ¯ 1 α 1 2 3 g 1,3 4 ω ¯ 1 α 2 2 3 g 1,2 8 ω ¯ 1 α 1 α 2 cos γ 2 + f 1 2 ω ¯ 1 sin γ 1
γ ˙ 2 = 3 c 1,1 2 8 ω ¯ 1 c 2,2 2 8 ω ¯ 2 ε σ 1 9 g 1,1 8 ω ¯ 1 3 g 2,2 4 ω ¯ 1 α 1 2 9 g 1,3 4 ω ¯ 1 9 g 2,4 4 ω ¯ 1 α 2 2 + 3 f 1 2 ω ¯ 1 α 1 sin γ 1 9 g 1,2 α 1 α 2 8 ω ¯ 1 g 2,1 α 1 3 8 ω ¯ 1 α 2 cos γ 2
The steady solutions for these responses are determined through vanishing the right sides of Equations (28)–(31), which correspond to the fixed-point states of amplitudes of the first and second flapwise modes as well as their phases. The stability of the steady motion can be evaluated in terms of the eigenvalues of the following Jacobian matrix of the fixed points
J = α ˙ 1 α 1 α ˙ 1 α 2 α ˙ 1 γ ˙ 1 α ˙ 1 γ ˙ 2 α ˙ 2 α 1 α ˙ 2 α 2 α ˙ 2 γ ˙ 1 α ˙ 2 γ ˙ 2 γ ˙ 1 α 1 γ ˙ 1 α 2 γ ˙ 1 γ ˙ 1 γ ˙ 1 γ ˙ 2 γ ˙ 2 α 1 γ ˙ 2 α 2 γ ˙ 2 γ ˙ 1 γ ˙ 2 γ ˙ 2 ,
where the components of J are provided in Equations (A16)–(A19), Appendix A. The response is said to be asymptotically stable if all of the eigenvalues possess negative real parts.
Afterwards, the approximations of  q ¯ 1  and  q ¯ 2  can be obtained as
q ¯ 1 = α 1 cos l Ω t γ 1 + O ε 3
q ¯ 2 = α 2 cos 3 l Ω t 3 γ 1 γ 2 + O ε 3
When the combination resonance happens, the steady responses are determined through conditions
[ ( ε σ 1 + 3 ε σ 2 + c 2,2 2 8 ω ¯ 2 ) 6 g 1,2 ω ¯ 1 α 2 2 ( ε σ 1 + c 1,1 2 8 ω ¯ 1 ) 2 g 2,1 ω ¯ 2 α 1 2 + ( 3 g 1,3 g 2,1 2 ω ¯ 1 ω ¯ 2 9 g 1,2 g 2,2 2 ω ¯ 1 ω ¯ 2 ) α 1 2 α 2 2 + 3 g 1,1 g 2,1 4 ω ¯ 1 ω ¯ 2 α 1 4 9 g 1,2 g 2,4 4 ω ¯ 1 ω ¯ 2 α 2 4 2 ] = ( f 1 g 2,1 α 1 ω ¯ 1 ω ¯ 2 ) 2 ( g 2,1 c 1,1 ω ¯ 2 α 1 2 + 3 g 1,2 c 2,2 ω ¯ 1 α 2 2 ) 2
ε σ 1 + 3 ε σ 2 + c 2,2 2 8 ω ¯ 2 2 g 2,2 α 1 2 4 ω ¯ 2 2 g 2,4 α 2 2 8 ω ¯ 2 = g 2,1 α 1 3 8 ω ¯ 2 α 2 2 c 2,2 2 4
Nonetheless,  α 1 , 2  are constrained by the inequalities
f 1 g 2,1 α 1 ω ¯ 1 ω ¯ 2 2 f 1 g 2,1 c 1,1 ω ¯ 2 α 1 2 + 3 g 1,2 c 2,2 ω ¯ 1 α 2 2 2 0 ,   g 2,1 α 1 3 4 ω ¯ 2 α 2 2 c 2,2 2 4 0
in order to assure the absolutes of  sin γ 1,2  and  cos γ 1,2  in Equations (28)–(31) are not greater than 1.
In the solution procedure, a simulation is implemented using the Runge–Kutta method with variable step size to obtain numerical results of Equation (18) and the analytical results of the transient system (28)–(31), whereas the analytical steady Equations (35)–(36) are solved through Newton’s iteration method. Nevertheless, it is necessary to solve the amplitude of the steady response (i.e., B in Equation (2)) via iterations before the solution of the governing equations begins. Such an iterative process for B is described in Appendix B.

4. Nonlinear Dynamics and Combination Resonance of Shrouded Blade

A flexible shrouded blade is used for examples and discussions. The major parameters involving the blade design and operation are provided in Table 1. The natural frequencies are  ω ¯ 1 = 0.7580  and  ω ¯ 2 = 2.9238  for the first-order and second-order modes in the flapwise direction. However, these frequencies become  ω ¯ 1 = 0.9993  and  ω ¯ 2 = 2.9954  when one considers the action of contact and friction on the shroud. The condition for concurrence of an internal resonance of type 1:3 between the two flapwise modes and the primary resonance of the first flapwise mode is expected when considering the aforementioned contact and friction.
To validate the proposed method, the steady responses of the original Equation (18) are solved with the analytical, multiple-scale results  q ¯ k  through Equations (33)–(34), which are obtained by computing Equations (28)–(31). As shown in Figure 5, these analytical solutions are found in good agreement with their numerical counterparts. Their discrepancies are inevitable, since higher order terms have been excluded from multiple-scale-based expansions, and they appear more noticeable for the second flapwise mode than the first. In summary, the proposed modeling and approximate approach are applicable for steady-response solutions of the shrouded blade in the situation of combination resonance.

4.1. Steady Responses of Flexible Shrouded Blade

For flexible blades, the tilting angle and contact stiffness of the shroud are key to the contact status of adjacent shrouds. Further, the blade rotating induces centrifugal load and the subsequent spinning softening that reduces the axial stiffness of the blade. In the following, the original governing Equation (18) is solved with various tilting angle, contact stiffness and rotation speed values to reflect their influence on the response.

4.1.1. Effect of Tilting Shroud

The bifurcation diagram of the blade vibration is demonstrated using the tilting angle of shroud varying between 0 and 90°, Figure 6. Based on the characteristic of the steady response of the first-order and second-order flapwise modes, the range of the tilting angle is divided into three stages, denoted by stage-I, -II and -III, respectively, with  α 1   = 38.96° and  α 2   = 83.65° being the dividing points that separate the stages.
The steady responses of the first-order mode  q ¯ 1  under variable  α  are shown in Figure 6a. In stage-I,  q ¯ 1  is quasi-periodic motion. Then,  q ¯ 1  becomes periodic-1, and the resonance peak of the primary resonance,  ω ¯ 1 = l , can be observed in stage-II. In addition, the jumping of the amplitude and the phenomenon of so-called “hardening” can be seen for stage-II, owing to the geometric nonlinearity originating from Equation (13). Finally, for stage-III,  q ¯ 1  is converted back to quasi-periodic motion. It is worthwhile to mention the curving of the resonance summit towards the side where the ratio between the frequency of external excitation and natural frequency,  l / ω ¯ 1 , goes up is defined as the hardening spring. In these cases, the jumping of the amplitude and the unstable motion appear in the steady response. However, the steady response is obtained through the Runge–Kutta method in Section 4.1, by which only the stable motion can be observed.
Corresponding to the first flapwise mode, the periodicity of  q ¯ 2  is the same as that of  q ¯ 1  in stage-I and -III, but  q ¯ 2  represents 3-periodic motion in stage-II, Figure 6b. In addition, the jumping of amplitude and the hardening phenomenon caused by the primary resonance can also be observed in stage-II. This is because the primary and internal resonances concur as  ω ¯ 1 = l  and  ω ¯ 2 = 3 ω ¯ 1 , which leads to the desired combination resonance. Therefore, the primary resonance is also observed in  q ¯ 2  due to the exchanges in energy between the first-order and second-order modes in the flapwise direction.
Further examinations are made focusing on the blade vibration carrying four different tilts angle of shrouds: 4.07° (in stage-I), 40.28°, 40.28° (in stage-II) and 85.48° (in stage-III), and the results are provided through Figure 7, Figure 8, Figure 9 and Figure 10. As found in Figure 7, both  q ¯ 1  and  q ¯ 2  are quasi-periodic motions with spectra containing the aerodynamic frequency  f Q , the first flapwise frequency  f q 1  and the combined frequencies of  f q 1 + 2 f Q 2 f q 1 + f Q  and  3 f Q . Further, the amplitude of  q ¯ 1  is synchronous and π is out-of-phase against the amplitude of  q ¯ 2 . This implies that it is likely to have internal resonance when  α  is small, since the curving of the resonance summit occurs as the blade appears hardened in this scenario. The energy carried by the first order flapwise mode is able to transmit to the second flapwise mode and vice versa.
The steady responses of the blade with  α =  40.28° and 44.81° are shown in Figure 8 and Figure 9, respectively.  q ¯ 1  and  q ¯ 2  are identified as periodic-1 motion and periodic-3 motion, respectively, and both  q ¯ 1  and  q ¯ 2  are governed by  f Q  and its odd-fold multiples (i.e.,  3 f Q ,   5 f Q  and so forth). The exchanges in energy between the modes can be observed clearly in Figure 8 but not in Figure 9. This is explained by the fact that the primary resonance is subjected to a strong enhancement on the steady response, which happens when  α =  44.81°. Although the amplitude increases under primary resonance, the change in motion caused solely by internal resonance is minor and thus can hardly be observed in Figure 9.
The steady response of the blade when  α =  85.48° is shown in Figure 10. Both  q ¯ 1  and  q ¯ 2  are found to be quasi-periodic with dominant spectra of  f Q , the first flapwise frequency  f q 1  and the combined frequencies of  f q 1 + 2 f Q 2 f q 1 + f Q  and  3 f Q . With the current section of α, the aforementioned exchanges in energy between the aforementioned modes cannot be observed. In fact, the frequency difference between these modes grows with the increasing α (c.f. Figure 6b), which weakens the effect of internal resonance on these modes.

4.1.2. Effect of Contact Stiffness

In the following, the vibration response of the blade is studied using different stiffness values of normal contact. The bifurcation of the blade vibration is shown in Figure 11, where the contact stiffness ranges from 1 MN/m to 10 MN/m. Similar to Section 4.1.1, the range of  k s  is divided into three stages, where  k s 1 = 1.72  and  k s 2 = 6.97  MN/m are the dividing points of these stages.
In Figure 11a, the steady  q ¯ 1  is identified as a periodic-1 motion in the entire range of  k s . Again, jumping of amplitude and hardening spring due to the primary resonance can be observed in stage-II. As for the steady second flapwise mode  q ¯ 2 , it is periodic-1 in stages-I and III but periodic-3 in stage-II, Figure 11b. Moreover, the energy exchanges caused by internal resonance occur in stage-II; thereby, the jumping of the amplitude and the hardening spring due to the primary resonance can also be observed.
Further explorations are carried out using three different contact stiffness values, i.e., 1.407 MN/m (in stage-I), 2.764 MN/m (in stage-II) and 8.191 MN/m (in stage-III), and the response is provided through Figure 12, Figure 13 and Figure 14. As depicted in Figure 12, both  q ¯ 1  and  q ¯ 2  are periodic-1 motions and are spectrally governed by  f Q  and its odd multiples ( 3 f Q  and  5 f Q , etc.). Similar to previous cases, the amplitudes of  q ¯ 1  and  q ¯ 2  are both periodic, while the amplitude of  q ¯ 1  is synchronous and π is out-of-phase against the amplitude of  q ¯ 2 . Comparing the results presented in the three figures, it is found that the first flapwise modal response is always period-1. Nonetheless, the response of the second flapwise mode changes from a period-1 to period-3 motion and then returns to period-1 when the stiffness of normal contact increases.

4.1.3. The Effect of Rotation Speed

In the following, the contribution of the rotation speed to the blade vibration is studied. The bifurcation of the modal responses is shown in Figure 15, where the range of  Ω  is from 2000 r/min to 12,000 r/min. This range of speed can also be divided into three stages:  2.0 × 1 0 3 ~ Ω 1  (stage-I),  Ω 1 ~ Ω 2  (stage-II) and  Ω 2 ~ 1.2 × 1 0 4  r/min (stage-III), respectively, where  Ω 1 =  4.46 × 103 r/min and  Ω 2 =  8.08 × 103 r/min.
As shown in Figure 15a, the steady responses of  q ¯ 1  remains in a periodic-3 motion for the entire range of  Ω . Similar to the case of contact stiffness,  q ¯ 2  remains periodic-1 in stages-I and III, while  q ¯ 2  is periodic-3 in stage-II, Figure 15b. Moreover, the jumping of the amplitude and the hardening spring can also be observed in stage-II due to the combination resonance.
Next, three values of rotational speed, i.e., 2000 r/min (in stage-I), 7480 r/min (in stage-II) and 11,000 r/min (in stage-III), are chosen to examine the modal responses, Figure 16, Figure 17 and Figure 18. It is found that the steady first flapwise response is always period-1. However, the second flapwise modal response turns from the period-1 to period-3 motion and then goes back to period-1 with the increasing rotation speed. In addition, exchanges in energy can also be observed in stage-III, a sign of the interaction between modes of the first-order and second-order flapwise motion due to internal resonance.

4.2. Parametric Analysis in the Case of Combination Resonance

Based on Section 4.1, both primary and internal resonances can be observed in the steady responses of the first-order and second-order modes in the flapwise motion under variable tilting angle of the shroud, stiffness of normal contact and rotation speed. Also, combination resonance can be identified with the parameters in Table 1. In this section, the parametric influences of detuning and shroud design ( δ s   and   k s ) are highlighted on the blade response in the case of combination resonance. It is worth noting that the forcing and damping parameters (i.e.,  f ¯ k  and  c ¯ k ) as well as the geometric nonlinearity coefficients  g ¯ k , j  are all related to  δ s  and  k s . For convenience, a coefficient of damping is defined for the flapwise modes of interest:
ζ k = c ¯ k , k / 2 ω ¯ k , k
Based on Figure 19, for a fixed  δ s , the flapwise damping ratios  ζ 1 , 2  rise as  k s  increases. Both  ζ 1,2  decrease with the increasing  δ s  when  k s  is constant. Further, the effect of  δ s  on the modal dampings of the flapwise motion is noticeably weak compared with that of  k s . In addition, the coefficient of the nondimensional equations in Appendix A shows that  g ¯ k , j  ascends as  δ s  increases in contrast to  f ¯ k . Both  g ¯ k , j  and  f ¯ k  are insensitive to the contact stiffness.

4.2.1. Effect of  ε σ 1  with Various Contact Stiffness and Gap Dimension

The amplitude–frequency relations of both  α 1  and  α 2  are presented in Figure 20 using different  ε σ 1  and constant  k s   = 3.55 × 106 N/m and  δ s   = 2 × 10−5 m for the case of combination resonance. It is observed that the blade resembles a hardened spring in vibration, which can be attributed to the presence of the  g ¯ k , j  terms. Based on Figure 20a, both stable and unstable motions can be identified within different bands of frequency, along with four turning points, i.e., A1, B1, C1 and D1. The real and imaginary parts of eigenvalues are presented in Figure 21, corresponding to the steady responses depicted in Figure 20.
As far as the fixed points are concerned, the branch left to summit A1 ( ε σ 1 = 1.52 ) are stable focuses of  α 1  before the occurrence of primary resonance. The connecting route from A1 to B1 ( ε σ 1 = 0.88 ) involves identified saddle-typed fixed points that destabilize  α 1 . The branch between points B1 and C1 ( ε σ 1 = 2.00 ) has stable focuses at  α 1 , as opposed to the saddle points represented by the branch from point C1 to D1 ( ε σ 1 = 0.72 ). Moreover, the branch right of the turning point D1 returns to a stable focus. As a matter of fact,  α 1  jumps downward from A1 and C1 once  ε σ 1  increases to critical values. The upward jumping happens at D1 as  ε σ 1  decreases to a certain critical value. However, with the decreasing  ε σ 1 α 1  will jump at the turning point B1. The responses either jumps up and flows following the left branch of point A1 or goes down to the right branch of point D1.
As for  α 2  in the case of combination resonance, a similar phenomenon of hardening can be observed in Figure 20b through three stable regions, two unstable regions and four turning points. Comparing the amplitudes in Figure 20a,b, the branches of response between points B2 ( α 2 = 0.206 ) and C2 ( α 2 = 0.209 ) are greater than that at A2 ( α 2 = 0.124 ), which makes a difference in the direction of jumping at points A2 ( ε σ 1 = 1.52 ) and B2 ( ε σ 1 = 0.88 ); however,  α 1  in the branch B1 ( α 1 = 0.516 ) and C1 ( α 1 = 0.853 ) are less than that at A1 ( α 1 = 1.20 ). This confirms the existence of energy exchange between the first and second modes in the flapwise direction.
To demonstrate the effect of contact stiffness,  k s  is reduced by half (denoted subsequently by  0.5 k s ) and doubled (by  2 k s ), respectively. The solutions of  α 1  and  α 2  are presented in Figure 22 and Figure 23, respectively. Unlike the results shown in Figure 20, the decreasing  k s  widens the range of frequency corresponding to instability of the vibration and delays the onset frequency for primary resonance. Specifically, the frequency interval for unstable motion expands from (0.88, 1.52) and (0.72, 2.00) to (0.84, 7.05) and (0.83, 10.01), respectively. The larger  k s  is able to remove the instability of  α 1  and  α 2  under the combination resonance alongside the turning points. This is because the damping of the blade rises with the increasing  k s  (see Figure 19) and thus refrains the vibration of the blade. When  k s  is large enough, the hardening spring-like reaction of the blade at the resonance point can been completely suppressed.
To demonstrate the effect of the shroud gap, the original  δ s  is halved (denoted by  0.5 δ s ) and doubled (by  2 δ s ), respectively. The solutions of  α 1  and  α 2  are presented in Figure 24 and Figure 25, accordingly. Similar to the previous finding, the blade vibrates in higher amplitude as  δ s  becomes larger, and the onset of primary resonance postpones. In this case, the regions of  ε σ 1  corresponding to instability expand from (0.88, 1.52) and (0.72, 2.00) to (1.14, 3.05) and (0.93, 3.97) and narrow to (0.60, 0.76) and (0.60, 0.76) in Figure 25 when the gap doubles its original value.
Nevertheless, the effect of  δ s  is noticeably weaker than that of  k s . The rising  δ s  influences the level of vibration through changing the aerodynamic force ( f k ), damping ( ζ k ) and geometric nonlinear terms ( g ¯ k , j ). It is worth mentioning that changing  g ¯ k , j  has negligible contribution to the level of vibration [32,33]. An increasing  δ s  has a double effect on the response: it does reduce the blade damping (see Figure 19) and hence increases the amplitudes of  α 1  and  α 2 ; on the other hand, it prevents the response by decreasing the aerodynamic force. Comparing Figure 24 with Figure 25, it is found that the level of vibration decreases with the increasing  δ s ; thereby, the strengthening of increasing  δ s  on the level of vibration overwhelms its weakening.

4.2.2. Effect of  ε σ 2  with Various Contact Stiffness and Gap Dimension

In this subsection, the contribution of parameter  ε σ 2  to the steady response of the blade is investigated. Given  k s  = 3.55 × 106 N/m and  δ s  = 2.0 × 10−5 m, the amplitudes of the first- and second-order modal response of the flapwise motion are plotted in Figure 26. Unlike the previous discussion, the phenomenon of softening spring are discovered due to the geometric nonlinearity associated with  g ¯ k , j .
There are two regions of stable motion and one region of unstable motion along with two turning points, i.e., A1 and B1, as depicted in Figure 26a. The stability of the steady response is demonstrated through eigenvalue analysis at the fixed points of  α 1  and  α 2  in Figure 27. In viewing the fixed points, the branch left of point A1 ( ε σ 2 = 1.504 ) is found to be a stable focus of  α 1 . The connection between A1 and the lowest point B1 ( ε σ 2 = 1.408 ) corresponds to the saddle points of  α 1 , demonstrating an unstable motion of the blade. Indeed,  α 1  jumps downward from point A1 once  ε σ 2  increases to be critical, while an upward jumping is located at point B1. As shown in Figure 26b, the hardening spring phenomenon can be observed alongside two stable regions and one unstable region. With the increasing  ε σ 2 α 2  rises in the branch left of the turning point A2 and drops to the region left of the other turning point B2, which is opposite to the trend of  α 1  in Figure 26a. This suggests that  ε σ 2  is the key influence on the internal resonance.
The steady responses of the blade vibration versus  k s  are presented in Figure 28 and Figure 29 for the combination resonance. As shown in Figure 28, both the level of vibration and the interval of frequency for unstable motion are enlarged. There exists an overlap in the unstable region and the right branch of the turning point B1. The unstable region expands from  ε σ 2  (1.408, 1.504) to (1.423, 1.679). Similar to the solutions in Figure 23, the unstable region as well as the turning points can be removed by increasing  k s  for the same reason presented in Section 4.2.1.
Finally, the steady responses in the situation of combination resonance are presented in Figure 30 and Figure 31 with various shroud gap  δ s . It is found that a decreasing  δ s  escalates the vibration level and broadens the interval of instability of motion in the frequency domain. The unstable region now expands from (1.408, 1.504) to (1.783, 1.986) in the frequency domain at a halved  δ s  and narrows to (1.097, 1.116) in Figure 31. The effect of  k s  on the responses is significantly weaker than that of  k s  on  α 1,2 . In addition, the overlap in the unstable region and the right branch of the turning point B1 can also be observed in Figure 30a.

5. Conclusions

In this paper, the nonlinear dynamics of a flexible turbine blade are investigated considering the contact and friction of shrouds, focusing on the case of combination resonance of flapwise motion induced by a 1:3 internal resonance between the first- and second-order modes and a primary resonance of the first-order mode. The stiffness and damping properties are expressed by linearizing the normal and tangential forces of shrouds using the method of harmonic balance. The steady responses of the blade in the situation of combination resonance are obtained with the multiple-scale method. The contributions of detuning and shroud parameters to the steady responses of the blade in the case combination resonance are analyzed. Several conclusions are drawn as follows:
(1)
The combination resonance involving the flapwise modes can be triggered by co-existing internal and primary resonances. Nevertheless, such resonance will not likely occur for separated blades where contact and friction on their shroud interfaces are removed.
(2)
Motion of the first-order flapwise mode can change from quasi-periodic to period-1 and then back to quasi-periodic with an increasing tilting angle of shroud. As for the second-order mode, it can change from quasi-periodic to period-3 and finally quasi-periodic.
(3)
With increasing contact stiffness and rotation speed, the first-order flapwise motion is period-1, while the second modal response can change between period-1 and period-1.
(4)
For primary resonance, the hardening spring behavior of the blade is attributed to geometric nonlinearity in the blade motion. On the other hand, such a nonlinearity can lead to a softening spring behavior of the blade in the scenario of internal resonance.
(5)
Less shroud contact stiffness and gap dimension can result in a stronger vibration response, wider region of unstable motion in the frequency domain and delayed onset of combination resonance.

Author Contributions

Conceptualization, H.L. and Y.W.; methodology, G.Y. and Y.W.; validation, H.L. and Y.W.; formal analysis, G.Y.; investigation, H.L. and Y.W.; resources, D.M. and Z.Y.; writing—original draft preparation, H.L. and G.Y.; writing—review and editing, Y.W.; funding acquisition, Y.W. and D.M. All authors have read and agreed to the published version of the manuscript.

Funding

The International Cooperation Fund Project of DBJI, Dalian University of Technology (ICR2109); the Free Exploration Project of the State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment; and the National Science and Technology Major Project (2019-IV-0019-0087).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available to protect the information of key design parameters.

Acknowledgments

The authors are grateful to DBJI (the Joint Institute of Dalian University of Technology and Belarus State University) and the State Key Laboratory of Structural Analysis, Optimization and CAE Software for Industrial Equipment for their sponsorships.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. The authors declare that there are not any potential commercial interests.

Appendix A

The coefficients  c k ω k g k , j  and  f k  of Equation (16) are defined as
c k = 1 m k 0 L c v γ k 2 v d x c n sin α + c f cos α γ k 2 v L
ω k 2 = 1 m k 0 L E I γ k 2 v d x + 0 L x L m l Ω 2 e + y d y γ k 2 v d x 0 L Ω 2 m l γ k 2 v d x i Ω 2 m s + k n sin α + k f cos α γ k 2 v L
g k , 1 = 1 m k 0 L 1 2 E A γ k v γ 1 3 v d x
g k , 2 = 1 m k 0 L 1 2 E A γ k v γ 1 2 v γ 2 v d x
g k , 3 = 1 2 m k 0 L E A γ k v γ 1 v γ 2 2 v d x
g k , 4 = 1 2 m k 0 L E A γ k v γ 2 3 v d x
f k = 1 m k Q 0 γ k v L ,
where
m k = 0 L m l γ k 2 v d x + m s γ k 2 v
The coefficients  c ¯ k ω ¯ k g ¯ k , j  and  f ¯ k  of Equation (18) are defined as
c ¯ k = 1 Ω m k 0 L c v γ k 2 v d x c n sin α + c f cos α γ k 2 v L
ω ¯ k 2 = 1 Ω 2 m k 0 L E I γ k 2 v d x + 0 L x L m l Ω 2 e + y d y γ k 2 v d x 0 L Ω 2 m l γ k 2 v d x i Ω 2 m s + k n sin α + k f cos α γ k 2 v L
g ¯ k , 1 = δ s 2 Ω 2 m k 0 L 1 2 E A γ k v γ 1 3 v d x
g ¯ k , 2 = δ s 2 2 Ω 2 m k 0 L E A γ k v γ 1 2 v γ 2 v d x
g ¯ k , 3 = δ s 2 2 Ω 2 m k 0 L E A γ k v γ 1 v γ 2 2 v d x
g ¯ k , 4 = δ s 2 2 Ω 2 m k 0 L E A γ k v γ 2 3 v d x
f ¯ k = 1 δ s Ω 2 m k Q 0 γ k v L
The components of J in Equation (32) are defined as
J 1,1 = c 1,1 2 3 g 1,2 4 ω 1 α 1 α 2 sin γ 2 ,   J 1,2 = 3 g 1,2 8 ω 1 α 1 2 sin γ 2 ,   J 1,3 = f 1 2 ω 1 sin γ 1 ,   J 1,4 = 3 g 1,2 8 ω 1 α 1 2 α 2 cos γ 2 ,
J 2,1 = 3 g 2,1 8 ω 2 α 1 2 sin γ 2 ,   J 2,2 = c 2,2 2 ,   J 2,3 = 0 ,   J 2,4 = g 2,1 8 ω 2 α 1 3 cos γ 2
J 3,1 = 3 g 1,1 α 1 4 ω 1 3 g 1,2 α 2 8 ω 1 cos γ 2 f 1 2 ω 1 α 1 2 sin γ 1 , J 3,2 = 3 g 1,3 α 2 2 ω 1 3 g 1,2 α 1 8 ω 1 cos γ 2 ,   J 3,3 = f 1 2 ω 1 α 1 cos γ 1 ,   J 3,4 = 3 g 1,2 α 1 α 2 8 ω 1 sin γ 2
J 4,1 = 9 g 1,1 4 ω 1 3 g 2,2 2 ω 2 α 1 9 g 1,2 α 2 8 ω 1 3 g 2,1 α 1 2 8 ω 1 α 2 cos γ 2 3 f 1 2 ω 1 α 1 2 sin γ 1 , J 4,2 = 9 g 1,3 2 ω 1 3 g 2,4 4 ω 2 α 2 9 g 1,2 α 1 8 ω 1 + g 2,1 α 1 3 8 ω 1 α 2 2 cos γ 2 ,   J 4,3 = 3 f 1 2 ω 1 α 1 cos γ 1 ,   J 4,4 = 9 g 1,2 α 1 α 2 8 ω 1 g 2,1 α 1 3 8 ω 1 α 2 sin γ 2

Appendix B

An iterative can be applied to solve the amplitude of the steady vibration responses (i.e., B in Equation (2)) before the governing equations in the paper are solved. Let the initial value of B be denoted as B0. The properties of stiffness and damping can be determined by substituting B0 into Equations (8), (9), (11) and (12), which represented as  k n , 0 k f , 0 c n , 0  and  c f , 0 , respectively. Then, the steady responses of Equation (18) along with the transient systems (28)–(31) are solved using  k n , 0 k f , 0 c n , 0  and  c f , 0 , and the amplitude of the steady response is denoted by B1. The above process is repeated with B1 as the initial value until the error between Bi and Bi−1 is below the preassigned tolerance. At last, Bi is accepted as the amplitude of the desired steady responses.

References

  1. Cigeroglu, E. Development of Microslip Friction Models and Forced Response Prediction Methods for Frictionally Constrained Turbine Blades. Ph.D. Thesis, The Ohio State University, Columbus, OH, USA, 2007. [Google Scholar]
  2. Xie, F.; Ma, H.; Cui, C.; Wen, B. Vibration response comparison of twisted shrouded blades using different impact models. J. Sound Vib. 2017, 397, 171–191. [Google Scholar] [CrossRef]
  3. Ju, D.; Sun, Q. Modeling of a Wind Turbine Rotor Blade System. J. Vib. Acoust. 2017, 139, 051013. [Google Scholar] [CrossRef]
  4. Chiu, Y.J.; Chen, D.Z.; Yang, C.H. Influence on Coupling Vibration of Rotor System with Grouped Blades due to Mistuned Lacing Wire. Appl. Mech. Mater. 2011, 101–102, 1119–1125. [Google Scholar] [CrossRef]
  5. Ma, H.; Xie, F.; Nai, H.; Wen, B. Vibration characteristics analysis of rotating shrouded blades with impacts. J. Sound Vib. 2016, 378, 92–108. [Google Scholar] [CrossRef]
  6. Chatterjee, A.; Kotambkar, M.S. Modal characteristics of turbine blade packets under lacing wire damage induced mistuning. J. Sound Vib. 2015, 343, 49–70. [Google Scholar] [CrossRef]
  7. Shadmani, M.; Tikani, R.; Ziaei-Rad, S. On using a distributed-parameter model for modal analysis of a mistuned bladed disk rotor and extracting the statistical properties of its in-plane natural frequencies. J. Sound Vib. 2019, 438, 324–343. [Google Scholar] [CrossRef]
  8. Zhou, X.; Huang, K.; Li, Z. Geometrically nonlinear beam analysis of composite wind turbine blades based on quadrature element method. Int. J. Non-Linear Mech. 2018, 104, 87–99. [Google Scholar] [CrossRef]
  9. Karimi, A.H.; Shadmani, M. Nonlinear vibration analysis of a beam subjected to a random axial force. Arch. Appl. Mech. 2019, 89, 385–402. [Google Scholar] [CrossRef]
  10. Griffin, J.H. Friction Damping of Resonant Stresses in Gas Turbine Engine Airfoils. J. Eng. Power 1980, 102, 329–333. [Google Scholar] [CrossRef]
  11. Menq, C.-H.; Bielak, J.; Griffin, J. The influence of microslip on vibratory response, part I: A new microslip model. J. Sound Vib. 1986, 107, 279–293. [Google Scholar] [CrossRef]
  12. Iwan, W.D. On a Class of Models for the Yielding Behavior of Continuous and Composite Systems. J. Appl. Mech. 1967, 34, 612–617. [Google Scholar] [CrossRef]
  13. Yang, B.; Menq, C. Characterization of 3D contact kinematics and prediction of resonant response of structures having 3D frictional constraint. J. Sound Vib. 1998, 217, 909–925. [Google Scholar] [CrossRef]
  14. Nan, G.; Lou, J.; Song, C.; Tang, M. A New Approach for Solving Rub-Impact Dynamic Characteristics of Shrouded Blades Based on Macroslip Friction Model. Shock. Vib. 2020, 2020, 8147143. [Google Scholar] [CrossRef]
  15. Brach, R.M. Mechanical Impact Dynamics: Rigid Body Collisions; John Wiley & Sons: New York, NY, USA, 1991. [Google Scholar]
  16. He, B.; Ouyang, H.; He, S.; Ren, X.; Mei, Y. Dynamic analysis of integrally shrouded group blades with rubbing and impact. Nonlinear Dyn. 2018, 92, 2159–2175. [Google Scholar] [CrossRef]
  17. He, S.; Si, K.; He, B.; Yang, Z.; Wang, Y. Rub-Impact Dynamics of Shrouded Blades under Bending-Torsion Coupling Vibration. Symmetry 2021, 13, 1073. [Google Scholar] [CrossRef]
  18. Mashayekhi, F.; Nobari, A.; Zucca, S. Hybrid reduction of mistuned bladed disks for nonlinear forced response analysis with dry friction. Int. J. Non-Linear Mech. 2019, 116, 73–84. [Google Scholar] [CrossRef]
  19. She, H.; Li, C.; Tang, Q.; Wen, B. Veering and merging analysis of nonlinear resonance frequencies of an assembly bladed disk system. J. Sound Vib. 2021, 493, 115818. [Google Scholar] [CrossRef]
  20. Wei, S.-T.; Pierre, C. Localization Phenomena in Mistuned Assemblies with Cyclic Symmetry Part II: Forced Vibrations. J. Vib. Acoust. 1988, 110, 439–449. [Google Scholar] [CrossRef]
  21. Fang, X.; Tang, J.; Jordan, E.; Murphy, K. Crack induced vibration localization in simplified bladed-disk structures. J. Sound Vib. 2006, 291, 395–418. [Google Scholar] [CrossRef]
  22. Picou, A.; Capiez-Lernout, E.; Soize, C.; Mbaye, M. Robust dynamic analysis of detuned-mistuned rotating bladed disks with geometric nonlinearities. Comput. Mech. 2020, 65, 711–730. [Google Scholar] [CrossRef]
  23. Zhao, W.; Zhang, D.; Xie, Y. Vibration analysis of mistuned damped blades with nonlinear friction and contact. J. Low Freq. Noise Vib. Act. Control. 2019, 38, 1505–1521. [Google Scholar] [CrossRef]
  24. Larsen, J.W.; Nielsen, S.R.K. Nonlinear parametric instability of wind turbine wings. J. Sound Vib. 2007, 299, 64–82. [Google Scholar] [CrossRef]
  25. Karimi, B.; Moradi, H. Nonlinear kinematics analysis and internal resonance of wind turbine blade with coupled flapwise and edgewise vibration modes. J. Sound Vib. 2018, 435, 390–408. [Google Scholar] [CrossRef]
  26. Zhang, W.; Liu, G.; Siriguleng, B. Saturation phenomena and nonlinear resonances of rotating pretwisted laminated composite blade under subsonic air flow excitation. J. Sound Vib. 2020, 478, 115353. [Google Scholar] [CrossRef]
  27. Chu, S.; Cao, D.; Sun, S.; Pan, J.; Wang, L. Impact vibration characteristics of a shrouded blade with asymmetric gaps under wake flow excitations. Nonlinear Dyn. 2013, 72, 539–554. [Google Scholar] [CrossRef]
  28. Allara, M. A model for the characterization of friction contacts in turbine blades. J. Sound Vib. 2009, 320, 527–544. [Google Scholar] [CrossRef]
  29. Pai, P.; Rommel, B.; Schulz, M.J. Non-linear vibration absorbers using higher order internal resonances. J. Sound Vib. 2000, 234, 799–817. [Google Scholar] [CrossRef]
  30. Sayed, M.; Kamel, M. 1:2 and 1:3 internal resonance active absorber for non-linear vibrating system. Appl. Math. Model. 2012, 36, 310–332. [Google Scholar] [CrossRef]
  31. Eftekhari, M.; Ziaei-Rad, S.; Mahzoon, M. Vibration suppression of a symmetrically cantilever composite beam using internal resonance under chordwise base excitation. Int. J. Non-Linear Mech. 2013, 48, 86–100. [Google Scholar] [CrossRef]
  32. Li, L.; Li, Y.; Liu, Q.; Lv, H. Flapwise non-linear dynamics of wind turbine blades with both external and internal resonances. Int. J. Non-linear Mech. 2014, 61, 1–14. [Google Scholar] [CrossRef]
  33. Yuan, G.; Wang, Y. Internal, primary and combination resonances of a wind turbine blade with coupled flapwise and edgewise motions. J. Sound Vib. 2021, 514, 116439. [Google Scholar] [CrossRef]
  34. She, H.; Li, C. Effects of centrifugal stiffening and spin softening on nonlinear modal characteristics of cyclic blades with impact–friction coupling. Nonlinear Dyn. 2022, 110, 3229–3254. [Google Scholar] [CrossRef]
  35. Nayfeh, A.H.; Mook, D.T. Nonlinear Oscillations; John Wiley & Sons: New York, NY, USA, 1995. [Google Scholar]
Figure 1. Mechanical model of shrouded rotating disk.
Figure 1. Mechanical model of shrouded rotating disk.
Machines 12 00059 g001
Figure 2. Sketch of a flexible shrouded blade.
Figure 2. Sketch of a flexible shrouded blade.
Machines 12 00059 g002
Figure 3. The sketch of aerodynamic, normal and friction forces of the shroud.
Figure 3. The sketch of aerodynamic, normal and friction forces of the shroud.
Machines 12 00059 g003
Figure 4. Hysteretic loop of friction versus tangential relative displacement.
Figure 4. Hysteretic loop of friction versus tangential relative displacement.
Machines 12 00059 g004
Figure 5. Analytical and numerical steady responses. (a) First flapwise mode and (b) second flapwise mode.
Figure 5. Analytical and numerical steady responses. (a) First flapwise mode and (b) second flapwise mode.
Machines 12 00059 g005
Figure 6. Bifurcation of the blade under variable  α . (a q ¯ 1  and (b q ¯ 2 .
Figure 6. Bifurcation of the blade under variable  α . (a q ¯ 1  and (b q ¯ 2 .
Machines 12 00059 g006
Figure 7. Steady responses of the blade with α = 4.07°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 7. Steady responses of the blade with α = 4.07°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g007aMachines 12 00059 g007b
Figure 8. Steady responses of the blade with α = 40.28°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 8. Steady responses of the blade with α = 40.28°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g008aMachines 12 00059 g008b
Figure 9. Steady responses of the blade with α = 44.81°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 9. Steady responses of the blade with α = 44.81°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g009aMachines 12 00059 g009b
Figure 10. Steady responses of the blade with  α =  85.48°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 10. Steady responses of the blade with  α =  85.48°. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g010aMachines 12 00059 g010b
Figure 11. Bifurcation of the blade with variable  k s . (a q ¯ 1  and (b q ¯ 2 .
Figure 11. Bifurcation of the blade with variable  k s . (a q ¯ 1  and (b q ¯ 2 .
Machines 12 00059 g011
Figure 12. Steady responses of the blade with ks = 1.407 MN/m. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 12. Steady responses of the blade with ks = 1.407 MN/m. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g012
Figure 13. Steady responses of the blade with ks = 2.764 MN/m. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 13. Steady responses of the blade with ks = 2.764 MN/m. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g013aMachines 12 00059 g013b
Figure 14. Steady responses of the blade with ks = 8.191 MN/m. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 14. Steady responses of the blade with ks = 8.191 MN/m. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g014aMachines 12 00059 g014b
Figure 15. Bifurcation diagram of the blade under variable W. (a q ¯ 1  and (b q ¯ 2 . .
Figure 15. Bifurcation diagram of the blade under variable W. (a q ¯ 1  and (b q ¯ 2 . .
Machines 12 00059 g015
Figure 16. Steady responses with = 2000 r/min. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 16. Steady responses with = 2000 r/min. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g016
Figure 17. Steady responses with = 7480 r/min. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 17. Steady responses with = 7480 r/min. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g017aMachines 12 00059 g017b
Figure 18. Steady responses with = 11,000 r/min. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Figure 18. Steady responses with = 11,000 r/min. (a) Time history, (b) frequency spectra and (c) phase diagram of  q ¯ 1 ; (d) time history, (e) frequency spectra and (f) phase diagram of  q ¯ 2 .
Machines 12 00059 g018aMachines 12 00059 g018b
Figure 19. Damping ratio of the first and second modes versus  δ s  and  k s . (a ζ 1  and (b ζ 2 . The colors represent the magnitudes of  ζ 1  and  ζ 2 .
Figure 19. Damping ratio of the first and second modes versus  δ s  and  k s . (a ζ 1  and (b ζ 2 . The colors represent the magnitudes of  ζ 1  and  ζ 2 .
Machines 12 00059 g019
Figure 20. Effect of  ε σ 1  in the case of combination resonance. (a   α 1  and (b α 2 .
Figure 20. Effect of  ε σ 1  in the case of combination resonance. (a   α 1  and (b α 2 .
Machines 12 00059 g020
Figure 21. Effect of  ε σ 1  on stability. (a) First flapwise mode and (b) second flapwise mode.
Figure 21. Effect of  ε σ 1  on stability. (a) First flapwise mode and (b) second flapwise mode.
Machines 12 00059 g021
Figure 22. Effect of  ε σ 1  in the case of combination resonance with  0.5 k s   c o n t a c t   s t i f f n e s s . (a α 1  and (b α 2 .
Figure 22. Effect of  ε σ 1  in the case of combination resonance with  0.5 k s   c o n t a c t   s t i f f n e s s . (a α 1  and (b α 2 .
Machines 12 00059 g022
Figure 23. Effect of  ε σ 1  in the case of combination resonance with a  2 k s  contact stiffness. (a α 1  and (b α 2 .
Figure 23. Effect of  ε σ 1  in the case of combination resonance with a  2 k s  contact stiffness. (a α 1  and (b α 2 .
Machines 12 00059 g023
Figure 24. Effect of  ε σ 1  in the case of combination resonance with a  0.5 δ s  shroud gap. (a α 1  and (b α 2 .
Figure 24. Effect of  ε σ 1  in the case of combination resonance with a  0.5 δ s  shroud gap. (a α 1  and (b α 2 .
Machines 12 00059 g024
Figure 25. Effect of  ε σ 1  in the case of combination resonance with a  2 δ s  shroud gap. (a α 1  and (b α 2 .
Figure 25. Effect of  ε σ 1  in the case of combination resonance with a  2 δ s  shroud gap. (a α 1  and (b α 2 .
Machines 12 00059 g025
Figure 26. Effect of  ε σ 2  in the case of combination resonance. (a   α 1  and (b   α 2 .
Figure 26. Effect of  ε σ 2  in the case of combination resonance. (a   α 1  and (b   α 2 .
Machines 12 00059 g026
Figure 27. Effect of  ε σ 2  on stability  α 1,2 . (a) The first and (b) second flapwise modes.
Figure 27. Effect of  ε σ 2  on stability  α 1,2 . (a) The first and (b) second flapwise modes.
Machines 12 00059 g027
Figure 28. Effect of  ε σ 2  in the case of combination resonance with  0.5 k s  contact stiffness. (a α 1  and (b α 2 .
Figure 28. Effect of  ε σ 2  in the case of combination resonance with  0.5 k s  contact stiffness. (a α 1  and (b α 2 .
Machines 12 00059 g028
Figure 29. Effect of  ε σ 2  in the case of combination resonance with  2 k s  contact stiffness. (a α 1  and (b α 2 .
Figure 29. Effect of  ε σ 2  in the case of combination resonance with  2 k s  contact stiffness. (a α 1  and (b α 2 .
Machines 12 00059 g029
Figure 30. Effect of  ε σ 2  in the case of combination resonance with a  0.5 δ s    shroud gap. (a α 1  and (b α 2 .
Figure 30. Effect of  ε σ 2  in the case of combination resonance with a  0.5 δ s    shroud gap. (a α 1  and (b α 2 .
Machines 12 00059 g030
Figure 31. Effect of  ε σ 2  in the case of combination resonance with a  2 δ s  shroud gap. (a α 1  and (b α 2 .
Figure 31. Effect of  ε σ 2  in the case of combination resonance with a  2 δ s  shroud gap. (a α 1  and (b α 2 .
Machines 12 00059 g031
Table 1. Major parameters of the blade.
Table 1. Major parameters of the blade.
ParameterValueParameterValue
c v 10 N·s/me0.25 m
EI5.41 × 104 N·m2EA1.08 × 109 N
k s 3.55 MN/mL1.05 m
m l 40.13 kg/m m s 0.2 kg
Q 0 50α60°
δ s 2.0 × 10−5 ml1
μ0.36900 r/min
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, H.; Yuan, G.; Yu, Z.; Wang, Y.; Marmysh, D. Nonlinear Dynamics and Combination Resonance of a Flexible Turbine Blade with Contact and Friction of Shrouds. Machines 2024, 12, 59. https://doi.org/10.3390/machines12010059

AMA Style

Li H, Yuan G, Yu Z, Wang Y, Marmysh D. Nonlinear Dynamics and Combination Resonance of a Flexible Turbine Blade with Contact and Friction of Shrouds. Machines. 2024; 12(1):59. https://doi.org/10.3390/machines12010059

Chicago/Turabian Style

Li, Hua, Gaofei Yuan, Zifeng Yu, Yuefang Wang, and Dzianis Marmysh. 2024. "Nonlinear Dynamics and Combination Resonance of a Flexible Turbine Blade with Contact and Friction of Shrouds" Machines 12, no. 1: 59. https://doi.org/10.3390/machines12010059

APA Style

Li, H., Yuan, G., Yu, Z., Wang, Y., & Marmysh, D. (2024). Nonlinear Dynamics and Combination Resonance of a Flexible Turbine Blade with Contact and Friction of Shrouds. Machines, 12(1), 59. https://doi.org/10.3390/machines12010059

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

Article Metrics

Back to TopTop