Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 3275750, 12 pages
http://dx.doi.org/10.1155/2016/3275750

Research Article

Global Dynamics of a Compressor Blade with Resonances

Xiaoxia Bian,1,2,3 Fangqi Chen,² and Fengxian An1,2

¹Department of Mechanics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
²Department of Mathematics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
³School of Mathematics and Physics, Yancheng Institute of Technology, Yancheng 224051, China

Received 24 March 2016; Revised 16 June 2016; Accepted 12 July 2016

Academic Editor: Jaromir Horacek

Copyright © 2016 Xiaoxia Bian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The global bifurcations and chaotic dynamics of a thin-walled compressor blade for the resonant case of 2 : 1 internal resonance and primary resonance are investigated. With the aid of the normal theory, the desired form associated with a double zero and a pair of pure imaginary eigenvalues for the global perturbation method is obtained. Based on the simpler form, the method developed by Kovacic and Wiggins is used to find the existence of a Shilnikov-type homoclinic orbit. The results obtained here indicate that the orbit homoclinic to certain invariant sets for the resonance case which may lead to chaos in the sense of Smale horseshoes for the system. The chaotic motions of the rotating compressor blade are also found by using numerical simulation.

1. Introduction

Compressor blades are widely used in many fields of aerospace, aeronautic engineering, and mechanical industry due to their excellent mechanical properties. The problem of nonlinear dynamics of the rotating blades had attracted lots of research interest during the past decade. Various strategies and approaches have been proposed for nonlinear dynamics of rotating blades (see, e.g., [1–15]). However, theoretical analysis of global dynamics of the rotating blades has not been concerned in the current available literature. Several researchers have examined the global behaviors of plates, beams, and belt (see, e.g., [16–22]), but the results cannot be directly extended to the case of rotating blades.

Yang and Tsao investigated the vibration and stability of a pretwisted blade under nonconstant rotating speed in [1], and they also predicted the time-dependent rotating speed leads to a system with six parametric instability regions in primary and combination resonances. Surace et al. [2] dealt with the coupled bending-bending-torsion vibration of rotating pretwisted blades. Şakar and Sabuncu [3] presented the static stability and the dynamic stability of an aerofoil cross section rotating blade subjected to an axial periodic force and took into account the effects of coupling due to the center of flexure distance from the centroid, rotational speed, disk radius, and stagger angle. Al-Bedoor and Al-Qaisia [4] used a reduced-order nonlinear dynamic model to research the steady-state response of the rotating blade under the main shaft torsional vibration. Tang and Dowell [5] analyzed the nonlinear response of a nonrotating flexible rotor blade subjected to periodic gust excitations theoretically and experimentally. They reported that there exists a periodic or possibly chaotic behavior in the blade. Choi and Chou [6] studied the dynamic response of turbomachinery blades with general end restraints by applying the modified differential quadrature method. A Monte Carlo approach was employed to explore a supercritical Hopf bifurcation and random bifurcation of a two-dimensional nonlinear airfoil in turbulent flow by Poirel and Price [7]. Lacarbonara et al. [8, 9] established the governing equations of the blades under the centrifugal forces and discussed linear modal properties and the nonlinear modes of vibration away from internal resonances, respectively. Yao et al. [10] performed a nonlinear dynamic analysis of the rotating blade with varying rotating speed under high-temperature supersonic gas flow; furthermore, they [11] explored the contributions of nonlinearity, damping, and rotating speed to the steady-state nonlinear responses of the rotating blade, and they also investigated the effects of the rotating speed on nonlinear oscillations of the blade. Wang and Zhang discussed the stability of a spinning blade having periodically time varying coefficients for both linear model and geometric nonlinear model and obtained the stability boundary of linear model and stability of steady-state solutions of nonlinear model in [12].

In many cases, blades are usually modeled as a pretwisted, presetting, thin-walled rotating cantilever beam because the shape of the blade is very complex. Many researchers carried out studies on the dynamic behavior of the beam of this kind and obtained a lot of valuable results (see, e.g., [13–15]). Several methods have been developed to research the global bifurcation behaviors and chaotic dynamics in nonlinear systems that possess homoclinic or heteroclinic orbits. There are three methods: Melnikov method, global perturbation method, and energy-phase method. Melnikov gave the condition under which a homoclinic orbit in the unperturbation system would break under perturbation and at last lead to chaos in the system. Based on Melnikov method, Wiggins studied the global behaviors of the three basic systems [23]. Then, Kovacic and Wiggins [24] developed the global perturbation method to present Shilnikov-type homoclinic orbit for resonant system. The energy-phase method proposed by Haller and Wiggins [25, 26] detected the existence of single-pulse and multipulse homoclinic orbits in a class of near Hamilton systems. Applying the latter two methods, there were many applications to investigate the global behaviors (see, e.g., [16–22]).

In this paper, we obtain a sufficient condition for the existence of Shilnikov-type homoclinic orbit of a compressor blade with 2 : 1 internal resonance and primary resonance using normal form theory and global perturbation method. Firstly, the formulas of the simpler normal form associated with a double zero and a pair of pure imaginary eigenvalues are derived by normal form theory in Section 2. Then, the dynamics of unperturbed system and perturbed system are analyzed in Sections 3 and 4 in detail, respectively. The analysis indicate that Shilnikov-type homoclinic orbit exists in these cases. Finally, numerical simulations are given to confirm the result in Section 5 and the work ends in Section 6 with a short summary.

2. Formulation of the Problem

A thin-walled compressor blade of gas turbine engines with varying speed under high-temperature supersonic gas flow is considered in [11]. It is modeled as a pretwisted, presetting, thin-walled rotating cantilever beam, considering the geometric nonlinearity, centrifugal force, the aerodynamic load, and the perturbed angular speed.

The pretwisted flexible cantilever blade, with length mounted on a rigid hub with radius , is considered [11]. It rotates at a varying rotating speed around its polar axis where , where is the rotating speed at the steady-state and is a periodic perturbation. It is also allowed to vibrate flexurally in the plane making an angle , as shown in Figure 1(a). The rotating blade is treated as a pretwisted, presetting, thin-walled rotating cantilever beam. The length and width of the cross section of the beam in the and directions are and , respectively, and the thickness of the thin-walled beam is . For the purpose of describing the motion of the rotating blade, different coordinate systems are needed. The origin of the rotating coordinate system is located at the blade root, and are the principal axes of an arbitrary beam cross section in the local coordinates (Figure 1(b)), and the transformations between two coordinate systems are shown as and , . is denoted as the pretwist at the beam tip; then, is the pretwist angle of a current beam cross section. The local coordinate system is defined on the cross section of the beam to describe the geometric configuration and the cross section, where and are the circumferential and thickness coordinate variables in Figure 1(c); the notion represents the points off the middle surface; it is different from the notion ; the relationship is and . Assume that and represent the displacements of an arbitrary point and a point in the middle surface of the rotating blades on the , , and directions, respectively. and represent the rotations about the - and -axis, respectively.

(a)

(b)

(c)

Figure 1: The rotating cantilever beam mode: (a) general view with the global coordinate system, (b) cross section of the rotating beam, and (c) local coordinate systems of the cross section.

Based on the isotropic constitutive law, the nonlinear partial differential governing equations of motion for the pretwist, presetting, thin-walled rotating cantilever beam were derived by using Hamilton’s principle in [11]. Then, Galerkin procedure was applied to obtain the dimensionless governing differential equations of nonlinear vibration for the rotating blade by Yao et al. as follows:

where and are the amplitudes of normal modes, and are normal frequencies, and , , , and are damping parameters. plays the role of the nonlinearity, is the amplitude of excitation, and all the expressions of the coefficients can be found in [11].

We study the case of 2 : 1 internal resonance and primary resonance; the resonant relations are represented as and , , where and and are two detuning parameters.

Using the method of multiple scales, the averaged equations were obtained as follows [11]:

Equations (2) have a zero solution . Without the perturbation parameter , the Jacobian matrix at the origin is

The characteristic equation corresponding to the zero solution is

For convenience of the following analysis, let , , , , , and . When , , , and are simultaneously satisfied, the eigenvalues of system (2) without parameter have a nonsemisimple double zero and a pair of pure imaginary eigenvalues and , where . Assume , , , and , considering , , , and as the perturbation parameters; then, (2) without the perturbation parameters becomes

In this case, we have

Using the method in [27], a third-order normal form of (5) is obtained as

Normal form with perturbation parameters of system (2) is

We need to transform system (8) to a desired form in order to apply the global perturbation method. Let and , and use the transformations

And substituting (9) into the normal form (8) yields

where the Hamiltonian functions and are of the following form:

where .

3. Dynamics of the Unperturbed System

Setting in system (10), we obtain the unperturbed system. Obviously, the variable is a constant since , and the first three equations are completely independent of . Thus, we obtain two uncoupled single-degree-of-freedom nonlinear systems:

All possible fixed points in phase space can be classified as

where ; that is, as or . When , the only solution of system (13) is , and from the Jacobian matrix evaluated at the trivial solution, is a saddle point. At , the trivial solution may bifurcate into three solutions through a pitchfork bifurcation. From the Jacobian matrices evaluated at and , it is known that is a center and are two saddle points. The phase portrait is illustrated in Figure 2.

Figure 2: Phase portrait of the unperturbed system in the space.

From transformation (9), the variables and may actually represent the amplitude and phase of nonlinear oscillations. Therefore, assume that variable and put and , such that, for all , system (13) has two saddle points and one center , which is connected by heteroclinic orbits . In four-dimensional space , the set defined by

is a two-dimensional invariant manifold and it is normally hyperbolic [24].

has a three-dimensional stable manifold and an unstable manifold . Then, from [24], the existence of the heteroclinic orbits implies that and intersect nontransversally along a three-dimensional heteroclinic manifold denoted by , which can be written as

where is a constant determined by the initial conditions. The geometric structures of the stable and unstable manifolds of and are shown in Figure 3. It is seen that (13) is a Hamilton system with Hamiltonian

where . Then, we get the expressions of the pair of heteroclinic orbits as follows:

The dynamics restricted to the invariant manifold are described by the following equations:

where . From (19b), we have periodic orbits which are circles for each when , and the corresponding circle is a circle of fixed points when ; that is, As homoclinic orbit in plane is a heteroclinic connection in the four-dimensional space which is shown in Figure 4, is called resonance due to the vanishing frequency of rotation along the direction, and when , is not defined.

Figure 3: The geometric structures of and .

Figure 4: The geometry of trajectories homoclinic to the periodic orbits on and orbits heteroclinic to fixed points on the resonances.

Now, we consider the phase shift:

Substituting (18) into (19b) yields

Integrating (21) yields

Thus, the phase shift is expressed as

is a function of . It is illustrated in Figure 5.

Figure 5: The phase shift defined in (23).

4. Dynamics of the Perturbed System

As the manifold along with its stable manifold and unstable manifold is invariant under sufficiently small perturbations [24], under perturbation (when ), becomes a locally invariant two-dimensional manifold described as follows:

The flow on is obtained by substituting into (10):

Introduce the scale transformations

Substituting transformations (26) into (25) yields

When , (27) is reduced to

which is a Hamilton system with Hamiltonian

The fixed points of Hamilton system (28) are given by

The Jacobian matrix of (28) evaluated at these fixed points is

It is easy to find that is a saddle point and is a center. Therefore, there exists a homoclinic orbit connecting to itself. The phase portrait of system (28) is shown in Figure 6. By the analysis of Kovacic and Wiggins [24], we can obtain that, for sufficiently small , remains a saddle point and becomes a hyperbolic sink .

Figure 6: Dynamics of unperturbed system (28) on .

The phase portrait of perturbed system (27) is given in Figure 7; Hamilton function remains constant on the homoclinic orbit; that is, ; then, we have

Figure 7: Dynamics of perturbed system (27) on .

In order to consider the dynamics on in the neighborhood of , an annulus is defined as

where is a constant, which is chosen sufficient large so that the unperturbed homoclinic orbits are enclosed within the annulus. Denote and as the three-dimensional stable and unstable manifolds of , which are subsets of and , respectively. According to the analysis of [24], the existence of an orbit homoclinic to a saddle-focus point can lead to chaos. This type of homoclinic orbit is called Shilnikov-type homoclinic orbit. The point on has an orbit that comes out of in the four-dimensional space and may return to the annulus; it may approach asymptotically as and eventually complete a Shilnikov-type homoclinic orbit as shown in Figure 8.

Figure 8: Shilnikov-type homoclinic orbit to .

We need to confirm the existence of a Shilnikov type homoclinic orbit in two steps. First, we show by using Melnikov theory when Melnikov function has a simple zero. Second, we determine whether or not the trajectory in comes back in the domain of attraction of . Based on [24], the higher dimensional Melnikov function is given as

Using the division of integral method and the abovementioned analysis, (34) can be expressed as

With the aforementioned analysis, the first term can be evaluated as

The second term can be simplified as

The third term is changed into

By (36), (37), and (38), the Melnikov function may be expressed as

Now, we can require that the Melnikov function has a simple zero. That is, we require

Next, we examine whether the orbit on returns to the domain of attraction of . The condition is given by

where is an integer, , , and are given by (30) and (31), and is the change of angle. According to [24], when conditions (40) and (41) are satisfied simultaneously, there exists the Shilnikov-type chaos in the sense of Smale horseshoes in system (2).

5. Numerical Simulation of Chaotic Motions

Now fixed parameters are used in the abovementioned theory to simplify the calculation. Letting

condition (41) becomes

