CN104462624A - Motor temperature field data processing method based on multi-field coupling - Google Patents

Motor temperature field data processing method based on multi-field coupling Download PDF

Info

Publication number
CN104462624A
CN104462624A CN201310422756.0A CN201310422756A CN104462624A CN 104462624 A CN104462624 A CN 104462624A CN 201310422756 A CN201310422756 A CN 201310422756A CN 104462624 A CN104462624 A CN 104462624A
Authority
CN
China
Prior art keywords
motor
formula
loss
rotor
data processing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201310422756.0A
Other languages
Chinese (zh)
Inventor
马淑梅
刘雪梅
李爱平
宋守泰
谢楠
徐立云
陈杨杨
倪泽浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN201310422756.0A priority Critical patent/CN104462624A/en
Publication of CN104462624A publication Critical patent/CN104462624A/en
Pending legal-status Critical Current

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
    • Y02T10/82—Elements for improving aerodynamics

Landscapes

  • Iron Core Of Rotating Electric Machines (AREA)

Abstract

The invention relates to a motor temperature field data processing method based on multi-field coupling. The motor temperature field data processing method includes 1) establishing a two-dimensional electromagnetic field model and a three-dimensional fluid field model; 2) calculating electromagnetic field distribution of a motor for the two-dimensional electromagnetic field model, computing stator coil copper loss, rotor bar copper loss and stator and rotor core loss; 3) for the three-dimensional fluid field model, selecting a computing area to simplify the three-dimensional fluid field model and obtain the fluid field distribution of the motor; 4) computing heat radiation coefficient of the surface of the motor; 5) loading various loss values during motor operating into the three-dimensional fluid field model as a heat source to obtain various boundary conditions according to the heat radiation coefficient, and adopting the finite volume method to obtain the temperature field distribution of the motor. Compared with the prior art, the motor temperature field data processing method has the advantages that influence of the electromagnetic field on the temperature field is considered, the fluid field and the temperature field are coupled, and consideration is more complete, so that computing of the motor temperature field is higher in accuracy and analysis efficiency.

Description

