CN101769797A - Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor - Google Patents
Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor Download PDFInfo
- Publication number
- CN101769797A CN101769797A CN200910076036A CN200910076036A CN101769797A CN 101769797 A CN101769797 A CN 101769797A CN 200910076036 A CN200910076036 A CN 200910076036A CN 200910076036 A CN200910076036 A CN 200910076036A CN 101769797 A CN101769797 A CN 101769797A
- Authority
- CN
- China
- Prior art keywords
- temperature
- permanent magnet
- thermal
- model
- heat
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
The invention relates to a temperature rise analytical method for predicting temperature of a permanent magnet in a permanent magnet synchronous motor (PMSM), belonging to the application electrical engineering design field; the method is characterized in that: distributed heat source of a motor is analyzed by a filed-circuit compact coupling method, comprising eddy current loss in the permanent magnet, iron loss in an iron core and copper loss in armature; on the consideration of precision requirements, the coupling analysis of a magnetic field and a temperature filed can be realized by single-way coupling mode. A thermal model of the permanent magnet synchronous motor is built based on a mixing method of a novel equivalent heat network and a finite element, heat parameters are rationally selected by adopting a combining mode of experimental measurement and empirical formula, the heat transferring coefficient and cooling condition of the motor are described completely, a stator and a rotor can be systematically combined by adopting air gap joints in the heat network, the stator and rotor unified temperature rise model is formed, the difficulty of measuring air gap temperature is avoided, material parameters are adopted at the practical working temperature, so as to lead the analysis to be rational; the accurate and optical method for predicting the temperature of the permanent magnet is realized by special correction processing in experimental links; in addition, the design method is used to give suggestions for model selection of the permanent magnet material in the motor.
Description
Technical field
The invention belongs to the application electrical technology field, relate in particular to the temperature prediction analysis of permasyn morot (PMSM).
Background technology
Permasyn morot (PMSM) has efficient energy-saving, characteristics such as simple in structure, be that permagnetic synchronous motor need not dc excitation power and contact devices such as brush that attaches and slip ring specifically, can reduce machinery and electrical loss, the power factor height, starting moment is big etc.In order to design high efficiency motor, the temperature prediction of permanent magnet is absolutely necessary in the permasyn morot, particularly because the rotor radiating condition is poor, and the eddy current loss of permanent magnet might cause higher temperature rise, and then causes the permanent magnet loss of excitation, reduces motor performance.The main considerations of its temperature prediction has:
(1) permanent magnetic material is to the influence of motor performance.
(2) different carrier frequencies and modulation ratio are to the influence of eddy current loss in the permanent magnet.
(3) Design of ventilation of optimal design motor guarantees that good cooling condition has extremely important effect to reducing temperature rise.
In the design of the permasyn morot of low capacity, its temperature rise is not caused enough attention in early days.In fact,, in permanent magnet, have eddy current, and along with the raising of power of motor, this loss probably can cause big temperature rise owing to be subjected to the effect of magnetic field space harmonic wave and time harmonic.Neodymium iron boron is compared ferrite as using maximum permanent magnet materials, though coercive force and remanent magnetism height, conductivity height and poor heat resistance, because the rotor heat-sinking capability is poor, eddy current can make the permanent magnet heating, can produce the danger of degaussing.
The coupling analysis that temperature rise research of permagnetic synchronous motor is related to magnetic field and temperature field, general way utilizes actual measurement air gap temperature curve to carry out heat analysis earlier, rotor is handled as finite element model independently, the loss source obtains by simplified model, error on the thermal source will be the inaccurate key factor of thermal model, and will accurate in advance measurement air gap temperature curve be the difficult point problem.Aspect temperature rise calculating, Traditional calculating methods mostly adopts experimental formula to obtain thermal parameter, thermal model has been taked more approximate, in accurate general calculation method not very analytically.
Pass through on the basis of emi analysis eddy current loss among the present invention, increase other aspect loss analysis to obtain the complete thermal source distribution of this motor, adopt the unidirectional couplings mode to realize the coupling analysis in magnetic field and temperature field, set up thermal model based on the mixed method that a kind of novel equivalent thermal network and finite element combine, measure by experiment and mode that experimental formula combines is rationally chosen thermal parameter, the coefficient of heat transfer of complete description motor and cooling condition.By the air gap node in the heat supply network network rotor is organically combined, form the unified temperature rise model of rotor, can avoid difficulties such as air gap temperature measurement, the material parameter when adopting actual work temperature, make to analyze to have more applicability, realized predicting more exactly the requirement of permanent magnet temperature.
Summary of the invention
The object of the present invention is to provide a kind of mixing temperature rise analytical approach of predicting permanent magnet temperature in the permasyn morot exactly, its step is as follows:
Step (1) adopts a tight road coupling process to find the solution magnetic field model, and thermal source distributes to obtain accurately:
The field is described with the Maxwell equation, and the road is described with Kirchhoff's equation, and the equation of motion is described with torque balance, at last magnetic vector position A and current i is found the solution as variable simultaneously.Above-mentioned equation forms the discretization equation of non-linear matrix form with Galerkin method finite element variation form discretize.Afterwards, the time is discrete by implicit Euler method, obtains recursion equation.Non-linear matrix equation mountain Newton-Raphson method is found the solution, and along with the motor rotation, the magnetic vector position and the electric current in each step can obtain.
Motor runs on rated speed, at first carries out zero load and analyzes, obtain back-emf and with the initial phase of rotor reference position.
The rated voltage that will have certain phasing degree is imported the rated power that obtains motor in the road coupling model, obtains the eddy current in its rotating magnetic field, electric current and the permanent magnet.
Eddy current in the permanent magnet is described with finite element method, calculates its eddy current loss with Ohm law then, and ignores magnetic hysteresis loss.Iron loss unshakable in one's determination is tried to achieve by the experimental formula that material database provides, and comprises that the magnetic flux density data are carried out FFT decomposes joule loss that obtains and the magnetic hysteresis loss that obtains according to area of hysteresis loop.Armature ohmic loss is calculated by Ohm law according to measured resistance value.
Described magnetic field model analysis comprises the voltage source input at least.
Step (2), heat of mixing analytical approach is taked the method that combines with suitable heat supply network network, utilizes finite element method can estimate high-precision motor internal Temperature Distribution.
Step (2.1) is being ignored under the situation of radiation, and during steady-state operation, the temperature partial derivative is zero, and boundary value problem is expressed as:
Wherein, S
1, S
2Be respectively adiabatic boundary surface and convection heat transfer boundary surface, k is the coefficient of heat conductivity of boundary surface vertical direction, and α is a coefficient of heat transfer, T
fIt is coolant temperature.
Then, list and said temperature field partial differential equation equivalent conditions variational problem, carry out element subdivision, set up the linear algebraic equation group and can find the solution the temperature value that obtains each point.
Step (2.2), the equivalent thermal network method combines with finite element model, has considered thermal resistance and the casing and the extraneous hot-fluid exchange of motor internal air and casing; Having simplified some influences less thermal parameter to the result, comprising: ignore the heat interchange of bearing and casing, the thermal resistance of bearing and rotor core thinks that it directly contacts each other; But the thermal boundary of stator core and casing but is one of the key element that will accurately analyze, also has heat interchange between air-gap and the stator core with to streamed realization.
Eddy current loss thermal source major part will be passed to stator core and flow to casing again by carrying out forced convertion with air gap in the permanent magnet, and it is through passing to casing with the convection current of inner air that fraction is wherein arranged; Iron loss in the coil in copper loss and the stator core is carried out heat interchange with direct and casing.
At first, the method that adopts experiment measuring and experimental formula to combine is chosen suitable thermal parameter.When rotor rotated, the coefficient of heat transfer of rotor outer surface (generally was taken as 10W/ (m greater than the Natural Heat Convection coefficient
2K)), can calculate by following formula:
h
1＝6.6/10
5×V
0.67/lg
0.33W/(cm
2·K)
Wherein, V is a rotor speed, and cm/s, lg are gas lengths, and cm is converted into system international again.When the Reynolds number greater than 200 the time, this equation is set up; In this model, through calculating, the Reynolds number is 1984, and the air in the air gap is a turbulent flow, shows that its coefficient of heat conductivity is much larger than still air.
Then, calculate thermal contact resistance R
1Heat interchange between expression permanent magnet and the rotor core, R
2Heat interchange between expression coil and the stator core, it is relevant with coefficient of heat conductivity, thickness and the contact area of B-class insulation material.Between inner air and casing, use equivalent thermal resistance R
3Weigh, its computing formula is:
R
3＝1/(A
l×h
4)，
Wherein, A
lBe the casing inner surface area, h
4It is the thermal boundary condition between inner air and rotor core or the permanent magnet.
Heat interchange between computer shell and the environment is determined by experiment again, and casing is forced cooling through coaxial mounted fan, in order to increase the precision of calculating, the cooling of the conduction between casing and the ground is also considered in the model go.
Step (3), magnetic thermal coupling model, the order of magnitude of thermal time constant is big more a lot of than electromagnetic time constant in the permasyn morot, ignore the influence of temperature rise to electromagnetic parameter, consider the accuracy requirement of practical problems, adopt folk prescription to realize that to coupling scheme the coupling analysis in magnetic field and temperature field can significantly shorten computing time, the counting yield height.Complete the importing in the temperature rise model that matches of thermal source distribution results of previous magnetic model analysis, can carry out follow-up temperature field analysis.Its step is as follows:
Step (3.1) to computing machine input material magnetic attribute and the hot attribute of material, comprising: permanent magnet conductivity, core material BH curve, rotating shaft magnetic permeability, the coefficient of heat conductivity of various materials, specific heat capacity and density.
Step (3.2) is set up geometric model, loads boundary condition.In the magnetic model, utilize the periodic boundary condition of motor, choose an electric cycle modeling, upwards be 1/2 model therefore in week.Considering the axial stratification structure of permanent magnet, set up the three-dimensional element motor model on the two dimensional model basis, utilize axial symmetry, only need set up the motor section model of half block magnet thickness, is 1/30 model in the axial direction therefore.In thermal model, import lossy data for convenience, take and the model consistent element subdivision mode of magnetic; Side surface at model is provided with periodic boundary condition equally, to show that temperature is identical on this surface.In the magnetic model, on the branch aspect of permanent magnets, the vertical component of A is zero, and eddy current is parallel to this plane; On the plane of symmetry, the parallel component of A is zero, and φ is zero, and magnetic field is parallel to this plane, and eddy current is perpendicular to this plane.In thermal model, be provided with between permanent magnet, stator core and the air gap that thermal boundary condition shows and air-flow between convection heat transfer, it is a node in the heat supply network network.Set initial temperature for then whole model and surrounding environment, and radiation heat transfer is included in the convection heat transfer to improve computational accuracy.
Step (3.3), coefficient of heat transfer determine the computation process of thermal boundary condition, h
1Thermal boundary between expression permanent magnet and the air gap; h
2The thermal boundary of expression stator core or coil and air gap is with h
1Corresponding, the turbulence characteristic of air gap air has been described; h
3Thermal boundary between expression stator core and the casing, wherein k is the coefficient of heat conductivity (in the time of 20～120 ℃, span is 0.0259～0.0334, and the working temperature of casing is approximately 46 ℃, therefore this value is chosen as 0.027) of air; h
4Thermal boundary between expression permanent magnet or rotor core and the inner air, this value role in the heat supply network network is limited, takes estimated value among the present invention, and it should be to loose between natural convection to flow coefficient and h
1Between number.Pass through h among the present invention
1And h
2Etc. parameter stator and rotor are organically united, do the as a whole temperature rise analysis of carrying out, this helps accurately describing the conduction of motor internal thermal source, needing to have avoided a series of problems such as measurement air gap temperature.
Step (3.4), the computation process of thermal contact resistance, R
1What describe is thermal resistance between permanent magnet and the rotor core, and A is a contact area; R
2What describe is thermal resistance between coil and the stator core, and wherein k is the coefficient of heat conductivity (span is 0.12～0.16) of B-class insulation material, and A is the surface area of coil, promptly with the contact area of stator slot.
Step (3.5), the determining of cooling condition.R
3What describe is equivalent thermal resistance between inner air and the casing; T describes is hot-fluid between casing and the external environment condition, and promptly the cooling condition in this motor model will adopt the lumped parameter of equivalence to simulate h
5The thermal convection border of expression casing and extraneous air, it comprises two parts heat interchange factor, at first is the heat loss through convection coefficient of the base-plates surface that caused by fan cooled, and rule of thumb formula calculates, and specifically is expressed as follows
α＝11+3v
Wherein, v is a cross-ventilation speed, can be obtained by experiment measuring, and α equals 44W/ (m
2K).Another part is the heat interchange between casing and the ground, and the coefficient of heat conductivity of iron is 38W/ (mK), and sectional area is 2cm
2, length is 5cm, so the conductive heat flow of iron block is 0.152W/K, equivalence is 14W/ (m for coefficient of heat transfer
2K), A
2Be the surface area of casing, so h
5Equal 58W/ (m
2K).This value can oppositely check according to the casing temperature through infrared heat image instrument measuring and obtain, because this parametric description is the integral body heat transfer of motor, only influences the surface temperature of casing, can not influence the relative value of internal temperature.
Method for designing among the present invention, carry out the temperature rise analysis by the mixed method that novel equivalent thermal network and finite element combine, rotor is organically combined, consider the heat interchange activity between the rotor, whole motor is become a complete temperature rise analytical model, can avoid difficulties such as air gap temperature measurement, have more applicability, its analysis result also will have more with reference to property.Simultaneously, the mode that combines by experimental formula and experiment measuring provides thinking for the thermal parameter of the type motor, and particularly the correction that is treated to model about the casing cooling condition provides new method.Use this analysis method can predict the permanent magnet temperature more exactly, can provide some of the recommendations, and provide safety guarantee for design of electrical motor for the type selecting of permanent magnetic material in the motor.
Description of drawings
Fig. 1 is an equivalent thermal network of using the present invention's design.
Fig. 2 is the analysis process figure that uses the magnetic thermal coupling model of the present invention's design.
Fig. 3 is the boundary condition of using in the thermal model of the present invention's design.
Fig. 4 uses that motor temperature distributes under a kind of PWM voltage of the present invention's design.
Fig. 5 uses stator core Temperature Distribution under a kind of PWM voltage of the present invention's design.
Table 1 is a kind of motor master data of using the present invention's design.
Table 2 is hot attributes of a kind of material of using the present invention's design.
Table 3 is a computation process of using a kind of thermal boundary condition of the present invention's design.
Table 4 is a computation process of using a kind of casing thermal parameter of the present invention's design.
Embodiment
Below in conjunction with accompanying drawing principle of the present invention and concrete embodiment are described once.
Table 1 is the parameter and the material properties of simulating electric machine used in the embodiment of the invention, is a bench teat electricity-testing machine.
Fig. 1 is the equivalent thermal network that combines with finite element model selected in the embodiment of the invention.Each module all is the pairing finite element part of each parts in the motor, interknits by various thermal elements between them, and be an organic whole.Main hot-fluid approach has: eddy current loss thermal source major part will be passed to stator core and flow to casing again by carrying out forced convertion with air gap in the permanent magnet, and it is through passing to casing with the convection current of inner air that fraction is wherein arranged; Iron loss in the coil in copper loss and the stator core is carried out heat interchange with direct and casing.
Fig. 2 is the analysis process figure of magnetic thermal coupling model in the embodiment of the invention.Consider the accuracy requirement of practical problems, adopt folk prescription to realize that to coupling scheme the coupling analysis in magnetic field and temperature field can significantly shorten computing time,, can carry out the temperature rise analysis complete the importing in the thermal model of thermal source distributed data of magnetic analysis.
Table 1
Table 2 is hot attributes of material used in the embodiment of the invention, and the coefficient of heat transfer of precise material coefficient of heat conductivity and boundary surface is accurately to analyze the electric machine temperature rise model based, and the material parameter when adopting actual work temperature makes analysis more accurate.
Table 2
Fig. 3 is the boundary condition in the thermal model selected in the embodiment of the invention, and table 3 is the computation process of used thermal boundary condition in the embodiment of the invention, comprises the numerical value of its concrete intermediate variable, and the Reynolds number satisfies the formula demand in this model.Considered the thermal convection effect that the motor internal air is played simultaneously, erected with the heat interchange between the casing for inner some structure and be connected bridge.
Table 3
Table 4 is the computation process of casing thermal parameter in the embodiment of the invention, and concrete effect is to be expressed as determining of cooling condition.From heat loss through convection and and the earth between heat loss through conduction two aspect factors described this motor and be in cooling condition under the nominal operation state.Because this cooling condition can measure by experiment, thereby and proofread and correct concrete numerical value in the temperature rise model, so the temperature rise that the method for being taked among the present invention can be used for the type motor is effectively accurately analyzed.
Table 4
As a sample calculation, carry out the Research on Temperature Rise under this test motor PWM inverter power supply among the present invention.Fig. 4 is motor temperature distribution under the PWM voltage in the embodiment of the invention.The lossy data of each unit is assigned to coupling model respectively as distributed thermal source, and thermal model is provided with thermal boundary condition description and extraneous heat interchange.As can be seen, the temperature difference between rotor and the stator is very big, has 40 ℃ approximately, if do not cooled off sufficiently, permanent magnet burns possibly.As the high spot reviews object, the maximum temperature of permanent magnet is 123.1 ℃; Because permanent magnet is an even short and small heat conductor, the highest minimum temperature difference of whole magnetic patch has only about 3 ℃.On rotor surface Temperature Distribution in rotational direction, at the every utmost point magnet in positive veer back a higher temperature peak is arranged all as can be known, reverse directions has a small peak in addition, and this phenomenon is corresponding with the characteristic distributions of eddy current loss density, and just Temperature Distribution is more level and smooth.
Fig. 5 is a stator core Temperature Distribution under the PWM voltage in the embodiment of the invention because and a little less than the heat interchange between the rotor, so the temperature of stator inboard is little more a lot of than permanent magnet, Temperature Distribution presents low inside and high outside gradient, has shown the pipeline of heat.
Claims (2)
1. a temperature rise analytical approach that is used for predicting permasyn morot permanent magnet temperature is characterized in that, contains following steps successively:
Step (1) adopts a tight road coupling process to find the solution magnetic field model, and thermal source distributes to obtain accurately:
The field is described with the Maxwell equation, and the road is described with Kirchhoff's equation, and the equation of motion is described with torque balance, at last magnetic vector position A and current i is found the solution as variable simultaneously; Above-mentioned equation forms the discretization equation of non-linear matrix form with Galerkin method finite element variation form discretize; Afterwards, the time is discrete by implicit Euler method, obtains recursion equation; Non-linear matrix equation is found the solution by the Newton-Raphson method, and along with the motor rotation, the magnetic vector position and the electric current in each step can obtain;
Motor runs on rated speed, at first carries out zero load and analyzes, obtain back-emf and with the initial phase of rotor reference position;
The rated voltage that will have certain phasing degree is imported the rated power that obtains motor in the road coupling model, obtains the eddy current in its rotating magnetic field, electric current and the permanent magnet;
Eddy current in the permanent magnet is described with finite element method, calculates its eddy current loss with Ohm law then, and ignores magnetic hysteresis loss; Iron loss unshakable in one's determination is tried to achieve by the experimental formula that material database provides, and comprises that the magnetic flux density data are carried out FFT decomposes joule loss that obtains and the magnetic hysteresis loss that obtains according to area of hysteresis loop; Armature ohmic loss is calculated by Ohm law according to measured resistance value;
Described magnetic field model analysis comprises the voltage source input at least;
Step (2), heat of mixing analytical approach is taked the method that combines with suitable heat supply network network, utilizes finite element method can estimate high-precision motor internal Temperature Distribution;
Step (2.1) is being ignored under the situation of radiation, and during steady-state operation, the temperature partial derivative is zero, and boundary value problem is expressed as:
Wherein, S
1, S
2Be respectively adiabatic boundary surface and convection heat transfer boundary surface, k is the coefficient of heat conductivity of boundary surface vertical direction, and α is a coefficient of heat transfer, T
fIt is coolant temperature;
Then, list and said temperature field partial differential equation equivalent conditions variational problem, carry out element subdivision, set up the linear algebraic equation group and can find the solution the temperature value that obtains each point;
Step (2.2), the equivalent thermal network method combines with finite element model, has considered thermal resistance and the casing and the extraneous hot-fluid exchange of motor internal air and casing; Having simplified some influences less thermal parameter to the result, comprising: ignore the heat interchange of bearing and casing, the thermal resistance of bearing and rotor core thinks that it directly contacts each other; But the thermal boundary of stator core and casing but is one of the key element that will accurately analyze, also has heat interchange between air-gap and the stator core with to streamed realization;
Eddy current loss thermal source major part will be passed to stator core and flow to casing again by carrying out forced convertion with air gap in the permanent magnet, and it is through passing to casing with the convection current of inner air that fraction is wherein arranged; Iron loss in the coil in copper loss and the stator core is carried out heat interchange with direct and casing;
At first, the method that adopts experiment measuring and experimental formula to combine is chosen suitable thermal parameter; When rotor rotated, the coefficient of heat transfer of rotor outer surface (generally was taken as 10W/ (m greater than the Natural Heat Convection coefficient
2K)), can calculate by following formula:
h
1＝6.6/10
5×V
0.67/lg
0.33W/(cm
2·K)
Wherein, V is a rotor speed, and cm/s, lg are gas lengths, and cm is converted into system international again.When the Reynolds number greater than 200 the time, this equation is set up; In this model, through calculating, the Reynolds number is 1984, and the air in the air gap is a turbulent flow, shows that its coefficient of heat conductivity is much larger than still air;
Then, calculate thermal contact resistance R
1Heat interchange between expression permanent magnet and the rotor core, R
2Heat interchange between expression coil and the stator core, it is relevant with coefficient of heat conductivity, thickness and the contact area of B-class insulation material; Between inner air and casing, use equivalent thermal resistance R
3Weigh, its computing formula is:
R
3＝1/(A
1×h
4)，
Wherein, A
1Be the casing inner surface area, h
4It is the thermal boundary condition between inner air and rotor core or the permanent magnet;
Heat interchange between computer shell and the environment is determined by experiment again, and casing is forced cooling through coaxial mounted fan, in order to increase the precision of calculating, the cooling of the conduction between casing and the ground is also considered in the model go;
Step (3), magnetic thermal coupling model, the order of magnitude of thermal time constant is big more a lot of than electromagnetic time constant in the permasyn morot, ignore the influence of temperature rise to electromagnetic parameter, consider the accuracy requirement of practical problems, adopt folk prescription to realize that to coupling scheme the coupling analysis in magnetic field and temperature field can significantly shorten computing time, the counting yield height; Complete the importing in the temperature rise model that matches of thermal source distribution results of previous magnetic model analysis, can carry out follow-up temperature field analysis; Its step is as follows:
Step (3.1) to computing machine input material magnetic attribute and the hot attribute of material, comprising: permanent magnet conductivity, core material BH curve, rotating shaft magnetic permeability, the coefficient of heat conductivity of various materials, specific heat capacity and density;
Step (3.2) is set up geometric model, loads boundary condition; In the magnetic model, utilize the periodic boundary condition of motor, choose an electric cycle modeling, upwards be 1/2 model therefore in week; Considering the axial stratification structure of permanent magnet, set up the three-dimensional element motor model on the two dimensional model basis, utilize axial symmetry, only need set up the motor section model of half block magnet thickness, is 1/30 model in the axial direction therefore; In thermal model, import lossy data for convenience, take and the model consistent element subdivision mode of magnetic; Side surface at model is provided with periodic boundary condition equally, to show that temperature is identical on this surface; In the magnetic model, on the branch aspect of permanent magnets, the vertical component of A is zero, and eddy current is parallel to this plane; On the plane of symmetry, the parallel component of A is zero, and φ is zero, and magnetic field is parallel to this plane, and eddy current is perpendicular to this plane; In thermal model, be provided with between permanent magnet, stator core and the air gap that thermal boundary condition shows and air-flow between convection heat transfer, it is a node in the heat supply network network; Set initial temperature for then whole model and surrounding environment, and radiation heat transfer is included in the convection heat transfer to improve computational accuracy;
Step (3.3), coefficient of heat transfer determine the computation process of thermal boundary condition, h
1Thermal boundary between expression permanent magnet and the air gap; h
2The thermal boundary of expression stator core or coil and air gap is with h
1Corresponding, the turbulence characteristic of air gap air has been described; h
3Thermal boundary between expression stator core and the casing, wherein k is the coefficient of heat conductivity (in the time of 20～120 ℃, span is 0.0259～0.0334, and the working temperature of casing is approximately 46 ℃, therefore this value is chosen as 0.027) of air; h
4Thermal boundary between expression permanent magnet or rotor core and the inner air, this value role in the heat supply network network is limited, takes estimated value among the present invention, and it should be to loose between natural convection to flow coefficient and h
1Between several the present invention in pass through h
1And h
2Etc. parameter stator and rotor are organically united, do the as a whole temperature rise analysis of carrying out, this helps accurately describing the conduction of motor internal thermal source, needing to have avoided a series of problems such as measurement air gap temperature;
Step (3.4), the computation process of thermal contact resistance, R
1What describe is thermal resistance between permanent magnet and the rotor core, and A is a contact area; R
2What describe is thermal resistance between coil and the stator core, and wherein k is the coefficient of heat conductivity (span is 0.12～0.16) of B-class insulation material, and A is the surface area of coil, promptly with the contact area of stator slot;
Step (3.5), the determining of cooling condition.R
3What describe is equivalent thermal resistance between inner air and the casing; T describes is hot-fluid between casing and the external environment condition, and promptly the cooling condition in this motor model will adopt the lumped parameter of equivalence to simulate h
5The thermal convection border of expression casing and extraneous air, it comprises two parts heat interchange factor, at first is the heat loss through convection coefficient of the base-plates surface that caused by fan cooled, and rule of thumb formula calculates, and specifically is expressed as follows
α＝11+3v
Wherein, v is a cross-ventilation speed, can be obtained by experiment measuring, and α equals 44W/ (m
2K).Another part is the heat interchange between casing and the ground, and the coefficient of heat conductivity of iron is 38W/ (mK), and sectional area is 2cm
2, length is 5cm, so the conductive heat flow of iron block is 0.152W/K, equivalence is 14W/ (m for coefficient of heat transfer
2K), A
2Be the surface area of casing, so h
5Equal 58W/ (m
2K); This value can oppositely check according to the casing temperature through infrared heat image instrument measuring and obtain, because this parametric description is the integral body heat transfer of motor, only influences the surface temperature of casing, can not influence the relative value of internal temperature.
2. a kind of temperature rise analytical approach that is used for predicting permasyn morot permanent magnet temperature according to claim 1, it is characterized in that, under a kind of PWM inverter power supply of prediction in the permasyn morot during permanent magnet temperature, analyze that to differ with experiment value be 6.5 ℃.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200910076036A CN101769797A (en) | 2009-01-06 | 2009-01-06 | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200910076036A CN101769797A (en) | 2009-01-06 | 2009-01-06 | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101769797A true CN101769797A (en) | 2010-07-07 |
Family
ID=42502795
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200910076036A Pending CN101769797A (en) | 2009-01-06 | 2009-01-06 | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101769797A (en) |
Cited By (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156000A (en) * | 2011-03-17 | 2011-08-17 | 北京经纬恒润科技有限公司 | Electric motor, electric motor winding temperature detection method and device as well as electric motor winding thermal protection method and device |
| CN103076562A (en) * | 2011-10-25 | 2013-05-01 | 北京大豪科技股份有限公司 | Motor temperature rise measurement method and measurement device |
| CN103344915A (en) * | 2013-06-27 | 2013-10-09 | 深圳市汇川技术股份有限公司 | Synchronous motor temperature rise test system and method |
| CN103400010A (en) * | 2013-08-08 | 2013-11-20 | 英特工程仿真技术（大连）有限公司 | Permanent magnet synchronous motor temperature rise heat dissipation analysis method based on multi-field coupling technology |
| CN103776556A (en) * | 2012-10-18 | 2014-05-07 | 发那科株式会社 | Temperature estimation apparatus for estimating temperature of motor |
| CN103970947A (en) * | 2014-05-05 | 2014-08-06 | 同济大学 | Real-time temperature rise measurement method of electromagnetic system of low-voltage apparatus based on finite element analysis |
| CN104008228A (en) * | 2014-04-28 | 2014-08-27 | 北京交通大学 | Novel method for solving generator stator end winding surface radiating coefficient |
| CN104034445A (en) * | 2014-06-30 | 2014-09-10 | 南车株洲电力机车研究所有限公司 | Permanent magnet temperature on-line detection method and device |
| CN104158463A (en) * | 2014-09-05 | 2014-11-19 | 南车株洲电力机车研究所有限公司 | Rotor temperature monitoring method for permanent magnet synchronous motor and system therefor |
| US8967857B2 (en) * | 2011-07-29 | 2015-03-03 | Fanuc Corporation | Temperature detection device that detects temperature of rotor of motor |
| CN104537154A (en) * | 2014-12-09 | 2015-04-22 | 天津大学 | Permanent magnetic spherical motor temperature rise forecasting method based on three-dimensional equivalent heat network model |
| CN104793141A (en) * | 2015-05-07 | 2015-07-22 | 哈尔滨工业大学 | Method for testing temperature rise of three-phase alternating-current permanent magnet motor |
| CN105676026A (en) * | 2016-01-11 | 2016-06-15 | 山东理工大学 | Analysis method for analyzing temperature rise and loss characteristics of electromagnetic actuator |
| CN105890811A (en) * | 2016-04-08 | 2016-08-24 | 南京深科博业电气股份有限公司 | Motor temperature rise measurement method |
| CN106368814A (en) * | 2016-08-31 | 2017-02-01 | 华中科技大学 | Real-time monitoring method for highest temperature in cylinder head of internal combustion engine |
| CN106446364A (en) * | 2016-09-08 | 2017-02-22 | 东南大学 | Temperature field-thermal circuit direct coupling-based motor heat analysis method |
| CN106570274A (en) * | 2016-11-07 | 2017-04-19 | 广东电网有限责任公司电力科学研究院 | High overload transformer temperature rise calculation method and apparatus |
| 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 |
| CN107490736A (en) * | 2017-08-07 | 2017-12-19 | 北京工业大学 | The method and device that a kind of nondestructive measurement electronic functional module internal temperature and thermal resistance are formed |
| CN108875255A (en) * | 2018-07-04 | 2018-11-23 | 黑龙江科技大学 | Permanent magnet drive motor temperature rise analytical method based on electric car Real-road Driving Cycle |
| 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 |
| CN109490775A (en) * | 2018-10-31 | 2019-03-19 | 中车大连机车车辆有限公司 | Temperature rise model detection method and frame control locomotive based on traction electric machine in frame control locomotive |
| CN109738704A (en) * | 2018-12-15 | 2019-05-10 | 上海工程技术大学 | A kind of external rotor permanent magnet hub motor permanent magnet loss acquisition methods |
| CN109752651A (en) * | 2017-11-03 | 2019-05-14 | 株洲中车时代电气股份有限公司 | A kind of method and system of traction electric machine overtemperature failure predication |
| CN110287581A (en) * | 2019-06-21 | 2019-09-27 | 天津大学 | Based on the comprehensive RC heat transfer model of architectural design parameter flooring radiation energy supplying system |
| CN110309535A (en) * | 2019-05-23 | 2019-10-08 | 中南大学 | A kind of permanent magnet part loss of excitation failure simulation method and fault injector |
| CN110323994A (en) * | 2019-07-24 | 2019-10-11 | 重庆长安新能源汽车科技有限公司 | Method, system, vehicle and the computer readable storage medium of real-time online estimation motor rotor temperature |
| CN110362890A (en) * | 2019-06-28 | 2019-10-22 | 西安交通大学 | A kind of calculation method of meter and the variable-frequency motor iron loss resistance under PWM harmonic condition |
| CN110427664A (en) * | 2019-07-18 | 2019-11-08 | 大连理工大学 | A kind of permanent magnetic coupling temperature field analysis method based on magnetic-thermo coupling |
| CN110532633A (en) * | 2019-08-02 | 2019-12-03 | 大连理工大学 | A kind of permanent magnetic coupling heat analysis method based on equivalent thermal network |
| CN111950797A (en) * | 2020-08-21 | 2020-11-17 | 中国科学院合肥物质科学研究院 | Local temperature prediction method for high-power water-cooling bus with connector |
| CN112134410A (en) * | 2020-09-28 | 2020-12-25 | 揭阳市汇宝昌电器有限公司 | Optimized heat dissipation design method for steady-state temperature field model of permanent magnet motor |
| CN112182869A (en) * | 2020-09-21 | 2021-01-05 | 盖耀辉 | Motor winding equivalent model, establishment method and motor temperature field analysis method |
| CN116848386A (en) * | 2021-02-17 | 2023-10-03 | Abb瑞士股份有限公司 | Method of generating a thermal model of a system comprising an electric machine |
-
2009
- 2009-01-06 CN CN200910076036A patent/CN101769797A/en active Pending
Cited By (59)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156000B (en) * | 2011-03-17 | 2012-12-12 | 北京经纬恒润科技有限公司 | Electric motor, electric motor winding temperature detection method and device as well as electric motor winding thermal protection method and device |
| US9234803B2 (en) | 2011-03-17 | 2016-01-12 | Beijing Jingwei Hirain Technologies Co., Ltd. | Method and device for detecting winding temperature, method and device for thermal protection of a motor |
| CN102156000A (en) * | 2011-03-17 | 2011-08-17 | 北京经纬恒润科技有限公司 | Electric motor, electric motor winding temperature detection method and device as well as electric motor winding thermal protection method and device |
| US8967857B2 (en) * | 2011-07-29 | 2015-03-03 | Fanuc Corporation | Temperature detection device that detects temperature of rotor of motor |
| CN103076562A (en) * | 2011-10-25 | 2013-05-01 | 北京大豪科技股份有限公司 | Motor temperature rise measurement method and measurement device |
| CN103076562B (en) * | 2011-10-25 | 2016-01-20 | 北京大豪科技股份有限公司 | Electric machine temperature rise measuring method and measurement mechanism |
| CN103776556A (en) * | 2012-10-18 | 2014-05-07 | 发那科株式会社 | Temperature estimation apparatus for estimating temperature of motor |
| CN103344915A (en) * | 2013-06-27 | 2013-10-09 | 深圳市汇川技术股份有限公司 | Synchronous motor temperature rise test system and method |
| CN103344915B (en) * | 2013-06-27 | 2016-03-16 | 深圳市汇川技术股份有限公司 | synchronous motor temperature rise test system and method |
| CN103400010A (en) * | 2013-08-08 | 2013-11-20 | 英特工程仿真技术（大连）有限公司 | Permanent magnet synchronous motor temperature rise heat dissipation analysis method based on multi-field coupling technology |
| CN103400010B (en) * | 2013-08-08 | 2016-08-17 | 英特工程仿真技术（大连）有限公司 | A kind of permagnetic synchronous motor temperature rise analysis of heat transmission method based on multi-scenarios method technology |
| CN104008228A (en) * | 2014-04-28 | 2014-08-27 | 北京交通大学 | Novel method for solving generator stator end winding surface radiating coefficient |
| CN103970947A (en) * | 2014-05-05 | 2014-08-06 | 同济大学 | Real-time temperature rise measurement method of electromagnetic system of low-voltage apparatus based on finite element analysis |
| CN104034445B (en) * | 2014-06-30 | 2016-02-17 | 南车株洲电力机车研究所有限公司 | A kind of permanent magnet temperature on-line testing method and device |
| CN104034445A (en) * | 2014-06-30 | 2014-09-10 | 南车株洲电力机车研究所有限公司 | Permanent magnet temperature on-line detection method and device |
| US10374545B2 (en) | 2014-09-05 | 2019-08-06 | Crrc Zhuzhou Electric Locomotive Research Institute Co., Ltd. | Rotor temperature monitoring method and system for permanent magnet synchronous motor |
| CN104158463A (en) * | 2014-09-05 | 2014-11-19 | 南车株洲电力机车研究所有限公司 | Rotor temperature monitoring method for permanent magnet synchronous motor and system therefor |
| CN104158463B (en) * | 2014-09-05 | 2016-03-30 | 南车株洲电力机车研究所有限公司 | A kind of temperature monitoring method and system of permagnetic synchronous motor |
| CN104537154A (en) * | 2014-12-09 | 2015-04-22 | 天津大学 | Permanent magnetic spherical motor temperature rise forecasting method based on three-dimensional equivalent heat network model |
| CN104537154B (en) * | 2014-12-09 | 2017-12-05 | 天津大学 | Permanent magnet spherical motor temperature prediction method based on three-dimensional equivalent ther mal network model |
| CN104793141A (en) * | 2015-05-07 | 2015-07-22 | 哈尔滨工业大学 | Method for testing temperature rise of three-phase alternating-current permanent magnet motor |
| CN104793141B (en) * | 2015-05-07 | 2018-05-22 | 哈尔滨工业大学 | A kind of Thee-phase alternating current permanent-magnetic motor temperature rise test method |
| CN105676026A (en) * | 2016-01-11 | 2016-06-15 | 山东理工大学 | Analysis method for analyzing temperature rise and loss characteristics of electromagnetic actuator |
| CN105890811A (en) * | 2016-04-08 | 2016-08-24 | 南京深科博业电气股份有限公司 | Motor temperature rise measurement method |
| CN105890811B (en) * | 2016-04-08 | 2018-08-21 | 南京深科博业电气股份有限公司 | A kind of motor temperature measurement method |
| CN106368814A (en) * | 2016-08-31 | 2017-02-01 | 华中科技大学 | Real-time monitoring method for highest temperature in cylinder head of internal combustion engine |
| CN106368814B (en) * | 2016-08-31 | 2019-03-05 | 华中科技大学 | The method of real-time of maximum temperature in combustion engine cylinder head |
| WO2018045721A1 (en) * | 2016-09-08 | 2018-03-15 | 东南大学 | Thermal analysis method for motors having directly coupled temperature field and thermal circuit |
| CN106446364A (en) * | 2016-09-08 | 2017-02-22 | 东南大学 | Temperature field-thermal circuit direct coupling-based motor heat analysis method |
| CN106446364B (en) * | 2016-09-08 | 2019-06-25 | 东南大学 | A kind of direct-coupled motor heat analysis method of temperature field-Re Lu |
| CN106570274A (en) * | 2016-11-07 | 2017-04-19 | 广东电网有限责任公司电力科学研究院 | High overload transformer temperature rise calculation method and apparatus |
| CN106570274B (en) * | 2016-11-07 | 2020-11-17 | 广东电网有限责任公司电力科学研究院 | Temperature rise calculation method and device for high-overload transformer |
| CN107291970A (en) * | 2016-12-19 | 2017-10-24 | 沈阳工业大学 | A kind of computational methods in air to air cooling middle size motor temperature field |
| CN107490736A (en) * | 2017-08-07 | 2017-12-19 | 北京工业大学 | The method and device that a kind of nondestructive measurement electronic functional module internal temperature and thermal resistance are formed |
| 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 |
| CN109752651B (en) * | 2017-11-03 | 2021-08-17 | 株洲中车时代电气股份有限公司 | Method and system for predicting over-temperature fault of traction motor |
| CN109752651A (en) * | 2017-11-03 | 2019-05-14 | 株洲中车时代电气股份有限公司 | A kind of method and system of traction electric machine overtemperature failure predication |
| CN108875255A (en) * | 2018-07-04 | 2018-11-23 | 黑龙江科技大学 | Permanent magnet drive motor temperature rise analytical method based on electric car Real-road Driving Cycle |
| CN109490775A (en) * | 2018-10-31 | 2019-03-19 | 中车大连机车车辆有限公司 | Temperature rise model detection method and frame control locomotive based on traction electric machine in frame control locomotive |
| CN109490775B (en) * | 2018-10-31 | 2021-07-06 | 中车大连机车车辆有限公司 | Temperature rise performance detection method based on traction motor in rack-controlled locomotive and rack-controlled locomotive |
| 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 |
| CN109376475B (en) * | 2018-11-27 | 2022-02-08 | 哈尔滨工业大学 | Multi-turn close-wound coil temperature field calculation method combining thermal resistance network method and finite difference method |
| CN109738704A (en) * | 2018-12-15 | 2019-05-10 | 上海工程技术大学 | A kind of external rotor permanent magnet hub motor permanent magnet loss acquisition methods |
| CN110309535A (en) * | 2019-05-23 | 2019-10-08 | 中南大学 | A kind of permanent magnet part loss of excitation failure simulation method and fault injector |
| CN110287581B (en) * | 2019-06-21 | 2023-01-24 | 天津大学 | Floor radiation energy supply system comprehensive RC heat transfer system based on building design parameters |
| CN110287581A (en) * | 2019-06-21 | 2019-09-27 | 天津大学 | Based on the comprehensive RC heat transfer model of architectural design parameter flooring radiation energy supplying system |
| CN110362890A (en) * | 2019-06-28 | 2019-10-22 | 西安交通大学 | A kind of calculation method of meter and the variable-frequency motor iron loss resistance under PWM harmonic condition |
| CN110362890B (en) * | 2019-06-28 | 2021-03-16 | 西安交通大学 | Method for calculating iron loss resistance of variable frequency motor under PWM harmonic condition |
| CN110427664A (en) * | 2019-07-18 | 2019-11-08 | 大连理工大学 | A kind of permanent magnetic coupling temperature field analysis method based on magnetic-thermo coupling |
| CN110323994A (en) * | 2019-07-24 | 2019-10-11 | 重庆长安新能源汽车科技有限公司 | Method, system, vehicle and the computer readable storage medium of real-time online estimation motor rotor temperature |
| CN110532633B (en) * | 2019-08-02 | 2021-06-18 | 大连理工大学 | Permanent magnet coupler thermal analysis method based on equivalent thermal network |
| CN110532633A (en) * | 2019-08-02 | 2019-12-03 | 大连理工大学 | A kind of permanent magnetic coupling heat analysis method based on equivalent thermal network |
| CN111950797A (en) * | 2020-08-21 | 2020-11-17 | 中国科学院合肥物质科学研究院 | Local temperature prediction method for high-power water-cooling bus with connector |
| CN111950797B (en) * | 2020-08-21 | 2023-03-10 | 中国科学院合肥物质科学研究院 | Local temperature prediction method for high-power water-cooling bus with connector |
| CN112182869A (en) * | 2020-09-21 | 2021-01-05 | 盖耀辉 | Motor winding equivalent model, establishment method and motor temperature field analysis method |
| CN112182869B (en) * | 2020-09-21 | 2024-05-31 | 盖耀辉 | Motor winding equivalent model and establishment method and motor temperature field analysis method |
| CN112134410A (en) * | 2020-09-28 | 2020-12-25 | 揭阳市汇宝昌电器有限公司 | Optimized heat dissipation design method for steady-state temperature field model of permanent magnet motor |
| CN116848386A (en) * | 2021-02-17 | 2023-10-03 | Abb瑞士股份有限公司 | Method of generating a thermal model of a system comprising an electric machine |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101769797A (en) | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor | |
| CN109063337B (en) | Multi-field-road coupling simulation calculation method for permanent magnet wind driven generator | |
| CN104182580B (en) | Eddy-current loss analysis method for permanent magnet wind generators | |
| CN104537154B (en) | Permanent magnet spherical motor temperature prediction method based on three-dimensional equivalent ther mal network model | |
| CN103353926B (en) | A kind of motor temperature distribution method of real-time | |
| Le et al. | A novel stator cooling structure for yokeless and segmented armature axial flux machine with heat pipe | |
| Zhu et al. | Cooling system design of a high-speed PMSM based on a coupled fluidic–thermal model | |
| Nasab et al. | Analytical thermal model of natural-convection cooling in axial flux machines | |
| Weili et al. | Calculation of ventilation cooling, three-dimensional electromagnetic fields, and temperature fields of the end region in a large water–hydrogen–hydrogen-cooled turbogenerator | |
| Tong et al. | Research on the airflow and thermal performance in a large forced air-cooled permanent magnet synchronous machine | |
| Li et al. | Research on temperature and braking performance of water-cooled eddy current retarder | |
| Xu et al. | Research on thermal capacity of a high-torque-density direct drive permanent magnet synchronous machine based on a temperature cycling module | |
| Wrobel et al. | Performance analysis and thermal modeling of a high-energy-density prebiased inductor | |
| Li et al. | Loss calculation and thermal analysis of ultra-high speed permanent magnet motor | |
| Wen et al. | Thermal analysis of coupled resonant coils for an electric vehicle wireless charging system | |
| Wang et al. | A review of recent developments in permanent magnet eddy current couplers technology | |
| Li et al. | Numerical calculation of multicoupled fields in large salient synchronous generator | |
| Plait et al. | Low computational cost semi-analytical magnetostatic model for magnetocaloric refrigeration systems | |
| Zheng et al. | Analytical prediction and analysis of electromagnetic-thermal fields in PM eddy current couplings with injected harmonics into magnet shape for torque improvement | |
| Chen et al. | Thermal field and stress analysis of induction motor with stator inter-turn fault | |
| Abubakar et al. | Electromagnetic and thermal analysis of 225 kW high-speed PMSM for centrifugal blower applications | |
| Han et al. | Numerical analysis of end part temperature in the turbogenerator end region with magnetic shield structure under the different operation conditions | |
| Jia et al. | Modeling and analysis of electromagnetic field and temperature field of permanent-magnet synchronous motor for automobiles | |
| Xu et al. | Effect of helium on temperature rise of helium blower drive motor in high-temperature gas-cooled reactor | |
| Cao et al. | Influence of single/dual ventilation path on fluid field and temperature field of HVLSSR-PMSM with air-cooled hybrid ventilation systems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| DD01 | Delivery of document by public notice |
Addressee: Li Hu Document name: Notification of Publication of the Application for Invention |
|
| DD01 | Delivery of document by public notice |
Addressee: Li Hu Document name: Notification of before Expiration of Request of Examination as to Substance |
|
| DD01 | Delivery of document by public notice |
Addressee: Li Hu Document name: Notification that Application Deemed to be Withdrawn |
|
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20100707 |