CN113107606A - Thermodynamic calculation and design algorithm for transverse stage of steam turbine - Google Patents

Thermodynamic calculation and design algorithm for transverse stage of steam turbine Download PDF

Info

Publication number
CN113107606A
CN113107606A CN202110505471.8A CN202110505471A CN113107606A CN 113107606 A CN113107606 A CN 113107606A CN 202110505471 A CN202110505471 A CN 202110505471A CN 113107606 A CN113107606 A CN 113107606A
Authority
CN
China
Prior art keywords
blade
inlet
velocity
outlet
movable
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.)
Granted
Application number
CN202110505471.8A
Other languages
Chinese (zh)
Other versions
CN113107606B (en
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.)
Harbin Turbine Co Ltd
Hadian Power Equipment National Engineering Research Center Co Ltd
Original Assignee
Harbin Turbine 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 Harbin Turbine Co Ltd filed Critical Harbin Turbine Co Ltd
Priority to CN202110505471.8A priority Critical patent/CN113107606B/en
Publication of CN113107606A publication Critical patent/CN113107606A/en
Application granted granted Critical
Publication of CN113107606B publication Critical patent/CN113107606B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • F01D21/003Arrangements for testing or measuring
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D21/00Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D25/00Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D5/00Blades; Blade-carrying members; Heating, heat-insulating, cooling or antivibration means on the blades or the members
    • F01D5/12Blades
    • F01D5/14Form or construction
    • F01D5/141Shape, i.e. outer, aerodynamic form
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D5/00Blades; Blade-carrying members; Heating, heat-insulating, cooling or antivibration means on the blades or the members
    • F01D5/12Blades
    • F01D5/14Form or construction
    • F01D5/147Construction, i.e. structural features, e.g. of weight-saving hollow blades
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01DNON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
    • F01D9/00Stators
    • F01D9/02Nozzles; Nozzle boxes; Stator blades; Guide conduits, e.g. individual nozzles
    • F01D9/04Nozzles; Nozzle boxes; Stator blades; Guide conduits, e.g. individual nozzles forming ring or sector
    • F01D9/041Nozzles; Nozzle boxes; Stator blades; Guide conduits, e.g. individual nozzles forming ring or sector using blades
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Architecture (AREA)
  • Fluid Mechanics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Control Of Turbines (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)

Abstract

A calculation and design algorithm for the thermodynamic power of a transverse stage of a steam turbine belongs to the technical field of steam turbines. The method aims to solve the problems that a three-dimensional CFD is adopted for design in a traditional design algorithm, the design process is complex, the design accuracy cannot be guaranteed, and the working efficiency is low. The method of iterative iteration by adopting the three-dimensional CFD method has complex process, one scheme needs about 20-30 working hours, and the time is longer if the model selection problem is considered. By adopting the method, one scheme can be controlled within 3 working hours even if the type selection is considered, and the period is shortened by about 90 percent. After a plurality of schemes are needed in consideration of engineering design, the final construction scheme can be determined, and the time is saved remarkably. The method is suitable for designing and calculating the transverse stage of the steam turbine.

Description

Thermodynamic calculation and design algorithm for transverse stage of steam turbine
Technical Field
The invention relates to the technical field of steam turbines, in particular to a thermodynamic calculation and design algorithm for a transverse stage of a steam turbine.
Background
In recent years, with the increase in the capacity of a steam turbine unit, a structure of full-cycle steam admission and a steam supplement valve has appeared. In order to adapt to the whole-cycle steam admission, the arrangement mode of a volute and a transverse stage is required. The transverse stage is different from the traditional axial flow stage and the traditional centripetal stage, no mature one-dimensional design method exists at present, and a three-dimensional CFD method is required. The traditional axial flow stage and centripetal stage have the same parameters of a static blade outlet and a movable blade inlet, while the transverse stage and the centripetal stage have different positions and greatly change thermal parameters. The transverse stage design not only needs to consider one-dimensional thermodynamic calculation, movable and static blade matching and the like, but also influences the through-flow capacity of the unit, and the adoption of a three-dimensional CFD process is complicated and seriously restricts the development of a steam turbine.
In summary, the conventional design algorithm is designed by using three-dimensional CFD, which is a complicated design process, and cannot ensure the design accuracy, and the working efficiency is low.
Disclosure of Invention
The invention provides a thermal calculation and design algorithm for a transverse stage of a steam turbine, aiming at solving the problems that the traditional design algorithm is designed by adopting three-dimensional CFD (computational fluid dynamics), the design process is complex, the design accuracy cannot be ensured, and the working efficiency is low.
The invention relates to a thermodynamic calculation method for a transverse stage of a steam turbine, which comprises the following specific algorithms:
step one, simplifying a volute and a transverse stage into a traditional axial flow stage according to a traditional design idea, carrying out through-flow design, and determining thermal and geometric parameters;
step two, extracting boundary conditions from the result of the step one, and using the boundary conditions as input data of the method of the invention, wherein the input data comprises static pressure p behind the movable blade4Static enthalpy i4A flow rate G; total enthalpy of volute inlet
Figure BDA0003058207030000011
Inlet total pressure
Figure BDA0003058207030000012
Height L of moving blade4Diameter of pitch circle of moving blade D4Angle beta of outlet of moving blade relative to steam flow4(ii) a The rotation speed omega; effective enthalpy drop
Figure BDA0003058207030000013
Step three, according to the static pressure p behind the movable blades4Static enthalpy i4Determining the density rho by means of a water vapor calculation program4
Step four, calculating the linear velocity of the outlet of the movable blade
Figure BDA0003058207030000014
According to the conservation of mass, the axial velocity C behind the movable blade is obtained4z
Step five, solving the tangential velocity C of the outlet of the movable blade according to the velocity triangle4u
Sixthly, the outlet relative speed of the movable blade
Figure BDA0003058207030000015
Step seven, setting a movable blade speed coefficient psi;
step eight, assuming the axial speed C of the inlet of the movable blade3Z
Nine steps, inlet tangential velocity of movable vane
Figure BDA0003058207030000021
Step ten, moving blade inlet absolute speed C3
Eleven, according to the conservation of energy, the total enthalpy of the inlet of the movable vane
Figure BDA0003058207030000022
Total enthalpy of inlet and volute
Figure BDA0003058207030000023
Equality, calculating static enthalpy of inlet of moving blade
Figure BDA0003058207030000024
Is prepared from (i)3,ρ3) Investigating the water vapour table3
Step twelve, inlet density of moving blade
Figure BDA0003058207030000025
Is prepared from (i)3,ρ3) Obtaining the inlet entropy s of the movable vane by adopting a water vapor calculation program3(ii) a Is prepared from (S)3,p4) Adopting a water vapor calculation program to obtain the isentropic enthalpy drop i of the movable blades4s(ii) a By Δ h4=i4-i4sTo obtain the isentropic relative velocity of the outlet of the movable vane
Figure BDA0003058207030000026
Relative velocity of bucket outlet
Figure BDA0003058207030000027
Step thirteen, comparing W calculated in step twelve4W with step six4If the relative error is less than 10-4Continuing the next step; otherwise, revision C3z, repeating the steps nine to twelve until stepW of step twelve and step six4Relative error less than 10-4(ii) a Through the steps, all thermodynamic parameters, speed components and geometric parameters before and after the movable blade can be obtained;
fourteen, setting a speed coefficient epsilon from a section 2-2 to a section 3-3 of the computing station; according to the law of circular quantity, the tangential velocity of the outlet of the stator blade 3
Figure BDA0003058207030000028
Fifthly, obtaining all thermal parameters, velocity components and geometric parameters of the outlet of the stationary blade 3 by imitating the processes from the seventh step to the thirteenth step;
sixthly, obtaining the tangential speed of the volute inlet according to the thermal parameters of the volute inlet in the step one
Figure BDA0003058207030000029
Then, the total energy loss coefficient xi and the volute speed coefficient from the 0-0 section to the 4-4 section of the calculation station are obtained
Figure BDA00030582070300000210
According to the law of circular quantity, stator blade inlet tangential velocity
Figure BDA00030582070300000211
Seventhly, obtaining all thermal parameters, speed components and geometric parameters of the inlet of the static blade 3 by imitating the process from the step seven to the step thirteen;
further, C calculated in the fourth step4zIs calculated by the formula C4z=G/(C4z*π*D4*L44);
Further, in the fifth step, according to the velocity triangle, the formula for obtaining the tangential velocity of the outlet of the movable blade is C4u=C4z*tan(β4)-U4
Further, the absolute speed C of the inlet of the movable blade is calculated in the step ten3Is of the formula
Figure BDA0003058207030000031
The invention relates to a design algorithm of a turbine transverse stage, which comprises the following specific steps:
step one, transversely placing a fixed blade inlet steam flow angle
Figure BDA0003058207030000032
Step two, according to a1Selecting proper transverse stator blade molded line to ensure the geometric angle of molded line inlet and ainRespectively meet the design standard;
and step three, determining the relative grid distance T/b corresponding to the highest efficiency point by combining the molded line loss library according to the selected stator blade molded line. Note that T here is the pitch corresponding to the pitch circle diameter;
step four, transversely arranging the steam flow angle of the stationary blade outlet
Figure BDA0003058207030000033
Note that t here is the pitch of the steam outlet side pitch circle diameter pair;
step five, according to aoutAnd determining the installation angle gamma of the transverse static blade so as to complete transverse static blade selection.
Compared with the prior art, the invention has the following beneficial effects:
the invention overcomes the defects of the prior art and thoroughly solves the problem that the transverse stage of the steam turbine is lack of an effective one-dimensional program. In the traditional method, overall parameters such as transverse enthalpy drop and loss are distributed according to axial flow calculation in one dimension, and a three-dimensional CFD method is adopted to iterate repeatedly to obtain a result meeting the requirement, even the problem that the one-dimensional requirement cannot be met and the one-dimensional calculation needs to be modified again occurs;
the design period is greatly shortened, and the working efficiency is improved. The method of iterative iteration by adopting the three-dimensional CFD method has complex process, one scheme needs about 20-30 working hours, and the time is longer if the model selection problem is considered. By adopting the method, one scheme can be controlled within 3 working hours even if the type selection is considered, and the period is shortened by about 90 percent. After a plurality of schemes are needed in consideration of engineering design, the final construction scheme can be determined, and the time is saved remarkably;
the precision completely meets the design requirements. The calculation method of the invention is based on the basic principles of fluid mechanics and impeller machinery, and corrects the speed coefficient (namely loss) which only affects the precision based on a large amount of three-dimensional CFD results, and the error between the speed coefficient and the three-dimensional CFD results can be controlled within 3%.
The invention carries out program embedding aiming at the particularity of the selection of the transverse stationary blade, and can realize the combined design of one-dimensional calculation and selection. The invention is suitable for all the horizontal stage structures of the steam turbine with a plurality of inlets of the volute, double split flow of the movable blades and the like, and has wide adaptability.
Drawings
FIG. 1 is a schematic diagram of a prior art horizontal stage computing station;
FIG. 2 is a schematic diagram of a prior art axial flow computing station;
FIG. 3 is a schematic diagram of a prior art centrifugal stage computing station;
FIG. 4 is a schematic diagram of a turbine cross-stage thermodynamic computing plant according to the present invention;
FIG. 5 is a cross-stage vane sizing graph of a turbine cross-stage design algorithm in accordance with the present invention;
FIG. 6 is a volute dual-inlet configuration of the present invention;
FIG. 7 is a cross-level dual-split architecture of the present invention;
the figure includes a transverse vane 1, a blade 2, and a vane 3.
Detailed Description
The first embodiment is as follows: the present embodiment is described with reference to fig. 1 to 7, and a specific algorithm of the thermodynamic calculation of the turbine transverse stage according to the present embodiment is as follows:
step one, simplifying a volute and a transverse stage into a traditional axial flow stage according to a traditional design idea, carrying out through-flow design, and determining thermal and geometric parameters;
step two, extracting boundary conditions from the result of the step one, and using the boundary conditions as input data of the method of the invention, wherein the input data comprises static pressure p behind the movable blade4Static enthalpy i4A flow rate G; total enthalpy of volute inlet
Figure BDA0003058207030000041
Inlet total pressure
Figure BDA0003058207030000042
Height L of moving blade4Diameter of pitch circle of moving blade D4Angle beta of outlet of moving blade relative to steam flow4(ii) a The rotation speed omega; effective enthalpy drop
Figure BDA0003058207030000043
Step three, according to the static pressure p behind the movable blades4Static enthalpy i4Determining the density rho by means of a water vapor calculation program4
Step four, calculating the linear velocity of the outlet of the movable blade
Figure BDA0003058207030000044
According to the conservation of mass, the axial velocity C behind the movable blade is obtained4z
Step five, solving the tangential velocity C of the outlet of the movable blade according to the velocity triangle4u
Sixthly, the outlet relative speed of the movable blade
Figure BDA0003058207030000045
Step seven, setting a movable blade speed coefficient psi;
step eight, assuming the axial speed C of the inlet of the movable blade3Z
Nine steps, inlet tangential velocity of movable vane
Figure BDA0003058207030000046
Step ten, moving blade inlet absolute speed C3
Eleven, according to the conservation of energy, the total enthalpy of the inlet of the movable vane
Figure BDA0003058207030000047
Total enthalpy of inlet and volute
Figure BDA0003058207030000048
Equality, calculating static enthalpy of inlet of moving blade
Figure BDA0003058207030000049
Is prepared from (i)3,ρ3) Investigating the water vapour table3
Step twelve, inlet density of moving blade
Figure BDA0003058207030000051
Is prepared from (i)3,ρ3) Obtaining the inlet entropy s of the movable vane by adopting a water vapor calculation program3(ii) a Is prepared from (S)3,p4) Adopting a water vapor calculation program to obtain the isentropic enthalpy drop i of the movable blades4s(ii) a By Δ h4=i4-i4sTo obtain the isentropic relative velocity of the outlet of the movable vane
Figure BDA0003058207030000052
Relative velocity of bucket outlet
Figure BDA0003058207030000053
Step thirteen, comparing W calculated in step twelve4W with step six4If the relative error is less than 10-4Continuing the next step; otherwise, revision C3z, repeating the steps nine to twelve until the W of the step twelve and the step six4Relative error less than 10-4(ii) a Through the steps, all thermodynamic parameters, speed components and geometric parameters before and after the movable blade can be obtained;
fourteen, setting a speed coefficient epsilon from a section 2-2 to a section 3-3 of the computing station; according to the law of circular quantity, the tangential velocity of the outlet of the stator blade 3
Figure BDA0003058207030000054
Fifthly, obtaining all thermal parameters, velocity components and geometric parameters of the outlet of the stationary blade 3 by imitating the processes from the seventh step to the thirteenth step;
sixthly, obtaining the tangential speed of the volute inlet according to the thermal parameters of the volute inlet in the step one
Figure BDA0003058207030000055
Then, the total energy loss coefficient xi and the volute speed coefficient from the 0-0 section to the 4-4 section of the calculation station are obtained
Figure BDA0003058207030000056
According to the law of circular quantity, stator blade inlet tangential velocity
Figure BDA0003058207030000057
Seventhly, obtaining all thermal parameters, speed components and geometric parameters of the inlet of the static blade 3 by imitating the process from the step seven to the step thirteen.
The second embodiment is as follows: the present embodiment will be described with reference to fig. 1, and the present embodiment is a further limitation of the thermodynamic calculation described in the first embodiment, the thermodynamic calculation of the turbine horizontal stage in the present embodiment, and C calculated in the fourth step4zIs calculated by the formula C4z=G/(C4z*π*D4*L44)。
The third concrete implementation mode: the present embodiment will be described with reference to fig. 1, which is a further limitation of the thermodynamic calculation described in the first embodiment, and in the thermodynamic calculation of a turbine horizontal stage described in the present embodiment, the equation for determining the tangential velocity of the outlet of the rotor blade from the velocity triangle in the fifth step is C4u=C4z*tan(β4)-U4
The fourth concrete implementation mode: the present embodiment will be described with reference to fig. 1, and the present embodiment is a further limitation of the thermodynamic calculation described in the first embodiment, and in the thermodynamic calculation of the turbine horizontal stage described in the present embodiment, the absolute speed C of the inlet of the rotor blade is calculated in the tenth step3Is of the formula
Figure BDA0003058207030000061
The fifth concrete implementation mode: the present embodiment is described with reference to fig. 1 to 7, and the design algorithm of the turbine transverse stage according to the present embodiment specifically includes the following steps:
step one, transversely placing a fixed blade inlet steam flow angle
Figure BDA0003058207030000062
Step two, according to a1Selecting proper transverse stator blade molded line to ensure the geometric angle of molded line inlet and ainRespectively meet the design standard;
and step three, determining the relative grid distance T/b corresponding to the highest efficiency point by combining the molded line loss library according to the selected stator blade molded line. Note that T here is the pitch corresponding to the pitch circle diameter;
step four, transversely arranging the steam flow angle of the stationary blade outlet
Figure BDA0003058207030000063
Note that t here is the pitch of the steam outlet side pitch circle diameter pair;
step five, according to aoutAnd determining the installation angle gamma of the transverse static blade so as to complete transverse static blade selection.
The sixth specific implementation mode: the present embodiment is described with reference to fig. 1 to 7, and a specific algorithm of the thermodynamic calculation of the turbine transverse stage according to the present embodiment is as follows: taking an N660-31/600/620/620 type ultra-supercritical double-reheat extraction condensing steam turbine as an example, the method of the invention is used for calculating the parameters of a high-pressure single-flow transverse stage, and the calculation results of typical parameters are shown in the following table:
parameter name Unit of Three dimensional CFD results The invention calculates the result Error from three-dimensional CFD
Flow rate t/h 1489.92 1501.1 +0.75%
Static pressure behind the moving blade MPa 8.51 8.51 0.0%
Static pressure in front of moving blade MPa 8.76 8.85 +1.03%
Static pressure behind stationary blade MPa 8.87 8.78 -1.01%
Static pressure before stationary blade MPa 9.09 9.14 +0.56%
Rear velocity of moving blade m/s 60.21 60.5 +0.48%
Forward speed of moving blade m/s 173.39 175.33 +1.12%
Stator blade rear velocity m/s 144.99 143.12 -1.29%
Stator blade front velocity m/s 48.87 49.03 +0.32%
The seventh embodiment: the present embodiment is described with reference to fig. 1 to 7, and a specific algorithm of the thermodynamic calculation of the turbine transverse stage according to the present embodiment is as follows: taking an N660-31/600/620/620 type ultra-supercritical double-reheat extraction condensing steam turbine as an example, the method of the invention is used for calculating the parameters of the low-pressure double-flow-direction transverse stage, and the calculation results of typical parameters are shown in the following table:
parameter name Unit of Three dimensional CFD results The invention calculates the result Error from three-dimensional CFD
Flow rate t/h 601.44 598.12 -0.55%
Static pressure behind the moving blade MPa 0.53 0.54 +1.88%
Static pressure in front of moving blade MPa 0.64 0.63 -1.56%
Static pressure behind stationary blade MPa 0.68 0.66 -2.94%
Static pressure before stationary blade MPa 0.80 0.82 +2.50%
Rear velocity of moving blade m/s 81.50 82.47 +1.19%
Forward speed of moving blade m/s 359.58 371.29 +3.26%
Stator blade rear velocity m/s 324.20 314.33 -3.04%
Stator blade front velocity m/s 131.80 134.37 +1.95%
The physical meaning of the variables in the invention is as follows:
variables are as follows: subscripts:
p: static pressure 0: calculate value 0
i: static enthalpy 1: computing station 1
G: and (3) flow rate 2: computing station 2
L: leaf height 3: computing station 3
D: pitch circle diameter 4: computing station 4
Beta: relative steam flow angle of the outlet of the movable blade 5: computing station 5
Omega: rotating speed s: isentropic
z: axial component
Hu: effective enthalpy drop
ρ: density u: tangential component
U: linear velocity r: radial component
Psi: moving blade velocity coefficient in: stationary blade steam inlet side
C: resultant velocity out: stationary blade steam outlet side
s: entropy of the entropy
Epsilon: speed coefficient of computing station 2-computing station 3
W: relative velocity
Δ h: static-isentropic enthalpy
r: radius of
O: size of throat
t: pitch of
T: pitch of
γ: mounting angle
Xi: calculating total energy loss coefficients for stations 0-0 to 4-4
μ: a volute velocity coefficient.

Claims (5)

1. A thermal calculation of a turbine transverse stage is characterized in that: the specific algorithm is as follows:
step one, simplifying a volute and a transverse stage into a traditional axial flow stage according to a traditional design idea, carrying out through-flow design, and determining thermal and geometric parameters;
step two, extracting boundary conditions from the result of the step one to be used as the basic bookInput data for the inventive method, including static pressure p behind the rotor blade4Static enthalpy i4A flow rate G; total enthalpy of volute inlet
Figure FDA0003058207020000011
Inlet total pressure
Figure FDA0003058207020000012
Height L of moving blade4Diameter of pitch circle of moving blade D4Angle beta of outlet of moving blade relative to steam flow4(ii) a The rotation speed omega; effective enthalpy drop
Figure FDA0003058207020000013
Step three, according to the static pressure p behind the movable blades4Static enthalpy i4Determining the density rho by means of a water vapor calculation program4
Step four, calculating the linear velocity of the outlet of the movable blade
Figure FDA0003058207020000014
According to the conservation of mass, the axial velocity C behind the movable blade is obtained4z
Step five, solving the tangential velocity C of the outlet of the movable blade according to the velocity triangle4u
Sixthly, the outlet relative speed of the movable blade
Figure FDA0003058207020000015
Step seven, setting a movable blade speed coefficient psi;
step eight, assuming the axial speed C of the inlet of the movable blade3Z
Nine steps, inlet tangential velocity of movable vane
Figure FDA0003058207020000016
Step ten, moving blade inlet absolute speed C3
Eleven, according to the conservation of energy, the total enthalpy of the inlet of the movable vane
Figure FDA0003058207020000017
Total enthalpy of inlet and volute
Figure FDA0003058207020000018
Equality, calculating static enthalpy of inlet of moving blade
Figure FDA0003058207020000019
Is prepared from (i)3,ρ3) Investigating the water vapour table3
Step twelve, inlet density of moving blade
Figure FDA00030582070200000110
Is prepared from (i)3,ρ3) Obtaining the inlet entropy s of the movable vane by adopting a water vapor calculation program3(ii) a Is prepared from (S)3,p4) Adopting a water vapor calculation program to obtain the isentropic enthalpy drop i of the movable blades4s(ii) a By Δ h4=i4-i4sTo obtain the isentropic relative velocity of the outlet of the movable vane
Figure FDA00030582070200000111
Relative velocity of bucket outlet
Figure FDA00030582070200000112
Step thirteen, comparing W calculated in step twelve4W with step six4If the relative error is less than 10-4Continuing the next step; otherwise, revision C3z, repeating the steps nine to twelve until the W of the step twelve and the step six4Relative error less than 10-4(ii) a Through the steps, all thermodynamic parameters, speed components and geometric parameters before and after the movable blade can be obtained;
fourteen, setting a speed coefficient epsilon from a section 2-2 to a section 3-3 of the computing station; according to the law of circular quantity, the outlet tangential velocity of the stator blade (3)
Figure FDA0003058207020000021
Fifthly, obtaining all thermal parameters, velocity components and geometric parameters of the outlet of the stationary blade (3) by imitating the processes from the seventh step to the thirteenth step;
sixthly, obtaining the tangential speed of the volute inlet according to the thermal parameters of the volute inlet in the step one
Figure FDA0003058207020000022
Then, the total energy loss coefficient xi and the volute speed coefficient from the 0-0 section to the 4-4 section of the calculation station are obtained
Figure FDA0003058207020000023
According to the law of circular quantity, stator blade inlet tangential velocity
Figure FDA0003058207020000024
Seventhly, obtaining all thermodynamic parameters, velocity components and geometric parameters of the inlet of the static blade (3) by imitating the process from the step seven to the step thirteen.
2. A turbine transverse stage thermodynamic calculation as claimed in claim 1 in which: c calculated in the step four4zIs calculated by the formula C4z=G/(C4z*π*D4*L44)。
3. A turbine transverse stage thermodynamic calculation as claimed in claim 1 in which: in the fifth step, according to the velocity triangle, the formula for obtaining the tangential velocity of the outlet of the movable blade is C4u=C4z*tan(β4)-U4
4. A turbine transverse stage thermodynamic calculation as claimed in claim 1 in which: in the step ten, the absolute speed C of the inlet of the movable blade is calculated3Is of the formula
Figure FDA0003058207020000025
5. A design algorithm for a turbine transverse stage is characterized in that: the specific algorithm is as follows:
step one, transversely placing a fixed blade inlet steam flow angle
Figure FDA0003058207020000026
Step two, according to a1Selecting proper transverse stator blade molded line to ensure the geometric angle of molded line inlet and ainRespectively meet the design standard;
and step three, determining the relative grid distance T/b corresponding to the highest efficiency point by combining the molded line loss library according to the selected stator blade molded line. Note that T here is the pitch corresponding to the pitch circle diameter;
step four, transversely arranging the steam flow angle of the stationary blade outlet
Figure FDA0003058207020000027
Note that t here is the pitch of the steam outlet side pitch circle diameter pair;
step five, according to aoutAnd determining the installation angle gamma of the transverse static blade so as to complete transverse static blade selection.
CN202110505471.8A 2021-05-10 2021-05-10 Thermodynamic calculation and design algorithm for transverse stage of steam turbine Active CN113107606B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110505471.8A CN113107606B (en) 2021-05-10 2021-05-10 Thermodynamic calculation and design algorithm for transverse stage of steam turbine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110505471.8A CN113107606B (en) 2021-05-10 2021-05-10 Thermodynamic calculation and design algorithm for transverse stage of steam turbine

Publications (2)

Publication Number Publication Date
CN113107606A true CN113107606A (en) 2021-07-13
CN113107606B CN113107606B (en) 2023-03-24

Family

ID=76721408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110505471.8A Active CN113107606B (en) 2021-05-10 2021-05-10 Thermodynamic calculation and design algorithm for transverse stage of steam turbine

Country Status (1)

Country Link
CN (1) CN113107606B (en)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19755982A1 (en) * 1997-12-17 1999-07-08 Abb Patent Gmbh Steam turbine with wet steam regime
JP2002349203A (en) * 2001-05-24 2002-12-04 National Institute For Materials Science Virtual turbine calculation method
CN1526915A (en) * 2003-02-05 2004-09-08 西门子公司 Steam turbine and method for running steam turbine
CN101050710A (en) * 2006-04-07 2007-10-10 孙敏超 Mixed flow type turbine vane
JP2010265826A (en) * 2009-05-15 2010-11-25 Toshiba Corp Nozzle box for steam turbine and steam turbine
CN102608914A (en) * 2011-12-22 2012-07-25 西安交通大学 Optimization design method of radial-flow-type hydraulic turbine
US20130195641A1 (en) * 2010-07-14 2013-08-01 Isis Innovation Ltd Vane assembly for an axial flow turbine
CN103541774A (en) * 2013-11-14 2014-01-29 上海汽轮机厂有限公司 Method for designing turbine blades
CN105332952A (en) * 2015-11-02 2016-02-17 南京航空航天大学 Small-bend adjustable stator design method
CN106227967A (en) * 2016-08-01 2016-12-14 杭州汽轮机股份有限公司 Industrial steam turbine low-pressure stage group vane type line optimization method
CN109026206A (en) * 2018-07-19 2018-12-18 哈尔滨汽轮机厂有限责任公司 It is a kind of with spiral case into the integral high pressure nuclear steam turbine of vapour formula high-pressure inner cylinder
US20190136704A1 (en) * 2016-05-10 2019-05-09 TURBODEN S. p. A. Mixed flow optimized turbine
CN109812300A (en) * 2019-02-12 2019-05-28 哈尔滨汽轮机厂有限责任公司 It is a kind of consider pretwist amount the small enthalpy drop blade dimensions of steam turbine determine method
CN110863861A (en) * 2019-11-21 2020-03-06 哈尔滨汽轮机厂有限责任公司 Multi-stage small enthalpy drop blade pneumatic design method for improving through-flow efficiency of steam turbine
CN112412546A (en) * 2020-11-23 2021-02-26 东方电气集团东方汽轮机有限公司 Nozzle chamber of industrial steam turbine without median plane
CN112483202A (en) * 2020-12-14 2021-03-12 哈尔滨汽轮机厂有限责任公司 360-degree volute tangential double-side steam inlet high-medium pressure inner cylinder for steam turbine

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19755982A1 (en) * 1997-12-17 1999-07-08 Abb Patent Gmbh Steam turbine with wet steam regime
JP2002349203A (en) * 2001-05-24 2002-12-04 National Institute For Materials Science Virtual turbine calculation method
US20040184908A1 (en) * 2002-02-05 2004-09-23 Detlef Haje Steam turbine and method for operating a steam turbine
CN1526915A (en) * 2003-02-05 2004-09-08 西门子公司 Steam turbine and method for running steam turbine
CN101050710A (en) * 2006-04-07 2007-10-10 孙敏超 Mixed flow type turbine vane
JP2010265826A (en) * 2009-05-15 2010-11-25 Toshiba Corp Nozzle box for steam turbine and steam turbine
US20130195641A1 (en) * 2010-07-14 2013-08-01 Isis Innovation Ltd Vane assembly for an axial flow turbine
CN102608914A (en) * 2011-12-22 2012-07-25 西安交通大学 Optimization design method of radial-flow-type hydraulic turbine
CN103541774A (en) * 2013-11-14 2014-01-29 上海汽轮机厂有限公司 Method for designing turbine blades
CN105332952A (en) * 2015-11-02 2016-02-17 南京航空航天大学 Small-bend adjustable stator design method
US20190136704A1 (en) * 2016-05-10 2019-05-09 TURBODEN S. p. A. Mixed flow optimized turbine
CN106227967A (en) * 2016-08-01 2016-12-14 杭州汽轮机股份有限公司 Industrial steam turbine low-pressure stage group vane type line optimization method
CN109026206A (en) * 2018-07-19 2018-12-18 哈尔滨汽轮机厂有限责任公司 It is a kind of with spiral case into the integral high pressure nuclear steam turbine of vapour formula high-pressure inner cylinder
CN109812300A (en) * 2019-02-12 2019-05-28 哈尔滨汽轮机厂有限责任公司 It is a kind of consider pretwist amount the small enthalpy drop blade dimensions of steam turbine determine method
CN110863861A (en) * 2019-11-21 2020-03-06 哈尔滨汽轮机厂有限责任公司 Multi-stage small enthalpy drop blade pneumatic design method for improving through-flow efficiency of steam turbine
CN112412546A (en) * 2020-11-23 2021-02-26 东方电气集团东方汽轮机有限公司 Nozzle chamber of industrial steam turbine without median plane
CN112483202A (en) * 2020-12-14 2021-03-12 哈尔滨汽轮机厂有限责任公司 360-degree volute tangential double-side steam inlet high-medium pressure inner cylinder for steam turbine

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
A. ROMAGNOLI ETC: "Performance prediction of a nozzled and nozzleless mixed-flow turbine in steady conditions", 《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》 *
A. ROMAGNOLI ETC: "Performance prediction of a nozzled and nozzleless mixed-flow turbine in steady conditions", 《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》, vol. 53, no. 8, 31 August 2011 (2011-08-31), pages 557 - 574, XP028096871, DOI: 10.1016/j.ijmecsci.2011.05.003 *
刘云锋: "某大型涡扇发动机低压涡轮气动设计与分析", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *
刘云锋: "某大型涡扇发动机低压涡轮气动设计与分析", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》, 15 February 2009 (2009-02-15) *
奚忠等: "小型向心透平一维和三维设计分析", 《航空动力学报》, no. 07 *
李文福: "跨音速向心汽轮机气动设计分析", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》, 15 October 2007 (2007-10-15) *
李斌: "先进涡轮气动设计规律的研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *
李斌: "先进涡轮气动设计规律的研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》, 15 February 2012 (2012-02-15) *
潘春雨等: "汽轮机横置静叶级涡轮静叶造型技术研究", 《汽轮机技术》 *
潘春雨等: "汽轮机横置静叶级涡轮静叶造型技术研究", 《汽轮机技术》, vol. 61, no. 01, 28 February 2019 (2019-02-28), pages 26 - 28 *

Also Published As

Publication number Publication date
CN113107606B (en) 2023-03-24

Similar Documents

Publication Publication Date Title
CN111914362B (en) Self-adaptive method for turbofan engine model in research and development stage
CN110929357A (en) Pneumatic design method for high-performance ship gas turbine compressor
CN109598081B (en) Radial-flow turbine pneumatic optimization method based on data dimension reduction and multi-two-dimensional flow surface
Wang et al. Effect of the tip clearance variation on the performance of a centrifugal compressor with considering impeller deformation
CN112685968B (en) Axial flow compressor pneumatic design method based on space load customization thought
CN111027148B (en) Automatic calibration and industrial axial flow compressor performance calculation method for loss lag angle model
CN115186398B (en) Method for determining key angle parameters of inlet guide vane modeling of axial flow compressor
CN114444331A (en) Stage characteristic matching method of multistage axial flow compressor
Aminjan et al. Design and Simulation of Radial Flow Turbine Impeller and Investigation Thermodynamic Properties of Flow in LE and TE.
Goryachkin et al. Three-stage low pressure compressor modernization by means of optimization methods
Popov et al. Multicriteria optimization of axial low pressure compressor of gas turbine power plant
CN109670244B (en) Area adjustment method for guide device of gas turbine after repairing turboshaft engine
CN113107606B (en) Thermodynamic calculation and design algorithm for transverse stage of steam turbine
CN101158991A (en) Pneumatic power design method for big and small blade compressor
Rusanov et al. The influence of stator-rotor interspace overlap of meridional contours on the efficiency of high-pressure steam turbine stages
Lei et al. A rapid and automatic optimal design method for six-stage axial-flow industry compressor
CN116595874A (en) Impeller mechanical performance prediction model parameter optimization method and device and storage medium
Knierim et al. Design process for an advanced reactor coolant pump for a 1400 MW nuclear power plant
Carcasci et al. Heavy duty gas turbine simulation: a compressor IGV airfoil off-design characterization
Sammak et al. Reduced-order through-flow design code for highly loaded, cooled axial turbines
CN112761969A (en) Optimization method for structural parameters of centrifugal pump
Senoo et al. Development and verification tests of titanium 50-inch and 60-inch last stage blades for steam turbines
Krivosheev et al. Optimization of the Stage Number and Parameter Distribution in the Flow Passage at GTE Compressor and Turbine Design
Ziegler et al. Development of a Novel Axial Compressor Generation for Industrial Applications: Part 1—Compressor Design and Performance
Xia et al. Improving Gas Turbine Performance Through Reassembling Degraded Components: An Experimental and Computational Study

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Liu Yunfeng

Inventor after: Ye Dongping

Inventor after: Ma Tianyin

Inventor after: Zhao Hongyu

Inventor after: Liu Zhiqiang

Inventor after: Li Yufeng

Inventor after: Wang Jian

Inventor after: Liu Changchun

Inventor after: Guan Jiwei

Inventor after: Pan Chunyu

Inventor after: Ma Yiliang

Inventor after: Guan Chun

Inventor before: Liu Yunfeng

Inventor before: Ma Tianyin

Inventor before: Zhao Hongyu

Inventor before: Li Yufeng

Inventor before: Wang Jian

Inventor before: Liu Changchun

Inventor before: Guan Jiwei

Inventor before: Pan Chunyu

Inventor before: Ma Yiliang

Inventor before: Guan Chun

Inventor before: Ye Dongping

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230117

Address after: 150000 building 3, high tech production base, Nangang District, Harbin City, Heilongjiang Province

Applicant after: HARBIN TURBINE Co.,Ltd.

Applicant after: HADIAN POWER EQUIPMENT NATIONAL ENGINEERING RESEARCH CENTER CO.,LTD.

Address before: 150046 No. three power road 345, Xiangfang District, Heilongjiang, Harbin

Applicant before: HARBIN TURBINE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant