CN114707360A - Novel stator permanent magnet motor iron loss calculation method - Google Patents

Novel stator permanent magnet motor iron loss calculation method Download PDF

Info

Publication number
CN114707360A
CN114707360A CN202210510042.4A CN202210510042A CN114707360A CN 114707360 A CN114707360 A CN 114707360A CN 202210510042 A CN202210510042 A CN 202210510042A CN 114707360 A CN114707360 A CN 114707360A
Authority
CN
China
Prior art keywords
iron loss
harmonic
calculation method
permanent magnet
elliptical
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210510042.4A
Other languages
Chinese (zh)
Inventor
朱晓锋
金子辰
丁树业
李光志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN202210510042.4A priority Critical patent/CN114707360A/en
Publication of CN114707360A publication Critical patent/CN114707360A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • 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/60—Other road transportation technologies with climate change mitigation effect
    • Y02T10/64—Electric machine technologies in electromobility

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Iron Core Of Rotating Electric Machines (AREA)

Abstract

The invention provides a novel method for calculating iron loss of a stator permanent magnet motor, which comprises the following steps: firstly, selecting representative points for analyzing the magnetic density change rule in the iron core; then obtaining an elliptical magnetization trajectory diagram through simulation; decomposing an elliptical harmonic magnetic field; solving magnetic flux density amplitude and direct current magnetic bias amplitude on the major axis and the minor axis of each subellipse harmonic; the iron loss is finally calculated by the proposed formula. The novel method for calculating the iron loss of the stator permanent magnet motor, provided by the invention, takes the influence of various factors into consideration, and comprises the following steps: the method has the advantages of high accuracy, simple and convenient calculation process, capability of calculating iron loss of different motors under different working conditions and the like, and provides guidance for design and optimization of the motors.

Description

