CN103400010A - Permanent magnet synchronous motor temperature rise heat dissipation analysis method based on multi-field coupling technology - Google Patents
Permanent magnet synchronous motor temperature rise heat dissipation analysis method based on multi-field coupling technology Download PDFInfo
- Publication number
- CN103400010A CN103400010A CN2013103438189A CN201310343818A CN103400010A CN 103400010 A CN103400010 A CN 103400010A CN 2013103438189 A CN2013103438189 A CN 2013103438189A CN 201310343818 A CN201310343818 A CN 201310343818A CN 103400010 A CN103400010 A CN 103400010A
- Authority
- CN
- China
- Prior art keywords
- analysis
- model
- heat
- coupling
- conjugation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
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
Abstract
The invention discloses a permanent magnet synchronous motor temperature rise heat dissipation analysis method based on the multi-field coupling technology. The method comprises the first step of carrying out the analysis of electromagnetic loss by the establishment of an electromagnetic field model of an electric machine to obtain the copper loss on a coil, the iron loss of an iron core and the eddy-current loss of permanent magnet materials, the second step of establishing a fluid solid strong coupling conjugate heat conduction model of the structural temperature field analysis and the fluid heat dissipation analysis of the electric machine, the third step of establishing a bothway weak coupling model of the electromagnetic loss and the fluid solid strong coupling conjugate heat conduction analysis, and the fourth step of achieving the convergence of a single physical model and the convergence of the load transferring between the electromagnetic field model and the fluid solid strong coupling conjugate heat conduction model by the adaptation of a repetitive cycle iterative method, wherein if the condition of the convergence meets the condition set by a user, the coupling is balanced, and the computing process is finished. The permanent magnet synchronous motor temperature rise heat dissipation analysis method integrates the advantages of the strong coupling and the weak coupling, the temperature rise and the heat dissipation process when the electric machine works are precisely simulated in a unified software platform, the stability of the convergence in the process of computing is enhanced, and computing precision and computing efficiency are improved.
Description
Technical field
The present invention relates to simulation technical field, a kind of coupling analysis technology that is applicable to permanent magnet synchronous motor for new energy automobile or other industrial permagnetic synchronous motor electromagnetism heat-dissipating-Re conduction-fluid for radiating heats specifically, the present invention can effectively improve the analysis precision of electric machine temperature rise heat radiation, and optimizes the design of electrical motor analysis.
Background technology
Along with the increasingly severe and environmental problem of non-renewable energy resources crisis to increasingly sharpening that human survival threatens, design and produce the hybrid vehicle that replaces oil or focus that electric automobile has become various countries' concern.Motor driver, as the key distinction of automobile with novel power and traditional fuel-engined vehicle, becomes the key that new-energy automobile successfully designs undoubtedly.Temperature Distribution under motor different operating state not only affects the characteristic of electromagnetic material, also there is vital impact while on the loss of excitation characteristic of permanent magnet, so temperature is controlled and the High Precision Simulation analysis becomes the key that permanent magnet synchronous motor for new energy automobile designs.
More common machine analysis software JMAG and ANSOFT-Maxwell self do not have the CFD software module at present, traditionally, when motor is carried out to temperature and diagnosis refrigeration, generally adopts third party CFD software, and the fluid boundary of liquid coolant is simplified to processing.Concrete, at first liquid coolant is solved under the solid coupled interface temperature conditions of the stream of supposition, obtain approximate coefficient of heat transfer and the Temperature Distribution of the solid coupled interface of stream, then thermoanalytical input boundary condition carries out quasi-static analysis using this result as motor, can't guarantee coupled interface temperature continuity and energy conservation like this, affect the precision of temperature rise analysis of heat transmission, thereby affect the reliability of design of electrical motor.
Summary of the invention
Defect in view of the prior art existence, the objective of the invention is to provide a kind of analytical approach of analyzing permagnetic synchronous motor temperature rise heat radiation based on many coupling techniques, the method can be carried out many coupling analysis to electromagnetic field distribution, electromagnetism heat-dissipating, the conduction of rotor stator solid thermal, the fluid for radiating heat of motor according to actual geometry designs, material parameter, the condition of work of motor, obtain high-precision Temperature Distribution, guarantee the normal operating conditions of motor.
To achieve these goals, technical scheme of the present invention:
Electromagnetism-temperature-fluid field three field analysis functions that a kind of method based on many coupling techniques analysis permagnetic synchronous motor temperature rise heat radiations is integrated, implementation step is as follows:
Step 1: carry out the electromagnetic consumable analysis by setting up the motor electromagnetic field model, obtain copper loss, the iron loss in iron core and the eddy current loss on permanent magnetic material on coil;
Step 2: set up the solid strong coupling conjugation heat conduction model of stream that electric machine structure temperature field analysis and fluid for radiating heat are analyzed, adopting the multi-point constraint method is that the MPC technology is guaranteed the continuity of motor runner solid interface and liquid coolant interface temperature and the conservativeness of hot-fluid, with dirty solid strong coupling conjugation heat, conducts referred to as the conduction of the solid conjugation heat of stream;
Step 3: the two-way weak coupling model of setting up electromagnetic consumable analysis and the solid conjugation heat conduction analysis of stream: the loss that the electromagnetic consumable analysis is obtained shine upon by volume as the heat-dissipating amount and interpolation method is loaded in the solid conjugation heat conduction analysis of stream, and the Temperature Distribution result that will flow solid conjugation heat conduction analysis is shone upon and is delivered in the electromagnetic consumable analysis;
Step 4: utilize the repetitive cycling alternative manner to realize the convergence that the load between the solid conjugation heat conduction model of the convergence of independent physical model and electromagnetic field model and stream is transmitted, when the condition of convergence met the user and imposes a condition, the coupling balance, calculated and finish.Described independent physical model is the solid conjugation heat conduction model of electromagnetic field model and stream.
Compared with prior art, beneficial effect of the present invention: the advantage of the present invention is integrated strong coupling method and weak coupling method, propose many bidirectional coupled technology and the method for brand-new motor electromagnetic heat-dissipating-Re conduction-fluid for radiating heat, realized the multiple physical field high precision temperature rise analysis of heat transmission of motor; The present invention has realized, under unified software platform, the temperature rise in machine operation and heat radiation process are carried out to accurate emulation, and the mode that adopts strong coupling to combine with weak coupling has increased the stability that restrains while calculating, and has improved computational accuracy and counting yield.
The accompanying drawing explanation
Fig. 1: the flow chart of steps that the present invention is based on the permagnetic synchronous motor temperature rise analysis of heat transmission method of many coupling techniques;
Fig. 2: the integrated figure of permagnetic synchronous motor temperature rise analysis of heat transmission methodological function that the present invention is based on many coupling techniques;
Fig. 3: the permagnetic synchronous motor temperature rise analysis of heat transmission method model graph of a relation that the present invention is based on many coupling techniques;
Fig. 4: embodiment of the present invention current of electric distributes and the direction of magnetization schematic diagram;
Fig. 5: the finite element model schematic diagram of embodiment of the present invention permagnetic synchronous motor;
Fig. 6: the magnetic intensity vector figure of embodiment of the present invention permanent magnet;
Fig. 7: the magnetic intensity vector figure of embodiment of the present invention stator;
In figure: 1, fluid, 2, air, 3, air gap, 4, stator, 5, coil, 6, permanent magnet, 7, rotor.
Embodiment
Below in conjunction with accompanying drawing and specific embodiment, further illustrate technical scheme of the present invention:
As shown in Figure 2, electromagnetism-temperature-fluid field three field analysis functions that the present invention is integrated, realized the integrated of under unified software platform electromagnetic consumable analytic function, structure heat conduction analysis function and fluid for radiating heat analytic function.Concrete, as shown in Figure 3, by the multi-point constraint method, realized the strong coupling of the non-matching grid that structure temperature field analysis and fluid for radiating heat are analyzed; By the load that solves a coupling interface, transmit the two-way weak coupling that coupling technique has been realized electromagnetic consumable analysis and the solid conjugation heat conduction analysis of stream.
As shown in Figure 1, to relate generally to step as follows in the present invention:
Step 1: carry out the electromagnetic consumable analysis by setting up the electromagnetic field model, obtain copper loss, the iron loss in iron core and the eddy current loss on permanent magnet material on coil;
Described electromagnetic consumable analysis comprises adopts Edge Finite Element Method to analyze to the electromagnetic potential vector, to Maxwell equation, carries out the finite element discretization processing and adopts method of Lagrange multipliers to guarantee uniqueness of solution.
Described iron loss comprises magnetic hysteresis loss and eddy current loss.
Step 2: set up the solid strong coupling conjugation heat conduction model of stream that electric machine structure temperature field analysis and fluid for radiating heat are analyzed, adopting the multi-point constraint method is that the MPC technology is guaranteed the continuity of motor runner solid interface and liquid coolant interface temperature and the conservativeness of hot-fluid, with dirty solid strong coupling conjugation heat, conducts referred to as the conduction of the solid conjugation heat of stream;
The solid strong coupling conjugation heat conduction model of described foundation stream comprises:
1, set up electric machine structure temperature field analysis create thermal conduction module: according to the real work characteristics of permagnetic synchronous motor, create temperature field analysis and in motor model, between different electric machine assemblies, set up in a different manner the thermo-contact relation, creating the equivalent heat loop;
Described thermo-contact relation comprises:
A, for the heat transfer between stator, rotor, consider the flow effect of air gap inner fluid, adopt experimental formula formula (1) to carry out the setting of the coefficient of heat transfer:
In formula, v is rotor speed, and unit is cm/s, and d is the air-gap separation between rotor, and unit is cm, h
1For the coefficient of heat transfer between rotor, be the amount relevant with rotating speed, rotating speed is larger, and coefficient of heat transfer value also correspondence is larger;
The accounting equation of the contact heat exchange coefficient between b, stator and coil and rotor and permanent magnet all adopts as the formula (2), this contact heat exchange coefficient be with said stator and coil or rotor and permanent magnet between the thermal conductivity λ of the material amount relevant with the gap d between stator and coil or rotor and permanent magnet;
C, for other contact positions in motor model, adopt the multi-point constraint method to realize the effect that the heat free of losses is transmitted; For the position that contacts with air, all adopt and the heat radiation of cross-ventilation heat exchange mode, finally set up complete hot loop.
In above-mentioned thermo-contact relation, set and all adopt penalty function method to realize the setting of the coefficient of heat transfer.
2, create the fluid for radiating heat analysis: adopt Navier-Stokes equation in Finite Element convection cell analysis of heat transmission and energy conservation equation to carry out discrete, obtain the flow equation take pressure, speed, temperature as degree of freedom;
3,, according to fluid and electric machine structure coupled boundary condition, set up the strong coupling condition that electric machine structure temperature field analysis and fluid for radiating heat are analyzed.The strong coupling of described electric machine structure temperature field analysis and fluid for radiating heat analysis adopts the multi-point constraint method to realize the continuity of coupled interface temperature and the balance of hot-fluid, and the subdivision of the non-coordination grid of permission coupled interface, the mesh mapping technology in coupled interface between service routine during described non-coordination grid coupling.
Step 3: the two-way weak coupling model of setting up electromagnetic consumable analysis and the solid conjugation heat conduction analysis of stream: the loss that the electromagnetic consumable analysis is obtained shine upon by volume as the heat-dissipating amount and interpolation method is loaded in the solid conjugation heat conduction analysis of stream, and the Temperature Distribution result that will flow solid conjugation heat conduction analysis is shone upon and is delivered in the electromagnetic consumable analysis;
The two-way weak coupling model of the analysis of described establishment electromagnetic consumable and the solid conjugation heat conduction analysis of stream comprises two aspects: the one, and electromagnetic consumable analyzes the transmission of the solid conjugation heat conduction analysis heat-dissipating amount of stream: electromagnetism heat-dissipating amount is analyzed mesh mapping to the solid conjugation heat conduction analysis grid of stream and flows solid conjugation heat conduction analysis from electromagnetic consumable; The 2nd, the stream impact of conjugation heat conduction analysis result on the electromagnetic material characteristic admittedly: be delivered in the electromagnetic consumable analysis flowing solid conjugation heat conduction analysis result and analyzing grid by mapping techniques, thereby the consideration temperature is on the impact of electromagnetic material characteristic.
Concrete comprises:
1, when relating to non-matching interface coupling, adopt mesh mapping and interpolation technique
Wherein, adopt the map operation of Bucket Search and based on the linear interpolation mode of unit local coordinate during non-matching interface coupling.Linear interpolation mode based on the unit local coordinate comprises the following steps again: find unit and the local coordinate value of impact point on the source item grid, according to the value of each node of source item grid cell, and impact point adopts the method for linear interpolation to try to achieve the numerical value of impact point in the local coordinate value of source item grid cell.
2, from the electromagnetic field model, to the thermal losses data of the solid conjugation heat conduction model of stream, transmit, and to the temperature transmission of electromagnetic field model, be all to carry out the data transmission by network Socket technology from flowing solid conjugation heat conduction model, guaranteed the real time data transmission in operational process between physical field.
Step 4: utilize the repetitive cycling alternative manner to realize the convergence that the load between the solid conjugation heat conduction model of the convergence of independent physical model and electromagnetic field model and stream is transmitted, when the condition of convergence met the user and imposes a condition, the coupling balance, calculated and finish.
The convergence of calculating in order to ensure coupling realizes the solid conjugation heat conduction model coupling of electromagnetic field-stream iteration balance, maneuvering load relaxes factor technique the load mitigation factor is applied on the load vectors that passes between the analytic target model, and realizes by enough iterationses the balance that load is transmitted.
F
I+1=F
i+ α * (F
ext-F
i) formula (3)
In formula (3), F
I+1The load that the i+1 iteration step is delivered to the stand under load model, F
iThe load that last iteration step i is delivered to the stand under load model, F
extBe that current up-to-date should applying outward carried, α relaxes the factor, and numerical value is between (0,1.0); Load transmission convergence control criterion is expressed as follows:
|| F
I+1-F
i||/|| F
i||<ε formula (4)
In formula (4), ε is the decimal of realizing appointment, || * || be vectorial mould.
In order to further illustrate technical solution hereinafter of the present invention, provide a typical embodiment:
Present case is carried out permagnetic synchronous motor temperature rise analysis of heat transmission by electromagnetism-temperature-three of fluid fields coupling technique.In the course of work of motor, heating is inevitable.If the temperature rise of motor is too high, can affect the insulating property of motor, threaten the normal operation of motor, guarantee that therefore the temperature rise of motor reasonably is being necessary in scope.Usually motor has cooling system, such as ventilate, installation refrigeratory etc., carrys out cooling motor, and the temperature of motor internal is controlled in certain scope.Present case uses many coupling techniques to analyze accurately the temperature rise heat radiation process of motor.
Motor is mainly by rotor, stator, and permanent magnet, coil, air forms.Be assumed to be the permagnetic synchronous motor of 4 utmost points, 24 grooves, the material of motor is as shown in table 1, and distribution of current and permanent magnet direction of magnetization are as shown in Figure 4.
Table 1 list of materials
| Assembly | Stator, rotor | Coil | Permanent magnet | Air | Fluid |
| Relative permeability | Vary with |
1 | 1.03 | 1 | — |
| Resistivity (Ω m) | — | 0.001 | — | — | — |
| Coercive force (Am) | — | — | 920000 | — | — |
| Density (Kg/m 3) | 7850 | 8950 | 7500 | 1.293 | 998.2 |
| Thermal conductivity (W/ (mK)) | 23 | 380 | 20 | 0.02624 | 0.6 |
| Specific heat (J/ (KKg)) | 460 | 380 | 460 | 1012 | 4182 |
| Coefficient of viscosity (Kg/ms) | — | — | — | — | 0.001003 |
Be actuated to three-phase current, the concrete numerical value of electric current is as follows:
I
U=A*sin (2 π ft+ θ π/180) formula (5)
I
V=A*sin (2 π ft+ (θ-120) π/180) formula (6)
I
W=A*sin (2 π ft+ (θ+120) π/180) formula (7)
Wherein, I
U, I
V, I
WRespectively U, V, the electric current of W phase; A is the amplitude of electric current, is worth for 4A; F is the frequency of motor, is worth for 60HZ; T is the time; θ is starting phase angle, according to motor characteristic, is made as 240 °.
In present case, to use cooling jacket to carry out cooling for motor, and water jacket has used simplified model to represent with fluid section, and as shown in Figure 5, material is as shown in table 1 for model.
While setting up finite element model, for whole system, done following hypothesis: (1) fluid is incompressible Newtonian fluid, and is in laminar condition; (2) material is considered as isotropy; (3) the motor each several part is considered as the homogeneous heater, and negligible axial is conducted heat, and thinks that motor shaft equates to the temperature of each xsect.Motor finite element model after simplification as shown in Figure 5.
In this case, result is analyzed and obtained to the heat-dissipating analysis of heat transmission of motor by INTESIM software.INTESIM software, as the general emulation of engineering multiple physical field and Optimization Platform, has abundant physical field model, comprises electromagnetic field, temperature field, fluid field etc.; Can solve complicated coupling field problem, comprise hot fluid field etc.
Below by the concise and to the point flow process of telling about motor heat-dissipating analysis of heat transmission:
Step 1: set up the electromagnetic field model
1, import grid assembly.
2, select method of Lagrange multipliers to guarantee uniqueness of solution.
3, according in list of materials, creating material and giving the assembly of appointment.Relative permeability, resistivity and the iron loss attribute of material are set.By material, give assembly, and be assembly designating unit type.
4, the load boundary condition is set.Top and bottom loading seamed edge flux restrained boundary condition, voltage constraint for assembly; Rotating part arranges angular velocity of rotation; Coil loads three-phase current.
5, iron loss being set calculates.The basic frequency of electromagnetic field analysis is set, the computing method of iron loss.
Step 2: set up the solid conjugation heat conduction model of stream that temperature field analysis and fluid for radiating heat are analyzed
1, set up models for temperature field create thermal conduction module
(1) import grid assembly.
(2) create material, thermal conductivity, specific heat and density are set.By material, give corresponding assembly, and be assembly designating unit type.
(3) the load boundary condition is set.Assembly surface is arranged to the surface heat concurrent condition.
(4) set up the thermo-contact relation.According to formula (1), formula (2), can obtain, the coefficient of heat transfer between stator, rotor is 130W/ (m
2K), the coefficient of heat transfer between stator and coil and permanent magnet and rotor is 100W/ (m
2K), by penalty function method, set up the thermo-contact relation.
2, set up the hot fluid field and carry out the fluid for radiating heat analysis
(1) import grid assembly.
(2) create material, thermal conductivity, specific heat, density, coefficient of viscosity are set.By material, give assembly, and be assembly designating unit type.
(3) the load boundary condition is set.For the assembly entrance arranges the constraint of velocity boundary condition; The interior outside of convection cell arranges constraint vx=vy=vz=0m/s; The top and bottom of fluid arrange vz=0m/s; The temperature of inlet porting and the pressure of outlet.
3, the strong coupling condition is set
Strong coupling condition between set temperature field and hot fluid field, realize the transmission of temperature field and hot fluid field coupling interface temperature by the multi-point constraint method.
Step 3: the two-way weak coupling model that creates electromagnetic consumable and the solid conjugation heat conduction analysis of stream
Set up the solid conjugation heat conduction model thermal_fluid of electromagnetic field model emag and stream.Electromagnetic field model emag comprises the electromagnetic field model of setting up in step 1; The solid conjugation heat conduction model thermal_fluid of stream comprises models for temperature field and the hot fluid field model of setting up in step 2.
Next set the weak coupling of electromagnetic field model emag and the solid conjugation heat conduction model thermal_fluid data of stream.Electromagnetic field model emag passes to the solid average heat-dissipating amount of conjugation heat conduction model thermal_fluid of stream, and the solid conjugation heat conduction model of stream is returned to the electromagnetic field model by temperature, thereby sets up two-way weak coupling model.
Step 4: the weak coupling condition of convergence and solve flow process and control to set between program
Electromagnetic field model emag is made as transient analysis, and T.T. is 0.02s, and the time sub-step is 0.001s.The solid conjugation heat conduction model thermal_fluid of stream is steady-state analysis.The coupling time of electromagnetic field model and the solid conjugation heat conduction model of stream is 0.02s; The condition of convergence is: tolerance is made as 1E-5 relatively.
After above setting completes, file is solved.Utilize the repetitive cycling alternative manner, when the condition of convergence met the user and imposes a condition, the coupling balance, calculated and finish.
Aftertreatment
Calculate and finish the corresponding result data of rear importing and check result.The data such as the exportable magnetic field intensity of aftertreatment, magnetic are close, current density, temperature, flow velocity.Fig. 6, Fig. 7 are respectively the magnetic intensity vector figure of permanent magnet, stator, as can be seen from Figure the direction of magnetic field intensity and distribution.
The above; it is only preferably embodiment of the present invention; but protection scope of the present invention is not limited to this; anyly be familiar with those skilled in the art in the technical scope that the present invention discloses; according to technical scheme of the present invention and inventive concept thereof, be equal to replacement or changed, within all should being encompassed in protection scope of the present invention.
Claims (10)
1. the temperature rise of the permagnetic synchronous motor based on many coupling techniques analysis of heat transmission method is characterized in that: comprise following steps:
Step 1: carry out the electromagnetic consumable analysis by setting up the motor electromagnetic field model, obtain copper loss, the iron loss in iron core and the eddy current loss on permanent magnetic material on coil;
Step 2: set up the solid strong coupling conjugation heat conduction model of stream that electric machine structure temperature field analysis and fluid for radiating heat are analyzed, adopting the multi-point constraint method is that the MPC technology is guaranteed the continuity of motor runner solid interface and liquid coolant interface temperature and the conservativeness of hot-fluid, with dirty solid strong coupling conjugation heat, conducts referred to as the conduction of the solid conjugation heat of stream;
Step 3: the two-way weak coupling model of setting up electromagnetic consumable analysis and the solid conjugation heat conduction analysis of stream: the loss that the electromagnetic consumable analysis is obtained shine upon by volume as the heat-dissipating amount and interpolation method is loaded in the solid conjugation heat conduction analysis of stream, and the Temperature Distribution result that will flow solid conjugation heat conduction analysis is shone upon and is delivered in the electromagnetic consumable analysis;
Step 4: utilize the repetitive cycling alternative manner to realize the convergence of independent physical model and the convergence that the load between electromagnetic field model and the solid conjugation heat conduction model of stream is transmitted, when the condition of convergence meets the user and imposes a condition, the coupling balance, calculate and finish, wherein said independent physical model is the solid conjugation heat conduction model of electromagnetic field model and stream.
2. analytical approach according to claim 1, it is characterized in that: described electromagnetic consumable analysis comprises adopts Edge Finite Element Method to analyze to the electromagnetic potential vector, to Maxwell equation, carries out the finite element discretization processing and adopts method of Lagrange multipliers to guarantee uniqueness of solution.
3. analytical approach according to claim 1 is characterized in that: the solid strong coupling conjugation heat conduction model of described foundation stream comprises:
ⅰ, set up the electric machine structure temperature field analysis and create thermal conduction module: according to the real work characteristics of permagnetic synchronous motor, create temperature field analysis and in motor model, between different electric machine assemblies, set up in a different manner the thermo-contact relation, creating the equivalent heat loop;
ⅱ, create the fluid for radiating heat analysis: adopt Navier-Stokes equation in Finite Element convection cell analysis of heat transmission and energy conservation equation to carry out discrete, obtain the flow equation take pressure, speed, temperature as degree of freedom;
ⅲ, according to fluid and electric machine structure coupled boundary condition, set up the strong coupling condition that electric machine structure temperature field analysis and fluid for radiating heat are analyzed.
4. analytical approach according to claim 3, it is characterized in that: described thermo-contact relation comprises:
A, for the heat transfer between stator, rotor, consider the flow effect of air gap inner fluid, adopt experimental formula formula (1) to carry out the setting of the coefficient of heat transfer:
In formula, v is rotor speed, and unit is cm/s, and d is the air-gap separation between rotor, and its unit is cm, h
1For the coefficient of heat transfer between stator rotor;
Contact heat exchange coefficient between b, stator and coil and rotor and permanent magnet all uses accounting equation as the formula (2):
In formula, λ is the thermal conductivity of material between stator and coil or rotor and permanent magnet, and d is the gap between stator and coil or rotor and permanent magnet;
C, for other contact positions in motor model, adopt the multi-point constraint method to realize the effect that the heat free of losses is transmitted; For the position that contacts with air, all adopt and the heat radiation of cross-ventilation heat exchange mode, set up complete hot loop.
5. analytical approach according to claim 4, is characterized in that: in described thermo-contact relation, adopt penalty function method to realize the setting of the coefficient of heat transfer.
6. analytical approach according to claim 3, it is characterized in that: the strong coupling of described electric machine structure temperature field analysis and fluid for radiating heat analysis adopts the multi-point constraint method to realize the continuity of coupled interface temperature and the balance of hot-fluid, and the subdivision of the non-coordination grid of permission coupled interface, the mesh mapping technology in coupled interface between service routine during described non-coordination grid coupling.
7. analytical approach according to claim 1 is characterized in that: the two-way weak coupling model of the analysis of described establishment electromagnetic consumable and the solid conjugation heat conduction analysis of stream adopts mesh mapping and interpolation technique when relating to non-matching interface coupling.
8. analytical approach according to claim 7 is characterized in that: during described non-matching interface coupling, mesh mapping and interpolation technique are the map operation of Bucket Search and based on the linear interpolation mode of unit local coordinate; Wherein the linear interpolation mode based on the unit local coordinate comprises the following steps again: find unit and the local coordinate value of impact point on the source item grid, according to the value of each node of source item grid cell, and impact point adopts the method for linear interpolation to try to achieve the numerical value of impact point in the local coordinate value of source item grid cell.
9. analytical approach according to claim 1, it is characterized in that: described from the electromagnetic field model to the solid conjugation heat conduction model of stream the transmission of thermal losses data and from flowing solid conjugation heat conduction model, to the temperature transmission of electromagnetic field model, be all to carry out the data transmission by network Socket technology, to have guaranteed the real time data transmission in operational process between physical field.
10. analytical approach according to claim 1, it is characterized in that: described realization convergence: maneuvering load relaxes factor technique the load mitigation factor is applied on the load vectors that passes between the analytic target model, and realizes by enough iterationses the balance that load is transmitted:
F
I+1=F
i+ α * (F
ext-F
i) formula (3)
In formula (3), F
I+1The load that the i+1 iteration step is delivered to the stand under load model, F
iThe load that last iteration step i is delivered to the stand under load model, F
extBe that current up-to-date should applying outward carried, α relaxes the factor, and numerical value is between (0,1.0); Load transmission convergence control criterion is expressed as follows:
|| F
I+1-F
i||/|| F
i||<ε formula (4)
In formula (4), ε is the decimal of realizing appointment, || * || be vectorial mould.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310343818.9A CN103400010B (en) | 2013-08-08 | 2013-08-08 | A kind of permagnetic synchronous motor temperature rise analysis of heat transmission method based on multi-scenarios method technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310343818.9A CN103400010B (en) | 2013-08-08 | 2013-08-08 | A kind of permagnetic synchronous motor temperature rise analysis of heat transmission method based on multi-scenarios method technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103400010A true CN103400010A (en) | 2013-11-20 |
| CN103400010B CN103400010B (en) | 2016-08-17 |
Family
ID=49563638
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310343818.9A Active CN103400010B (en) | 2013-08-08 | 2013-08-08 | A kind of permagnetic synchronous motor temperature rise analysis of heat transmission method based on multi-scenarios method technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103400010B (en) |
Cited By (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103617367A (en) * | 2013-12-06 | 2014-03-05 | 三峡大学 | Irregular mesh mapping method used in electromagnetic field-flow field-temperature field coupling calculation |
| CN104182580A (en) * | 2014-08-21 | 2014-12-03 | 内蒙古工业大学 | Eddy-current loss analysis method for permanent magnet wind generators |
| CN104537193A (en) * | 2015-01-21 | 2015-04-22 | 英特工程仿真技术（大连）有限公司 | System for realizing coupling analysis of strength of multi-physic field under unified platform |
| CN104537154A (en) * | 2014-12-09 | 2015-04-22 | 天津大学 | Permanent magnetic spherical motor temperature rise forecasting method based on three-dimensional equivalent heat network model |
| CN105354354A (en) * | 2015-09-28 | 2016-02-24 | 沈阳航空航天大学 | Method for calculating temperature field of main driving motor of electric aircraft |
| CN108509726A (en) * | 2018-03-30 | 2018-09-07 | 河北工业大学 | Electro spindle heat/optimizing thermal solution analysis method based on heat flow piercement emulation |
| CN108875235A (en) * | 2018-06-28 | 2018-11-23 | 奇瑞汽车股份有限公司 | The three-dimensional steady state heat transfer property analysis method and device of the water jacket of automobile engine |
| CN109063311A (en) * | 2018-07-26 | 2018-12-21 | 黑龙江科技大学 | A kind of permanent magnet synchronous motor temperature rise analysis of heat transmission method |
| 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 |
| CN109344557A (en) * | 2018-11-30 | 2019-02-15 | 安徽江淮汽车集团股份有限公司 | A kind of analysis method and system of automobile hot property |
| CN109861172A (en) * | 2019-03-29 | 2019-06-07 | 北京经纬恒润科技有限公司 | Motor overheat protection method and device |
| CN105354354B (en) * | 2015-09-28 | 2019-07-16 | 沈阳航空航天大学 | A kind of electric airplane main drive motor Calculation Method of Temperature Field |
| CN110032810A (en) * | 2019-04-17 | 2019-07-19 | 燕山大学 | Copper pipe rolls cooling procedure emulation mode, equipment and computer readable storage medium |
| CN110555289A (en) * | 2019-09-27 | 2019-12-10 | 哈尔滨理工大学 | Cloud computing-based motor stator winding multi-physical field coupling computing platform and computing method |
| CN110702272A (en) * | 2019-09-24 | 2020-01-17 | 中国第一汽车股份有限公司 | Method and device for estimating temperature of motor rotor, motor and storage medium |
| CN111104764A (en) * | 2020-01-21 | 2020-05-05 | 湖南科技大学 | Structured grid optimization division method for alternating current motor rotor conducting bar thermal analysis model |
| CN111737779A (en) * | 2020-06-23 | 2020-10-02 | 哈尔滨工业大学 | Simulation method for barium titanate ceramic PTC thermal management structure of vehicle-mounted air conditioner |
| CN111881611A (en) * | 2020-07-31 | 2020-11-03 | 珠海格力电器股份有限公司 | Simulation processing method and device for physical field of motor |
| CN113591264A (en) * | 2021-04-27 | 2021-11-02 | 平高集团有限公司 | Temperature calculation method and device for high-voltage switch |
| CN115329624A (en) * | 2022-07-20 | 2022-11-11 | 重庆金康动力新能源有限公司 | Motor rotor temperature rise simulation method |
| CN116432484A (en) * | 2023-03-23 | 2023-07-14 | 中国科学院高能物理研究所 | Multi-physical field coupling numerical simulation method for cooling quench process of liquid helium-free superconducting cavity |
| CN117669337A (en) * | 2024-02-01 | 2024-03-08 | 苏州保邦电气有限公司 | Design method of ultra-high-speed permanent magnet magnetic suspension motor |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101769797A (en) * | 2009-01-06 | 2010-07-07 | 李虎 | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor |
| CN102194034A (en) * | 2011-06-29 | 2011-09-21 | 福州大学 | Electromagnetic valve temperature rise computational analysis method |
| CN102592016A (en) * | 2011-12-31 | 2012-07-18 | 英特工程仿真技术（大连）有限公司 | Engineering multiphysics coupling analysis method |
| US20120239245A1 (en) * | 2009-12-04 | 2012-09-20 | Hiroshi Tamura | Numerical analysis method of a vehicle drive |
-
2013
- 2013-08-08 CN CN201310343818.9A patent/CN103400010B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101769797A (en) * | 2009-01-06 | 2010-07-07 | 李虎 | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor |
| US20120239245A1 (en) * | 2009-12-04 | 2012-09-20 | Hiroshi Tamura | Numerical analysis method of a vehicle drive |
| CN102194034A (en) * | 2011-06-29 | 2011-09-21 | 福州大学 | Electromagnetic valve temperature rise computational analysis method |
| CN102592016A (en) * | 2011-12-31 | 2012-07-18 | 英特工程仿真技术（大连）有限公司 | Engineering multiphysics coupling analysis method |
Non-Patent Citations (3)
| Title |
|---|
| 孔晓光: "高速永磁电机定子损耗和温升研究", <工程科技Ⅱ辑> * |
| 张琪: "高密度车用永磁电机流固耦合传热仿真分析", 《电机与控制应用》 * |
| 饶明忠 等: "电磁场计算中的棱边有限元法", 《中国电机工程学报》 * |
Cited By (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103617367A (en) * | 2013-12-06 | 2014-03-05 | 三峡大学 | Irregular mesh mapping method used in electromagnetic field-flow field-temperature field coupling calculation |
| CN103617367B (en) * | 2013-12-06 | 2016-08-17 | 三峡大学 | Special-shaped mesh mapping method in electromagnetic field-flow field-temperature field coupling calculating |
| CN104182580A (en) * | 2014-08-21 | 2014-12-03 | 内蒙古工业大学 | Eddy-current loss analysis method for permanent magnet wind generators |
| CN104182580B (en) * | 2014-08-21 | 2017-01-25 | 内蒙古工业大学 | Eddy-current loss analysis method for permanent magnet wind generators |
| 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 |
| CN104537193B (en) * | 2015-01-21 | 2017-12-29 | 英特工程仿真技术（大连）有限公司 | The system that multiple physical field power coupling analysis is realized under the unified platform |
| CN104537193A (en) * | 2015-01-21 | 2015-04-22 | 英特工程仿真技术（大连）有限公司 | System for realizing coupling analysis of strength of multi-physic field under unified platform |
| 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 |
| CN108509726A (en) * | 2018-03-30 | 2018-09-07 | 河北工业大学 | Electro spindle heat/optimizing thermal solution analysis method based on heat flow piercement emulation |
| CN108509726B (en) * | 2018-03-30 | 2021-08-10 | 河北工业大学 | Electric spindle heat generation/heat dissipation optimization analysis method based on heat-fluid-solid coupling simulation |
| CN108875235A (en) * | 2018-06-28 | 2018-11-23 | 奇瑞汽车股份有限公司 | The three-dimensional steady state heat transfer property analysis method and device of the water jacket of automobile engine |
| CN109063311A (en) * | 2018-07-26 | 2018-12-21 | 黑龙江科技大学 | A kind of permanent magnet synchronous motor temperature rise analysis of heat transmission method |
| 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 |
| CN109063337B (en) * | 2018-08-03 | 2022-03-29 | 内蒙古工业大学 | Multi-field-road coupling simulation calculation method for permanent magnet wind driven generator |
| CN109344557A (en) * | 2018-11-30 | 2019-02-15 | 安徽江淮汽车集团股份有限公司 | A kind of analysis method and system of automobile hot property |
| CN109861172A (en) * | 2019-03-29 | 2019-06-07 | 北京经纬恒润科技有限公司 | Motor overheat protection method and device |
| CN110032810A (en) * | 2019-04-17 | 2019-07-19 | 燕山大学 | Copper pipe rolls cooling procedure emulation mode, equipment and computer readable storage medium |
| CN110702272A (en) * | 2019-09-24 | 2020-01-17 | 中国第一汽车股份有限公司 | Method and device for estimating temperature of motor rotor, motor and storage medium |
| CN110702272B (en) * | 2019-09-24 | 2021-11-05 | 中国第一汽车股份有限公司 | Method and device for estimating temperature of motor rotor, motor and storage medium |
| 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 |
| CN111104764A (en) * | 2020-01-21 | 2020-05-05 | 湖南科技大学 | Structured grid optimization division method for alternating current motor rotor conducting bar thermal analysis model |
| CN111737779A (en) * | 2020-06-23 | 2020-10-02 | 哈尔滨工业大学 | Simulation method for barium titanate ceramic PTC thermal management structure of vehicle-mounted air conditioner |
| CN111881611B (en) * | 2020-07-31 | 2023-07-14 | 珠海格力电器股份有限公司 | Simulation processing method and device for motor physical field |
| CN111881611A (en) * | 2020-07-31 | 2020-11-03 | 珠海格力电器股份有限公司 | Simulation processing method and device for physical field of motor |
| CN113591264A (en) * | 2021-04-27 | 2021-11-02 | 平高集团有限公司 | Temperature calculation method and device for high-voltage switch |
| CN113591264B (en) * | 2021-04-27 | 2024-03-26 | 平高集团有限公司 | Temperature calculation method and device for high-voltage switch |
| CN115329624B (en) * | 2022-07-20 | 2023-08-25 | 重庆金康动力新能源有限公司 | Motor rotor temperature rise simulation method |
| CN115329624A (en) * | 2022-07-20 | 2022-11-11 | 重庆金康动力新能源有限公司 | Motor rotor temperature rise simulation method |
| CN116432484A (en) * | 2023-03-23 | 2023-07-14 | 中国科学院高能物理研究所 | Multi-physical field coupling numerical simulation method for cooling quench process of liquid helium-free superconducting cavity |
| CN116432484B (en) * | 2023-03-23 | 2024-04-16 | 中国科学院高能物理研究所 | Multi-physical field coupling numerical simulation method for cooling quench process of liquid helium-free superconducting cavity |
| CN117669337A (en) * | 2024-02-01 | 2024-03-08 | 苏州保邦电气有限公司 | Design method of ultra-high-speed permanent magnet magnetic suspension motor |
| CN117669337B (en) * | 2024-02-01 | 2024-04-09 | 苏州保邦电气有限公司 | Design method of ultra-high-speed permanent magnet magnetic suspension motor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103400010B (en) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103400010A (en) | Permanent magnet synchronous motor temperature rise heat dissipation analysis method based on multi-field coupling technology | |
| CN109063311A (en) | A kind of permanent magnet synchronous motor temperature rise analysis of heat transmission method | |
| CN106446364B (en) | A kind of direct-coupled motor heat analysis method of temperature field-Re Lu | |
| CN110412462B (en) | Transient temperature prediction method for permanent magnet synchronous motor for aviation | |
| Chen et al. | Analysis of temperature field and water cooling of outer rotor in-wheel motor for electric vehicle | |
| CN105577069B (en) | The active thermal optimal control method and device of a kind of drive system of electric automobile | |
| CN101769797A (en) | Temperature rise analytical method for predicting temperature of permanent magnet in permanent magnet synchronous motor | |
| Wang et al. | Ultrafast steady-state multiphysics model for PM and synchronous reluctance machines | |
| CN103413013B (en) | A kind of in-situ heat simulating analysis of battery cell in lithium ion battery pack system | |
| CN104915472B (en) | Engine-cooling system optimization Simulation computational methods | |
| CN104537154A (en) | Permanent magnetic spherical motor temperature rise forecasting method based on three-dimensional equivalent heat network model | |
| CN104182580A (en) | Eddy-current loss analysis method for permanent magnet wind generators | |
| Liu et al. | Exploration of interactive thermal influence characteristics of power and air conditioning system based on 1D/3D coupling calculation in electric vehicle underhood | |
| Hao et al. | Transient thermal model of a permanent magnet synchronous planar motor considering spreading thermal resistance | |
| Sequeira et al. | Validation and parametric investigations of an internal permanent magnet motor using a lumped parameter thermal model | |
| He et al. | Research on temperature rise characteristics of vehicle motors under bench working condition | |
| Wahl et al. | Active temperature control of electric drivetrains for efficiency increase | |
| Wang et al. | Design and optimization of driving motor cooling water pipeline structure based on a comprehensive evaluation method and CNN-PSO | |
| CN106096157A (en) | A kind of axial ventilation permagnetic synchronous motor Calculation Method of Temperature Field | |
| Xie et al. | An efficient multidisciplinary design research for the integrated low speed permanent magnet motor system based on analytical and numerical hybrid analysis | |
| Salar Elahi et al. | Analytical technique for determination of toroidal plasma displacement | |
| Tan et al. | CFD simulation and investigation on flow field characteristics and temperature predictions of an electric generator unit | |
| Li et al. | Bidirectional electromagnetic–thermal coupling analysis for permanent magnet traction motors under complex operating conditions | |
| Nituca | Thermal Analisys for a Double Sided Linear Induction Motor | |
| He et al. | Methods for temperature estimation and monitoring of permanent magnet: a technology review and future trends |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |