CN112560367A - Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator - Google Patents

Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator Download PDF

Info

Publication number
CN112560367A
CN112560367A CN202011361490.XA CN202011361490A CN112560367A CN 112560367 A CN112560367 A CN 112560367A CN 202011361490 A CN202011361490 A CN 202011361490A CN 112560367 A CN112560367 A CN 112560367A
Authority
CN
China
Prior art keywords
hole
axial
air
wind
coil
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
CN202011361490.XA
Other languages
Chinese (zh)
Other versions
CN112560367B (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.)
Shanghai Electric Power Generation Equipment Co Ltd
Original Assignee
Shanghai Electric Power Generation Equipment 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 Shanghai Electric Power Generation Equipment Co Ltd filed Critical Shanghai Electric Power Generation Equipment Co Ltd
Priority to CN202011361490.XA priority Critical patent/CN112560367B/en
Publication of CN112560367A publication Critical patent/CN112560367A/en
Application granted granted Critical
Publication of CN112560367B publication Critical patent/CN112560367B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Windings For Motors And Generators (AREA)
  • Motor Or Generator Cooling System (AREA)

Abstract

According to the calculation method for the hot spot temperature of the axial double-hole copper bar rotor coil of the generator, the calculation process is modularly designed according to functions, the whole calculation process is divided into the input data module, the wind path iterative calculation module and the loss and temperature rise calculation module, the calculation method has the advantages that the function modules are clear, the calculation program is easy to use and upgrade, the pertinence is strong, and the calculation precision is high aiming at the temperature rise calculation of the double-hole copper bar structure. Meanwhile, the calculation of the temperature rise of the surface of the outlet of the radial wind hole surface heat dissipation and the axial wind hole surface heat dissipation is considered, the method is suitable for the temperature rise design of the existing two-pole 800 MW-level to 1260 MW-level generator, can be used for developing generators adopting axial double-hole copper bar rotor coils in more capacity ranges, and can meet the new calculation method of the rotor hot spot temperature, which has higher application degree and can ensure the precision, in the product optimization and improvement, so as to meet the calculation requirement of the generator adopting the axial double-hole copper bar rotor.

Description

Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator
Technical Field
The invention relates to a method for calculating the hot spot temperature of an axial double-hole copper bar rotor coil of a generator, which can be used for calculating and researching the hot spot temperature of the axial double-hole copper bar rotor coil of the generator and developing and researching a generator product and belongs to the technical field of electromagnetic design of the generator.
Background
The hot spot temperature of the rotor coil of the generator is one of main design data of the rotor coil of the generator, the accuracy of the design value not only influences the determination of the main size of the rotor of the generator, but also influences the operation reliability and the service life of the generator, and the accurate calculation of the hot spot temperature of the rotor coil becomes an important link in the design process of a generator product.
The prior art does not have a special calculation method for the double-hole copper bar, and originally adopts the introduced Siemens program TGS3248 to carry out approximate calculation, and the calculation method is a calculation method for calculating a single-hole copper bar which is formed by combining an upper half copper bar and a lower half copper bar into a circle. Because the structure of the axial double-hole copper bar is inconsistent with the structure, the axial double-hole copper bar is not completely applicable in practice. First, the air volume in the ventilation hole cannot be accurately calculated because the air path structures of the two are different. Secondly, the copper loss of the rotor coil at the air outlet cannot be calculated. The heat dissipation coefficient at the air outlet cannot be calculated. Therefore, the hot-spot temperature rise of the rotor coil cannot be accurately calculated. Further, when the TGS3248 is used for calculation, various assumptions are made, and the use is complicated.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the prior art does not have a special calculation method for the hot spot temperature of the double-hole copper bar rotor coil.
In order to solve the problems, the technical scheme of the invention is to provide a method for calculating the hot spot temperature of a generator axial double-hole copper bar rotor coil, which is characterized by comprising the following steps of:
step 1, inputting calculation parameters and calculating the structure size:
the calculation parameters comprise:
number of one pole coil: n _ coil;
number of turns per slot of rotor: NRS;
the rotor body is long: lp;
the arc part of the No. 1 coil is long from the body: l iscoil1_to_body；
Coil end arc portion width: b iscoil_end；
Coil pitch: l iscoil_to_coil；
The distance between the near body side air inlet hole and the body is as follows: l isinvent_to_body；
Distance between the near-end-portion-side air inlet hole and the end coil: l isinvent_to_endturn；
Air outlet hole spacing: l isoutvent_to_outvent；
The axial vent hole is high: h isv；
Axial vent hole width: bv；
The copper wire is high: h isn；
Half the width of the double-hole copper wire: bn；
Copper line outer chamfer radius: r1；
Copper line inside chamfer radius: r3；
Wall thickness in the copper wire width direction: thick _ cir;
circumferential width of the air outlet hole: bot；
Axial length of the air outlet: lot；
Air outlet hole center distance: t is tot；
Circumferential width of the air inlet hole: bin；
Axial length of air inlet hole: lin；
Rated rotating speed of the generator: RPM;
depth of rotor groove: h isrs；
Rotor outer diameter: dr；
Pressure difference between the air inlet and the air outlet of the rotor coil: pENTR；
Half value of generator rated excitation: i isf；
Cold air temperature: t isamb；
Area of vent hole on rotor wedge under backing strip: a. the12；
Area of rotor slot wedge exhaust vent: a. the13；
Cooling gas pressure: p;
calculating the length of the linear end distance of the maximum coil from the body:
Lcoilmax_to_body＝Lcoil_to_body+(N_coil-1)*(Bcoil_end+Lcoil_to_coil)；
wall thickness in the copper wire height direction: WEB ═ h (h)n-hv)/2；
The height of the axial air hole side wall is protruded in the circumferential direction of the upper wall and the lower wall of the copper wire after the radial air outlet is formed:
bdiff＝thick_cir+bv-bot；
step 2, to the ventilation structure of diplopore copper bar, calculate the rotor wind speed size in the ventilation hole, judge whether satisfy the calculation accuracy requirement, specifically include:
step 2.1, average wind path length
Figure BDA0002804100780000031
Figure BDA0002804100780000032
Step 2.2, the sectional area A of the air duct of the rotor copper wirev＝hv·bv-(4-π)R3 2；
Wetted perimeter S of axial wind holev＝2·(hv+bv)-R3·(8-2π)；
Hydraulic diameter of axial wind hole
Figure BDA0002804100780000033
Step 2.3, air duct inlet area Ain＝Av；
Step 2.4, air duct outlet area Aov＝(lot·bot)；
Wetted perimeter S of radial wind holesov＝2(bot+lot)-lot；
Hydraulic diameter of radial wind hole
Figure BDA0002804100780000034
Step 2.5, area A of vent hole on filler strip at top turn conductor11＝Aov；
Step 2.6, air pressure drop coefficient of air duct inlet
Figure BDA0002804100780000035
Step 2.7, wind pressure drop coefficient K at the secondary turningbe＝2.584；
Step 2.8, the air pressure drop coefficient of the axial air duct
Figure BDA0002804100780000036
frIteratively determined by the results of the calculation of step 2.15.3;
step 2.9, carve the wind pressure drop coefficient of the bottom filler strip
Figure BDA0002804100780000037
Step 2.10, wind pressure drop coefficient of air outlet of slot wedge
Figure BDA0002804100780000041
Step 2.11, considering resistance coefficient increased by orifice plate effect of air outlet
Figure BDA0002804100780000042
Step 2.12, total wind pressure drop coefficient of rotor wind path
Kv＝KOutlet _ Orifice plate Effect _ Orifice _4+Kfr+(Kin+Kbe+K12+K13)；
Step 2.13 wind speed of rotor surface
Figure BDA0002804100780000043
Wind speed at rotor slot bottom
Figure BDA0002804100780000044
Centrifugal wind pressure
Figure BDA0002804100780000045
Step 2.14, total pressure drop Δ P in the air path is equal to PENTR + Pcen；
Step 2.15, a wind speed and wind amount calculation part comprises the following steps
Step 2.15.1, setting the wind speed of the axial vent hole as V1；
Figure BDA0002804100780000046
Step 2.15.2 Reynolds number Rev1＝P·V1·Dv；
Step 2.15.3, coefficient of friction
Figure BDA0002804100780000047
Figure BDA0002804100780000048
The deviation is less than 0.0001, if not, V in formula 2.15.1 needs to be adjusted1A value;
total wind pressure drop of step 2.15.4 wind path
Figure BDA0002804100780000049
Figure BDA00028041007800000410
The deviation is less than 0.0001, if not, V in formula 2.15.1 needs to be adjusted1A value;
step 2.15.5, air volume of each circle of copper wire
Figure BDA00028041007800000411
Step 3, loss and temperature rise calculation: calculating the loss of each circle of copper wire, the loss of the copper wire in unit length, the surface temperature rise of the axial vent hole and the gas temperature rise; and aiming at the air outlet area where the hot spot temperature is located, the gas speed and the surface temperature rise at the outlet are respectively calculated, and the method specifically comprises the following steps:
step 3.1, take away long loss of copper
Figure BDA0002804100780000051
Step 3.2 axial conductive sectional area A of copper wirecu＝bn·hn-(4-π)R1 2-Av；
Step 3.3 No radial hole section conductive section area A at air outletcuh1＝Acu；
Step 3.4, the air outlet part is provided with a radial hole section conductive sectional area
Acuh2＝Acu-2WEB·(bot-thick_cir)-thick_cir·hn；
Step 3.5, an air hole pitch internal resistance of a radial hole section is not arranged at the air outlet part:
temperature at rotor coil outlet is Tf_out_initial，Tf_out_initialIteratively determined from the results of the 3.33 calculations,
Figure BDA0002804100780000052
step 3.6, an air outlet part is provided with an inner resistor with a radial hole section and an air hole pitch
Figure BDA0002804100780000053
Step 3.7 air outlet one air hole pitch internal resistance Rov＝Rot1+Rot2；
Step 3.8 equivalent area of copper wire in air-out area
Figure BDA0002804100780000054
Step 3.9, multiplying power reduction of equivalent area of copper wire in air outlet area
Figure BDA0002804100780000055
Step 3.10, copper consumption of each circle of copper wire:
average temperature T of straight section of rotor coilcoil_avg_initial，Tcoil_avg_initialIteratively determined from the results of the 3.14 calculations,
Figure BDA0002804100780000056
step 3.11 copper per turn per inch copper wire
Figure BDA0002804100780000057
Step 3.12, axially cooled copper surface heat dissipation coefficient
Figure BDA0002804100780000061
Step 3.13, temperature rise on the surface of the axial vent hole
Figure BDA0002804100780000062
Step 3.14 temperature rise of gas in straight-line segment
Figure BDA0002804100780000063
Average temperature of straight line segment of coil
Figure BDA0002804100780000064
Figure BDA0002804100780000065
The deviation is less than 0.0001, if not, T in 3.10 needs to be adjustedcoil_avg_initialA value of (d);
step 3.15, the copper consumption W of each circle of copper wire in one air outlet pitch at the air outletot＝If 2·Rov；
Step 3.16, gas velocity at outlet
Figure BDA0002804100780000066
Step 3.17 Reynolds number R of outlet gasot＝P·Vov·Dov；
Step 3.18, coefficient of heat dissipation of copper surface at outlet
Figure BDA0002804100780000067
Step 3.19, surface temperature rise at the outlet
Figure BDA0002804100780000068
Step 3.20, estimating the surface wind speed of the axial wind hole at the outlet
Figure BDA0002804100780000069
Step 3.21, Reynolds number R of surface wind of axial wind hole at outletev_o_ax＝P·Vo_ax·Dv；
Step 3.22, axial air hole surface heat dissipation coefficient at outlet
Figure BDA00028041007800000610
Step 3.23, exit surface temperature rise (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure BDA00028041007800000611
The position where the hot spot temperature exists is the central line of the rotor body, radial air outlet holes are formed in both sides of the position, and the radial air outlet channels of the holes 1 share the same pitch temperature;
in the pitch, the axial air hole speed of the hole-free side of the radial air duct is Hax, and the axial air hole speed of the hole-close side 2 is Hax _ o;
step 3.24, taking away copper loss of hole 1 axial air hole
Figure BDA0002804100780000071
Step 3.25 air temperature rise of axial air hole of hole 1
Figure BDA0002804100780000072
Step 3.26, surface temperature rise at the hole 1 (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure BDA0002804100780000073
Step 3.27 temperature T of well 1Root _ hole _1＝θsot _ plux _ ax _ Aperture _1+θg _ Aperture _1+Tamb；
The temperature of the radial air outlet duct at the hole 2 is within a same pitch at two sides, and in the pitch, the axial air hole speed of one side of the radial air outlet duct is Hax, and the axial air hole speed of the other side of the radial air outlet duct is Hax;
step 3.28, taking away copper loss of axial air hole of hole 2
Figure BDA0002804100780000074
Step 3.29, wind temperature rise of axial wind hole of hole 2
Figure BDA0002804100780000075
Step 3.30, surface temperature rise at the hole 2 (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure BDA0002804100780000076
Step 3.31, temperature T at hole 2Root _ hole _2＝θg _ Aperture _2+θsot _ plux _ ax _ Aperture _2+Tamb；
Step 3.32 average temperature on both sides of the well 1
Figure BDA0002804100780000077
Step 3.33 average temperature on both sides of well 2
Figure BDA0002804100780000078
Figure BDA0002804100780000081
The deviation is less than 0.0001, if not, T in 3.5 is adjustedf_out_initialThe value is obtained.
Preferably, this is achieved using MATHCAD.
Compared with the prior art, the invention has the beneficial effects that:
the calculation process is modularly designed according to functions, the whole calculation process is divided into an input data module, an air path iterative calculation module and a loss and temperature rise calculation module, the calculation process is clear, the calculation program is easy to use and upgrade, the pertinence is strong and the calculation accuracy is high aiming at the temperature rise calculation of the double-hole copper bar structure.
The method aims at the ventilation calculation of the structure of the double-hole copper bar, the calculation of the resistance and the loss of the copper bar in the air outlet area and the calculation of the surface wind speed of the axial wind hole at the air outlet, simultaneously considers the surface heat dissipation of the radial wind hole and the surface temperature rise calculation of the outlet of the surface heat dissipation of the axial wind hole, is suitable for the temperature rise design of the existing two-pole 800 MW-level to 1260 MW-level generator, can be used for developing generators adopting axial double-hole copper bar rotor coils in more capacity ranges, can meet the requirement of a new rotor hot point temperature calculation method which is higher in application degree and can ensure the precision in the optimization and improvement of products, and meets the calculation requirement of the generator adopting the axial double-hole.
Drawings
FIG. 1 is a schematic view of a typical two-hole copper bar rotor coil;
FIG. 2 is a schematic view of typical double-hole copper bar rotor coil air outlet;
FIG. 3 is a cross-sectional schematic view of a typical double-hole copper bar rotor coil copper wire;
FIG. 4 is a schematic view of an exemplary dual hole copper bar rotor coil outlet;
FIG. 5 is a schematic diagram of a hot spot position of a rotor coil of a double-hole copper bar rotor;
Detailed Description
In order to make the invention more comprehensible, preferred embodiments are described in detail below with reference to the accompanying drawings.
According to the method for calculating the hot spot temperature of the axial double-hole copper bar rotor coil of the generator, the modular design is carried out according to functions, the whole calculation process is divided into an input data module, an air path iteration calculation module and a loss and temperature rise calculation module, and the method has the advantages of clear functional modules and easy use and upgrade of calculation programs. As can be seen from figure 1, the copper bars on the two sides are designed in a rotational symmetry mode, so that one side can be taken for calculation, then the two sides are subjected to synthesis calculation, and the hot point temperature is located in the No. 7 coil. The calculation is carried out in English units with length in.
The method comprises the following specific steps:
step 1, inputting calculation parameters and calculating the structure size:
the calculation parameters comprise:
number of one pole coil: n _ coil ═ 7;
number of turns per slot of rotor: NRS ═ 7;
the rotor body is long: lp is 7830/25.4;
the arc part of the No. 1 coil is long from the body: l iscoil1_to_body＝80/25.4；
Coil end arc portion width: b iscoil_end＝41/25.4；
Coil pitch: l iscoil_to_coil＝20/25.4；
The distance between the near body side air inlet hole and the body is as follows: l isinvent_to_body＝42.5/25.4；
Distance between the near-end-portion-side air inlet hole and the end coil: l isinvent_to_endturn＝37.5/25.4；
Air outlet hole spacing: l isoutvent_to_outvent＝40/25.4；
The axial vent hole is high: h isv＝10/25.4；
Axial vent hole width: bv＝12.5/25.4；
The copper wire is high: h isn＝15/25.4；
Half the width of the double-hole copper wire: bn＝24.4/25.4；
Copper line outer chamfer radius: r1＝1.8/25.4；
Copper line inside chamfer radius: r3＝2.5/25.4；
Wall thickness in the copper wire width direction: thick _ cir is 3/25.4;
circumferential width of the air outlet hole: bot＝13.5/25.4；
Axial length of the air outlet: lot＝12/25.4；
Air outlet hole center distance: t is tot＝40/25.4；
Circumferential width of the air inlet hole: bin＝bv；
Axial length of air inlet hole: lin＝hv；
Rated rotating speed of the generator: RPM is 3000;
depth of rotor groove: h isrs＝6.102；
Rotor outer diameter: dr＝49.21；
Pressure difference between the air inlet and the air outlet of the rotor coil: pENTR＝175；
Half value of generator rated excitation: i isf6474/2, the rotor is a double-row hole, and half of the rotor coil is taken in the embodiment;
cold air temperature: t isamb＝40；
Area of vent hole on rotor wedge under backing strip: a. the12＝17*17.8/25.42；
Area of rotor slot wedge exhaust vent: a. the13＝π162/4*25.42；
Cooling gas pressure: p ═ (79.77+ 14.7)/14.7;
calculating the length of the linear end distance of the maximum coil from the body:
Lcoilmax_to_body＝Lcoil_to_body+(N_coil-1)*(Bcoil_end+Lcoil_to_coil)；
Lcoilmax_to_body*25.4＝446；
wall thickness in the copper wire height direction: WEB ═ h (h)n-hv)/2；WEB*25.4＝2.5；
The height of the axial air hole side wall is protruded in the circumferential direction of the upper wall and the lower wall of the copper wire after the radial air outlet is formed:
bdiff＝thick_cir+bv-bot；bdiff*25.4＝2；
step 2, aiming at the ventilation structure of the double-hole copper bar, calculating the wind speed of the rotor in the ventilation hole, performing control of a circulation strategy by using the MATHCAD, and judging whether the calculation precision requirement is met (see steps 2.15.3 and 2.15.4), wherein the method specifically comprises the following steps:
step 2.1, average wind path length
Figure BDA0002804100780000101
Figure BDA0002804100780000102
Step 2.2, the sectional area A of the air duct of the rotor copper wirev＝hv·bv-(4-π)R3 2＝0.185；
Wetted perimeter S of axial wind holev＝2·(hv+bv)-R3·(8-2π)＝1.603；
Hydraulic diameter of axial wind hole
Figure BDA0002804100780000103
Step 2.3, air duct inlet area Ain＝Av＝0.185；
Step 2.4, air duct outlet area Aov＝(lot·bot)＝0.251；
Wetted perimeter S of radial wind holesov＝2(bot+lot)-lot＝1.535；
Hydraulic diameter of radial wind hole
Figure BDA0002804100780000104
Step 2.5, area A of vent hole on filler strip at top turn conductor11＝Aov＝0.251；
Step 2.6, air pressure drop coefficient of air duct inlet
Figure BDA0002804100780000111
Step 2.7, wind pressure drop coefficient K at the secondary turningbe＝2.584；
Step 2.8, the air pressure drop coefficient of the axial air duct
Figure BDA0002804100780000112
frIteratively determined by the calculation of step 2.15.3, fr＝0.00526；
Step 2.9, carve the wind pressure drop coefficient of the bottom filler strip
Figure BDA0002804100780000113
Step 2.10, wind pressure drop coefficient of air outlet of slot wedge
Figure BDA0002804100780000114
Step 2.11, considering resistance coefficient increased by orifice plate effect of air outlet
Figure BDA0002804100780000115
Step 2.12, total wind pressure drop coefficient of rotor wind path
Kv＝KOutlet _ Orifice plate Effect _ Orifice _4+Kfr+(Kin+Kbe+K12+K13)＝11.832；
Step 2.13 wind speed of rotor surface
Figure BDA0002804100780000116
Wind speed at rotor slot bottom
Figure BDA0002804100780000117
Centrifugal wind pressure
Figure BDA0002804100780000118
Step 2.14, total pressure drop Δ P in the air path is equal to PENTR + Pcen＝212.504；
Step 2.15, a wind speed and wind amount calculation part comprises the following steps
Step 2.15.1, setting the wind speed of the axial vent hole as V1＝17630；
Figure BDA0002804100780000119
Step 2.15.2 Reynolds number Rev1＝P·V1·Dv＝5.244×104；
Step 2.15.3, coefficient of friction
Figure BDA00028041007800001110
Figure BDA0002804100780000121
The deviation is less than 0.0001, if not, V in formula 2.15.1 needs to be adjusted1A value;
total wind pressure drop of step 2.15.4 wind path
Figure BDA0002804100780000122
Figure BDA0002804100780000123
The deviation is less than 0.0001, if not, V in formula 2.15.1 needs to be adjusted1A value;
step 2.15.5, air volume of each circle of copper wire
Figure BDA0002804100780000124
Step 3, loss and temperature rise calculation: calculating the loss of each circle of copper wire, the loss of the copper wire in unit length, the surface temperature rise of the axial vent hole and the gas temperature rise; and aiming at the air outlet area where the hot spot temperature is located, the gas speed and the surface temperature rise at the outlet are respectively calculated. Aiming at the structure of the outlet of the double-hole copper bar, the wind speed, the heat dissipation coefficient, the loss and the like of the axial wind hole and the radial wind port at different positions of the outlet are specially analyzed, and the heat conduction at two sides of the double-hole copper bar is considered. The method specifically comprises the following steps:
step 3.1, take away long loss of copper
Figure BDA0002804100780000125
Step 3.2 axial conductive sectional area A of copper wirecu＝bn·hn-(4-π)R1 2-Av＝0.378；
Step 3.3 No radial hole section conductive section area A at air outletcuh1＝Acu＝0.378；
Step 3.4, the air outlet part is provided with a radial hole section conductive sectional area
Acuh2＝Acu-2WEB·(bot-thick_cir)-thick_cir·hn＝0.226；
Step 3.5, an air hole pitch internal resistance of a radial hole section is not arranged at the air outlet part:
temperature at rotor coil outlet is Tf_out_initial＝114.705，Tf_out_initialIteratively determined from the results of the 3.33 calculations,
Figure BDA0002804100780000126
step 3.6, an air outlet part is provided with an inner resistor with a radial hole section and an air hole pitch
Figure BDA0002804100780000127
Step 3.7 air outlet one air hole pitch internal resistance Rov＝Rot1+Rot2＝4.66×10-6；
Step 3.8 equivalent area of copper wire in air-out area
Figure BDA0002804100780000131
Step 3.9, multiplying power reduction of equivalent area of copper wire in air outlet area
Figure BDA0002804100780000132
Step 3.10, copper consumption of each circle of copper wire:
average temperature T of straight section of rotor coilcoil_avg_initial76.949 degrees Celsius, Tcoil_avg_initialIteratively determined from the results of the 3.14 calculations,
Figure BDA0002804100780000133
step 3.11 copper per turn per inch copper wire
Figure BDA0002804100780000134
Step 3.12, axially cooled copper surface heat dissipation coefficient
Figure BDA0002804100780000135
Step 3.13, temperature rise on the surface of the axial vent hole
Figure BDA0002804100780000136
Step 3.14,Gas temperature rise of straight line segment
Figure BDA0002804100780000137
Average temperature of straight line segment of coil
Figure BDA0002804100780000138
Figure BDA0002804100780000139
Deviation of requirement is less than 0.0001, if not, T in 3.10 needs to be adjustedcoil_avg_initialA value;
step 3.15, the copper consumption W of each circle of copper wire in one air outlet pitch at the air outletot＝If 2·Rov＝48.828；
Step 3.16, gas velocity at outlet
Figure BDA0002804100780000141
Step 3.17 Reynolds number R of outlet gasot＝P·Vov·Dov＝5.473×104；
Step 3.18, coefficient of heat dissipation of copper surface at outlet
Figure BDA0002804100780000142
Step 3.19, surface temperature rise at the outlet
Figure BDA0002804100780000143
Step 3.20, estimating the surface wind speed of the axial wind hole at the outlet
Figure BDA0002804100780000144
Step 3.21, Reynolds number R of surface wind of axial wind hole at outletev_o_ax＝P·Vo_ax·Dv＝1.732×104；
Step 3.22, axial wind at the outletCoefficient of heat dissipation of pore surface
Figure BDA0002804100780000145
Step 3.23, exit surface temperature rise (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure BDA0002804100780000146
The position where the hot spot temperature exists is the central line of the rotor body, radial air outlet holes are formed in both sides of the position, and the radial air outlet channels of the holes 1 share the same pitch temperature;
in the pitch, the axial air hole speed of the hole-free side of the radial air duct is Hax, and the axial air hole speed of the hole-close side 2 is Hax _ o;
step 3.24, taking away copper loss of hole 1 axial air hole
Figure BDA0002804100780000147
Step 3.25 air temperature rise of axial air hole of hole 1
Figure BDA0002804100780000148
Step 3.26, surface temperature rise at the hole 1 (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure BDA0002804100780000149
Step 3.27 temperature T of well 1Root _ hole _1＝θsot _ plux _ ax _ Aperture _1+θg _ Aperture _1+Tamb＝110.651；
The temperature of the radial air outlet duct at the hole 2 is within a same pitch at two sides, and in the pitch, the axial air hole speed of one side of the radial air outlet duct is Hax, and the axial air hole speed of the other side of the radial air outlet duct is Hax;
step 3.28, taking away copper loss of axial air hole of hole 2
Figure BDA0002804100780000151
Step 3.29, wind temperature rise of axial wind hole of hole 2
Figure BDA0002804100780000152
Step 3.30, surface temperature rise at the hole 2 (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure BDA0002804100780000153
Step 3.31, temperature T at hole 2Root _ hole _2＝θg _ Aperture _2+θsot _ plux _ ax _ Aperture _2+Tamb＝118.755；
Step 3.32 average temperature on both sides of the well 1
Figure BDA0002804100780000154
Step 3.33 average temperature on both sides of well 2
Figure BDA0002804100780000155
Figure BDA0002804100780000156
The deviation is less than 0.0001, if not, T in 3.5 is adjustedf_out_initialThe value is obtained.

Claims (2)

1. A method for calculating the hot spot temperature of a generator axial double-hole copper bar rotor coil is characterized by comprising the following steps:
step 1, inputting calculation parameters and calculating the structure size:
the calculation parameters comprise:
number of one pole coil: n _ coil;
NRS is the number of turns of each slot of the rotor;
the rotor body is long: lp;
the arc part of the No. 1 coil is long from the body: l iscoil1_to_body；
Coil end arc portion width: b iscoil_end；
Coil pitch: l iscoil_to_coil；
The distance between the near body side air inlet hole and the body is as follows: l isinvent_to_body；
Distance between the near-end-portion-side air inlet hole and the end coil: l isinvent_to_endturn；
Air outlet hole spacing: l isoutvent_to_outvent；
The axial vent hole is high: h isv；
Axial vent hole width: bv；
The copper wire is high: h isn；
Half the width of the double-hole copper wire: bn；
Copper line outer chamfer radius: r1；
Copper line inside chamfer radius: r3；
Wall thickness in the copper wire width direction: thick _ cir;
circumferential width of the air outlet hole: bot；
Axial length of the air outlet: lot；
Air outlet hole center distance: t is tot；
Circumferential width of the air inlet hole: bin；
Axial length of air inlet hole: lin；
Rated rotating speed of the generator: RPM;
depth of rotor groove: h isrs；
Rotor outer diameter: dr；
Pressure difference between the air inlet and the air outlet of the rotor coil: pENTR；
Half value of generator rated excitation: i isf；
Cold air temperature: t isamb；
Area of vent hole on rotor wedge under backing strip: a. the12；
Area of rotor slot wedge exhaust vent: a. the13；
Cooling gas pressure: p;
calculating the length of the linear end distance of the maximum coil from the body:
Lcoilmax_to_body＝Lcoil_to_body+(N_coil-1)*(Bcoil_end+Lcoil_to_coil)；
wall thickness in the copper wire height direction: WEB ═ h (h)n-hv)/2；
The height of the axial air hole side wall is protruded in the circumferential direction of the upper wall and the lower wall of the copper wire after the radial air outlet is formed:
bdiff＝thick_cir+bv-bot；
step 2, to the ventilation structure of diplopore copper bar, calculate the rotor wind speed size in the ventilation hole, judge whether satisfy the calculation accuracy requirement, specifically include:
step 2.1, average wind path length
Figure FDA0002804100770000021
Figure FDA0002804100770000022
Step 2.2, the sectional area A of the air duct of the rotor copper wirev＝hv·bv-(4-π)R3 2；
Wetted perimeter S of axial wind holev＝2·(hv+bv)-R3·(8-2π)；
Hydraulic diameter of axial wind hole
Figure FDA0002804100770000023
Step 2.3, air duct inlet area Ain＝Av；
Step 2.4, air duct outlet area Aov＝(lot·bot)；
Wetted perimeter S of radial wind holesov＝2(bot+lot)-lot；
Hydraulic diameter of radial wind hole
Figure FDA0002804100770000024
Step 2.5, area A of vent hole on filler strip at top turn conductor11＝Aov；
Step 2.6, air pressure drop coefficient of air duct inlet
Figure FDA0002804100770000025
Step 2.7, wind pressure drop coefficient K at the secondary turningbe＝2.584；
Step 2.8, the air pressure drop coefficient of the axial air duct
Figure FDA0002804100770000031
frIteratively determined by the results of the calculation of step 2.15.3;
step 2.9, carve the wind pressure drop coefficient of the bottom filler strip
Figure FDA0002804100770000032
Step 2.10, wind pressure drop coefficient of air outlet of slot wedge
Figure FDA0002804100770000033
Step 2.11, considering resistance coefficient increased by orifice plate effect of air outlet
Figure FDA0002804100770000034
Step 2.12, total wind pressure drop coefficient of rotor wind path
Kv＝KOutlet _ Orifice plate Effect _ Orifice _4+Kfr+(Kin+Kbe+K12+K13)；
Step 2.13 wind speed of rotor surface
Figure FDA0002804100770000035
Wind speed at rotor slot bottom
Figure FDA0002804100770000036
Centrifugal wind pressure
Figure FDA0002804100770000037
Step 2.14, total pressure drop Δ P in the air path is equal to PENTR + Pcen；
Step 2.15, a wind speed and wind amount calculation part comprises the following steps
Step 2.15.1, setting the wind speed of the axial vent hole as V1；
Figure FDA0002804100770000038
Step 2.15.2 Reynolds number Rev1＝P·V1·Dv；
Step 2.15.3, coefficient of friction
Figure FDA0002804100770000039
Figure FDA00028041007700000310
The deviation is less than 0.0001, if not, V in formula 2.15.1 needs to be adjusted1A value;
total wind pressure drop of step 2.15.4 wind path
Figure FDA0002804100770000041
Figure FDA0002804100770000042
The deviation is less than 0.0001, if not, V in formula 2.15.1 needs to be adjusted1A value;
step 2.15.5, air volume of each circle of copper wire
Figure FDA0002804100770000043
Step 3, loss and temperature rise calculation: calculating the loss of each circle of copper wire, the loss of the copper wire in unit length, the surface temperature rise of the axial vent hole and the gas temperature rise; and aiming at the air outlet area where the hot spot temperature is located, the gas speed and the surface temperature rise at the outlet are respectively calculated, and the method specifically comprises the following steps:
step 3.1, take away long loss of copper
Figure FDA0002804100770000044
Step 3.2 axial conductive sectional area A of copper wirecu＝bn·hn-(4-π)R1 2-Av；
Step 3.3 No radial hole section conductive section area A at air outletcuh1＝Acu；
Step 3.4, the air outlet part is provided with a radial hole section conductive sectional area
Acuh2＝Acu-2WEB·(bot-thick_cir)-thick_cir·hn；
Step 3.5, an air hole pitch internal resistance of a radial hole section is not arranged at the air outlet part:
temperature at rotor coil outlet is Tf_out_initial，Tf_out_initialIteratively determined from the results of the 3.33 calculations,
Figure FDA0002804100770000045
step 3.6, an air outlet part is provided with an inner resistor with a radial hole section and an air hole pitch
Figure FDA0002804100770000046
Step 3.7 air outlet one air hole pitch internal resistance Rov＝Rot1+Rot2；
Step 3.8 equivalent area of copper wire in air-out area
Figure FDA0002804100770000047
Step 3.9, multiplying power reduction of equivalent area of copper wire in air outlet area
Figure FDA0002804100770000048
Step 3.10, copper consumption of each circle of copper wire:
average temperature T of straight section of rotor coilcoil_avg_initial，Tcoil_avg_initialIteratively determined from the results of the 3.14 calculations,
Figure FDA0002804100770000051
step 3.11 copper per turn per inch copper wire
Figure FDA0002804100770000052
Step 3.12, axially cooled copper surface heat dissipation coefficient
Figure FDA0002804100770000053
Step 3.13, temperature rise on the surface of the axial vent hole
Figure FDA0002804100770000054
Step 3.14 temperature rise of gas in straight-line segment
Figure FDA0002804100770000055
Average temperature of straight line segment of coil
Figure FDA0002804100770000056
Figure FDA0002804100770000057
The deviation is less than 0.0001, if not, T in 3.10 needs to be adjustedcoil_avg_initialA value of (d);
step 3.15, the copper consumption W of each circle of copper wire in one air outlet pitch at the air outletot＝If 2·Rov；
Step 3.16, gas velocity at outlet
Figure FDA0002804100770000058
Step 3.17 Reynolds number R of outlet gasot＝P·Vov·Dov；
Step 3.18, coefficient of heat dissipation of copper surface at outlet
Figure FDA0002804100770000059
Step 3.19, surface temperature rise at the outlet
Figure FDA00028041007700000510
Step 3.20, estimating the surface wind speed of the axial wind hole at the outlet
Figure FDA00028041007700000511
Step 3.21, Reynolds number R of surface wind of axial wind hole at outletev_o_ax＝P·Vo_ax·Dv；
Step 3.22, axial air hole surface heat dissipation coefficient at outlet
Figure FDA0002804100770000061
Step 3.23, exit surface temperature rise (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure FDA0002804100770000062
The position where the hot spot temperature exists is the central line of the rotor body, radial air outlet holes are formed in both sides of the position, and the radial air outlet channels of the holes 1 share the same pitch temperature;
in the pitch, the axial air hole speed of the hole-free side of the radial air duct is Hax, and the axial air hole speed of the hole-close side 2 is Hax _ o;
step 3.24, taking away copper loss of hole 1 axial air hole
Figure FDA0002804100770000063
Step 3.25 air temperature rise of axial air hole of hole 1
Figure FDA0002804100770000064
Step 3.26, surface temperature rise at the hole 1 (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure FDA0002804100770000065
Step 3.27 temperature T of well 1Root _ hole _1＝θsot _ plux _ ax _ Aperture _1+θg _ Aperture _1+Tamb；
The temperature of the radial air outlet duct at the hole 2 is within a same pitch at two sides, and in the pitch, the axial air hole speed of one side of the radial air outlet duct is Hax, and the axial air hole speed of the other side of the radial air outlet duct is Hax;
step 3.28, taking away copper loss of axial air hole of hole 2
Figure FDA0002804100770000066
Step 3.29, wind temperature rise of axial wind hole of hole 2
Figure FDA0002804100770000067
Step 3.30, surface temperature rise at the hole 2 (considering radial wind hole surface heat radiation and axial wind hole surface heat radiation)
Figure FDA0002804100770000068
Step 3.31, temperature T at hole 2Root _ hole _2＝θg _ Aperture _2+θsot _ plux _ ax _ Aperture _2+Tamb；
Step 3.32 average temperature on both sides of the well 1
Figure FDA0002804100770000071
Step 3.33 average temperature on both sides of well 2
Figure FDA0002804100770000072
Figure FDA0002804100770000073
The deviation is less than 0.0001, if not, T in 3.5 is adjustedf_out_initialThe value is obtained.
2. The method for calculating the hot spot temperature of the generator axial double-hole copper bar rotor coil as claimed in claim 1, wherein the method comprises the following steps: it is realized by MATHCAD.
CN202011361490.XA 2020-11-27 2020-11-27 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator Active CN112560367B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011361490.XA CN112560367B (en) 2020-11-27 2020-11-27 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011361490.XA CN112560367B (en) 2020-11-27 2020-11-27 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator

Publications (2)

Publication Number Publication Date
CN112560367A true CN112560367A (en) 2021-03-26
CN112560367B CN112560367B (en) 2022-05-27

Family

ID=75046470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011361490.XA Active CN112560367B (en) 2020-11-27 2020-11-27 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator

Country Status (1)

Country Link
CN (1) CN112560367B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113593752A (en) * 2021-08-07 2021-11-02 济南宝世达新材料有限公司 Hollow conductor and processing technology thereof

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462624A (en) * 2013-09-16 2015-03-25 同济大学 Motor temperature field data processing method based on multi-field coupling
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
CN105450122A (en) * 2016-01-14 2016-03-30 重庆大学 IGBT device junction temperature fluctuation inhibition method of doubly-fed wind turbine generator system machine-side current transformer
CN106777459A (en) * 2016-11-10 2017-05-31 北京交通大学 The computational methods in generator shaft radial rotor temperature field
US20190173414A1 (en) * 2018-01-29 2019-06-06 Yangzhou University Prediction Method for Reliability Degree of Running Temperature Rise of a Large and Medium-sized Motor
CN110688738A (en) * 2019-09-10 2020-01-14 北京交通大学 Method for calculating temperature field of nuclear power turbonator with axial ventilation system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462624A (en) * 2013-09-16 2015-03-25 同济大学 Motor temperature field data processing method based on multi-field coupling
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
CN105450122A (en) * 2016-01-14 2016-03-30 重庆大学 IGBT device junction temperature fluctuation inhibition method of doubly-fed wind turbine generator system machine-side current transformer
CN106777459A (en) * 2016-11-10 2017-05-31 北京交通大学 The computational methods in generator shaft radial rotor temperature field
US20190173414A1 (en) * 2018-01-29 2019-06-06 Yangzhou University Prediction Method for Reliability Degree of Running Temperature Rise of a Large and Medium-sized Motor
CN110688738A (en) * 2019-09-10 2020-01-14 北京交通大学 Method for calculating temperature field of nuclear power turbonator with axial ventilation system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113593752A (en) * 2021-08-07 2021-11-02 济南宝世达新材料有限公司 Hollow conductor and processing technology thereof
CN113593752B (en) * 2021-08-07 2023-01-13 孙福学 Hollow conductor and processing technology thereof

Also Published As

Publication number Publication date
CN112560367B (en) 2022-05-27

Similar Documents

Publication Publication Date Title
EP2676000B1 (en) Integrated axial and tangential serpentine cooling circuit in a turbine airfoil
US8864438B1 (en) Flow control insert in cooling passage for turbine vane
WO2016201877A1 (en) Radial ventilation cooling structure for motor
JP4486114B2 (en) Rotating electric machine
US9188016B2 (en) Multi-orifice plate for cooling flow control in vane cooling passage
US8827632B1 (en) Integrated TBC and cooling flow metering plate in turbine vane
CN112560367B (en) Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator
JP6570734B2 (en) Cooling system for electric machine
RU2619664C2 (en) Cooled composite sheets for gas turbine
EP2944762B1 (en) Airfoil with improved cooling
CN216146188U (en) System of electric motor cooling jackets and system of coaxially surrounding electric motor cooling jackets
CN104135114A (en) High-voltage air and water cooling doubly-fed wind driven generator
Zhou et al. Ventilation cooling design for a novel 350-MW air-cooled turbo generator
WO2023169063A1 (en) Outer rotor permanent magnet roller stator and stator assembly thereof
CN104539071B (en) Ventilation channel steel, manufacturing method thereof, ventilation structure and motor
JP6272525B1 (en) Rotating electric machine
KR101719061B1 (en) Ventilative channel steel as well as manufacturing method, ventilating structure and motor thereof
WO2015107419A1 (en) A dendritic tube circular fin heat exchanger
CN207868915U (en) A kind of fluid-cooled electrical machine shell
CN111983404B (en) Hotspot calculation and positioning method of lightweight vehicle-mounted traction transformer
CN113890234B (en) Closed motor cooling device with circulation convection between rotor holes
Liu et al. Thermal analysis and cooling structure design of axial flux permanent magnet synchronous motor for electrical vehicle
JP5978800B2 (en) Rotating electric machine and cooling method of rotating electric machine stator
CN216873021U (en) Structure for reducing temperature rise of submersible motor rotor by adopting oblique heat conducting grooves
CN111578750B (en) Secondary flow heat exchange device and system

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
GR01 Patent grant
GR01 Patent grant