Novel stator permanent magnet motor iron loss calculation method
Technical Field
The invention relates to the field of motor iron loss calculation, in particular to a novel stator permanent magnet motor iron loss calculation method considering the influences of an elliptic rotation magnetization phenomenon, a direct-current magnetic bias phenomenon, a small magnetic hysteresis loop phenomenon and PWM harmonic current.
Background
With global warming, shortage of fossil fuels and environmental pollution, how to utilize energy more efficiently becomes a problem that people must face in the development process. The permanent magnet synchronous motor has the advantages of small volume, high power density, high efficiency, high output torque, reliable operation and the like, is widely applied in many fields, and particularly plays an important role in the development of new energy automobiles in the field of automobile industry.
The traditional permanent magnet motor mostly adopts the surface mounting type, plug-in type or embedded type and other processes to fix the permanent magnet on the rotor, so that the problems of difficult heat dissipation, irreversible demagnetization of the permanent magnet caused by overhigh temperature rise and the like can be generated in the running process of the motor, and the output of the motor is limited. To solve this problem, in the 50 s of the last century, researchers have developed studies on stator permanent magnet motors. Different from a rotor permanent magnet motor, the stator permanent magnet motor arranges the permanent magnets and the armature windings in the stator, the simple rotor structure not only improves the running reliability of the motor, but also can play a role of a fan when the rotor rotates, so that the heat radiation performance of the motor is enhanced.
Stator permanent magnet motors are receiving more and more attention due to their substantial advantages of solid structure, high power (torque) density, good heat dissipation, and high fault tolerance. It is considered to be a combination of switched reluctance and conventional permanent magnet brushless motors and is becoming more and more widely used in many electrical industry applications. However, few researchers have studied iron loss calculation models for stator permanent magnet motors. Since the structure of the motor is completely different from that of a conventional permanent magnet motor, and the loss distribution is completely different, it is very important to research the loss of the motor.
Through retrieval, the Chinese patent with the application number of CN202110629757.7 discloses a method and a model for rapidly calculating the iron loss of a high-speed motor based on the rotating speed. The motor is a high-speed motor, the method and the model are based on a classical iron loss calculation model, a variable coefficient iron loss calculation model considering factors such as higher harmonics, rotating magnetization, skin effect and small hysteresis loop is obtained, the model can fully reflect the change of iron loss along with the flux density waveform distortion rate, the flux density amplitude and the motor frequency, and the relationship between the iron loss and the rotating speed is obtained by carrying out nonlinear curve fitting on the motor iron loss calculation results at different rotating speeds, so that an iron loss calculation model related to the rotating speed is obtained, and the model can quickly obtain the corresponding iron loss at different motor rotating speeds and has higher calculation accuracy. The method for establishing the iron loss rapid calculation model has wide application range, is not only suitable for calculating the iron loss of the high-speed permanent magnet motor, but also can be used for other types of high-speed motors.
The invention provides a novel method for calculating iron loss of a stator permanent magnet motor. The method comprises the following steps: firstly, selecting a representative point for analyzing the change rule of the magnetic density in the iron core; then obtaining an elliptical magnetization trajectory diagram through simulation; decomposing an elliptical harmonic magnetic field; solving the magnetic flux density amplitude and the direct current magnetic bias amplitude on the major axis and the minor axis of each subellipse harmonic; the iron loss is finally calculated by the proposed formula. The novel method for calculating the iron loss of the stator permanent magnet motor, provided by the invention, takes the influence of various factors into consideration, and comprises the following steps: the method has the advantages of high accuracy, simple and convenient calculation process, capability of calculating iron loss of different motors under different working conditions and the like, and provides guidance for design and optimization of the motors.
Compared with the above patent documents, the solution proposed by the present invention is different in that:
1. the influence factors considered by the invention include elliptic rotation magnetization phenomenon, direct current magnetic bias phenomenon, small magnetic hysteresis loop phenomenon and PWM harmonic current.
2. The research object of the invention is a stator permanent magnet motor;
the research method of the invention is as follows: firstly, selecting a representative point for analyzing the change rule of the magnetic density in the iron core; then obtaining an elliptical magnetization trajectory diagram through simulation; decomposing an elliptical harmonic magnetic field; solving the magnetic flux density amplitude and the direct current magnetic bias amplitude on the major axis and the minor axis of each subellipse harmonic; the iron loss is finally calculated by the proposed formula.
Disclosure of Invention
In order to overcome the defects of few consideration factors and inaccurate result of the traditional motor iron loss calculation method, the invention provides a novel iron loss calculation method which particularly aims at a stator permanent magnet motor and considers the influence of elliptic rotation magnetization phenomenon, direct current magnetic bias phenomenon, small magnetic hysteresis loop phenomenon and PWM harmonic current
In order to achieve the purpose, the invention adopts the technical scheme that:
a novel stator permanent magnet motor iron loss calculation method is characterized by comprising the following steps: comprises the following steps:
step 1: selecting a representative point, and selecting a representative point,
selecting representative points on a stator iron core and a rotor iron core of the motor;
step 2: obtaining an elliptical magnetization trace diagram,
solving a oscillogram of each representative point, which changes along with the tangential and radial flux densities, by an Ansys Maxwell software simulation, and obtaining an elliptical magnetization trajectory graph of the change relation between the radial flux density and the tangential flux density of each representative point;
and step 3: the decomposition of the elliptical harmonic magnetic field is carried out,
decomposing the irregular elliptical magnetic field of each representative point into a series of elliptical harmonic flux density vectors by means of a harmonic analysis principle;
and 4, step 4: find BmAnd the sum of the delta B and the delta B,
2 mutually orthogonal alternating magnetizations are adopted for each subelliptic harmonic magnetic field to be equivalent, namely, the elliptic magnetic flux density is decomposed into the directions of a long axis and a short axis to calculate the iron loss, and then the magnetic flux density amplitude B on the long axis and the short axis of each subelliptic harmonic is sequentially solvedmAnd a DC bias amplitude Δ B;
and 5: the iron loss is calculated and calculated,
b of each sub-elliptical harmonic of each representative pointmAnd substituting the calculated values and the delta B into an iron loss calculation formula in sequence to obtain an iron loss density value of each representative point, taking the average value of each representative point as the iron loss density of the motor, and finally multiplying the iron loss density by the iron loss volume to obtain the iron loss, wherein the iron loss calculation formula is as follows:
Figure BDA0003639051690000031
wherein, PFeThe total iron loss density is obtained; phHysteresis loss density; peIs the eddy current loss density; k is a radical ofhIs a hysteresis loss coefficient; k is a radical ofeIs the eddy current loss coefficient; f is fundamental frequency, vf is v-order elliptic harmonic frequency; b ismMvAnd BmNvMagnetic flux density amplitudes of the v-th order elliptic harmonic waves in the major axis direction and the minor axis direction respectively; delta BMvAnd Δ BNvThe magnitudes of the direct current magnetic biases of the v-th order elliptical harmonic waves in the major axis direction and the minor axis direction are respectively.
As a preferred technical scheme of the invention: seven representative points are selected and respectively represent a stator tooth tip, a stator tooth middle part, a stator tooth connecting yoke part, a stator yoke, a rotor tooth tip, a rotor tooth connecting yoke part and a rotor yoke.
As a preferred technical scheme of the invention: the stator tooth point has a direct current magnetic bias induction phenomenon, the rotor tooth point has a small magnetic hysteresis loop phenomenon, and an elliptical rotation magnetization phenomenon exists in the iron core.
As a preferred technical scheme of the invention: the iron loss calculation method takes into account the influence of the elliptical rotation magnetization phenomenon on the iron loss.
As a preferred technical scheme of the invention: the iron loss calculation method takes the influence of the direct current magnetic biasing phenomenon on the iron loss into account.
As a preferred technical scheme of the invention: when the DC bias only affects the hysteresis loss, the hysteresis loss density PhAnd dc bias induction Δ B is expressed as:
Ph(ΔB)＝Ph0ε(ΔB)
ε(ΔB)＝1+kdcΔBβ+klΔB2
in the formula, Ph(Δ B) and Ph0Hysteresis losses, k, with and without influence of dc-biased magnetization induction, respectivelydc、klAnd β is a coefficient of no physical significance, equal to 0.4575, 0.376 and 5.5975 by nonlinear fitting, respectively, as can be seen by the above equationThe effect on the iron loss calculation can be accounted for by ε (Δ B).
As a preferred technical scheme of the invention: the iron loss calculation method takes into account the influence of the small hysteresis loop phenomenon on the iron loss.
As a preferred technical scheme of the invention: the iron loss calculation method takes into account the influence of PWM harmonic current on iron loss.
The novel method for calculating the iron loss of the stator permanent magnet motor, provided by the invention, takes the influence of various factors into consideration, and comprises the following steps: the elliptic rotation magnetization phenomenon, the direct current magnetic biasing phenomenon, the small magnetic hysteresis loop phenomenon and the PWM harmonic current improve the accuracy of the calculation of the iron loss model and provide guidance for the design and optimization of the motor.
Compared with the prior art, the invention has the beneficial effects that:
(1) compared with the traditional iron loss calculation method, the iron loss calculation method provided by the invention considers the influence of various factors, and the calculation accuracy is improved.
(2) The iron loss calculation method provided by the invention is easy to understand, convenient to operate, simple and convenient in calculation process and less in calculation time.
(3) The iron loss calculation method provided by the invention can be used for calculating the iron loss of all motors and not only can be used for stator permanent magnet motors.
(4) The iron loss calculation method provided by the invention can calculate the iron loss values of the motor under no-load and load working conditions.
Drawings
FIG. 1 is a schematic diagram of the logic structure of the present invention;
FIG. 2 is a graph of typical spot locations and flux density variation waveforms for an example of the present invention;
FIG. 3 is a typical trace of spot flux density variation for an embodiment of the present invention;
FIG. 4 is an elliptical harmonic exploded view of a rotor tooth tip of an embodiment of the present invention;
FIG. 5 is a schematic illustration of an elliptical rotating magnetization;
FIG. 6 is a diagram illustrating DC magnetic bias and small hysteresis loop;
FIG. 7 is a schematic diagram showing the effect of PWM harmonic current on the magnetic flux density of an iron core;
FIG. 8 is a diagram illustrating sub-elliptical magnetic fields and their harmonic decomposition.
List of reference numerals: 1. a stator tooth tip; 2. in the stator teeth; 3. a stator tooth yoke portion; 4. a stator yoke; 5. a rotor tooth tip; 6. a rotor tooth yoke portion; 7. a rotor yoke.
Detailed Description
The invention is described in further detail below with reference to the following detailed description and accompanying drawings:
the invention provides a novel method for calculating iron loss of a stator permanent magnet motor, which comprises the following steps:
step 1: selecting a representative point, and selecting a representative point,
referring to fig. 2, taking a flux-reversing motor as a representative stator permanent magnet motor as an example, some representative points are selected on a stator core and a rotor core of the motor to analyze a flux density variation law in the core. Points 1-7 in fig. 2 represent the stator tooth tip 1, the stator tooth middle 2, the stator tooth yoke 3, the stator yoke 4, the rotor tooth tip 5, the rotor tooth yoke 6 and the rotor yoke 7, respectively.
Step 2: obtaining an elliptical magnetization trace diagram,
the oscillogram of each representative point, namely the tangential and radial flux densities, which change along with the position of the rotor is solved through Ansys Maxwell software simulation, as shown in FIG. 2, and an elliptical magnetization trajectory diagram of the change relation between the radial flux density and the tangential flux density of each representative point is obtained, as shown in FIG. 3.
From fig. 2 and 3, it can be concluded that:
(1) the direct current magnetic bias induction phenomenon obviously exists on the stator tooth tip 1;
(2) the rotor tooth tip 5 has some minor hysteresis loops;
(3) the elliptical rotation magnetization phenomenon widely exists in the iron core.
And step 3: the decomposition of the elliptical harmonic magnetic field is carried out,
for an irregular elliptical magnetic field, it is decomposed into a series of elliptical harmonic flux density vectors by the principle of harmonic analysis, as shown in fig. 4.
And 4, step 4: find BmAnd the sum of the delta B and the delta B,
2 mutually orthogonal alternating magnetizations are adopted for each subelliptic harmonic magnetic field for equivalence, namely, the elliptic magnetic density is decomposed into the directions of a long axis and a short axis for iron loss calculation, and the magnetic density amplitude B on the long axis and the short axis of each subelliptic harmonic is sequentially solvedmAnd a dc bias amplitude Δ B.
And 5: the iron loss is calculated and calculated,
b of each elliptical harmonic of each pointmAnd substituting the delta B into an iron loss calculation formula in sequence to obtain an iron loss density value of each point, taking the average value of each point as the iron loss density of the motor, and finally multiplying the iron loss density by the iron loss volume to obtain the iron loss, wherein the iron loss calculation formula is as follows:
Figure BDA0003639051690000051
wherein, PFeThe total iron loss density is obtained; phHysteresis loss density; p iseIs the eddy current loss density; k is a radical of formulahIs a hysteresis loss coefficient; k is a radical ofeIs the eddy current loss coefficient; f is fundamental frequency, vf is v-order elliptic harmonic frequency; bmMvAnd BmNvThe magnetic flux density amplitudes of the v-order elliptic harmonic waves in the major axis direction and the minor axis direction respectively; delta BMvAnd Δ BNvThe magnitudes of the direct current magnetic biases of the v-th order elliptical harmonic waves in the major axis direction and the minor axis direction are respectively.
In this embodiment: the iron loss calculation method takes into account the influence of the elliptical rotation magnetization phenomenon on the iron loss.
The ferromagnetic material has 2 magnetization modes of alternating magnetization and rotating magnetization, and the rotating magnetization is divided into circular rotating magnetization and elliptical rotating magnetization. The direction of the alternating magnetization magnetic field is unchanged, and the magnitude of the alternating magnetization magnetic field changes along with time; the size of the magnetic field of the circular rotating magnetization is unchanged, and the magnetization direction changes along with time; for elliptical magnetization, both magnitude and magnetization direction change with time. A schematic diagram of the elliptical rotational magnetization is shown in fig. 5. it can be seen from fig. 5 that the major and minor axes of the ellipse do not coincide with the radial and tangential directions. In order to ensure the accuracy of the calculation of the iron loss, the elliptical magnetization track is decomposed along the major axis and the minor axis, and the iron is calculated respectivelyThe iron consumption is more reasonable. The major axis direction of the magnetization ellipse is defined as a straight line direction between a point of magnetic flux density having the largest absolute value and the origin on the rotating magnetization trajectory, and the minor axis direction is a direction perpendicular to the major axis. B ismajAnd BminThe magnetic flux density amplitudes in the major and minor axis directions, respectively.
In the calculation of iron loss in the prior art, k.yamazaki proposes a hysteresis loop superposition calculation method, and the expression is as follows:
Figure BDA0003639051690000061
Figure BDA0003639051690000062
wherein B isrmAnd BtmThe magnitudes of the flux densities in the radial and tangential directions of each hysteresis loop, respectively. N is a radical ofrAnd NtThe total number of hysteresis loops for the radial and tangential components, respectively. In the Yamazaki model, the flux density trajectory is divided radially and tangentially for calculation, as shown in fig. 4, which reduces the accuracy of the core loss calculation. Therefore, the invention calculates by decomposing along the long axis and the short axis, and the expression is converted into
Figure BDA0003639051690000063
Figure BDA0003639051690000064
The solution of the iron loss is carried out by decomposing the elliptical rotation magnetization track into the directions of the long axis and the short axis, and the solution is more accurate than the traditional solution of decomposing the magnetic flux density into the radial direction and the tangential direction.
In this embodiment: the iron loss calculation method takes the influence of the direct current magnetic biasing phenomenon on the iron loss into account.
The DC magnetic bias phenomenon is shown in FIG. 6 and exists due to the alternating magnetic densityIn the dc magnetic biasing, the corresponding hysteresis loop will be distorted, and the corresponding hysteresis loss will be increased. It is generally believed that the DC bias affects only the hysteresis loss, and the hysteresis loss density PhAnd dc bias induction Δ B is expressed as:
Ph(ΔB)＝Ph0ε(ΔB)
ε(ΔB)＝1+kdcΔBβ+klΔB2
in the formula, Ph(Δ B) and Ph0With and without hysteresis losses due to dc bias magnetization induction, respectively. k is a radical ofdc、klAnd β are coefficients of no physical significance, equal to 0.4575, 0.376, and 5.5975 by non-linear fit, respectively. From the above equation it can be seen that the effect of dc bias magnetization on the iron loss calculation can be accounted for by epsilon (deltab).
In this embodiment: the iron loss calculation method takes into account the influence of the small hysteresis loop phenomenon on the iron loss.
Because the stator permanent magnet motor adopts a double salient pole structure, the cogging effect is obvious, and the magnetic flux density waveform in an iron core (particularly a rotor tooth) is seriously distorted, so that a small magnetic hysteresis loop appears, as shown in fig. 6.
In this embodiment: the iron loss calculation method takes into account the influence of PWM harmonic current on iron loss.
The high-speed permanent magnet brushless motor powered by the inverter contains high-frequency current harmonic waves in the stator current due to PWM modulation, and the high-frequency harmonic wave current can increase the core loss.
As shown in fig. 7, the distortion current injected by the inverter may create more and smaller magnetic hysteresis loops, and thus additional hysteresis losses, than the no-load condition and the sinusoidal current input. The magnetic flux density harmonics caused by the PWM carrier also generate additional eddy current losses.
In order to rapidly, conveniently and accurately obtain the number and corresponding amplitude of each small hysteresis loop, an irregular elliptical magnetic field can be decomposed into a series of elliptical harmonic flux density vectors by means of a harmonic analysis principle, as shown in fig. 8. 2 mutually orthogonal alternating magnetizations are adopted for each subelliptical harmonic magnetic field along the directions of the long axis and the short axis to achieve equivalence, and in consideration of the direct current magnetic bias influence epsilon (delta B), the iron loss calculation formula is as follows:
Figure BDA0003639051690000071
in summary, the novel method for calculating the iron loss of the stator permanent magnet motor provided by the invention takes the influence of various factors into consideration, and comprises the following steps: the method has the advantages that the elliptic rotation magnetization phenomenon, the direct current magnetic biasing phenomenon, the small magnetic hysteresis loop phenomenon and the PWM harmonic current are adopted, the calculation accuracy of the iron loss model is improved, and guidance is provided for the design and optimization of the motor.
Compared with the prior art, the invention has the following advantages:
(1) compared with the traditional iron loss calculation method, the iron loss calculation method provided by the invention considers the influence of various factors, and the calculation accuracy is improved.
(2) The iron loss calculation method provided by the invention is easy to understand, convenient to operate, simple and convenient in calculation process and less in calculation time.
(3) The iron loss calculation method provided by the invention can be used for calculating the iron loss of all motors and not only can be used for stator permanent magnet motors.
(4) The iron loss calculation method provided by the invention can calculate the iron loss values of the motor under no-load and load working conditions.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art should be considered to be within the technical scope of the present invention, and the technical solutions and the inventive concepts thereof according to the present invention should be equivalent or changed within the scope of the present invention.

Claims (8)

1. A novel stator permanent magnet motor iron loss calculation method is characterized by comprising the following steps: comprises the following steps:
step 1: selecting a representative point, and selecting a representative point,
selecting representative points on a stator iron core and a rotor iron core of the motor;
step 2: obtaining an elliptical magnetization trace diagram,
solving a oscillogram of each representative point, which changes along with the tangential and radial flux densities, by an Ansys Maxwell software simulation, and obtaining an elliptical magnetization trajectory graph of the change relation between the radial flux density and the tangential flux density of each representative point;
and step 3: the decomposition of the elliptical harmonic magnetic field is carried out,
decomposing the irregular elliptical magnetic field of each representative point into a series of elliptical harmonic flux density vectors by means of a harmonic analysis principle;
and 4, step 4: find BmAnd the sum of the delta B and the delta B,
2 mutually orthogonal alternating magnetizations are adopted for each subelliptic harmonic magnetic field to be equivalent, namely, the elliptic magnetic flux density is decomposed into the directions of a long axis and a short axis to calculate the iron loss, and then the magnetic flux density amplitude B on the long axis and the short axis of each subelliptic harmonic is sequentially solvedmAnd a DC bias amplitude Δ B;
and 5: the iron loss is calculated and calculated,
b of each sub-elliptical harmonic of each representative pointmAnd substituting the calculated values and the delta B into an iron loss calculation formula in sequence to obtain an iron loss density value of each representative point, taking the average value of each representative point as the iron loss density of the motor, and finally multiplying the iron loss density by the iron loss volume to obtain the iron loss, wherein the iron loss calculation formula is as follows:
Figure FDA0003639051680000011
wherein, PFeThe total iron loss density is obtained; phHysteresis loss density; peIs the eddy current loss density; k is a radical ofhIs a hysteresis loss coefficient; k is a radical ofeIs the eddy current loss coefficient; f is fundamental frequency, vf is v-order elliptic harmonic frequency; b ismMvAnd BmNvMagnetic flux density amplitudes of the v-th order elliptic harmonic waves in the major axis direction and the minor axis direction respectively; delta BMvAnd Δ BNvThe magnitudes of the direct current magnetic biases of the v-th order elliptical harmonic waves in the major axis direction and the minor axis direction are respectively.
2. The novel stator permanent magnet motor iron loss calculation method according to claim 1, characterized in that: seven representative points are selected and respectively represent a stator tooth tip (1), a stator tooth middle (2), a stator tooth connecting yoke portion (3), a stator yoke (4), a rotor tooth tip (5), a rotor tooth connecting yoke portion (6) and a rotor yoke (7).
3. The novel stator permanent magnet motor iron loss calculation method according to claim 2, characterized in that: there is direct current magnetic biasing induction phenomenon on stator prong (1), rotor prong (5) have little magnetic hysteresis loop phenomenon, have oval rotatory magnetization phenomenon in the iron core.
4. The novel stator permanent magnet motor iron loss calculation method according to claim 1 or 3, characterized in that: the iron loss calculation method takes into account the influence of the elliptical rotation magnetization phenomenon on the iron loss.
5. The novel stator permanent magnet motor iron loss calculation method according to claim 1 or 3, characterized in that: the iron loss calculation method takes the influence of the direct current magnetic biasing phenomenon on the iron loss into account.
6. The novel stator permanent magnet motor iron loss calculation method according to claim 5, characterized in that: when the DC bias only affects the hysteresis loss, the hysteresis loss density PhAnd dc bias induction Δ B is expressed as:
Ph(ΔB)＝Ph0ε(ΔB)
ε(ΔB)＝1+kdcΔBβ+klΔB2
in the formula, Ph(Δ B) and Ph0Hysteresis losses, k, with and without influence of DC-biased magnetization induction, respectivelydc、klAnd beta is a coefficient without physical meaning, and is respectively equal to 0.4575, 0.376 and 5.5975 through nonlinear fitting, and the influence of direct current bias magnetization on iron loss calculation can be seen through the formulaTo account for by epsilon (deltab).
7. The novel stator permanent magnet motor iron loss calculation method according to claim 1 or 3, characterized in that: the iron loss calculation method takes into account the influence of the small hysteresis loop phenomenon on the iron loss.
8. The novel stator permanent magnet motor iron loss calculation method according to claim 1 or 3, characterized in that: the iron loss calculation method takes into account the influence of PWM harmonic current on iron loss.
CN202210510042.4A 2022-05-11 2022-05-11 Novel stator permanent magnet motor iron loss calculation method Pending CN114707360A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210510042.4A CN114707360A (en) 2022-05-11 2022-05-11 Novel stator permanent magnet motor iron loss calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210510042.4A CN114707360A (en) 2022-05-11 2022-05-11 Novel stator permanent magnet motor iron loss calculation method

Publications (1)

Publication Number Publication Date
CN114707360A true CN114707360A (en) 2022-07-05

Family

ID=82176470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210510042.4A Pending CN114707360A (en) 2022-05-11 2022-05-11 Novel stator permanent magnet motor iron loss calculation method

Country Status (1)

Country Link
CN (1) CN114707360A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116956666A (en) * 2023-07-11 2023-10-27 青岛大学 Axial asymmetric motor iron loss calculation method considering three-dimensional magnetic density distribution
WO2023234013A1 (en) * 2022-05-30 2023-12-07 Ｊｆｅスチール株式会社 Data set generation method, electromagnetic field analysis method, and computer program

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023234013A1 (en) * 2022-05-30 2023-12-07 Ｊｆｅスチール株式会社 Data set generation method, electromagnetic field analysis method, and computer program
CN116956666A (en) * 2023-07-11 2023-10-27 青岛大学 Axial asymmetric motor iron loss calculation method considering three-dimensional magnetic density distribution
CN116956666B (en) * 2023-07-11 2024-04-26 青岛大学 Axial asymmetric motor iron loss calculation method considering three-dimensional magnetic density distribution

Similar Documents

Publication Publication Date Title
Kang et al. A novel cogging torque reduction method for interior-type permanent-magnet motor
Li et al. Analysis of torque capability and quality in Vernier permanent-magnet machines
Zhou et al. Influence of magnet shape on the cogging torque of a surface-mounted permanent magnet motor
CN114707360A (en) Novel stator permanent magnet motor iron loss calculation method
Li et al. Flux-focusing permanent magnet machines with modular consequent-pole rotor
Zhang et al. Loss calculation and thermal analysis for high-speed permanent magnet synchronous machines
Jian et al. Analytical calculation of magnetic field in surface-inset permanent magnet motors
Qiu et al. Performance analysis and comparison of PMSM with concentrated winding and distributed winding
Yang et al. Analysis of the electromagnetic performance of homopolar inductor machine through nonlinear magnetic equivalent circuit and air-gap permeance function
Li et al. Electromagnetic performance analysis of consequent-pole PM machine with asymmetric magnetic pole
Duan et al. Influence of different Halbach arrays on performance of permanent magnet synchronous motors
Tavana et al. Analysis and design of magnetic pole shape in linear permanent-magnet machine
Woo et al. Electromagnetic characteristic analysis of high-speed motors with rare-earth and ferrite permanent magnets considering current harmonics
Xiang et al. Development of a novel radial-flux machine with enhanced torque profile employing quasi-cylindrical pm pattern
Zhao et al. Harmonic suppression and torque ripple reduction of a high-speed permanent magnet spindle motor
Hua et al. General power equation of switched reluctance machines and power density comparison
Geng et al. Influence of harmonic voltage coupling on torque ripple of permanent magnet synchronous motor
Qiu et al. Research on the influence of driving harmonic on electromagnetic field and temperature field of permanent magnet synchronous motor
CN114513106B (en) Permanent magnet synchronous linear motor magnetic pole structure design method based on PWM (pulse-Width modulation) concept
CN114741933A (en) Mathematical modulation model and magnetic field analysis method for double-permanent-magnet vernier motor
Renyuan Research and development of rare earth permanent magnet machines
Wang et al. Performance analysis and comparison of interior permanent magnet machine with FSCW and distributed windings for a hybrid electric vehicle
Qiu et al. Influence of different rotor magnetic circuit structure on the performance of permanent magnet synchronous motor
Chen et al. Optimisation analysis of cogging torque of in-wheel motor based on Ansoft parametric modelling for electric vehicles
Zhang et al. Accurate Calculation of Iron Loss Based on Frequency Division and Piecewise Variable Coefficient for Magnetic Field Modulation Permanent Magnet Machine

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination