Method for profiling blades of an axial turbomachine

Abstract

A method for profiling blades of an axial turbomachine includes preparing a geometric model of a blade profile; determining an oscillation mode of the geometric model; calculating a time profile of a position-dependent disruptive pressure in a channel between two adjacent blade profiles over an oscillation period of the oscillation belonging to the oscillation mode, changing the geometric model and determining a different oscillation mode for the modified geometric model; and determining the damping of the oscillation using the disruptive pressure profile calculated previously and accepting the modified geometric model for the case that the damping of the oscillation turns out to be greater than calculated, otherwise repeating the last two steps with another modified geometric model.

Claims

1. A method for profiling blades of an axial turbomachine comprising: a) preparing a geometric model of a blade profile (R, −1, +1); b) determining an oscillation mode of the geometric model; c) calculating a time profile of a position-dependent disruptive pressure p({right arrow over (x)}, t)=p.sub.ges({right arrow over (x)}, t)−p.sub.0({right arrow over (x)} ) in a channel between two adjacent blade profiles over an oscillation period of an oscillation belonging to the oscillation mode, wherein in order to calculate the profile it is assumed that the position-dependent disruptive pressure p({right arrow over (x)}, t) in the channel is proportional to a magnitude of a surface ({right arrow over (x)}, t) which extends in the channel and is arranged perpendicular to the flow lines between blade surfaces bounding the channel, and perpendicular to these blade surfaces, wherein {right arrow over (x)} is a coordinate on one of the two blade surfaces, t is time, p.sub.0({right arrow over (x)}) is a position-dependent pressure when the blade profile is stationary, and p.sub.ges({right arrow over (x)}, t) is a position-dependent total pressure; d) determining a damping of an oscillation caused by the position-dependent disruptive; e) changing the geometric model and determining a different oscillation mode for the changed geometric model, as well as carrying out step c) with the changed geometric model and the different oscillation mode; and f) determining a damping of an oscillation using the disruptive pressure profile p({right arrow over (x)}, t) calculated in step e) and accepting the changed geometric model for the case that the damping of the oscillation turns out to be greater than damping of the oscillationin step d), otherwise repeating steps e) and f) with another changed geometric model.

2. The method as claimed in claim 1, wherein an integration of the position-dependent disruptive pressure is performed in step d) along the complete channel.

3. The method as claimed in claim 1, further comprising: g) numerically calculating an oscillation, engendered by a flow, of the geometric model accepted in step f) and calculating the damping of an oscillation.

4. The method as claimed in claim 3, wherein a convection and an inertia of the flow are taken into account in step g).

5. The method as claimed in claim 1, wherein the oscillation mode is represented in complex form: {right arrow over (99)}={right arrow over (a)}.sub.c+i{right arrow over (a)}.sub.s, wherein {right arrow over (ϕ)} is the oscillation mode, {right arrow over (a)}.sub.c and {right arrow over (a)}.sub.s are local deformations.

6. The method as claimed in claim 5, wherein {right arrow over (a)}.sub.c and {right arrow over (a)}.sub.s are prepared by means of a finite element structure solver.

7. The method as claimed in claim 5, wherein ({right arrow over (x)}, t) is illustrated according to ({right arrow over (x)}, t)=.sub.0({right arrow over (x)})+Re└(Â.sub.1e.sup.iβ.sup.1e.sup.iσ+Â.sub.Re.sup.iβ.sup.R)e.sup.iωt┘, wherein .sub.0({right arrow over (x)}) is the area ({right arrow over (x)}, t) when the blade profiles are stationary, .sub.R=|{right arrow over (n)}.sub.R.Math.{right arrow over (ϕ)}|.sub.R is a projection of the oscillation mode onto the normal {right arrow over (n)}.sub.R to the blade surface at the coordinate {right arrow over (x)} for the one blade profile, Â.sub.1=|{right arrow over (n)}.sub.1.Math.{right arrow over (ϕ)}|.sub.1 is a projection of the oscillation mode onto the normal {right arrow over (n)}.sub.1 to the blade surface at a point on the blade surface of the other blade profile, β.sub.1 and β.sub.R are the relative phase offset between an imaginary component and a real component, σ is the phase offset of the oscillations of the two blade profiles, and ω is the angular frequency.

8. The method as claimed in claim 7, wherein the damping of the oscillation is determined in step d) for different values of β.sub.1, β.sub.R and σ.

9. The method as claimed in claim 7, wherein .sub.0({right arrow over (x)}) is determined in that a circle is found on the blade surface of the one blade profile at the coordinate {right arrow over (x)} which is tangential to the blade surface of the one blade profile at the coordinate {right arrow over (x)} and the blade surface of the other blade profile at the point, and .sub.0({right arrow over (x)}) is selected in such a way that .sub.0({right arrow over (x)}) is tangential to the normal {right arrow over (n)}.sub.R at the coordinate {right arrow over (x)} and to the normal {right arrow over (n)}.sub.1 at that point.

