VIBRATION NUMERICAL ANALYSISOF COUNTER-ROTATING TURBINE WITH WAKE-FLOWUSING FLUID-STRUCTURE INTERACTION METHOD

Zhao Zhenhua,Lu Wenliang,Chen Wei,Wu Tieying

(1.College of Energy and Power Engineering,NUAA,29 Yudao Street,Nanjing,210016,P.R.China;2.The Third Academy,China Aerospace Science and Technology Corporation,Beijing,100074,P.R.China)

INTRODUCTION

The counter-rotating turbine is that,in a twin spool turbo-engine,two rotors rotate in the opposite directions,and there may be no nozzle guide vane(NGV)between the two rotors in order to save weight,to increase thrust-weight ratio and probably to improve engine performance.In recent years,great progresses have been made in technology of the counter-rotating turbine.Studies show that the use of counter-rotating turbine structure can reduce the number of components,lighten the structural quality,improve aerodynamic performance,reduce gyroscopic moment and so on.Counter-rotating turbines have great many advantages over conventional turbines[1-3].

1 METHOD FOR FLUID-STRUCTURE COUPLING CACULATION

Fluid-solid coupling mechanics is a branch of the fluid and solid mechanics interaction,and many researchers have made great contributions in this area.A K Slone[5]solved the problem in the fluid-solid coupling of transient dynamics with the momentum conservation equation of thestructure and the Euler equation of fluid mass,momentum and energy conservation.Zhang Qun and Toshiaki Hisada[6]proposed two kinds of methods for the fluid-solid coupling calculation. V Carstens[7]derived the control equations of gassolid coupled system with the introduction of modal coordinates of the vector aeroelastic equations and 3-DN-Sequations simultaneously.

There are generally two ways to solve the fluid-solid coupling problems.One is the alternate solution method,i.e.,the weak-coupling method;another is the overall solution method,i.e.,the strong-coupling method.

In this paper,the weak-coupling method is applied.The procedure can be described as follows:at moment ti,firstly to solvetheflow field,then to transfer the calculated pressure of the flow field into the solid model for the calculation of the solid domain to find displacement which modifies fluid grid and to re-compute flow field repeatedly using the new grid,and so on.Such a superposition,when the calculation reaches local convergence,the coupling calculation goes to the next step.Calculation process is shown in Fig.1.After getting convergence at ti,computation goes to next step for ti+1 until the final step.

Fig.1 Calculation process

2 NUMARICAL ANALYSIS

2.1 Physical model and boundary conditions

An finite element model(FEM)of a counterrotating turbine is taken as an example.It has a high pressure(HP)NGV stator,a HP rotor and a LP rotor.There are 80 high pressure turbine(HPT)blades in HP rotor and 86 LPT blades in LP rotor(Fig.2).Rotation speed is 14 978 r/min for HPT and 9 426 r/min for LPT.The study is focused on the LPT rotor blade,to analyze the structural strength and vibration properties in the design situation.

Fig.2 Calculation model of CFX

At theinlet of HPNGV,total pressureis P*3 and the distribution of total temperaturevs blade height is shown in Fig.3.Theaverageoutlet static pressure of LPT is P4.

This unsteady coupling calculation takes ideal gas and unsteady Reynolds-averaged N-Sequations as the control equations,using high-precision convection format.The second order backward Euler transient scheme is used in the transient format,and dynamic grid and k-X model are used.

Fig.3 Distribution of inlet temperature vs bladeheight

2.2 Calculation model

The CFX model of fluid calculation has 134 268 elements(Fig.2).Transient coupling calculations and FEM of blade and bound form(the root of solid sticks)are shown in Fig.4.It con-sists 3 648 elements and 18 984 nodes in total.The surfaces of LPT rotor blades are coupling surface,the coupling parameter is pneumatic pressure.The dynamic grid is applied in the region of the LPT rotor blades(Fig.2).

按說這類人應該混得很好，然而現實中往往難有作為。原因很簡單，他們的奉承話在“受眾”看來不值錢。《圍城》中有個褚慎明，常寫信給各國知名哲學家，說自己如何愛讀他們的書，把雜志書評欄里贊美他們的話，改頭換面算自己的意見。哲學家并非歌星、影星，平日里沒什么人“粉”他們，忽然發現在遙遠的東方還有自己的忠實粉絲，高興得險些忘掉了哲學。那些虛榮的老頭子拿著信向同行賣弄，不料彼此都收到這樣一封信，頓時都掃興生氣了……泛濫的奉承話如同嚴重超發的貨幣，一文不值。表哥是經濟學碩士，不知道為什么不明白這個道理。也許受市井家庭影響太大，相對封閉的成長環境多少也限制了他的格局。

Fig.4 FEM of LPT blade

2.3 Approach of temperature load

ANSYS/CFX does not provide function of coupling temperature between FEM and computational fluid dynamics(CFD)model.A program is developed to transfer the temperature from CFD model to FEM.

In one work state of the engine,the temperature in turbine is changing at the different moment since computation is unsteady.However,unsteady temperatures will not have a significant impact on turbine blades because of their thermal latency in the solid.In this case,temperatures in blades can be considered as steady.Considering the most dangerous situation,the temperature field at the moment of top temperature gradient will be exerted on the LPT rotor blade for analysis.

The approach steps are as follows:

(1)Using CFX to calculate unsteady temperaturedistributions of LPT rotor blade without deformation;

(2)Taking the one which has the greatest temperature gradient as example 1;

(3)Transforming this temperature field into the surface nodes of FEM by loading conversion developed by the authors;

(4)Computing the whole temperature distribution for all the nodes of FEM;

(5)Applying temperature results to FEM for further unsteady iteration analysis.

3-D linear interpolation method with the weighted average theory is applied to transform the temperature load.In Fig.5,point P is a node of FEM and the grid is CFD model.Centered on P,a circlewith r0 is defined.Every nodein CFD model in that circle is taken into account for weighted average for point P.

Fig.5 Blade surface

The formula used in conversion program is

where riis the distance between CFD node i and P,Ti the corresponding temperature in point i,and index i the CFD node number of points located in that circle. The program allows user to choose suitable r 0 to get more accurate interpolation results.In addition,this program can be used for other parameters such as pressure,and so on.

Basically,material properties vary with temperatures,like Young′s module.These variations are considered in the computation.

2.4 Calculation result and analysis

Aerodynamic pressure,temperature and centrifugal load areacted on the turbine rotor blade.The LPT rotor blade is calculated used the fluidsolid coupling numerical calculation method under two kinds of blade temperatures(examples 1 and 2)in order to analyze forced vibrations of the LPT rotor blade under aerodynamic loads.The temperatures of the example 2 are globally proportionally lower than thoseof the example 1 con-sidering the cooling effect.

Three points on pressure surface of the blade are selected for analysis,points A and B(Fig.6)for stress response and points B and C for displacement response.Point A is near the hub of the blade and point C near the tip.

Fig.6 Feature point position on LPT blade

Extracting the pressure fluctuation of the place approach point B(Fig.7).In order to analyze the characteristics of wake-flow,the pressure curveis transformed in frequency domain using fast Fourier transform(FFT)method.Fluid excitation force shows broadband characteristics,and contains seven major frequency ranged from 6 006.5 Hz to 39 448.1 Hz.The frequency with maximum amplitudeis 33 459.6 Hz,and frequencies 6 006.5,21 500.7,and 27 489.2 Hz follow(Fig.8).

Fig.7 Pressure curves of place approaching point B

Fig.8 Frequency spectrum of pressure at place approaching point B

Forced response of the LPT rotor blade is obtained with transient coupling method under different temperature conditions, the stress curves of points A and B are shown in Fig.9 and the displacement curves of points B and C are shown in Fig.10.

Fig.9 Stress curves of points A and B under different temperature

Fig.10 Displacement curves of points B and C under different temperature

Form Fig.9 and Fig.10,we can find that:(1)The curves of the displacement are smoother than the curves of stress.Displacement and stress are not proportional.That means the vibration is not in its single natural modal.(2)Temperature has effect on the amplitude of stress response.The amplitude of stress response increases while the temperature decreases,but the mean stress increases with temperature. Here,the mean stress includes centrifugal effect and thermal stress because of temperature gradient in the blade.(3)The amplitude of displacement decreases when temperatureincreases,but the mean displacement increases.

3 SPECTRUM ANALYSIS OF BLADE FORCED RESPONSE

In order to analyze the forced response char-acteristics of LPT blade,the stress and displacement response signals obtained from the above computation of the selected points on the blade under different temperature field are transformed in frequency domain using FFT method.Spectrograms of the stress curves of points A,B and the displacement curves of points B,C of the examples 1,2 areshown in Figs.11-14.The frequency statistics of the selected points is listed in Table 1.

Table 1 Frequency statistics under diff erent temperaturefields Hz

The frequencies of point a in the examples 1,2 are the first natural frequency of the blade under the respectively state(Fig.11).The frequencies of point d in Fig.11(a),point c in Fig.12(a),point d in Fig.13(a)and point d in Fig.14(a)are the seventh natural frequency of blade.The frequencies of point b in Fig.12(b),point c in Fig.13(b),and point b in Fig.14(b)are the second natural frequency of the blade.2 196,2 201 and 3 418,3 420,3 196 Hzin example 1 are frequency-doubled and tripled harmonic.2 321 and 3 760,3 789 Hz in example 2 are also frequency-doubled and tripled harmonic.

Fig.11 Frequency spectrum of stress at point A under different temperature

Fig.12 Frequency spectrum of stress at point B under different temperature

In summary,the forced response of the LPT rotor blade in example 1 is mainly composed by the first and the seventh natural modes.It can be seen through the frequency spectrum that the first natural mode dominates the response.The forced response in example 2 is composed by the first and the second natural modes,and the first mode dominates the response too.

Fig.13 Frequency spectrum of displacement at point B under different temperature

Fig.14 Frequency spectrum of displacement at point C under different temperature

Comparison of the frequency characteristics of wake aerodynamic force and blade vibration shows that:(1)Due to the wide frequency band of exciting force with high energy,the first order vibration modal is easily aroused.(2)The higher order modals may also appear and they are not ignored when the amplitude is big enough.

4 CONCLUSIONS

By using ANSYS/CFX software,the process of computation of fluid-structure coupling is demonstrated successfully.Through this process the counter-rotating turbine is analyzed aerodynamically and structurally. After temperature conversion,vibration is analyzed in the LPT blade under two different temperature distributions.Based on above results,the main conclusions are:

(1)The curves of the displacement are smoother than that of stress.They are not proportional.That means the vibration is not in single natural mode in this operation rating of the turbine.

(2)Under wake aerodynamic force excitation,the first order modal vibrationis more easily aroused and high order response cannot be ignored.

(3)With different temperature fields,natural modals of the blade are different and vibration responses are also different.

(4)The blade vibration response under wake-flow excitation conditions is complex and needs further studies.

[1] Louis J F.Axial flow contra-rotating turbines[R].ASME 85-GB-218,1985.

[2] Liu Siyong,Wang Ping,Fang Xiangjun,et al.Research on aerodynamics design of counter-rotating turbine without intermediate vanes[J].Gas Turbine Exprement and Research,2002,15(1):20-27.(in Chinese)

[3] Ji Lucheng.Review and prospect on research of counter-rotating turbine[J].Aeroengine,2006,32(4):49-53.(in Chinese)

[4] Ji Lucheng.Analysis for technique challenges on counter-rotating turbomachinery[J]. Journal of Propulsion Technology,2007,28(1): 40-44.(in Chinese)

[5] Slone A K.Dynamic fluid-structureinteraction using finite volume unstructured mesh procedures[J].Computer and Structures,2002,80:371-390.

[6] Zhang Qun,Hisada Toshiaki.Studies of the strong coupling and weak coupling methods in FSI analysis[J].Int J Numer Meth Enging,2004(60):2013-2029.

[7] Carstens V,Kemme R,Schmitts S.Coupled simulation of fluid-structure interaction in turbo machinery[J].Aerospace and Technology,2003,7:298-306.