CN107766598A - Turbine optimum timing location determining method and device - Google Patents

Turbine optimum timing location determining method and device Download PDF

Info

Publication number
CN107766598A
CN107766598A CN201610692257.7A CN201610692257A CN107766598A CN 107766598 A CN107766598 A CN 107766598A CN 201610692257 A CN201610692257 A CN 201610692257A CN 107766598 A CN107766598 A CN 107766598A
Authority
CN
China
Prior art keywords
tail
movable vane
vane
circumferential
upstream
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610692257.7A
Other languages
Chinese (zh)
Inventor
张宝
严红明
谭智勇
郑建弘
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AECC Commercial Aircraft Engine Co Ltd
Original Assignee
AECC Commercial Aircraft Engine Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by AECC Commercial Aircraft Engine Co Ltd filed Critical AECC Commercial Aircraft Engine Co Ltd
Priority to CN201610692257.7A priority Critical patent/CN107766598A/en
Publication of CN107766598A publication Critical patent/CN107766598A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)

Abstract

The present invention discloses a kind of turbine optimum timing location determining method and device.This method includes：Upstream vane tail is obtained in downstream transport process along journey VELOCITY DISTRIBUTION by steady numerical simulation result；According to described circumferential total drift distance of the upstream vane tail in downstream stator porch is determined along journey VELOCITY DISTRIBUTION；Turbine optimum timing position is obtained according to the circumferential total drift distance.The present invention by permanent numerical computations inverting tail downstream transport process along journey VELOCITY DISTRIBUTION, and combine the engineering model proposed and calculate tail in downstream stator porch circumferential offset distance, turbine optimum timing position can be obtained by the offset distance, so as to avoid substantial amounts of unsteady numerical simulations process.

Description

