CN109359339A - Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct - Google Patents
Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct Download PDFInfo
- Publication number
- CN109359339A CN109359339A CN201811081735.6A CN201811081735A CN109359339A CN 109359339 A CN109359339 A CN 109359339A CN 201811081735 A CN201811081735 A CN 201811081735A CN 109359339 A CN109359339 A CN 109359339A
- Authority
- CN
- China
- Prior art keywords
- rotor
- generator
- temperature
- loss
- indicate
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000004364 calculation method Methods 0.000 title claims abstract description 34
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 80
- 229910052802 copper Inorganic materials 0.000 claims abstract description 65
- 239000010949 copper Substances 0.000 claims abstract description 65
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims abstract description 64
- 229910052742 iron Inorganic materials 0.000 claims abstract description 39
- 238000009413 insulation Methods 0.000 claims description 34
- 238000004804 winding Methods 0.000 claims description 27
- 239000012530 fluid Substances 0.000 claims description 25
- 238000012546 transfer Methods 0.000 claims description 20
- 239000013598 vector Substances 0.000 claims description 12
- 208000011580 syndromic disease Diseases 0.000 claims description 7
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 6
- 239000001257 hydrogen Substances 0.000 claims description 6
- 229910052739 hydrogen Inorganic materials 0.000 claims description 6
- 230000005672 electromagnetic field Effects 0.000 claims description 4
- 101710198693 Invasin Proteins 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000004134 energy conservation Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 claims description 2
- 230000003993 interaction Effects 0.000 claims description 2
- 238000010248 power generation Methods 0.000 claims description 2
- 238000009423 ventilation Methods 0.000 abstract description 16
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 208000008312 Tooth Loss Diseases 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000112 cooling gas Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
Abstract
The nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct that the present invention provides a kind of, belongs to steam turbine generator performance detection technical field.Establish Generator Rotor Temperature field computation model；Calculate the copper loss and iron loss of generator amature；The temperature field of the generator amature operation is determined using finite volume method, obtains the draught temperature of mean temperature and air duct in conjunction with the Generator Rotor Temperature field computation model according to the copper loss and iron loss；The copper loss of generator amature is recalculated according to the mean temperature, and is iterated calculating using finite volume method in conjunction with the draught temperature, until the difference of the temperature value of mean temperature and calculating copper loss is less than given threshold.The present invention greatly improves the computational accuracy of the three-dimensional temperature field of the rotor of steam turbo generator with axial ventilation ditch, is more met the thermo parameters method of actual motion.
Description
Technical field
The present invention relates to rotor of steam turbo generator performance detection technical fields, and in particular to a kind of with axial ventilation ditch
Large nuclear-power rotor of steam turbo generator thermo parameters method calculation method.
Background technique
The rotor ventilation system of large turbo-type generator mainly has two class of ventilation from sub-slot, gap-pickup and axial ventilation, and
Ventilation from sub-slot and gap-pickup ventilating patterns have the characteristics that structure is complicated but rotor temperature rise is uniform, and axial ventilation pattern has
The feature that structure is simple but thermal unbalance coefficient is high.Since temperature distribution is non-uniform for the rotor windings with axial ventilation pattern,
Hot(test)-spot temperature and mean temperature are affected by copper loss calculated value.
Steam turbine generator is one of main equipment of power grid, is the immediate producer of electric energy, about 80% electric energy is by steamer
Generator generates.And the up to more than one hundred million members of the cost of nuclear power steam turbine generator, safe power generation are also directly related to a very big section
The stabilization electricity consumption in domain, thus it is guaranteed that the normal operation of steam turbine generator, the direct or indirect health hair for being related to national economy
Exhibition.
If ignoring the influence of the temperature value of winding copper loss, under the influence of biggish thermal unbalance coefficient, it is difficult to quasi-
Really calculate the Temperature Distribution of rotor.Since the rotor hot(test)-spot temperature of large turbo-type generator is difficult to survey, rotor is accurately calculated
Temperature field is one of the method for assessing engine health operation.
Summary of the invention
The purpose of the present invention is to provide a kind of influences for considering winding copper loss to temperature, in larger thermal unbalance coefficient
Under the influence of can accurately calculate the nuclear power rotor of steam turbo generator thermo parameters method calculation method of temperature of rotor field distribution, to solve
Technical problem present in above-mentioned background technique.
To achieve the goals above, this invention takes following technical solutions:
A kind of nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct provided by the invention, the party
Method includes following process step:
Step S110: Generator Rotor Temperature field computation model is established；
Step S120: the copper loss and iron loss of generator amature are calculated；
Step S130: according to the copper loss and iron loss, in conjunction with the Generator Rotor Temperature field computation model, utilization is limited
Volumetric method determines the temperature field of the generator amature operation, obtains the draught temperature of mean temperature and air duct；
Step S140: the copper loss of generator amature is recalculated according to the mean temperature, and in conjunction with the draught temperature
It is iterated calculating using finite volume method, until the difference of the temperature value of mean temperature and calculating copper loss is less than given threshold.
Further, the step S110 is specifically included:
According to the practical structures and size of steam turbine generator, the three-dimensional entity model of Generator Rotor Temperature field is established；It is right
The three-dimensional entity model carries out grid dividing, establishes corresponding Generator Rotor Temperature field computation model；
The Generator Rotor Temperature field computation model include rotor core, winding inner copper line, turn-to-turn insulation, major insulation,
Slot wedge, air duct and cushion block；
The rotor core includes rotor tooth, and the adjacent rotor tooth is the arcuate structure of unequal spacing, and two is adjacent
There is cuboid slot between arcuate structure, winding is equipped in the cuboid slot, the winding includes copper bar, turn-to-turn insulation, master
Insulation and slot wedge, the copper bar and the turn-to-turn insulation are arranged alternately in the cuboid slot, and the slot wedge is arranged in close
On the upside of the copper bar of the notch of cuboid slot, major insulation is arranged in the outside of copper bar, turn-to-turn insulation and slot wedge, offsets with rotor tooth
It touches, circulation has hydrogen-cooled fluid in the air duct, and the air duct has cushion block in exit region to make the hydrogen-cooled fluid from diameter
It is flowed out to air outlet.
Further, the step S120 is specifically included:
The iron loss includes basic iron loss and additional iron loss, sets the initial calculation temperature of copper loss and iron losses computation, according to
Loss calculation equation calculates copper loss and basic iron loss based on Theory of Electromagnetic Field method；
The additional iron loss includes that stator field the higher hamonic wave short circuit loss, the stator slot ripples that generate in rotor surface exist
No-load loss and rotor that short circuit loss that rotor surface generates, stator slot ripples are generated in rotor surface and hydrogen-cooled fluid
Frictional dissipation.
Further, which is characterized in that the calculation method of the additional iron loss specifically includes:
Wherein, Q2SHIndicate additional iron loss, Q2vKIndicate the short circuit loss that stator field higher hamonic wave is generated in rotor surface,
Q2zKIndicate the short circuit loss that stator slot ripples is generated in rotor surface, Q2zoIndicate the sky that stator slot ripples is generated in rotor surface
Load-loss, Q2HIndicate the frictional dissipation of rotor and hydrogen-cooled fluid；
Wherein,Indicate that stator field higher hamonic wave generates the design factor of loss, AS in rotor surface1Indicate stator
Specific electric load, KδIndicate that Carter's coefficient, δ indicate unilateral air gap, DiIndicate that stator core internal diameter, p are magnetic field number of pole-pairs, l2It indicates to turn
Sub- body length, f indicate rotor rotation frequency；
Wherein,Indicate that stator slot ripples generates the design factor of loss, Z in rotor surface1Indicate number of stator slots；
Wherein, BδIndicate the magnetic flux density in air gap, KδnIndicate the Carter's coefficient of stator slot, lefIndicate that stator core is grown only
Degree；
Wherein, D2Indicate rotor diameter, PHIndicate the pressure of hydrogen-cooled fluid, PHOIndicate the normal hydrogen pressure in motor.
Further, the step S130 is specifically included:
Generator Rotor Temperature field computation model is calculated, the Generator Rotor Temperature field computation is calculated separately out
Heat transfer equation, boundary face and the inlet and outlet of model.
Further, the calculating of the heat transfer equation is as follows:
To the Generator Rotor Temperature field computation model, using standard fluid structurecoupling k- ε modeling Equations of Turbulence, benefit
It is as follows with standard fluid structurecoupling k- ε modeling Equations of Turbulence:
In formula, k is Turbulent Kinetic, and ε is invasin, and ρ is fluid density, and V is fluid velocity vectors, and t is time, GkFor
Turbulent generation rate, utFor turbulent flow viscosity, G1ε、G2εFor constant, σkAnd σεFor turbulent Planck's constant；
It is described to meet the following conditions using standard fluid structurecoupling k- ε modeling Equations of Turbulence:
Mass-conservation equation:
Momentum conservation equation:
Energy conservation equation:
Wherein, vrFor relative velocity vector, r is the position vector of the micro unit in rotating coordinate system, and p is to act on air
Static pressure on micro unit, τ are the viscous stress for acting on micro unit surface that generates due to molecular viscosity effect, ρ (2 Ω ×
vr+ Ω × Ω × r) it is Ke Liaolili, F is the body force on micro unit, and T is temperature, and v is absolute velocity, and λ is thermally conductive system
Number, c is specific heat at constant pressure, SrThe ratio of the heat and c that are generated for unit volume inner heat source；
The Generator Rotor Temperature field is calculated using fluid-wall interaction method, then the heat transfer side of Generator Rotor Temperature field
Journey are as follows:
In formula, λx、λy、λzIt is the heat transfer coefficient on different directions respectively；qVFor internal heat resource density.
Further, the boundary face includes radiating surface and insulation face；
The boundary face and inlet and outlet meet in the boundary condition of the Generator Rotor Temperature field computation model: described
In Generator Rotor Temperature field computation model, the insulation face meets:
In formula, T is the temperature of object, and n is boundary normal vector；
The radiating surface meets:
In formula, λ is thermal coefficient, and α is coefficient of heat transfer, TfFor ambient temperature；
The inlet and outlet meet the specified criteria on inlet pressure boundary and outlet pressure boundary, the inlet pressure and described
Outlet pressure is determined according to the productive experiment value of the steam turbine generator.
Further, the step S140 is specifically included:
According to the heat transfer equation, boundary face and inlet and outlet, the temperature of rotor regularity of distribution is obtained, according to the flat of rotor windings
Equal temperature calculates the close value of heat of copper loss of rotor consumption, is attached to rotor windings, and additional calculated turn for the first time again
The wind-warm syndrome of sub- air outlet carries out Temperature calculating, until calculate the close value of heat of copper loss of rotor consumption and rotor temperature field calculate around
The mean temperature difference of group is less than preset threshold value, to obtain the higher rotor temperature field of precision.
The invention has the advantages that: the calculation model for temperature field by building rotor, and using finite volume method, analytical calculation turns
Temperature field when son operates normally, has obtained mean temperature, the maximum temperature of the winding of the Temperature Distribution of rotor, especially rotor
Meter is utilized since rotor mean temperature and the temperature value of calculating heat flow density use are inconsistent with the wind-warm syndrome of air duct air outlet
The mean temperature of calculating recalculates heat flow density, and adds the wind of new heat flow density and rotor ventilation ditch air outlet
Temperature, carries out the first time iteration in temperature field, and iteration is to using the calculated temperature value of Finite Volume Method and calculate heat flow density
Temperature value it is almost the same, effectively raise the computational accuracy in temperature field, make the calculating of the biggish motor of thermal unbalance coefficient
As a result it is more in line with reality, practical to engineering have directiveness.
The additional aspect of the present invention and advantage will be set forth in part in the description, these will become from the following description
Obviously, or practice through the invention is recognized.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill of field, without creative efforts, it can also be obtained according to these attached drawings others
Attached drawing.
Fig. 1 is the nuclear power rotor of steam turbo generator thermo parameters method calculating side that air duct is had described in the embodiment of the present invention
Method flow chart.
Fig. 2 is the nuclear power rotor of steam turbo generator calculation model for temperature field knot that air duct is had described in the embodiment of the present invention
Composition.
Fig. 3 is that the nuclear power rotor of steam turbo generator air duct flow direction structure described in the embodiment of the present invention with air duct shows
It is intended to.
Fig. 4 is the nuclear power rotor of steam turbo generator calculation model for temperature field structure that air duct is had described in the embodiment of the present invention
Construction method flow chart.
Wherein: 1- rotor tooth, 2- rotor yoke, 3- copper bar rotor, 4- turn-to-turn insulation, 5- slot wedge, 6- major insulation, 7- ventilation
Ditch, 8- cushion block.
Specific embodiment
Embodiments of the present invention are described below in detail, the example of the embodiment is shown in the accompanying drawings, wherein from beginning
Same or similar element or module with the same or similar functions are indicated to same or similar label eventually.Below by ginseng
The embodiment for examining attached drawing description is exemplary, and for explaining only the invention, and is not construed as limiting the claims.
Those skilled in the art of the present technique are appreciated that unless expressly stated, singular " one " used herein, " one
It is a ", " described " and "the" may also comprise plural form.It is to be further understood that being arranged used in specification of the invention
Diction " comprising " refer to that there are the feature, integer, step, operation, element and/or modules, but it is not excluded that in the presence of or addition
Other one or more features, integer, step, operation, element, module and/or their group.
Those skilled in the art of the present technique are appreciated that unless otherwise defined, all terms used herein (including technology art
Language and scientific term) there is meaning identical with the general understanding of those of ordinary skill in fields of the present invention.Should also
Understand, those terms such as defined in the general dictionary, which should be understood that, to be had and the meaning in the context of the prior art
The consistent meaning of justice, and unless defined as here, it will not be explained in an idealized or overly formal meaning.
In order to facilitate understanding of embodiments of the present invention, further by taking specific embodiment as an example below in conjunction with attached drawing to be solved
Explanation is released, and embodiment does not constitute the restriction to the embodiment of the present invention.
Those of ordinary skill in the art are it should be understood that attached drawing is the schematic diagram of one embodiment, the portion in attached drawing
Part or device are not necessarily implemented necessary to the present invention.
Embodiment
As shown in Figure 1, the embodiment of the present invention provides a kind of vast capacity nuclear power rotor of steam turbo generator axial ventilation system
Multiple physical field coupling calculation, this method make to calculate by the wind-warm syndrome of the heat close value and air outlet of iteration difference winding temperature
As a result more accurately, rotor of steam turbo generator higher for thermal unbalance coefficient and that actual measurement difficulty is big, can preferably meet
Engineering actual demand.
Specifically, building nuclear power turbine generator rotor 1/8 is circumferential first, and half of three-dimensional entity model axially grown, and it is right
The three-dimensional entity model carries out grid dividing, establishes corresponding Generator Rotor Temperature field computation model；
Copper bar copper loss, the basic iron loss of teeth portion and the additional iron loss of nuclear power steam turbine generator are calculated, and by copper bar copper loss, teeth portion
It is close that basic iron loss and additional iron loss are separately converted to heat, then by the close corresponding position for being attached to the calculation model for temperature field of heat；
Generator Rotor Temperature field computation model is calculated, the flat of the Temperature Distribution of rotor, especially winding is obtained
Equal temperature；
According to the mean temperature of rotor windings, the copper loss of copper bar is recalculated, and be converted into heat it is close, be attached to rotor temperature
Field computation model, while the wind-warm syndrome of additional first calculated rotor ventilation ditch air outlet are spent, to rotor ventilation ditch air outlet
Wind-warm syndrome value is added for first time calculating as a result, carrying out Finite Volume Method calculating, until the average temperature of calculated rotor windings
Angle value differs very little with the temperature value for calculating copper loss, and iteration finishes.
As shown in Fig. 2, coupling meter according to the multiple physical field of the axial complicated ventilating system of vast capacity rotor of steam turbo generator
The three-dimensional entity model in calculation method building nuclear power rotor of steam turbo generator temperature field, comprising:
According to the practical structures and size of nuclear power steam turbine generator, the 3D solid mould of Generator Rotor Temperature field is established
Type；
Three-dimensional entity model respectively include: rotor core, winding inner copper line, turn-to-turn insulation, major insulation, slot wedge and ventilation
Ditch.
Grid dividing is carried out to the three-dimensional entity model, establishes corresponding Generator Rotor Temperature field computation model, is wrapped
It includes:
Establish the Generator Rotor Temperature field computation model for corresponding to the three-dimensional entity model of the generator amature；
Generator Rotor Temperature field computation model, comprising: rotor core, winding inner copper line, turn-to-turn insulation, major insulation, slot
Wedge and air duct；
Rotor core includes 2 two parts of rotor tooth 1 and rotor yoke, constitutes the main part of the three-dimensional entity model,
Rotor tooth is connect with the rotor yoke, and rotor tooth is the arcuate structure of unequal spacing, is had among two adjacent arcuate structures rectangular
The slot of body structure placed winding in the slot of rectangular parallelepiped structure, include copper bar 3, turn-to-turn insulation 4, major insulation 6 and slot wedge in winding
5, copper bar and turn-to-turn insulation are arranged alternately in slot, and slot wedge is arranged on the upside of the copper bar of notch, and major insulation is arranged in neatly
The outside of the copper bar of arrangement, turn-to-turn insulation and slot wedge, closely rotor tooth have air duct 7 in the copper bar, have hydrogen in air duct
Gas flows through, and has cushion block 8 in exit region, flows out hydrogen from radial air outlet, air duct schematic diagram is as shown in Figure 3.
Calculate nuclear power steam turbine generator copper bar copper loss, the basic iron loss of teeth portion and additional iron loss, and by the copper bar copper loss,
It is close that the basic iron loss of teeth portion and additional iron loss are separately converted to heat, then by the close correspondence portion for being attached to the calculation model for temperature field of heat
Position, comprising:
The copper bar of nuclear power rotor of steam turbo generator is direct current generation, only basic copper loss, the initial calculation of copper loss
Temperature is 75 DEG C general in engineering.
According to the expression formula of loss calculation, copper loss and basic iron loss are calculated based on Theory of Electromagnetic Field method.
The iron core additional iron losses of rotor consumes, and is the relevant knowledge based on electromagnetic theory and mechanical theory, using following formula
It is calculated:
The added losses Q on rotor core surface2SHMainly have: 1) stator field higher hamonic wave generates short in rotor surface
Path loss consumes Q2vK；2) the short circuit loss Q that stator slot ripples is generated in rotor surface2zk；3) stator slot ripples is generated in rotor surface
No-load loss Q2zo；4) the frictional dissipation Q of rotor and cooling gas2H, unit of account kW:
In formula,The design factor of loss is generated in rotor surface for stator field higher hamonic wave；AS1It is negative for stator line
Lotus (A/cm)；KδFor Carter's coefficient；δ is unilateral air gap (mm)；DiFor stator core internal diameter；P is grade logarithm；l2For rotor body
Length (mm)；F is frequency (Hz).
In formula,The design factor of loss is generated in rotor surface for stator slot ripples；Z1For number of stator slots.
In formula, BδFor the magnetic flux density (Gs) in air gap；KδnFor the Carter's coefficient of stator slot；lefFor stator core clear length
(mm)。
In formula, D2For rotor diameter (mm)；PHFor used hydrogen pressure；PHOFor the normal hydrogen pressure in motor.
Copper bar copper loss, the iron of the nuclear power steam turbine generator are calculated based on Theory of Electromagnetic Field method and machine theoretical method
The basic iron loss of core and additional iron loss；
By the basic iron loss of copper bar copper loss, iron core for the nuclear power steam turbine generator being calculated and additional iron loss, add respectively
To Generator Rotor Temperature field computation model.
Generator Rotor Temperature field computation model is calculated, obtain corresponding rotor is operated normally and overloaded three
Tie up temperature field, comprising:
Calculate separately out heat transfer equation, boundary face and the inlet and outlet of the Generator Rotor Temperature field computation model.
Generator Rotor Temperature field computation model is calculated, the Generator Rotor Temperature field computation is calculated separately out
Heat transfer equation, boundary face and the inlet and outlet of model, comprising:
To the Generator Rotor Temperature field computation model, using the fluid structurecoupling k- ε modeling Equations of Turbulence of standard,
Specifically meet following formula (6)~(8):
Mass-conservation equation:
Momentum conservation equation:
Energy conservation equation:
Wherein, ρ is density, vrFor relative velocity vector, r is the position vector of the micro unit in rotating coordinate system, and p is to make
For the static pressure on air micro unit, τ is the viscous stress for acting on micro unit surface generated by molecular viscosity effect,
ρ(2Ω×vr+ Ω × Ω × r) it is Ke Liaolili, F is the body force on micro unit, and T is temperature, and v is absolute velocity, and λ is
Thermal coefficient, c are specific heat at constant pressure, SrThe ratio of the heat and c that are generated for unit volume inner heat source；
It is as follows using the fluid structurecoupling k- ε modeling Equations of Turbulence of standard:
In formula, k is Turbulent Kinetic, and ε is invasin, and ρ is fluid density, and V is fluid velocity vectors, and t is time, GkFor
Turbulent generation rate, utFor turbulent flow viscosity, G1ε、G2εFor constant, σkAnd σεFor turbulent Planck's constant；
Generator Rotor Temperature field is calculated using fluid and heat transfer coupling process, then the heat transfer side of Generator Rotor Temperature field
Journey are as follows:
In formula, λx、λy、λzIt is the heat transfer coefficient on different directions respectively；qVFor internal heat resource density.
Generator Rotor Temperature field computation model is calculated, Generator Rotor Temperature field computation model is calculated separately out
Heat transfer equation, boundary face and inlet and outlet, further includes:
Boundary face includes: radiating surface and insulation face；
The boundary condition of boundary face and inlet and outlet in Generator Rotor Temperature field is as follows:
In Generator Rotor Temperature field computation model, the insulation face is four sides, comprising: S2, S3, S4, S6, insulation face point
Do not meet the following conditions:
In formula, T is the temperature of object, and n is boundary normal vector；
Radiating surface is two sides, comprising: S1And S5, meet the following conditions respectively:
In formula, λ is thermal coefficient, and α is coefficient of heat transfer, TfFor ambient temperature；
Inlet and outlet are respectively as follows: S7And S8, import and export the given item for meeting inlet pressure boundary and outlet pressure boundary respectively
Part, inlet pressure and outlet pressure are provided according to nuclear power Testing of Turbo-generator value.
As shown in figure 4, the inventive embodiments, which provide Finite Volume Method, calculates rotor of steam turbo generator temperature field, place
Managing process includes following processing step: according to the practical structures and size of the nuclear power rotor of steam turbo generator of a 1100MW,
1/8 circumference of rotor of foundation, the three-dimensional entity model of half of axial direction part are established as shown in Fig. 2, dividing to physical model
Corresponding to the Generator Rotor Temperature field computation model of Generator Rotor Temperature field physical model, as shown in Figure 2.
In Fig. 2, the computation model include: rotor core, winding inner copper line, turn-to-turn insulation, major insulation, slot wedge,
Air duct and cushion block.
According to electromagnetic theory and mechanical theory, copper loss of rotor consumption, basic iron loss and additional iron losses consumption are calculated, is converted into
Heat is close, is attached to corresponding position.
According to heat transfer equation, boundary face and the inlet and outlet of the rotor three-dimensional temperature field, the temperature of rotor regularity of distribution is obtained,
According to the mean temperature of rotor windings, the close value of heat of copper loss of rotor consumption is calculated, is attached to rotor windings again, and add
The wind-warm syndrome of calculated rotor air outlet for the first time carries out Temperature calculating, until calculating the hot close value of copper loss of rotor consumption and turning
The mean temperature that sub- Temperature calculating goes out winding is very nearly the same, to obtain the higher rotor temperature field of precision.
The embodiment of the invention provides the multiple physical field couplings of the axial complicated ventilating system of vast capacity rotor of steam turbo generator
Calculation method is closed, is by taking the new model of large nuclear-power rotor of steam turbo generator as an example；But this method is not limited solely to nuclear power vapour
The rotor portion of turbine generator, while being suitable for the biggish motor temperature field computation of any kind thermal unbalance coefficient, have wide
General applicability.
The multiple physical field of the axial complicated ventilating system of vast capacity rotor of steam turbo generator provided by the invention, which couples, to be calculated
Method has the advantage that compared with traditional Steady-State Thermal Field calculation method
The present invention has the three-dimensional temperature field computation model of the rotor of steam turbo generator of axial ventilation system by establishing, and obtains
Out to draw a conclusion: the iterative calculation in temperature field is carried out by the close calculated value of heat generated to the copper loss under different temperatures, it can be with
The accuracy for greatly improving temperature of rotor field computation can preferably refer to the higher steam turbine generator of thermal unbalance coefficient
Lead its safe operation.
In conclusion the three-dimensional temperature of rotor of steam turbo generator of the embodiment of the present invention by building with axial ventilation system
Field computation model is spent, using finite volume method, the temperature that analytical calculation rotor calculates the close value of heat is consistent with Calculated Results of Temperature Field,
A kind of new approach is provided for the accurate grasp Temperature Distribution of steam turbine generator in operation, and greatly improves judgement hair
The accuracy of motor safe operation.
As seen through the above description of the embodiments, those skilled in the art can be understood that the present invention can
It realizes by means of software and necessary general hardware platform.Based on this understanding, technical solution of the present invention essence
On in other words the part that contributes to existing technology can be embodied in the form of software products, the computer software product
It can store in storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions are used so that a computer equipment
(can be personal computer, server or the network equipment etc.) executes the certain of each embodiment or embodiment of the invention
Method described in part.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
In the technical scope disclosed by the present invention, any changes or substitutions that can be easily thought of by anyone skilled in the art,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with scope of protection of the claims
Subject to.
Claims (8)
1. a kind of nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct, which is characterized in that including such as
Lower process step:
Step S110: Generator Rotor Temperature field computation model is established；
Step S120: the copper loss and iron loss of generator amature are calculated；
Step S130: limited bulk is utilized in conjunction with the Generator Rotor Temperature field computation model according to the copper loss and iron loss
Method determines the temperature field of the generator amature operation, obtains the draught temperature of mean temperature and air duct；
Step S140: the copper loss of generator amature is recalculated according to the mean temperature, and is utilized in conjunction with the draught temperature
Finite volume method is iterated calculating, until the difference of the temperature value of mean temperature and calculating copper loss is less than given threshold.
2. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 1 with air duct,
It is characterized in that, the step S110 is specifically included:
According to the practical structures and size of steam turbine generator, the three-dimensional entity model of Generator Rotor Temperature field is established；To described
Three-dimensional entity model carries out grid dividing, establishes corresponding Generator Rotor Temperature field computation model；
The Generator Rotor Temperature field computation model include rotor core, winding inner copper line, turn-to-turn insulation, major insulation, slot wedge,
Air duct and cushion block；
The rotor core includes rotor tooth, and the adjacent rotor tooth is the arcuate structure of unequal spacing, two adjacent arcs
There is cuboid slot between structure, winding is equipped in the cuboid slot, the winding includes copper bar, turn-to-turn insulation, major insulation
And slot wedge, the copper bar and the turn-to-turn insulation are arranged alternately in the cuboid slot, the slot wedge is arranged in close to rectangular
On the upside of the copper bar of the notch of body slot, major insulation is arranged in the outside of copper bar, turn-to-turn insulation and slot wedge, institute inconsistent with rotor tooth
Stating circulation in air duct has hydrogen-cooled fluid, and the air duct has cushion block in exit region to make the hydrogen-cooled fluid from radial outlet air
Mouth outflow.
3. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 2 with air duct,
It is characterized in that, the step S120 is specifically included:
The iron loss includes basic iron loss and additional iron loss, the initial calculation temperature of copper loss and iron losses computation is set, according to loss
Accounting equation calculates copper loss and basic iron loss based on Theory of Electromagnetic Field method；
The additional iron loss include stator field higher hamonic wave generated in rotor surface short circuit loss, stator slot ripples is in rotor
The friction of no-load loss and rotor and hydrogen-cooled fluid that short circuit loss, the stator slot ripples of surface generation are generated in rotor surface
Loss.
4. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 3 with air duct,
It is characterized in that, the calculation method of the additional iron loss specifically includes:
Wherein, Q2SHIndicate additional iron loss, Q2vKIndicate the short circuit loss that stator field higher hamonic wave is generated in rotor surface, Q2zK
Indicate the short circuit loss that stator slot ripples is generated in rotor surface, Q2zoIndicate the zero load that stator slot ripples is generated in rotor surface
Loss, Q2HIndicate the frictional dissipation of rotor and hydrogen-cooled fluid；
Wherein,Indicate that stator field higher hamonic wave generates the design factor of loss, AS in rotor surface1Indicate that stator line is negative
Lotus, KδIndicate that Carter's coefficient, δ indicate unilateral air gap, DiIndicate that stator core internal diameter, p are magnetic field number of pole-pairs, l2Indicate rotor sheet
Body length, f indicate rotor rotation frequency；
Wherein,Indicate that stator slot ripples generates the design factor of loss, Z in rotor surface1Indicate number of stator slots；
Wherein, BδIndicate the magnetic flux density in air gap, KδnIndicate the Carter's coefficient of stator slot, lefIndicate stator core clear length；
Wherein, D2Indicate rotor diameter, PHIndicate the pressure of hydrogen-cooled fluid, PHOIndicate the normal hydrogen pressure in motor.
5. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 4 with air duct,
It is characterized in that, the step S130 is specifically included:
Generator Rotor Temperature field computation model is calculated, the Generator Rotor Temperature field computation model is calculated separately out
Heat transfer equation, boundary face and inlet and outlet.
6. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 5 with air duct,
It is characterized in that, the calculating of the heat transfer equation is as follows:
Mark is utilized using standard fluid structurecoupling k- ε modeling Equations of Turbulence to the Generator Rotor Temperature field computation model
Quasi- fluid structurecoupling k- ε modeling Equations of Turbulence is as follows:
In formula, k is Turbulent Kinetic, and ε is invasin, and ρ is fluid density, and V is fluid velocity vectors, and t is time, GkFor turbulent flow
Generation rate, utFor turbulent flow viscosity, G1ε、G2εFor constant, σkAnd σεFor turbulent Planck's constant；
It is described to meet the following conditions using standard fluid structurecoupling k- ε modeling Equations of Turbulence:
Mass-conservation equation:
▽(ρvr)=0；
Momentum conservation equation:
Energy conservation equation:
Wherein, vrFor relative velocity vector, r is the position vector of the micro unit in rotating coordinate system, and p is to act on air infinitesimal
Static pressure on body, τ are the viscous stress for acting on micro unit surface generated by molecular viscosity effect, ρ (2 Ω × vr+Ω
× Ω × r) it is Ke Liaolili, F is the body force on micro unit, and T is temperature, and v is absolute velocity, and λ is thermal coefficient, and c is
Specific heat at constant pressure, SrThe ratio of the heat and c that are generated for unit volume inner heat source；
The Generator Rotor Temperature field is calculated using fluid-wall interaction method, then the heat transfer equation of Generator Rotor Temperature field are as follows:
In formula, λx、λy、λzIt is the heat transfer coefficient on different directions respectively；qVFor internal heat resource density.
7. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 5 with air duct,
It is characterized in that,
The boundary face includes radiating surface and insulation face；
The boundary face and inlet and outlet meet in the boundary condition of the Generator Rotor Temperature field computation model: in the power generation
In machine rotor calculation model for temperature field, the insulation face meets:
In formula, T is the temperature of object, and n is boundary normal vector；
The radiating surface meets:
In formula, λ is thermal coefficient, and α is coefficient of heat transfer, TfFor ambient temperature；
The inlet and outlet meet the specified criteria on inlet pressure boundary and outlet pressure boundary, the inlet pressure and the outlet
Pressure is determined according to the productive experiment value of the steam turbine generator.
8. the nuclear power rotor of steam turbo generator thermo parameters method calculation method according to claim 7 with air duct,
It is characterized in that, the step S140 is specifically included:
According to the heat transfer equation, boundary face and inlet and outlet, the temperature of rotor regularity of distribution is obtained, according to the average temperature of rotor windings
Degree calculates the close value of heat of copper loss of rotor consumption, is attached to rotor windings again, and additional rotor calculated for the first time goes out
The wind-warm syndrome in air port carries out Temperature calculating, until the close value of heat of calculating copper loss of rotor consumption and rotor temperature field calculate winding
Mean temperature difference is less than preset threshold value, to obtain the higher rotor temperature field of precision.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811081735.6A CN109359339A (en) | 2018-09-17 | 2018-09-17 | Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811081735.6A CN109359339A (en) | 2018-09-17 | 2018-09-17 | Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109359339A true CN109359339A (en) | 2019-02-19 |
Family
ID=65350915
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811081735.6A Pending CN109359339A (en) | 2018-09-17 | 2018-09-17 | Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109359339A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110031506A (en) * | 2019-04-25 | 2019-07-19 | 北京交通大学 | The calculation method of magneto air gap thermal coefficient |
| CN110688738A (en) * | 2019-09-10 | 2020-01-14 | 北京交通大学 | Method for calculating temperature field of nuclear power turbonator with axial ventilation system |
| CN111928965A (en) * | 2020-07-13 | 2020-11-13 | 中广核核电运营有限公司 | Rotor temperature monitoring method and device, computer equipment and storage medium |
| CN112327603A (en) * | 2020-10-13 | 2021-02-05 | 南京航空航天大学 | Method for predicting thermal bending vibration in magnetic suspension bearing rotor system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105245041A (en) * | 2015-10-15 | 2016-01-13 | 北京交通大学 | Method for determining surface heat dissipation coefficients of ventilating ducts of stator and rotor of large motor |
| CN106126861A (en) * | 2016-07-18 | 2016-11-16 | 北京交通大学 | The computational methods of the coefficient of heat transfer in air-cooled steam turbine generator teeth portion and back of the yoke portion |
| CN107301302A (en) * | 2017-07-12 | 2017-10-27 | 北京交通大学 | A kind of computational methods of Stator Coil added losses |
-
2018
- 2018-09-17 CN CN201811081735.6A patent/CN109359339A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105245041A (en) * | 2015-10-15 | 2016-01-13 | 北京交通大学 | Method for determining surface heat dissipation coefficients of ventilating ducts of stator and rotor of large motor |
| CN106126861A (en) * | 2016-07-18 | 2016-11-16 | 北京交通大学 | The computational methods of the coefficient of heat transfer in air-cooled steam turbine generator teeth portion and back of the yoke portion |
| CN107301302A (en) * | 2017-07-12 | 2017-10-27 | 北京交通大学 | A kind of computational methods of Stator Coil added losses |
Non-Patent Citations (1)
| Title |
|---|
| 苏营等: "1100 MW半速汽轮发电机转子绕组电阻率对转子温度场的影响", 《电机技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110031506A (en) * | 2019-04-25 | 2019-07-19 | 北京交通大学 | The calculation method of magneto air gap thermal coefficient |
| CN110688738A (en) * | 2019-09-10 | 2020-01-14 | 北京交通大学 | Method for calculating temperature field of nuclear power turbonator with axial ventilation system |
| CN111928965A (en) * | 2020-07-13 | 2020-11-13 | 中广核核电运营有限公司 | Rotor temperature monitoring method and device, computer equipment and storage medium |
| CN112327603A (en) * | 2020-10-13 | 2021-02-05 | 南京航空航天大学 | Method for predicting thermal bending vibration in magnetic suspension bearing rotor system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106126861B (en) | The calculation method of the coefficient of heat transfer in air-cooled steam turbine generator teeth portion and back of the yoke portion | |
| CN109359339A (en) | Nuclear power rotor of steam turbo generator thermo parameters method calculation method with air duct | |
| CN109101693B (en) | Method for calculating rotor multi-physical field based on air path-fluid-heat transfer cooperative coupling | |
| CN107301302B (en) | A kind of calculation method of Stator Coil added losses | |
| CN105245041B (en) | A method of solving large-size machine rotor air duct surface coefficient of heat transfer | |
| Tong et al. | Totally enclosed self-circulation axial ventilation system design and thermal analysis of a 1.65-MW direct-drive PMSM | |
| Li et al. | Influence of rotor radial ventilation ducts number on temperature distribution of rotor excitation winding and fluid flow state between two poles of a fully air-cooled hydro-generator | |
| Weili et al. | Influence of rotation on rotor fluid and temperature distribution in a large air-cooled hydrogenerator | |
| Lu et al. | Simulation and analysis of thermal fields of rotor multislots for nonsalient-pole motor | |
| CN109829220A (en) | The multi- scenarios method calculation method of the Rotor's Transient Temperature of air-cooled steam turbine generator | |
| Yi et al. | Self-pumped air-cooling design for a high-speed high-specific-power motor | |
| Xia et al. | Analyzing temperature rise and fluid flow of high-power-density and high-voltage induction motor in the starting process | |
| Weili et al. | Influence of the end ventilation structure change on the temperature distribution in the end region of large water–hydrogen–hydrogen cooled turbogenerator | |
| CN110688738A (en) | Method for calculating temperature field of nuclear power turbonator with axial ventilation system | |
| Xu et al. | Heat transfer characteristics of external ventilated path in compact high-voltage motor | |
| Wang et al. | Coupled fluid-solid heat transfer of a gas and liquid cooling PMSM including rotor rotation | |
| Jichao et al. | Calculation and analysis of complex fluid flow and thermal fields in a fully air-cooled hydrogenerator | |
| Su et al. | Influence of incident angle of hydrogen in the ventilation ducts on multi-physical fields of 1100 MW turbine-generator rotor | |
| Xu et al. | Research on heat transfer of submersible motor based on fluid network decoupling | |
| Li et al. | Multi-physical fields of rotor windings with axial-radial ventilation system for 1100 MW nuclear half-speed turbine generator | |
| Kang et al. | Thermal analysis of the cooling system with the circulation between rotor holes of enclosed PMSMs based on modified models | |
| Ai et al. | Research on thermal characteristics of internal ventilated paths in compact medium high-voltage motor based fluid network decoupling | |
| Ge et al. | Thermal analysis of axial‐radial hybrid ventilation motor and stator ventilation channel improvement | |
| Ryzhov et al. | Three-dimensional thermal stator model of a fully air-cooled turbogenerator | |
| Lee et al. | Thermofluidic analysis of interior permanent magnet synchronous motors with internal air circulation by protrusion-shaped flow inducers for effective thermal management |
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 |
Application publication date: 20190219 |
|
| RJ01 | Rejection of invention patent application after publication |