Study on temperature influence factors and electromagnetic heat coupling of in-wheel motor for micro-electric vehicle

With the energy drying up and worsening environment, the new energy vehicles have obtained more attentions, and electric vehicles have an extraordinary development in many countries of the world because of their advantages of emission-free, low noise, energy saving and high efficiency [1]. Permanent magnet brushless direct current motor (PMBLDCM) are often used to drive electric vehicles, which are regarded as the in-wheel motors of electric vehicles, especially the short distance transportation micro-electric vehicle. Because PMBLDCM possesses the merits of high power density, light weight, simple structure, small volume, high starting torque, strong overload capacity, good maneuverability and controllability, high reliability, low noise etc [2, 3]. In this paper, the inwheel motor of micro-electric vehicle is PMBLDCM, which is installed in the wheel to improve the utilization rate of the body space and abandon the traditional clutch, gear, Transmission Bridge and other mechanical components. In-wheel motor has the advantages of reducing the vehicle weight, decreasing the mechanical losses and improving the driving sensitivity. In the micro-electric vehicle, the in-wheel motor possesses the narrow installation space, complex working condition, high power densities, no cooling fan and bad cooling environment and so on, which not only easily causes the temperature too high, but also reduces the work performance and service life of the in-wheel motor. Therefore, it is vital significance to study thermal loss mechanism and temperature influence factors of in-wheel motor, in order to meet the development requirements of micro-electric vehicle in efficiency, reliability and safety. Owing to complex thermal model structure, difficulty of dealing with windings and determining of convection heat transfer coefficient and limitation of ventilation cooling system, thermal loss mechanism and temperature field influence factors of in-wheel motor in micro-electric vehicle are calculated and analyzed difficultly. Furthermore, it is more difficult to master the inner temperature change rules when there is thermal breakdown in the inwheel motor. In order to accurately calculate inner thermal loss and temperature field distribution of in-wheel motor, a lot of researches have been conducted and many practical feasible methods were given. Referring to related articles, the methods of thermal analysis of motor mainly involve thermal parameters method, network topologic method, equivalent thermal network, the finite element method and the equivalent heat road method, etc [4]. Li and Yang [5] adopted the finite element method to study the static and transient static temperature field distribution of the pressure level panel ultrasonic motor, where mechanical vibration and friction loss were calculated, and the influence of the high temperature was analyzed by the equivalent insulation, but it lacked the analysis of load influence on temperature field. Based on coupling physical field in generator, the mathematical model of temperature field was established, and the electromagnetic field and temperature field distribution were calculated and analyzed by the numerical method, but this method lacked the test [6]. Wang et al [7] applied the equivalent heat road method to calculate the thermal field distribution of motor, but the lack of test. All-around three-dimensional finite element model of temperature field was established in the claw motor of vehicle, the boundary condition problems of temperature field were handled effectively, All-around static and transient state temperature field under rated load were calculated, and the test results were in accord with the calculation results, but this method lacked the analysis of speed influence on temperature field [8]. The mathematical model of thermal loss was established, the temperature field distribution was calculated, the influence of the different voltages and frequencies on temperature distribution was simulated and analyzed, and the measurement results were in accord with the calculation results, but the lack of analysis of transient temperature field [9]. According to above analysis, the object of the present study is thermal loss and temperature field of inwheel motor in micro-electric vehicle. In the light of the complex structure, complex heat exchange, high power density, high electromagnetic load and high thermal load in the in-wheel motor, some key factors are studied for simplicities calculation, saving time and improving accuracy. A finite element method is used to simulate and analyze thermal loss and temperature field influence factors in this paper. The mathematical model of temperature field is established, thermal loss is calculated, and the transient state temperature fields are studied under its nominal rated load. The influence of the load and revolving speed on the temperature rise of the in-wheel motor is calculated and discussed. Test bench of in-wheel motor is established, which is used to verify the rationality and effectiveness of simulation results.


Introduction
With the energy drying up and worsening environment, the new energy vehicles have obtained more attentions, and electric vehicles have an extraordinary development in many countries of the world because of their advantages of emission-free, low noise, energy saving and high efficiency [1].Permanent magnet brushless direct current motor (PMBLDCM) are often used to drive electric vehicles, which are regarded as the in-wheel motors of electric vehicles, especially the short distance transportation micro-electric vehicle.Because PMBLDCM possesses the merits of high power density, light weight, simple structure, small volume, high starting torque, strong overload capacity, good maneuverability and controllability, high reliability, low noise etc [2,3].In this paper, the inwheel motor of micro-electric vehicle is PMBLDCM, which is installed in the wheel to improve the utilization rate of the body space and abandon the traditional clutch, gear, Transmission Bridge and other mechanical components.In-wheel motor has the advantages of reducing the vehicle weight, decreasing the mechanical losses and improving the driving sensitivity.In the micro-electric vehicle, the in-wheel motor possesses the narrow installation space, complex working condition, high power densities, no cooling fan and bad cooling environment and so on, which not only easily causes the temperature too high, but also reduces the work performance and service life of the in-wheel motor.Therefore, it is vital significance to study thermal loss mechanism and temperature influence factors of in-wheel motor, in order to meet the development requirements of micro-electric vehicle in efficiency, reliability and safety.
Owing to complex thermal model structure, difficulty of dealing with windings and determining of convection heat transfer coefficient and limitation of ventilation cooling system, thermal loss mechanism and temperature field influence factors of in-wheel motor in micro-electric vehicle are calculated and analyzed difficultly.Furthermore, it is more difficult to master the inner temperature change rules when there is thermal breakdown in the inwheel motor.In order to accurately calculate inner thermal loss and temperature field distribution of in-wheel motor, a lot of researches have been conducted and many practical feasible methods were given.Referring to related articles, the methods of thermal analysis of motor mainly involve thermal parameters method, network topologic method, equivalent thermal network, the finite element method and the equivalent heat road method, etc [4].Li and Yang [5] adopted the finite element method to study the static and transient static temperature field distribution of the pres-sure level panel ultrasonic motor, where mechanical vibration and friction loss were calculated, and the influence of the high temperature was analyzed by the equivalent insulation, but it lacked the analysis of load influence on temperature field.Based on coupling physical field in generator, the mathematical model of temperature field was established, and the electromagnetic field and temperature field distribution were calculated and analyzed by the numerical method, but this method lacked the test [6].Wang et al [7] applied the equivalent heat road method to calculate the thermal field distribution of motor, but the lack of test.All-around three-dimensional finite element model of temperature field was established in the claw motor of vehicle, the boundary condition problems of temperature field were handled effectively, All-around static and transient state temperature field under rated load were calculated, and the test results were in accord with the calculation results, but this method lacked the analysis of speed influence on temperature field [8].The mathematical model of thermal loss was established, the temperature field distribution was calculated, the influence of the different voltages and frequencies on temperature distribution was simulated and analyzed, and the measurement results were in accord with the calculation results, but the lack of analysis of transient temperature field [9].
According to above analysis, the object of the present study is thermal loss and temperature field of inwheel motor in micro-electric vehicle.In the light of the complex structure, complex heat exchange, high power density, high electromagnetic load and high thermal load in the in-wheel motor, some key factors are studied for simplicities calculation, saving time and improving accuracy.A finite element method is used to simulate and analyze thermal loss and temperature field influence factors in this paper.The mathematical model of temperature field is established, thermal loss is calculated, and the transient state temperature fields are studied under its nominal rated load.The influence of the load and revolving speed on the temperature rise of the in-wheel motor is calculated and discussed.Test bench of in-wheel motor is established, which is used to verify the rationality and effectiveness of simulation results.

Method and analysis
When thermal loss is analyzed, it is necessary to calculate the thermal loss of the motor, which is regarded as heat source of the motor.The total loss of the in-wheel motor is P z which includes winding copper loss, iron core loss, permanent magnet eddy loss, mechanical loss and stray loss [8][9], as shown in Eq. ( 1): here P Cu is winding copper loss, P Fe is iron core loss, P me is permanent magnet eddy loss, P m is mechanical loss, and P s is stray loss.In engineering field, mechanical loss and stray loss are generally determined by the motor structure, engineering experience and computational formula.

Winding copper loss
Losses caused by leakage magnetic field and leakage electric field are very small, which are easily neglected.According to the Joule's and Lenz's law, winding copper loss is equal to the product of the resistance and square of the winding current.The in-wheel motor adopts three phase winding in this paper, and the calculation formula is expressed in Eq. (2).When the motor is the transient thermal model, the inner temperature of the motor changes constantly, and the resistance value also changes constantly, therefore winding copper loss of the motor varies with time [10].The resistance temperature formula is expressed in Eq. ( 3): ( ) where I is the current of each winding, R τ is the resistance of each winding at τ, R 0 is the resistance of each winding at 0, λ is the temperature coefficient of the resistance.