Turbine optimum timing location determining method and device
Technical field
The present invention relates to turbine field, more particularly to a kind of turbine optimum timing location determining method and device.
Background technology
High aeroperformance is typically pursued in turbine pneumatic design, and with the development of design method, increasingly show permanent The limitation of design.It is understood that turbine interior flow field substantially has height non-stationarity, and fully excavate unsteady The pneumatic income that effect is brought is new breakthrough mouth.Utilize the unsteady aerodynamic effect in flow field can by being laid out the timing position of turbine Further to improve the isentropic efficiency of turbine.
Different timing positions influence difference to turbine isentropic efficiency, an optimum timing position be present.It is optimal to find this Timing position, prior art typically carry out unsteady numerical simulations to turbine difference timing position and then obtain turbine optimum timing Position.But there are following several respects in this method：First, current unsteady computation may relate to blade reduction problem, mistake Journey is cumbersome；2nd, the unsteady computation of multiple examples requires higher to computer resource；3rd, numerical simulation, which iterates, easily makes Grown into calculating cycle.
The content of the invention
In view of above technical problem, can the invention provides a kind of turbine optimum timing location determining method and device Turbine optimum timing position is searched out with simple and fast, substantial amounts of unsteady numerical simulations process can be avoided.
According to an aspect of the present invention, there is provided a kind of turbine optimum timing location determining method, including：
Upstream vane tail is obtained in downstream transport process along journey VELOCITY DISTRIBUTION by steady numerical simulation result；
According to described circumferential total drift distance of the upstream vane tail in downstream stator porch is determined along journey VELOCITY DISTRIBUTION；
Turbine optimum timing position is obtained according to the circumferential total drift distance.
In one embodiment of the invention, it is described to determine upstream vane tail in downstream along journey VELOCITY DISTRIBUTION according to described The circumferential total drift distance of stator porch includes：
Determine that upstream vane tail is led by the very first time and upstream of movable vane passage along journey VELOCITY DISTRIBUTION according to described Leaf tail mark passes through movable vane and the second time of downstream stator axial direction spacing；
The circumferential total drift distance is determined according to the very first time and second time.
It is in one embodiment of the invention, described that along journey VELOCITY DISTRIBUTION to determine that upstream vane tail passes through dynamic according to described The very first time of leaf passage includes：
In axial direction take in movable vane passage at least 3 points, wherein include movable vane porch, movable vane at described at least 3 points Exit and at least one movable vane passage intermediate point；
Obtain tail axial velocity value at described from steady numerical simulation result at least 3 points；
According in tail axial velocity value, movable vane axial direction chord length, at least one movable vane passage at least three point Between point to movable vane leading edge axial line distance, determine the very first time.
It is in one embodiment of the invention, described that along journey VELOCITY DISTRIBUTION to determine that upstream vane tail passes through dynamic according to described Leaf and the second time of downstream stator axial direction spacing include：
In axial direction take in movable vane passage at 2 points, include movable vane exit and downstream stator entrance at described 2 points Place；
Tail axial velocity value at described 2 points is obtained from steady numerical simulation result；
According to the axial spacing of tail axial velocity value, movable vane and downstream stator at 2 point, when determining described second Between.
In one embodiment of the invention, it is described that the circumference is determined according to the very first time and second time Total drift distance includes：
The first of upstream vane tail is determined according to the axis spacing of upstream vane and movable vane, upstream vane outlet geometry angle Stage circumferential offset distance, wherein, the first stage circumferential offset distance reaches movable vane porch for upstream vane tail Circumferential offset distance；
The of upstream vane tail is determined according to movable vane botanical origin radial direction radius, turbine rotating speed and the very first time Two-stage circumferential offset distance, wherein, second stage circumferential offset distance is upstream vane tail by movable vane passage, is arrived Up to the circumferential offset distance in movable vane exit；
According to the axis of movable vane botanical origin radial direction radius, turbine rotating speed, movable vane outlet geometry angle, movable vane and downstream stator Spacing and second time determine the phase III circumferential offset distance of upstream vane tail, wherein, week phase III It is circumferential offset distance of the upstream vane tail from movable vane exit arrival downstream stator entrance to offset distance；
By first stage circumferential offset distance, second stage circumferential offset distance and phase III circumferential offset distance And be used as the circumferential total drift distance.
In one embodiment of the invention, it is described that turbine optimum timing position is obtained according to the circumferential total drift distance Put including：
Determined according to circumferential total drift distance, downstream stator botanical origin radial direction radius and the downstream stator lobe numbers Turbine optimum timing position, wherein the turbine optimum timing position is the circumferentially opposed position of downstream stator and upstream vane Put.
According to another aspect of the present invention, there is provided a kind of turbine optimum timing position determining means, including speed obtain Module, total drift acquisition module and optimum timing acquisition module, wherein：
Speed acquiring module, for obtaining upstream vane tail in downstream transport process by steady numerical simulation result Along journey VELOCITY DISTRIBUTION；
Total drift acquisition module, for determining upstream vane tail in downstream stator entrance along journey VELOCITY DISTRIBUTION according to described The circumferential total drift distance at place；
Optimum timing acquisition module, for obtaining turbine optimum timing position according to the circumferential total drift distance.
In one embodiment of the invention, the total drift acquisition module includes time acquisition unit and total drift obtains Unit, wherein：
Time acquisition unit, for determining upstream vane tail by movable vane passage along journey VELOCITY DISTRIBUTION according to described One time and upstream vane tail pass through movable vane and the second time of downstream stator axial direction spacing；
Total drift acquiring unit, for according to the very first time and second time determine the circumferential total drift away from From.
In one embodiment of the invention, the time acquisition unit includes very first time acquisition submodule, wherein：
Very first time acquisition submodule, in axial direction taking in movable vane passage at least 3 points, wherein it is described at least 3 points include movable vane porch, movable vane exit and at least one movable vane passage intermediate point；From steady numerical simulation result Tail axial velocity value at least three points described in obtaining；According to tail axial velocity value, movable vane shaft orientation string at least three point Long, described at least one movable vane passage intermediate point determines the very first time to the axial line distance of movable vane leading edge.
In one embodiment of the invention, the time acquisition unit includes the second time acquisition submodule, wherein：
Second time acquisition submodule, in axial direction taking in movable vane passage at 2 points, described 2 points include movable vane Exit and downstream stator porch；Tail axial velocity value at described 2 points is obtained from steady numerical simulation result；Root According to the axial spacing of tail axial velocity value, movable vane and downstream stator at 2 point, second time is determined.
In one embodiment of the invention, the total drift acquiring unit includes the first skew acquisition submodule, second Acquisition submodule, the 3rd skew acquisition submodule and total drift acquisition submodule are offset, wherein：
First skew acquisition submodule, for the axis spacing according to upstream vane and movable vane, upstream vane outlet geometry Angle determines the first stage circumferential offset distance of upstream vane tail, wherein, the first stage circumferential offset distance is upstream Stator tail reaches the circumferential offset distance of movable vane porch；
Second skew acquisition submodule, for according to movable vane botanical origin radial direction radius, turbine rotating speed and it is described first when Between determine the second stage circumferential offset distance of upstream vane tail, wherein, second stage circumferential offset distance is upstream Stator tail passes through movable vane passage, the circumferential offset distance in arrival movable vane exit；
3rd skew acquisition submodule, for exporting geometry according to movable vane botanical origin radial direction radius, turbine rotating speed, movable vane The axis spacing at angle, movable vane and downstream stator and second time determine the phase III circumferential offset of upstream vane tail Distance, wherein, the phase III circumferential offset distance reaches downstream stator entrance for upstream vane tail from movable vane exit Circumferential offset distance；
Total drift acquisition submodule, for by the first stage circumferential offset distance, second stage circumferential offset distance With phase III circumferential offset distance and be used as the circumferential total drift distance.
In one embodiment of the invention, the optimum timing acquisition module be used for according to the circumferential total drift away from Turbine optimum timing position is determined from, downstream stator botanical origin radial direction radius and downstream stator lobe numbers, wherein described Turbine optimum timing position is the circumferential relative position of downstream stator and upstream vane.
The present invention can by permanent numerical computations inverting tail downstream transport process along journey VELOCITY DISTRIBUTION, and combine The engineering model of it is proposed calculates tail in downstream stator porch circumferential offset distance, and turbine can be obtained by the offset distance Optimum timing position, so as to avoid substantial amounts of unsteady numerical simulations process.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the schematic diagram of turbine optimum timing location determining method one embodiment of the present invention.
Fig. 2 is the schematic diagram of the unsteady matching engineering model of turbine edge line in one embodiment of the invention.
Fig. 3 is according to the schematic diagram that circumferential total drift distance is determined along journey VELOCITY DISTRIBUTION in one embodiment of the invention.
Fig. 4 be one embodiment of the invention in tail movable vane passage interior edge axial velocity profile schematic diagram.
Fig. 5 be another embodiment of the present invention in tail movable vane passage interior edge axial velocity profile schematic diagram.
Fig. 6 is the signal for determining circumferential total drift distance in one embodiment of the invention according to the very first time and the second time Figure.
Fig. 7 is the schematic diagram of turbine optimum timing position determining means one embodiment of the present invention.
Fig. 8 is the schematic diagram of total drift acquisition module in one embodiment of the invention.
Fig. 9 is the schematic diagram of total drift acquisition module in another embodiment of the present invention.
Figure 10 is the schematic diagram of total drift acquiring unit in one embodiment of the invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete Site preparation describes, it is clear that described embodiment is only part of the embodiment of the present invention, rather than whole embodiments.Below Description only actually at least one exemplary embodiment is illustrative, is never used as to the present invention and its application or makes Any restrictions.Based on the embodiment in the present invention, those of ordinary skill in the art are not making creative work premise Lower obtained every other embodiment, belongs to the scope of protection of the invention.
Unless specifically stated otherwise, the part and positioned opposite, the digital table of step otherwise illustrated in these embodiments Do not limited the scope of the invention up to formula and numerical value.
Simultaneously, it should be appreciated that for the ease of description, the size of the various pieces shown in accompanying drawing is not according to reality Proportionate relationship draw.
It may be not discussed in detail for technology, method and apparatus known to person of ordinary skill in the relevant, but suitable In the case of, the technology, method and apparatus should be considered as authorizing part for specification.
In shown here and discussion all examples, any occurrence should be construed as merely exemplary, without It is as limitation.Therefore, the other examples of exemplary embodiment can have different values.
It should be noted that：Similar label and letter represents similar terms in following accompanying drawing, therefore, once a certain Xiang Yi It is defined, then it need not be further discussed in subsequent accompanying drawing in individual accompanying drawing.
Time (clocking effect), i.e. turbine isentropic efficiency arrange circumferentially opposed position with adjacent to and opposite static leaf Put change and change.
It is found by the applicant that：The essence of time is exactly to allow low energy fluid to develop close proximity to blade surface so as to improve leaf Turbine isentropic efficiency.Based on above-mentioned essence, if can obtain when tail reaches downstream stator porch with respect to upstream vane trailing edge Circumferential offset distance, downstream stator blade can be displaced in relevant position, realize and improve turbine efficiency using time Purpose.
Fig. 1 is the schematic diagram of turbine optimum timing location determining method one embodiment of the present invention.Preferably, this implementation Example can be performed by turbine optimum timing position determining means of the present invention.This method comprises the following steps：
Step 1, in engineering model as shown in Figure 2, upstream vane tail is obtained by the inverting of steady numerical simulation result Mark is in downstream transport process along journey VELOCITY DISTRIBUTION.
Step 2, determine upstream vane tail in the circumferential total inclined of downstream stator porch along journey VELOCITY DISTRIBUTION according to described Move distance.
In one embodiment of the invention, as shown in figure 3, step 2 can include：
Step 21, very first time t of the upstream vane tail by movable vane passage is determined along journey VELOCITY DISTRIBUTION according to described1。
In one embodiment of the invention, step 21 can include：At least three are in axial direction taken in movable vane passage Point, wherein described at least 3 points include movable vane porch, movable vane exit and at least one movable vane passage intermediate point；From calmly Tail axial velocity value at least three points described in being obtained in constant value analog result；According to tail (upstream at least three point Stator tail) axial velocity value, movable vane axial direction chord length, at least one movable vane passage intermediate point to movable vane leading edge axis away from From determining the very first time t1。
In a specific embodiment of the invention, step 21 can include：In axial direction take in movable vane passage at 3 points, Wherein described 3 points are included at movable vane porch, movable vane exit and the convex back of the body of suction surface；Obtained from steady numerical simulation result Obtain tail axial velocity value at 3 point；And assume tail in movable vane passage interior edge axial velocity profile as shown in figure 4, then The very first time t is determined according to formula 11。
Wherein, L2Represent movable vane axial direction chord length, LmRepresent the convex back of the body of movable vane passage suction surface to movable vane leading edge axis away from From VZ1Represent tail in movable vane porch axial velocity, VZ2Represent tail in movable vane exit axial velocity, VZmRepresent tail Axial velocity at the convex back of the body of movable vane passage suction surface.
In another specific embodiment of the present invention, step 21 can include：In axial direction take in movable vane passage at 5 points, Wherein described 5 points include movable vane porch, movable vane exit, at the convex back of the body of suction surface, at movable vane porch and the convex back of the body of suction surface it Between a points and the convex back of the body of suction surface at b points between movable vane exit；Described five are obtained from steady numerical simulation result Tail axial velocity value at point；And assume tail in movable vane passage interior edge axial velocity profile as shown in figure 5, then according to formula 2 determine the very first time t1。
Wherein, L2Represent movable vane axial direction chord length, Lm、La、LbRepresent respectively at the convex back of the body of movable vane passage suction surface, a points, b points arrive The axial line distance of movable vane leading edge, VZ1Represent tail in movable vane porch axial velocity, VZ2Represent tail in movable vane exit axial direction Speed, VZm、VZa、VZbRepresent tail at the convex back of the body of movable vane passage suction surface, the axial velocity of a points, b points respectively.
Step 22, determine that upstream vane tail passes through movable vane and downstream stator axial direction spacing along journey VELOCITY DISTRIBUTION according to described The second time t2。
In one embodiment of the invention, step 21 can include：In axial direction take in movable vane passage at 2 points, institute Stating at 2 points includes movable vane exit and downstream stator porch；Movable vane exit tail is obtained from steady numerical simulation result Mark axial velocity value VZ2, downstream stator porch tail axial velocity value VZ3；According to tail axial velocity value at 2 point, The axial space D of movable vane and downstream stator2, second time t is determined according to formula 32。
Step 23, according to the very first time t1With second time t2Determine the circumferential total drift distance.
In one embodiment of the invention, as shown in fig. 6, step 23 can include：
Step 231, as shown in Fig. 2 engineering model, according to the axis space D of upstream vane and movable vane1, upstream vane goes out Mouth geometry angle α, the first stage circumferential offset distance S of upstream vane tail is determined according to formula 41, wherein, the first stage Circumferential offset distance S1For the circumferential offset distance of upstream vane tail in-position A (movable vane porch).
S1=D1×tanα (4)
Step 232, according to movable vane botanical origin radial direction radius R1, turbine rotating speed N and the very first time t1, according to formula 5 determine the second stage circumferential offset distance S of upstream vane tail2, wherein, the second stage circumferential offset distance S2To be upper Swim stator tail and pass through movable vane passage, the circumferential offset distance of arrival movable vane trailing edge exit (position B as shown in Figure 2).
Step 233, according to movable vane botanical origin radial direction radius R1, turbine rotating speed N, movable vane outlet geometry angle beta, movable vane and under Swim the axis space D of stator2And second time t2, determine that the phase III of upstream vane tail is circumferentially inclined according to formula 6 Move distance S3, wherein, the phase III circumferential offset distance S3It is (as shown in Figure 2 from movable vane exit for upstream vane tail Position B) reach downstream stator entrance (position C as shown in Figure 2) circumferential offset distance.
Step 234, according to formula 7 by the first stage circumferential offset distance S1, second stage circumferential offset distance S2With Phase III circumferential offset distance S3And be used as the circumferential total drift distance S, wherein what the circumferential total drift distance S referred to Be tail produces to when being transported to downstream stator entrance from upstream vane, circumferential migration distance altogether.
S=S1+S2+S3 (7)
Step 3, turbine optimum timing position Δ L is obtained according to the circumferential total drift distance S.
In one embodiment of the invention, step 3 can include：According to the circumferential total drift distance S, downstream stator Botanical origin radial direction radius R2And downstream stator lobe numbers Num, turbine optimum timing position Δ L is determined according to formula 8, its Described in turbine optimum timing position Δ L be downstream stator and upstream vane circumferential relative position, K is integer.
The turbine optimum timing location determining method provided based on the above embodiment of the present invention, reached by calculating tail Circumferential offset distance obtains the engineering model of optimum timing position during the position of downstream, and specific implementation can pass through permanent numerical computations Inverting tail downstream transport process along journey VELOCITY DISTRIBUTION, and combine the engineering model proposed and calculate tail to enter in downstream stator Circumferential offset distance at mouthful, can obtain turbine optimum timing position, the above embodiment of the present invention avoids by the offset distance Substantial amounts of unsteady numerical simulations process.
The deficiency for aiming to overcome that the permanent design method of existing turbine of the above embodiment of the present invention, is determined with turbine Based on normal design method, the UNSTEADY FLOW phase in each section of adjustment blade, before being finally presented as adjustment blade, trailing edge line Relative tertiary location, so as to provide a kind of engineering model of the unsteady matching of fast and effeciently impeller good luck line, it is possible thereby to Turbine optimum timing position is fast and effeciently searched out during aerodynamic design.
Thus the above embodiment of the present invention has the advantages that：(1) the above embodiment of the present invention can utilize sequential Effect improves turbine isentropic efficiency；(2) the above embodiment of the present invention can search out turbine optimum timing position with simple and fast； (3) the numerical computations elapsed time that the above embodiment of the present invention is used is few, and it is few to take computer resource；(4) the above-mentioned reality of the present invention Unsteady numerical computations can be avoided by applying example, improve designer's operating efficiency.
Fig. 7 is the schematic diagram of turbine optimum timing position determining means one embodiment of the present invention.Impeller shown in Fig. 7 Machine optimum timing position determining means can include speed acquiring module 700, total drift acquisition module 800 and optimum timing and obtain Module 900, wherein：
Speed acquiring module 700, transported for obtaining upstream vane tail by steady numerical simulation result downstream Process along journey VELOCITY DISTRIBUTION.
Total drift acquisition module 800, for determining upstream vane tail in downstream stator along journey VELOCITY DISTRIBUTION according to described The circumferential total drift distance of porch.
Optimum timing acquisition module 900, for obtaining turbine optimum timing position according to the circumferential total drift distance.
In one embodiment of the invention, the optimum timing acquisition module 900 can be used for total according to the circumference Offset distance, downstream stator botanical origin radial direction radius and downstream stator lobe numbers, when determining that turbine is optimal according to formula 8 Tagmeme is put, wherein the turbine optimum timing position is the circumferential relative position of downstream stator and upstream vane.
Fig. 8 is the schematic diagram of total drift acquisition module in one embodiment of the invention.As shown in figure 8, in Fig. 7 embodiments Total drift acquisition module 800 can include time acquisition unit 810 and total drift acquiring unit 820, wherein：
Time acquisition unit 810, for determining that upstream vane tail passes through movable vane passage along journey VELOCITY DISTRIBUTION according to described The very first time and upstream vane tail pass through the second time of movable vane and downstream stator axial direction spacing.
Total drift acquiring unit 820, for determining that the circumference is always inclined according to the very first time and second time Move distance.
Fig. 9 is the schematic diagram of total drift acquisition module in another embodiment of the present invention.Compared with Fig. 8 embodiments, Fig. 9 is implemented In example, the time acquisition unit 810 can include the time acquisition submodule 812 of very first time acquisition submodule 811 and second, Wherein：
Very first time acquisition submodule 811, in axial direction taking in movable vane passage at least 3 points, wherein it is described extremely Few 3 points include movable vane porch, movable vane exit and at least one movable vane passage intermediate point；From steady numerical simulation result Tail axial velocity value at least three points described in middle acquisition；According to tail axial velocity value, movable vane axial direction at least three point The axial line distance of chord length, at least one movable vane passage intermediate point to movable vane leading edge, according to described in formula 1 or the determination of formula 2 The very first time.
Second time acquisition submodule 812, in axial direction taking in movable vane passage at 2 points, described 2 points include moving Leaf exit and downstream stator porch；Tail axial velocity value at described 2 points is obtained from steady numerical simulation result； According to the axial spacing of tail axial velocity value, movable vane and downstream stator at 2 point, when determining described second according to formula 3 Between.
Figure 10 is the schematic diagram of total drift acquiring unit in one embodiment of the invention.As shown in Figure 10, Fig. 8 or Fig. 9 is real Apply total drift acquiring unit 820 in example can include the first skew acquisition submodule 821, second offset acquisition submodule 822, 3rd skew acquisition submodule 823 and total drift acquisition submodule 824, wherein：
First skew acquisition submodule 821, it is several for the axis spacing according to upstream vane and movable vane, upstream vane outlet What angle, the first stage circumferential offset distance of upstream vane tail is determined according to formula 4, wherein, the first stage is circumferentially inclined Move the circumferential offset distance that distance reaches movable vane porch for upstream vane tail.
Second skew acquisition submodule 822, for according to movable vane botanical origin radial direction radius, turbine rotating speed and described the One time, the second stage circumferential offset distance of upstream vane tail is determined according to formula 5, wherein, the second stage is circumferential Offset distance is that upstream vane tail passes through movable vane passage, the circumferential offset distance in arrival movable vane exit.
3rd skew acquisition submodule 823, for being exported according to movable vane botanical origin radial direction radius, turbine rotating speed, movable vane The axis spacing at geometry angle, movable vane and downstream stator and second time, the of upstream vane tail is determined according to formula 6 Three stage circumferential offset distances, wherein, phase III circumferential offset distance for upstream vane tail from movable vane exit to Up to the circumferential offset distance of downstream stator entrance.
Total drift acquisition submodule 824, for according to formula 7 by the first stage circumferential offset distance, second stage Circumferential offset distance and phase III circumferential offset distance and be used as the circumferential total drift distance.
The turbine optimum timing position determining means provided based on the above embodiment of the present invention, reached by calculating tail Circumferential offset distance obtains the engineering model of optimum timing position during the position of downstream, and specific implementation can pass through permanent numerical computations Inverting tail downstream transport process along journey VELOCITY DISTRIBUTION, and combine the engineering model proposed and calculate tail to enter in downstream stator Circumferential offset distance at mouthful, can obtain turbine optimum timing position, the above embodiment of the present invention avoids by the offset distance Substantial amounts of unsteady numerical simulations process.
The deficiency for aiming to overcome that the permanent design method of existing turbine of the above embodiment of the present invention, is determined with turbine Based on normal design method, the UNSTEADY FLOW phase in each section of adjustment blade, before being finally presented as adjustment blade, trailing edge line Relative tertiary location, so as to provide a kind of engineering model of the unsteady matching of fast and effeciently impeller good luck line, it is possible thereby to Turbine optimum timing position is fast and effeciently searched out during aerodynamic design.
Thus the above embodiment of the present invention has the advantages that：(1) the above embodiment of the present invention can utilize sequential Effect improves turbine isentropic efficiency；(2) the above embodiment of the present invention can search out turbine optimum timing position with simple and fast； (3) the numerical computations elapsed time that the above embodiment of the present invention is used is few, and it is few to take computer resource；(4) the above-mentioned reality of the present invention Unsteady numerical computations can be avoided by applying example, improve designer's operating efficiency.
Speed acquiring module 700 described above, total drift acquisition module 800, optimum timing acquisition module 900 etc. Functional unit can be implemented as perform function described herein general processor, programmable logic controller (PLC) (PLC), Digital signal processor (DSP), application specific integrated circuit (ASIC), field programmable gate array (FPGA) or other programmable patrol Volume device, discrete gate either transistor logic, discrete hardware components or it is any appropriately combined.
So far, the present invention is described in detail.In order to avoid the design of the masking present invention, it is public that this area institute is not described Some details known.Those skilled in the art as described above, can be appreciated how to implement technology disclosed herein completely Scheme.
One of ordinary skill in the art will appreciate that hardware can be passed through by realizing all or part of step of above-described embodiment To complete, by program the hardware of correlation can also be instructed to complete, described program can be stored in a kind of computer-readable In storage medium, storage medium mentioned above can be read-only storage, disk or CD etc..
Description of the invention provides for the sake of example and description, and is not exhaustively or by the present invention It is limited to disclosed form.Many modifications and variations are obvious for the ordinary skill in the art.Select and retouch State embodiment and be to more preferably illustrate the principle and practical application of the present invention, and one of ordinary skill in the art is managed The present invention is solved so as to design the various embodiments with various modifications suitable for special-purpose.