From (23), is a function of ; then, is a function of . The figure of shows that exists when , so Melnikov function has a simple zero (Figure 9). , , and are presented in Figure 10; we can see ; that is, condition (41) is satisfied. Then, has a Shilnikov homoclinic orbit for sufficiently small . We choose (1) and (2) to do numerical simulations. We use numerical approach to explore the existence of chaotic motions of the rotating thin-walled blade. In Figure 10, we show the existence of the chaotic responses of the thin-walled blade to the forcing excitation. , and other parameters and initial conditions were chosen as , , , , , , , , , , , and . Figure 11 shows the phase portraits on the planes and the wave forms on plane based on (2). With the same parameters, we get the portraits on the planes and the wave forms on plane based on (1). They are shown in Figure 12; the chaotic motion demonstrated in Figures 11 and 12 is Shilnikov-type multipulse chaotic motion. Therefore, the numerical results agree with the theoretical predictions qualitatively.

Figure 9: The zeros of Melnikov’s function.

Figure 10: Graphs of , , and .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11: The chaotic motions of the compressor blade based on (2): (a) the phase portrait on plane (), (b) the waveform on plane (), (c) the phase portrait on plane (), (d) the waveform on plane (), (e) the phase portrait in three-dimensional space (), and (f) the phase portrait in three-dimensional space ().

(a)

(b)

(c)

(d)

(e)

(f)

Figure 12: The chaotic motions of the compressor blade based on (1): (a) the phase portrait on plane (), (b) the waveform on plane (), (c) the phase portrait on plane (), (d) the waveform on plane (), (e) the phase portrait in three-dimensional space (), and (f) the phase portrait in three-dimensional space ().

6. Conclusions

The global bifurcations and chaotic dynamics of the thin-walled compressor blade with varying speed are investigated for the first time by using the analytical and numerical approaches simultaneously when the averaged equations have one nonsemisimple double zero and a pair of pure imaginary eigenvalues. The study is focused on coexistence of 2 : 1 internal resonance and primary resonance. Normal theory is utilized to find the explicit expressions of the simpler normal form of the averaged equations with a double zero and a pair of pure imaginary eigenvalues. Based on the Melnikov method and its extensions to resonance cases developed by Kovacic and Wiggins, the thin-walled compressor blade can undergo homoclinic bifurcation and the Shilnikov-type homoclinic orbit; that is, there exists chaotic motion in full four-dimensional averaged system. Finally, the Dynamics software is used to perform numerical simulation. The numerical results show the existence of chaotic motions in the averaged equations, which illustrate the predictions obtained by the theoretical analysis. The chaotic motions in averaged equations can lead to the amplitude modulated chaotic oscillations in the original system under certain conditions. Therefore, there are Shilnikov-type single-pulse chaotic motions for the thin-walled rotating compressor blade. This is the extension of the results obtained by Yao et al. [11]. We believe that our results give a direct explanation for the jumping behaviors observed in this class of the compressor blade under in-plane and moment excitations.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (11572148).

References

S. M. Yang and S. M. Tsao, “Dynamics of a pretwisted blade under nonconstant rotating speed,” Computers & Structures, vol. 62, no. 4, pp. 643–651, 1997.

G. Surace, V. Anghel, and C. Mares, “Coupled bending-bending-torsion vibration analysis of rotating pretwisted blades: an integral formulation and numerical examples,” Journal of Sound and Vibration, vol. 206, no. 4, pp. 473–486, 1997.

G. Şakar and M. Sabuncu, “Dynamic stability of a rotating asymmetric cross-section blade subjected to an axial periodic force,” International Journal of Mechanical Sciences, vol. 45, no. 9, pp. 1467–1482, 2003.

B. O. Al-Bedoor and A. A. Al-Qaisia, “Stability analysis of rotating blade bending vibration due to torsional excitation,” Journal of Sound and Vibration, vol. 282, no. 3–5, pp. 1065–1083, 2005.

D. M. Tang and E. H. Dowell, “Nonlinear response of a non-rotating rotor blade to a periodic gust,” Journal of Fluids and Structures, vol. 10, no. 7, pp. 721–742, 1996.

S.-T. Choi and Y.-T. Chou, “Vibration analysis of elastically supported turbomachinery blades by the modified differential quadrature method,” Journal of Sound and Vibration, vol. 240, no. 5, pp. 937–953, 2001.

D. Poirel and S. J. Price, “Bifurcation characteristics of a two-dimensional structurally non-linear airfoil in turbulent flow,” Nonlinear Dynamics, vol. 48, no. 4, pp. 423–435, 2007.