A kind of electric motor temperature field data processing method based on multi-scenarios method
Technical field
The present invention relates to simulation calculation field, especially relate to a kind of electric motor temperature field data processing method based on multi-scenarios method.
Background technology
Motor is the important bridge of being come from the conversion of initial energy form by electric energy, is the device most of electric energy being converted to mechanical energy.It is simple that asynchronous machine has structure, manufactures, operation and maintenance is convenient, reliable, and the advantage such as lighter in weight, cost be lower.Asynchronous machine has higher operational efficiency and good operating characteristic, close to constant-speed operation in from zero load to full-load range, can meet the transmission requirement of most of industrial and agricultural production machinery.Along with the development of electric manufacturing, the improvement of the type of cooling, material and progress, the single-machine capacity of motor is in continuous increase, motor electromagnetic load and thermal load also improve thereupon, the unit volume loss produced when motor runs constantly increases, it is serious to make that the motor feels hot, and this will directly have influence on serviceable life and the operational reliability of motor.Large-size asynchronous motor is due to complex structure, and cost is high, usually all can cause huge economic loss after breaking down.
Also be an importance in motor research to the research of the heat problem of motor always, heat problem directly affects serviceable life and the reliability of operation of motor, such as, when stator winding is short-circuited fault, it is overheated that the short-circuit current that numerical value is very large can produce, and burns winding and iron core.For most of unusual service condition, the principal element finally constituted a threat to motor safety is overheated, this is determined by following two aspects: first, temperature rise is the key factor affecting insulation life, and the temperature of large-size machine regional area and parts may reach very high level.Secondly, because the intensity of metal material in motor and hardness decline with the rising of temperature, in temperature too high or temperature is not too high but the duration is long time all can cause the decline of metal strength, thus have influence on the safe operation of motor.And large-size machine due to single-machine capacity large, motor load also improves thereupon, and the unit volume loss produced when motor runs is larger, makes large-size machine fever phenomenon more serious.Therefore, calculating exactly and study large-size asynchronous motor temperature rise is manufacture and the department of operation questions of common interest.
Traditional electric machine temperature rise computing method have formula of reduction method, Equivalent heat path method.Formula of reduction method is that the whole core loss of supposition and copper loss are only shed by stator (or rotor) cylindrical cooling surface, does not have heat interchange between the live part of armature winding copper and contact portion.This method advantage calculates simple, and shortcoming is accurate, sufficiently complete not, approximate treatment can only go out the average temperature rising of motor.Equivalent heat path method forms Equivalent heat path according to thermal conduction study and Circuit theory, and the thermal source in thermal circuits is the copper loss of winding, iron loss, and loss heat, by various corresponding thermal resistance, by thermal source to heat eliminating medium transmission, forms a complicated ther mal network.Some theorems in circuit network and method is adopted to calculate the average temperature rising of each live part of motor.The method computational accuracy is higher than formula of reduction, can obtain the overall temperature rise of motor and average temperature rising.But also have and can not calculate local temperature rise and the large shortcoming of workload.And use the temperature field method of modern numerical method to be the temperature rise calculation method risen recently, namely domain is separated into many junior units, sets up heat-conduction equation in each cell, more overall system of equations is solved.Temperature field method transfers research object to microcosmic from macroscopic view, from totally forwarding local unit to, tries to achieve temperature and the temperature rise of every bit.
Separately motor stator temperature field or rotor temperature field are solved being also mainly limited in the Calculation Method of Temperature Field of motor at present, and also less to the research work of motor stator and rotor temperature field simultaneous solution.Some model only considers stator temperature field, and some model only considers rotor temperature field, so just have ignored the conduction of heat of air-gap between rotor.In addition, temperature field and fluid field and electromagnetic field in close relations, some researchs are now mainly for the coupling in fluid field and temperature field, and the coupling in electromagnetic field and temperature field, not to the research that many, motor considers.
Summary of the invention
Object of the present invention be exactly in order to overcome above-mentioned prior art exist defect and a kind of electric motor temperature field data processing method based on multi-scenarios method is provided.
Object of the present invention can be achieved through the following technical solutions: a kind of electric motor temperature field data processing method based on multi-scenarios method, is characterized in that, comprise the following steps:
1) use finite element software to set up physical model of electrical machine, comprise and set up two-dimensional electromagnetic field model according to Theory of Electromagnetic Field and set up three dimensional fluid field model according to Hydrodynamics Theory;
2) for two-dimensional electromagnetic field model, Finite Element Method is utilized to set up the field circuit method finite element equation solved in territory, calculate the magnetic distribution of motor, comprise magnetic flux distribution and distribution of current, and calculate stator coil copper loss, rotor bar copper loss and stator and rotor cores loss;
3) for three dimensional fluid field model, choose zoning and simplify three dimensional fluid field model, model after simplification comprises air gap flow field and Temperature calculating physical model and stator-rotor iron core, winding solid temperature field computational physics model between rotor, three dimensional fluid field model after simplification is imported in finite element software and carries out stress and strain model, adopt finite volume method to carry out iterative computation to zoning, obtain the fluid field distribution of motor;
4) according to the wind speed everywhere in ventilating system and electric machine rotor structural parameters, the coefficient of heat transfer on motor surface is calculated by hydrodynamics method formula and law of conservation of mass;
5) using step 2) the every loss value of the motor that draws of electromagnetic field analysis when running be loaded in three dimensional fluid field model as thermal source, according to step 4) coefficient of heat transfer that calculates draws all kinds of boundary condition, adopt finite volume method to carry out Coupled Numerical Simulation to the zoning in electric machine rotor temperature field under different operating mode, finally obtain the thermo parameters method of motor.
Step 1) described in two-dimensional electromagnetic field model be motor shaft to two-dimensional section.
Step 2) described in calculating stator coil copper loss P cucomputing formula be
P cu=∑I x 2R x，
In formula, I xfor the electric current in stator winding x, R xduring for being converted to benchmark job temperature, the resistance of winding x, carries out aftertreatment and obtains coil of stator of motor copper loss in finite element software.
Step 2) described in calculating rotor bar copper loss P ' cucomputing formula be
P′ cu＝∑I y 2R y，
In formula, I yfor the electric current in rotor windings y, R yduring for being converted to benchmark job temperature, the resistance of winding y, carries out aftertreatment and obtains the loss of rotor conductive copper in finite element software.
Step 2) described in stator and rotor cores loss comprise the basic loss of stator core and stator core added losses P fep, the basic loss of described stator core comprises yoke portion basic loss P unshakable in one's determination fejbasic loss P unshakable in one's determination with teeth portion fet.
Described calculating yoke portion basic loss P unshakable in one's determination fejexpression formula be P fej=K αp fejg j,
In formula, G jfor yoke portion clean iron unshakable in one's determination; K αfor experience factor, for alternating current generator, as capacity P nduring <100kVA, K α=1.5, as capacity P nduring>=100kVA, K α=1.3, for direct current generator, K α=3.6; Specific loss p fejcomputing formula be
p Fej = p 10 / 50 B j 2 ( f f 0 ) 1.3 ,
In formula, B jfor peakflux density value in yoke, p 10/50for working as B=1T, during F=50Hz, the siliconized plate loss in unit weight;
Described calculating teeth portion basic loss P unshakable in one's determination fetexpression formula be P fet=K βp fetg t,
In formula, G tfor teeth portion clean iron unshakable in one's determination; K βfor experience factor, for asynchronous machine, as capacity P nduring <100kVA, K β=2.0, as capacity P nduring>=100kVA, K β=1.7, for direct current generator, K β=4.0; Specific loss p fetcomputing formula be
p Fet = p 10 / 50 B t 2 ( f f 0 ) 1.3 ,
In formula, B tfor the mean value of teeth portion magnetic flux density.
Described calculating stator core added losses P fepexperimental formula be
P Fep=K 0(B 0t) 2(zn) 1.5A p，
In formula, K 0for the coefficient relevant with magnetic pole material therefor, B 0for slot ripples magnetic flux density maximal value, t is armature tooth distance, and z is the number of teeth, and n is magnetic pole rotating speed, A pfor magnetic pole area.
Step 3) described in zoning concrete grammar of choosing be in electric machine rotor circumference, choose 1/12nd of rotor as zoning.
Step 4) described in the coefficient of heat transfer on motor surface specifically comprise:
A) yoke portion coefficient of heat transfer α in motor Radial ventilation duct e,
in formula, v efor the mean wind speed in yoke portion in Radial ventilation duct;
B) teeth portion coefficient of heat transfer α in motor Radial ventilation duct c,
in formula, v cfor the mean wind speed of teeth portion in Radial ventilation duct;
C) the coefficient of heat transfer α of insulation sides in air duct b,
in formula, v bfor the wind speed of insulation sides in air duct;
D) the convection coefficient α of air gap surface δ,
α δ=28 (1+ ω δ 0.5), in formula, ω δfor mean wind speed, ω δ=v/2,
V is rotor peripheral speed, wherein n nfor motor speed, D 2for rotor outer circle diameter;
E) the coefficient of heat transfer α of stator core cylindrical and rotor core inner circle w,
in formula, D hfor equivalent diameter, λ is Measurement of Gas Thermal Conductivity .N ufor nusselt number,
Nusselt number N ucalculate with following formula:
N u = 0.152 R e 0.9 P r 0.833 [ 2.25 ln ( 0.144 R e 0.9 ) + 13.2 P r - 5.8 ] ,
In formula, P rfor Prandtl number, R efor Reynolds number.
Described motor is large-scale squirrel cage asynchronous motor.
The invention provides a kind of easily across the coupling calculation between physical field, carry out motor multiple physical field coupling analysis across different software, and calculate the thermo parameters method of motor.
Compared with prior art, the inventive method applies the function of existing CAE software to greatest extent, and improve electric machine temperature rise computational accuracy, save the time, for the optimal design of large-size asynchronous motor, the enforcement of control strategy, safe operation provide theoretical foundation.
Accompanying drawing explanation
Fig. 1 is schematic flow sheet of the present invention.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.
As shown in Figure 1, a kind of electric motor temperature field data processing method based on multi-scenarios method, is characterized in that, comprise the following steps:
Step 101, uses CAD software to set up physical model of electrical machine according to motor entity M, comprises and sets up two-dimensional electromagnetic field model according to Theory of Electromagnetic Field and set up three dimensional fluid field model according to Hydrodynamics Theory;
According to Theory of Electromagnetic Field, motor shaft is equal everywhere to can be regarded as, and therefore chooses motor shaft and solves physical model to any one two-dimensional section as electromagnetic field, set up two-dimensional physical model in CAD software, carry out stress and strain model.
Step 102, for two-dimensional electromagnetic field model, corresponding load loading is carried out under different operating mode to motor, the electromagnetic field of finite element method to motor is utilized to calculate, obtain corresponding motor magnetic flux Density Distribution and distribution of current, consider the impact of end on electromagnetic field, adopt coupled field-circuit method to revise, stator coil end, rotor conductor end are taken into account;
Step 103, and calculate stator coil copper loss, rotor bar copper loss and stator and rotor cores loss;
Core loss is decided by the alternative frequency of motor-field to a great extent, and magnetic field alternative frequency is f in asynchronous machine rotor iron core 2=sf 1, the revolutional slip s of general motor is about about 1.5%, therefore the close alternative frequency of rotor core magnetic is very low, and rotor iron loss is negligible relative to stator, so only calculate stator iron loss in general design.
(1) stator coil copper loss P is calculated cucomputing formula be
P cu = &Sigma; I x 2 R x ,
In formula, I xfor the electric current in stator winding x, R xduring for being converted to benchmark job temperature, the resistance of winding x, carries out aftertreatment and obtains coil of stator of motor copper loss in finite element software.
(2) rotor bar copper loss P ' is calculated cucomputing formula be
P cu &prime; = &Sigma; I y 2 R y ,
In formula, I yfor the electric current in rotor windings y, R yduring for being converted to benchmark job temperature, the resistance of winding y, carries out aftertreatment and obtains the loss of rotor conductive copper in finite element software.
(3) the stator and rotor cores loss described in comprises the basic loss of stator core and stator core added losses P fepbasic loss is the loss that main field produces during alternation in iron core, added losses to be slotted the loss (i.e. additional idle loss) causing air-gap permeance harmonic field to cause in the other side's iron core due to rotor, and the basic loss of described stator core comprises yoke portion basic loss P unshakable in one's determination fejbasic loss P unshakable in one's determination with teeth portion fet.
The basic loss P of electric machine iron core feexpression formula be P fe=K αp feg fe, in formula, G fefor iron core uses iron only, K αfor due to siliconized plate processing, magnetic flux distribution uneven, and its not reason such as sinusoidal variations and coefficient of causing loss to increase in time, p fefor the loss of unit quality, also claim specific loss.Due in iron core, yoke portion is different with teeth portion magnetic flux distribution, and thus in iron core, the iron loss of yoke portion and teeth portion should calculate respectively, therefore calculates the basic loss of stator core and comprises loss and the teeth portion basic loss unshakable in one's determination substantially of the yoke portion iron core of calculating.
Described calculating yoke portion basic loss P unshakable in one's determination fejexpression formula be P fej=K αp fejg j, in formula, G jfor yoke portion clean iron unshakable in one's determination, K αfor experience factor, for alternating current generator, as capacity P nduring <100kVA, K α=1.5; As capacity P nduring>=100kVA, K α=1.3, for direct current generator, K α=3.6; Specific loss p fejcomputing formula be in formula, B jfor peakflux density value in yoke, p 10/50for working as B=1T, during F=50Hz, the loss in unit weight, its value can be looked into and get from siliconized plate damage curve.
Described calculating teeth portion basic loss P unshakable in one's determination fetexpression formula be P fet=K βp fetg t, in formula, G tfor teeth portion clean iron unshakable in one's determination, K βfor experience factor, for asynchronous machine, as capacity P nduring <100kVA, K β=2.0; As capacity P nduring>=100kVA, K β=1.7, for direct current generator, K β=4.0; Specific loss p fetcomputing formula be in formula, B tfor the mean value of teeth portion magnetic flux density.
Described calculating stator core added losses P fepconcrete grammar is: added losses unshakable in one's determination add iron loss also known as zero load, is to be caused by the higher hamonic wave magnetic field in air gap, and it calculates experimental formula is P fep=K 0(B 0t) 2(zn) 1.5a p, in formula, K 0for the coefficient relevant with magnetic pole material therefor, B 0for slot ripples magnetic flux density maximal value, t is armature tooth distance, and z is the number of teeth, and n is magnetic pole rotating speed, A pfor magnetic pole area.
Step 104, for three dimensional fluid field model, because electric machine rotor has symmetry in the circumferential, choose 1/12nd of electric machine rotor in the circumferential as zoning, under the prerequisite not affecting computational accuracy, consider follow-up stress and strain model and numerical simulation, three dimensional fluid field model is reasonably simplified, final mask comprises air gap flow field and Temperature calculating physical model and stator-rotor iron core between rotor, winding solid temperature field computational physics model, three dimensional fluid field model after simplification is imported in finite element software and carries out stress and strain model, finite volume method is adopted to carry out iterative computation to zoning, obtain the fluid field distribution of motor, obtain velocity distribution situation in ventilation slot and air gap,
Step 105, according to the wind speed everywhere in ventilating system and electric machine rotor structural parameters, calculates the coefficient of heat transfer on motor surface by hydrodynamics method formula and law of conservation of mass;
The coefficient of heat transfer on motor surface specifically comprises:
A) yoke portion coefficient of heat transfer α in motor Radial ventilation duct e,
in formula, v efor the mean wind speed in yoke portion in Radial ventilation duct;
B) teeth portion coefficient of heat transfer α in motor Radial ventilation duct c,
in formula, v cfor the mean wind speed of teeth portion in Radial ventilation duct;
C) the coefficient of heat transfer α of insulation sides in air duct b,
in formula, v bfor the wind speed of insulation sides in air duct;
D) the convection coefficient α of air gap surface δ,
in formula, ω δfor mean wind speed, ω δ=v/2,
V is rotor peripheral speed, wherein n nfor motor speed, D 2for rotor outer circle diameter;
E) the coefficient of heat transfer α of stator core cylindrical and rotor core inner circle w,
in formula, D hfor equivalent diameter, λ is Measurement of Gas Thermal Conductivity, N ufor nusselt number,
Nusselt number N ucalculate with following formula:
N u = 0.152 R e 0.9 P r 0.833 [ 2.25 ln ( 0.144 R e 0.9 ) + 13.2 P r - 5.8 ] ,
In formula, P rfor Prandtl number, R efor Reynolds number.
Above wind speed calculating by fluid field, the coefficient of heat transfer at other positions also can use the same method and obtain.
Step 106, every loss value when motor step 103 drawn runs is loaded in three dimensional fluid field model as thermal source, coefficient of heat transfer step 105 calculated is as boundary condition, and other boundary condition is set: environment temperature, inlet boundary, outlet border and periodic boundary, adopt finite volume method to carry out Coupled Numerical Simulation to the zoning in electric machine rotor temperature field under different operating mode, finally obtain the thermo parameters method of motor.
Described motor is large-scale squirrel cage asynchronous motor, and described finite element software is ANSYS software.
The invention provides a kind of easily across the coupling calculation between physical field, carry out motor multiple physical field coupling analysis across different software, and calculate the thermo parameters method of motor.Compared with traditional electric machine temperature rise computing method, the inventive method applies the function of existing CAE software to greatest extent, and improve electric machine temperature rise computational accuracy, save the time, for the optimal design of large-size asynchronous motor, the enforcement of control strategy, safe operation provide theoretical foundation.

Claims (10)

1., based on an electric motor temperature field data processing method for multi-scenarios method, it is characterized in that, comprise the following steps:
1) use finite element software to set up physical model of electrical machine, comprise and set up two-dimensional electromagnetic field model according to Theory of Electromagnetic Field and set up three dimensional fluid field model according to Hydrodynamics Theory;
2) for two-dimensional electromagnetic field model, Finite Element Method is utilized to set up the field circuit method finite element equation solved in territory, calculate the magnetic distribution of motor, comprise magnetic flux distribution and distribution of current, and calculate stator coil copper loss, rotor bar copper loss and stator and rotor cores loss;
3) for three dimensional fluid field model, choose zoning and simplify three dimensional fluid field model, model after simplification comprises air gap flow field and Temperature calculating physical model and stator-rotor iron core, winding solid temperature field computational physics model between rotor, three dimensional fluid field model after simplification is imported in finite element software and carries out stress and strain model, adopt finite volume method to carry out iterative computation to zoning, obtain the fluid field distribution of motor;
4) according to the wind speed everywhere in ventilating system and electric machine rotor structural parameters, the coefficient of heat transfer on motor surface is calculated by hydrodynamics method formula and law of conservation of mass;
5) using step 2) the every loss value of the motor that draws of electromagnetic field analysis when running be loaded in three dimensional fluid field model as thermal source, according to step 4) coefficient of heat transfer that calculates draws all kinds of boundary condition, adopt finite volume method to carry out Coupled Numerical Simulation to the zoning in electric machine rotor temperature field under different operating mode, finally obtain the thermo parameters method of motor.
2. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, is characterized in that, step 1) described in two-dimensional electromagnetic field model be motor shaft to two-dimensional section.
3. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, is characterized in that, step 2) described in calculating stator coil copper loss P cucomputing formula be
P cu = &Sigma; I x 2 R x ,
In formula, I xfor the electric current in stator winding x, R xduring for being converted to benchmark job temperature, the resistance of winding x, carries out aftertreatment and obtains coil of stator of motor copper loss in finite element software.
4. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, is characterized in that, step 2) described in calculating rotor bar copper loss computing formula be
P cu &prime; = &Sigma; I y 2 R y ,
In formula, I yfor the electric current in rotor windings y, R yduring for being converted to benchmark job temperature, the resistance of winding y, carries out aftertreatment and obtains the loss of rotor conductive copper in finite element software.
5. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, is characterized in that, step 2) described in stator and rotor cores loss comprise the basic loss of stator core and stator core added losses P fep, the basic loss of described stator core comprises yoke portion basic loss P unshakable in one's determination fejbasic loss P unshakable in one's determination with teeth portion fet.
6. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 5, is characterized in that, described calculating yoke portion basic loss P unshakable in one's determination fejexpression formula be P fej=K αp fejg j, in formula, G jfor yoke portion clean iron unshakable in one's determination; K αfor experience factor, for alternating current generator, as capacity P nduring < 100kVA, K α=1.5, as capacity P nduring>=100kVA, K α=1.3, for direct current generator, K α=3.6; Specific loss p fejcomputing formula be
p Fej = p 10 / 50 B j 2 ( f f 0 ) 1.3 ,
In formula, B jfor peakflux density value in yoke, P 10/50for working as B=1T, during F=50Hz, the siliconized plate loss in unit weight;
Described calculating teeth portion basic loss P unshakable in one's determination fetexpression formula be P fet=K βp fetg t, in formula, G tfor teeth portion clean iron unshakable in one's determination; K βfor experience factor, for asynchronous machine, as capacity P nduring < 100kVA, K β=2.0, as capacity P nduring>=100kVA, K β=1.7, for direct current generator, K β=4.0; Specific loss p fetcomputing formula be
p Fet = p 10 / 50 B t 2 ( f f 0 ) 1.3 ,
In formula, B tfor the mean value of teeth portion magnetic flux density.
7. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 5, is characterized in that, described calculating stator core added losses P fepexperimental formula be
P Fep=K 0(B 0t) 2(zn) 1.5A p，
In formula, K 0for the coefficient relevant with magnetic pole material therefor, B 0for slot ripples magnetic flux density maximal value, t is armature tooth distance, and z is the number of teeth, and n is magnetic pole rotating speed, A pfor magnetic pole area.
8. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, it is characterized in that, step 3) described in zoning concrete grammar of choosing be in electric machine rotor circumference, choose 1/12nd of rotor as zoning.
9. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, is characterized in that, step 4) described in the coefficient of heat transfer on motor surface specifically comprise:
A) yoke portion coefficient of heat transfer α in motor Radial ventilation duct e,
in formula, v efor the mean wind speed in yoke portion in Radial ventilation duct;
B) teeth portion coefficient of heat transfer α in motor Radial ventilation duct c,
in formula, v cfor the mean wind speed of teeth portion in Radial ventilation duct;
C) the coefficient of heat transfer α of insulation sides in air duct b,
in formula, v bfor the wind speed of insulation sides in air duct;
D) the convection coefficient α of air gap surface δ,
α δ=28 (1+ ω δ 0.5), in formula, ω δfor mean wind speed, ω δ=v/2,
V is rotor peripheral speed, wherein n nfor motor speed, D 2for rotor outer circle diameter;
E) the coefficient of heat transfer α of stator core cylindrical and rotor core inner circle w,
in formula, D hfor equivalent diameter, λ is Measurement of Gas Thermal Conductivity, N ufor nusselt number,
Nusselt number N ucalculate with following formula:
N u = 0.152 R e 0.9 P r 0.833 [ 2.25 ln ( 0.144 R e 0.9 ) + 13.2 P r - 5.8 ] ,
In formula, P rfor Prandtl number, R efor Reynolds number.
10. a kind of electric motor temperature field data processing method based on multi-scenarios method according to claim 1, it is characterized in that, described motor is large-scale squirrel cage asynchronous motor.
CN201310422756.0A 2013-09-16 2013-09-16 Motor temperature field data processing method based on multi-field coupling Pending CN104462624A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310422756.0A CN104462624A (en) 2013-09-16 2013-09-16 Motor temperature field data processing method based on multi-field coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310422756.0A CN104462624A (en) 2013-09-16 2013-09-16 Motor temperature field data processing method based on multi-field coupling

Publications (1)

Publication Number Publication Date
CN104462624A true CN104462624A (en) 2015-03-25

Family

ID=52908659

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310422756.0A Pending CN104462624A (en) 2013-09-16 2013-09-16 Motor temperature field data processing method based on multi-field coupling

Country Status (1)

Country Link
CN (1) CN104462624A (en)

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105116336A (en) * 2015-10-09 2015-12-02 中国农业大学 Testing device and testing method for online loading of three-phase asynchronous motor online loading
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
CN105354354A (en) * 2015-09-28 2016-02-24 沈阳航空航天大学 Method for calculating temperature field of main driving motor of electric aircraft
CN105404751A (en) * 2015-12-07 2016-03-16 中国电子科技集团公司第三十二研究所 Method for realizing network unification based on heat-force-electromagnetic field
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
CN106383971A (en) * 2016-10-28 2017-02-08 沈阳工业大学 Improved analytical model of motor stator core vibration induced by magnetostriction
CN106599429A (en) * 2016-12-06 2017-04-26 中南大学 Fault injection method and fault injector of squirrel cage asynchronous motor rotor conducting bar
CN106777459A (en) * 2016-11-10 2017-05-31 北京交通大学 The computational methods in generator shaft radial rotor temperature field
CN107025366A (en) * 2017-05-23 2017-08-08 中南大学 Composite autoclave molding temperature field interactive mode approach of coupled numerical simulation
CN107291970A (en) * 2016-12-19 2017-10-24 沈阳工业大学 A kind of computational methods in air to air cooling middle size motor temperature field
CN107391884A (en) * 2017-08-28 2017-11-24 天津大学 Double remaining permagnetic synchronous motor temperature rise calculation methods based on equivalent thermal network model
CN107633143A (en) * 2017-09-25 2018-01-26 湖南大学 A kind of large-scale directly driven wind-powered machine stator radiator structure Optimization Design
CN108052724A (en) * 2017-12-08 2018-05-18 长沙学院 A kind of unit cost annual energy production optimised method of permanent-magnetic wind driven generator
CN109063337A (en) * 2018-08-03 2018-12-21 内蒙古工业大学 A kind of more coupling magnetic field with electric circuit simulation methods of permanent-magnetic wind driven generator
CN109067073A (en) * 2018-08-31 2018-12-21 温州大学 The method of totally-enclosed forced air-cooled structure permanent magnet traction motor radiator structure optimization
CN109190172A (en) * 2018-08-02 2019-01-11 北京交通大学 The calculation method of high speed rotor of motor air gap temperature field based on air gap rotating flow
CN109347405A (en) * 2018-10-08 2019-02-15 重庆长安汽车股份有限公司 A kind of evaluation method and estimating system of motor rotor temperature
CN109376475A (en) * 2018-11-27 2019-02-22 哈尔滨工业大学 A kind of multiturn densely packed coil Calculation Method of Temperature Field of Thermal self-ignition in conjunction with finite difference calculus
CN109900477A (en) * 2019-02-25 2019-06-18 河海大学 Double winding bearing-free switch reluctance motor models for temperature field based on heat source subdivision
CN105354354B (en) * 2015-09-28 2019-07-16 沈阳航空航天大学 A kind of electric airplane main drive motor Calculation Method of Temperature Field
CN110031506A (en) * 2019-04-25 2019-07-19 北京交通大学 The calculation method of magneto air gap thermal coefficient
CN110224376A (en) * 2019-07-11 2019-09-10 哈尔滨理工大学 The determination method of more duty three phase induction motor protection system thresholds
CN110412462A (en) * 2019-07-06 2019-11-05 陕西航空电气有限责任公司 A kind of aviation permasyn morot transient temperature prediction technique
CN110427664A (en) * 2019-07-18 2019-11-08 大连理工大学 A kind of permanent magnetic coupling temperature field analysis method based on magnetic-thermo coupling
CN110555289A (en) * 2019-09-27 2019-12-10 哈尔滨理工大学 Cloud computing-based motor stator winding multi-physical field coupling computing platform and computing method
CN110990973A (en) * 2019-12-08 2020-04-10 国网江苏省电力有限公司检修分公司 Electromagnetic field simulation calculation method and system for large synchronous phase modulator
CN111159860A (en) * 2019-12-16 2020-05-15 凡登（江苏）新型材料有限公司 Method for calculating temperature distribution in annealing process of copper strip of annealing device
CN111159953A (en) * 2019-12-31 2020-05-15 中国科学院电工研究所 Physical field analysis method, device and system of low-temperature motor
CN111222251A (en) * 2019-12-08 2020-06-02 国网江苏省电力有限公司检修分公司 Large synchronous phase modulator stator temperature field calculation method and system
CN112446150A (en) * 2020-11-25 2021-03-05 国网江苏省电力有限公司检修分公司 Three-dimensional visualization method and system for temperature field of double-water internal cooling synchronous phase modifier
CN112560367A (en) * 2020-11-27 2021-03-26 上海电气电站设备有限公司 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator
CN112597628A (en) * 2020-11-30 2021-04-02 广西电网有限责任公司电力科学研究院 Method and system for estimating temperature rise of transformer
CN112800685A (en) * 2021-03-26 2021-05-14 湖南大学 High-precision motor multi-physical field coupling simulation calculation method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080029506A1 (en) * 2006-07-19 2008-02-07 Encap Technologies Inc. Electromagnetic device with integrated fluid flow path
CN103279627A (en) * 2013-06-17 2013-09-04 清华大学 Heat-machinery-abrasion coupling analysis numerical simulation method based on finite element
CN103886206A (en) * 2014-03-25 2014-06-25 河海大学 Method and system for analyzing influence of joints on shaft voltage of power generator

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080029506A1 (en) * 2006-07-19 2008-02-07 Encap Technologies Inc. Electromagnetic device with integrated fluid flow path
CN103279627A (en) * 2013-06-17 2013-09-04 清华大学 Heat-machinery-abrasion coupling analysis numerical simulation method based on finite element
CN103886206A (en) * 2014-03-25 2014-06-25 河海大学 Method and system for analyzing influence of joints on shaft voltage of power generator

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
丁舜年: "《大型电机的发热与冷却》", 31 July 1992 *
周封: "大型同步电机内流场及温度场耦合的研究", 《万方学位论文全文数据库》 *
孔晓光等: "高速永磁电机的损耗计算与温度场分析", 《电工技术学报》 *
李伟力等: "感应电动机定转子全域温度场数值计算及相关因素敏感性分析", 《中国电机工程学报》 *
温嘉斌等: "中型高压异步电动机三维温度场耦合计算与分析", 《电机与控制学报》 *
魏永田等: "《电机内热交换》", 31 August 1998 *

Cited By (48)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105354354A (en) * 2015-09-28 2016-02-24 沈阳航空航天大学 Method for calculating temperature field of main driving motor of electric aircraft
CN105354354B (en) * 2015-09-28 2019-07-16 沈阳航空航天大学 A kind of electric airplane main drive motor Calculation Method of Temperature Field
CN105116336A (en) * 2015-10-09 2015-12-02 中国农业大学 Testing device and testing method for online loading of three-phase asynchronous motor online loading
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
CN105404751B (en) * 2015-12-07 2018-10-23 中国电子科技集团公司第三十二研究所 Method for realizing network unification based on heat-force-electromagnetic field
CN105404751A (en) * 2015-12-07 2016-03-16 中国电子科技集团公司第三十二研究所 Method for realizing network unification based on heat-force-electromagnetic field
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
CN106126861B (en) * 2016-07-18 2019-04-30 北京交通大学 The calculation method of the coefficient of heat transfer in air-cooled steam turbine generator teeth portion and back of the yoke portion
CN106383971A (en) * 2016-10-28 2017-02-08 沈阳工业大学 Improved analytical model of motor stator core vibration induced by magnetostriction
CN106383971B (en) * 2016-10-28 2019-12-27 沈阳工业大学 Improved motor stator core vibration analysis model caused by magnetostriction
CN106777459A (en) * 2016-11-10 2017-05-31 北京交通大学 The computational methods in generator shaft radial rotor temperature field
CN106599429A (en) * 2016-12-06 2017-04-26 中南大学 Fault injection method and fault injector of squirrel cage asynchronous motor rotor conducting bar
CN107291970A (en) * 2016-12-19 2017-10-24 沈阳工业大学 A kind of computational methods in air to air cooling middle size motor temperature field
CN107025366A (en) * 2017-05-23 2017-08-08 中南大学 Composite autoclave molding temperature field interactive mode approach of coupled numerical simulation
CN107025366B (en) * 2017-05-23 2020-05-15 中南大学 Interactive coupling numerical simulation method for composite material autoclave molding temperature field
CN107391884A (en) * 2017-08-28 2017-11-24 天津大学 Double remaining permagnetic synchronous motor temperature rise calculation methods based on equivalent thermal network model
CN107391884B (en) * 2017-08-28 2020-06-26 天津大学 Temperature rise calculation method of dual-redundancy permanent magnet synchronous motor based on equivalent thermal network model
CN107633143B (en) * 2017-09-25 2019-11-19 湖南大学 A kind of directly driven wind-powered machine stator radiator structure optimum design method of large size
CN107633143A (en) * 2017-09-25 2018-01-26 湖南大学 A kind of large-scale directly driven wind-powered machine stator radiator structure Optimization Design
CN108052724A (en) * 2017-12-08 2018-05-18 长沙学院 A kind of unit cost annual energy production optimised method of permanent-magnetic wind driven generator
CN109190172A (en) * 2018-08-02 2019-01-11 北京交通大学 The calculation method of high speed rotor of motor air gap temperature field based on air gap rotating flow
CN109063337B (en) * 2018-08-03 2022-03-29 内蒙古工业大学 Multi-field-road coupling simulation calculation method for permanent magnet wind driven generator
CN109063337A (en) * 2018-08-03 2018-12-21 内蒙古工业大学 A kind of more coupling magnetic field with electric circuit simulation methods of permanent-magnetic wind driven generator
CN109067073B (en) * 2018-08-31 2020-07-28 温州大学 Method for optimizing heat dissipation structure of permanent magnet traction motor with totally-enclosed forced air cooling structure
CN109067073A (en) * 2018-08-31 2018-12-21 温州大学 The method of totally-enclosed forced air-cooled structure permanent magnet traction motor radiator structure optimization
CN109347405A (en) * 2018-10-08 2019-02-15 重庆长安汽车股份有限公司 A kind of evaluation method and estimating system of motor rotor temperature
CN109376475A (en) * 2018-11-27 2019-02-22 哈尔滨工业大学 A kind of multiturn densely packed coil Calculation Method of Temperature Field of Thermal self-ignition in conjunction with finite difference calculus
CN109900477A (en) * 2019-02-25 2019-06-18 河海大学 Double winding bearing-free switch reluctance motor models for temperature field based on heat source subdivision
CN110031506A (en) * 2019-04-25 2019-07-19 北京交通大学 The calculation method of magneto air gap thermal coefficient
CN110412462A (en) * 2019-07-06 2019-11-05 陕西航空电气有限责任公司 A kind of aviation permasyn morot transient temperature prediction technique
CN110412462B (en) * 2019-07-06 2021-06-08 陕西航空电气有限责任公司 Transient temperature prediction method for permanent magnet synchronous motor for aviation
CN110224376A (en) * 2019-07-11 2019-09-10 哈尔滨理工大学 The determination method of more duty three phase induction motor protection system thresholds
CN110224376B (en) * 2019-07-11 2021-05-28 哈尔滨理工大学 Method for determining threshold value of multi-working-system three-phase induction motor protection system
CN110427664A (en) * 2019-07-18 2019-11-08 大连理工大学 A kind of permanent magnetic coupling temperature field analysis method based on magnetic-thermo coupling
CN110555289A (en) * 2019-09-27 2019-12-10 哈尔滨理工大学 Cloud computing-based motor stator winding multi-physical field coupling computing platform and computing method
CN110555289B (en) * 2019-09-27 2023-12-15 哈尔滨理工大学 Motor stator winding multi-physical field coupling computing platform and computing method based on cloud computing
CN111222251A (en) * 2019-12-08 2020-06-02 国网江苏省电力有限公司检修分公司 Large synchronous phase modulator stator temperature field calculation method and system
CN111222251B (en) * 2019-12-08 2023-09-26 国网江苏省电力有限公司检修分公司 Large synchronous camera stator temperature field calculation method and system
CN110990973A (en) * 2019-12-08 2020-04-10 国网江苏省电力有限公司检修分公司 Electromagnetic field simulation calculation method and system for large synchronous phase modulator
CN111159860A (en) * 2019-12-16 2020-05-15 凡登（江苏）新型材料有限公司 Method for calculating temperature distribution in annealing process of copper strip of annealing device
CN111159953A (en) * 2019-12-31 2020-05-15 中国科学院电工研究所 Physical field analysis method, device and system of low-temperature motor
CN111159953B (en) * 2019-12-31 2024-04-05 中国科学院电工研究所 Physical field analysis method, device and system of low-temperature motor
CN112446150A (en) * 2020-11-25 2021-03-05 国网江苏省电力有限公司检修分公司 Three-dimensional visualization method and system for temperature field of double-water internal cooling synchronous phase modifier
CN112446150B (en) * 2020-11-25 2024-06-04 国网江苏省电力有限公司检修分公司 Three-dimensional visualization method and system for temperature field of double-water internal cooling synchronous camera
CN112560367B (en) * 2020-11-27 2022-05-27 上海电气电站设备有限公司 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator
CN112560367A (en) * 2020-11-27 2021-03-26 上海电气电站设备有限公司 Method for calculating hot spot temperature of axial double-hole copper bar rotor coil of generator
CN112597628A (en) * 2020-11-30 2021-04-02 广西电网有限责任公司电力科学研究院 Method and system for estimating temperature rise of transformer
CN112800685A (en) * 2021-03-26 2021-05-14 湖南大学 High-precision motor multi-physical field coupling simulation calculation method

Similar Documents

Publication Publication Date Title
CN104462624A (en) Motor temperature field data processing method based on multi-field coupling
Fan et al. Water cold plates for efficient cooling: Verified on a permanent-magnet machine with concentrated winding
Fan et al. Ventilation and thermal improvement of radial forced air-cooled FSCW permanent magnet synchronous wind generators
Kang et al. Self-circulation cooling structure design of permanent magnet machines for electric vehicle
Popescu et al. Study of the thermal aspects in brushless permanent magnet machines performance
Fan et al. A dynamic multilayer winding thermal model for electrical machines with concentrated windings
Liang et al. Thermal design and optimization of a water-cooling permanent magnet synchronous in-wheel motor
Zhu et al. Cooling system design of a high-speed PMSM based on a coupled fluidic–thermal model
Xu et al. Research on thermal capacity of a high-torque-density direct drive permanent magnet synchronous machine based on a temperature cycling module
Tong et al. Research on the airflow and thermal performance in a large forced air-cooled permanent magnet synchronous machine
Yu et al. Cooling analysis of high-speed stator-permanent magnet flux-switching machines for fuel-cell electric vehicle compressor
Cheng et al. Research on induction motor for mini electric vehicles
Li et al. Numerical calculation of multicoupled fields in large salient synchronous generator
CN111324974B (en) Optimization method and device of air-cooled generator based on stator tooth internal cooling ventilation structure
Nategh et al. Thermal analysis of permanent-magnet synchronous reluctance machines
CN105243249A (en) Three-phase induction motor transient temperature rise calculation model and measuring calculation method
Xia et al. Coupled airflow-thermal network model for large-scale offshore wind generators with modular fractional-slot concentrated winding
Feiyang et al. Influence of copper shield structure on 3-D electromagnetic field, fluid and temperature fields in end region of large turbogenerator
Ni et al. Thermal analysis of a hybride excitation flux-switching motor with water-cooling system
Guo et al. Thermal Design and Simulation of winding cooling for permanent magnet synchronous motor of electric vehicle
Li et al. Design and analysis of different cooling schemes of a flux‐modulated permanent magnet in‐wheel motor for electric vehicle applications
Liang et al. Combined lumped-parameter and simplified 2-D analytical thermal model of totally enclosed water cooled PM machine
Zhao et al. Design and optimization of water cooling system for high power density permanent magnet motor based on multi-physics field simulation
Sim et al. Temperature prediction of oil-cooled IPMSM for in-wheel direct-drive through lumped parameter thermal model
Tong et al. Cooling system design and thermal analysis of multibrid permanent magnet wind generator

Legal Events

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

Application publication date: 20150325