Claims (12)

  1. A kind of 1. turbine optimum timing location determining method, it is characterised in that including：
    Upstream vane tail is obtained in downstream transport process along journey VELOCITY DISTRIBUTION by steady numerical simulation result；
    According to described circumferential total drift distance of the upstream vane tail in downstream stator porch is determined along journey VELOCITY DISTRIBUTION；
    Turbine optimum timing position is obtained according to the circumferential total drift distance.
  2. 2. according to the method for claim 1, it is characterised in that described to determine upstream vane along journey VELOCITY DISTRIBUTION according to described Circumferential total drift distance of the tail in downstream stator porch includes：
    According to described the very first time and upstream vane tail of the upstream vane tail by movable vane passage are determined along journey VELOCITY DISTRIBUTION Mark passes through movable vane and the second time of downstream stator axial direction spacing；
    The circumferential total drift distance is determined according to the very first time and second time.
  3. 3. according to the method for claim 2, it is characterised in that described to determine upstream vane along journey VELOCITY DISTRIBUTION according to described Tail is included by the very first time of movable vane passage：
    In axial direction take in movable vane passage at least 3 points, wherein include movable vane porch, movable vane outlet at described at least 3 points Place and at least one movable vane passage intermediate point；
    Obtain tail axial velocity value at described from steady numerical simulation result at least 3 points；
    According to tail axial velocity value, movable vane axial direction chord length, at least one movable vane passage intermediate point at least three point To the axial line distance of movable vane leading edge, the very first time is determined.
  4. 4. according to the method for claim 2, it is characterised in that described to determine upstream vane along journey VELOCITY DISTRIBUTION according to described Tail is included by movable vane and the second time of downstream stator axial direction spacing：
    In axial direction take in movable vane passage at 2 points, include movable vane exit and downstream stator porch at described 2 points；
    Tail axial velocity value at described 2 points is obtained from steady numerical simulation result；
    According to the axial spacing of tail axial velocity value, movable vane and downstream stator at 2 point, second time is determined.
  5. 5. according to the method any one of claim 2-4, it is characterised in that described according to the very first time and described Second time determined that the circumferential total drift distance included：
    The first stage of upstream vane tail is determined according to the axis spacing of upstream vane and movable vane, upstream vane outlet geometry angle Circumferential offset distance, wherein, the first stage circumferential offset distance reaches the circumference of movable vane porch for upstream vane tail Offset distance；
    The second-order of upstream vane tail is determined according to movable vane botanical origin radial direction radius, turbine rotating speed and the very first time Section circumferential offset distance, wherein, the second stage circumferential offset distance is dynamic by movable vane passage, arrival for upstream vane tail The circumferential offset distance in leaf exit；
    According between the axis of movable vane botanical origin radial direction radius, turbine rotating speed, movable vane outlet geometry angle, movable vane and downstream stator Away from and second time determine the phase III circumferential offset distance of upstream vane tail, wherein, the phase III is circumferential Offset distance is the circumferential offset distance that upstream vane tail reaches downstream stator entrance from movable vane exit；
    By the sum of first stage circumferential offset distance, second stage circumferential offset distance and phase III circumferential offset distance As the circumferential total drift distance.
  6. 6. according to the method any one of claim 1-4, it is characterised in that described according to the circumferential total drift distance Obtaining turbine optimum timing position includes：
    Impeller is determined according to circumferential total drift distance, downstream stator botanical origin radial direction radius and the downstream stator lobe numbers Machine optimum timing position, wherein the turbine optimum timing position is the circumferential relative position of downstream stator and upstream vane.
  7. 7. a kind of turbine optimum timing position determining means, it is characterised in that obtain mould including speed acquiring module, total drift Block and optimum timing acquisition module, wherein：
    Speed acquiring module, for obtaining upstream vane tail on the edge of downstream transport process by steady numerical simulation result Journey VELOCITY DISTRIBUTION；
    Total drift acquisition module, for determining upstream vane tail in downstream stator porch along journey VELOCITY DISTRIBUTION according to described Circumferential total drift distance；
    Optimum timing acquisition module, for obtaining turbine optimum timing position according to the circumferential total drift distance.
  8. 8. device according to claim 7, it is characterised in that the total drift acquisition module include time acquisition unit and Total drift acquiring unit, wherein：
    Time acquisition unit, for according to it is described along journey VELOCITY DISTRIBUTION determine upstream vane tail by movable vane passage first when Between and upstream vane tail pass through the second time of movable vane and downstream stator axial direction spacing；
    Total drift acquiring unit, for determining the circumferential total drift distance according to the very first time and second time.
  9. 9. device according to claim 8, it is characterised in that the time acquisition unit includes very first time acquisition submodule Block, wherein：
    Very first time acquisition submodule, in axial direction taking in movable vane passage at least 3 points, wherein described at least 3 points Including movable vane porch, movable vane exit and at least one movable vane passage intermediate point；Obtained from steady numerical simulation result Tail axial velocity value at least three point；According to tail axial velocity value, movable vane axial direction chord length, institute at least three point At least one movable vane passage intermediate point is stated to the axial line distance of movable vane leading edge, determines the very first time.
  10. 10. device according to claim 8, it is characterised in that the time acquisition unit includes the second time acquisition Module, wherein：
    Second time acquisition submodule, in axial direction taking in movable vane passage at 2 points, described 2 points include movable vane and export Place and downstream stator porch；Tail axial velocity value at described 2 points is obtained from steady numerical simulation result；According to institute The axial spacing of tail axial velocity value, movable vane and downstream stator at 2 points is stated, determines second time.
  11. 11. according to the device any one of claim 8-10, it is characterised in that the total drift acquiring unit includes the One skew acquisition submodule, the second skew acquisition submodule, the 3rd skew acquisition submodule and total drift acquisition submodule, its In：
    First skew acquisition submodule, it is true for the axis spacing according to upstream vane and movable vane, upstream vane outlet geometry angle Determine the first stage circumferential offset distance of upstream vane tail, wherein, the first stage circumferential offset distance is upstream vane Tail reaches the circumferential offset distance of movable vane porch；
    Second skew acquisition submodule, for true according to movable vane botanical origin radial direction radius, turbine rotating speed and the very first time Determine the second stage circumferential offset distance of upstream vane tail, wherein, the second stage circumferential offset distance is upstream vane Tail passes through movable vane passage, the circumferential offset distance in arrival movable vane exit；
    3rd skew acquisition submodule, for according to movable vane botanical origin radial direction radius, turbine rotating speed, movable vane export geometry angle, The axis spacing and second time of movable vane and downstream stator determine the phase III circumferential offset of upstream vane tail away from From, wherein, the phase III circumferential offset distance reaches downstream stator entrance for upstream vane tail from movable vane exit Circumferential offset distance；
    Total drift acquisition submodule, for by first stage circumferential offset distance, second stage circumferential offset distance and the Three stage circumferential offset distances and be used as the circumferential total drift distance.
  12. 12. according to the device any one of claim 7-10, it is characterised in that
    The optimum timing acquisition module be used for according to the circumferential total drift distance, downstream stator botanical origin radial direction radius and Downstream stator lobe numbers determine turbine optimum timing position, wherein the turbine optimum timing position be downstream stator and The circumferential relative position of upstream vane.
CN201610692257.7A 2016-08-19 2016-08-19 Turbine optimum timing location determining method and device Pending CN107766598A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610692257.7A CN107766598A (en) 2016-08-19 2016-08-19 Turbine optimum timing location determining method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610692257.7A CN107766598A (en) 2016-08-19 2016-08-19 Turbine optimum timing location determining method and device

Publications (1)

Publication Number Publication Date
CN107766598A true CN107766598A (en) 2018-03-06

Family

ID=61262029

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610692257.7A Pending CN107766598A (en) 2016-08-19 2016-08-19 Turbine optimum timing location determining method and device

Country Status (1)

Country Link
CN (1) CN107766598A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109578085A (en) * 2018-12-26 2019-04-05 中国船舶重工集团公司第七0三研究所 A method of it is tilted by guide vane and weakens the unsteady active force of turbine rotor blade
CN109737069A (en) * 2019-01-31 2019-05-10 浙江理工大学 For studying the adjustable experimental provision of multistage centrifugal pump time

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5486091A (en) * 1994-04-19 1996-01-23 United Technologies Corporation Gas turbine airfoil clocking
CN1955440A (en) * 2005-10-28 2007-05-02 中国科学院工程热物理研究所 Three-D sequential effect maximization method for multi-stage turbomachine
CN103670526A (en) * 2012-09-10 2014-03-26 通用电气公司 Method of clocking a turbine by reshaping the turbine's downstream airfoils

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5486091A (en) * 1994-04-19 1996-01-23 United Technologies Corporation Gas turbine airfoil clocking
CN1955440A (en) * 2005-10-28 2007-05-02 中国科学院工程热物理研究所 Three-D sequential effect maximization method for multi-stage turbomachine
CN103670526A (en) * 2012-09-10 2014-03-26 通用电气公司 Method of clocking a turbine by reshaping the turbine's downstream airfoils

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张宝: "《1+1对转涡轮时序效应数值研究及预测模型》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109578085A (en) * 2018-12-26 2019-04-05 中国船舶重工集团公司第七0三研究所 A method of it is tilted by guide vane and weakens the unsteady active force of turbine rotor blade
CN109578085B (en) * 2018-12-26 2021-06-22 中国船舶重工集团公司第七0三研究所 Method for weakening unsteady acting force of turbine movable blade through guide blade inclination
CN109737069A (en) * 2019-01-31 2019-05-10 浙江理工大学 For studying the adjustable experimental provision of multistage centrifugal pump time
CN109737069B (en) * 2019-01-31 2023-10-20 浙江理工大学 Adjustable experimental device for researching sequential effect of guide vanes of multistage centrifugal pump

Similar Documents

Publication Publication Date Title
Gier et al. Interaction of shroud leakage flow and main flow in a three-stage LP turbine
Gao Advanced numerical simulation of corner separation in a linear compressor cascade
Hoghooghi et al. A new low-cost swirler for axial micro hydro turbines of low head potential
CN107766598A (en) Turbine optimum timing location determining method and device
Gunn et al. Non-axisymmetric stator design for boundary layer ingesting fans
CN105197255B (en) A kind of air intake duct curve design method
Wisler Shock wave and flow velocity measurements in a high speed fan rotor using the laser velocimeter
CN109386381A (en) Flow splitter design method
CN108194412A (en) The uniaxial CO2 compressor latter end Stage of discharge coefficient 0.0056 and method for designing impeller
CN104791172B (en) The mixed flow pump turbine blade space method for determining position of bidirectional shunt
CN110610034A (en) Generation method of hydraulic characteristics of mixed-flow water turbine
Ko¨ nig et al. Clocking Effects in a 1.5-Stage Axial Turbine: Boundary Layer Behaviour at Midspan
Blaszczak Efficiency improvement and noise reduction through stator-stator clocking effect of a two-stage turbine
Xiaoxiong et al. Development of an improved streamline curvature approach for transonic axial compressor
Schneider et al. On the interaction of streamwise vorticity with a rotating turbine blade
Simpson et al. Design, validation and application of a radial cascade for centrifugal compressors
Le Guyader et al. Effects of tip shroud geometries on low-pressure turbine performance
CN108121838A (en) Turbine setting out curbs method and device
Bayeul-Lainé et al. Investigations inside a vaned diffuser of a centrifugal pump at low flowrates.
Moore et al. On three-dimensional flow in centrifugal impellers
Gier et al. Modeling and Analysis of Main Flow–Shroud Leakage Flow Interaction in LP Turbines
CN108561331A (en) 0.0264 light substance high energy head compressor model grade of discharge coefficient and design method
Hall et al. Non-axisymmetric Stator Design for Boundary Layer Ingesting Fans
Lichtfuss Customized profiles—The beginning of an era: a short history of blade design
CN106687684A (en) Method for determining the operating point of a hydraulic machine and installation for converting hydraulic energy

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20180306