Iron core loss
Alternating iron core loss and rotating iron core loss are produced because of the existence of the alternating magnetic field and rotating magnetic field.In order to calculate the iron core loss precisely, flux density change law of the iron core at any given time must be ascertained.Currently, a theory presented by Bertotti is wildly used in the engineering [7,11].Based on the different heat mechanism, iron core loss can be solved by separation and superposition, iron core loss consists of eddy current loss, magnetic hysteresis loss and additional losses, which coexist in the iron core.The machining accuracy of the yoke and tooth and the difference of the flux density are ignored, and iron core loss can be expressed as Eq. ( 4): Fe h e ex where P Fe is total iron core loss of unit mass, P h is magnetic hysteresis loss of unit mass, P e is eddy current loss of unit mass, and P ex is additional loss of unit mass.
According to the harmonic analysis principle, any given point flux density wave in the in-wheel motor can be decomposed into a series of elliptic harmonic flux density vector.The revolving iron core losses can be equivalent to the losses produced by two quadrature magnetic fields, for each harmonic component.Magnetic hysteresis loss of unit mass is expressed as Eq. ( 5), eddy-current loss of unit mass is expressed as Eq. ( 6), and additional losses of unit mass is expressed as Eq.(7).Based on the finite element method to calculate the flux density vector waveform of each unit within a cycle, the unit mass core loss of each unit was calculated, which then multiplied by the mass to obtain the core loss of this unit.The total core loss equals to the sum of the unit core loss of all parts: ( ) where k h is hysteresis loss coefficient, f is actual core flux frequency, B i is the long axis flux density, B j is the short axis flux density, k e is eddy current loss coefficient, k ex is additional loss coefficient, and T ex is flux density cycle.

Eddy current loss of permanent magnet
Permanent magnet possesses higher conductivity, which can induce the eddy.When the outer magnetic field changes, the eddy can produce the loss defined as eddy current loss of permanent magnet.Permanent magnet of inwheel motor adopts Nd-Fe-B, which possesses the satisfactory magnetic flux density performances, but its temperature coefficient is high and its heat resistance is bad.The permanent magnet is easy to weaken magnetism when the temperature of permanent magnet becomes higher to a certain extent, which can reduce the overall performance of the motor.For these reasons, it has practical significance to calculate and analyze the eddy current loss of permanent magnet accurately.According to the magnetic induction law, the induced electromotive force and induced current which are in eddy shape can be generated in the permanent magnet.Eddy current loss of permanent magnet is expressed as Eq. ( 8): where L a is the axial length of permanent magnet, L b is average width of permanent magnet, V is volume of permanent magnet, k me is proportion constant of electromotive force, f me is alternating frequency of magnetic field, B me is the biggest flux density of permanent magnet, and 1 ρ is resistivity of permanent magnet.

Mechanical loss
Mechanical loss consists of bearing friction loss and ventilation loss.Bearing friction loss is expressed as Eq. ( 9), which is mainly affected by the machining accuracy, assembly quality, bearings quality, lubricating grease and other aspects [10,11].Ventilation loss is expressed as Eq.(10), which is mainly related to the motor structure, revolving speed, air gap, and stacked thick: where P f is the bearing friction loss, F is the bearing loads, d is the center diameter of bearing, v is the peripheral speed of bearing center, P air is the ventilation loss, κ f is ventilation friction loss coefficient, ρ air is the air density, l is the cylinder length, ω is the revolving angular velocity, and r is the cylinder radius.

Stray loss
It is very difficult to calculate the stray loss of the in-wheel motor accurately, which is generally determined by the motor structure, engineering experience and computational formula [12,13].As the load increases, the wheel motor current increases, and the stray loss increases approximately with the square of the current.When the stator phase current is I 1 , stray loss of the in-wheel motor can be approximately expressed as Eq. ( 11): where P sN is the stray loss when the output power of the inwheel motor is in rated power, and I N is the rated current of the in-wheel motor.

Mathematical model of temperature field
According to energy conservation law and the basic principle of heat transfer, coefficient of thermal conductivity is constant for isotropic media [10].The corresponding boundary conditions of the solution domain are established, which can be divided into three categories, which include the first boundary condition, the second boundary condition and the third boundary condition.In the Cartesian coordinate system, the transient temperature field of the in-wheel motor can be obtained by the heat conduction controlling differential equations, as shown in Eq. (12).Eq. ( 12) expresses the mathematical model of transient temperature field for in-wheel motor, which includes the heat conduction controlling differential equations and the boundary conditions choice equations: ( ) where T is the temperature varying with time t, κ x , κ y and κ z are respectively the thermal conductivity of the material along the x, y, and z directions, and x y z κ κ κ = = , q s is the inner generated heat in per unit volume and per unit time, ρ is density of the motor material, and c is specific heat of the material.In static state temperature field, the temperature T does not vary with time t is time.S 1 is the boundary surface of the first class, T 0 is a known temperature during the static state heat conduction process which can also be regarded as the temperature during transient heat conduction process.S 2 is the boundary surface of the second class, n is the unit vector of boundary surface outer normal direction, κ n is the thermal conductivity of boundary surface outer normal direction, and q 0 is boundary surface heat flux.S 3 is the third boundary surface that is convective heat transfer surface, T ev is surrounding fluid temperature at the convective heat transfer boundary surface S 3 , and β is the surface heat convection coefficient at the boundary surface.

Equivalent winding
Thermal model of the stator equivalent winding was established.The whole winding was regarded as an equivalent part, and thermal conductivity of the equivalent winding was calculated.Some presumptions were given [10,11]: (a) the slot wire arrangement was welldistributed, and temperature difference was neglected; (b) insulating paint distribution of copper wire was welldistributed; (c) copper wire insulation and fill paint temperature change was linear; (d) winding impregnating paint was completely filled.Therefore, the stator slot insulation is expressed as Eq. ( 13), which was constituted of copper wire paint, varnish and slot insulation: where ε is the equivalent thermal conductivity coefficient of slot insulation, δ i is the equivalent thickness of various insulating materials, α i is the thermal conductivity coefficient of various insulating materials.

Convection heat transfer coefficient of air gap
The in-wheel motor of micro-electric vehicle adopted a closed construction, and there was no interchange of air between inner and outer environment.The air gap between stator and rotor was dealt with, and the revolving rotor made the air of air gap flow when the inwheel motor was running.On the one hand, the flowing air is influenced by tangential motion of rotor.On the other hand, the flowing air is blocked by the inner circle surface of stator.Therefore, the complex coupling field was formed by fluid field and temperature field in the air gap.In order to simplify the calculation, the thermal conductivity coefficient of static fluid was used to describe the heat transfer ability of the flowing air in the air gap, which brought to an effective heat exchange coefficient [12,13].If stator inner surface and rotor outer surface are smooth, the Renault coefficient in the air gap is expressed as Eq. ( 14), and the Critical Renault coefficient is expressed as Eq. ( 15): where n is the revolving speed of motor, δ air is the length of air gap, d out is the outer diameter of rotor, τ is the air motion viscosity coefficient, and d in is the inner diameter of stator.When Re < Re cr , the flowing air in the air gap is laminar flow, and convection heat transfer coefficient of air gap is expressed as Eq. ( 16).When Re > Re cr , the flowing air in the air gap is turbulent flow, and convection heat transfer coefficient of air gap is expressed as Eq. ( 17): Here, the convection heat transfer coefficient of air gap ( eff α ) is 33 W/(m 2 ℃), the thermal conductivity of air ( air α ) is 0.026 W/(m℃) and the anger XieErTe number ( Nu ) is tow.

Establishment of finite element model
The united simulation model of in-wheel motor was established by Ansoft Maxwell and Ansys Workbench software, and thermal loss and temperature field were simulated and analyzed.According to the control requirements of micro-electric vehicle system, the main technical parameters of motor were optimized and designed by Genetic Algorithm.The structure of poles was tile shape and radial magnetization.The main technical parameters of motor were shown in Table 1.Through the symmetric simplified treatment, the three-dimensional finite element model was established, as shown in Fig. 1.To ensure the correctness of the magnetic circuit calculation and magnetic field analysis, the finite element meshes were manually divided by the inside selection and surface approximation method, and the solution magnetic field was divided into 197588 units.The overall mesh subdivision model is shown in Fig. 2. Fig. 2 shows that the subdivision of finite element grid is relatively uniform distribution on the whole，and the grid subdivision has a higher density in the field where the magnetic field is rela-tively strong and changes greatly，such as the air gap part.Therefore, the subdivision of finite element grid can reach the desired effect and satisfy the accuracy of finite element calculation.
Fig. 1 Three-dimensional finite element model Fig. 2 Overall mesh subdivision model

Analysis of thermal loss
To reduce core loss, the stator and rotor adopted the core material parameters of Toyota Prius, which was non-linear B-H characteristic.To accurately obtain the core loss, loss curves of materials with different frequency was calculated, as shown in Fig. 3. Fig. 3 shows that iron loss factor can be easily worked out.The hysteresis loss factor k h is 184.234, the eddy loss factor k e is 0.27, and the additional loss factor k ex is 0.386.

Fig. 3 Loss curves of materials with different frequency
The winding copper loss, core loss and eddy current loss of permanent magnet were calculated by Ansoft Maxwell software.The heat source loss curve is shown in Fig. 4. Fig. 4 shows that the fluctuation of winding copper loss is very large before the in-wheel motor operates at 26 ms, and the maximum amplitude is 110 W, but the fluctuations of the core loss and the eddy current loss of permanent magnet are low.As the in-wheel motor can gener-ated a big fluctuant starting current.After the in-wheel motor operates at 26 ms, it is in stable operation, and the heat source loss tends to cyclical stability.The core loss of stator and rotor is greater than the eddy current loss of permanent magnet, but less than the winding copper loss, and the winding copper loss accounts for the main loss.In addition, Fig. 4 also shows that the winding copper loss is stability in 33 W, the core loss is stability in 3.4 W, and the eddy current loss of permanent magnet is stability in 1.6 W. According to the previous formula, it can be calculated that the mechanical loss is about 1.8 W and stray loss is about 1.5 W. The efficiency of in-wheel motor is up to 88.2%, which can fully meet the normal operation requirements of the micro-electric vehicle.Firstly, the thermal loss was calculated by Ansoft Maxwell software, which was changed into the heat flux by Ansys Workbench software.Secondly, the heat flux was loaded in the temperature model to calculate the rise of temperature.Finally, the temperature distribution of the heat source loss and each part of in-wheel motor was obtained.When the in-wheel motor is running, all heats come from its loss.The stator and rotor, equivalent winding and permanent magnet are not only the heat source components, but also the heat transfer components, while other parts just are the heat transfer components.
When the in-wheel motor was running at the rated load, the transient state temperature field distribution was calculated, and transient temperature curve of various components was simulated, as shown in Fig. 5.  5 shows that, in the entire operation, the temperature of equivalent winding is always higher than stator temperature, and the temperature of stator is always higher than the temperature of rotor and permanent magnets.The major reasons are given: on the one hand, the copper loss of equivalent winding is located in the stator; on the other hand, stator iron core loss is more than rotor iron core loss and the eddy current loss of permanent magnet.Temperature of the rotor and the permanent magnet is basically the same change trend and value.The major reasons are that the rotor iron core loss and the eddy current loss of permanent magnet are relatively small and the cooling environment is bad.Before 1800 s the winding temperature approximately increases linearly, temperature rising trend slows between 1800 and 4000 s, after 4000 s the temperature tends towards basic stability; Before 1650 s the stator temperature approximately increases linearly, temperature rising trend slows between 1650 and 3800 s, after 3800 s the temperature tends towards basic stability; Before 1500 s the temperature of the rotor and the permanent magnet approximately increases linearly, temperature rising trend slows between 1500 and 3600 s, after 3600 s the temperature tends towards basic stability.

Analysis of temperature influence factors
When the in-wheel motor was running stably for 1.5 hours at the rated speed of 1300 r min -1 , the influence of load on temperature rise is simulated, as shown in Fig. 6.Fig. 6 shows when the load of in-wheel motor is less than the rated load, the temperature rise rate of the winding and stator increases slowly.However, the load of in-wheel motor is more than 85 kg, the temperature rise rate of the winding and stator increases rapidly.At this time, the output power of the in-wheel motor increases, power loss also increases, which brings to lower efficiency.When the in-wheel motor was running stably for 1.5 hours at the rated load, the influence of revolving speed on temperature rise is simulated, as shown in Fig. 7. Fig. 7 shows the revolving speed of in-wheel motor is less than 750 r/min, and the temperature rise rate of the winding and stator increases slowly.The revolving speed of in-wheel motor is at 750 and 1500 r/min, and the temperature rise rate of the winding and stator increases linearly.However, the revolving speed of in-wheel motor is more than 1500 r/min, the temperature rise rate of the winding and stator increases rapidly.At this time, the output power of the in-wheel motor increases, power loss also increases, which brings to lower efficiency.Fig. 6 and Fig. 7 show that the load of in-wheel motor is more than 110 kg or the revolving speed of inwheel motor is more than 1800 r/min, and the temperature of winding and stator can be more than 100°C, which can lead to the overheating of in-wheel motor, the burning out of winding insulation layer, the failure of Hall sensor and the demagnetization of permanent magnet, and even seriously the in-wheel motor to stall.Therefore, the load and revolving speed of the in-wheel motor have an effect on the temperature rise of in-wheel motor, while over a certain value, rise of temperature increases rapidly, which can seriously reduce the efficiency and life of in-wheel motor.

Analysis of test results
Temperature rise test of in-wheel motors was executed at rated load and ambient temperature was 22°C.The surface temperatures of in-wheel motor shell and stator were measured by the thermal infrared imager，and the temperatures of stator winding and air gap were measured by the thermistor embedded [13][14][15].The comparison of the calculated results with the test average data is shown in Table 2. Table 2 shows that the calculated results have a slightly smaller than the test average data and the maximal relative error which can satisfy the project requirements is within 4.6%.There are two reasons for the relative error: one is system measuring error and the other is related to boundary conditions, assumptions of finite element model, the mesh quality and equivalent treatment processes method.In addition, Table also shows that the range of temperature is small among the parts, which indicates the inwheel motor has a good heat transfer effect.The maximum measure temperature is about 77.8°C which is inferior to the maximum working temperature 100°C of Nd-Fe-B magnet and the maximum allowable temperature 120°C of B-class insulation material as well as the limit temperature 80°C of winding.Therefore, loss of excitation of permanent magnet and insulation aging are not easy to reduce, which can meet the requirements of the normal work of the in-wheel motor at the rated conditions.In a word, the simulation and test results are basically consistent.

Conclusions
The results of study in summary can be classified as follows: 1.The total loss of the in-wheel motor is analyzed and calculated, in order to establish the mathematical model of heat loss.Winding copper loss, stator and rotor iron core loss and eddy current loss of permanent magnet are calculated, which are coupled to the temperature field as the heat sources.
2. Mathematical model of temperature field is established.The equivalent model of stator winding is adequately handled, convection heat transfer coefficients are calculated, and the heat distribution of in-wheel motor is analyzed.
3. The united simulation model of in-wheel motor was established by Ansoft Maxwell and Ansys Workbench software, and thermal loss and temperature field were simulated and analyzed.The fluctuation of winding copper loss is very large before the in-wheel motor operates at 26 ms, and the maximum amplitude is 110 W, and the stator winding loss is the largest loss.However the fluctuations of the core loss and the eddy current loss of permanent magnet are low.
4. Transient temperature curve of various components was simulated.The temperature of equivalent winding is always higher than stator temperature, and the temperature of stator is always higher than the temperature of rotor and permanent magnets.Temperature of the rotor and the permanent magnet is basically the same change trend and value.
5. The load and revolving speed of the in-wheel motor have an effect on the temperature rise of in-wheel motor, while over a certain value, rise of temperature increases rapidly, which can seriously reduce the efficiency and life of in-wheel motor.
6.The simulation and test results are a certain relative error, but there are basically consistent.

Q. P. Chen STUDY ON TEMPERATURE INFLUENCE FACTORS AND ELECTROMAGNETIC HEAT COUPLING ON IN-WHEEL MOTOR FOR MICRO-ELECTRIC VEHICLE S u m m a r y
In order to meet the development requirements of micro-electric vehicle in efficiency, reliability and safety, thermal loss mechanism and temperature influence factors of in-wheel motor are study.The finite element method is used to simulate and analyze internal complex thermal loss and temperature field.The mathematical model of thermal loss and temperature field was established.The united simulation model of in-wheel motor was established by Ansoft Maxwell and Ansys Workbench software, transient state temperature were simulated and analyzed.The influences of load and revolving speed on temperature rise were also studied.The conformity of the test results with the final simulation results indicates that this method can be used to provide a theoretical basis to make further optimal design of the new driving in-wheel motor.

Fig. 5
Fig.5Transient temperature curve of various components Fig.5shows that, in the entire operation, the temperature of equivalent winding is always higher than stator temperature, and the temperature of stator is always higher than the temperature of rotor and permanent magnets.The major reasons are given: on the one hand, the copper loss

Fig. 6 Fig. 7
Fig. 6 Influence of load on temperature rise

Table 1
Main technical parameters

Table 2
Comparison of the calculated results with the test average data