CN116822095A - Magnetic circuit modeling method for double-stator single-rotor axial permanent magnet motor - Google Patents
Magnetic circuit modeling method for double-stator single-rotor axial permanent magnet motor Download PDFInfo
- Publication number
- CN116822095A CN116822095A CN202310860251.6A CN202310860251A CN116822095A CN 116822095 A CN116822095 A CN 116822095A CN 202310860251 A CN202310860251 A CN 202310860251A CN 116822095 A CN116822095 A CN 116822095A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- motor
- magnetic circuit
- stator
- permanent magnet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000035699 permeability Effects 0.000 claims abstract description 49
- 239000000463 material Substances 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000004804 winding Methods 0.000 claims abstract description 21
- 230000004907 flux Effects 0.000 claims description 40
- 230000004323 axial length Effects 0.000 claims description 4
- 230000005415 magnetization Effects 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 claims 2
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 15
- 238000010586 diagram Methods 0.000 description 10
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Permanent Field Magnets Of Synchronous Machinery (AREA)
Abstract
The invention discloses a double-stator single-rotor axial permanent magnet motor magnetic circuit modeling method, and belongs to the field of motor optimization design; a modeling method for a magnetic circuit of a double-stator single-rotor axial permanent magnet motor comprises the following steps: defining geometrical parameters and material parameters of the motor parts; dividing nodes of motor components, and expressing the magnetic resistance and magnetomotive force of each component according to the geometric parameters and the material parameters of the components so as to construct a magnetic network model; based on a general magnetic network model, considering magnetic leakage in a stator slot, magnetic leakage of an adjacent permanent magnet through a long pole shoe and magnetic leakage of the end part position of a stator winding; defining a connection mode among units of the motor; obtaining complete magnetic circuit models at different rotor positions according to different rotor positions corresponding to different moments; and carrying out nonlinear iterative calculation on the magnetic permeability of the stator core based on the complete magnetic circuit model to obtain a final convergence result of the magnetic circuit.
Description
Technical Field
The invention belongs to the field of motor optimization design, and particularly relates to a double-stator single-rotor axial permanent magnet motor magnetic circuit modeling method.
Background
Permanent magnet synchronous motors have evolved greatly over the past several decades. Radial flux permanent magnet machines are of many types and have found wide application in various fields. However, radial flux permanent magnet machines may not be usable in some cases due to their excessive axial length. Such as wind turbines, flywheel energy storage, automotive applications, and other applications with axial space limitations. The magnetic flux direction of the axial permanent magnet motor is along the axial direction, so that the magnetic energy density is high, and the energy exchange space is also large, and therefore, the torque density of the axial permanent magnet motor is improved compared with that of the radial permanent magnet motor.
The most accurate method for calculating the performance of the axial permanent magnet motor is three-dimensional finite element analysis, and the three-dimensional finite element method is based on an accurate geometric model of the motor, can calculate a space magnetic field with a complex structure, but requires too much calculation time. The calculation time can be greatly reduced by expanding the axial motor to a two-dimensional linear motor model, however, the method generally ignores the influence of the magnetic flux leakage at the end position of the stator winding, and as the load increases, the average electromagnetic torque error from the three-dimensional finite element analysis increases.
The magnetic circuit method can design a special parameter magnetic circuit for a certain type of motor, and can rapidly predict the performance of the motor in an early stage of motor design. In addition, the preliminary optimization design of the motor can be realized by a magnetic circuit method. The existing magnetic circuit method for the double-stator single-rotor axial permanent magnet motor generally does not consider magnetic leakage of the end part of a stator winding, and does not predict the performance of the motor when the load changes, and has certain defects. Therefore, the modeling of the magnetic circuit method for the double-stator single-rotor axial permanent magnet motor still has higher research value.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a magnetic circuit modeling method for a double-stator single-rotor axial permanent magnet motor.
The aim of the invention can be achieved by the following technical scheme:
a modeling method for a magnetic circuit of a double-stator single-rotor axial permanent magnet motor comprises the following steps:
defining geometrical parameters and material parameters of the motor parts;
dividing nodes of motor components, and expressing the magnetic resistance and magnetomotive force of each component according to the geometric parameters and the material parameters of the components so as to construct a magnetic network model;
based on a magnetic network model, considering magnetic leakage in a stator slot, magnetic leakage of an adjacent permanent magnet through a long pole shoe and magnetic leakage of the end part position of a stator winding; defining a connection mode among units of the motor; obtaining complete magnetic circuit models at different rotor positions according to different rotor positions corresponding to different moments;
and carrying out nonlinear iterative calculation on the magnetic permeability of the stator core based on the complete magnetic circuit model to obtain a final convergence result of the magnetic circuit.
Further, the geometric parameters of the motor component include: geometric parameters of stator, rotor and windings, geometric position, and air gap;
the material parameters of the motor component include: the stator and rotor materials have electromagnetic properties including permeability and magnetizing strength.
Further, the node division of the motor component includes dividing the stator, the rotor, and the air gap.
Further, the magnetic resistance is divided into rectangular columnar magnetic resistance, rectangular and semicircular combined columnar magnetic resistance and annular columnar magnetic resistance;
wherein, for rectangular columnar magnetic resistance, magnetic resistance R in the magnetic flux direction thereof 1 It can be expressed as that,
wherein L represents the magnetic path length of a unit of the motor in the magnetic flux direction, mu 0 Is vacuum permeability, mu r For a certain motor component relative permeability, S represents the area through which the flux path passes;
for a rectangular and semicircular combined columnar reluctance, the reluctance R in the direction of the magnetic flux thereof 2 It can be expressed as that,
wherein R is h The radius of the semicircle is indicated,L depth representing the thickness of the linear motor; x represents an integral variable of the radius of the arc;
L depth can be expressed as:
L depth ＝R o -R i (3)
wherein R is i ，R o The inner diameter and the outer diameter of the axial motor are respectively;
for annular columnar magnetoresistance, magnetoresistance R in the direction of magnetic flux thereof 3 It can be expressed as that,
wherein R is m ，R n Respectively the inner diameter and the outer diameter of the circular ring, and theta represents the radian swept by the circular ring.
Further, an axial magnetomotive force F of the three-phase armature magnetomotive force abc Can be expressed as:
wherein N is c Indicating the number of turns, i a 、i b 、i c Respectively representing instantaneous values of three-phase currents;
its axial magnetomotive force F of permanent magnetic magnetomotive force pm Can be expressed as:
F pm ＝H c H pm (6)
wherein H is c Indicating the magnetization intensity of the permanent magnet, H pm Representing the axial length of the permanent magnet.
Further, the step of nonlinear iterative calculation of the magnetic permeability of the stator core includes:
s41, initializing parameters of an equivalent magnetic circuit, giving initial magnetic permeability, and calculating initial magnetic resistance of a stator core; it should be noted that;
s42, performing magnetic circuit calculation to obtain magnetic flux phi of each position;
s43, leading out the magnetic flux density B of each position to a working area, and obtaining new magnetic permeability through a B-H curve of a common material;
s44, calculating the relative error between the current magnetic permeability and the magnetic permeability of the last iteration, if the relative error is larger than a given error delta, updating the magnetic permeability and continuing iteration until the relative error is smaller than delta, namely, the calculation result is converged, recording the current magnetic flux and the magnetic flux density at the moment, and finishing the calculation of the network iteration model of the motor.
Further, the iterative expression of the permeability is:
μ k ＝p 1 μ k +p 2 μ k-1 (7)
wherein mu is k Sum mu k-1 The current permeability and the permeability of the last iteration, p 2 And p 1 Is the current magnetic permeability coefficient and the magnetic permeability coefficient after the last iteration, p 1 +p 2 ＝1。
Further, the relative error convergence judgment expression is:
a magnetic circuit modeling system of a double-stator single-rotor axial permanent magnet motor comprises
Parameter definition module: defining geometrical parameters and material parameters of the motor parts;
the magnetic network model building module: dividing nodes of motor components, and expressing the magnetic resistance and magnetomotive force of each component according to the geometric parameters and the material parameters of the components so as to construct a magnetic network model;
the magnetic circuit model building module: based on a magnetic network model, considering magnetic leakage in a stator slot, magnetic leakage of an adjacent permanent magnet through a long pole shoe and magnetic leakage of the end part position of a stator winding; defining a connection mode among units of the motor; obtaining complete magnetic circuit models at different rotor positions according to different rotor positions corresponding to different moments;
and an iterative calculation module: and carrying out nonlinear iterative calculation on the magnetic permeability of the stator core based on the complete magnetic circuit model to obtain a final convergence result of the magnetic circuit.
A computer storage medium storing a readable program capable of executing the modeling method described above when the program is run.
The invention has the beneficial effects that:
1. the magnetic circuit modeling method of the double-stator single-rotor axial permanent magnet motor can be used for efficiently establishing a magnetic network model of the motor;
2. by using the magnetic circuit modeling method of the double-stator single-rotor axial permanent magnet motor, provided by the invention, a magnetic network model of the motor can be established in the initial stage of motor design so as to quickly obtain the main performance of the motor, and the method has the characteristics of intuitiveness and easiness in implementation;
3. the motor magnetic network model established by the method for modeling the magnetic circuit of the double-stator single-rotor axial permanent magnet motor provided by the invention considers the leakage magnetic circuit of the stator winding end part position, can obtain the average electromagnetic torque of the motor under different load conditions, and is more accurate compared with a two-dimensional finite element model equivalent to a linear motor.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described, and it will be obvious to those skilled in the art that other drawings can be obtained according to these drawings without inventive effort.
FIG. 1 is a flow chart of a modeling method of the present invention;
FIG. 2 is a topological structure diagram of a 24-slot 20-pole double-stator single-rotor axial permanent magnet motor in the invention;
FIG. 3 is a diagram of a stator magnetic network model of the present invention;
FIG. 4 is a diagram of a rotor magnetic network model of the present invention;
FIG. 5 is a schematic diagram of an air gap division model at a first rotor position of the present invention;
FIG. 6 is a diagram of an air gap division model at a second rotor position of the present invention;
FIG. 7 is a diagram of an air gap division model at a third rotor position of the present invention;
FIG. 8 is a magnetic flux leakage model diagram of stator winding end position;
FIG. 9 is a diagram of a complete motor magnetic network model of the present invention;
FIG. 10 is a graph of the magnetic circuit method of the rotor at 0 electrical angle compared with the equivalent two-dimensional finite element axial air gap flux density of a 1-layer linear motor under no-load condition, and the harmonic analysis result;
FIG. 11 is a graph of the magnetic circuit method of the rotor at 15 electrical angle positions compared with the equivalent two-dimensional finite element axial air gap flux density of a 1-layer linear motor under no-load condition and the harmonic analysis result;
FIG. 12 is a graph comparing the magnetic circuit method, the equivalent two-dimensional finite element and three-dimensional finite element analysis of the 1-layer linear motor, and the empty-load phase flux linkage curve and the harmonic analysis result thereof under the condition that the motor operates at 750 rpm;
FIG. 13 is a graph comparing the results of the magnetic circuit method, the equivalent two-dimensional finite element and three-dimensional finite element analysis of the 1-layer linear motor, and the no-load opposite potential curve and the harmonic analysis thereof under the condition that the motor operates at 750 rpm;
fig. 14 is a graph of the average electromagnetic torque contrast curve of the equivalent two-dimensional finite element and three-dimensional finite element of a magnetic circuit method and a 1-layer linear motor, and the relative error between the magnetic circuit method and the three-dimensional finite element result when the slot current density is less than 30 amperes per square millimeter.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
As shown in fig. 2, the components of the analysis motor are typical common components in the motor, such as a stator, a rotor and windings, wherein the rotor is formed by embedding permanent magnets into non-magnetic conductive materials, an air gap exists between the stator and the rotor, and the motor uses three-phase symmetrical alternating current windings. In addition, the technology related to the invention is not only applicable to windings of this type, but also to windings of single-phase, multiphase, asymmetric motors.
In addition, the technology related to the invention is not only suitable for the motor, but also suitable for the topology of all the double-stator single-rotor axial permanent magnet motors with the same permanent magnet and stator width.
As shown in fig. 1, the method for modeling the magnetic circuit of the double-stator single-rotor axial permanent magnet motor comprises the following steps:
s1, defining geometric parameters and material parameters of a motor part;
the geometrical parameters of the motor components include: geometric parameters of stator, rotor, windings, geometric position, air gap, etc.;
the material parameters of the motor component include: the stator and rotor materials have electromagnetic properties, such as magnetic permeability, magnetizing strength and the like;
taking the middle section of the axial motor to be unfolded into a linear motor, and obtaining the geometric parameters and the material parameters of the linear motor; and calculating the length and the area of the linear motor in the magnetic flux direction according to the geometric parameters of the stator, the rotor and the air gap of the linear motor.
S2, dividing nodes of the motor parts, and expressing the magnetic resistance and magnetomotive force of each part according to the geometric parameters and the material parameters of the parts so as to construct a magnetic network model;
the motor part node division comprises nodes for dividing stators, rotors and air gaps, is used for defining the number of nodes of the magnetic network model, and divides the motor into proper units.
The magnetic resistance is divided into rectangular columnar magnetic resistance, columnar magnetic resistance combined by rectangle and semicircle and annular columnar magnetic resistance;
wherein, for rectangular columnar magnetic resistance, magnetic resistance R in the magnetic flux direction thereof 1 It can be expressed as that,
wherein L represents the magnetic path length of a unit of the motor in the magnetic flux direction, mu 0 Is vacuum permeability, mu r For a certain motor component relative permeability, S represents the area through which the flux path passes;
for a rectangular and semicircular combined columnar reluctance, the reluctance R in the direction of the magnetic flux thereof 2 It can be expressed as that,
wherein R is h Represents the radius of a semicircle, L depth Representing the thickness of the linear motor; x represents an integral variable of the radius of the arc;
L depth can be expressed as:
L depth ＝R o -R i (3)
wherein R is i ，R o The inner diameter and the outer diameter of the axial motor are respectively;
for annular columnar magnetoresistance, magnetoresistance R in the direction of magnetic flux thereof 3 It can be expressed as that,
wherein R is m ，R n Respectively the inner diameter and the outer diameter of the circular ring, wherein theta represents the radian swept by the circular ring;
for the three-phase armature magnetomotive force, the axial magnetomotive force F thereof abc Can be expressed as:
wherein N is c Indicating the number of turns, i a 、i b 、i c Respectively representing instantaneous values of three-phase currents;
for the permanent magnetic magnetomotive force, its axial magnetomotive force F pm Can be expressed as:
F pm ＝H c H pm (6)
wherein H is c Indicating the magnetization intensity of the permanent magnet, H pm Representing the axial length of the permanent magnet.
In the embodiment, the motor has symmetry in the axial direction, so that a unilateral stator and a unilateral rotor can be taken for modeling; the magnetic circuit model of the stator teeth is shown in fig. 3; wherein R is y Is the reluctance of the stator yoke, R t Is the reluctance of the stator teeth, R ss Is the magnetic resistance between the two pole shoes, R s Is the leakage magnetic resistance in the groove. R is R sh Is pole shoe magnetic resistance, the bottom of tooth is formed from R tt Modeling, F abc Is armature magnetomotive force F pm Is the magnetomotive force of a permanent magnet.
The magnetic network model of the permanent magnet is shown in fig. 4, and end leakage exists at the inner diameter and the outer diameter of the permanent magnet. Air gap division is a key part of magnetic circuit design, and the utilization range of the air gap is increased due to the existence of pole shoes; firstly, constructing a basic air gap reluctance frame according to the relative positions of the permanent magnets and the stator and the opposite areas; as shown in fig. 5 to 7, R is divided ap (annular air gap leakage reluctance between adjacent permanent magnets) and R pl Outside (parallelogram air gap leakage reluctance between adjacent permanent magnets), the air gap flux path is considered to be axial and has a high degree of accuracy. The stator winding end position magnetic leakage model is shown in fig. 8, and according to the electromagnetic field theoretical analysis and the finite element magnetic field calculation result, the stator winding end position magnetic leakage resistance is defined as two annular columnar magnetic resistances, and can be regarded as being connected in parallel with the in-slot magnetic leakage resistance.
S3, based on the magnetic circuit division of the S2 and the magnetic resistance or magnetomotive force of each component, on the basis of considering the magnetic leakage in the stator slot, the magnetic leakage of the adjacent permanent magnets through the long pole shoes and the magnetic leakage of the end parts of the permanent magnets, modeling the magnetic leakage circuit at the end parts of the stator winding so as to simulate the actual magnetic circuit under heavy load to define the connection mode among all units of the motor; obtaining complete magnetic circuit models at different rotor positions according to different rotor positions corresponding to different moments;
in this embodiment, the complete magnetic circuit model is shown in fig. 9; the model magnetic resistance sizeThe number of nodes and the connection mode have a great relation with the rotor position, and the magnetic circuit is different according to the rotor position. Since the example of the axial permanent magnet motor with the double stators and the single rotor is of a 24-slot/20-pole structure, the magnetic circuit topology structure of the rotor is the same (armature magnetomotive force is different) when the rotor rotates for 3 mechanical angles (30 electrical angles); in this embodiment, a magnetic circuit topology structure is established at two positions, R el The presence or absence of (stator winding end position leakage reluctance) reflects whether the effect of end flux leakage is considered.
S4, based on the complete magnetic circuit model obtained in the step S3, nonlinear iterative computation is carried out on the magnetic permeability of the stator core in the magnetic circuit, and a final convergence result of the magnetic circuit is obtained;
in a magnetic circuit model of the double-stator axial permanent magnet motor, the magnetic permeability of a stator core is nonlinear; in order to obtain the final convergence result of the magnetic circuit, nonlinear iterative computation is required;
the step of nonlinear iterative computation includes:
s41, initializing parameters of an equivalent magnetic circuit, giving initial magnetic permeability, and calculating initial magnetic resistance of a stator core; it should be noted that the initial permeability should be as close to the actual value as possible to reduce the number of iterations;
s42, performing magnetic circuit calculation to obtain magnetic flux phi of each position;
s43, leading out the magnetic flux density B of each position to a working area, and obtaining new magnetic permeability through a B-H curve of a common material;
s44, calculating the relative error between the current magnetic permeability and the magnetic permeability of the last iteration, if the relative error is larger than a given error delta, updating the magnetic permeability and continuing iteration until the relative error is smaller than delta, namely, the calculation result is converged, recording the current magnetic flux and the magnetic flux density at the moment, and finishing the calculation of the network iteration model of the motor.
Permeability may be iterated by:
μ k ＝p 1 μ k +p 2 μ k-1 (7)
wherein mu is k Sum mu k-1 The current permeability and the permeability of the last iteration, p 2 And p 1 Is the current magnetic permeability coefficient and the magnetic permeability coefficient after the last iteration, p 1 +p 2 ＝1；p 1 The larger the convergence, the better the convergence, the slower the convergence speed;
the relative error convergence judgment expression is:
according to the established double-stator single-rotor axial permanent magnet motor magnetic circuit model which takes magnetic leakage at the end part of the stator winding into account, referring to the operation working condition of the motor, calculating the average electromagnetic torque of the motor corresponding to different slot current densities, and comparing with a finite element result;
as shown in fig. 10 and 11 (a) and (b), in the no-load condition, the magnetic circuit method of the rotor at the positions of 0 electric angle and 15 electric angles is compared with the axial air gap flux density of the two-dimensional finite element equivalent to the 1-layer linear motor, and according to waveforms, it can be seen that the magnetic circuit method and the finite element result have certain errors in the no-load air gap flux density at certain positions, but the amplitude accuracy of the harmonic wave of each order is higher.
As shown in fig. 12, fig. 12 (a) and (b) are respectively a magnetic circuit method of the motor running at 750 rpm, equivalent two-dimensional finite element and three-dimensional finite element curves of the 1-layer linear motor, and comparison diagrams of harmonic analysis results of the curves, the fundamental wave amplitude error of which is about 4.3%.
As shown in fig. 13, (a) and (b) are respectively a magnetic circuit method under the condition that the motor operates at 750 rpm, an equivalent two-dimensional finite element and three-dimensional finite element curve between the two-dimensional finite element analysis and the three-dimensional finite element analysis, and a harmonic analysis result comparison diagram thereof, it can be obtained that the fundamental wave amplitude error between the magnetic circuit method and the three-dimensional finite element analysis result is about 4.0%, and the back electromotive force harmonic amplitude calculated by the magnetic circuit method is slightly larger. From fig. 12 and 13, it can be seen that the present invention predicts the no-load performance of the motor more accurately.
As shown in fig. 14, fig. 14 is a graph of magnetic circuit method, two-dimensional finite element equivalent to 1-layer linear motor, three-dimensional finite element average electromagnetic torque contrast curve, and relative error between the former two and three-dimensional finite element results. As can be seen from fig. 14, the average electromagnetic torque error between the magnetic circuit method and the three-dimensional finite element is controlled to be substantially within 5% when the cell current density is within 30 amperes per square millimeter. As can be seen from fig. 14, the present invention still has a more accurate torque prediction capability under the condition of high load of the axial motor with double stators and single rotor.
In the description of the present specification, the descriptions of the terms "one embodiment," "example," "specific example," and the like, mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing has shown and described the basic principles, principal features and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and that the above embodiments and descriptions are merely illustrative of the principles of the present invention, and various changes and modifications may be made without departing from the spirit and scope of the invention, which is defined in the appended claims.
Claims (10)
1. The modeling method for the magnetic circuit of the double-stator single-rotor axial permanent magnet motor is characterized by comprising the following steps of:
defining geometrical parameters and material parameters of the motor parts;
dividing nodes of motor components, and expressing the magnetic resistance and magnetomotive force of each component according to the geometric parameters and the material parameters of the components so as to construct a magnetic network model;
based on a magnetic network model, the connection mode among all units of the motor is defined by considering the magnetic leakage in a stator slot, the magnetic leakage of an adjacent permanent magnet through a long pole shoe and the magnetic leakage of the end part position of a stator winding; obtaining complete magnetic circuit models at different rotor positions according to different rotor positions corresponding to different moments;
based on the complete magnetic circuit model, nonlinear iterative calculation is carried out on the magnetic permeability of the stator core in the magnetic circuit, and a final convergence result of the magnetic circuit is obtained.
2. The method for modeling a magnetic circuit of a dual stator single rotor axial permanent magnet machine of claim 1, wherein the geometric parameters of the machine component include: geometric parameters of stator, rotor and windings, geometric position, and air gap;
the material parameters of the motor component include: the stator and rotor materials have electromagnetic properties including permeability and magnetizing strength.
3. The method of modeling a magnetic circuit of a dual stator single rotor axial permanent magnet machine of claim 1, wherein the node division of the machine component comprises nodes dividing the stator, rotor, and air gap.
4. The modeling method for the magnetic circuit of the double-stator single-rotor axial permanent magnet motor according to claim 1, wherein the magnetic resistances are divided into rectangular columnar magnetic resistances, columnar magnetic resistances combined by rectangular and semicircular magnetic resistances and annular columnar magnetic resistances;
wherein the magnetic resistance R in the rectangular columnar magnetic resistance magnetic flux direction 1 It can be expressed as that,
wherein L represents the magnetic path length of a unit of the motor in the magnetic flux direction, mu 0 Is vacuum permeability, mu r For a certain motor component relative permeability, S represents the area through which the flux path passes;
reluctance R in columnar reluctance magnetic flux direction of rectangular and semicircular combination 2 It can be expressed as that,
wherein R is h Represents the radius of a semicircle, L depth The thickness of the linear motor is represented, and x represents an integral variable of the radius of the arc;
L depth can be expressed as:
L depth ＝R o -R i (3)
wherein R is i ，R o The inner diameter and the outer diameter of the axial motor are respectively;
reluctance R in annular columnar reluctance flux direction 3 Can be expressed as:
wherein R is m ，R n Respectively the inner diameter and the outer diameter of the circular ring, and theta represents the radian swept by the circular ring.
5. The method for modeling a magnetic circuit of a double-stator single-rotor axial permanent magnet motor according to claim 1, wherein the axial magnetomotive force F of the magnetomotive force of the three-phase armature abc Can be expressed as:
wherein N is c Indicating the number of turns, i a 、i b 、i c Respectively representing instantaneous values of three-phase currents;
its axial magnetomotive force F of permanent magnetic magnetomotive force pm Can be expressed as:
F pm ＝H c H pm (6)
wherein H is c Indicating the magnetization intensity of the permanent magnet, H pm Representing the axial length of the permanent magnet.
6. The method for modeling a magnetic circuit of a dual-stator single-rotor axial permanent magnet motor according to claim 1, wherein the step of nonlinear iterative calculation of the magnetic permeability of the stator core comprises:
s41, initializing parameters of an equivalent magnetic circuit, giving initial magnetic permeability, and calculating initial magnetic resistance of a stator core; it should be noted that;
s42, performing magnetic circuit calculation to obtain magnetic flux phi of each position;
s43, leading out the magnetic flux density B of each position to a working area, and obtaining new magnetic permeability through a B-H curve of a common material;
s44, calculating the relative error between the current magnetic permeability and the magnetic permeability of the last iteration, if the relative error is larger than a given error delta, updating the magnetic permeability and continuing iteration until the relative error is smaller than delta, namely, the calculation result is converged, recording the current magnetic flux and the magnetic flux density at the moment, and finishing the calculation of the network iteration model of the motor.
7. The method for modeling a magnetic circuit of a dual-stator single-rotor axial permanent magnet motor according to claim 6, wherein the iterative expression of the magnetic permeability is:
μ k ＝p 1 μ k +p 2 μ k-1 (7)
wherein mu is k Sum mu k-1 The current permeability and the permeability of the last iteration, p 2 And p 1 Is the current magnetic permeability coefficient and the magnetic permeability coefficient after the last iteration, p 1 +p 2 ＝1。
8. The method for modeling a magnetic circuit of a dual-stator single-rotor axial permanent magnet motor according to claim 7, wherein the relative error convergence judgment expression is:
9. the utility model provides a two stator single rotor axial permanent magnet machine magnetic circuit modeling system which characterized in that includes
Parameter definition module: defining geometrical parameters and material parameters of the motor parts;
the magnetic network model building module: dividing nodes of motor components, and expressing the magnetic resistance and magnetomotive force of each component according to the geometric parameters and the material parameters of the components so as to construct a magnetic network model;
the magnetic circuit model building module: based on a magnetic network model, considering magnetic leakage in a stator slot, magnetic leakage of an adjacent permanent magnet through a long pole shoe and magnetic leakage of the end part position of a stator winding; defining a connection mode among units of the motor; obtaining complete magnetic circuit models at different rotor positions according to different rotor positions corresponding to different moments;
and an iterative calculation module: and carrying out nonlinear iterative calculation on the magnetic permeability of the stator core based on the complete magnetic circuit model to obtain a final convergence result of the magnetic circuit.
10. A computer storage medium storing a readable program, characterized in that the modeling method of any of claims 1-8 can be performed when the program is run.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310860251.6A CN116822095A (en) | 2023-07-13 | 2023-07-13 | Magnetic circuit modeling method for double-stator single-rotor axial permanent magnet motor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310860251.6A CN116822095A (en) | 2023-07-13 | 2023-07-13 | Magnetic circuit modeling method for double-stator single-rotor axial permanent magnet motor |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116822095A true CN116822095A (en) | 2023-09-29 |
Family
ID=88118478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310860251.6A Pending CN116822095A (en) | 2023-07-13 | 2023-07-13 | Magnetic circuit modeling method for double-stator single-rotor axial permanent magnet motor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116822095A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117634397A (en) * | 2023-12-01 | 2024-03-01 | 安徽工程大学 | Multi-objective optimization method and system based on two-dimensional equivalent model of axial flux permanent magnet motor |
-
2023
- 2023-07-13 CN CN202310860251.6A patent/CN116822095A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117634397A (en) * | 2023-12-01 | 2024-03-01 | 安徽工程大学 | Multi-objective optimization method and system based on two-dimensional equivalent model of axial flux permanent magnet motor |
| CN117634397B (en) * | 2023-12-01 | 2024-05-28 | 安徽工程大学 | Multi-objective optimization method and system based on two-dimensional equivalent model of axial flux permanent magnet motor |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bonthu et al. | Optimal torque ripple reduction technique for outer rotor permanent magnet synchronous reluctance motors | |
| Moghaddam et al. | Novel high-performance SynRM design method: An easy approach for a complicated rotor topology | |
| CN113094911B (en) | High-power factor design method for magnetic field modulation permanent magnet fault-tolerant motor | |
| Wu et al. | A hybrid field model for open-circuit field prediction in surface-mounted PM machines considering saturation | |
| Ji et al. | A novel linear permanent-magnet vernier machine with improved force performance | |
| Sun et al. | Electromagnetic analysis on novel rotor-segmented axial-field SRM based on dynamic magnetic equivalent circuit | |
| Liu et al. | A new modeling approach for permanent magnet vernier machine with modulation effect consideration | |
| Xu et al. | Analysis of a hybrid rotor permanent magnet motor based on equivalent magnetic network | |
| Zhao et al. | Performance analysis of a coreless axial-flux PMSM by an improved magnetic equivalent circuit model | |
| Jiang et al. | Modeling and analysis of spoke-type permanent magnet vernier machine based on equivalent magnetic network method | |
| CN116822095A (en) | Magnetic circuit modeling method for double-stator single-rotor axial permanent magnet motor | |
| Somesan et al. | Design of a permanent magnet flux-switching machine | |
| He et al. | Multi-layer quasi three-dimensional equivalent model of axial-flux permanent magnet synchronous machine | |
| Candelo-Zuluaga et al. | Customized PMSM design and optimization methodology for water pumping applications | |
| Caruso et al. | An IPMSM torque/weight and torque/moment of inertia ratio optimization | |
| CN114123698A (en) | Method for establishing equivalent magnetic network model of hybrid excitation double-stator magnetic suspension switched reluctance motor | |
| CN117195633B (en) | Magnetic field modeling and analysis method for stator-rotor double-permanent-magnet fault-tolerant motor | |
| Chai et al. | Design of an asymmetric rotor pole for wound field synchronous machines | |
| Mohammadi et al. | A comparison of different models for permanent magnet synchronous machines: Finite element analysis, DQ lumped parameter modeling, and magnetic equivalent circuit | |
| Somesan et al. | Sizing-designing procedure of the permanent magnet flux-switching machine based on a simplified analytical model | |
| CN111654124A (en) | Design method of five-phase permanent magnet fault-tolerant motor with high reluctance torque and high salient pole rate | |
| CN116467890A (en) | Modeling and analysis method for permanent magnet vernier rim propulsion motor based on continuous variable magnetic network | |
| Zhu et al. | Magnetic circuit modeling of dual-armature flux-switching permanent magnet machine | |
| CN115842496A (en) | High-precision magnetic network modeling method suitable for hybrid stator permanent magnet vernier motor | |
| Baek et al. | Optimal design and comparison of stator winding configurations in permanent magnet assisted synchronous reluctance generator |
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 |