10. The method as claimed in claim 9, wherein .sub.0({right arrow over (x)}) determined through an interpolation of a function, or through a second-order polynomial, between the coordinate {right arrow over (x)} and the point.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The invention is explained in more detail below with reference to the schematic drawings. Here:

(2) FIG. 1 shows a cross section through three adjacent blades,

(3) FIG. 2 shows a detail from FIG. 1,

(4) FIG. 3 shows an illustration of a channel between two adjacent blades,

(5) FIG. 4 shows an illustration of an oscillation mode of a blade, and

(6) FIG. 5 shows a damping profile over an oscillation period for different reduced frequencies.

DETAILED DESCRIPTION OF INVENTION

(7) An axial turbomachine, such as for example a gas turbine or a steam turbine, comprises arrays of blades. The blades may be guide blades and/or rotor blades which are arranged in a compressor and/or in a turbine. FIG. 1 shows three adjacent, neighboring blades. A geometric model of a blade profile for each of the three blades is illustrated. The middle blade profile is shown as a reference blade profile R, and the two blade profiles neighboring the reference blade profile R are shown as a first neighboring blade profile +1 and as a second neighboring blade profile −1. Each of the blade profiles R, −1, +1 comprises a respective front edge 1 and a rear edge 2. A channel 3 is arranged between each two adjacent blade profiles in which a working fluid flows when the axial turbomachine is operating. The solid lines in FIG. 1 represent the blade profiles R, +1, −1 in their stationary state. The dotted lines represent the oscillation of the blade profiles R, +1, −1, wherein the first neighboring blade profile +1 and the second neighboring blade profile −1 have an identical oscillation phase that is offset from the oscillation phase of the reference blade profile R.

(8) In the method for profiling the blade, in a step a) the geometric model of the blade profile R, −1, +1 is prepared, wherein the blade profiles R, −1, +1 in a row of blades are identical.

(9) The oscillation mode {right arrow over (ϕ)} of the geometric model is determined in a step b). As can be seen from FIG. 4, the oscillation mode {right arrow over (ϕ)} can be represented in complex form according to {right arrow over (ϕ)}={right arrow over (a)}.sub.c+i{right arrow over (a)}.sub.s, wherein {right arrow over (a)}.sub.c and {right arrow over (a)}.sub.s are local deformations. The oscillation mode then results as an overlay of these two modes. This overlay is displayed in FIG. 4 in a vector diagram 4, wherein the phase offset between the real and imaginary components is identified with β. {right arrow over (a)}.sub.c and {right arrow over (a)}.sub.s can be determined by means of a finite element structure solver, wherein the finite element structure solver receives the blade profile as an input.

(10) A temporal progression of a position-dependent disruptive pressure p({right arrow over (x)}, t) in the channel 3 over an oscillation period of the oscillation belonging to the oscillation mode is calculated in a step c). {right arrow over (x)} here is a coordinate on the surface of one of the two adjacent blade profiles, and t is the time. In order to calculate the disruptive pressure p({right arrow over (x)}, t) it is assumed that the position-dependent disruptive pressure p{right arrow over (x)}, t) in the channel 3 is proportional to the magnitude of a surface ({right arrow over (x)}, t) that extends in the channel and is arranged on the blade surfaces bounding the channel 3, perpendicular to these blade surfaces, and perpendicular to the flow lines in between.

(11) The approximation that the disruptive pressure p({right arrow over (x)}, t) is assumed to be proportional to the surface ({right arrow over (x)}, t) is illustrated in FIG. 3. The solid lines in FIG. 3 represent the two blade surfaces bounding the channel 3 in their stationary state. The dotted lines represent the two blade surfaces that bound the channel 3 with a deflection during their oscillation. With a small area ({right arrow over (x)}, t) the speed of the working fluid when the axial turbomachine is operating is high, whereby a low disruptive pressure p({right arrow over (x)}, t) is found. Conversely, it is similarly true that with a large area ({right arrow over (x)}, t) the speed of the working fluid when the axial turbomachine is operating is low, whereby a high disruptive pressure p({right arrow over (x)}, t) is found.

(12) ({right arrow over (x)}, t) is illustrated according to

(13) ({right arrow over (x)}, t)=.sub.0({right arrow over (x)}, t)+Re └(Â.sub.1e.sup.iβ.sup.1e.sup.iσ+Â.sub.Re.sup.iβ.sup.R)e.sup.iωt┘, wherein .sub.0({right arrow over (x)}) is the surface ({right arrow over (x)}, t) when the blade profile is stationary. As shown in FIG. 1, .sub.0({right arrow over (x)}) can be determined in that a circle 12 is found on the blade surface of the one blade profile at the coordinate {right arrow over (x)}, which circle is tangential to the blade surface of the one blade profile at the coordinate {right arrow over (x)} and the blade surface of the other blade profile at a point. For the case that a plurality of the circles 12 could be found, the circle 12 with the shortest diameter is selected. .sub.0({right arrow over (x)}) is then chosen in such a way that .sub.0({right arrow over (x)}) is tangential to the normal {right arrow over (n)}.sub.R to the surface of the one blade profile at the coordinate {right arrow over (x)} and to the normal {right arrow over (n)}.sub.1 to the surface of the other blade profile at the point. .sub.0({right arrow over (x)}) is determined for this purpose by interpolating a function between the coordinate {right arrow over (x)} and the point. The function can, for example, be a second -degree polynomial. As can be seen from FIG. 1, both normals {right arrow over (n)}.sub.R and {right arrow over (n)}.sub.1 point towards the circle center 11 of the circle 12, which is a consequence of the fact that the circle 12 is tangential to the two blade surfaces at the coordinate {right arrow over (x)} and the point.

(14) Â.sub.R=|{right arrow over (n)}.sub.R.Math.{right arrow over (ϕ)}|.sub.R is the projection of the oscillation mode {right arrow over (ϕ)} onto the normal {right arrow over (n)}.sub.R to the blade surface at the coordinate {right arrow over (x)} for the one blade profile. Â.sub.1=|{right arrow over (n)}.sub.1.Math.{right arrow over (ϕ)}|.sub.1 is the projection of the oscillation mode {right arrow over (ϕ)} onto the normal {right arrow over (n)}.sub.1 to the blade surface at the point on the blade surface of the other blade profile. The respective projections Â.sub.R and Â.sub.1 are illustrated in FIG. 2 as the lines from the respective point to the respective end of the lines. B.sub.1 and β.sub.R are the respective phase offset between the imaginary component and the real component. σ is the phase offset of the oscillations of the two blade profiles, and ω is the angular frequency.

(15) The damping of the oscillation caused by the disruptive pressure profile is determined in a step d). The disruptive pressure is multiplied for this purpose with the projection of the oscillation mode onto the normal to the blade surface. The product obtained as a result is integrated over the full channel 3, whereby the damping results.

(16) The geometric model is changed in a step e), and a different oscillation mode determined for the changed geometric model. Step c) is also carried out in step e) with the changed geometric model and the different oscillation mode.

(17) The damping of the oscillation due to the disruptive pressure profile p({right arrow over (x)}, t) calculated in step e) is determined in a step f). The changed geometric model is accepted for the case in which the damping of the oscillation is greater than is found in step d), otherwise the steps e) and f) are repeated with another changed geometric model.

(18) In step d) the damping of the oscillation for different values of β.sub.1, β.sub.R and σ are determined, so that the time-dependency of the surface ({right arrow over (x)}, t) is obtained.

(19) A numerical calculation of the oscillation, engendered by a flow, of the geometric model accepted in step f) and calculation of the damping of the oscillation take place in a step g).

(20) The disruptive pressure p({right arrow over (x)}, t) calculated in step c) and the damping calculated in step d) represent an approximation for the case in which k->0, where k is the reduced frequency. The reduced frequency can be determined according to

(21) $k=\frac{2*\Pi *c*f}{U},$
wherein f is the frequency of the oscillation, c is the chord length of the blade profile, and U is the speed of the working fluid in an axial position of the flow machine in which the front edges 1 of the blade profiles lie. The damping profile over an oscillation period for different reduced frequencies k calculated in step g) is illustrated in FIG. 5. The geometric model on which step g) is based is a geometric model accepted in step f). The phase angle σ is plotted on the abscissa 5, and the damping is plotted on the ordinate 6. A zero line 7 of the damping is also drawn on FIG. 5, wherein the negative values for the damping signify an excitation of the oscillation, and positive values of the damping signify a damping of the oscillation.

(22) FIG. 5 shows a damping value profile 8 at k=0.001, a damping value profile 9 at k=0.1 and a damping value profile 10 at k=0.21. Due to its lower reduced frequency, the damping value profile 8 is comparable to that determined in step d). It can be seen that the damping of the oscillation increases with increasingly reduced frequency. The geometric model accepted in step f) thus represents a good starting point for the subsequent numerical simulation of the oscillation.

(23) Although the invention has been more closely illustrated and described in more detail through the preferred exemplary embodiment, the invention is not restricted by the disclosed examples, and other variations can be derived from this by the person skilled in the art without leaving the scope of protection of the invention.