W. Lacarbonara, H. Arvin, and F. Bakhtiari-Nejad, “A geometrically exact approach to the overall dynamics of elastic rotating blades—part 1: linear modal properties,” Nonlinear Dynamics, vol. 70, no. 1, pp. 659–675, 2012.

H. Arvin, W. Lacarbonara, and F. B. Nejad, “A geometrically exact approach to the overall dynamics of elastic rotating blades part II: flapping nonlinear normal modes,” Nonlinear Dynamics, vol. 70, no. 3, pp. 2279–2301, 2012.

M. H. Yao, Y. P. Chen, and W. Zhang, “Nonlinear vibrations of blade with varying rotating speed,” Nonlinear Dynamics, vol. 68, no. 4, pp. 487–504, 2012.

M. H. Yao, W. Zhang, and Y. P. Chen, “Analysis on nonlinear oscillations and resonant responses of a compressor blade,” Acta Mechanica, vol. 225, no. 12, pp. 3483–3510, 2014.

F. Wang and W. Zhang, “Stability analysis of a nonlinear rotating blade with torsional vibrations,” Journal of Sound and Vibration, vol. 331, no. 26, pp. 5755–5773, 2012.

S.-Y. Oh, O. Song, and L. Librescu, “Effects of pretwist and presetting on coupled bending vibrations of rotating thin-walled composite beams,” International Journal of Solids & Structures, vol. 40, no. 5, pp. 1203–1224, 2003.

F. S. M. Jarrar and M. N. Hamdan, “Nonlinear vibrations and buckling of a flexible rotating beam: a prescribed torque approach,” Mechanism & Machine Theory, vol. 42, no. 8, pp. 919–939, 2007.

A. Warminska, E. Manoach, J. Warminski, and S. Samborski, “Regular and chaotic oscillations of a Timoshenko beam subjected to mechanical and thermal loadings,” Continuum Mechanics and Thermodynamics, vol. 27, no. 4-5, pp. 719–737, 2015.

Z. C. Feng and P. R. Sethna, “Global bifurcations in the motion of parametrically excited thin plates,” Nonlinear Dynamics, vol. 4, no. 4, pp. 389–408, 1993.

W. Zhang and M. H. Yao, “Theories of multi-pulse global bifurcations for high-dimensional systems and application to cantilever beam,” International Journal of Modern Physics B, vol. 22, no. 24, pp. 4089–4141, 2008.

W. Q. Yu and F. Q. Chen, “Global bifurcations and chaos in externally excited cyclic systems,” Communications in Nonlinear Science & Numerical Simulation, vol. 15, no. 12, pp. 4007–4019, 2010.

W. Zhang, M. H. Yao, and J. H. Zhang, “Using the extended Melnikov method to study the multi-pulse global bifurcations and chaos of a cantilever beam,” Journal of Sound and Vibration, vol. 319, no. 1-2, pp. 541–569, 2009.

M. H. Yao, W. Zhang, and J. W. Zu, “Multi-pulse chaotic dynamics in non-planar motion of parametrically excited viscoelastic moving belt,” Journal of Sound and Vibration, vol. 331, no. 11, pp. 2624–2653, 2012.

M. H. Yao and W. Zhang, “Multi-pulse chaotic motions of high-dimension nonlinear system for a laminated composite piezoelectric rectangular plate,” Meccanica, vol. 49, no. 2, pp. 365–392, 2014.

M. H. Yao, W. Zhang, and D. M. Wang, “Modeling and chaotic dynamics of the laminated composite piezoelectric rectangular plate,” Mathematical Problems in Engineering, vol. 2014, Article ID 345072, 19 pages, 2014.

S. Wiggins, Global Bifurcations and Chaos-Analytial Methods, vol. 73 of Applied Mathematical Sciences, Springer, 1988.

G. Kovacic and S. Wiggins, “Orbits homoclinic to resonances, with an application to chaos in a model of the forced and damped sine-Gordon equation,” Physica D: Nonlinear Phenomena, vol. 57, no. 1-2, pp. 185–225, 1992.

G. Haller and S. Wiggins, “N-pulse homoclinic orbits in perturbations of resonant hamiltonian systems,” Archive for Rational Mechanics & Analysis, vol. 130, no. 1, pp. 25–101, 1995.

G. Haller, Chaos Near Resonance, vol. 138 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1999.

W. Zhang, F. X. Wang, and J. W. Zu, “Computation of normal forms for high dimensional non-linear systems and application to non-planar non-linear oscillations of a cantilever beam,” Journal of Sound and Vibration, vol. 278, no. 4-5, pp. 949–974, 2004.